526 Mathematische Geografie
Refine
Has Fulltext
- yes (119)
Is part of the Bibliography
- yes (119) (remove)
Year of publication
Document Type
- Journal article (117)
- Doctoral Thesis (2)
Language
- English (119) (remove)
Keywords
- remote sensing (26)
- Sentinel-2 (14)
- time series (10)
- Google Earth Engine (9)
- MODIS (9)
- Sentinel-1 (8)
- earth observation (8)
- machine learning (8)
- Earth observation (6)
- Landsat (6)
Institute
EU-Project number / Contract (GA) number
- 227159 (1)
- 818182 (1)
- LIFE12 BIO/AT/000143 (1)
- LIFE20 NAT/AT/000049 (1)
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.
Contemporary climate change leads to earlier spring phenological events in Europe. In forests, in which overstory strongly regulates the microclimate beneath, it is not clear if further change equally shifts the timing of leaf unfolding for the over- and understory of main deciduous forest species, such as Fagus sylvatica L. (European beech). Furthermore, it is not known yet how this vertical phenological (mis)match — the phenological difference between overstory and understory — affects the remotely sensed satellite signal. To investigate this, we disentangled the start of season (SOS) of overstory F.sylvatica foliage from understory F. sylvatica foliage in forests, within nine quadrants of 5.8 × 5.8 km, stratified over a temperature gradient of 2.5 °C in Bavaria, southeast Germany, in the spring seasons of 2019 and 2020 using time lapse cameras and visual ground observations. We explained SOS dates and vertical phenological (mis)match by canopy temperature and compared these to Sentinel-2 derived SOS in response to canopy temperature. We found that overstory SOS advanced with higher mean April canopy temperature (visual ground observations: −2.86 days per °C; cameras: −2.57 days per °C). However, understory SOS was not significantly affected by canopy temperature. This led to an increase of vertical phenological mismatch with increased canopy temperature (visual ground observations: +3.90 days per °C; cameras: +2.52 days per °C). These results matched Sentinel-2-derived SOS responses, as pixels of higher canopy height advanced more by increased canopy temperature than pixels of lower canopy height. The results may indicate that, with further climate change, spring phenology of F. sylvatica overstory will advance more than F. sylvatica understory, leading to increased vertical phenological mismatch in temperate deciduous forests. This may have major ecological effects, but also methodological consequences for the field of remote sensing, as what the signal senses highly depends on the pixel mean canopy height and the vertical (mis)match.
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.
Recently, locust outbreaks around the world have destroyed agricultural and natural vegetation and caused massive damage endangering food security. Unusual heavy rainfalls in habitats of the desert locust (Schistocerca gregaria) and lack of monitoring due to political conflicts or inaccessibility of those habitats lead to massive desert locust outbreaks and swarms migrating over the Arabian Peninsula, East Africa, India and Pakistan. At the same time, swarms of the Moroccan locust (Dociostaurus maroccanus) in some Central Asian countries and swarms of the Italian locust (Calliptamus italicus) in Russia and China destroyed crops despite developed and ongoing monitoring and control measurements. These recent events underline that the risk and damage caused by locust pests is as present as ever and affects 100 million of human lives despite technical progress in locust monitoring, prediction and control approaches. Remote sensing has become one of the most important data sources in locust management. Since the 1980s, remote sensing data and applications have accompanied many locust management activities and contributed to an improved and more effective control of locust outbreaks and plagues. Recently, open-access remote sensing data archives as well as progress in cloud computing provide unprecedented opportunity for remote sensing-based locust management and research. Additionally, unmanned aerial vehicle (UAV) systems bring up new prospects for a more effective and faster locust control. Nevertheless, the full capacity of available remote sensing applications and possibilities have not been exploited yet. This review paper provides a comprehensive and quantitative overview of international research articles focusing on remote sensing application for locust management and research. We reviewed 110 articles published over the last four decades, and categorized them into different aspects and main research topics to summarize achievements and gaps for further research and application development. The results reveal a strong focus on three species — the desert locust, the migratory locust (Locusta migratoria), and the Australian plague locust (Chortoicetes terminifera) — and corresponding regions of interest. There is still a lack of international studies for other pest species such as the Italian locust, the Moroccan locust, the Central American locust (Schistocerca piceifrons), the South American locust (Schistocerca cancellata), the brown locust (Locustana pardalina) and the red locust (Nomadacris septemfasciata). In terms of applied sensors, most studies utilized Advanced Very-High-Resolution Radiometer (AVHRR), Satellite Pour l’Observation de la Terre VEGETATION (SPOT-VGT), Moderate-Resolution Imaging Spectroradiometer (MODIS) as well as Landsat data focusing mainly on vegetation monitoring or land cover mapping. Application of geomorphological metrics as well as radar-based soil moisture data is comparably rare despite previous acknowledgement of their importance for locust outbreaks. Despite great advance and usage of available remote sensing resources, we identify several gaps and potential for future research to further improve the understanding and capacities of the use of remote sensing in supporting locust outbreak- research and management.
By 2050, two-third of the world’s population will live in cities. In this study, we develop a framework for analyzing urban growth-related imperviousness in North Rhine-Westphalia (NRW) from the 1980s to date using Landsat data. For the baseline 2017-time step, official geodata was extracted to generate labelled data for ten classes, including three classes representing low, middle, and high level of imperviousness. We used the output of the 2017 classification and information based on radiometric bi-temporal change detection for retrospective classification. Besides spectral bands, we calculated several indices and various temporal composites, which were used as an input for Random Forest classification. The results provide information on three imperviousness classes with accuracies exceeding 75%. According to our results, the imperviousness areas grew continuously from 1985 to 2017, with a high imperviousness area growth of more than 167,000 ha, comprising around 30% increase. The information on the expansion of urban areas was integrated with population dynamics data to estimate the progress towards SDG 11. With the intensity analysis and the integration of population data, the spatial heterogeneity of urban expansion and population growth was analysed, showing that the urban expansion rates considerably excelled population growth rates in some regions in NRW. The study highlights the applicability of earth observation data for accurately quantifying spatio-temporal urban dynamics for sustainable urbanization and targeted planning.
Public safety and socio-economic development of the Jharia coalfield (JCF) in India is critically dependent on precise monitoring and comprehensive understanding of coal fires, which have been burning underneath for more than a century. This study utilizes New-Small BAseline Subset (N-SBAS) technique to compute surface deformation time series for 2017–2020 to characterize the spatiotemporal dynamics of coal fires in JCF. The line-of-sight (LOS) surface deformation estimated from ascending and descending Sentinel-1 SAR data are subsequently decomposed to derive precise vertical subsidence estimates. The most prominent subsidence (~22 cm) is observed in Kusunda colliery. The subsidence regions also correspond well with the Landsat-8 based thermal anomaly map and field evidence. Subsequently, the vertical surface deformation time-series is analyzed to characterize temporal variations within the 9.5 km\(^2\) area of coal fires. Results reveal that nearly 10% of the coal fire area is newly formed, while 73% persisted throughout the study period. Vulnerability analyses performed in terms of the susceptibility of the population to land surface collapse demonstrate that Tisra, Chhatatanr, and Sijua are the most vulnerable towns. Our results provide critical information for developing early warning systems and remediation strategies.
Supraglacial meltwater accumulation on ice sheets can be a main driver for accelerated ice discharge, mass loss, and global sea-level-rise. With further increasing surface air temperatures, meltwater-induced hydrofracturing, basal sliding, or surface thinning will cumulate and most likely trigger unprecedented ice mass loss on the Greenland and Antarctic ice sheets. While the Greenland surface hydrological network as well as its impacts on ice dynamics and mass balance has been studied in much detail, Antarctic supraglacial lakes remain understudied with a circum-Antarctic record of their spatio-temporal development entirely lacking. This study provides the first automated supraglacial lake extent mapping method using Sentinel-1 synthetic aperture radar (SAR) imagery over Antarctica and complements the developed optical Sentinel-2 supraglacial lake detection algorithm presented in our companion paper. In detail, we propose the use of a modified U-Net for semantic segmentation of supraglacial lakes in single-polarized Sentinel-1 imagery. The convolutional neural network (CNN) is implemented with residual connections for optimized performance as well as an Atrous Spatial Pyramid Pooling (ASPP) module for multiscale feature extraction. The algorithm is trained on 21,200 Sentinel-1 image patches and evaluated in ten spatially or temporally independent test acquisitions. In addition, George VI Ice Shelf is analyzed for intra-annual lake dynamics throughout austral summer 2019/2020 and a decision-level fused Sentinel-1 and Sentinel-2 maximum lake extent mapping product is presented for January 2020 revealing a more complete supraglacial lake coverage (~770 km\(^2\)) than the individual single-sensor products. Classification results confirm the reliability of the proposed workflow with an average Kappa coefficient of 0.925 and a F\(_1\)-score of 93.0% for the supraglacial water class across all test regions. Furthermore, the algorithm is applied in an additional test region covering supraglacial lakes on the Greenland ice sheet which further highlights the potential for spatio-temporal transferability. Future work involves the integration of more training data as well as intra-annual analyses of supraglacial lake occurrence across the whole continent and with focus on supraglacial lake development throughout a summer melt season and into Antarctic winter.
The monitoring of land cover and land use change is critical for assessing the provision of ecosystem services. One of the sources for long-term land cover change quantification is through the classification of historical and/or current maps. Little research has been done on historical maps using Object-Based Image Analysis (OBIA). This study applied an object-based classification using eCognition tool for analyzing the land cover based on historical maps in the Main river catchment, Upper Franconia, Germany. This allowed land use change analysis between the 1850s and 2015, a time span which covers the phase of industrialization of landscapes in central Europe. The results show a strong increase in urban area by 2600%, a severe loss of cropland (−24%), a moderate reduction in meadows (−4%), and a small gain in forests (+4%). The method proved useful for the application on historical maps due to the ability of the software to create semantic objects. The confusion matrix shows an overall accuracy of 82% for the automatic classification compared to manual reclassification considering all 17 sample tiles. The minimum overall accuracy was 65% for historical maps of poor quality and the maximum was 91% for very high-quality ones. Although accuracy is between high and moderate, coarse land cover patterns in the past and trends in land cover change can be analyzed. We conclude that such long-term analysis of land cover is a prerequisite for quantifying long-term changes in ecosystem services.
Forests in Germany cover around 11.4 million hectares and, thus, a share of 32% of Germany's surface area. Therefore, forests shape the character of the country's cultural landscape. Germany's forests fulfil a variety of functions for nature and society, and also play an important role in the context of climate levelling. Climate change, manifested via rising temperatures and current weather extremes, has a negative impact on the health and development of forests. Within the last five years, severe storms, extreme drought, and heat waves, and the subsequent mass reproduction of bark beetles have all seriously affected Germany’s forests. Facing the current dramatic extent of forest damage and the emerging long-term consequences, the effort to preserve forests in Germany, along with their diversity and productivity, is an indispensable task for the government. Several German ministries have and plan to initiate measures supporting forest health. Quantitative data is one means for sound decision-making to ensure the monitoring of the forest and to improve the monitoring of forest damage. In addition to existing forest monitoring systems, such as the federal forest inventory, the national crown condition survey, and the national forest soil inventory, systematic surveys of forest condition and vulnerability at the national scale can be expanded with the help of a satellite-based earth observation. In this review, we analysed and categorized all research studies published in the last 20 years that focus on the remote sensing of forests in Germany. For this study, 166 citation indexed research publications have been thoroughly analysed with respect to publication frequency, location of studies undertaken, spatial and temporal scale, coverage of the studies, satellite sensors employed, thematic foci of the studies, and overall outcomes, allowing us to identify major research and geoinformation product gaps.
Forecasting spatio-temporal dynamics on the land surface using Earth Observation data — a review
(2020)
Reliable forecasts on the impacts of global change on the land surface are vital to inform the actions of policy and decision makers to mitigate consequences and secure livelihoods. Geospatial Earth Observation (EO) data from remote sensing satellites has been collected continuously for 40 years and has the potential to facilitate the spatio-temporal forecasting of land surface dynamics. In this review we compiled 143 papers on EO-based forecasting of all aspects of the land surface published in 16 high-ranking remote sensing journals within the past decade. We analyzed the literature regarding research focus, the spatial scope of the study, the forecasting method applied, as well as the temporal and technical properties of the input data. We categorized the identified forecasting methods according to their temporal forecasting mechanism and the type of input data. Time-lagged regressions which are predominantly used for crop yield forecasting and approaches based on Markov Chains for future land use and land cover simulation are the most established methods. The use of external climate projections allows the forecasting of numerical land surface parameters up to one hundred years into the future, while auto-regressive time series modeling can account for intra-annual variances. Machine learning methods have been increasingly used in all categories and multivariate modeling that integrates multiple data sources appears to be more popular than univariate auto-regressive modeling despite the availability of continuously expanding time series data. Regardless of the method, reliable EO-based forecasting requires high-level remote sensing data products and the resulting computational demand appears to be the main reason that most forecasts are conducted only on a local scale. In the upcoming years, however, we expect this to change with further advances in the field of machine learning, the publication of new global datasets, and the further establishment of cloud computing for data processing.
Global Navigation Satellite System (GNSS) provides accurate positioning data for vehicular navigation in open outdoor environment. In an indoor environment, Light Detection and Ranging (LIDAR) Simultaneous Localization and Mapping (SLAM) establishes a two-dimensional map and provides positioning data. However, LIDAR can only provide relative positioning data and it cannot directly provide the latitude and longitude of the current position. As a consequence, GNSS/Inertial Navigation System (INS) integrated navigation could be employed in outdoors, while the indoors part makes use of INS/LIDAR integrated navigation and the corresponding switching navigation will make the indoor and outdoor positioning consistent. In addition, when the vehicle enters the garage, the GNSS signal will be blurred for a while and then disappeared. Ambiguous GNSS satellite signals will lead to the continuous distortion or overall drift of the positioning trajectory in the indoor condition. Therefore, an INS/LIDAR seamless integrated navigation algorithm and a switching algorithm based on vehicle navigation system are designed. According to the experimental data, the positioning accuracy of the INS/LIDAR navigation algorithm in the simulated environmental experiment is 50% higher than that of the Dead Reckoning (DR) algorithm. Besides, the switching algorithm developed based on the INS/LIDAR integrated navigation algorithm can achieve 80% success rate in navigation mode switching.
In Earth observation (EO), large-scale land-surface dynamics are traditionally analyzed by investigating aggregated classes. The increase in data with a very high spatial resolution enables investigations on a fine-grained feature level which can help us to better understand the dynamics of land surfaces by taking object dynamics into account. To extract fine-grained features and objects, the most popular deep-learning model for image analysis is commonly used: the convolutional neural network (CNN). In this review, we provide a comprehensive overview of the impact of deep learning on EO applications by reviewing 429 studies on image segmentation and object detection with CNNs. We extensively examine the spatial distribution of study sites, employed sensors, used datasets and CNN architectures, and give a thorough overview of applications in EO which used CNNs. Our main finding is that CNNs are in an advanced transition phase from computer vision to EO. Upon this, we argue that in the near future, investigations which analyze object dynamics with CNNs will have a significant impact on EO research. With a focus on EO applications in this Part II, we complete the methodological review provided in Part I.
Protection and recovery of natural resource and biodiversity requires accurate monitoring at multiple scales. Airborne Laser Scanning (ALS) provides high-resolution imagery that is valuable for monitoring structural changes to vegetation, providing a reliable reference for ecological analyses and comparison purposes, especially if used in conjunction with other remote-sensing and field products. However, the potential of ALS data has not been fully exploited, due to limits in data availability and validation. To bridge this gap, the global network for airborne laser scanner data (GlobALS) has been established as a worldwide network of ALS data providers that aims at linking those interested in research and applications related to natural resources and biodiversity monitoring. The network does not collect data itself but collects metadata and facilitates networking and collaborative research amongst the end-users and data providers. This letter describes this facility, with the aim of broadening participation in GlobALS.
In China, freshwater is an increasingly scarce resource and wetlands are under great pressure. This study focuses on China's second largest freshwater lake in the middle reaches of the Yangtze River — the Dongting Lake — and its surrounding wetlands, which are declared a protected Ramsar site. The Dongting Lake area is also a research region of focus within the Sino-European Dragon Programme, aiming for the international collaboration of Earth Observation researchers. ESA's Copernicus Programme enables comprehensive monitoring with area-wide coverage, which is especially advantageous for large wetlands that are difficult to access during floods. The first year completely covered by Sentinel-1 SAR satellite data was 2016, which is used here to focus on Dongting Lake's wetland dynamics. The well-established, threshold-based approach and the high spatio-temporal resolution of Sentinel-1 imagery enabled the generation of monthly surface water maps and the analysis of the inundation frequency at a 10 m resolution. The maximum extent of the Dongting Lake derived from Sentinel-1 occurred in July 2016, at 2465 km\(^2\), indicating an extreme flood year. The minimum size of the lake was detected in October, at 1331 km\(^2\). Time series analysis reveals detailed inundation patterns and small-scale structures within the lake that were not known from previous studies. Sentinel-1 also proves to be capable of mapping the wetland management practices for Dongting Lake polders and dykes. For validation, the lake extent and inundation duration derived from the Sentinel-1 data were compared with excerpts from the Global WaterPack (frequently derived by the German Aerospace Center, DLR), high-resolution optical data, and in situ water level data, which showed very good agreement for the period studied. The mean monthly extent of the lake in 2016 from Sentinel-1 was 1798 km\(^2\), which is consistent with the Global WaterPack, deviating by only 4%. In summary, the presented analysis of the complete annual time series of the Sentinel-1 data provides information on the monthly behavior of water expansion, which is of interest and relevance to local authorities involved in water resource management tasks in the region, as well as to wetland conservationists concerned with the Ramsar site wetlands of Dongting Lake and to local researchers.
Deep learning (DL) has great influence on large parts of science and increasingly established itself as an adaptive method for new challenges in the field of Earth observation (EO). Nevertheless, the entry barriers for EO researchers are high due to the dense and rapidly developing field mainly driven by advances in computer vision (CV). To lower the barriers for researchers in EO, this review gives an overview of the evolution of DL with a focus on image segmentation and object detection in convolutional neural networks (CNN). The survey starts in 2012, when a CNN set new standards in image recognition, and lasts until late 2019. Thereby, we highlight the connections between the most important CNN architectures and cornerstones coming from CV in order to alleviate the evaluation of modern DL models. Furthermore, we briefly outline the evolution of the most popular DL frameworks and provide a summary of datasets in EO. By discussing well performing DL architectures on these datasets as well as reflecting on advances made in CV and their impact on future research in EO, we narrow the gap between the reviewed, theoretical concepts from CV and practical application in EO.
Supraglacial lakes can have considerable impact on ice sheet mass balance and global sea-level-rise through ice shelf fracturing and subsequent glacier speedup. In Antarctica, the distribution and temporal development of supraglacial lakes as well as their potential contribution to increased ice mass loss remains largely unknown, requiring a detailed mapping of the Antarctic surface hydrological network. In this study, we employ a Machine Learning algorithm trained on Sentinel-2 and auxiliary TanDEM-X topographic data for automated mapping of Antarctic supraglacial lakes. To ensure the spatio-temporal transferability of our method, a Random Forest was trained on 14 training regions and applied over eight spatially independent test regions distributed across the whole Antarctic continent. In addition, we employed our workflow for large-scale application over Amery Ice Shelf where we calculated interannual supraglacial lake dynamics between 2017 and 2020 at full ice shelf coverage. To validate our supraglacial lake detection algorithm, we randomly created point samples over our classification results and compared them to Sentinel-2 imagery. The point comparisons were evaluated using a confusion matrix for calculation of selected accuracy metrics. Our analysis revealed wide-spread supraglacial lake occurrence in all three Antarctic regions. For the first time, we identified supraglacial meltwater features on Abbott, Hull and Cosgrove Ice Shelves in West Antarctica as well as for the entire Amery Ice Shelf for years 2017–2020. Over Amery Ice Shelf, maximum lake extent varied strongly between the years with the 2019 melt season characterized by the largest areal coverage of supraglacial lakes (~763 km\(^2\)). The accuracy assessment over the test regions revealed an average Kappa coefficient of 0.86 where the largest value of Kappa reached 0.98 over George VI Ice Shelf. Future developments will involve the generation of circum-Antarctic supraglacial lake mapping products as well as their use for further methodological developments using Sentinel-1 SAR data in order to characterize intraannual supraglacial meltwater dynamics also during polar night and independent of meteorological conditions. In summary, the implementation of the Random Forest classifier enabled the development of the first automated mapping method applied to Sentinel-2 data distributed across all three Antarctic regions.
Land surface temperature (LST) is a fundamental parameter within the system of the Earth’s surface and atmosphere, which can be used to describe the inherent physical processes of energy and water exchange. The need for LST has been increasingly recognised in agriculture, as it affects the growth phases of crops and crop yields. However, challenges in overcoming the large discrepancies between the retrieved LST and ground truth data still exist. Precise LST measurement depends mainly on accurately deriving the surface emissivity, which is very dynamic due to changing states of land cover and plant development. In this study, we present an LST retrieval algorithm for the combined use of multispectral optical and thermal UAV images, which has been optimised for operational applications in agriculture to map the heterogeneous and diverse agricultural crop systems of a research campus in Germany (April 2018). We constrain the emissivity using certain NDVI thresholds to distinguish different land surface types. The algorithm includes atmospheric corrections and environmental thermal emissions to minimise the uncertainties. In the analysis, we emphasise that the omission of crucial meteorological parameters and inaccurately determined emissivities can lead to a considerably underestimated LST; however, if the emissivity is underestimated, the LST can be overestimated. The retrieved LST is validated by reference temperatures from nearby ponds and weather stations. The validation of the thermal measurements indicates a mean absolute error of about 0.5 K. The novelty of the dual sensor system is that it simultaneously captures highly spatially resolved optical and thermal images, in order to construct the precise LST ortho-mosaics required to monitor plant diseases and drought stress and validate airborne and satellite data.
Information about land use/land cover (LULC) and their changes is useful for different stakeholders to assess future pathways of sustainable land use for food production as well as for nature conservation. In this study, we assess LULC changes in the Kilombero catchment in Tanzania, an important area of recent development in East Africa. LULC change is assessed in two ways: first, post-classification comparison (PCC) which allows us to directly assess changes from one LULC class to another, and second, spectral change detection. We perform LULC classification by applying random forests (RF) on sets of multitemporal metrics that account for seasonal within-class dynamics. For the spectral change detection, we make use of the robust change vector analysis (RCVA) and determine those changes that do not necessarily lead to another class. The combination of the two approaches enables us to distinguish areas that show (a) only PCC changes, (b) only spectral changes that do not affect the classification of a pixel, (c) both types of change, or (d) no changes at all. Our results reveal that only one-quarter of the catchment has not experienced any change. One-third shows both, spectral changes and LULC conversion. Changes detected with both methods predominantly occur in two major regions, one in the West of the catchment, one in the Kilombero floodplain. Both regions are important areas of food production and economic development in Tanzania. The Kilombero floodplain is a Ramsar protected area, half of which was converted to agricultural land in the past decades. Therefore, LULC monitoring is required to support sustainable land management. Relatively poor classification performances revealed several challenges during the classification process. The combined approach of PCC and RCVA allows us to detect spatial patterns of LULC change at distinct dimensions and intensities. With the assessment of additional classifier output, namely class-specific per-pixel classification probabilities and derived parameters, we account for classification uncertainty across space. We overlay the LULC change results and the spatial assessment of classification reliability to provide a thorough picture of the LULC changes taking place in the Kilombero catchment.
Illegal small-scale mining (galamsey) in South-Western Ghana has grown tremendously in the last decade and caused significant environmental degradation. Excessive cloud cover in the area has limited the use of optical remote sensing data to map and monitor the extent of these activities. This study investigated the use of annual time-series Sentinel-1 data to map and monitor illegal mining activities along major rivers in South-Western Ghana between 2015 and 2019. A change detection approach, based on three time-series features — minimum, mean, maximum — was used to compute a backscatter threshold value suitable to identify/detect mining-induced land cover changes in the study area. Compared to the mean and maximum, the minimum time-series feature (in both VH and VV polarization) was found to be more sensitive to changes in backscattering within the period of investigation. Our approach permitted the detection of new illegal mining areas on an annual basis. A backscatter threshold value of +1.65 dB was found suitable for detecting illegal mining activities in the study area. Application of this threshold revealed illegal mining area extents of 102 km\(^2\), 60 km\(^2\) and 33 km\(^2\) for periods 2015/2016–2016/2017, 2016/2017–2017/2018 and 2017/2018–2018/2019, respectively. The observed decreasing trend in new illegal mining areas suggests that efforts at stopping illegal mining yielded positive results in the period investigated. Despite the advantages of Synthetic Aperture Radar data in monitoring phenomena in cloud-prone areas, our analysis revealed that about 25% of the Sentinel-1 data, mostly acquired in March and October (beginning and end of rainy season respectively), were unusable due to atmospheric effects from high intensity rainfall events. Further investigation in other geographies and climatic regions is needed to ascertain the susceptibility of Sentinel-1 data to atmospheric conditions.
Many parts of sub-Saharan Africa (SSA) are prone to land use and land cover change (LULCC). In many cases, natural systems are converted into agricultural land to feed the growing population. However, despite climate change being a major focus nowadays, the impacts of these conversions on water resources, which are essential for agricultural production, is still often neglected, jeopardizing the sustainability of the socio-ecological system. This study investigates historic land use/land cover (LULC) patterns as well as potential future LULCC and its effect on water quantities in a complex tropical catchment in Tanzania. It then compares the results using two climate change scenarios. The Land Change Modeler (LCM) is used to analyze and to project LULC patterns until 2030 and the Soil and Water Assessment Tool (SWAT) is utilized to simulate the water balance under various LULC conditions. Results show decreasing low flows by 6–8% for the LULC scenarios, whereas high flows increase by up to 84% for the combined LULC and climate change scenarios. The effect of climate change is stronger compared to the effect of LULCC, but also contains higher uncertainties. The effects of LULCC are more distinct, although crop specific effects show diverging effects on water balance components. This study develops a methodology for quantifying the impact of land use and climate change and therefore contributes to the sustainable management of the investigated catchment, as it shows the impact of environmental change on hydrological extremes (low flow and floods) and determines hot spots, which are critical for environmental development.
Large-area remote sensing time-series offer unique features for the extensive investigation of our environment. Since various error sources in the acquisition chain of datasets exist, only properly validated results can be of value for research and downstream decision processes. This review presents an overview of validation approaches concerning temporally dense time-series of land surface geo-information products that cover the continental to global scale. Categorization according to utilized validation data revealed that product intercomparisons and comparison to reference data are the conventional validation methods. The reviewed studies are mainly based on optical sensors and orientated towards global coverage, with vegetation-related variables as the focus. Trends indicate an increase in remote sensing-based studies that feature long-term datasets of land surface variables. The hereby corresponding validation efforts show only minor methodological diversification in the past two decades. To sustain comprehensive and standardized validation efforts, the provision of spatiotemporally dense validation data in order to estimate actual differences between measurement and the true state has to be maintained. The promotion of novel approaches can, on the other hand, prove beneficial for various downstream applications, although typically only theoretical uncertainties are provided.
Advances in remote inventory and analysis of forest resources during the last decade have reached a level to be now considered as a crucial complement, if not a surrogate, to the long-existing field-based methods. This is mostly reflected in not only the use of multiple-band new active and passive remote sensing data for forest inventory, but also in the methodic and algorithmic developments and/or adoptions that aim at maximizing the predictive or calibration performances, thereby minimizing both random and systematic errors, in particular for multi-scale spatial domains. With this in mind, this editorial note wraps up the recently-published Remote Sensing special issue “Remote Sensing-Based Forest Inventories from Landscape to Global Scale”, which hosted a set of state-of-the-art experiments on remotely sensed inventory of forest resources conducted by a number of prominent researchers worldwide.
Estimating penetration-related X-band InSAR elevation bias: a study over the Greenland ice sheet
(2019)
Accelerating melt on the Greenland ice sheet leads to dramatic changes at a global scale. Especially in the last decades, not only the monitoring, but also the quantification of these changes has gained considerably in importance. In this context, Interferometric Synthetic Aperture Radar (InSAR) systems complement existing data sources by their capability to acquire 3D information at high spatial resolution over large areas independent of weather conditions and illumination. However, penetration of the SAR signals into the snow and ice surface leads to a bias in measured height, which has to be corrected to obtain accurate elevation data. Therefore, this study purposes an easy transferable pixel-based approach for X-band penetration-related elevation bias estimation based on single-pass interferometric coherence and backscatter intensity which was performed at two test sites on the Northern Greenland ice sheet. In particular, the penetration bias was estimated using a multiple linear regression model based on TanDEM-X InSAR data and IceBridge laser-altimeter measurements to correct TanDEM-X Digital Elevation Model (DEM) scenes. Validation efforts yielded good agreement between observations and estimations with a coefficient of determination of R\(^2\) = 68% and an RMSE of 0.68 m. Furthermore, the study demonstrates the benefits of X-band penetration bias estimation within the application context of ice sheet elevation change detection.
Reliable near-surface soil moisture (θ) information is crucial for supporting risk assessment of future water usage, particularly considering the vulnerability of agroforestry systems of Mediterranean environments to climate change. We propose a simple empirical model by integrating dual-polarimetric Sentinel-1 (S1) Synthetic Aperture Radar (SAR) C-band single-look complex data and topographic information together with in-situ measurements of θ into a random forest (RF) regression approach (10-fold cross-validation). Firstly, we compare two RF models' estimation performances using either 43 SAR parameters (θNov\(^{SAR}\)) or the combination of 43 SAR and 10 terrain parameters (θNov\(^{SAR+Terrain}\)). Secondly, we analyze the essential parameters in estimating and mapping θ for S1 overpasses twice a day (at 5 a.m. and 5 p.m.) in a high spatiotemporal (17 × 17 m; 6 days) resolution. The developed site-specific calibration-dependent model was tested for a short period in November 2018 in a field-scale agroforestry environment belonging to the “Alento” hydrological observatory in southern Italy. Our results show that the combined SAR + terrain model slightly outperforms the SAR-based model (θNov\(^{SAR+Terrain}\) with 0.025 and 0.020 m3 m\(^{−3}\), and 89% compared to θNov\(^{SAR}\) with 0.028 and 0.022 m\(^3\) m\(^{−3}\, and 86% in terms of RMSE, MAE, and R2). The higher explanatory power for θNov\(^{SAR+Terrain}\) is assessed with time-variant SAR phase information-dependent elements of the C2 covariance and Kennaugh matrix (i.e., K1, K6, and K1S) and with local (e.g., altitude above channel network) and compound topographic attributes (e.g., wetness index). Our proposed methodological approach constitutes a simple empirical model aiming at estimating θ for rapid surveys with high accuracy. It emphasizes potentials for further improvement (e.g., higher spatiotemporal coverage of ground-truthing) by identifying differences of SAR measurements between S1 overpasses in the morning and afternoon.
Forest systems provide crucial ecosystem functions to our environment, such as balancing carbon stocks and influencing the local, regional and global climate. A trend towards an increasing frequency of climate change induced extreme weather events, including drought, is hereby a major challenge for forest management. Within this context, the application of remote sensing data provides a powerful means for fast, operational and inexpensive investigations over large spatial scales and time. This study was dedicated to explore the potential of satellite data in combination with harmonic analyses for quantifying the vegetation response to drought events in German forests. The harmonic modelling method was compared with a z-score standardization approach and correlated against both, meteorological and topographical data. Optical satellite imagery from Landsat and the Moderate Resolution Imaging Spectroradiometer (MODIS) was used in combination with three commonly applied vegetation indices. Highest correlation scores based on the harmonic modelling technique were computed for the 6th harmonic degree. MODIS imagery in combination with the Normalized Difference Vegetation Index (NDVI) generated hereby best results for measuring spectral response to drought conditions. Strongest correlation between remote sensing data and meteorological measures were observed for soil moisture and the self-calibrated Palmer Drought Severity Index (scPDSI). Furthermore, forests regions over sandy soils with pine as the dominant tree type were identified to be particularly vulnerable to drought. In addition, topographical analyses suggested mitigated drought affects along hill slopes. While the proposed approaches provide valuable information about vegetation dynamics as a response to meteorological weather conditions, standardized in-situ measurements over larger spatial scales and related to drought quantification are required for further in-depth quality assessment of the used methods and data.
Biological soil crusts (BSCs) are thin microbiological vegetation layers that naturally develop in unfavorable higher plant conditions (i.e., low precipitation rates and high temperatures) in global drylands. They consist of poikilohydric organisms capable of adjusting their metabolic activities depending on the water availability. However, they, and with them, their ecosystem functions, are endangered by climate change and land-use intensification. Remote sensing (RS)-based studies estimated the BSC cover in global drylands through various multispectral indices, and few of them correlated the BSCs’ activity response to rainfall. However, the allocation of BSCs is not limited to drylands only as there are areas beyond where smaller patches have developed under intense human impact and frequent disturbance. Yet, those areas were not addressed in RS-based studies, raising the question of whether the methods developed in extensive drylands can be transferred easily. Our temperate climate study area, the ‘Lieberoser Heide’ in northeastern Germany, is home to the country’s largest BSC-covered area. We applied a Random Forest (RF) classification model incorporating multispectral Sentinel-2 (S2) data, indices derived from them, and topographic information to spatiotemporally map the BSC cover for the first time in Central Europe. We further monitored the BSC response to rainfall events over a period of around five years (June 2015 to end of December 2020). Therefore, we combined datasets of gridded NDVI as a measure of photosynthetic activity with daily precipitation data and conducted a change detection analysis. With an overall accuracy of 98.9%, our classification proved satisfactory. Detected changes in BSC activity between dry and wet conditions were found to be significant. Our study emphasizes a high transferability of established methods from extensive drylands to BSC-covered areas in the temperate climate. Therefore, we consider our study to provide essential impulses so that RS-based biocrust mapping in the future will be applied beyond the global drylands.
Grasslands cover one third of the earth’s terrestrial surface and are mainly used for livestock production. The usage type, use intensity and condition of grasslands are often unclear. Remote sensing enables the analysis of grassland production and management on large spatial scales and with high temporal resolution. Despite growing numbers of studies in the field, remote sensing applications in grassland biomes are underrepresented in literature and less streamlined compared to other vegetation types. By reviewing articles within research on satellite-based remote sensing of grassland production traits and management, we describe and evaluate methods and results and reveal spatial and temporal patterns of existing work. In addition, we highlight research gaps and suggest research opportunities. The focus is on managed grasslands and pastures and special emphasize is given to the assessment of studies on grazing intensity and mowing detection based on earth observation data. Grazing and mowing highly influence the production and ecology of grassland and are major grassland management types. In total, 253 research articles were reviewed. The majority of these studies focused on grassland production traits and only 80 articles were about grassland management and use intensity. While the remote sensing-based analysis of grassland production heavily relied on empirical relationships between ground-truth and satellite data or radiation transfer models, the used methods to detect and investigate grassland management differed. In addition, this review identified that studies on grassland production traits with satellite data often lacked including spatial management information into the analyses. Studies focusing on grassland management and use intensity mostly investigated rather small study areas with homogeneous intensity levels among the grassland parcels. Combining grassland production estimations with management information, while accounting for the variability among grasslands, is recommended to facilitate the development of large-scale continuous monitoring and remote sensing grassland products, which have been rare thus far.
This study compares the performance of the five widely used crop growth models (CGMs): World Food Studies (WOFOST), Coalition for Environmentally Responsible Economies (CERES)-Wheat, AquaCrop, cropping systems simulation model (CropSyst), and the semi-empiric light use efficiency approach (LUE) for the prediction of winter wheat biomass on the Durable Environmental Multidisciplinary Monitoring Information Network (DEMMIN) test site, Germany. The study focuses on the use of remote sensing (RS) data, acquired in 2015, in CGMs, as they offer spatial information on the actual conditions of the vegetation. Along with this, the study investigates the data fusion of Landsat (30 m) and Moderate Resolution Imaging Spectroradiometer (MODIS) (500 m) data using the spatial and temporal reflectance adaptive reflectance fusion model (STARFM) fusion algorithm. These synthetic RS data offer a 30-m spatial and one-day temporal resolution. The dataset therefore provides the necessary information to run CGMs and it is possible to examine the fine-scale spatial and temporal changes in crop phenology for specific fields, or sub sections of them, and to monitor crop growth daily, considering the impact of daily climate variability. The analysis includes a detailed comparison of the simulated and measured crop biomass. The modelled crop biomass using synthetic RS data is compared to the model outputs using the original MODIS time series as well. On comparison with the MODIS product, the study finds the performance of CGMs more reliable, precise, and significant with synthetic time series. Using synthetic RS data, the models AquaCrop and LUE, in contrast to other models, simulate the winter wheat biomass best, with an output of high R2 (>0.82), low RMSE (<600 g/m\(^2\)) and significant p-value (<0.05) during the study period. However, inputting MODIS data makes the models underperform, with low R2 (<0.68) and high RMSE (>600 g/m\(^2\)). The study shows that the models requiring fewer input parameters (AquaCrop and LUE) to simulate crop biomass are highly applicable and precise. At the same time, they are easier to implement than models, which need more input parameters (WOFOST and CERES-Wheat).
Maize cropping systems mapping using RapidEye observations in agro-ecological landscapes in Kenya
(2017)
Cropping systems information on explicit scales is an important but rarely available variable in many crops modeling routines and of utmost importance for understanding pests and disease propagation mechanisms in agro-ecological landscapes. In this study, high spatial and temporal resolution RapidEye bio-temporal data were utilized within a novel 2-step hierarchical random forest (RF) classification approach to map areas of mono- and mixed maize cropping systems. A small-scale maize farming site in Machakos County, Kenya was used as a study site. Within the study site, field data was collected during the satellite acquisition period on general land use/land cover (LULC) and the two cropping systems. Firstly, non-cropland areas were masked out from other land use/land cover using the LULC mapping result. Subsequently an optimized RF model was applied to the cropland layer to map the two cropping systems (2nd classification step). An overall accuracy of 93% was attained for the LULC classification, while the class accuracies (PA: producer’s accuracy and UA: user’s accuracy) for the two cropping systems were consistently above 85%. We concluded that explicit mapping of different cropping systems is feasible in complex and highly fragmented agro-ecological landscapes if high resolution and multi-temporal satellite data such as 5 m RapidEye data is employed. Further research is needed on the feasibility of using freely available 10–20 m Sentinel-2 data for wide-area assessment of cropping systems as an important variable in numerous crop productivity models.
Accurate and detailed spatial soil information is essential for environmental modelling, risk assessment and decision making. The use of Remote Sensing data as secondary sources of information in digital soil mapping has been found to be cost effective and less time consuming compared to traditional soil mapping approaches. But the potentials of Remote Sensing data in improving knowledge of local scale soil information in West Africa have not been fully explored. This study investigated the use of high spatial resolution satellite data (RapidEye and Landsat), terrain/climatic data and laboratory analysed soil samples to map the spatial distribution of six soil properties–sand, silt, clay, cation exchange capacity (CEC), soil organic carbon (SOC) and nitrogen–in a 580 km2 agricultural watershed in south-western Burkina Faso. Four statistical prediction models–multiple linear regression (MLR), random forest regression (RFR), support vector machine (SVM), stochastic gradient boosting (SGB)–were tested and compared. Internal validation was conducted by cross validation while the predictions were validated against an independent set of soil samples considering the modelling area and an extrapolation area. Model performance statistics revealed that the machine learning techniques performed marginally better than the MLR, with the RFR providing in most cases the highest accuracy. The inability of MLR to handle non-linear relationships between dependent and independent variables was found to be a limitation in accurately predicting soil properties at unsampled locations. Satellite data acquired during ploughing or early crop development stages (e.g. May, June) were found to be the most important spectral predictors while elevation, temperature and precipitation came up as prominent terrain/climatic variables in predicting soil properties. The results further showed that shortwave infrared and near infrared channels of Landsat8 as well as soil specific indices of redness, coloration and saturation were prominent predictors in digital soil mapping. Considering the increased availability of freely available Remote Sensing data (e.g. Landsat, SRTM, Sentinels), soil information at local and regional scales in data poor regions such as West Africa can be improved with relatively little financial and human resources.
Detailed information on the land cover types present and the horizontal position of the land–water interface is needed for sensitive coastal ecosystems throughout the Arctic, both to establish baselines against which the impacts of climate change can be assessed and to inform response operations in the event of environmental emergencies such as oil spills. Previous work has demonstrated potential for accurate classification via fusion of optical and SAR data, though what contribution either makes to model accuracy is not well established, nor is it clear what shorelines can be classified using optical or SAR data alone. In this research, we evaluate the relative value of quad pol RADARSAT-2 and Landsat 5 data for shoreline mapping by individually excluding both datasets from Random Forest models used to classify images acquired over Nunavut, Canada. In anticipation of the RADARSAT Constellation Mission (RCM), we also simulate and evaluate dual and compact polarimetric imagery for shoreline mapping. Results show that SAR data is needed for accurate discrimination of substrates as user’s and producer’s accuracies were 5–24% higher for models constructed with quad pol RADARSAT-2 and DEM data than models constructed with Landsat 5 and DEM data. Models based on simulated RCM and DEM data achieved significantly lower overall accuracies (71–77%) than models based on quad pol RADARSAT-2 and DEM data (80%), with Wetland and Tundra being most adversely affected. When classified together with Landsat 5 and DEM data, however, model accuracy was less affected by the SAR data type, with multiple polarizations and modes achieving independent overall accuracies within a range acceptable for operational mapping, at 89–91%. RCM is expected to contribute positively to ongoing efforts to monitor change and improve emergency preparedness throughout the Arctic.
Monitoring the spatio-temporal development of vegetation is a challenging task in heterogeneous and cloud-prone landscapes. No single satellite sensor has thus far been able to provide consistent time series of high temporal and spatial resolution for such areas. In order to overcome this problem, data fusion algorithms such as the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM) have been established and frequently used in recent years to generate high-resolution time series. In order to make it applicable to larger scales and to increase the input data availability especially in cloud-prone areas, an ESTARFM framework was developed in this study introducing several enhancements. An automatic filling of cloud gaps was included in the framework to make best use of available, even partly cloud-covered Landsat images. Furthermore, the ESTARFM algorithm was enhanced to automatically account for regional differences in the heterogeneity of the study area. The generation of time series was automated and the processing speed was accelerated significantly by parallelization. To test the performance of the developed ESTARFM framework, MODIS and Landsat-8 data were fused for generating an 8-day NDVI time series for a study area of approximately 98,000 km\(^{2}\) in West Africa. The results show that the ESTARFM framework can accurately produce high temporal resolution time series (average MAE (mean absolute error) of 0.02 for the dry season and 0.05 for the vegetative season) while keeping the spatial detail in such a heterogeneous, cloud-prone region. The developments introduced within the ESTARFM framework establish the basis for large-scale research on various geoscientific questions related to land degradation, changes in land surface phenology or agriculture
Burkina Faso ranges amongst the fastest growing countries in the world with an annual population growth rate of more than three percent. This trend has consequences for food security since agricultural productivity is still on a comparatively low level in Burkina Faso. In order to compensate for the low productivity, the agricultural areas are expanding quickly. The mapping and monitoring of this expansion is difficult, even on the basis of remote sensing imagery, since the extensive farming practices and frequent cloud coverage in the area make the delineation of cultivated land from other land cover and land use types a challenging task. However, as the rapidly increasing population could have considerable effects on the natural resources and on the regional development of the country, methods for improved mapping of LULCC (land use and land cover change) are needed. For this study, we applied the newly developed ESTARFM (Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model) framework to generate high temporal (8-day) and high spatial (30 m) resolution NDVI time series for all of Burkina Faso for the years 2001, 2007, and 2014. For this purpose, more than 500 Landsat scenes and 3000 MODIS scenes were processed with this automated framework. The generated ESTARFM NDVI time series enabled extraction of per-pixel phenological features that all together served as input for the delineation of agricultural areas via random forest classification at 30 m spatial resolution for entire Burkina Faso and the three years. For training and validation, a randomly sampled reference dataset was generated from Google Earth images and based on expert knowledge. The overall accuracies of 92% (2001), 91% (2007), and 91% (2014) indicate the well-functioning of the applied methodology. The results show an expansion of agricultural area of 91% between 2001 and 2014 to a total of 116,900 km\(^2\). While rainfed agricultural areas account for the major part of this trend, irrigated areas and plantations also increased considerably, primarily promoted by specific development projects. This expansion goes in line with the rapid population growth in most provinces of Burkina Faso where land was still available for an expansion of agricultural area. The analysis of agricultural encroachment into protected areas and their surroundings highlights the increased human pressure on these areas and the challenges of environmental protection for the future.
Rice is an important food crop and a large producer of green-house relevant methane. Accurate and timely maps of paddy fields are most important in the context of food security and greenhouse gas emission modelling. During their life-cycle, rice plants undergo a phenological development that influences their interaction with waves in the visible light and infrared spectrum. Rice growth has a distinctive signature in time series of remotely-sensed data. We used time series of MODIS (Moderate Resolution Imaging Spectroradiometer) products MOD13Q1 and MYD13Q1 and a one-class support vector machine to detect these signatures and classify paddy rice areas in continental China. Based on these classifications, we present a novel product for continental China that shows rice areas for the years 2002, 2005, 2010 and 2014 at 250-m resolution. Our classification has an overall accuracy of 0.90 and a kappa coefficient of 0.77 compared to our own reference dataset for 2014 and correlates highly with rice area statistics from China’s Statistical Yearbooks (R2 of 0.92 for 2010, 0.92 for 2005 and 0.90 for 2002). Moderate resolution time series analysis allows accurate and timely mapping of rice paddies over large areas with diverse cropping schemes.
Sea level rise contribution from the Antarctic ice sheet is influenced by changes in glacier and ice shelf front position. Still, little is known about seasonal glacier and ice shelf front fluctuations as the manual delineation of calving fronts from remote sensing imagery is very time-consuming. The major challenge of automatic calving front extraction is the low contrast between floating glacier and ice shelf fronts and the surrounding sea ice. Additionally, in previous decades, remote sensing imagery over the often cloud-covered Antarctic coastline was limited. Nowadays, an abundance of Sentinel-1 imagery over the Antarctic coastline exists and could be used for tracking glacier and ice shelf front movement. To exploit the available Sentinel-1 data, we developed a processing chain allowing automatic extraction of the Antarctic coastline from Seninel-1 imagery and the creation of dense time series to assess calving front change. The core of the proposed workflow is a modified version of the deep learning architecture U-Net. This convolutional neural network (CNN) performs a semantic segmentation on dual-pol Sentinel-1 data and the Antarctic TanDEM-X digital elevation model (DEM). The proposed method is tested for four training and test areas along the Antarctic coastline. The automatically extracted fronts deviate on average 78 m in training and 108 m test areas. Spatial and temporal transferability is demonstrated on an automatically extracted 15-month time series along the Getz Ice Shelf. Between May 2017 and July 2018, the fronts along the Getz Ice Shelf show mostly an advancing tendency with the fastest moving front of DeVicq Glacier with 726 ± 20 m/yr.
This study investigates synthetic aperture radar (SAR) time series of the Sentinel-1 mission acquired over the Atacama Desert, Chile, between March 2015 and December 2018. The contribution analyzes temporal and spatial variations of Sentinel-1 interferometric SAR (InSAR) coherence and exemplarily illustrates factors that are responsible for observed signal differences. The analyses are based on long temporal baselines (365–1090 days) and temporally dense time series constructed with short temporal baselines (12–24 days). Results are compared to multispectral data of Sentinel-2, morphometric features of the digital elevation model (DEM) TanDEM-X WorldDEM™, and to a detailed governmental geographic information system (GIS) dataset of the local hydrography. Sentinel-1 datasets are suited for generating extensive, nearly seamless InSAR coherence mosaics covering the entire Atacama Desert (>450 × 1100 km) at a spatial resolution of 20 × 20 meter per pixel. Temporal baselines over several years lead only to very minor decorrelation, indicating a very high signal stability of C-Band in this region, especially in the hyperarid uplands between the Coastal Cordillera and the Central Depression. Signal decorrelation was associated with certain types of surface cover (e.g., water or aeolian deposits) or with actual surface dynamics (e.g., anthropogenic disturbance (mining) or fluvial activity and overland flow). Strong rainfall events and fluvial activity in the periods 2015 to 2016 and 2017 to 2018 caused spatial patterns with significant signal decorrelation; observed linear coherence anomalies matched the reference channel network and indicated actual episodic and sporadic discharge events. In the period 2015–2016, area-wide loss of coherence appeared as strip-like patterns of more than 80 km length that matched the prevailing wind direction. These anomalies, and others observed in that period and in the period 2017–2018, were interpreted to be caused by overland flow of high magnitude, as their spatial location matched well with documented heavy rainfall events that showed cumulative precipitation amounts of more than 20 mm.
Air temperatures in the Arctic have increased substantially over the last decades, which has extensively altered the properties of the land surface. Capturing the state and dynamics of Land Surface Temperatures (LSTs) at high spatial detail is of high interest as LST is dependent on a variety of surficial properties and characterizes the land–atmosphere exchange of energy. Accordingly, this study analyses the influence of different physical surface properties on the long-term mean of the summer LST in the Arctic Mackenzie Delta Region (MDR) using Landsat 30 m-resolution imagery between 1985 and 2018 by taking advantage of the cloud computing capabilities of the Google Earth Engine. Multispectral indices, including the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) and Tasseled Cap greenness (TCG), brightness (TCB), and wetness (TCW) as well as topographic features derived from the TanDEM-X digital elevation model are used in correlation and multiple linear regression analyses to reveal their influence on the LST. Furthermore, surface alteration trends of the LST, NDVI, and NDWI are revealed using the Theil-Sen (T-S) regression method. The results indicate that the mean summer LST appears to be mostly influenced by the topographic exposition as well as the prevalent moisture regime where higher evapotranspiration rates increase the latent heat flux and cause a cooling of the surface, as the variance is best explained by the TCW and northness of the terrain. However, fairly diverse model outcomes for different regions of the MDR (R2 from 0.31 to 0.74 and RMSE from 0.51 °C to 1.73 °C) highlight the heterogeneity of the landscape in terms of influential factors and suggests accounting for a broad spectrum of different factors when modeling mean LSTs. The T-S analysis revealed large-scale wetting and greening trends with a mean decadal increase of the NDVI/NDWI of approximately +0.03 between 1985 and 2018, which was mostly accompanied by a cooling of the land surface given the inverse relationship between mean LSTs and vegetation and moisture conditions. Disturbance through wildfires intensifies the surface alterations locally and lead to significantly cooler LSTs in the long-term compared to the undisturbed surroundings.
Forest biodiversity conservation requires precise, area-wide information on the abundance and distribution of key habitat structures at multiple spatial scales. We combined airborne laser scanning (ALS) data with color-infrared (CIR) aerial imagery for identifying individual tree characteristics and quantifying multi-scale habitat requirements using the example of the three-toed woodpecker (Picoides tridactylus) (TTW) in the Bavarian Forest National Park (Germany). This bird, a keystone species of boreal and mountainous forests, is highly reliant on bark beetles dwelling in dead or dying trees. While previous studies showed a positive relationship between the TTW presence and the amount of deadwood as a limiting resource, we hypothesized a unimodal response with a negative effect of very high deadwood amounts and tested for effects of substrate quality. Based on 104 woodpecker presence or absence locations, habitat selection was modelled at four spatial scales reflecting different woodpecker home range sizes. The abundance of standing dead trees was the most important predictor, with an increase in the probability of TTW occurrence up to a threshold of 44–50 dead trees per hectare, followed by a decrease in the probability of occurrence. A positive relationship with the deadwood crown size indicated the importance of fresh deadwood. Remote sensing data allowed both an area-wide prediction of species occurrence and the derivation of ecological threshold values for deadwood quality and quantity for more informed conservation management.
Inadequate land management and agricultural activities have largely resulted in land degradation in Burkina Faso. The nationwide governmental and institutional driven implementation and adoption of soil and water conservation measures (SWCM) since the early 1960s, however, is expected to successively slow down the degradation process and to increase the agricultural output. Even though relevant measures have been taken, only a few studies have been conducted to quantify their effect, for instance, on soil erosion and environmental restoration. In addition, a comprehensive summary of initiatives, implementation strategies, and eventually region-specific requirements for adopting different SWCM is missing. The present study therefore aims to review the different SWCM in Burkina Faso and implementation programs, as well as to provide information on their effects on environmental restoration and agricultural productivity. This was achieved by considering over 143 studies focusing on Burkina Faso’s experience and research progress in areas of SWCM and soil erosion. SWCM in Burkina Faso have largely resulted in an increase in agricultural productivity and improvement in food security. Finally, this study aims at supporting the country’s informed decision-making for extending already existing SWCM and for deriving further implementation strategies.
Forest ecosystems fulfill a whole host of ecosystem functions that are essential for life on our planet. However, an unprecedented level of anthropogenic influences is reducing the resilience and stability of our forest ecosystems as well as their ecosystem functions. The relationships between drivers, stress, and ecosystem functions in forest ecosystems are complex, multi-faceted, and often non-linear, and yet forest managers, decision makers, and politicians need to be able to make rapid decisions that are data-driven and based on short and long-term monitoring information, complex modeling, and analysis approaches. A huge number of long-standing and standardized forest health inventory approaches already exist, and are increasingly integrating remote-sensing based monitoring approaches. Unfortunately, these approaches in monitoring, data storage, analysis, prognosis, and assessment still do not satisfy the future requirements of information and digital knowledge processing of the 21st century. Therefore, this paper discusses and presents in detail five sets of requirements, including their relevance, necessity, and the possible solutions that would be necessary for establishing a feasible multi-source forest health monitoring network for the 21st century. Namely, these requirements are: (1) understanding the effects of multiple stressors on forest health; (2) using remote sensing (RS) approaches to monitor forest health; (3) coupling different monitoring approaches; (4) using data science as a bridge between complex and multidimensional big forest health (FH) data; and (5) a future multi-source forest health monitoring network. It became apparent that no existing monitoring approach, technique, model, or platform is sufficient on its own to monitor, model, forecast, or assess forest health and its resilience. In order to advance the development of a multi-source forest health monitoring network, we argue that in order to gain a better understanding of forest health in our complex world, it would be conducive to implement the concepts of data science with the components: (i) digitalization; (ii) standardization with metadata management after the FAIR (Findability, Accessibility, Interoperability, and Reusability) principles; (iii) Semantic Web; (iv) proof, trust, and uncertainties; (v) tools for data science analysis; and (vi) easy tools for scientists, data managers, and stakeholders for decision-making support.
Impervious surface areas (ISA) are heavily influenced by urban structure and related structural features. We examined the effects of object-based impervious surface spatial pattern analysis on land surface temperature and population density in Guangzhou, China, in comparison to classic per-pixel analyses. An object-based support vector machine (SVM) and a linear spectral mixture analysis (LSMA) were integrated to estimate ISA fraction using images from the Chinese HJ-1B satellite for 2009 to 2011. The results revealed that the integrated object-based SVM-LSMA algorithm outperformed the traditional pixel-wise LSMA algorithm in classifying ISA fraction. More specifically, the object-based ISA spatial patterns extracted were more suitable than pixel-wise patterns for urban heat island (UHI) studies, in which the UHI areas (landscape surface temperature >37 °C) generally feature high ISA fraction values (ISA fraction >50%). In addition, the object-based spatial patterns enable us to quantify the relationship of ISA with population density (correlation coefficient >0.2 in general), with global human settlement density (correlation coefficient >0.2), and with night-time light map (correlation coefficient >0.4), and, whereas pixel-wise ISA did not yield significant correlations. These results indicate that object-based spatial patterns have a high potential for UHI detection and urbanization monitoring. Planning measures that aim to reduce the urbanization impacts and UHI intensities can be better supported.
Regardless of political boundaries, river basins are a functional unit of the Earth’s land surface and provide an abundance of resources for the environment and humans. They supply livelihoods supported by the typical characteristics of large river basins, such as the provision of freshwater, irrigation water, and transport opportunities. At the same time, they are impacted i.e., by human-induced environmental changes, boundary conflicts, and upstream–downstream inequalities. In the framework of water resource management, monitoring of river basins is therefore of high importance, in particular for researchers, stake-holders and decision-makers. However, land surface and surface water properties of many major river basins remain largely unmonitored at basin scale. Several inventories exist, yet consistent spatial databases describing the status of major river basins at global scale are lacking. Here, Earth observation (EO) is a potential source of spatial information providing large-scale data on the status of land surface properties. This review provides a comprehensive overview of existing research articles analyzing major river basins primarily using EO. Furthermore, this review proposes to exploit EO data together with relevant open global-scale geodata to establish a database and to enable consistent spatial analyses and evaluate past and current states of major river basins.
Exploring the potential of C-Band SAR in contributing to burn severity mapping in tropical savanna
(2019)
The ability to map burn severity and to understand how it varies as a function of time of year and return frequency is an important tool for landscape management and carbon accounting in tropical savannas. Different indices based on optical satellite imagery are typically used for mapping fire scars and for estimating burn severity. However, cloud cover is a major limitation for analyses using optical data over tropical landscapes. To address this pitfall, we explored the suitability of C-band Synthetic Aperture Radar (SAR) data for detecting vegetation response to fire, using experimental fires in northern Australia. Pre- and post-fire results from Sentinel-1 C-band backscatter intensity data were compared to those of optical satellite imagery and were corroborated against structural changes on the ground that we documented through terrestrial laser scanning (TLS). Sentinel-1 C-band backscatter (VH) proved sensitive to the structural changes imparted by fire and was correlated with the Normalised Burn Ratio (NBR) derived from Sentinel-2 optical data. Our results suggest that C-band SAR holds potential to inform the mapping of burn severity in savannas, but further research is required over larger spatial scales and across a broader spectrum of fire regime conditions before automated products can be developed. Combining both Sentinel-1 SAR and Sentinel-2 multi-spectral data will likely yield the best results for mapping burn severity under a range of weather conditions.
The worldwide demand for food has been increasing due to the rapidly growing global population, and agricultural lands have increased in extent to produce more food crops. The pattern of cropland varies among different regions depending on the traditional knowledge of farmers and availability of uncultivated land. Satellite images can be used to map cropland in open areas but have limitations for detecting undergrowth inside forests. Classification results are often biased and need to be supplemented with field observations. Undercover cropland inside forests in the Bale Mountains of Ethiopia was assessed using field observed percentage cover of land use/land cover classes, and topographic and location parameters. The most influential factors were identified using Boosted Regression Trees and used to map undercover cropland area. Elevation, slope, easterly aspect, distance to settlements, and distance to national park were found to be the most influential factors determining undercover cropland area. When there is very high demand for growing food crops, constrained under restricted rights for clearing forest, cultivation could take place within forests as an undercover. Further research on the impact of undercover cropland on ecosystem services and challenges in sustainable management is thus essential.
Background:
Schistosomiasis is a water-based disease that affects an estimated 250 million people, mainly in sub-Saharan Africa. The transmission of schistosomiasis is spatially and temporally restricted to freshwater bodies that contain schistosome cercariae released from specific snails that act as intermediate hosts. Our objective was to assess the contribution of remote sensing applications and to identify remaining challenges in its optimal application for schistosomiasis risk profiling in order to support public health authorities to better target control interventions.
Methods:
We reviewed the literature (i) to deepen our understanding of the ecology and the epidemiology of schistosomiasis, placing particular emphasis on remote sensing; and (ii) to fill an identified gap, namely interdisciplinary research that bridges different strands of scientific inquiry to enhance spatially explicit risk profiling. As a first step, we reviewed key factors that govern schistosomiasis risk. Secondly, we examined remote sensing data and variables that have been used for risk profiling of schistosomiasis. Thirdly, the linkage between the ecological consequence of environmental conditions and the respective measure of remote sensing data were synthesised.
Results:
We found that the potential of remote sensing data for spatial risk profiling of schistosomiasis is - in principle - far greater than explored thus far. Importantly though, the application of remote sensing data requires a tailored approach that must be optimised by selecting specific remote sensing variables, considering the appropriate scale of observation and modelling within ecozones. Interestingly, prior studies that linked prevalence of Schistosoma infection to remotely sensed data did not reflect that there is a spatial gap between the parasite and intermediate host snail habitats where disease transmission occurs, and the location (community or school) where prevalence measures are usually derived from.
Conclusions:
Our findings imply that the potential of remote sensing data for risk profiling of schistosomiasis and other neglected tropical diseases has yet to be fully exploited.
Disentangling the relative effects of bushmeat availability on human nutrition in central Africa
(2015)
We studied links between human malnutrition and wild meat availability within the Rainforest Biotic Zone in central Africa. We distinguished two distinct hunted mammalian diversity distributions, one in the rainforest areas (Deep Rainforest Diversity, DRD) containing taxa of lower hunting sustainability, the other in the northern rainforest-savanna mosaic, with species of greater hunting potential (Marginal Rainforest Diversity, MRD). Wild meat availability, assessed by standing crop mammalian biomass, was greater in MRD than in DRD areas. Predicted bushmeat extraction was also higher in MRD areas. Despite this, stunting of children, a measure of human malnutrition, was greater in MRD areas. Structural equation modeling identified that, in MRD areas, mammal diversity fell away from urban areas, but proximity to these positively influenced higher stunting incidence. In DRD areas, remoteness and distance from dense human settlements and infrastructures explained lower stunting levels. Moreover, stunting was higher away from protected areas. Our results suggest that in MRD areas, forest wildlife rational use for better human nutrition is possible. By contrast, the relatively low human populations in DRD areas currently offer abundant opportunities for the continued protection of more vulnerable mammals and allow dietary needs of local populations to be met.
The focus of this analysis is on the early detection of forest health changes, specifically that of Norway spruce (Picea abies L. Karst.). In this analysis, we planned to examine the time (degree of early detection), spectral wavelengths and appropriate method for detecting vitality changes. To accomplish this, a ring-barking experiment with seven subsequent laboratory needle measurements was carried out in 2013 and 2014 in an area in southeastern Germany near Altötting. The experiment was also accompanied by visual crown condition assessment. In total, 140 spruce trees in groups of five were ring-barked with the same number of control trees in groups of five that were selected as reference trees in order to compare their development. The laboratory measurements were analysed regarding the separability of ring-barked and control samples using spectral reflectance, vegetation indices and derivative analysis. Subsequently, a random forest classifier for determining important spectral wavelength regions was applied. Results from the methods are consistent and showed a high importance of the visible (VIS) spectral region, very low importance of the near-infrared (NIR) and minor importance of the shortwave infrared (SWIR) spectral region. Using spectral reflectance data as well as indices, the earliest separation time was found to be 292 days after ring-barking. The derivative analysis showed that a significant separation was observed 152 days after ring-barking for six spectral features spread through VIS and SWIR. A significant separation was detected using a random forest classifier 292 days after ring-barking with 58% separability. The visual crown condition assessment was analysed regarding obvious changes of vitality and the first indication was observed 302 days after ring-barking as bark beetle infestation and yellowing of foliage in the ring-barked trees only. This experiment shows that an early detection, compared with visual crown assessment, is possible using the proposed methods for this specific data set. This study will contribute to ongoing research for early detection of vitality changes that will support foresters and decision makers.
Beyond maximum independent set: an extended integer programming formulation for point labeling
(2017)
Map labeling is a classical problem of cartography that has frequently been approached by combinatorial optimization. Given a set of features in a map and for each feature a set of label candidates, a common problem is to select an independent set of labels (that is, a labeling without label–label intersections) that contains as many labels as possible and at most one label for each feature. To obtain solutions of high cartographic quality, the labels can be weighted and one can maximize the total weight (rather than the number) of the selected labels. We argue, however, that when maximizing the weight of the labeling, the influences of labels on other labels are insufficiently addressed. Furthermore, in a maximum-weight labeling, the labels tend to be densely packed and thus the map background can be occluded too much. We propose extensions of an existing model to overcome these limitations. Since even without our extensions the problem is NP-hard, we cannot hope for an efficient exact algorithm for the problem. Therefore, we present a formalization of our model as an integer linear program (ILP). This allows us to compute optimal solutions in reasonable time, which we demonstrate both for randomly generated point sets and an existing data set of cities. Moreover, a relaxation of our ILP allows for a simple and efficient heuristic, which yielded near-optimal solutions for our instances.
In this study, polarimetric Synthetic Aperture Radar (PolSAR) data at X-, C- and L-Bands, acquired by the satellites: TerraSAR-X (2011), Radarsat-2 (2011), ALOS (2010) and ALOS-2 (2016), were used to characterize the tundra land cover of a test site located close to the town of Tuktoyaktuk, NWT, Canada. Using available in situ ground data collected in 2010 and 2012, we investigate PolSAR scattering characteristics of common tundra land cover classes at X-, C- and L-Bands. Several decomposition features of quad-, co-, and cross-polarized data were compared, the correlation between them was investigated, and the class separability offered by their different feature spaces was analyzed. Certain PolSAR features at each wavelength were sensitive to the land cover and exhibited distinct scattering characteristics. Use of shorter wavelength imagery (X and C) was beneficial for the characterization of wetland and tundra vegetation, while L-Band data highlighted differences of the bare ground classes better. The Kennaugh Matrix decomposition applied in this study provided a unified framework to store, process, and analyze all data consistently, and the matrix offered a favorable feature space for class separation. Of all elements of the quad-polarized Kennaugh Matrix, the intensity based elements K0, K1, K2, K3 and K4 were found to be most valuable for class discrimination. These elements contributed to better class separation as indicated by an increase of the separability metrics squared Jefferys Matusita Distance and Transformed Divergence. The increase in separability was up to 57% for Radarsat-2 and up to 18% for ALOS-2 data.
Interactions between different formative processes are reflected in the internal structure of rock glaciers. Therefore, the detection of subsurface conditions can help to enhance our understanding of landform development. For an assessment of subsurface conditions, we present an analysis of the spatial variability of active layer thickness, ground ice content and frost table topography for two different rock glaciers in the Eastern Swiss Alps by means of quasi-3-D electrical resistivity imaging (ERI). This approach enables an extensive mapping of subsurface structures and a spatial overlay between site-specific surface and subsurface characteristics. At Nair rock glacier, we discovered a gradual descent of the frost table in a downslope direction and a constant decrease of ice content which follows the observed surface topography. This is attributed to ice formation by refreezing meltwater from an embedded snow bank or from a subsurface ice patch which reshapes the permafrost layer. The heterogeneous ground ice distribution at Uertsch rock glacier indicates that multiple processes on different time domains were involved in the development. Resistivity values which represent frozen conditions vary within a wide range and indicate a successive formation which includes several advances, past glacial overrides and creep processes on the rock glacier surface. In combination with the observed topography, quasi-3-D ERI enables us to delimit areas of extensive and compressive flow in close proximity. Excellent data quality was provided by a good coupling of electrodes to the ground in the pebbly material of the investigated rock glaciers. Results show the value of the quasi-3-D ERI approach but advise the application of complementary geophysical methods for interpreting the results.
Home on the Range: Factors Explaining Partial Migration of African Buffalo in a Tropical Environment
(2012)
Partial migration (when only some individuals in a population undertake seasonal migrations) is common in many species and geographical contexts. Despite the development of modern statistical methods for analyzing partial migration, there have been no studies on what influences partial migration in tropical environments. We present research on factors affecting partial migration in African buffalo (Syncerus caffer) in northeastern Namibia. Our dataset is derived from 32 satellite tracking collars, spans 4 years and contains over 35,000 locations. We used remotely sensed data to quantify various factors that buffalo experience in the dry season when making decisions on whether and how far to migrate, including potential man-made and natural barriers, as well as spatial and temporal heterogeneity in environmental conditions. Using an information-theoretic, non-linear regression approach, our analyses showed that buffalo in this area can be divided into 4 migratory classes: migrants, non-migrants, dispersers, and a new class that we call "expanders". Multimodel inference from least-squares regressions of wet season movements showed that environmental conditions (rainfall, fires, woodland cover, vegetation biomass), distance to the nearest barrier (river, fence, cultivated area) and social factors (age, size of herd at capture) were all important in explaining variation in migratory behaviour. The relative contributions of these variables to partial migration have not previously been assessed for ungulates in the tropics. Understanding the factors driving migratory decisions of wildlife will lead to better-informed conservation and land-use decisions in this area.
Robust risk assessment requires accurate flood intensity area mapping to allow for the identification of populations and elements at risk. However, available flood maps in West Africa lack spatial variability while global datasets have resolutions too coarse to be relevant for local scale risk assessment. Consequently, local disaster managers are forced to use traditional methods such as watermarks on buildings and media reports to identify flood hazard areas. In this study, remote sensing and Geographic Information System (GIS) techniques were combined with hydrological and statistical models to delineate the spatial limits of flood hazard zones in selected communities in Ghana, Burkina Faso and Benin. The approach involves estimating peak runoff concentrations at different elevations and then applying statistical methods to develop a Flood Hazard Index (FHI). Results show that about half of the study areas fall into high intensity flood zones. Empirical validation using statistical confusion matrix and the principles of Participatory GIS show that flood hazard areas could be mapped at an accuracy ranging from 77% to 81%. This was supported with local expert knowledge which accurately classified 79% of communities deemed to be highly susceptible to flood hazard. The results will assist disaster managers to reduce the risk to flood disasters at the community level where risk outcomes are first materialized.
River deltas belong to the most densely settled places on earth. Although they only account for 5% of the global land surface, over 550 million people live in deltas. These preferred livelihood locations, which feature flat terrain, fertile alluvial soils, access to fluvial and marine resources, a rich wetland biodiversity and other advantages are, however, threatened by numerous internal and external processes. Socio-economic development, urbanization, climate change induced sea level rise, as well as flood pulse changes due to upstream water diversion all lead to changes in these highly dynamic systems. A thorough understanding of a river delta's general setting and intra-annual as well as long-term dynamic is therefore crucial for an informed management of natural resources. Here, remote sensing can play a key role in analyzing and monitoring these vast areas at a global scale. The goal of this study is to demonstrate the potential of intra-annual time series analyses at dense temporal, but coarse spatial resolution for inundation characterization in five river deltas located in four different countries. Based on 250 m MODIS reflectance data we analyze inundation dynamics in four densely populated Asian river deltas-namely the Yellow River Delta (China), the Mekong Delta (Vietnam), the Irrawaddy Delta (Myanmar), and the Ganges-Brahmaputra (Bangladesh, India)-as well as one very contrasting delta: the nearly uninhabited polar Mackenzie Delta Region in northwestern Canada for the complete time span of one year (2013). A complex processing chain of water surface derivation on a daily basis allows the generation of intra-annual time series, which indicate inundation duration in each of the deltas. Our analyses depict distinct inundation patterns within each of the deltas, which can be attributed to processes such as overland flooding, irrigation agriculture, aquaculture, or snowmelt and thermokarst processes. Clear differences between mid-latitude, subtropical, and polar deltas are illustrated, and the advantages and limitations of the approach for inundation derivation are discussed.
Due to their negative water budget most recent semi-/arid regions are characterized by vast evaporates (salt lakes and salty soils). We recently identified those hyper-saline environments as additional sources for a multitude of volatile halogenated organohalogens (VOX). These compounds can affect the ozone layer of the stratosphere and play a key role in the production of aerosols. A remote sensing based analysis was performed in the Southern Aral Sea basin, providing information of major soil types as well as their extent and spatial and temporal evolution. VOX production has been determined in dry and moist soil samples after 24 h. Several C1- and C2 organohalogens have been found in hyper-saline topsoil profiles, including CH3Cl, CH3Br, CHBr3 and CHCl3. The range of organohalogens also includes trans-1,2-dichloroethene (DCE), which is reported here to be produced naturally for the first time. Using MODIS time series and supervised image classification a daily production rate for DCE has been calculated for the 15 000 km\(^2\) ranging research area in the southern Aralkum. The applied laboratory setup simulates a short-term change in climatic conditions, starting from dried-out saline soil that is instantly humidified during rain events or flooding. It describes the general VOX production potential, but allows only for a rough estimation of resulting emission loads. VOX emissions are expected to increase in the future since the area of salt affected soils is expanding due to the regressing Aral Sea. Opportunities, limits and requirements of satellite based rapid change detection and salt classification are discussed.
This study investigates a two component decomposition technique for HH/VV-polarized PolSAR (Polarimetric Synthetic Aperture Radar) data. The approach is a straight forward adaption of the Yamaguchi decomposition and decomposes the data into two scattering contributions: surface and double bounce under the assumption of a negligible vegetation scattering component in Tundra environments. The dependencies between the features of this two and the classical three component Yamaguchi decomposition were investigated for Radarsat-2 (quad) and TerraSAR-X (HH/VV) data for the Mackenzie Delta Region, Canada. In situ data on land cover were used to derive the scattering characteristics and to analyze the correlation among the PolSAR features. The double bounce and surface scattering features of the two and three component scattering model (derived from pseudo-HH/VV- and quad-polarized data) showed similar scattering characteristics and positively correlated-R2 values of 0.60 (double bounce) and 0.88 (surface scattering) were observed. The presence of volume scattering led to differences between the features and these were minimized for land cover classes of low vegetation height that showed little volume scattering contribution. In terms of separability, the quad-polarized Radarsat-2 data offered the best separation of the examined tundra land cover types and will be best suited for the classification. This is anticipated as it represents the largest feature space of all tested ones. However; the classes “wetland” and “bare ground” showed clear positions in the feature spaces of the C- and X-Band HH/VV-polarized data and an accurate classification of these land cover types is promising. Among the possible dual-polarization modes of Radarsat-2 the HH/VV was found to be the favorable mode for the characterization of the aforementioned tundra land cover classes due to the coherent acquisition and the preserved co-pol. phase. Contrary, HH/HV-polarized and VV/VH-polarized data were found to be best suited for the characterization of mixed and shrub dominated tundra.
The Sentinel-1 Satellite (S-1) of ESA's Copernicus Mission delivers freely available C-Band Synthetic Aperture Radar (SAR) data that are suited for interferometric applications (InSAR). The high geometric resolution of less than fifteen meter and the large coverage offered by the Interferometric Wide Swath mode (IW) point to new perspectives on the comprehension and understanding of surface changes, the quantification and monitoring of dynamic processes, especially in arid regions. The contribution shows the application of S-1 intensities and InSAR coherences in time series analysis for the delineation of changes related to fluvial morphodynamics in Damghan, Iran. The investigations were carried out for the period from April to October 2015 and exhibit the potential of the S-1 data for the identification of surface disturbances, mass movements and fluvial channel activity in the surroundings of the Damghan Playa. The Amplitude Change Detection highlighted extensive material movement and accumulation - up to sizes of more than 4,000 m in width - in the east of the Playa via changes in intensity. Further, the Coherence Change Detection technique was capable to indicate small-scale channel activity of the drainage system that was neither recognizable in the S-1 intensity nor the multispectral Landsat-8 data. The run off caused a decorrelation of the SAR signals and a drop in coherence. Seen from a morphodynamic point of view, the results indicated a highly dynamic system and complex tempo-spatial patterns were observed that will be subject of future analysis. Additionally, the study revealed the necessity to collect independent reference data on fluvial activity in order to train and adjust the change detector.
Cropping Intensity in the Aral Sea Basin and Its Dependency from the Runoff Formation 2000–2012
(2016)
This study is aimed at a better understanding of how upstream runoff formation affected the cropping intensity (CI: number of harvests) in the Aral Sea Basin (ASB) between 2000 and 2012. MODIS 250 m NDVI time series and knowledge-based pixel masking that included settlement layers and topography features enabled to map the irrigated cropland extent (iCE). Random forest models supported the classification of cropland vegetation phenology (CVP: winter/summer crops, double cropping, etc.). CI and the percentage of fallow cropland (PF) were derived from CVP. Spearman’s rho was selected for assessing the statistical relation of CI and PF to runoff formation in the Amu Darya and Syr Darya catchments per hydrological year. Validation in 12 reference sites using multi-annual Landsat-7 ETM+ images revealed an average overall accuracy of 0.85 for the iCE maps. MODIS maps overestimated that based on Landsat by an average factor of ~1.15 (MODIS iCE/Landsat iCE). Exceptional overestimations occurred in case of inaccurate settlement layers. The CVP and CI maps achieved overall accuracies of 0.91 and 0.96, respectively. The Amu Darya catchment disclosed significant positive (negative) relations between upstream runoff with CI (PF) and a high pressure on the river water resources in 2000–2012. Along the Syr Darya, reduced dependencies could be observed, which is potentially linked to the high number of water constructions in that catchment. Intensified double cropping after drought years occurred in Uzbekistan. However, a 10 km × 10 km grid of Spearman’s rho (CI and PF vs. upstream runoff) emphasized locations at different CI levels that are directly affected by runoff fluctuations in both river systems. The resulting maps may thus be supportive on the way to achieve long-term sustainability of crop production and to simultaneously protect the severely threatened environment in the ASB. The gained knowledge can be further used for investigating climatic impacts of irrigation in the region.
This study analyzed the spatiotemporal pattern of settlement expansion in Abuja, Nigeria, one of West Africa’s fastest developing cities, using geoinformation and ancillary datasets. Three epochs of Land-use Land-cover (LULC) maps for 1986, 2001 and 2014 were derived from Landsat images using support vector machines (SVM). Accuracy assessment (AA) of the LULC maps based on the pixel count resulted in overall accuracy of 82%, 92% and 92%, while the AA derived from the error adjusted area (EAA) method stood at 69%, 91% and 91% for 1986, 2001 and 2014, respectively. Two major techniques for detecting changes in the LULC epochs involved the use of binary maps as well as a post-classification comparison approach. Quantitative spatiotemporal analysis was conducted to detect LULC changes with specific focus on the settlement development pattern of Abuja, the federal capital city (FCC) of Nigeria. Logical transitions to the urban category were modelled for predicting future scenarios for the year 2050 using the embedded land change modeler (LCM) in the IDRISI package. Based on the EAA, the result showed that urban areas increased by more than 11% between 1986 and 2001. In contrast, this value rose to 17% between 2001 and 2014. The LCM model projected LULC changes that showed a growing trend in settlement expansion, which might take over allotted spaces for green areas and agricultural land if stringent development policies and enforcement measures are not implemented. In conclusion, integrating geospatial technologies with ancillary datasets offered improved understanding of how urbanization processes such as increased imperviousness of such a magnitude could influence the urban microclimate through the alteration of natural land surface temperature. Urban expansion could also lead to increased surface runoff as well as changes in drainage geography leading to urban floods.
Advancing land degradation in the irrigated areas of Central Asia hinders sustainable development of this predominantly agricultural region. To support decisions on mitigating cropland degradation, this study combines linear trend analysis and spatial logistic regression modeling to expose a land degradation trend in the Khorezm region, Uzbekistan, and to analyze the causes. Time series of the 250-m MODIS NDVI, summed over the growing seasons of 2000–2010, were used to derive areas with an apparent negative vegetation trend; this was interpreted as an indicator of land degradation. About one third (161,000 ha) of the region’s area experienced negative trends of different magnitude. The vegetation decline was particularly evident on the low-fertility lands bordering on the natural sandy desert, suggesting that these areas should be prioritized in mitigation planning. The results of logistic modeling indicate that the spatial pattern of the observed trend is mainly associated with the level of the groundwater table (odds = 330 %), land-use intensity (odds = 103 %), low soil quality (odds = 49 %), slope (odds = 29 %), and salinity of the groundwater (odds = 26 %). Areas, threatened by land degradation, were mapped by fitting the estimated model parameters to available data. The elaborated approach, combining remote-sensing and GIS, can form the basis for developing a common tool for monitoring land degradation trends in irrigated croplands of Central Asia.
Rice is the most important food crop in Asia, and the timely mapping and monitoring of paddy rice fields subsequently emerged as an important task in the context of food security and modelling of greenhouse gas emissions. Rice growth has a distinct influence on Synthetic Aperture Radar (SAR) backscatter images, and time-series analysis of C-band images has been successfully employed to map rice fields. The poor data availability on regional scales is a major drawback of this method. We devised an approach to classify paddy rice with the use of all available Envisat ASAR WSM (Advanced Synthetic Aperture Radar Wide Swath Mode) data for our study area, the Mekong Delta in Vietnam. We used regression-based incidence angle normalization and temporal averaging to combine acquisitions from multiple tracks and years. A crop phenology-based classifier has been applied to this time series to detect single-, double- and triple-cropped rice areas (one to three harvests per year), as well as dates and lengths of growing seasons. Our classification has an overall accuracy of 85.3% and a kappa coefficient of 0.74 compared to a reference dataset and correlates highly with official rice area statistics at the provincial level (R-2 of 0.98). SAR-based time-series analysis allows accurate mapping and monitoring of rice areas even under adverse atmospheric conditions.
An Overview of the Regional Experiments for Land-atmosphere Exchanges 2012 (REFLEX 2012) Campaign
(2015)
The REFLEX 2012 campaign was initiated as part of a training course on the organization of an airborne campaign to support advancement of the understanding of land-atmosphere interaction processes. This article describes the campaign, its objectives and observations, remote as well as in situ. The observations took place at the experimental Las Tiesas farm in an agricultural area in the south of Spain. During the period of ten days, measurements were made to capture the main processes controlling the local and regional land-atmosphere exchanges. Apart from multi-temporal, multi-directional and multi-spatial space-borne and airborne observations, measurements of the local meteorology, energy fluxes, soil temperature profiles, soil moisture profiles, surface temperature, canopy structure as well as leaf-level measurements were carried out. Additional thermo-dynamical monitoring took place at selected sites. After presenting the different types of measurements, some examples are given to illustrate the potential of the observations made.
This study aimed to optimise the application, efficiency and interpretability of quasi-3D resistivity imaging for investigating the heterogeneous permafrost distribution at mountain sites by a systematic forward modelling approach. A three dimensional geocryologic model, representative for most mountain permafrost settings, was developed. Based on this geocryologic model quasi-3D models were generated by collating synthetic orthogonal 2D arrays, demonstrating the effects of array types and electrode spacing on resolution and interpretability of the inversion results. The effects of minimising the number of 2D arrays per quasi-3D grid were tested by enlarging the spacing between adjacent lines and by reducing the number of perpendicular tie lines with regard to model resolution and loss of information value. Synthetic and measured quasi-3D models were investigated with regard to the lateral and vertical resolution, reliability of inverted resistivity values, the possibility of a quantitative interpretation of resistivities and the response of the inversion process on the validity of quasi-3D models. Results show that setups using orthogonal 2D arrays with electrode spacings of 2 m and 3 m are capable of delineating lateral heterogeneity with high accuracy and also deliver reliable data on active layer thickness. Detection of permafrost thickness, especially if the permafrost base is close to the penetration depth of the setups, and the reliability of absolute resistivity values emerged to be a weakness of the method. Quasi-3D imaging has proven to be a promising tool for investigating permafrost in mountain environments especially for delineating the often small-scale permafrost heterogeneity, and therefore provides an enhanced possibility for aligning permafrost distribution with site specific surface properties and morphological settings.
The natural environment and livelihoods in the Lower Mekong Basin (LMB) are significantly affected by the annual hydrological cycle. Monitoring of soil moisture as a key variable in the hydrological cycle is of great interest in a number of Hydrological and agricultural applications. In this study we evaluated the quality and spatiotemporal variability of the soil moisture product retrieved from C-band scatterometers data across the LMB sub-catchments. The soil moisture retrieval algorithm showed reasonable performance in most areas of the LMB with the exception of a few sub-catchments in the eastern parts of Laos, where the land cover is characterized by dense vegetation. The best performance of the retrieval algorithm was obtained in agricultural regions. Comparison of the available in situ evaporation data in the LMB and the Basin Water Index (BWI), an indicator of the basin soil moisture condition, showed significant negative correlations up to R = −0.85. The inter-annual variation of the calculated BWI was also found corresponding to the reported extreme hydro-meteorological events in the Mekong region. The retrieved soil moisture data show high correlation (up to R = 0.92) with monthly anomalies of precipitation in non-irrigated regions. In general, the seasonal variability of soil moisture in the LMB was well captured by the retrieval method. The results of analysis also showed significant correlation between El Niño events and the monthly BWI anomaly measurements particularly for the month May with the maximum correlation of R = 0.88.
Mobile laser scanning puts high requirements on the accuracy of the positioning systems and the calibration of the measurement system. We present a novel algorithmic approach for calibration with the goal of improving the measurement accuracy of mobile laser scanners. We describe a general framework for calibrating mobile sensor platforms that estimates all configuration parameters for any arrangement of positioning sensors, including odometry. In addition, we present a novel semi-rigid Simultaneous Localization and Mapping (SLAM) algorithm that corrects the vehicle position at every point in time along its trajectory, while simultaneously improving the quality and precision of the entire acquired point cloud. Using this algorithm, the temporary failure of accurate external positioning systems or the lack thereof can be compensated for. We demonstrate the capabilities of the two newly proposed algorithms on a wide variety of datasets.
Recently, several backpack-mounted systems, also known as personal laser scanning systems, have been developed. They consist of laser scanners or cameras that are carried by a human operator to acquire measurements of the environment while walking. These systems were first designed to overcome the challenges of mapping indoor environments with doors and stairs. While the human operator inherently has the ability to open doors and to climb stairs, the flexible movements introduce irregularities of the trajectory to the system. To compete with other mapping systems, the accuracy of these systems has to be evaluated. In this paper, we present an extensive evaluation of our backpack mobile mapping system in indoor environments. It is shown that the system can deal with the normal human walking motion, but has problems with irregular jittering. Moreover, we demonstrate the applicability of the backpack in a suitable urban scenario.
Background
Schistosomiasis is the most widespread water-based disease in sub-Saharan Africa. Transmission is governed by the spatial distribution of specific freshwater snails that act as intermediate hosts and human water contact patterns. Remote sensing data have been utilized for spatially explicit risk profiling of schistosomiasis. We investigated the potential of remote sensing to characterize habitat conditions of parasite and intermediate host snails and discuss the relevance for public health.
Methodology
We employed high-resolution remote sensing data, environmental field measurements, and ecological data to model environmental suitability for schistosomiasis-related parasite and snail species. The model was developed for Burkina Faso using a habitat suitability index (HSI). The plausibility of remote sensing habitat variables was validated using field measurements. The established model was transferred to different ecological settings in Côte d’Ivoire and validated against readily available survey data from school-aged children.
Principal Findings
Environmental suitability for schistosomiasis transmission was spatially delineated and quantified by seven habitat variables derived from remote sensing data. The strengths and weaknesses highlighted by the plausibility analysis showed that temporal dynamic water and vegetation measures were particularly useful to model parasite and snail habitat suitability, whereas the measurement of water surface temperature and topographic variables did not perform appropriately. The transferability of the model showed significant relations between the HSI and infection prevalence in study sites of Côte d’Ivoire.
Conclusions/Significance
A predictive map of environmental suitability for schistosomiasis transmission can support measures to gain and sustain control. This is particularly relevant as emphasis is shifting from morbidity control to interrupting transmission. Further validation of our mechanistic model needs to be complemented by field data of parasite- and snail-related fitness. Our model provides a useful tool to monitor the development of new hotspots of potential schistosomiasis transmission based on regularly updated remote sensing data.
Accurate quantification of land use/cover change (LULCC) is important for efficient environmental management, especially in regions that are extremely affected by climate variability and continuous population growth such as West Africa. In this context, accurate LULC classification and statistically sound change area estimates are essential for a better understanding of LULCC processes. This study aimed at comparing mono-temporal and multi-temporal LULC classifications as well as their combination with ancillary data and to determine LULCC across the heterogeneous landscape of southwest Burkina Faso using accurate classification results. Landsat data (1999, 2006 and 2011) and ancillary data served as input features for the random forest classifier algorithm. Five LULC classes were identified: woodland, mixed vegetation, bare surface, water and agricultural area. A reference database was established using different sources including high-resolution images, aerial photo and field data. LULCC and LULC classification accuracies, area and area uncertainty were computed based on the method of adjusted error matrices. The results revealed that multi-temporal classification significantly outperformed those solely based on mono-temporal data in the study area. However, combining mono-temporal imagery and ancillary data for LULC classification had the same accuracy level as multi-temporal classification which is an indication that this combination is an efficient alternative to multi-temporal classification in the study region, where cloud free images are rare. The LULCC map obtained had an overall accuracy of 92%. Natural vegetation loss was estimated to be 17.9% ± 2.5% between 1999 and 2011. The study area experienced an increase in agricultural area and bare surface at the expense of woodland and mixed vegetation, which attests to the ongoing deforestation. These results can serve as means of regional and global land cover products validation, as they provide a new validated data set with uncertainty estimates in heterogeneous ecosystems prone to classification errors.
Schistosomiasis is a widespread water-based disease that puts close to 800 million people at risk of infection with more than 250 million infected, mainly in sub-Saharan Africa. Transmission is governed by the spatial distribution of specific freshwater snails that act as intermediate hosts and the frequency, duration and extent of human bodies exposed to infested water sources during human water contact. Remote sensing data have been utilized for spatially explicit risk profiling of schistosomiasis. Since schistosomiasis risk profiling based on remote sensing data inherits a conceptual drawback if school-based disease prevalence data are directly related to the remote sensing measurements extracted at the location of the school, because the disease transmission usually does not exactly occur at the school, we took the local environment around the schools into account by explicitly linking ecologically relevant environmental information of potential disease transmission sites to survey measurements of disease prevalence. Our models were validated at two sites with different landscapes in Côte d’Ivoire using high- and moderateresolution remote sensing data based on random forest and partial least squares regression. We found that the ecologically relevant modelling approach explained up to 70% of the variation in Schistosoma infection prevalence and performed better compared to a purely pixelbased modelling approach. Furthermore, our study showed that model performance increased as a function of enlarging the school catchment area, confirming the hypothesis that suitable environments for schistosomiasis transmission rarely occur at the location of survey measurements.
Remote sensing for disease risk profiling: a spatial analysis of schistosomiasis in West Africa
(2014)
Global environmental change leads to the emergence of new human health risks. As a consequence, transmission opportunities of environment-related diseases are transformed and human infection with new emerging pathogens increase. The main motivation for this study is the considerable demand for disease surveillance and monitoring in relation to dynamic environmental drivers. Remote sensing (RS) data belong to the key data sources for environmental modelling due to their capabilities to deliver spatially continuous information repeatedly for large areas with an ecologically adequate spatial resolution.
A major research gap as identified by this study is the disregard of the spatial mismatch inherent in current modelling approaches of profiling disease risk using remote sensing data. Typically, epidemiological data are aggregated at school or village level. However, these point data do neither represent the spatial distribution of habitats, where disease-related species find their suitable environmental conditions, nor the place, where infection has occurred. As a consequence, the prevalence data and remotely sensed environmental variables, which aim to characterise the habitat of disease-related species, are spatially disjunct.
The main objective of this study is to improve RS-based disease risk models by incorporating the ecological and spatial context of disease transmission. Exemplified by the analysis of the human schistosomiasis disease in West Africa, this objective includes the quantification of the impact of scales and ecological regions on model performance.
In this study, the conditions that modify the transmission of schistosomiasis are reviewed in detail. A conceptual underpinning of the linkages between geographical RS measures, disease transmission ecology, and epidemiological survey data is developed. During a field-based analysis, environmental suitability for schistosomiasis transmission was assessed on the ground, which is then quantified by a habitat suitability index (HSI) and applied to RS data. This conceptual model of environmental suitability is refined by the development of a hierarchical model approach that statistically links school-based disease prevalence with the ecologically relevant measurements of RS data. The statistical models of schistosomiasis risk are derived from two different algorithms; the Random Forest and the partial least squares regression (PLSR). Scale impact is analysed based on different spatial resolutions of RS data. Furthermore, varying buffer extents are analysed around school-based measurements. Three distinctive sites of Burkina Faso and Côte d’Ivoire are specifically modelled to represent a gradient of ecozones from dry savannah to tropical rainforest including flat and mountainous regions.
The model results reveal the applicability of RS data to spatially delineate and quantitatively evaluate environmental suitability for the transmission of schistosomiasis. In specific, the multi-temporal derivation of water bodies and the assessment of their riparian vegetation coverage based on high-resolution RapidEye and Landsat data proofed relevant. In contrast, elevation data and water surface temperature are constraint in their ability to characterise habitat conditions for disease-related parasites and freshwater snail species. With increasing buffer extent observed around the school location, the performance of statistical models increases, improving the prediction of transmission risk. The most important RS variables identified to model schistosomiasis risk are the measure of distance to water bodies, topographic variables, and land surface temperature (LST). However, each ecological region requires a different set of RS variables to optimise the modelling of schistosomiasis risk. A key result of the hierarchical model approach is its superior performance to explain the spatial risk of schistosomiasis.
Overall, this study stresses the key importance of considering the ecological and spatial context for disease risk profiling and demonstrates the potential of RS data. The methodological approach of this study contributes substantially to provide more accurate and relevant geoinformation, which supports an efficient planning and decision-making within the public health sector.
Agriculture is mankind’s primary source of food production and plays the key role for cereal supply to humanity. One of the future challenges will be to feed a constantly growing population, which is expected to reach more than nine billion by 2050. The potential to expand cropland is limited, and enhancing agricultural production efficiency is one important means to meet the future food demand. Hence, there is an increasing demand for dependable, accurate and comprehensive agricultural intelligence on crop production. The value of satellite earth observation (EO) data for agricultural monitoring is well recognized. One fundamental requirement for agricultural monitoring is routinely updated information on crop acreage and the spatial distribution of crops. With the technical advancement of satellite sensor systems, imagery with higher temporal and finer spatial resolution became available. The classification of such multi-temporal data sets is an effective and accurate means to produce crop maps, but methods must be developed that can handle such large and complex data sets. Furthermore, to properly use satellite EO for agricultural production monitoring a high temporal revisit frequency over vast geographic areas is often necessary. However, this often limits the spatial resolution that can be used. The challenge of discriminating pixels that correspond to a particular crop type, a prerequisite for crop specific agricultural monitoring, remains daunting when the signal encoded in pixels stems from several land uses (mixed pixels), e.g. over heterogeneous landscapes where individual fields are often smaller than individual pixels.
The main purposes of the presented study were (i) to assess the influence of input dimensionality and feature selection on classification accuracy and uncertainty in object-based crop classification, (ii) to evaluate if combining classifier algorithms can improve the quality of crop maps (e.g. classification accuracy), (iii) to assess the spatial resolution requirements for crop identification via image classification.
Reporting on the map quality is traditionally done with measures that stem from the confusion matrix based on the hard classification result. Yet, these measures do not consider the spatial variation of errors in maps. Measures of classification uncertainty can be used for this purpose, but they have attained only little attention in remote sensing studies. Classifier algorithms like the support vector machine (SVM) can estimate class memberships (the so called soft output) for each classified pixel or object. Based on these estimations, measures of classification uncertainty can be calculated, but it has not been analysed in detail, yet, if these are reliable in predicting the spatial distribution of errors in maps. In this study, SVM was applied for the classification of agricultural crops in irrigated landscapes in Middle Asia at the object-level. Five different categories of features were calculated from RapidEye time series data as classification input. The reliability of classification uncertainty measures like entropy, derived from the soft output of SVM, with regard to predicting the spatial distribution of error was evaluated. Further, the impact of the type and dimensionality of the input data on classification uncertainty was analysed. The results revealed that SMVs applied to the five feature categories separately performed different in classifying different types of crops. Incorporating all five categories of features by concatenating them into one stacked vector did not lead to an increase in accuracy, and partly reduced the model performance most obviously because of the Hughes phenomena. Yet, applying the random forest (RF) algorithm to select a subset of features led to an increase of classification accuracy of the SVM. The feature group with red edge-based indices was the most important for general crop classification, and the red edge NDVI had an outstanding importance for classifying crops. Two measures of uncertainty were calculated based on the soft output from SVM: maximum a-posteriori probability and alpha quadratic entropy. Irrespective of the measure used, the results indicate a decline in classification uncertainty when a dimensionality reduction was performed. The two uncertainty measures were found to be reliable indicators to predict errors in maps. Correctly classified test cases were associated with low uncertainty, whilst incorrectly test cases tended to be associated with higher uncertainty.
The issue of combining the results of different classifier algorithms in order to increase classification accuracy was addressed. First, the SVM was compared with two other non-parametric classifier algorithms: multilayer perceptron neural network (MLP) and RF. Despite their comparatively high classification performance, each of the tested classifier algorithms tended to make errors in different parts of the input space, e.g. performed different in classifying crops. Hence, a combination of the complementary outputs was envisaged. To this end, a classifier combination scheme was proposed, which is based on existing algebraic operators. It combines the outputs of different classifier algorithms at the per-case (e.g. pixel or object) basis. The per-case class membership estimations of each classifier algorithm were compared, and the reliability of each classifier algorithm with respect to classifying a specific crop class was assessed based on the confusion matrix. In doing so, less reliable classifier algorithms were excluded at the per-class basis before the final combination. Emphasis was put on evaluating the selected classification algorithms under limiting conditions by applying them to small input datasets and to reduced training sample sets, respectively. Further, the applicability to datasets from another year was demonstrated to assess temporal transferability. Although the single classifier algorithms performed well in all test sites, the classifier combination scheme provided consistently higher classification accuracies over all test sites and in different years, respectively. This makes this approach distinct from the single classifier algorithms, which performed different and showed a higher variability in class-wise accuracies. Further, the proposed classifier combination scheme performed better when using small training set sizes or when applied to small input datasets, respectively.
A framework was proposed to quantitatively define pixel size requirements for crop identification via image classification. That framework is based on simulating how agricultural landscapes, and more specifically the fields covered by one crop of interest, are seen by instruments with increasingly coarser resolving power. The concept of crop specific pixel purity, defined as the degree of homogeneity of the signal encoded in a pixel with respect to the target crop type, is used to analyse how mixed the pixels can be (as they become coarser) without undermining their capacity to describe the desired surface properties (e.g. to distinguish crop classes via supervised or unsupervised image classification). This tool can be modulated using different parameterizations to explore trade-offs between pixel size and pixel purity when addressing the question of crop identification. Inputs to the experiments were eight multi-temporal images from the RapidEye sensor. Simulated pixel sizes ranged from 13 m to 747.5 m, in increments of 6.5 m. Constraining parameters for crop identification were defined by setting thresholds for classification accuracy and uncertainty. Results over irrigated agricultural landscapes in Middle Asia demonstrate that the task of finding the optimum pixel size did not have a “one-size-fits-all” solution. The resulting values for pixel size and purity that were suitable for crop identification proved to be specific to a given landscape, and for each crop they differed across different landscapes. Over the same time series, different crops were not identifiable simultaneously in the season and these requirements further changed over the years, reflecting the different agro-ecological conditions the investigated crops were growing in. Results further indicate that map quality (e.g. classification accuracy) was not homogeneously distributed in a landscape, but that it depended on the spatial structures and the pixel size, respectively. The proposed framework is generic and can be applied to any agricultural landscape, thereby potentially serving to guide recommendations for designing dedicated EO missions that can satisfy the requirements in terms of pixel size to identify and discriminate crop types.
Regarding the operationalization of EO-based techniques for agricultural monitoring and its application to a broader range of agricultural landscapes, it can be noted that, despite the high performance of existing methods (e.g. classifier algorithms), transferability and stability of such methods remain one important research issue. This means that methods developed and tested in one place might not necessarily be portable to another place or over several years, respectively. Specifically in Middle Asia, which was selected as study region in this thesis, classifier combination makes sense due to its easy implementation and because it enhanced classification accuracy for classes with insufficient training samples. This observation makes it interesting for operational contexts and when field reference data availability is limited. Similar to the transferability of methods, the application of only one certain kind of EO data (e.g. with one specific pixel size) over different landscapes needs to be revisited and the synergistic use of multi-scale data, e.g. combining remote sensing imagery of both fine and coarse spatial resolution, should be fostered. The necessity to predict and control the effects of spatial and temporal scale on crop classification is recognized here as a major goal to achieve in EO-based agricultural monitoring.