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.

Livestock plays an important economic role in Niger, especially in the semi-arid regions, while being highly vulnerable as a result of the large inter-annual variability of precipitation and, hence, rangeland production. This study aims to support effective rangeland management by developing an approach for mapping rangeland biomass production. The observed spatiotemporal variability of biomass production is utilised to build a model based on ground and remote sensing data for the period 2001 to 2015. Once established, the model can also be used to estimate herbaceous biomass for the current year at the end of the season without the need for new ground data. The phenology-based seasonal cumulative Normalised Difference Vegetation Index (cNDVI), computed from 10-day image composites of the Moderate-resolution Imaging Spectroradiometer (MODIS) NDVI data, was used as proxy for biomass production. A linear regression model was fitted with multi-annual field measurements of herbaceous biomass at the end of the growing season. In addition to a general model utilising all available sites for calibration, different aggregation schemes (i.e., grouping of sites into calibration units) of the study area with a varying number of calibration units and different biophysical meaning were tested. The sampling sites belonging to a specific calibration unit of a selected scheme were aggregated to compute the regression. The different aggregation schemes were evaluated with respect to their predictive power. The results gathered at the different aggregation levels were subjected to cross-validation (cv), applying a jackknife technique (leaving out one year at a time). In general, the model performance increased with increasing model parameterization, indicating the importance of additional unobserved and spatially heterogeneous agro-ecological effects (which might relate to grazing, species composition, optical soil properties, etc.) in modifying the relationship between cNDVI and herbaceous biomass at the end of the season. The biophysical aggregation scheme, the calibration units for which were derived from an unsupervised ISODATA classification utilising 10-day NDVI images taken between January 2001 and December 2015, showed the best performance in respect to the predictive power (R²_cv = 0.47) and the cross-validated root-mean-square error (398 kg·ha⁻¹) values, although it was not the model with the highest number of calibration units. The proposed approach can be applied for the timely production of maps of estimated biomass at the end of the growing season before field measurements are made available. These maps can be used for the improved management of rangeland resources, for decisions on fire prevention and aid allocation, and for the planning of more in-depth field missions.

Grasslands are an important part of pre-Alpine and Alpine landscapes. Despite the economic value and the significant role of grasslands in carbon and nitrogen (N) cycling, spatially explicit information on grassland biomass and quality is rarely available. Remotely sensed data from unmanned aircraft systems (UASs) and satellites might be an option to overcome this gap. Our study aims to investigate the potential of low-cost UAS-based multispectral sensors for estimating above-ground biomass (dry matter, DM) and plant N concentration. In our analysis, we compared two different sensors (Parrot Sequoia, SEQ; MicaSense RedEdge-M, REM), three statistical models (linear model; random forests, RFs; gradient-boosting machines, GBMs), and six predictor sets (i.e. different combinations of raw reflectance, vegetation indices, and canopy height). Canopy height information can be derived from UAS sensors but was not available in our study. Therefore, we tested the added value of this structural information with in situ measured bulk canopy height data. A combined field sampling and flight campaign was conducted in April 2018 at different grassland sites in southern Germany to obtain in situ and the corresponding spectral data. The hyper-parameters of the two machine learning (ML) approaches (RF, GBM) were optimized, and all model setups were run with a 6-fold cross-validation. Linear models were characterized by very low statistical performance measures, thus were not suitable to estimate DM and plant N concentration using UAS data. The non-linear ML algorithms showed an acceptable regression performance for all sensor–predictor set combinations with average (avg; cross-validated, cv) $R_{\text{cv}}^{\mathrm{2}}$ of 0.48, RMSE_cv,avg of 53.0 g m², and rRMSE_cv,avg (relative) of 15.9 % for DM and with $R_{\text{cv,avg}}^{\mathrm{2}}$ of 0.40, RMSE_cv,avg of 0.48 wt %, and rRMSE_{cv, avg} of 15.2 % for plant N concentration estimation. The optimal combination of sensors, ML algorithms, and predictor sets notably improved the model performance. The best model performance for the estimation of DM ($R_{\text{cv}}^{\mathrm{2}}=\mathrm{0.67}$ , RMSE_cv=41.9 g m², rRMSE_cv=12.6 %) was achieved with an RF model that utilizes all possible predictors and REM sensor data. The best model for plant N concentration was a combination of an RF model with all predictors and SEQ sensor data ($R_{\text{cv}}^{\mathrm{2}}=\mathrm{0.47}$ , RMSE_cv=0.45 wt %, rRMSE_cv=14.2 %). DM models with the spectral input of REM performed significantly better than those with SEQ data, while for N concentration models, it was the other way round. The choice of predictors was most influential on model performance, while the effect of the chosen ML algorithm was generally lower. The addition of canopy height to the spectral data in the predictor set significantly improved the DM models. In our study, calibrating the ML algorithm improved the model performance substantially, which shows the importance of this step.

The data set contains information on aboveground vegetation traits of > 100 georeferenced locations within ten temperate pre-Alpine grassland plots in southern Germany. The grasslands were sampled in April 2018 for the following traits: bulk canopy height; weight of fresh and dry biomass; dry weight percentage of the plant functional types (PFT) non-green vegetation, legumes, non-leguminous forbs, and graminoids; total green area index (GAI) and PFT-specific GAI; plant water content; plant carbon and nitrogen content (community values and PFT-specific values); as well as leaf mass per area (LMA) of PFT. In addition, a species specific inventory of the plots was conducted in June 2020 and provides plot-level information on grassland type and plant species composition. The data set was obtained within the framework of the SUSALPS project (“Sustainable use of alpine and pre-alpine grassland soils in a changing climate”;https://www.susalps.de/) to providein-situdata for the calibration and validation of remote sensing based models to estimate grassland traits.

Semi-natural, agriculturally used grasslands provide important ecologic and economic services, such as feed supply. In mountain regions, pastures are the dominant agricultural system and face more severe climate change impacts than lowlands. Climate change threatens ecosystem functions, such as aboveground net primary production [ANPP] and its nutrient content. It is necessary to understand the impacts of climate change and land-management on such ecosystems to develop management practices to sustainably maintain provision of ecosystem services under future climatic conditions. We studied the effect of climate change and different land-use intensities on plant-soil communities by the downslope translocation of plant-soil mesocosms along an elevation gradient in 2016, and the subsequent application of two management types (extensive vs. intensive). Communities’ response to ANPP and leaf carbon (C), nitrogen (N), and phosphorus (P) content was quantified over the subsequent two years after translocation. ANPP increased with warming in 2017 under both management intensities, but this effect was amplified by intensive land-use management. In 2018, ANPP of intensively managed communities decreased, in comparison to 2017, from 35% to 42%, while extensively managed communities maintained their production levels. The changes in ANPP are coupled with an exceptionally dry year in 2018, with up to 100 more days of drought conditions. The C:N of extensively managed communities was higher than those of intensively managed ones, and further increased in 2018, potentially indicating shifts in resource allocation strategies that may explain production stability. Our results revealed a low resistance of intensively managed communities’ ANPP under especially dry conditions. The ability to alter resource allocation likely enables a constant level of production under extensive management, but this ability is lost under intensive management. Thus, future drought events may leave intensive management as a non-sustainable farming practice, and ultimately threaten ecosystem services of montane pastures.

Grasslands are an important part of pre-Alpine and Alpine landscapes. Despite the economic value and the significant role of grasslands in carbon and nitrogen (N) cycling, spatially explicit information on grassland biomass and quality is rarely available. Remotely sensed data from unmanned aircraft systems (UAS) and satellites might be an option to overcome this gap. Our study aims to investigate the potential of low-cost UAS-based multispectral sensors for estimating above-ground biomass (dry matter, DM) and plant N concentration. In our analysis, we compared two different sensors (Parrot Sequoia, SEQ; MicaSense RedEdge-M, REM), three statistical models (Linear Model; Random Forests, RF; Gradient Boosting Machines, GBM) and six predictor sets (i.e. different combinations of raw reflectance, vegetation indices, and canopy height). Canopy height information can be derived from UAS sensors, but was not available in our study. Therefore, we tested the added value of this structural information with in-situ measured bulk canopy height data. A combined field sampling and flight campaign was conducted in April 2018 at different grassland sites in Southern Germany to obtain in-situ and the corresponding spectral data. The hyper-parameters of the two machine learning (ML) approaches (RF, GBM) were optimized and all model set-ups were run with a six-fold cross-validation. Linear models were characterized by very low statistical performance measures, thus were not suitable to estimate DM and plant N concentration using UAS data. The non-linear ML algorithms showed an acceptable regression performance for all sensor-predictor set combinations with average (avg) R²_cv of 0.48, RMSE_{cv, avg} of 53.0 g m² and rRMSE_{cv, avg} of 15.9 % for DM, and with R²_{cv, avg} of 0.40, RMSE_{cv, avg} of 0.48 wt.% and rRMSE_{cv, avg} of 15.2 % for plant N concentration estimation. The optimal combination of sensors, ML algorithms and predictor sets notably improved the model performance. The best model performance for the estimation of DM (R²_cv = 0.67, RMSE_cv = 41.9 g m², rRMSE_cv = 12.6 %) was achieved with a RF model that utilizes all possible predictors and REM sensor data. The best model for plant N concentration was a combination of a RF model with all predictors and SEQ sensor data (R²_cv = 0.47, RMSE_cv = 0.45 wt.%, rRMSE_cv = 14.2 %). DM models with the spectral input of REM performed significantly better than those with SEQ data, while for N concentration models it was the other way round. The choice of predictors was most influential on model performance, while the effect of the chosen ML algorithm was generally lower. The addition of canopy height to the spectral data in the predictor set significantly improved the DM models. In our study, calibrating ML algorithm improved the model performance substantially, which shows the importance of this step.