TY - JOUR A1 - Dhillon, Maninder Singh A1 - Dahms, Thorsten A1 - Kuebert-Flock, Carina A1 - Rummler, Thomas A1 - Arnault, Joel A1 - Steffan-Dewenter, Ingolf A1 - Ullmann, Tobias T1 - Integrating random forest and crop modeling improves the crop yield prediction of winter wheat and oil seed rape JF - Frontiers in Remote Sensing N2 - 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. KW - crop modeling KW - random forest KW - machine learning KW - NDVI KW - satellite KW - landsat KW - sentinel-2 KW - winter wheat Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-301462 SN - 2673-6187 VL - 3 ER - TY - JOUR A1 - Kanmegne Tamga, Dan A1 - Latifi, Hooman A1 - Ullmann, Tobias A1 - Baumhauer, Roland A1 - Thiel, Michael A1 - Bayala, Jules T1 - Modelling the spatial distribution of the classification error of remote sensing data in cocoa agroforestry systems JF - Agroforestry Systems N2 - 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. KW - cocoa mapping KW - geographically weighted regression KW - Sentinel-1 KW - Sentinel-2 KW - Shannon entropy KW - spatial error assessment Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-324139 SN - 0167-4366 VL - 97 IS - 1 ER - TY - JOUR A1 - Philipp, Marius A1 - Dietz, Andreas A1 - Ullmann, Tobias A1 - Kuenzer, Claudia T1 - A circum-Arctic monitoring framework for quantifying annual erosion rates of permafrost coasts JF - Remote Sensing N2 - 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. KW - permafrost KW - coastal erosion KW - circum-Arctic KW - deep learning KW - change vector analysis KW - Google Earth Engine KW - synthetic aperture RADAR Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-304447 SN - 2072-4292 VL - 15 IS - 3 ER - TY - JOUR A1 - Dhillon, Maninder Singh A1 - Kübert-Flock, Carina A1 - Dahms, Thorsten A1 - Rummler, Thomas A1 - Arnault, Joel A1 - Steffan-Dewenter, Ingolf A1 - Ullmann, Tobias T1 - Evaluation of MODIS, Landsat 8 and Sentinel-2 data for accurate crop yield predictions: a case study using STARFM NDVI in Bavaria, Germany JF - Remote Sensing N2 - 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. KW - MODIS KW - Sentinel-2 KW - Landsat 8 KW - sustainable agriculture KW - decision-making KW - winter wheat KW - oil seed rape KW - resolution Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-311132 SN - 2072-4292 VL - 15 IS - 7 ER - TY - JOUR A1 - Dhillon, Maninder Singh A1 - Dahms, Thorsten A1 - Kübert-Flock, Carina A1 - Liepa, Adomas A1 - Rummler, Thomas A1 - Arnault, Joel A1 - Steffan-Dewenter, Ingolf A1 - Ullmann, Tobias T1 - Impact of STARFM on crop yield predictions: fusing MODIS with Landsat 5, 7, and 8 NDVIs in Bavaria Germany JF - Remote Sensing N2 - 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. KW - MOD13Q1 KW - precision agriculture KW - fusion KW - sustainable agriculture KW - decision making KW - winter wheat KW - oil-seed rape KW - crop models Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-311092 SN - 2072-4292 VL - 15 IS - 6 ER - TY - JOUR A1 - Reinermann, Sophie A1 - Asam, Sarah A1 - Gessner, Ursula A1 - Ullmann, Tobias A1 - Kuenzer, Claudia T1 - Multi-annual grassland mowing dynamics in Germany BT - spatio-temporal patterns and the influence of climate, topographic and socio-political conditions JF - Frontiers in Environmental Science N2 - 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. KW - remote sensing KW - Sentinel-2 KW - time series KW - cutting KW - management KW - pasture KW - meadow KW - Earth observation Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-320700 SN - 2296-665X VL - 11 ER - TY - JOUR A1 - Schäfer, Christian A1 - Fäth, Julian A1 - Kneisel, Christof A1 - Baumhauer, Roland A1 - Ullmann, Tobias T1 - Multidimensional hydrological modeling of a forested catchment in a German low mountain range using a modular runoff and water balance model JF - Frontiers in Forests and Global Change N2 - 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. KW - forest ecology KW - forest hydrology KW - WaSiM-ETH KW - drought stress indicators KW - beech Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-357358 VL - 6 ER - TY - JOUR A1 - Buchelt, Sebastian A1 - Skov, Kirstine A1 - Rasmussen, Kerstin Krøier A1 - Ullmann, Tobias T1 - Sentinel-1 time series for mapping snow cover depletion and timing of snowmelt in Arctic periglacial environments: case study from Zackenberg and Kobbefjord, Greenland JF - The Cryosphere N2 - 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. KW - Greenland KW - Sentinel-1 (S-1) synthetic aperture radar (SAR) KW - snow cover depletion Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-300139 VL - 16 IS - 2 SP - 625 EP - 646 ER - TY - JOUR A1 - Philipp, Marius A1 - Dietz, Andreas A1 - Ullmann, Tobias A1 - Kuenzer, Claudia T1 - Automated extraction of annual erosion rates for Arctic permafrost coasts using Sentinel-1, Deep Learning, and Change Vector Analysis JF - Remote Sensing N2 - 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. KW - permafrost KW - coastal erosion KW - deep learning KW - change vector analysis KW - Google Earth Engine KW - synthetic aperture RADAR Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-281956 SN - 2072-4292 VL - 14 IS - 15 ER - TY - JOUR A1 - Rösch, Moritz A1 - Sonnenschein, Ruth A1 - Buchelt, Sebastian A1 - Ullmann, Tobias T1 - Comparing PlanetScope and Sentinel-2 imagery for mapping mountain pines in the Sarntal Alps, Italy JF - Remote Sensing N2 - 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. KW - mountain pines KW - PlanetScope KW - Sentinel-2 KW - gray level co-occurrence matrix Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-281945 SN - 2072-4292 VL - 14 IS - 13 ER -