526 Mathematische Geografie
Refine
Has Fulltext
- yes (114)
Is part of the Bibliography
- yes (114) (remove)
Year of publication
Document Type
- Journal article (111)
- Doctoral Thesis (3)
Keywords
- remote sensing (26)
- Sentinel-2 (13)
- time series (10)
- Google Earth Engine (9)
- MODIS (9)
- Sentinel-1 (8)
- earth observation (8)
- machine learning (8)
- Earth observation (6)
- Landsat (6)
Institute
- Institut für Geographie und Geologie (114) (remove)
EU-Project number / Contract (GA) number
- 227159 (1)
- 818182 (1)
- LIFE12 BIO/AT/000143 (1)
- LIFE20 NAT/AT/000049 (1)
The Essential Climate Variable (ECV) Permafrost is currently undergoing strong changes due to rising ground and air temperatures. Surface movement, forming characteristic landforms such as rock glaciers, is one key indicator for mountain permafrost. Monitoring this movement can indicate ongoing changes in permafrost; therefore, rock glacier velocity (RGV) has recently been added as an ECV product. Despite the increased understanding of rock glacier dynamics in recent years, most observations are either limited in terms of the spatial coverage or temporal resolution. According to recent studies, Sentinel-1 (C-band) Differential SAR Interferometry (DInSAR) has potential for monitoring RGVs at high spatial and temporal resolutions. However, the suitability of DInSAR for the detection of heterogeneous small-scale spatial patterns of rock glacier velocities was never at the center of these studies. We address this shortcoming by generating and analyzing Sentinel-1 DInSAR time series over five years to detect small-scale displacement patterns of five high alpine permafrost environments located in the Central European Alps on a weekly basis at a range of a few millimeters. Our approach is based on a semi-automated procedure using open-source programs (SNAP, pyrate) and provides East-West displacement and elevation change with a ground sampling distance of 5 m. Comparison with annual movement derived from orthophotos and unpiloted aerial vehicle (UAV) data shows that DInSAR covers about one third of the total movement, which represents the proportion of the year suited for DInSAR, and shows good spatial agreement (Pearson R: 0.42–0.74, RMSE: 4.7–11.6 cm/a) except for areas with phase unwrapping errors. Moreover, the DInSAR time series unveils spatio-temporal variations and distinct seasonal movement dynamics related to different drivers and processes as well as internal structures. Combining our approach with in situ observations could help to achieve a more holistic understanding of rock glacier dynamics and to assess the future evolution of permafrost under changing climatic conditions.
Introduction: Grasslands cover one third of the agricultural area in Germany and are mainly used for fodder production. However, grasslands fulfill many other ecosystem functions, like carbon storage, water filtration and the provision of habitats. In Germany, grasslands are mown and/or grazed multiple times during the year. The type and timing of management activities and the use intensity vary strongly, however co-determine grassland functions. Large-scale spatial information on grassland activities and use intensity in Germany is limited and not openly provided. In addition, the cause for patterns of varying mowing intensity are usually not known on a spatial scale as data on the incentives of farmers behind grassland management decisions is not available.
Methods: We applied an algorithm based on a thresholding approach utilizing Sentinel-2 time series to detect grassland mowing events to investigate mowing dynamics in Germany in 2018–2021. The detected mowing events were validated with an independent dataset based on the examination of public webcam images. We analyzed spatial and temporal patterns of the mowing dynamics and relationships to climatic, topographic, soil or socio-political conditions.
Results: We found that most intensively used grasslands can be found in southern/south-eastern Germany, followed by areas in northern Germany. This pattern stays the same among the investigated years, but we found variations on smaller scales. The mowing event detection shows higher accuracies in 2019 and 2020 (F1 = 0.64 and 0.63) compared to 2018 and 2021 (F1 = 0.52 and 0.50). We found a significant but weak (R2 of 0–0.13) relationship for a spatial correlation of mowing frequency and climate as well as topographic variables for the grassland areas in Germany. Further results indicate a clear value range of topographic and climatic conditions, characteristic for intensive grassland use. Extensive grassland use takes place everywhere in Germany and on the entire spectrum of topographic and climatic conditions in Germany. Natura 2000 grasslands are used less intensive but this pattern is not consistent among all sites.
Discussion: Our findings on mowing dynamics and relationships to abiotic and socio-political conditions in Germany reveal important aspects of grassland management, including incentives of farmers.
A circum-Arctic monitoring framework for quantifying annual erosion rates of permafrost coasts
(2023)
This study demonstrates a circum-Arctic monitoring framework for quantifying annual change of permafrost-affected coasts at a spatial resolution of 10 m. Frequent cloud coverage and challenging lighting conditions, including polar night, limit the usability of optical data in Arctic regions. For this reason, Synthetic Aperture RADAR (SAR) data in the form of annual median and standard deviation (sd) Sentinel-1 (S1) backscatter images covering the months June–September for the years 2017–2021 were computed. Annual composites for the year 2020 were hereby utilized as input for the generation of a high-quality coastline product via a Deep Learning (DL) workflow, covering 161,600 km of the Arctic coastline. The previously computed annual S1 composites for the years 2017 and 2021 were employed as input data for the Change Vector Analysis (CVA)-based coastal change investigation. The generated DL coastline product served hereby as a reference. Maximum erosion rates of up to 67 m per year could be observed based on 400 m coastline segments. Overall highest average annual erosion can be reported for the United States (Alaska) with 0.75 m per year, followed by Russia with 0.62 m per year. Out of all seas covered in this study, the Beaufort Sea featured the overall strongest average annual coastal erosion of 1.12 m. Several quality layers are provided for both the DL coastline product and the CVA-based coastal change analysis to assess the applicability and accuracy of the output products. The predicted coastal change rates show good agreement with findings published in previous literature. The proposed methods and data may act as a valuable tool for future analysis of permafrost loss and carbon emissions in Arctic coastal environments.
Monitoring forest conditions is an essential task in the context of global climate change to preserve biodiversity, protect carbon sinks and foster future forest resilience. Severe impacts of heatwaves and droughts triggering cascading effects such as insect infestation are challenging the semi-natural forests in Germany. As a consequence of repeated drought years since 2018, large-scale canopy cover loss has occurred calling for an improved disturbance monitoring and assessment of forest structure conditions. The present study demonstrates the potential of complementary remote sensing sensors to generate wall-to-wall products of forest structure for Germany. The combination of high spatial and temporal resolution imagery from Sentinel-1 (Synthetic Aperture Radar, SAR) and Sentinel-2 (multispectral) with novel samples on forest structure from the Global Ecosystem Dynamics Investigation (GEDI, LiDAR, Light detection and ranging) enables the analysis of forest structure dynamics. Modeling the three-dimensional structure of forests from GEDI samples in machine learning models reveals the recent changes in German forests due to disturbances (e.g., canopy cover degradation, salvage logging). This first consistent data set on forest structure for Germany from 2017 to 2022 provides information of forest canopy height, forest canopy cover and forest biomass and allows estimating recent forest conditions at 10 m spatial resolution. The wall-to-wall maps of the forest structure support a better understanding of post-disturbance forest structure and forest resilience.
Satellite-derived land surface temperature dynamics in the context of global change — a review
(2023)
Satellite-derived Land Surface Temperature (LST) dynamics have been increasingly used to study various geophysical processes. This review provides an extensive overview of the applications of LST in the context of global change. By filtering a selection of relevant keywords, a total of 164 articles from 14 international journals published during the last two decades were analyzed based on study location, research topic, applied sensor, spatio-temporal resolution and scale and employed analysis methods. It was revealed that China and the USA were the most studied countries and those that had the most first author affiliations. The most prominent research topic was the Surface Urban Heat Island (SUHI), while the research topics related to climate change were underrepresented. MODIS was by far the most used sensor system, followed by Landsat. A relatively small number of studies analyzed LST dynamics on a global or continental scale. The extensive use of MODIS highly determined the study periods: A majority of the studies started around the year 2000 and thus had a study period shorter than 25 years. The following suggestions were made to increase the utilization of LST time series in climate research: The prolongation of the time series by, e.g., using AVHRR LST, the better representation of LST under clouds, the comparison of LST to traditional climate change measures, such as air temperature and reanalysis variables, and the extension of the validation to heterogenous sites.
The increasing availability and variety of global satellite products and the rapid development of new algorithms has provided great potential to generate a new level of data with different spatial, temporal, and spectral resolutions. However, the ability of these synthetic spatiotemporal datasets to accurately map and monitor our planet on a field or regional scale remains underexplored. This study aimed to support future research efforts in estimating crop yields by identifying the optimal spatial (10 m, 30 m, or 250 m) and temporal (8 or 16 days) resolutions on a regional scale. The current study explored and discussed the suitability of four different synthetic (Landsat (L)-MOD13Q1 (30 m, 8 and 16 days) and Sentinel-2 (S)-MOD13Q1 (10 m, 8 and 16 days)) and two real (MOD13Q1 (250 m, 8 and 16 days)) NDVI products combined separately to two widely used crop growth models (CGMs) (World Food Studies (WOFOST), and the semi-empiric Light Use Efficiency approach (LUE)) for winter wheat (WW) and oil seed rape (OSR) yield forecasts in Bavaria (70,550 km\(^2\)) for the year 2019. For WW and OSR, the synthetic products’ high spatial and temporal resolution resulted in higher yield accuracies using LUE and WOFOST. The observations of high temporal resolution (8-day) products of both S-MOD13Q1 and L-MOD13Q1 played a significant role in accurately measuring the yield of WW and OSR. For example, L- and S-MOD13Q1 resulted in an R\(^2\) = 0.82 and 0.85, RMSE = 5.46 and 5.01 dt/ha for WW, R\(^2\) = 0.89 and 0.82, and RMSE = 2.23 and 2.11 dt/ha for OSR using the LUE model, respectively. Similarly, for the 8- and 16-day products, the simple LUE model (R\(^2\) = 0.77 and relative RMSE (RRMSE) = 8.17%) required fewer input parameters to simulate crop yield and was highly accurate, reliable, and more precise than the complex WOFOST model (R\(^2\) = 0.66 and RRMSE = 11.35%) with higher input parameters. Conclusively, both S-MOD13Q1 and L-MOD13Q1, in combination with LUE, were more prominent for predicting crop yields on a regional scale than the 16-day products; however, L-MOD13Q1 was advantageous for generating and exploring the long-term yield time series due to the availability of Landsat data since 1982, with a maximum resolution of 30 m. In addition, this study recommended the further use of its findings for implementing and validating the long-term crop yield time series in different regions of the world.
Rapid and accurate yield estimates at both field and regional levels remain the goal of sustainable agriculture and food security. Hereby, the identification of consistent and reliable methodologies providing accurate yield predictions is one of the hot topics in agricultural research. This study investigated the relationship of spatiotemporal fusion modelling using STRAFM on crop yield prediction for winter wheat (WW) and oil-seed rape (OSR) using a semi-empirical light use efficiency (LUE) model for the Free State of Bavaria (70,550 km\(^2\)), Germany, from 2001 to 2019. A synthetic normalised difference vegetation index (NDVI) time series was generated and validated by fusing the high spatial resolution (30 m, 16 days) Landsat 5 Thematic Mapper (TM) (2001 to 2012), Landsat 7 Enhanced Thematic Mapper Plus (ETM+) (2012), and Landsat 8 Operational Land Imager (OLI) (2013 to 2019) with the coarse resolution of MOD13Q1 (250 m, 16 days) from 2001 to 2019. Except for some temporal periods (i.e., 2001, 2002, and 2012), the study obtained an R\(^2\) of more than 0.65 and a RMSE of less than 0.11, which proves that the Landsat 8 OLI fused products are of higher accuracy than the Landsat 5 TM products. Moreover, the accuracies of the NDVI fusion data have been found to correlate with the total number of available Landsat scenes every year (N), with a correlation coefficient (R) of +0.83 (between R\(^2\) of yearly synthetic NDVIs and N) and −0.84 (between RMSEs and N). For crop yield prediction, the synthetic NDVI time series and climate elements (such as minimum temperature, maximum temperature, relative humidity, evaporation, transpiration, and solar radiation) are inputted to the LUE model, resulting in an average R\(^2\) of 0.75 (WW) and 0.73 (OSR), and RMSEs of 4.33 dt/ha and 2.19 dt/ha. The yield prediction results prove the consistency and stability of the LUE model for yield estimation. Using the LUE model, accurate crop yield predictions were obtained for WW (R\(^2\) = 0.88) and OSR (R\(^2\) = 0.74). Lastly, the study observed a high positive correlation of R = 0.81 and R = 0.77 between the yearly R\(^2\) of synthetic accuracy and modelled yield accuracy for WW and OSR, respectively.
On a daily basis, political decisions are made, often with their full extent of impact being unclear. Not seldom, the decisions and policy measures implemented result in direct or indirect unintended negative impacts, such as on the natural environment, which can vary in time, space, nature, and severity. To achieve a more sustainable world with equitable societies requires fundamental rethinking of our policymaking. It calls for informed decision making and a monitoring of political impact for which evidence-based knowledge is necessary. The most powerful tool to derive objective and systematic spatial information and, thus, add to transparent decisions is remote sensing (RS). This review analyses how spaceborne RS is used by the scientific community to provide evidence for the policymaking process. We reviewed 194 scientific publications from 2015 to 2020 and analysed them based on general insights (e.g., study area) and RS application-related information (e.g., RS data and products). Further, we classified the studies according to their degree of science–policy integration by determining their engagement with the political field and their potential contribution towards four stages of the policy cycle: problem identification/knowledge building, policy formulation, policy implementation, and policy monitoring and evaluation. Except for four studies, we found that studies had not directly involved or informed the policy field or policymaking process. Most studies contributed to the stage problem identification/knowledge building, followed by ex post policy impact assessment. To strengthen the use of RS for policy-relevant studies, the concept of the policy cycle is used to showcase opportunities of RS application for the policymaking process. Topics gaining importance and future requirements of RS at the science–policy interface are identified. If tackled, RS can be a powerful complement to provide policy-relevant evidence to shed light on the impact of political decisions and thus help promote sustainable development from the core.
The positive phase of the subtropical Indian Ocean dipole (SIOD) is one of the climatic modes in the subtropical southern Indian Ocean that influences the austral summer inter-annual rainfall variability in parts of southern Africa. This paper examines austral summer rain-bearing circulation types (CTs) in Africa south of the equator that are related to the positive SIOD and the dynamics through which specific rainfall regions in southern Africa can be influenced by this relationship. Four austral summer rain-bearing CTs were obtained. Among the four CTs, the CT that featured (i) enhanced cyclonic activity in the southwest Indian Ocean; (ii) positive widespread rainfall anomaly in the southwest Indian Ocean; and (iii) low-level convergence of moisture fluxes from the tropical South Atlantic Ocean, tropical Indian Ocean, and the southwest Indian Ocean, over the south-central landmass of Africa, was found to be related to the positive SIOD climatic mode. The relationship also implies that positive SIOD can be expected to increase the amplitude and frequency of occurrence of the aforementioned CT. The linkage between the CT related to the positive SIOD and austral summer homogeneous regions of rainfall anomalies in Africa south of the equator showed that it is the principal CT that is related to the inter-annual rainfall variability of the south-central regions of Africa, where the SIOD is already known to significantly influence its rainfall variability. Hence, through the large-scale patterns of atmospheric circulation associated with the CT, the SIOD can influence the spatial distribution and intensity of rainfall over the preferred landmass through enhanced moisture convergence.
The July 2021 heavy rainfall episode in parts of Western Europe caused devastating floods, specifically in Germany. This study examines circulation types (CTs) linked to extreme precipitation in Germany. It was investigated if the classified CTs can highlight the anomaly in synoptic patterns that contributed to the unusual July 2021 heavy rainfall in Germany. The North Atlantic Oscillation was found to be the major climatic mode related to the seasonal and inter-annual variations of most of the classified CTs. On average, wet (dry) conditions in large parts of Germany can be linked to westerly (northerly) moisture fluxes. During spring and summer seasons, the mid-latitude cyclone when located over the North Sea disrupts onshore moisture transport from the North Atlantic Ocean by westerlies driven by the North Atlantic subtropical anticyclone. The CT found to have the highest probability of being associated with above-average rainfall in large part of Germany features (i) enhancement and northward track of the cyclonic system over the Mediterranean; (ii) northward track of the North Atlantic anticyclone, further displacing poleward, the mid-latitude cyclone over the North Sea, enabling band of westerly moisture fluxes to penetrate Germany; (iii) cyclonic system over the Baltic Sea coupled with northeast fluxes of moisture to Germany; (iv) and unstable atmospheric conditions over Germany. In 2021, a spike was detected in the amplitude and frequency of occurrence of the aforementioned wet CT suggesting that in addition to the nearly stationary cut-off low over central Europe, during the July flood episode, anomalies in the CT contributed to the heavy rainfall event.
The occurrence of a likely graptolite in lowest Wuliuan strata of the Franconian Forest almost certainly records the oldest known graptolithoid hemichordate in West Gondwana and possibly the oldest graptolite presently known. The fossil is a delicate, erect, apparently unbranched rhabdosome with narrow thecae tentatively assigned to the poorly known genus Ovetograptus of the Dithecodendridae. This report includes an overview of pre-Furongian graptolithoids with slight corrections on the stratigraphic position of earlier reported species.
New U–Pb age and Hf isotope data obtained on detrital zircon grains from Au- and U-bearing Archaean quartz-pebble conglomerates in the Singhbhum Craton, eastern India, specifically the Upper Iron Ore Group in the Badampahar Greenstone Belt and the Phuljhari Formation below the Dhanjori Group provide insights into the zircon provenance and maximum age of sediment deposition. The most concordant, least disturbed \(^{207}\)Pb/\(^{206}\)Pb ages cover the entire range of known magmatic and higher grade metamorphic events in the craton from 3.48 to 3.06 Ga and show a broad maximum between 3.38 and 3.18 Ga. This overlap is also mimicked by Lu–Hf isotope analyses, which returned a wide range in \(_{εHf}\)(t) values from + 6 to − 5, in agreement with the range known from zircon grains in igneous and metamorphic rocks in the Singhbhum Craton. A smaller but distinct age peak centred at 3.06 Ga corresponds to the age of the last major magmatic intrusive event, the emplacement of the Mayurbhanj Granite and associated gabbro, picrite and anorthosite. Thus, these intrusive rocks must form a basement rather than being intrusive into the studied conglomerates as previously interpreted. The corresponding detrital zircon grains all have a subchondritic Hf isotopic composition. The youngest reliable zircon ages of 3.03 Ga in the case of the basal Upper Iron Ore Group in the east of the craton and 3.00 Ga for the Phuljhari Formation set an upper limit on the age of conglomerate sedimentation. Previously published detrital zircon age data from similarly Au-bearing conglomerates in the Mahagiri Quartzite in the Upper Iron Ore Group in the south of the craton gave a somewhat younger maximum age of sedimentation of 2.91 Ga. There, the lower limit on sedimentation is given by an intrusive relationship with a c. 2.8 Ga granite. The time window thus defined for conglomerate deposition on the Singhbhum Craton is almost identical to the age span established for the, in places, Au- and U-rich conglomerates in the Kaapvaal Craton of South Africa: the 2.98–2.78 Ga Dominion Group and Witwatersrand Supergroup in South Africa. Since the recognition of first major concentration of gold on Earth’s surface by microbial activity having taken place at around 2.9 Ga, independent of the nature of the hinterland, the above similarity in age substantially increases the potential for discovering Witwatersrand-type gold and/or uranium deposits on the Singhbhum Craton. Further age constraints are needed there, however, to distinguish between supposedly less fertile (with respect to Au) > 2.9 Ga and more fertile < 2.9 Ga successions.
The effects of drought on tree mortality at forest stands are not completely understood. For assessing their water supply, knowledge of the small-scale distribution of soil moisture as well as its temporal changes is a key issue in an era of climate change. However, traditional methods like taking soil samples or installing data loggers solely collect parameters of a single point or of a small soil volume. Electrical resistivity tomography (ERT) is a suitable method for monitoring soil moisture changes and has rarely been used in forests. This method was applied at two forest sites in Bavaria, Germany to obtain high-resolution data of temporal soil moisture variations. Geoelectrical measurements (2D and 3D) were conducted at both sites over several years (2015–2018/2020) and compared with soil moisture data (matric potential or volumetric water content) for the monitoring plots. The greatest variations in resistivity values that highly correlate with soil moisture data were found in the main rooting zone. Using the ERT data, temporal trends could be tracked in several dimensions, such as the interannual increase in the depth of influence from drought events and their duration, as well as rising resistivity values going along with decreasing soil moisture. The results reveal that resistivity changes are a good proxy for seasonal and interannual soil moisture variations. Therefore, 2D- and 3D-ERT are recommended as comparatively non-laborious methods for small-spatial scale monitoring of soil moisture changes in the main rooting zone and the underlying subsurface of forested sites. Higher spatial and temporal resolution allows a better understanding of the water supply for trees, especially in times of drought.
A fuzzy classification scheme that results in physically interpretable meteorological patterns associated with rainfall generation is applied to classify homogeneous regions of boreal summer rainfall anomalies in Germany. Four leading homogeneous regions are classified, representing the western, southeastern, eastern, and northern/northwestern parts of Germany with some overlap in the central parts of Germany. Variations of the sea level pressure gradient across Europe, e.g., between the continental and maritime regions, is the major phenomenon that triggers the time development of the rainfall regions by modulating wind patterns and moisture advection. Two regional climate models (REMO and CCLM4) were used to investigate the capability of climate models to reproduce the observed summer rainfall regions. Both regional climate models (RCMs) were once driven by the ERA-Interim reanalysis and once by the MPI-ESM general circulation model (GCM). Overall, the RCMs exhibit good performance in terms of the regionalization of summer rainfall in Germany; though the goodness-of-match with the rainfall regions/patterns from observational data is low in some cases and the REMO model driven by MPI-ESM fails to reproduce the western homogeneous rainfall region. Under future climate change, virtually the same leading modes of summer rainfall occur, suggesting that the basic synoptic processes associated with the regional patterns remain the same over Germany. We have also assessed the added value of bias-correcting the MPI-ESM driven RCMs using a simple linear scaling approach. The bias correction does not significantly alter the identification of homogeneous rainfall regions and, hence, does not improve their goodness-of-match compared to the observed patterns, except for the one case where the original RCM output completely fails to reproduce the observed pattern. While the linear scaling method improves the basic statistics of precipitation, it does not improve the simulated meteorological patterns represented by the precipitation regimes.
Cocoa growing is one of the main activities in humid West Africa, which is mainly grown in pure stands. It is the main driver of deforestation and encroachment in protected areas. Cocoa agroforestry systems which have been promoted to mitigate deforestation, needs to be accurately delineated to support a valid monitoring system. Therefore, the aim of this research is to model the spatial distribution of uncertainties in the classification cocoa agroforestry. The study was carried out in Côte d’Ivoire, close to the Taï National Park. The analysis followed three steps (i) image classification based on texture parameters and vegetation indices from Sentinel-1 and -2 data respectively, to train a random forest algorithm. A classified map with the associated probability maps was generated. (ii) Shannon entropy was calculated from the probability maps, to get the error maps at different thresholds (0.2, 0.3, 0.4 and 0.5). Then, (iii) the generated error maps were analysed using a Geographically Weighted Regression model to check for spatial autocorrelation. From the results, a producer accuracy (0.88) and a user’s accuracy (0.91) were obtained. A small threshold value overestimates the classification error, while a larger threshold will underestimate it. The optimal value was found to be between 0.3 and 0.4. There was no evidence of spatial autocorrelation except for a smaller threshold (0.2). The approach differentiated cocoa from other landcover and detected encroachment in forest. Even though some information was lost in the process, the method is effective for mapping cocoa plantations in Côte d’Ivoire.
Performance assessment of CORDEX regional climate models in wind speed simulations over Zambia
(2023)
There is no single solution to cutting emissions, however, renewable energy projects that are backed by rigorous ex-ante assessments play an important role in these efforts. An inspection of literature reveals critical knowledge gaps in the understanding of future wind speed variability across Zambia, thus leading to major uncertainties in the understanding of renewable wind energy potential over the country. Several model performance metrics, both statistical and graphical were used in this study to examine the performance of CORDEX Africa Regional Climate Models (RCMs) in simulating wind speed across Zambia. Results indicate that wind speed is increasing at the rate of 0.006 m s\(^{−1}\) per year. RCA4-GFDL-ESM2M, RCA4-HadGEM2-ES, RCA4-IPSL-CM5A-MR, and RCA4-CSIRO-MK3.6.0 were found to correctly simulate wind speed increase with varying magnitudes on the Sen’s estimator of slope. All the models sufficiently reproduce the annual cycle of wind speed with a steady increase being observed from April reaching its peak around August/September and beginning to drop in October. Apart from RegCM4-MPI-ESM and RegCM4-HadGEM2, the performance of RCMs in simulating spatial wind speed patterns is generally good although they overestimate it by ~ 1 m s\(^{−1}\) in the western and southern provinces of the country. Model performance metrics indicate that with a correlation coefficient of 0.5, a root mean square error of 0.4 m s\(^{−1}\), an RSR value of 7.7 and a bias of 19.9%, RCA4-GFDL-ESM2M outperforms all other models followed by RCA4-HadGEM2, and RCA4-CM5A-MR respectively. These results, therefore, suggest that studies that use an ensemble of RCA4-GFDL-ESM2M, RCA4-HadGEM2, and RCA4-CM5A-MR would yield useful results for informing future renewable wind energy potential in Zambia.
Performance of a regional climate model with interactive vegetation (REMO-iMOVE) over Central Asia
(2022)
The current study evaluates the regional climate model REMO (v2015) and its new version REMO-iMOVE, including interactive vegetation and plant functional types (PFTs), over two Central Asian domains for the period of 2000–2015 at two different horizontal resolutions (0.44° and 0.11°). Various statistical metrices along with mean bias patterns for precipitation, temperature, and leaf area index have been used for the model evaluation. A better representation of the spatial pattern of precipitation is found at 0.11° resolution over most of Central Asia. Regarding the mean temperature, both model versions show a high level of agreement with the validation data, especially at the higher resolution. This also reduces the biases in maximum and minimum temperature. Generally, REMO-iMOVE shows an improvement regarding the temperature bias but produces a larger precipitation bias compared to the REMO conventional version with interannually static vegetation. Since the coupled version is capable to simulate the mean climate of Central Asia like its parent version, both can be used for impact studies and future projections. However, regarding the new vegetation scheme and its spatiotemporal representation exemplified by the leaf area index, REMO-iMOVE shows a clear advantage over REMO. This better simulation is caused by the implementation of more realistic and interactive vegetation and related atmospheric processes which consequently add value to the regional climate model.
The Moroccan locust has been considered one of the most dangerous agricultural pests in the Mediterranean region. The economic importance of its outbreaks diminished during the second half of the 20th century due to a high degree of agricultural industrialization and other human-caused transformations of its habitat. Nevertheless, in Sardinia (Italy) from 2019 on, a growing invasion of this locust species is ongoing, being the worst in over three decades. Locust swarms destroyed crops and pasture lands of approximately 60,000 ha in 2022. Drought, in combination with increasing uncultivated land, contributed to forming the perfect conditions for a Moroccan locust population upsurge. The specific aim of this paper is the quantification of land cover land use (LCLU) influence with regard to the recent locust outbreak in Sardinia using remote sensing data. In particular, the role of untilled, fallow, or abandoned land in the locust population upsurge is the focus of this case study. To address this objective, LCLU was derived from Sentinel-2A/B Multispectral Instrument (MSI) data between 2017 and 2021 using time-series composites and a random forest (RF) classification model. Coordinates of infested locations, altitude, and locust development stages were collected during field observation campaigns between March and July 2022 and used in this study to assess actual and previous land cover situation of these locations. Findings show that 43% of detected locust locations were found on untilled, fallow, or uncultivated land and another 23% within a radius of 100 m to such areas. Furthermore, oviposition and breeding sites are mostly found in sparse vegetation (97%). This study demonstrates that up-to-date remote sensing data and target-oriented analyses can provide valuable information to contribute to early warning systems and decision support and thus to minimize the risk concerning this agricultural pest. This is of particular interest for all agricultural pests that are strictly related to changing human activities within transformed habitats.
A novel method for detecting and delineating coppice trees in UAV images to monitor tree decline
(2022)
Monitoring tree decline in arid and semi-arid zones requires methods that can provide up-to-date and accurate information on the health status of the trees at single-tree and sample plot levels. Unmanned Aerial Vehicles (UAVs) are considered as cost-effective and efficient tools to study tree structure and health at small scale, on which detecting and delineating tree crowns is the first step to extracting varied subsequent information. However, one of the major challenges in broadleaved tree cover is still detecting and delineating tree crowns in images. The frequent dominance of coppice structure in degraded semi-arid vegetation exacerbates this problem. Here, we present a new method based on edge detection for delineating tree crowns based on the features of oak trees in semi-arid coppice structures. The decline severity in individual stands can be analyzed by extracting relevant information such as texture from the crown area. Although the method presented in this study is not fully automated, it returned high performances including an F-score = 0.91. Associating the texture indices calculated in the canopy area with the phenotypic decline index suggested higher correlations of the GLCM texture indices with tree decline at the tree level and hence a high potential to be used for subsequent remote-sensing-assisted tree decline studies.
Ouagadougou and Bobo-Dioulasso remain the two major urban centers in Burkina Faso with an increasing trend in human footprint. The research aimed at analyzing the Land Use/Land Cover (LULC) dynamics in the two cities between 2003 and 2021 using intensity analysis, which decomposes LULC changes into interval, category and transition levels. The satellite data used for this research were composed of surface reflectance imagery from Landsat 5, Landsat 7 and Landsat 8 acquired from the Google Earth Engine Data Catalogue. The Random Forest, Support Vector Machine and Gradient Tree Boost algorithms were employed to run supervised image classifications for four selected years including 2003, 2009, 2015 and 2021. The results showed that the landscape is changing in both cities due to rapid urbanization. Ouagadougou experienced more rapid changes than Bobo-Dioulasso, with a maximum annual change intensity of 3.61% recorded between 2015 and 2021 against 2.22% in Bobo-Dioulasso for the period 2009–2015. The transition of change was mainly towards built-up areas, which gain targeted bare and agricultural lands in both cities. This situation has led to a 78.12% increase of built-up surfaces in Ouagadougou, while 42.24% of agricultural land area was lost. However, in Bobo-Dioulasso, the built class has increased far more by 140.67%, and the agricultural land areas experienced a gain of 1.38% compared with the 2003 baseline. The study demonstrates that the human footprint is increasing in both cities making the inhabitants vulnerable to environmental threats such as flooding and the effect of an Urban Heat Island, which is information that could serve as guide for sustainable urban land use planning.
Forests are essential for global environmental well-being because of their rich provision of ecosystem services and regulating factors. Global forests are under increasing pressure from climate change, resource extraction, and anthropologically-driven disturbances. The results are dramatic losses of habitats accompanied with the reduction of species diversity. There is the urgent need for forest biodiversity monitoring comprising analysis on α, β, and γ scale to identify hotspots of biodiversity. Remote sensing enables large-scale monitoring at multiple spatial and temporal resolutions. Concepts of remotely sensed spectral diversity have been identified as promising methodologies for the consistent and multi-temporal analysis of forest biodiversity. This review provides a first time focus on the three spectral diversity concepts “vegetation indices”, “spectral information content”, and “spectral species” for forest biodiversity monitoring based on airborne and spaceborne remote sensing. In addition, the reviewed articles are analyzed regarding the spatiotemporal distribution, remote sensing sensors, temporal scales and thematic foci. We identify multispectral sensors as primary data source which underlines the focus on optical diversity as a proxy for forest biodiversity. Moreover, there is a general conceptual focus on the analysis of spectral information content. In recent years, the spectral species concept has raised attention and has been applied to Sentinel-2 and MODIS data for the analysis from local spectral species to global spectral communities. Novel remote sensing processing capacities and the provision of complementary remote sensing data sets offer great potentials for large-scale biodiversity monitoring in the future.
Snow is a vital environmental parameter and dynamically responsive to climate change, particularly in mountainous regions. Snow cover can be monitored at variable spatial scales using Earth Observation (EO) data. Long-lasting remote sensing missions enable the generation of multi-decadal time series and thus the detection of long-term trends. However, there have been few attempts to use these to model future snow cover dynamics. In this study, we, therefore, explore the potential of such time series to forecast the Snow Line Elevation (SLE) in the European Alps. We generate monthly SLE time series from the entire Landsat archive (1985–2021) in 43 Alpine catchments. Positive long-term SLE change rates are detected, with the highest rates (5–8 m/y) in the Western and Central Alps. We utilize this SLE dataset to implement and evaluate seven uni-variate time series modeling and forecasting approaches. The best results were achieved by Random Forests, with a Nash–Sutcliffe efficiency (NSE) of 0.79 and a Mean Absolute Error (MAE) of 258 m, Telescope (0.76, 268 m), and seasonal ARIMA (0.75, 270 m). Since the model performance varies strongly with the input data, we developed a combined forecast based on the best-performing methods in each catchment. This approach was then used to forecast the SLE for the years 2022–2029. In the majority of the catchments, the shift of the forecast median SLE level retained the sign of the long-term trend. In cases where a deviating SLE dynamic is forecast, a discussion based on the unique properties of the catchment and past SLE dynamics is required. In the future, we expect major improvements in our SLE forecasting efforts by including external predictor variables in a multi-variate modeling approach.
Drought is a recurring natural climatic hazard event over terrestrial land; it poses devastating threats to human health, the economy, and the environment. Given the increasing climate crisis, it is likely that extreme drought phenomena will become more frequent, and their impacts will probably be more devastating. Drought observations from space, therefore, play a key role in dissimilating timely and accurate information to support early warning drought management and mitigation planning, particularly in sparse in-situ data regions. In this paper, we reviewed drought-related studies based on Earth observation (EO) products in Southeast Asia between 2000 and 2021. The results of this review indicated that drought publications in the region are on the increase, with a majority (70%) of the studies being undertaken in Vietnam, Thailand, Malaysia and Indonesia. These countries also accounted for nearly 97% of the economic losses due to drought extremes. Vegetation indices from multispectral optical remote sensing sensors remained a primary source of data for drought monitoring in the region. Many studies (~21%) did not provide accuracy assessment on drought mapping products, while precipitation was the main data source for validation. We observed a positive association between spatial extent and spatial resolution, suggesting that nearly 81% of the articles focused on the local and national scales. Although there was an increase in drought research interest in the region, challenges remain regarding large-area and long time-series drought measurements, the combined drought approach, machine learning-based drought prediction, and the integration of multi-sensor remote sensing products (e.g., Landsat and Sentinel-2). Satellite EO data could be a substantial part of the future efforts that are necessary for mitigating drought-related challenges, ensuring food security, establishing a more sustainable economy, and the preservation of the natural environment in the region.
Spatiotemporal Fusion Modelling Using STARFM: Examples of Landsat 8 and Sentinel-2 NDVI in Bavaria
(2022)
The increasing availability and variety of global satellite products provide a new level of data with different spatial, temporal, and spectral resolutions; however, identifying the most suited resolution for a specific application consumes increasingly more time and computation effort. The region’s cloud coverage additionally influences the choice of the best trade-off between spatial and temporal resolution, and different pixel sizes of remote sensing (RS) data may hinder the accurate monitoring of different land cover (LC) classes such as agriculture, forest, grassland, water, urban, and natural-seminatural. To investigate the importance of RS data for these LC classes, the present study fuses NDVIs of two high spatial resolution data (high pair) (Landsat (30 m, 16 days; L) and Sentinel-2 (10 m, 5–6 days; S), with four low spatial resolution data (low pair) (MOD13Q1 (250 m, 16 days), MCD43A4 (500 m, one day), MOD09GQ (250 m, one-day), and MOD09Q1 (250 m, eight day)) using the spatial and temporal adaptive reflectance fusion model (STARFM), which fills regions’ cloud or shadow gaps without losing spatial information. These eight synthetic NDVI STARFM products (2: high pair multiply 4: low pair) offer a spatial resolution of 10 or 30 m and temporal resolution of 1, 8, or 16 days for the entire state of Bavaria (Germany) in 2019. Due to their higher revisit frequency and more cloud and shadow-free scenes (S = 13, L = 9), Sentinel-2 (overall R\(^2\) = 0.71, and RMSE = 0.11) synthetic NDVI products provide more accurate results than Landsat (overall R\(^2\) = 0.61, and RMSE = 0.13). Likewise, for the agriculture class, synthetic products obtained using Sentinel-2 resulted in higher accuracy than Landsat except for L-MOD13Q1 (R\(^2\) = 0.62, RMSE = 0.11), resulting in similar accuracy preciseness as S-MOD13Q1 (R\(^2\) = 0.68, RMSE = 0.13). Similarly, comparing L-MOD13Q1 (R\(^2\) = 0.60, RMSE = 0.05) and S-MOD13Q1 (R\(^2\) = 0.52, RMSE = 0.09) for the forest class, the former resulted in higher accuracy and precision than the latter. Conclusively, both L-MOD13Q1 and S-MOD13Q1 are suitable for agricultural and forest monitoring; however, the spatial resolution of 30 m and low storage capacity makes L-MOD13Q1 more prominent and faster than that of S-MOD13Q1 with the 10-m spatial resolution.
Nationwide and consistent information on agricultural land use forms an important basis for sustainable land management maintaining food security, (agro)biodiversity, and soil fertility, especially as German agriculture has shown high vulnerability to climate change. Sentinel-1 and Sentinel-2 satellite data of the Copernicus program offer time series with temporal, spatial, radiometric, and spectral characteristics that have great potential for mapping and monitoring agricultural crops. This paper presents an approach which synergistically uses these multispectral and Synthetic Aperture Radar (SAR) time series for the classification of 17 crop classes at 10 m spatial resolution for Germany in the year 2018. Input data for the Random Forest (RF) classification are monthly statistics of Sentinel-1 and Sentinel-2 time series. This approach reduces the amount of input data and pre-processing steps while retaining phenological information, which is crucial for crop type discrimination. For training and validation, Land Parcel Identification System (LPIS) data were available covering 15 of the 16 German Federal States. An overall map accuracy of 75.5% was achieved, with class-specific F1-scores above 80% for winter wheat, maize, sugar beet, and rapeseed. By combining optical and SAR data, overall accuracies could be increased by 6% and 9%, respectively, compared to single sensor approaches. While no increase in overall accuracy could be achieved by stratifying the classification in natural landscape regions, the class-wise accuracies for all but the cereal classes could be improved, on average, by 7%. In comparison to census data, the crop areas could be approximated well with, on average, only 1% of deviation in class-specific acreages. Using this streamlined approach, similar accuracies for the most widespread crop types as well as for smaller permanent crop classes were reached as in other Germany-wide crop type studies, indicating its potential for repeated nationwide crop type mapping.
Inland surface water is often the most accessible freshwater source. As opposed to groundwater, surface water is replenished in a comparatively quick cycle, which makes this vital resource — if not overexploited — sustainable. From a global perspective, freshwater is plentiful. Still, depending on the region, surface water availability is severely limited. Additionally, climate change and human interventions act as large-scale drivers and cause dramatic changes in established surface water dynamics. Actions have to be taken to secure sustainable water availability and usage. This requires informed decision making based on reliable environmental data. Monitoring inland surface water dynamics is therefore more important than ever. Remote sensing is able to delineate surface water in a number of ways by using optical as well as active and passive microwave sensors. In this review, we look at the proceedings within this discipline by reviewing 233 scientific works. We provide an extensive overview of used sensors, the spatial and temporal resolution of studies, their thematic foci, and their spatial distribution. We observe that a wide array of available sensors and datasets, along with increasing computing capacities, have shaped the field over the last years. Multiple global analysis-ready products are available for investigating surface water area dynamics, but so far none offer high spatial and temporal resolution.
In most countries, freight is predominantly transported by road cargo trucks. We present a new satellite remote sensing method for detecting moving trucks on roads using Sentinel-2 data. The method exploits a temporal sensing offset of the Sentinel-2 multispectral instrument, causing spatially and spectrally distorted signatures of moving objects. A random forest classifier was trained (overall accuracy: 84%) on visual-near-infrared-spectra of 2500 globally labelled targets. Based on the classification, the target objects were extracted using a developed recursive neighbourhood search. The speed and the heading of the objects were approximated. Detections were validated by employing 350 globally labelled target boxes (mean F\(_1\) score: 0.74). The lowest F\(_1\) score was achieved in Kenya (0.36), the highest in Poland (0.88). Furthermore, validated at 26 traffic count stations in Germany on in sum 390 dates, the truck detections correlate spatio-temporally with station figures (Pearson r-value: 0.82, RMSE: 43.7). Absolute counts were underestimated on 81% of the dates. The detection performance may differ by season and road condition. Hence, the method is only suitable for approximating the relative truck traffic abundance rather than providing accurate absolute counts. However, existing road cargo monitoring methods that rely on traffic count stations or very high resolution remote sensing data have limited global availability. The proposed moving truck detection method could fill this gap, particularly where other information on road cargo traffic are sparse by employing globally and freely available Sentinel-2 data. It is inferior to the accuracy and the temporal detail of station counts, but superior in terms of spatial coverage.
Detection of grassland mowing events for Germany by combining Sentinel-1 and Sentinel-2 time series
(2022)
Grasslands cover one-third of the agricultural area in Germany and play an important economic role by providing fodder for livestock. In addition, they fulfill important ecosystem services, such as carbon storage, water purification, and the provision of habitats. These ecosystem services usually depend on the grassland management. In central Europe, grasslands are grazed and/or mown, whereby the management type and intensity vary in space and time. Spatial information on the mowing timing and frequency on larger scales are usually not available but would be required in order to assess the ecosystem services, species composition, and grassland yields. Time series of high-resolution satellite remote sensing data can be used to analyze the temporal and spatial dynamics of grasslands. Within this study, we aim to overcome the drawbacks identified by previous studies, such as optical data availability and the lack of comprehensive reference data, by testing the time series of various Sentinel-2 (S2) and Sentinal-1 (S1) parameters and combinations of them in order to detect mowing events in Germany in 2019. We developed a threshold-based algorithm by using information from a comprehensive reference dataset of heterogeneously managed grassland parcels in Germany, obtained by RGB cameras. The developed approach using the enhanced vegetation index (EVI) derived from S2 led to a successful mowing event detection in Germany (60.3% of mowing events detected, F1-Score = 0.64). However, events shortly before, during, or shortly after cloud gaps were missed and in regions with lower S2 orbit coverage fewer mowing events were detected. Therefore, S1-based backscatter, InSAR, and PolSAR features were investigated during S2 data gaps. From these, the PolSAR entropy detected mowing events most reliably. For a focus region, we tested an integrated approach by combining S2 and S1 parameters. This approach detected additional mowing events, but also led to many false positive events, resulting in a reduction in the F1-Score (from 0.65 of S2 to 0.61 of S2 + S1 for the focus region). According to our analysis, a majority of grasslands in Germany are only mown zero to two times (around 84%) and are probably additionally used for grazing. A small proportion is mown more often than four times (3%). Regions with a generally higher grassland mowing frequency are located in southern, south-eastern, and northern Germany.
The Northern Bald Ibis (Geronticus eremita, NBI) is an endangered migratory species, which went extinct in Europe in the 17th century. Currently, a translocation project in the frame of the European LIFE program is carried out, to reintroduce a migratory population with breeding colonies in the northern and southern Alpine foothills and a common wintering area in southern Tuscany. The population meanwhile consists of about 200 individuals, with about 90% of them carrying a GPS device on their back. We used biologging data from 2021 to model the habitat suitability for the species in the northern Alpine foothills. To set up a species distribution model, indices describing environmental conditions were calculated from satellite images of Landsat-8, and in addition to the well-proven use of optical remote sensing data, we also included Sentinel-1 actively sensed observation data, as well as climate and urbanization data. A random forest model was fitted on NBI GPS positions, which we used to identify regions with high predicted foraging suitability within the northern Alpine foothills. The model resulted in 84.5% overall accuracy. Elevation and slope had the highest predictive power, followed by grass cover and VV intensity of Sentinel-1 radar data. The map resulting from the model predicts the highest foraging suitability for valley floors, especially of Inn, Rhine, and Salzach-Valley as well as flatlands, like the Swiss Plateau and the agricultural areas surrounding Lake Constance. Areas with a high suitability index largely overlap with known historic breeding sites. This is particularly noteworthy because the model only refers to foraging habitats without considering the availability of suitable breeding cliffs. Detailed analyses identify the transition zone from extensive grassland management to intensive arable farming as the northern range limit. The modeling outcome allows for defining suitable areas for further translocation and management measures in the frame of the European NBI reintroduction program. Although required in the international IUCN translocation guidelines, the use of models in the context of translocation projects is still not common and in the case of the Northern Bald Ibis not considered in the present Single Species Action Plan of the African-Eurasian Migratory Water bird Agreement. Our species distribution model represents a contemporary snapshot, but sustainability is essential for conservation planning, especially in times of climate change. In this regard, a further model could be optimized by investigating sustainable land use, temporal dynamics, and climate change scenarios.
A first assessment of canopy cover loss in Germany's forests after the 2018–2020 drought years
(2022)
Central Europe was hit by several unusually strong periods of drought and heat between 2018 and 2020. These droughts affected forest ecosystems. Cascading effects with bark beetle infestations in spruce stands were fatal to vast forest areas in Germany. We present the first assessment of canopy cover loss in Germany for the period of January 2018–April 2021. Our approach makes use of dense Sentinel-2 and Landsat-8 time-series data. We computed the disturbance index (DI) from the tasseled cap components brightness, greenness, and wetness. Using quantiles, we generated monthly DI composites and calculated anomalies in a reference period (2017). From the resulting map, we calculated the canopy cover loss statistics for administrative entities. Our results show a canopy cover loss of 501,000 ha for Germany, with large regional differences. The losses were largest in central Germany and reached up to two-thirds of coniferous forest loss in some districts. Our map has high spatial (10 m) and temporal (monthly) resolution and can be updated at any time.
The analysis of the Earth system and interactions among its spheres is increasingly important to improve the understanding of global environmental change. In this regard, Earth observation (EO) is a valuable tool for monitoring of long term changes over the land surface and its features. Although investigations commonly study environmental change by means of a single EO-based land surface variable, a joint exploitation of multivariate land surface variables covering several spheres is still rarely performed. In this regard, we present a novel methodological framework for both, the automated processing of multisource time series to generate a unified multivariate feature space, as well as the application of statistical time series analysis techniques to quantify land surface change and driving variables. In particular, we unify multivariate time series over the last two decades including vegetation greenness, surface water area, snow cover area, and climatic, as well as hydrological variables. Furthermore, the statistical time series analyses include quantification of trends, changes in seasonality, and evaluation of drivers using the recently proposed causal discovery algorithm Peter and Clark Momentary Conditional Independence (PCMCI). We demonstrate the functionality of our methodological framework using Indo-Gangetic river basins in South Asia as a case study. The time series analyses reveal increasing trends in vegetation greenness being largely dependent on water availability, decreasing trends in snow cover area being mostly negatively coupled to temperature, and trends of surface water area to be spatially heterogeneous and linked to various driving variables. Overall, the obtained results highlight the value and suitability of this methodological framework with respect to global climate change research, enabling multivariate time series preparation, derivation of detailed information on significant trends and seasonality, as well as detection of causal links with minimal user intervention. This study is the first to use multivariate time series including several EO-based variables to analyze land surface dynamics over the last two decades using the causal discovery algorithm PCMCI.
Intercomparison of satellite-derived vegetation phenology is scarce in remote locations because of the limited coverage area and low temporal resolution of field observations. By their reliable near-ground observations and high-frequency data collection, PhenoCams can be a robust tool for intercomparison of land surface phenology derived from satellites. This study aims to investigate the transition dates of black spruce (Picea mariana (Mill.) B.S.P.) phenology by comparing fortnightly the MODIS normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI) extracted using the Google Earth Engine (GEE) platform with the daily PhenoCam-based green chromatic coordinate (GCC) index. Data were collected from 2016 to 2019 by PhenoCams installed in six mature stands along a latitudinal gradient of the boreal forests of Quebec, Canada. All time series were fitted by double-logistic functions, and the estimated parameters were compared between NDVI, EVI, and GCC. The onset of GCC occurred in the second week of May, whereas the ending of GCC occurred in the last week of September. We demonstrated that GCC was more correlated with EVI (R\(^2\) from 0.66 to 0.85) than NDVI (R\(^2\) from 0.52 to 0.68). In addition, the onset and ending of phenology were shown to differ by 3.5 and 5.4 days between EVI and GCC, respectively. Larger differences were detected between NDVI and GCC, 17.05 and 26.89 days for the onset and ending, respectively. EVI showed better estimations of the phenological dates than NDVI. This better performance is explained by the higher spectral sensitivity of EVI for multiple canopy leaf layers due to the presence of an additional blue band and an optimized soil factor value. Our study demonstrates that the phenological observations derived from PhenoCam are comparable with the EVI index. We conclude that EVI is more suitable than NDVI to assess phenology in evergreen species of the northern boreal region, where PhenoCam data are not available. The EVI index could be used as a reliable proxy of GCC for monitoring evergreen species phenology in areas with reduced access, or where repeated data collection from remote areas are logistically difficult due to the extreme weather.
Freely available satellite data at Google Earth Engine (GEE) cloud platform enables vegetation phenology analysis across different scales very efficiently. We evaluated seasonal and annual phenology of the old-growth Hyrcanian forests (HF) of northern Iran covering an area of ca. 1.9 million ha, and also focused on 15 UNESCO World Heritage Sites. We extracted bi-weekly MODIS-NDVI between 2017 and 2020 in GEE, which was used to identify the range of NDVI between two temporal stages. Then, changes in phenology and growth were analyzed by Sentinel 2-derived Temporal Normalized Phenology Index. We modelled between seasonal phenology and growth by additionally considering elevation, surface temperature, and monthly precipitation. Results indicated considerable difference in onset of forests along the longitudinal gradient of the HF. Faster growth was observed in low- and uplands of the western zone, whereas it was lower in both the mid-elevations and the western outskirts. Longitudinal range was a major driver of vegetation growth, to which environmental factors also differently but significantly contributed (p < 0.0001) along the west-east gradient. Our study developed at GEE provides a benchmark to examine the effects of environmental parameters on the vegetation growth of HF, which cover mountainous areas with partly no or limited accessibility.
Earth Observation satellite data allows for the monitoring of the surface of our planet at predefined intervals covering large areas. However, there is only one medium resolution sensor family in orbit that enables an observation time span of 40 and more years at a daily repeat interval. This is the AVHRR sensor family. If we want to investigate the long-term impacts of climate change on our environment, we can only do so based on data that remains available for several decades. If we then want to investigate processes with respect to climate change, we need very high temporal resolution enabling the generation of long-term time series and the derivation of related statistical parameters such as mean, variability, anomalies, and trends. The challenges to generating a well calibrated and harmonized 40-year-long time series based on AVHRR sensor data flown on 14 different platforms are enormous. However, only extremely thorough pre-processing and harmonization ensures that trends found in the data are real trends and not sensor-related (or other) artefacts. The generation of European-wide time series as a basis for the derivation of a multitude of parameters is therefore an extremely challenging task, the details of which are presented in this paper.
Land Surface Temperature (LST) is an important parameter for tracing the impact of changing climatic conditions on our environment. Describing the interface between long- and shortwave radiation fluxes, as well as between turbulent heat fluxes and the ground heat flux, LST plays a crucial role in the global heat balance. Satellite-derived LST is an indispensable tool for monitoring these changes consistently over large areas and for long time periods. Data from the AVHRR (Advanced Very High-Resolution Radiometer) sensors have been available since the early 1980s. In the TIMELINE project, LST is derived for the entire operating period of AVHRR sensors over Europe at a 1 km spatial resolution. In this study, we present the validation results for the TIMELINE AVHRR daytime LST. The validation approach consists of an assessment of the temporal consistency of the AVHRR LST time series, an inter-comparison between AVHRR LST and in situ LST, and a comparison of the AVHRR LST product with concurrent MODIS (Moderate Resolution Imaging Spectroradiometer) LST. The results indicate the successful derivation of stable LST time series from multi-decadal AVHRR data. The validation results were investigated regarding different LST, TCWV and VA, as well as land cover classes. The comparisons between the TIMELINE LST product and the reference datasets show seasonal and land cover-related patterns. The LST level was found to be the most determinative factor of the error. On average, an absolute deviation of the AVHRR LST by 1.83 K from in situ LST, as well as a difference of 2.34 K from the MODIS product, was observed.
Coal mining, an important human activity, disturbs soil organic carbon (SOC) accumulation and decomposition, eventually affecting terrestrial carbon cycling and the sustainability of human society. However, changes of SOC content and their relation with influential factors in coal mining areas remained unclear. In the study, predictive models of SOC content were developed based on field sampling and Landsat images for different land-use types (grassland, forest, farmland, and bare land) of the largest coal mining area in China (i.e., Shendong). The established models were employed to estimate SOC content across the Shendong mining area during 1990–2020, followed by an investigation into the impacts of climate change and human disturbance on SOC content by a Geo-detector. Results showed that the models produced satisfactory results (R\(^2\) > 0.69, p < 0.05), demonstrating that SOC content over a large coal mining area can be effectively assessed using remote sensing techniques. Results revealed that average SOC content in the study area rose from 5.67 gC·kg\(^{−1}\) in 1990 to 9.23 gC·kg\(^{−1}\) in 2010 and then declined to 5.31 gC·Kg\(^{−1}\) in 2020. This could be attributed to the interaction between the disturbance of soil caused by coal mining and the improvement of eco-environment by land reclamation. Spatially, the SOC content of farmland was the highest, followed by grassland, and that of bare land was the lowest. SOC accumulation was inhibited by coal mining activities, with the effect of high-intensity mining being lower than that of moderate- and low-intensity mining activities. Land use was found to be the strongest individual influencing factor for SOC content changes, while the interaction between vegetation coverage and precipitation exerted the most significant influence on the variability of SOC content. Furthermore, the influence of mining intensity combined with precipitation was 10 times higher than that of mining intensity alone.
Vietnam's 3260 km coastline is densely populated, experiences rapid urban and economic growth, and faces at the same time a high risk of coastal hazards. Satellite archives provide a free and powerful opportunity for long-term area-wide monitoring of the coastal zone. This paper presents an automated analysis of coastline dynamics from 1986 to 2021 for Vietnam's entire coastal zone using the Landsat archive. The proposed method is implemented within the cloud-computing platform Google Earth Engine to only involve publicly and globally available datasets and tools. We generated annual coastline composites representing the mean-high water level and extracted sub-pixel coastlines. We further quantified coastline change rates along shore-perpendicular transects, revealing that half of Vietnam's coast did not experience significant change, while the remaining half is classified as erosional (27.7%) and accretional (27.1%). A hotspot analysis shows that coastal segments with the highest change rates are concentrated in the low-lying deltas of the Mekong River in the south and the Red River in the north. Hotspots with the highest accretion rates of up to +47 m/year are mainly associated with the construction of artificial coastlines, while hotspots with the highest erosion rates of −28 m/year may be related to natural sediment redistribution and human activity.
An increasing amount of Brazilian rainforest is being lost or degraded for various reasons, both anthropogenic and natural, leading to a loss of biodiversity and further global consequences. Especially in the Brazilian state of Mato Grosso, soy production and large-scale cattle farms led to extensive losses of rainforest in recent years. We used a spectral mixture approach followed by a decision tree classification based on more than 30 years of Landsat data to quantify these losses. Research has shown that current methods for assessing forest degradation are lacking accuracy. Therefore, we generated classifications to determine land cover changes for each year, focusing on both cleared and degraded forest land. The analyses showed a decrease in forest area in Mato Grosso by 28.8% between 1986 and 2020. In order to measure changed forest structures for the selected period, fragmentation analyses based on diverse landscape metrics were carried out for the municipality of Colniza in Mato Grosso. It was found that forest areas experienced also a high degree of fragmentation over the study period, with an increase of 83.3% of the number of patches and a decrease of the mean patch area of 86.1% for the selected time period, resulting in altered habitats for flora and fauna.
Mapping of lava flows in unvegetated areas of active volcanoes using optical satellite data is challenging due to spectral similarities of volcanic deposits and the surrounding background. Using very high-resolution PlanetScope data, this study introduces a novel object-oriented classification approach for mapping lava flows in both vegetated and unvegetated areas during several eruptive phases of three Indonesian volcanoes (Karangetang 2018/2019, Agung 2017, Krakatau 2018/2019). For this, change detection analysis based on PlanetScope imagery for mapping loss of vegetation due to volcanic activity (e.g., lava flows) is combined with the analysis of changes in texture and brightness, with hydrological runoff modelling and with analysis of thermal anomalies derived from Sentinel-2 or Landsat-8. Qualitative comparison of the mapped lava flows showed good agreement with multispectral false color time series (Sentinel-2 and Landsat-8). Reports of the Global Volcanism Program support the findings, indicating the developed lava mapping approach produces valuable results for monitoring volcanic hazards. Despite the lack of bands in infrared wavelengths, PlanetScope proves beneficial for the assessment of risk and near-real-time monitoring of active volcanoes due to its high spatial (3 m) and temporal resolution (mapping of all subaerial volcanoes on a daily basis).
Comparing PlanetScope and Sentinel-2 imagery for mapping mountain pines in the Sarntal Alps, Italy
(2022)
The mountain pine (Pinus mugo ssp. Mugo Turra) is an important component of the alpine treeline ecotone and fulfills numerous ecosystem functions. To understand and quantify the impacts of increasing logging activities and climatic changes in the European Alps, accurate information on the occurrence and distribution of mountain pine stands is needed. While Earth observation provides up-to-date information on land cover, space-borne mapping of mountain pines is challenging as different coniferous species are spectrally similar, and small-structured patches may remain undetected due to the sensor’s spatial resolution. This study uses multi-temporal optical imagery from PlanetScope (3 m) and Sentinel-2 (10 m) and combines them with additional features (e.g., textural statistics (homogeneity, contrast, entropy, spatial mean and spatial variance) from gray level co-occurrence matrix (GLCM), topographic features (elevation, slope and aspect) and canopy height information) to overcome the present challenges in mapping mountain pine stands. Specifically, we assessed the influence of spatial resolution and feature space composition including the GLCM window size for textural features. The study site is covering the Sarntal Alps, Italy, a region known for large stands of mountain pine. Our results show that mountain pines can be accurately mapped (PlanetScope (90.96%) and Sentinel-2 (90.65%)) by combining all features. In general, Sentinel-2 can achieve comparable results to PlanetScope independent of the feature set composition, despite the lower spatial resolution. In particular, the inclusion of textural features improved the accuracy by +8% (PlanetScope) and +3% (Sentinel-2), whereas accuracy improvements of topographic features and canopy height were low. The derived map of mountain pines in the Sarntal Alps supports local forest management to monitor and assess recent and ongoing anthropogenic and climatic changes at the treeline. Furthermore, our study highlights the importance of freely available Sentinel-2 data and image-derived textural features to accurately map mountain pines in Alpine environments.
Arctic permafrost coasts become increasingly vulnerable due to environmental drivers such as the reduced sea-ice extent and duration as well as the thawing of permafrost itself. A continuous quantification of the erosion process on large to circum-Arctic scales is required to fully assess the extent and understand the consequences of eroding permafrost coastlines. This study presents a novel approach to quantify annual Arctic coastal erosion and build-up rates based on Sentinel-1 (S1) Synthetic Aperture RADAR (SAR) backscatter data, in combination with Deep Learning (DL) and Change Vector Analysis (CVA). The methodology includes the generation of a high-quality Arctic coastline product via DL, which acted as a reference for quantifying coastal erosion and build-up rates from annual median and standard deviation (sd) backscatter images via CVA. The analysis was applied on ten test sites distributed across the Arctic and covering about 1038 km of coastline. Results revealed maximum erosion rates of up to 160 m for some areas and an average erosion rate of 4.37 m across all test sites within a three-year temporal window from 2017 to 2020. The observed erosion rates within the framework of this study agree with findings published in the previous literature. The proposed methods and data can be applied on large scales and, prospectively, even for the entire Arctic. The generated products may be used for quantifying the loss of frozen ground, estimating the release of stored organic material, and can act as a basis for further related studies in Arctic coastal environments.
The monitoring of species and functional diversity is of increasing relevance for the development of strategies for the conservation and management of biodiversity. Therefore, reliable estimates of the performance of monitoring techniques across taxa become important. Using a unique dataset, this study investigates the potential of airborne LiDAR-derived variables characterizing vegetation structure as predictors for animal species richness at the southern slopes of Mount Kilimanjaro. To disentangle the structural LiDAR information from co-factors related to elevational vegetation zones, LiDAR-based models were compared to the predictive power of elevation models. 17 taxa and 4 feeding guilds were modeled and the standardized study design allowed for a comparison across the assemblages. Results show that most taxa (14) and feeding guilds (3) can be predicted best by elevation with normalized RMSE values but only for three of those taxa and two of those feeding guilds the difference to other models is significant. Generally, modeling performances between different models vary only slightly for each assemblage. For the remaining, structural information at most showed little additional contribution to the performance. In summary, LiDAR observations can be used for animal species prediction. However, the effort and cost of aerial surveys are not always in proportion with the prediction quality, especially when the species distribution follows zonal patterns, and elevation information yields similar results.
The fast and accurate yield estimates with the increasing availability and variety of global satellite products and the rapid development of new algorithms remain a goal for precision agriculture and food security. However, the consistency and reliability of suitable methodologies that provide accurate crop yield outcomes still need to be explored. The study investigates the coupling of crop modeling and machine learning (ML) to improve the yield prediction of winter wheat (WW) and oil seed rape (OSR) and provides examples for the Free State of Bavaria (70,550 km2), Germany, in 2019. The main objectives are to find whether a coupling approach [Light Use Efficiency (LUE) + Random Forest (RF)] would result in better and more accurate yield predictions compared to results provided with other models not using the LUE. Four different RF models [RF1 (input: Normalized Difference Vegetation Index (NDVI)), RF2 (input: climate variables), RF3 (input: NDVI + climate variables), RF4 (input: LUE generated biomass + climate variables)], and one semi-empiric LUE model were designed with different input requirements to find the best predictors of crop monitoring. The results indicate that the individual use of the NDVI (in RF1) and the climate variables (in RF2) could not be the most accurate, reliable, and precise solution for crop monitoring; however, their combined use (in RF3) resulted in higher accuracies. Notably, the study suggested the coupling of the LUE model variables to the RF4 model can reduce the relative root mean square error (RRMSE) from −8% (WW) and −1.6% (OSR) and increase the R
2 by 14.3% (for both WW and OSR), compared to results just relying on LUE. Moreover, the research compares models yield outputs by inputting three different spatial inputs: Sentinel-2(S)-MOD13Q1 (10 m), Landsat (L)-MOD13Q1 (30 m), and MOD13Q1 (MODIS) (250 m). The S-MOD13Q1 data has relatively improved the performance of models with higher mean R
2 [0.80 (WW), 0.69 (OSR)], and lower RRMSE (%) (9.18, 10.21) compared to L-MOD13Q1 (30 m) and MOD13Q1 (250 m). Satellite-based crop biomass, solar radiation, and temperature are found to be the most influential variables in the yield prediction of both crops.
Enhancing digital and precision agriculture is currently inevitable to overcome the economic and environmental challenges of the agriculture in the 21st century. The purpose of this study was to generate and compare management zones (MZ) based on the Sentinel-2 satellite data for variable rate application of mineral nitrogen in wheat production, calculated using different remote sensing (RS)-based models under varied soil, yield and crop data availability. Three models were applied, including (1) a modified “RS- and threshold-based clustering”, (2) a “hybrid-based, unsupervised clustering”, in which data from different sources were combined for MZ delineation, and (3) a “RS-based, unsupervised clustering”. Various data processing methods including machine learning were used in the model development. Statistical tests such as the Paired Sample T-test, Kruskal–Wallis H-test and Wilcoxon signed-rank test were applied to evaluate the final delineated MZ maps. Additionally, a procedure for improving models based on information about phenological phases and the occurrence of agricultural drought was implemented. The results showed that information on agronomy and climate enables improving and optimizing MZ delineation. The integration of prior knowledge on new climate conditions (drought) in image selection was tested for effective use of the models. Lack of this information led to the infeasibility of obtaining optimal results. Models that solely rely on remote sensing information are comparatively less expensive than hybrid models. Additionally, remote sensing-based models enable delineating MZ for fertilizer recommendations that are temporally closer to fertilization times.
Mapping aquaculture ponds for the coastal zone of Asia with Sentinel-1 and Sentinel-2 time series
(2021)
Asia dominates the world's aquaculture sector, generating almost 90 percent of its total annual global production. Fish, shrimp, and mollusks are mainly farmed in land-based pond aquaculture systems and serve as a primary protein source for millions of people. The total production and area occupied for pond aquaculture has expanded rapidly in coastal regions in Asia since the early 1990s. The growth of aquaculture was mainly boosted by an increasing demand for fish and seafood from a growing world population. The aquaculture sector generates income and employment, contributes to food security, and has become a billion-dollar industry with high socio-economic value, but has also led to severe environmental degradation. In this regard, geospatial information on aquaculture can support the management of this growing food sector for the sustainable development of coastal ecosystems, resources, and human health. With free and open access to the rapidly growing volume of data from the Copernicus Sentinel missions as well as machine learning algorithms and cloud computing services, we extracted coastal aquaculture at a continental scale. We present a multi-sensor approach that utilizes Earth observation time series data for the mapping of pond aquaculture within the entire Asian coastal zone, defined as the onshore area up to 200 km from the coastline. In this research, we developed an object-based framework to detect and extract aquaculture at a single-pond level based on temporal features derived from high-spatial-resolution SAR and optical satellite data acquired from the Sentinel-1 and Sentinel-2 satellites. In a second step, we performed spatial and statistical data analyses of the Earth-observation-derived aquaculture dataset to investigate spatial distribution and identify production hotspots at various administrative units at regional, national, and sub-national scale.
Semi-arid tree covers, in both high and coppice growth forms, play an essential role in protecting water and soil resources and provides multiple ecosystem services across fragile ecosystems. Thus, they require continuous inventories. Quantification of forest structure in these tree covers provides important measures for their management and biodiversity conservation. We present a framework, based on consumer-grade UAV photogrammetry, to separately estimate primary variables of tree height (H) and crown area (A) across diverse coppice and high stands dominated by Quercus brantii Lindl. along the latitudinal gradient of Zagros mountains of western Iran. Then, multivariate linear regressions were parametrized with H and A to estimate the diameter at breast height (DBH) of high trees because of its importance to accelerate the existing practical DBH inventories across Zagros Forests. The estimated variables were finally applied to a model tree aboveground biomass (AGB) for both vegetative growth forms by local allometric equations and Random Forest models. In each step, the estimated variables were evaluated against the field reference values, indicating practically high accuracies reaching root mean square error (RMSE) of 0.68 m and 4.74 cm for H and DBH, as well as relative RMSE < 10% for AGB estimates. The results generally suggest an effective framework for single tree-based attribute estimation over mountainous, semi-arid coppice, and high stands.
Wetlands are one of the most important ecosystems due to their critical services to both humans and the environment. Therefore, wetland mapping and monitoring are essential for their conservation. In this regard, remote sensing offers efficient solutions due to the availability of cost-efficient archived images over different spatial scales. However, a lack of sufficient consistent training samples at different times is a significant limitation of multi-temporal wetland monitoring. In this study, a new training sample migration method was developed to identify unchanged training samples to be used in wetland classification and change analyses over the International Shadegan Wetland (ISW) areas of southwestern Iran. To this end, we first produced the wetland map of a reference year (2020), for which we had training samples, by combining Sentinel-1 and Sentinel-2 images and the Random Forest (RF) classifier in Google Earth Engine (GEE). The Overall Accuracy (OA) and Kappa coefficient (KC) of this reference map were 97.93% and 0.97, respectively. Then, an automatic change detection method was developed to migrate unchanged training samples from the reference year to the target years of 2018, 2019, and 2021. Within the proposed method, three indices of the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), and the mean Standard Deviation (SD) of the spectral bands, along with two similarity measures of the Euclidean Distance (ED) and Spectral Angle Distance (SAD), were computed for each pair of reference–target years. The optimum threshold for unchanged samples was also derived using a histogram thresholding approach, which led to selecting the samples that were most likely unchanged based on the highest OA and KC for classifying the test dataset. The proposed migration sample method resulted in high OAs of 95.89%, 96.83%, and 97.06% and KCs of 0.95, 0.96, and 0.96 for the target years of 2018, 2019, and 2021, respectively. Finally, the migrated samples were used to generate the wetland map for the target years. Overall, our proposed method showed high potential for wetland mapping and monitoring when no training samples existed for a target year.
Earth observation time series are well suited to monitor global surface dynamics. However, data products that are aimed at assessing large-area dynamics with a high temporal resolution often face various error sources (e.g., retrieval errors, sampling errors) in their acquisition chain. Addressing uncertainties in a spatiotemporal consistent manner is challenging, as extensive high-quality validation data is typically scarce. Here we propose a new method that utilizes time series inherent information to assess the temporal interpolation uncertainty of time series datasets. For this, we utilized data from the DLR-DFD Global WaterPack (GWP), which provides daily information on global inland surface water. As the time series is primarily based on optical MODIS (Moderate Resolution Imaging Spectroradiometer) images, the requirement of data gap interpolation due to clouds constitutes the main uncertainty source of the product. With a focus on different temporal and spatial characteristics of surface water dynamics, seven auxiliary layers were derived. Each layer provides probability and reliability estimates regarding water observations at pixel-level. This enables the quantification of uncertainty corresponding to the full spatiotemporal range of the product. Furthermore, the ability of temporal layers to approximate unknown pixel states was evaluated for stratified artificial gaps, which were introduced into the original time series of four climatologic diverse test regions. Results show that uncertainty is quantified accurately (>90%), consequently enhancing the product's quality with respect to its use for modeling and the geoscientific community.
Fresh water is a vital natural resource. Earth observation time-series are well suited to monitor corresponding surface dynamics. The DLR-DFD Global WaterPack (GWP) provides daily information on globally distributed inland surface water based on MODIS (Moderate Resolution Imaging Spectroradiometer) images at 250 m spatial resolution. Operating on this spatiotemporal level comes with the drawback of moderate spatial resolution; only coarse pixel-based surface water quantification is possible. To enhance the quantitative capabilities of this dataset, we systematically access subpixel information on fractional water coverage. For this, a linear mixture model is employed, using classification probability and pure pixel reference information. Classification probability is derived from relative datapoint (pixel) locations in feature space. Pure water and non-water reference pixels are located by combining spatial and temporal information inherent to the time-series. Subsequently, the model is evaluated for different input sets to determine the optimal configuration for global processing and pixel coverage types. The performance of resulting water fraction estimates is evaluated on the pixel level in 32 regions of interest across the globe, by comparison to higher resolution reference data (Sentinel-2, Landsat 8). Results show that water fraction information is able to improve the product's performance regarding mixed water/non-water pixels by an average of 11.6% (RMSE). With a Nash-Sutcliffe efficiency of 0.61, the model shows good overall performance. The approach enables the systematic provision of water fraction estimates on a global and daily scale, using only the reflectance and temporal information contained in the input time-series.
A disease is non-communicable when it is not transferred from one person to another. Typical examples include all types of cancer, diabetes, stroke, or allergies, as well as mental diseases. Non-communicable diseases have at least two things in common — environmental impact and chronicity. These diseases are often associated with reduced quality of life, a higher rate of premature deaths, and negative impacts on a countries' economy due to healthcare costs and missing work force. Additionally, they affect the individual's immune system, which increases susceptibility toward communicable diseases, such as the flu or other viral and bacterial infections. Thus, mitigating the effects of non-communicable diseases is one of the most pressing issues of modern medicine, healthcare, and governments in general. Apart from the predisposition toward such diseases (the genome), their occurrence is associated with environmental parameters that people are exposed to (the exposome). Exposure to stressors such as bad air or water quality, noise, extreme heat, or an overall unnatural surrounding all impact the susceptibility to non-communicable diseases. In the identification of such environmental parameters, geoinformation products derived from Earth Observation data acquired by satellites play an increasingly important role. In this paper, we present a review on the joint use of Earth Observation data and public health data for research on non-communicable diseases. We analyzed 146 articles from peer-reviewed journals (Impact Factor ≥ 2) from all over the world that included Earth Observation data and public health data for their assessments. Our results show that this field of synergistic geohealth analyses is still relatively young, with most studies published within the last five years and within national boundaries. While the contribution of Earth Observation, and especially remote sensing-derived geoinformation products on land surface dynamics is on the rise, there is still a huge potential for transdisciplinary integration into studies. We see the necessity for future research and advocate for the increased incorporation of thematically profound remote sensing products with high spatial and temporal resolution into the mapping of exposomes and thus the vulnerability and resilience assessment of a population regarding non-communicable diseases.