Refine
Has Fulltext
- yes (32)
Is part of the Bibliography
- yes (32)
Year of publication
Document Type
- Journal article (31)
- Doctoral Thesis (1)
Language
- English (32)
Keywords
- Sentinel-2 (5)
- TerraSAR-X (5)
- Google Earth Engine (4)
- PolSAR (4)
- Sentinel-1 (4)
- winter wheat (4)
- InSAR (3)
- Radarsat-2 (3)
- SAR (3)
- arctic (3)
Institute
EU-Project number / Contract (GA) number
- 20-3044-2-11 (1)
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.
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 ecosystem of the high northern latitudes is affected by the recently changing environmental conditions. The Arctic has undergone a significant climatic change over the last decades. The land coverage is changing and a phenological response to the warming is apparent. Remotely sensed data can assist the monitoring and quantification of these changes. The remote sensing of the Arctic was predominantly carried out by the usage of optical sensors but these encounter problems in the Arctic environment, e.g. the frequent cloud cover or the solar geometry. In contrast, the imaging of Synthetic Aperture Radar is not affected by the cloud cover and the acquisition of radar imagery is independent of the solar illumination. The objective of this work was to explore how polarimetric Synthetic Aperture Radar (PolSAR) data of TerraSAR-X, TanDEM-X, Radarsat-2 and ALOS PALSAR and interferometric-derived digital elevation model data of the TanDEM-X Mission can contribute to collect meaningful information on the actual state of the Arctic Environment. The study was conducted for Canadian sites of the Mackenzie Delta Region and Banks Island and in situ reference data were available for the assessment. The up-to-date analysis of the PolSAR data made the application of the Non-Local Means filtering and of the decomposition of co-polarized data necessary.
The Non-Local Means filter showed a high capability to preserve the image values, to keep the edges and to reduce the speckle. This supported not only the suitability for the interpretation but also for the classification. The classification accuracies of Non-Local Means filtered data were in average +10% higher compared to unfiltered images. The correlation of the co- and quad-polarized decomposition features was high for classes with distinct surface or double bounce scattering and a usage of the co-polarized data is beneficial for regions of natural land coverage and for low vegetation formations with little volume scattering. The evaluation further revealed that the X- and C-Band were most sensitive to the generalized land cover classes. It was found that the X-Band data were sensitive to low vegetation formations with low shrub density, the C-Band data were sensitive to the shrub density and the shrub dominated tundra. In contrast, the L-Band data were less sensitive to the land cover. Among the different dual-polarized data the HH/VV-polarized data were identified to be most meaningful for the characterization and classification, followed by the HH/HV-polarized and the VV/VH-polarized data. The quad-polarized data showed highest sensitivity to the land cover but differences to the co-polarized data were small. The accuracy assessment showed that spectral information was required for accurate land cover classification. The best results were obtained when spectral and radar information was combined. The benefit of including radar data in the classification was up to +15% accuracy and most significant for the classes wetland and sparse vegetated tundra. The best classifications were realized with quad-polarized C-Band and multispectral data and with co-polarized X-Band and multispectral data. The overall accuracy was up to 80% for unsupervised and up to 90% for supervised classifications. The results indicated that the shortwave co-polarized data show promise for the classification of tundra land cover since the polarimetric information is sensitive to low vegetation and the wetlands. Furthermore, co-polarized data provide a higher spatial resolution than the quad-polarized data.
The analysis of the intermediate digital elevation model data of the TanDEM-X showed a high potential for the characterization of the surface morphology. The basic and relative topographic features were shown to be of high relevance for the quantification of the surface morphology and an area-wide application is feasible. In addition, these data were of value for the classification and delineation of landforms. Such classifications will assist the delineation of geomorphological units and have potential to identify locations of actual and future morphologic activity.
In this work the potential of polarimetric Synthetic Aperture Radar (PolSAR) data of dual-polarized TerraSAR-X (HH/VV) and quad-polarized Radarsat-2 was examined in combination with multispectral Landsat 8 data for unsupervised and supervised classification of tundra land cover types of Richards Island, Canada. The classification accuracies as well as the backscatter and reflectance characteristics were analyzed using reference data collected during three field work campaigns and include in situ data and high resolution airborne photography. The optical data offered an acceptable initial accuracy for the land cover classification. The overall accuracy was increased by the combination of PolSAR and optical data and was up to 71% for unsupervised (Landsat 8 and TerraSAR-X) and up to 87% for supervised classification (Landsat 8 and Radarsat-2) for five tundra land cover types. The decomposition features of the dual and quad-polarized data showed a high sensitivity for the non-vegetated substrate (dominant surface scattering) and wetland vegetation (dominant double bounce and volume scattering). These classes had high potential to be automatically detected with unsupervised classification techniques.
Crop mapping in West Africa is challenging, due to the unavailability of adequate satellite images (as a result of excessive cloud cover), small agricultural fields and a heterogeneous landscape. To address this challenge, we integrated high spatial resolution multi-temporal optical (RapidEye) and dual polarized (VV/VH) SAR (TerraSAR-X) data to map crops and crop groups in northwestern Benin using the random forest classification algorithm. The overall goal was to ascertain the contribution of the SAR data to crop mapping in the region. A per-pixel classification result was overlaid with vector field boundaries derived from image segmentation, and a crop type was determined for each field based on the modal class within the field. A per-field accuracy assessment was conducted by comparing the final classification result with reference data derived from a field campaign. Results indicate that the integration of RapidEye and TerraSAR-X data improved classification accuracy by 10%–15% over the use of RapidEye only. The VV polarization was found to better discriminate crop types than the VH polarization. The research has shown that if optical and SAR data are available for the whole cropping season, classification accuracies of up to 75% are achievable.
Snow cover (SC) and timing of snowmelt are key regulators of a wide range of Arctic ecosystem functions. Both are strongly influenced by the amplified Arctic warming and essential variables to understand environmental changes and their dynamics. This study evaluates the potential of Sentinel-1 (S-1) synthetic aperture radar (SAR) time series for monitoring SC depletion and snowmelt with high spatiotemporal resolution to capture their understudied small-scale heterogeneity. We use 97 dual-polarized S-1 SAR images acquired over northeastern Greenland and 94 over southwestern Greenland in the interferometric wide swath mode from the years 2017 and 2018. Comparison of S-1 intensity against SC fraction maps derived from orthorectified terrestrial time-lapse imagery indicates that SAR backscatter can increase before a decrease in SC fraction is observed. Hence, the increase in backscatter is related to changing snowpack properties during the runoff phase as well as decreasing SC fraction. We here present a novel empirical approach based on the temporal evolution of the SAR signal to identify start of runoff (SOR), end of snow cover (EOS) and SC extent for each S-1 observation date during melt using backscatter thresholds as well as the derivative. Comparison of SC with orthorectified time-lapse imagery indicates that HV polarization outperforms HH when using a global threshold. The derivative avoids manual selection of thresholds and adapts to different environmental settings and seasonal conditions. With a global configuration (threshold: 4 dB; polarization: HV) as well as with the derivative, the overall accuracy of SC maps was in all cases above 75 % and in more than half of cases above 90 %. Based on the physical principle of SAR backscatter during snowmelt, our approach is expected to work well in other low-vegetation areas and, hence, could support large-scale SC monitoring at high spatiotemporal resolution (20 m, 6 d) with high accuracy.
Cocoa growing is one of the main activities in humid West Africa, which is mainly grown in pure stands. It is the main driver of deforestation and encroachment in protected areas. Cocoa agroforestry systems which have been promoted to mitigate deforestation, needs to be accurately delineated to support a valid monitoring system. Therefore, the aim of this research is to model the spatial distribution of uncertainties in the classification cocoa agroforestry. The study was carried out in Côte d’Ivoire, close to the Taï National Park. The analysis followed three steps (i) image classification based on texture parameters and vegetation indices from Sentinel-1 and -2 data respectively, to train a random forest algorithm. A classified map with the associated probability maps was generated. (ii) Shannon entropy was calculated from the probability maps, to get the error maps at different thresholds (0.2, 0.3, 0.4 and 0.5). Then, (iii) the generated error maps were analysed using a Geographically Weighted Regression model to check for spatial autocorrelation. From the results, a producer accuracy (0.88) and a user’s accuracy (0.91) were obtained. A small threshold value overestimates the classification error, while a larger threshold will underestimate it. The optimal value was found to be between 0.3 and 0.4. There was no evidence of spatial autocorrelation except for a smaller threshold (0.2). The approach differentiated cocoa from other landcover and detected encroachment in forest. Even though some information was lost in the process, the method is effective for mapping cocoa plantations in Côte d’Ivoire.
A circum-Arctic monitoring framework for quantifying annual erosion rates of permafrost coasts
(2023)
This study demonstrates a circum-Arctic monitoring framework for quantifying annual change of permafrost-affected coasts at a spatial resolution of 10 m. Frequent cloud coverage and challenging lighting conditions, including polar night, limit the usability of optical data in Arctic regions. For this reason, Synthetic Aperture RADAR (SAR) data in the form of annual median and standard deviation (sd) Sentinel-1 (S1) backscatter images covering the months June–September for the years 2017–2021 were computed. Annual composites for the year 2020 were hereby utilized as input for the generation of a high-quality coastline product via a Deep Learning (DL) workflow, covering 161,600 km of the Arctic coastline. The previously computed annual S1 composites for the years 2017 and 2021 were employed as input data for the Change Vector Analysis (CVA)-based coastal change investigation. The generated DL coastline product served hereby as a reference. Maximum erosion rates of up to 67 m per year could be observed based on 400 m coastline segments. Overall highest average annual erosion can be reported for the United States (Alaska) with 0.75 m per year, followed by Russia with 0.62 m per year. Out of all seas covered in this study, the Beaufort Sea featured the overall strongest average annual coastal erosion of 1.12 m. Several quality layers are provided for both the DL coastline product and the CVA-based coastal change analysis to assess the applicability and accuracy of the output products. The predicted coastal change rates show good agreement with findings published in previous literature. The proposed methods and data may act as a valuable tool for future analysis of permafrost loss and carbon emissions in Arctic coastal environments.
The increasing availability and variety of global satellite products and the rapid development of new algorithms has provided great potential to generate a new level of data with different spatial, temporal, and spectral resolutions. However, the ability of these synthetic spatiotemporal datasets to accurately map and monitor our planet on a field or regional scale remains underexplored. This study aimed to support future research efforts in estimating crop yields by identifying the optimal spatial (10 m, 30 m, or 250 m) and temporal (8 or 16 days) resolutions on a regional scale. The current study explored and discussed the suitability of four different synthetic (Landsat (L)-MOD13Q1 (30 m, 8 and 16 days) and Sentinel-2 (S)-MOD13Q1 (10 m, 8 and 16 days)) and two real (MOD13Q1 (250 m, 8 and 16 days)) NDVI products combined separately to two widely used crop growth models (CGMs) (World Food Studies (WOFOST), and the semi-empiric Light Use Efficiency approach (LUE)) for winter wheat (WW) and oil seed rape (OSR) yield forecasts in Bavaria (70,550 km\(^2\)) for the year 2019. For WW and OSR, the synthetic products’ high spatial and temporal resolution resulted in higher yield accuracies using LUE and WOFOST. The observations of high temporal resolution (8-day) products of both S-MOD13Q1 and L-MOD13Q1 played a significant role in accurately measuring the yield of WW and OSR. For example, L- and S-MOD13Q1 resulted in an R\(^2\) = 0.82 and 0.85, RMSE = 5.46 and 5.01 dt/ha for WW, R\(^2\) = 0.89 and 0.82, and RMSE = 2.23 and 2.11 dt/ha for OSR using the LUE model, respectively. Similarly, for the 8- and 16-day products, the simple LUE model (R\(^2\) = 0.77 and relative RMSE (RRMSE) = 8.17%) required fewer input parameters to simulate crop yield and was highly accurate, reliable, and more precise than the complex WOFOST model (R\(^2\) = 0.66 and RRMSE = 11.35%) with higher input parameters. Conclusively, both S-MOD13Q1 and L-MOD13Q1, in combination with LUE, were more prominent for predicting crop yields on a regional scale than the 16-day products; however, L-MOD13Q1 was advantageous for generating and exploring the long-term yield time series due to the availability of Landsat data since 1982, with a maximum resolution of 30 m. In addition, this study recommended the further use of its findings for implementing and validating the long-term crop yield time series in different regions of the world.
Mapping buried paleogeographical features of the Nile Delta (Egypt) using the Landsat archive
(2020)
The contribution highlights the use of Landsat spectral-temporal metrics (STMs) for the detection of surface anomalies that are potentially related to buried near-surface paleogeomorphological deposits in the Nile Delta (Egypt), in particular for a buried river branch close to Buto. The processing was completed in the Google Earth Engine (GEE) for the entire Nile Delta and for selected seasons of the year (summer/winter) using Landsat data from 1985 to 2019. We derived the STMs of the tasseled cap transformation (TC), the Normalized Difference Wetness Index (NDWI), and the Normalized Difference Vegetation Index (NDVI). These features were compared to historical topographic maps of the Survey of Egypt, CORONA imagery, the digital elevation model of the TanDEM-X mission, and modern high-resolution satellite imagery. The results suggest that the extent of channels is best revealed when differencing the median NDWI between summer (July/August) and winter (January/February) seasons (ΔNDWI). The observed difference is likely due to lower soil/plant moisture during summer, which is potentially caused by coarser-grained deposits and the morphology of the former levee. Similar anomalies were found in the immediate surroundings of several Pleistocene sand hills (“geziras”) and settlement mounds (“tells”) of the eastern delta, which allowed some mapping of the potential near-surface continuation. Such anomalies were not observed for the surroundings of tells of the western Nile Delta. Additional linear and meandering ΔNDWI anomalies were found in the eastern Nile Delta in the immediate surroundings of the ancient site of Bubastis (Tell Basta), as well as several kilometers north of Zagazig. These anomalies might indicate former courses of Nile river branches. However, the ΔNDWI does not provide an unambiguous delineation.
Geoarchaeological information presented here pertains to a subsidiary Nile channel that once flowed west of the main Sebennitic distributary and discharged its water and sediments at Egypt’s then north-central deltaic coast. Periodical paleoclimatic episodes during the later Middle and Upper Holocene included decreased rainfall and increased aridity that reduced the Nile’s flow levels and thus likely disrupted nautical transport and anthropogenic activity along this channel. Such changes in this deltaic sector, positioned adjacent to the Levantine Basin in the Eastern Mediterranean, can be attributed to climatic shifts triggered as far as the North Atlantic to the west, and African highland source areas of the Egyptian Nile to the south. Of special interest in a study core recovered along the channel are several sediment sequences without anthropogenic material that are interbedded between strata comprising numerous potsherds. The former are interpreted here as markers of increased regional aridity and reduced Nile flow which could have periodically disrupted the regional distribution of goods and nautical activities. Such times occurred ~5000 years B.P., ~4200–4000 years B.P., ~3200–2800 years B.P., ~2300–2200 years B.P., and more recently. Periods comparable to these are also identified by altered proportions of pollen, isotopic and compositional components in different radiocarbon-dated Holocene cores recovered elsewhere in the Nile delta, the Levantine region to the east and north of Egypt, and in the Faiyum depression south of the delta.
Sufficient plant-available water is one of the most important requirements for vital, stable, and well-growing forest stands. In the face of climate change, there are various approaches to derive recommendations considering tree species selection based on plant-available water provided by measurements or simulations. Owing to the small-parcel management of Central European forests as well as small-spatial variation of soil and stand properties, in situ data collection for individual forest stands of large areas is not feasible, considering time and cost effort. This problem can be addressed using physically based modeling, aiming to numerically simulate the water balance. In this study, we parameterized, calibrated, and verified the hydrological multidimensional WaSiM-ETH model to assess the water balance at a spatial resolution of 30 m in a German forested catchment area (136.4 km2) for the period 2000–2021 using selected in situ data, remote sensing products, and total runoff. Based on the model output, drought-sensitive parameters, such as the difference between potential and effective stand transpiration (Tdiff) and the water balance, were deduced from the model, analyzed, and evaluated. Results show that the modeled evapotranspiration (ET) correlated significantly (R2 = 0.80) with the estimated ET using MODIS data (MOD16A2GFv006). Compared with observed daily, monthly, and annual runoff data, the model shows a good performance (R2: 0.70|0.77|0.73; Kling–Gupta efficiency: 0.59|0.62|0.83; volumetric efficiency: 0.52|0.60|0.83). The comparison with in situ data from a forest monitoring plot, established at the end of 2020, indicated good agreement between observed and simulated interception and soil water content. According to our results, WaSiM-ETH is a potential supplement for forest management, owing to its multidimensionality and the ability to model soil water balance for large areas at comparable high spatial resolution. The outputs offer, compared to non-distributed models (like LWF-Brook90), spatial differentiability, which is important for small-scale parceled forests, regarding stand structure and soil properties. Due to the spatial component offered, additional verification possibilities are feasible allowing a reliable and profound verification of the model and its parameterization.
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.
Rapid and accurate yield estimates at both field and regional levels remain the goal of sustainable agriculture and food security. Hereby, the identification of consistent and reliable methodologies providing accurate yield predictions is one of the hot topics in agricultural research. This study investigated the relationship of spatiotemporal fusion modelling using STRAFM on crop yield prediction for winter wheat (WW) and oil-seed rape (OSR) using a semi-empirical light use efficiency (LUE) model for the Free State of Bavaria (70,550 km\(^2\)), Germany, from 2001 to 2019. A synthetic normalised difference vegetation index (NDVI) time series was generated and validated by fusing the high spatial resolution (30 m, 16 days) Landsat 5 Thematic Mapper (TM) (2001 to 2012), Landsat 7 Enhanced Thematic Mapper Plus (ETM+) (2012), and Landsat 8 Operational Land Imager (OLI) (2013 to 2019) with the coarse resolution of MOD13Q1 (250 m, 16 days) from 2001 to 2019. Except for some temporal periods (i.e., 2001, 2002, and 2012), the study obtained an R\(^2\) of more than 0.65 and a RMSE of less than 0.11, which proves that the Landsat 8 OLI fused products are of higher accuracy than the Landsat 5 TM products. Moreover, the accuracies of the NDVI fusion data have been found to correlate with the total number of available Landsat scenes every year (N), with a correlation coefficient (R) of +0.83 (between R\(^2\) of yearly synthetic NDVIs and N) and −0.84 (between RMSEs and N). For crop yield prediction, the synthetic NDVI time series and climate elements (such as minimum temperature, maximum temperature, relative humidity, evaporation, transpiration, and solar radiation) are inputted to the LUE model, resulting in an average R\(^2\) of 0.75 (WW) and 0.73 (OSR), and RMSEs of 4.33 dt/ha and 2.19 dt/ha. The yield prediction results prove the consistency and stability of the LUE model for yield estimation. Using the LUE model, accurate crop yield predictions were obtained for WW (R\(^2\) = 0.88) and OSR (R\(^2\) = 0.74). Lastly, the study observed a high positive correlation of R = 0.81 and R = 0.77 between the yearly R\(^2\) of synthetic accuracy and modelled yield accuracy for WW and OSR, respectively.
Introduction: Grasslands cover one third of the agricultural area in Germany and are mainly used for fodder production. However, grasslands fulfill many other ecosystem functions, like carbon storage, water filtration and the provision of habitats. In Germany, grasslands are mown and/or grazed multiple times during the year. The type and timing of management activities and the use intensity vary strongly, however co-determine grassland functions. Large-scale spatial information on grassland activities and use intensity in Germany is limited and not openly provided. In addition, the cause for patterns of varying mowing intensity are usually not known on a spatial scale as data on the incentives of farmers behind grassland management decisions is not available.
Methods: We applied an algorithm based on a thresholding approach utilizing Sentinel-2 time series to detect grassland mowing events to investigate mowing dynamics in Germany in 2018–2021. The detected mowing events were validated with an independent dataset based on the examination of public webcam images. We analyzed spatial and temporal patterns of the mowing dynamics and relationships to climatic, topographic, soil or socio-political conditions.
Results: We found that most intensively used grasslands can be found in southern/south-eastern Germany, followed by areas in northern Germany. This pattern stays the same among the investigated years, but we found variations on smaller scales. The mowing event detection shows higher accuracies in 2019 and 2020 (F1 = 0.64 and 0.63) compared to 2018 and 2021 (F1 = 0.52 and 0.50). We found a significant but weak (R2 of 0–0.13) relationship for a spatial correlation of mowing frequency and climate as well as topographic variables for the grassland areas in Germany. Further results indicate a clear value range of topographic and climatic conditions, characteristic for intensive grassland use. Extensive grassland use takes place everywhere in Germany and on the entire spectrum of topographic and climatic conditions in Germany. Natura 2000 grasslands are used less intensive but this pattern is not consistent among all sites.
Discussion: Our findings on mowing dynamics and relationships to abiotic and socio-political conditions in Germany reveal important aspects of grassland management, including incentives of farmers.
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).
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.
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.