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 km²), 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² 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² [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.

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²) 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² = 0.82 and 0.85, RMSE = 5.46 and 5.01 dt/ha for WW, R² = 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² = 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² = 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.

Climate and land‐use change are key drivers of environmental degradation in the Anthropocene, but too little is known about their interactive effects on biodiversity and ecosystem services. Long‐term data on biodiversity trends are currently lacking. Furthermore, previous ecological studies have rarely considered climate and land use in a joint design, did not achieve variable independence or lost statistical power by not covering the full range of environmental gradients. Here, we introduce a multi‐scale space‐for‐time study design to disentangle effects of climate and land use on biodiversity and ecosystem services. The site selection approach coupled extensive GIS‐based exploration (i.e. using a Geographic information system) and correlation heatmaps with a crossed and nested design covering regional, landscape and local scales. Its implementation in Bavaria (Germany) resulted in a set of study plots that maximise the potential range and independence of environmental variables at different spatial scales. Stratifying the state of Bavaria into five climate zones (reference period 1981–2010) and three prevailing land‐use types, that is, near‐natural, agriculture and urban, resulted in 60 study regions (5.8 × 5.8 km quadrants) covering a mean annual temperature gradient of 5.6–9.8°C and a spatial extent of ~310 × 310 km. Within these regions, we nested 180 study plots located in contrasting local land‐use types, that is, forests, grasslands, arable land or settlement (local climate gradient 4.5–10°C). This approach achieved low correlations between climate and land use (proportional cover) at the regional and landscape scale with |r ≤ 0.33| and |r ≤ 0.29| respectively. Furthermore, using correlation heatmaps for local plot selection reduced potentially confounding relationships between landscape composition and configuration for plots located in forests, arable land and settlements. The suggested design expands upon previous research in covering a significant range of environmental gradients and including a diversity of dominant land‐use types at different scales within different climatic contexts. It allows independent assessment of the relative contribution of multi‐scale climate and land use on biodiversity and ecosystem services. Understanding potential interdependencies among global change drivers is essential to develop effective restoration and mitigation strategies against biodiversity decline, especially in expectation of future climatic changes. Importantly, this study also provides a baseline for long‐term ecological monitoring programs.

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² = 0.71, and RMSE = 0.11) synthetic NDVI products provide more accurate results than Landsat (overall R² = 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² = 0.62, RMSE = 0.11), resulting in similar accuracy preciseness as S-MOD13Q1 (R² = 0.68, RMSE = 0.13). Similarly, comparing L-MOD13Q1 (R² = 0.60, RMSE = 0.05) and S-MOD13Q1 (R² = 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.

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 R² (>0.82), low RMSE (<600 g/m²) and significant p-value (<0.05) during the study period. However, inputting MODIS data makes the models underperform, with low R² (600 g/m²). 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).

During the past decade, the carbon emission patterns due to groundwater pumping in central Punjab have changed enormously, and today the sector is one of the main contributors of C-emissions. This study identifies the critical blocks of central Punjab that should be aimed for reducing the carbon emissions of the groundwater economy. This study is based on the regional scale and highlights various factors that play crucial role in contributing to these emissions. Overlay analysis of the choropleth maps (tube well density, groundwater draft, groundwater depth and C-emissions) showed that SC-III zone was the most critical zone of central Punjab. Blocks Bhawanigarh, Dharmkot, Malerkotla I, Malerkotla II, Sangrur, Sherpur and Sunam were considered as threat to the environment and should be targeted to reduce C-emissions.