@article{PhilippDietzUllmannetal.2023,
  author    = {Philipp, Marius and Dietz, Andreas and Ullmann, Tobias and Kuenzer, Claudia},
  title     = {A circum-Arctic monitoring framework for quantifying annual erosion rates of permafrost coasts},
  series = {Remote Sensing},
  volume    = {15},
  journal   = {Remote Sensing},
  number    = {3},
  issn      = {2072-4292},
  doi       = {10.3390/rs15030818},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-304447},
  year      = {2023},
  abstract  = {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.},
  language  = {en}
}
@article{DhillonKuebertFlockDahmsetal.2023,
  author    = {Dhillon, Maninder Singh and K{\"u}bert-Flock, Carina and Dahms, Thorsten and Rummler, Thomas and Arnault, Joel and Steffan-Dewenter, Ingolf and Ullmann, Tobias},
  title     = {Evaluation of MODIS, Landsat 8 and Sentinel-2 data for accurate crop yield predictions: a case study using STARFM NDVI in Bavaria, Germany},
  series = {Remote Sensing},
  volume    = {15},
  journal   = {Remote Sensing},
  number    = {7},
  issn      = {2072-4292},
  doi       = {10.3390/rs15071830},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-311132},
  year      = {2023},
  abstract  = {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.},
  language  = {en}
}
@article{DhillonDahmsKuebertFlocketal.2023,
  author    = {Dhillon, Maninder Singh and Dahms, Thorsten and K{\"u}bert-Flock, Carina and Liepa, Adomas and Rummler, Thomas and Arnault, Joel and Steffan-Dewenter, Ingolf and Ullmann, Tobias},
  title     = {Impact of STARFM on crop yield predictions: fusing MODIS with Landsat 5, 7, and 8 NDVIs in Bavaria Germany},
  series = {Remote Sensing},
  volume    = {15},
  journal   = {Remote Sensing},
  number    = {6},
  issn      = {2072-4292},
  doi       = {10.3390/rs15061651},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-311092},
  year      = {2023},
  abstract  = {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.},
  language  = {en}
}
@article{KanmegneTamgaLatifiUllmannetal.2023,
  author    = {Kanmegne Tamga, Dan and Latifi, Hooman and Ullmann, Tobias and Baumhauer, Roland and Thiel, Michael and Bayala, Jules},
  title     = {Modelling the spatial distribution of the classification error of remote sensing data in cocoa agroforestry systems},
  series = {Agroforestry Systems},
  volume    = {97},
  journal   = {Agroforestry Systems},
  number    = {1},
  issn      = {0167-4366},
  doi       = {10.1007/s10457-022-00791-2},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-324139},
  pages     = {109-119},
  year      = {2023},
  abstract  = {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{\^o}te d'Ivoire, close to the Ta{\"i} 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{\^o}te d'Ivoire.},
  language  = {en}
}
@article{DhillonDahmsKuebertFlocketal.2022,
  author    = {Dhillon, Maninder Singh and Dahms, Thorsten and K{\"u}bert-Flock, Carina and Steffan-Dewenter, Ingolf and Zhang, Jie and Ullmann, Tobias},
  title     = {Spatiotemporal Fusion Modelling Using STARFM: Examples of Landsat 8 and Sentinel-2 NDVI in Bavaria},
  series = {Remote Sensing},
  volume    = {14},
  journal   = {Remote Sensing},
  number    = {3},
  issn      = {2072-4292},
  doi       = {10.3390/rs14030677},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-323471},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}
@article{ReinermannGessnerAsametal.2022,
  author    = {Reinermann, Sophie and Gessner, Ursula and Asam, Sarah and Ullmann, Tobias and Schucknecht, Anne and Kuenzer, Claudia},
  title     = {Detection of grassland mowing events for Germany by combining Sentinel-1 and Sentinel-2 time series},
  series = {Remote Sensing},
  volume    = {14},
  journal   = {Remote Sensing},
  number    = {7},
  issn      = {2072-4292},
  doi       = {10.3390/rs14071647},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-267164},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}
@article{UllmannMoellerBaumhaueretal.2022,
  author    = {Ullmann, Tobias and M{\"o}ller, Eric and Baumhauer, Roland and Lange-Athinodorou, Eva and Meister, Julia},
  title     = {A new Google Earth Engine tool for spaceborne detection of buried palaeogeographical features - examples from the Nile Delta (Egypt)},
  series = {E\&G Quaternary Science Journal},
  volume    = {71},
  journal   = {E\&G Quaternary Science Journal},
  number    = {2},
  doi       = {10.5194/egqsj-71-243-2022},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-300164},
  pages     = {243-247},
  year      = {2022},
  abstract  = {No abstract available.},
  language  = {en}
}
@article{LappeUllmannBachofer2022,
  author    = {Lappe, Ronja and Ullmann, Tobias and Bachofer, Felix},
  title     = {State of the Vietnamese coast — assessing three decades (1986 to 2021) of coastline dynamics using the Landsat archive},
  series = {Remote Sensing},
  volume    = {14},
  journal   = {Remote Sensing},
  number    = {10},
  issn      = {2072-4292},
  doi       = {10.3390/rs14102476},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-275281},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}
@article{RoeschSonnenscheinBucheltetal.2022,
  author    = {R{\"o}sch, Moritz and Sonnenschein, Ruth and Buchelt, Sebastian and Ullmann, Tobias},
  title     = {Comparing PlanetScope and Sentinel-2 imagery for mapping mountain pines in the Sarntal Alps, Italy},
  series = {Remote Sensing},
  volume    = {14},
  journal   = {Remote Sensing},
  number    = {13},
  issn      = {2072-4292},
  doi       = {10.3390/rs14133190},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-281945},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}
@article{PhilippDietzUllmannetal.2022,
  author    = {Philipp, Marius and Dietz, Andreas and Ullmann, Tobias and Kuenzer, Claudia},
  title     = {Automated extraction of annual erosion rates for Arctic permafrost coasts using Sentinel-1, Deep Learning, and Change Vector Analysis},
  series = {Remote Sensing},
  volume    = {14},
  journal   = {Remote Sensing},
  number    = {15},
  issn      = {2072-4292},
  doi       = {10.3390/rs14153656},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-281956},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}