On the uncertainty of stream networks derived from elevation data: the error propagation approach

DEM error propagation methodology is extended to the derivation of vector-based objects (stream networks) using geostatistical simulations. First, point sampled elevations are used to fit a variogram model. Next 100 DEM realizations are generated using conditional sequential Gaussian simulation; the stream network map is extracted for each of these realizations, and the collection of stream networks is analyzed to quantify the error propagation. At each grid cell, the probability of the occurrence of a stream and the propagated error are estimated. The method is illustrated using two small data sets: Baranja hill (30 m grid cell size; 16 512 pixels; 6367 sampled elevations), and Zlatibor (30 m grid cell size; 15 000 pixels; 2051 sampled elevations). All computations are run in the open source software for statistical computing R: package geoR is used to fit variogram; package gstat is used to run sequential Gaussian simulation; streams are extracted using the open source GIS SAGA via the RSAGA library. The resulting stream error map (Information entropy of a Bernoulli trial) clearly depicts areas where the extracted stream network is least precise – usually areas of low local relief and slightly convex (0–10 difference from the mean value). In both cases, significant parts of the study area (17.3% for Baranja Hill; 6.2% for Zlatibor) show high error (H>0.5) of locating streams. By correlating the propagated uncertainty of the derived stream network with various land surface parameters sampling of height measurements can be optimized so that delineated streams satisfy the required accuracy level. Such error propagation tool should become a standard functionality in any modern GIS. Remaining issue to be tackled is the computational burden of geostatistical simulations: this framework is at the moment limited to small data sets with several hundreds of points. Scripts and data sets used in this article are available on-line via the www.geomorphometry.org website and can be easily adopted/adjusted to any similar case study

Modelling sea level driven change of Macaronesian archipelago configurations since 120 kyr BP

The MacArthur and Wilson island biogeography theory relates species diversity on islands as the result of equilibrium between extinctions and colonization events which rates depend on island size and isolation. Although island size and isolation can be considered static on ecological timescales (<100 years) they are not static on longer time scales. Since the last million years sea levels fluctuate with a period of ca. 120 kyr between -120 m and up to +10 m MSL (Mean Sea Level). Due to these sea level changes islands have changed in size and ultimately may have drowned or emerged. The rate and degree of their drowning depends on island morphometry and the shape of the sea level change curve. We explore the effects of global sea level cycles on the configuration of archipelagos and volcanic islands of Macaronesia. The results indicate that the islands changed shape considerably during the last 120 kyr. Notably the period between 80 kyr and 15 kyr ago sea levels were at least 80 m lower than present and several islands now isolated were merged or were much larger than present. Recent shrinking of islands due to the sea level rise since the last glacial maximum period (20 kyr BP) led to more than 50% reductions in island size, significant loss of coastal habitat and a significant increase in isolation by the increase of distances between islands and island and continents. Island size reduction must have induced pressures especially on terrestrial insular ecosystems, inducing upward migrations and interspecies competitions, and probable extinctions. The splitting of merged islands must have led to separations of populations leading to gene flow losses for some biota. Present day islands are not representative for the mean island configurationsduring the last Myr but rather represent an anomaly. Islands at present are smallest and mostisolated and this configuration makes the insular biota even more vulnerable to human impact

Spatio-temporal prediction of daily temperatures using time-series of MODIS LST images

A computational framework to generate daily temperature maps using time-series of publicly available MODIS MOD11A2 product Land Surface Temperature (LST) images (1 km resolution; 8-day composites) is illustrated using temperature measurements from the national network of meteorological stations (159) in Croatia. The input data set contains 57,282 ground measurements of daily temperature for the year 2008. Temperature was modeled as a function of latitude, longitude, distance from the sea, elevation, time, insolation, and the MODIS LST images. The original rasters were first converted to principal components to reduce noise and filter missing pixels in the LST images. The residual were next analyzed for spatio-temporal auto-correlation; sum-metric separable variograms were fitted to account for zonal and geometric space-time anisotropy. The final predictions were generated for time-slices of a 3D space-time cube, constructed in the R environment for statistical computing. The results show that the space-time regression model can explain a significant part of the variation in station-data (84%). MODIS LST 8-day (cloud-free) images are unbiased estimator of the daily temperature, but with relatively low precision (±4.1°C); however their added value is that they systematically improve detection of local changes in land surface temperature due to local meteorological conditions and/or active heat sources (urban areas, land cover classes). The results of 10–fold cross-validation show that use of spatio-temporal regression-kriging and incorporation of time-series of remote sensing images leads to significantly more accurate maps of temperature than if plain spatial techniques were used. The average (global) accuracy of mapping temperature was ±2.4°C. The regression-kriging explained 91% of variability in daily temperatures, compared to 44% for ordinary kriging. Further software advancement—interactive space-time variogram exploration and automated retrieval, resampling and filtering of MODIS images—are anticipated

Soil nutrient maps of Sub-Saharan Africa: assessment of soil nutrient content at 250 m spatial resolution using machine learning

Spatial predictions of soil macro and micro-nutrient content across Sub-Saharan Africa at 250 m spatial resolution and for 0–30 cm depth interval are presented. Predictions were produced for 15 target nutrients: organic carbon (C) and total (organic) nitrogen (N), total phosphorus (P), and extractable—phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), sulfur (S), sodium (Na), iron (Fe), manganese (Mn), zinc (Zn), copper (Cu), aluminum (Al) and boron (B). Model training was performed using soil samples from ca. 59,000 locations (a compilation of soil samples from the AfSIS, EthioSIS, One Acre Fund, VitalSigns and legacy soil data) and an extensive stack of remote sensing covariates in addition to landform, lithologic and land cover maps. An ensemble model was then created for each nutrient from two machine learning algorithms

Development of global soil information facilities

ISRIC - World Soil Information has a mandate to serve the international community as custodian of global soil information and to increase awareness and understanding of the role of soils in major global issues. To adapt to the current demand for soil information, ISRIC is updating its enterprise data management system, including procedures for registering acquired data, such as lineage, versioning, quality assessment, and control. Data can be submitted, queried, and analysed using a growing range of web-based services - ultimately aiming at full and open exchange of data, metadata, and products - through the ICSU-accredited World Data Centre for Soils

OEMC D2.1 Report "Stakeholder Committee and Open- Earth-Monitor Design" workshop

This deliverable of the Open-Earth-Monitor project describes the approach taken to compile and categorize the project's stakeholders, and provides recommendations for the future stakeholder interactions based on the results of a survey from the OEMC design workshop, which took place during the project's kick-off meeting in July, 2022

Comparison of Different Methods of Automated Landform Classification at the Drainage Basin Scale: Examples from the Southern Italy

In this work, we tested the reliability of two different methods of automated landform classification (ACL) in three geological domains of the southern Italian chain with contrasting morphological features. ACL maps deriving from the TPI-based (topographic position index) algorithm are strictly dependent to the search input parameters and they are not able to fully capture landforms of different size. Geomorphons-based classification has shown a higher potential and can represent a powerful method of ACL, although it should be improved with the introduction of additional DEM-based parameters for the extraction of landform classe

Quo vadis, smallholder forest landscape?: An introduction to the LPB-RAP model

The impacts of the Anthropocene on climate and biodiversity pose societal and ecological problems that may only be solved by ecosystem restoration. Local to regional actions are required, which need to consider the prevailing present and future conditions of a certain landscape extent. Modeling approaches can be of help to support management efforts and to provide advice to policy making. We present stage one of the LaForeT-PLUC-BE model (Landscape Forestry in the Tropics–PCRaster Land Use Change–Biogeographic & Economic model; in short: LPB) and its thematic expansion module RAP (Restoration Areas Potentials). LPB-RAP is a high-resolution pixel-based scenario tool that relies on a range of explicit land use types (LUTs) to describe various forest types and the environment. It simulates and analyzes future landscape configurations under consideration of climate, population and land use change long-term. Simulated Land Use Land Cover Change (LULCC) builds on dynamic, probabilistic modeling incorporating climatic and anthropogenic determinants as well as restriction parameters to depict a sub-national regional smallholder-dominated forest landscape. The model delivers results for contrasting scenario settings by simulating without and with potential Forest and Landscape Restoration (FLR) measures. FLR potentials are depicted by up to five RAP-LUTs. The model builds on user-defined scenario inputs, such as the Shared Socioeconomic Pathways (SSP) and Representative Concentration Pathways (RCP). Model application is here exemplified for the SSP2-RCP4.5 scenario in the time frame 2018–2100 on the hectare scale in annual resolution using Esmeraldas province, Ecuador, as a case study area. The LPB-RAP model is a novel, heuristic Spatial Decision Support System (SDSS) tool for smallholder-dominated forest landscapes, supporting near-time top-down planning measures with long-term bottom-up modeling. Its application should be followed up by FLR on-site investigations and stakeholder participation across all involved scales