@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{NillUllmannKneiseletal.2019,
  author    = {Nill, Leon and Ullmann, Tobias and Kneisel, Christof and Sobiech-Wolf, Jennifer and Baumhauer, Roland},
  title     = {Assessing Spatiotemporal Variations of Landsat Land Surface Temperature and Multispectral Indices in the Arctic Mackenzie Delta Region between 1985 and 2018},
  series = {Remote Sensing},
  volume    = {11},
  journal   = {Remote Sensing},
  number    = {19},
  issn      = {2072-4292},
  doi       = {10.3390/rs11192329},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-193301},
  year      = {2019},
  abstract  = {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.},
  language  = {en}
}
@article{UllmannSauerbreyHoffmeisteretal.2019,
  author    = {Ullmann, Tobias and Sauerbrey, Julia and Hoffmeister, Dirk and May, Simon Matthias and Baumhauer, Roland and Bubenzer, Olaf},
  title     = {Assessing Spatiotemporal Variations of Sentinel-1 InSAR Coherence at Different Time Scales over the Atacama Desert (Chile) between 2015 and 2018},
  series = {Remote Sensing},
  volume    = {11},
  journal   = {Remote Sensing},
  number    = {24},
  issn      = {2072-4292},
  doi       = {10.3390/rs11242960},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-193836},
  pages     = {2960},
  year      = {2019},
  abstract  = {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.},
  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{DhillonDahmsKuebertFlocketal.2023,
  author    = {Dhillon, Maninder Singh and Dahms, Thorsten and Kuebert-Flock, Carina and Rummler, Thomas and Arnault, Joel and Steffan-Dewenter, Ingolf and Ullmann, Tobias},
  title     = {Integrating random forest and crop modeling improves the crop yield prediction of winter wheat and oil seed rape},
  series = {Frontiers in Remote Sensing},
  volume    = {3},
  journal   = {Frontiers in Remote Sensing},
  issn      = {2673-6187},
  doi       = {10.3389/frsen.2022.1010978},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-301462},
  year      = {2023},
  abstract  = {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.},
  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{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{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}
}