Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Vet. Sci., 19 December 2025

Sec. Animal Behavior and Welfare

Volume 12 - 2025 | https://doi.org/10.3389/fvets.2025.1678123

This article is part of the Research TopicAdvances in Precision Livestock Management for Grazing Ruminant SystemsView all 11 articles

Satellite-based monitoring of forage quality in grasslands of the United Kingdom using sentinel-2 data and random forest regression

J. G. N. Irisarri,,
J. G. N. Irisarri1,2,3*M. A. TexeiraM. A. Texeira2P. HarrisP. Harris3K. G. PembletonK. G. Pembleton4
  • 1Department of Ecosystem Science and Management, University of Wyoming, Laramie, WY, United States
  • 2Laboratorio de Análisis Regional y Teledetección (LART), Instituto de Investigaciones Fisiológicas y Ecológicas Vinculadas a la Agricultura (IFEVA), Departamento de Métodos Cuantitativos y Sistemas de Información, Facultad de Agronomía, Universidad de Buenos Aires, Buenos Aires, Argentina
  • 3Rothamsted Research North Wyke, Okehampton, United Kingdom
  • 4Centre for Sustainable Agricultural Systems, University of Southern Queensland, Toowoomba, QLD, Australia

Introduction: In the temperate grasslands of the UK, forage quality is a key factor influencing both animal performance and environmental impact. Because forage quality strongly affects rumen fermentation, improving it can reduce enteric methane emissions and mitigate animal nutritional stress. However, large-scale monitoring of forage quality remains limited due to the reliance on destructive, labor-intensive, and costly sampling methods. Remote sensing offers a promising alternative for scalable monitoring.

Methods: We explored an indicative approach combining optical remote sensing (Sentinel-2) with random forest regression (RFR) models to predict four critical forage quality attributes: crude protein (CP), water-soluble carbohydrates (WSC), neutral detergent fiber (NDF), and acid detergent fiber (ADF). Calibration and validation were performed using >9,500 georeferenced observations collected between 2020 and 2022 at the North Wyke Farm Platform in southwest UK. Forage quality was measured using near-infrared (NIR) sensors mounted on agricultural machinery across paddocks containing permanent and improved pastures. Sentinel-2 spectral predictors included visible, NIR, and red-edge bands, and model performance was evaluated using R² and RMSE metrics.

Results: Model performance was strong across all four forage quality attributes, with R² values ranging from 0.77 to 0.86 and consistently low RMSE values, indicating high predictive accuracy. Red-edge and NIR wavelengths were the most influential predictors. Improved pastures generally exhibited higher forage quality—characterized by lower ADF and higher WSC concentrations—than permanent pastures. Model-predicted seasonal changes were modest, whereas spatial contrasts between paddocks were much more pronounced.

Discussion: The calibrated models are suitable for forage systems with species composition and quality ranges similar to those represented in our dataset but should not be directly applied to other forage types without recalibration. Overall, this work demonstrates the potential of Sentinel-2 remote sensing combined with machine-learning approaches for tolerably accurate, large-scale forage quality monitoring. Such advancements could help improve grazing management, support nutritional planning, and contribute to efforts aimed at reducing methane emissions from livestock systems.

1 Introduction

In temperate regions worldwide, grazing-based livestock production supplies a substantial share of high-quality animal protein and dairy products (1). Forage quality is a key determinant of animal productivity and the environmental footprint (2, 3). In the face of climate change, temperate grasslands are emerging as crucial for food security. Their unique ability to utilize fibrous materials enables sustainable production where other systems may struggle. Detailed monitoring of forage quality dynamics is thus essential to enhance animal performance while minimizing environmental impacts (47). This study focuses on temperate grasslands in the United Kingdom, providing a context where pasture-based livestock production is predominant and forage quality monitoring is critical for sustainable management. This requires addressing a fundamental challenge in grasslands—capturing changes in key forage quality attributes across both time and space with high frequency.

At least four forage quality attributes can be considered critical for animal performance and environmental outcomes: crude protein (CP), water-soluble carbohydrates (WSC), and the fiber fractions neutral detergent fiber (NDF) and acid detergent fiber (ADF). Optimal CP levels enhance feed conversion efficiency via rumen microbial synthesis, while excess CP leads to nitrous oxide (N₂O) emissions and nitrate leaching; mismatches with WSC can also increase methane (CH₄) production through inefficient fermentation (812). Several authors (812) showed that high WSC levels enhance fiber digestibility, leading to greater weight gain and milk production. In contrast, high WSC levels enhance fiber digestibility, leading to greater weight gain and milk production. These authors (812) also found that WSC-rich diets reduce the acetate:propionate ratio in rumen fermentation, favoring propionate, a less methanogenic pathway and potentially lowering enteric CH₄ emissions per unit of animal product. NDF represents the structural fraction of the plant that limits voluntary dry matter (DM) intake. Its digestibility depends on lignin content and microbial activity in the rumen (13, 14). Moderate NDF levels are necessary to maintain rumen motility, but high levels reduce intake and passage rate, limiting animal performance. ADF includes the least digestible forage components (cellulose and lignin) and is a direct predictor of DM digestibility. High ADF levels reduce overall digestibility, lower energy availability, and increase CH4 emissions per unit of metabolizable energy (13). Indeed, digestibility metrics such as organic matter digestibility (OMD) or dry matter digestibility (DMD) are often derived from ADF and NDF values, highlighting the integrative role of these attributes in estimating overall forage quality.

Traditional prediction of these attributes relies on destructive, costly lab analysis with limited scale (15). Machinery-mounted NIR sensors provide high-resolution alternatives (≈1.5 m) during harvests. One such system is the John Deere HarvestLab™ 3,000, which performs real-time NIR-based measurements of forage constituents during harvest. It operates using calibration curves developed in laboratory settings, matched to the reflectance spectra of the crop material. The system is widely used in commercial agriculture and has been certified for accuracy in dry matter estimation by independent organizations such as the Deutsche Landwirtschafts-Gesellschaft, DLG. Sentinel-2 satellite observations, from the visible, red-edge and NIR bands, show strong potential as predictors of forage biochemistry (6, 1621). Spectral bands in the visible, red, red-edge, and NIR regions (bands 2, 3, 4, 5, 6, 7, 8, and 8A) capture variations in canopy structure and composition, which can influence livestock post-ingestive behavioral responses. However, calibrating and validating satellite observations requires reliable and representative reference data.

This study aims to calibrate and validate random forest regression (RFR) models for predicting four key forage quality attributes (CP, WSC, ADF, and NDF) using NIR-harvester data. It also identifies key Sentinel-2 spectral portions for predictions, enabling better grazing decisions that enhance animal welfare by reducing nutritional stress. Models were developed for the agroecological conditions of southwest UK, trained with 70% of the available data and validated with the remaining 30%, using an optimized RFR configuration. This integrated approach is designed to support accurate and cost-effective forage monitoring, with potential applications in pasture management, forage harvest planning, and precision livestock nutrition.

2 Materials and methods

2.1 Ethics statement

No animal studies are presented in this manuscript. No ethical clearance was required as the study did not involve direct animal experimentation.

2.2 Site

The study was conducted at the North Wyke Farm Platform (NWFP), located in southwest England, UK (50°46′12″N, 3°54′05″W), at an elevation ranging from 120 to 190 meters above sea level (22). The region has a temperate humid climate, with an average annual precipitation of 1,031 mm (range: 705–1,361 mm) and average daily minimum and maximum temperatures of 6.8 °C and 13.5 °C, respectively, for the period 1982–2019 (23). The predominant soils belong to the Hallsworth (Dystric Gleysol) and Halstow (Gleyic Cambisol) series, characterized by a moderately stony clay loam surface layer (~36% clay) over a denser mottled clay subsoil (~60% clay), derived from Carboniferous rocks (23).

The NWFP comprises 63 ha with contrasting pasture treatments: permanent (Green farmlet) and improved (Blue farmlet) pastures, representing typical temperate livestock systems (24). The study focused on data from 2020–2022 (24). In our analysis, the forage quality data used for RFR calibration and validation came exclusively from the two pasture types within the Green and Blue farmlets—permanent pastures dominated by Lolium perenne L. (Green) and improved pastures composed of Lolium perenne L. and Trifolium repens L. (Blue). These pasture systems represent typical livestock production systems in temperate environments (24).

2.3 Data

Ground reference data were obtained from a NIR sensor mounted on a John Deere forage harvester, while the remote spectral data were derived from Sentinel-2 satellite imagery. The NIRS sensor used in this study was the John Deere HarvestLab™ 3,000, which employs proprietary calibration models developed from thousands of laboratory-analyzed forage samples. These calibrations are embedded in the device and benchmarked against wet-chemical reference methods, as validated by independent testing (DLG, calibration model LKS 04/18).1 The sensor reports forage constituents including dry matter, crude protein (CP), neutral detergent fiber (NDF), acid detergent fiber (ADF), and water-soluble carbohydrates (WSC), primarily on a dry matter basis. Moisture content is also measured directly, allowing conversion to fresh matter values if needed. While the calibration specifications (e.g., precision, sensitivity) are proprietary, the sensor has been validated for field use and shows performance comparable to laboratory standards. These calibrations were used to generate reference values for model training and evaluation. The sensor operates in the NIRS range (1,100–2,500 nm) and is integrated into the harvester’s monitoring system, allowing data to be recorded at an approximate spatial resolution of 1.5 meters in the direction of travel. This enables detailed characterization of within-paddock variability.

The forage quality data for the NWFP’s Green and Blue farmlets were collected during the 2020, 2021, and 2022 harvest campaigns, at two points in the year: a first harvest in late May or early June, and a second in mid to late August. All sampled paddocks were managed for silage production during the study period. No grazing occurred prior to harvest, ensuring that biomass accumulation reflected ungrazed growth. This management context is relevant for interpreting forage composition, as it tends to result in higher biomass and more uniform canopy structure compared to grazed systems. Data were automatically recorded by the harvester’s integrated system, which includes high-precision GPS, and later downloaded from the John Deere online platform. Each harvester observation represented an approximate area of 4–5 m2, based on a 1.5 m resolution in the direction of travel and the width of the cutting head. Geolocation accuracy was ensured through triangulation with a local ground station at the North Wyke head office, providing sub-meter precision for all georeferenced points. For data analysis, the data were exported in shapefile format. A cleaning process was applied, including the removal of observations with missing values for any attribute and the exclusion of outliers identified through variable relationship inspection. NIRS measurements were performed post-harvest using a sensor mounted on a forage harvester. Therefore, the spectral data used for calibration represent mowed biomass rather than canopy-level estimates.

Sentinel-2 data were extracted via Google Earth Engine (25), selecting the closest low-cloud image (<20% cover, typically within 8 ± 5 days prior, range = 3–13 days) for each NIR point. Bands B2–B8A and B11 were used. The extraction was performed via point sampling on each selected image, and the spectral values were linked to each forage quality sampling point from the harvester. This procedure enabled the association of satellite spectral data with high-resolution ground reference data, facilitating the calibration of forage quality prediction models. Each Sentinel-2 pixel typically aggregated between 12 to 14 NIR observations.

2.4 Data analysis

To predict forage quality attributes from satellite data, RFR models were fitted using functions available in the `randomForest` package in R. Independent RFRs were developed for each of the four attributes: CP, WSC, NDF and ADF. The dataset included a total of 9,531 observations, which were randomly split into 70% for training and 30% for validation. The predictor variables included the optical bands from Sentinel-2 (B2, B3, B4, B5, B6, B7, B8, and B8A). Each RFR was trained using 2,000 trees and five randomly selected predictors (mtry = 5) at each split, values determined through manual tuning that optimized predictive accuracy and model stability. The relative importance of each predictor variable was assessed using the mean decrease in the accuracy metric and visualized through bar plots. Model validation was performed on the test dataset by comparing predicted values with observed values using scatter plots and calculating the Pearson correlation coefficient (r). Hyperparameter tuning was performed manually, without formal grid or random search. Multiple configurations were evaluated heuristically, and the final model parameters were selected based on performance on independent test datasets. In addition to Pearson correlation, we calculated R2, RMSE, and MAE to provide a more comprehensive evaluation of model accuracy and generalization, as shown in Figure 1.

Figure 1
Violin plots on the left display distributions of ADF, NDF, WSC, and CP percentages. Maps on the right show spatial distribution of these components across fields, using color gradients to denote concentration levels, with ADF and NDF in the top row and WSC and CP in the bottom row.

Figure 1. Non-spatial and spatial distributions of the four key forage quality attributes predicted using near-infrared (NIR) sensors mounted on farm machinery (harvesters) by the John Deere system, at six paddocks of the North Wyke Farm Platform (NWFP). The violin plots (left) show the distribution of values for: (A) acid detergent fiber (ADF), (B) neutral detergent fiber (NDF), (C) Sugar (water-soluble carbohydrates), and (D) crude protein (CP). The maps (right) display the spatial variation of each attribute both within and across paddocks. Warm colors indicate higher values. Maps represent the aggregation of all NIRS measurements collected between 2020 and 2023, across multiple sampling events. All values are expressed on a dry matter basis.

3 Results

The four forage quality attributes predicted by the John Deere system showed marked variability across the six NWFP paddocks (Figure 2). WSC ranged fourfold (4–16%), CP twofold (8–16%), NDF 1.5-fold (40–60%), and ADF 1.4-fold (26–36%). These datasets provided the target variables for the RFR models. These datasets, reflecting potential shifts in stem-to-leaf ratios and selectivity in grazing, provided the target variables for RFR models.

Figure 2
Four bar charts labeled A, B, C, and D display the relative importance of spectral bands. Chart A (ADF) shows S2_B8, S2_B6, and others with importance up to 4. Chart B (NDF) indicates S2_B6 and S2_B8 reaching up to 7.5. Chart C (WSC) highlights S2_B6 and S2_B4 up to 5. Chart D (CP) shows S2_B6 and S2_B8 up to 2.

Figure 2. Relative (variable) importance of different Sentinel 2 spectral bands for the random forest regression (RFR) prediction of four forage quality components: (A) Acid Detergent Fiber (ADF), (B) Neutral Detergent Fiber (NDF), (C) Water-Soluble Carbohydrates (WSC), and (D) Crude Protein (CP).

For the fitted Random Forest Regression (RFR) models, Sentinel-2 NIR (S2_B8), red-edge (S2_B6, B7, B8A), and visible (S2_B4) bands were most relevant. SWIR band B11 was not influential, possibly due to limited moisture variability in this system (despite literature support in drier regions), as shown by the variable importance plots (Figure 3). In the case of ADF, the red-edge bands (S2_B6, B7) and narrow NIR (S2_B8A) contributed most to the model’s performance. For NDF, a combination of red-edge and NIR bands were the most influential. For WSC (sugars), the red-edge band S2_B6 and the red band (S2_B4, in the visible spectrum) dominated in importance, indicating a strong relationship between these wavelengths and soluble carbohydrate content. Finally, for CP (crude protein), variable importance was more evenly distributed among visible (S2_B4), red-edge (S2_B6), and NIR (S2_B8) bands, suggesting greater spectral complexity in its prediction using RFR.

Figure 3
Scatter plots show relationships between predicted and observed values for ADF, NDF, WSC, and CP. Each plot includes a black regression line and a diagonal red line. The plots display high correlations, with R-squared values between 0.769 and 0.855, indicating strong predictive performance. RMSE and MAE values are also noted for each plot.

Figure 3. Relationship between observed (John Deere predicted) and random forest regression (RFR) (satellite-informed) predicted values for four forage quality attributes: (A) Acid Detergent Fiber (ADF), (B) Neutral Detergent Fiber (NDF), (C) Water-Soluble Carbohydrates (WSC), and (D) Crude Protein (CP). Each panel shows a scatter plot with the Pearson correlation coefficient (r) between observed values (y-axis) and values predicted by the RFRs (x-axis), using independent validation data.

The study RFRs, informed by satellite data, showed high predictive power for all four forage quality attributes as indicated by the performance metrics (Figure 4). Among traits, CP exhibited the best overall performance (R2 = 0.86, RMSE = 0.68, MAE = 0.44), followed by WSC (R2 = 0.81, RMSE = 1.11, MAE = 0.77) and NDF (R2 = 0.84, RMSE = 1.36, MAE = 0.87). ADF predictions were also robust (R2 = 0.77, RMSE = 1.10, MAE = 0.73), though slightly more dispersed at extreme values. Specifically, the RMSE values indicate typical prediction errors of approximately 0.7 percentage points for CP, 1.1 for WSC, 1.1 for ADF, and 1.4 for NDF. These error levels represent <10% of the observed range for each attribute, confirming the models’ strong predictive accuracy for practical forage quality monitoring. All Random Forest model files (.rds) and accompanying documentation are available through the Zenodo repository.2

Figure 4
Four line graphs depict predicted percentages of ADF, NDF, WSC, and CP over months for improved and permanent pastures. Colored lines and shaded areas represent the different pastures, showing monthly variability and trends.

Figure 4. Seasonal dynamics of forage quality attributes predicted by Random Forest Regression (RFR) models in two functional pasture communities: permanent and improved pastures. Median monthly values are shown for (A) Acid Detergent Fiber (ADF), (B) Neutral Detergent Fiber (NDF), (C) Water-Soluble Carbohydrates (WSC), and (D) Crude Protein (CP). Lines represent the monthly medians, and shaded areas denote the interquartile range, summarizing the spatial and temporal variability in model predictions across the study area. Values integrate the 2020–2022 period, and the RFR models were trained using bi-seasonal field data collected during those years, capturing spatial, seasonal, and interannual variability in forage quality relevant to ruminant nutrition.

Seasonal dynamics of forage quality attributes predicted by the RFR models revealed distinct patterns between the two functional communities, namely the permanent and improved pastures of the NWFP (Figure 1). ADF concentration followed a U-shaped pattern, with minimum values around mid-year and peaks at the beginning and end of the year. In improved pastures, lower ADF levels were maintained from May to August, indicating better forage quality for this structural component. NDF also showed noticeable seasonal fluctuations, with peaks in April, July, and October, and a gradual decline toward the end of the year. Permanent pastures exhibited slightly higher NDF values for most of the year, except during the October peak. In contrast, WSC and CP showed limited seasonal variation, remaining relatively stable throughout the year. Across all attributes, the shaded interquartile ranges revealed substantial temporal and spatial variability, particularly for ADF and NDF during spring and early summer, reflecting heterogeneous canopy structures and phenological stages among fields. Conversely, WSC and CP displayed narrower variability bands during mid-summer, suggesting more uniform forage quality at the peak of vegetative growth.

Spatial patterns of predicted forage quality attributes across the study area showed marked heterogeneity among fields and greater than seasonal shifts (Figure 5). Spatial variability was primarily expressed across fields rather than within them, suggesting that differences in management and vegetation composition were main drivers of spatial heterogeneity in predicted quality attributes. This field-level variability remained consistent throughout the months, underscoring the ability of Sentinel-2 spectral data combined with RFR modeling to detect management-driven contrasts in forage quality across the landscape. Seasonal variation suggested that in April, higher ADF and NDF values were concentrated in the central and northern sectors of the study area. By June, both structural components declined, reflecting the mid-season increase in forage quality. In contrast, WSC and CP exhibited spatially variable but generally higher values in June, coinciding with the peak of the growing season. In September, ADF and NDF increased again, while WSC and CP declined, indicating a late-season reduction in forage quality.

Figure 5
A series of four maps displaying predicted percentages for ADF, NDF, WSC, and CP over April, June, and September. Each map shows geographical variations with color gradients representing different percentage values. The maps have latitude and longitude coordinates, with color bars indicating percentage ranges. ADF ranges from 26 to 32 percent, NDF from 47.5 to 57.5 percent, WSC from 7.5 to 12.5 percent, and CP from 10 to 15 percent.

Figure 5. Spatial patterns of forage quality attributes predicted by Random Forest Regression (RFR) models derived from Sentinel-2 spectral data across the study area. Mean predictions for Acid Detergent Fiber (ADF), Neutral Detergent Fiber (NDF), Water-Soluble Carbohydrates (WSC), and Crude Protein (CP) are shown for April, June, and September, integrating the 2020–2022 period. Each map represents the spatial distribution of model-predicted values averaged per month and attribute, highlighting spatial heterogeneity and seasonal shifts in forage quality. The RFR models were trained using bi-seasonal field data collected during 2020–2022, allowing the predictions to capture both spatial and temporal variation relevant to ruminant nutrition.

4 Discussion

Sentinel-2 data combined with Random Forest regression (RFR) accurately predicted the four key forage quality attributes (CP, WSC, ADF, and NDF), achieving validation statistics comparable to those reported in previous studies (1721). The novelty of this work lies in the integration of satellite data with near-infrared reference information collected directly from a machinery-mounted NIRS system, enabling large-scale calibration under operational farm conditions. This highlights Sentinel-2’s robustness for estimating forage quality in grasslands of the United Kingdom, consistent with results observed in other temperate regions (20, 21). The approach is therefore applicable to forage systems with species composition and quality attribute ranges like those represented in our dataset. Extrapolation to environments with markedly different seasonal dynamics or forage types—such as those where crude protein regularly falls below 6%—would exceed the calibration range and is not recommended.

The variable importance analysis revealed that red-edge, near-infrared (NIR), and shortwave infrared (SWIR) bands (B5, B6, B7, B11, and B12) were the most influential for estimating the four forage quality attributes. This pattern is consistent with previous studies emphasizing the relevance of these spectral regions for assessing structural and chemical properties of vegetation, including lignin, cellulose, and soluble carbohydrates (20, 21). In particular, the red-edge band B7 showed the highest contribution, supporting earlier findings that identify it as a key predictor of dry matter and nutrient content (21). The combined role of red-edge and SWIR wavelengths enhances the prediction of both biomass and forage quality by reducing reflectance saturation in dense canopies (19, 20).

The seasonal dynamics of forage quality attributes showed very small amplitudes (<1 percentage point) compared to the variability observed during calibration (e.g., ADF ranging from 13 to 57% or CP from 6 to 36%). This low variability suggests that forage quality in this system is highly stable throughout the year, despite phenological changes. However, the spatial patterns revealed that fine-scale variability in forage quality across fields exceeded the average intra-annual variation observed within individual pastures. Such strong spatial contrasts likely reflect underlying heterogeneity in management and potentially factors associated with topographic positioning, clearly associated with SOC (26). This pattern is like that observed under similar pastures composition in Southern Australia (27). Moreover, the average values of ADF (~29%), NDF (~50%), and CP (~13%) fall within the upper range compared to global literature, where high-quality forages typically present ADF < 35%, NDF < 50%, and CP between 8–18% (28). Rather than implying broad management recommendations, we emphasize that spatially explicit information systems like this can support adaptive grazing strategies, for example by identifying areas of higher forage quality instead of relying on fixed rotations. It should be noted that our dataset does not include very low biomass levels typical of heavily grazed swards, such conditions often correspond to younger tissues with higher nutritive value, meaning that excluding them could lead to a slight underestimation of forage quality in contexts where quality is already high globally (28). Although improved pastures showed slightly higher forage quality values, the differences were modest. This reflects the fact that the primary goal of the improved pasture system is not only to enhance nutritional value, but to reduce nitrogen inputs and improve environmental sustainability (29).

Beyond decision-support for grazing management, the availability of spatially explicit nutritional information also offers significant value for improving greenhouse gas (GHG) estimation frameworks in pasture-based systems. First, remote estimates of crude protein can support a better understanding of individual dry matter intake as a percentage of live weight, a key determinant of enteric methane production under Tier 2 methodologies (30). Second, forage digestibility directly influences methane emissions per unit of dry matter intake, thus defining emission intensity (g CH₄ kg−1 DMI). While our models do not predict digestibility directly, digestibility (%) can be derived from the predicted ADF and CP using established nutritional frameworks, such as the Digestibility (Van Soest) model (31) or the Digestibility (Dairy model) tailored to temperate production systems (32, 33). By combining satellite-predicted forage constituents with physiologically grounded digestibility equations, this approach has the potential to provide more realistic and spatially refined estimates of enteric methane emissions, thereby supporting climate-smart grazing strategies.

Although the satellite-based models presented here performed well within the conditions of this study, their transferability to other pasture systems must be approached with caution. Machine learning approaches such as Random Forest rely on the coverage of the calibration domain: predictions can become unreliable when applied to forage conditions, botanical compositions, or management regimes that fall outside the biophysical space represented in the training data (34, 35). In our case, the observed ADF, NDF, and CP ranges align closely with those reported for temperate British pasture systems through national initiative known as GrassCheck GB (36), which supports confidence in our predictions. Nonetheless, broader adoption would require local calibration datasets to ensure validity under different livestock species, fertilization regimes, and climatic contexts. This highlights an important limitation of the current work but also a clear pathway toward operational scalability.

5 Conclusion

This study demonstrates the feasibility of integrating NIRS-derived forage quality data with Sentinel-2 spectral information to model and map key nutritional attributes of pastures at landscape scale. By combining multi-year and bi-seasonal field measurements with machine learning techniques, we provide a framework for predicting crude protein, fiber fractions, and sugar content with high spatial and temporal resolution. These findings offer a scalable and cost-effective approach for satellite-based forage monitoring, with direct relevance to temperate humid agricultural pastures. The models developed here can support decision-making in pasture management, forage harvest planning, and precision livestock nutrition. More broadly, this work contributes to the advancement of remote sensing applications in Rangeland and Grassland ecology, enabling researchers and producers to better assess forage dynamics and nutritional value across space and time. The applicability of this approach is limited to forage systems with species composition and quality ranges like those represented in our dataset; extrapolation to markedly different environments would require additional calibration data.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

JI: Validation, Conceptualization, Writing – review & editing, Supervision, Methodology, Data curation, Writing – original draft, Investigation, Formal analysis, Visualization, Software. MT: Investigation, Writing – review & editing, Validation, Methodology. PH: Resources, Writing – review & editing, Visualization, Investigation, Funding acquisition, Project administration, Methodology, Supervision. KP: Writing – review & editing, Investigation, Validation, Methodology, Visualization.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by UKRI-BBSRC (UK Research and Innovation – Biotechnology and Biological Sciences Research Council) through the Resilient Farming Futures Institute Strategic Programme (grant numbers BBS/E/RH/230004B and BBS/E/RH/230004C), the North Wyke Farm Platform National Bioscience Research Infrastructure (BBS/E/RH/23NB0008), and the Sensing Oats with APSIM project (BB/X018539/1). Additional support was provided by the University of Wyoming NASA EPSCoR Faculty Research Grant Program (Grant #80NSSC19M0061).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

References

1. Taube, F, Gierus, M, Hermann, A, Loges, R, and Schönbach, P. Grassland and globalization – challenges for north-west European grass and forage research. Grass Forage Sci. (2014) 69:2–16. doi: 10.1111/GFS.12043

Crossref Full Text | Google Scholar

2. Ban, Y, Neves, LA, Guan, LL, and McAllister, T. Modifying the rumen environment to reduce greenhouse gas emissions. Reduc. Greenhouse Gas Emissions Livestock Prod. (2021):287–330. doi: 10.1201/9781003048213-11

Crossref Full Text | Google Scholar

3. de Ruiter, H, Macdiarmid, JI, Matthews, RB, Kastner, T, Lynd, LR, and Smith, P. Total global agricultural land footprint associated with UK food supply 1986–2011. Glob Environ Change. (2017) 43:72–81. doi: 10.1016/J.GLOENVCHA.2017.01.007

Crossref Full Text | Google Scholar

4. Irisarri, JGN, Durante, M, Derner, JD, Oesterheld, M, and Augustine, DJ. Remotely sensed spatiotemporal variation in crude protein of shortgrass steppe forage. Remote Sens. (2022) 14:854. doi: 10.3390/rs14040854

Crossref Full Text | Google Scholar

5. Raynor, EJ, Derner, JD, Hartman, MD, Dorich, CD, Parton, WJ, Hendrickson, JR, et al. Secondary production of the central rangeland region of the United States. Ecol Appl. (2024) 34:e2978. doi: 10.1002/eap.2978,

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kearney, SP, Porensky, LM, Augustine, DJ, Gaffney, R, and Derner, JD. Monitoring standing herbaceous biomass and thresholds in semiarid rangelands from harmonized Landsat 8 and Sentinel-2 imagery to support within-season adaptive management. Remote Sens Environ. (2022) 271:112907. doi: 10.1016/J.RSE.2022.112907

Crossref Full Text | Google Scholar

7. Kearney, SP, Porensky, LM, Augustine, DJ, Derner, JD, and Gao, F. Predicting spatial-temporal patterns of diet quality and large herbivore performance using satellite time series. Ecol Appl. (2022) 32:e2503. doi: 10.1002/EAP.2503,

PubMed Abstract | Crossref Full Text | Google Scholar

8. Rivera, JE, and Chará, J. CH4 and N2O emissions from cattle excreta: a review of main drivers and mitigation strategies in grazing systems. Front Sustain Food Syst. (2021) 5:657936. doi: 10.3389/FSUFS.2021.657936/XML

Crossref Full Text | Google Scholar

9. Batista, ED, Detmann, E, Valadares Filho, SC, Titgemeyer, EC, and Valadares, RFD. The effect of CP concentration in the diet on urea kinetics and microbial usage of recycled urea in cattle: a meta-analysis. Animal. (2017) 11:1303–11. doi: 10.1017/S1751731116002822,

PubMed Abstract | Crossref Full Text | Google Scholar

10. Rivero, MJ, Keim, JP, Balocchi, OA, and Lee, MRF. In vitro fermentation patterns and methane output of perennial ryegrass differing in water-soluble carbohydrate and nitrogen concentrations. Animals. (2020) 10:1076. doi: 10.3390/ANI10061076,

PubMed Abstract | Crossref Full Text | Google Scholar

11. Lee, MRF, Brooks, AE, Moorby, JM, Humphreys, MO, Theodorou, MK, MacRae, JC, et al. In vitro investigation into the nutritive value of Lolium perenne bred for an elevated concentration of water-soluble carbohydrate and the added effect of sample processing: freeze-dried and ground vs. frozen and thawed. Anim Res. (2002) 51:269–77. doi: 10.1051/ANIMRES:2002023

Crossref Full Text | Google Scholar

12. Villalba, JJ, Ates, S, and MacAdam, JW. Non-fiber carbohydrates in forages and their influence on beef production systems. Front Sustain Food Syst. (2021) 5:566338. doi: 10.3389/fsufs.2021.566338/xml

Crossref Full Text | Google Scholar

13. Jung, HG, and Allen, MS. Characteristics of plant cell walls affecting intake and digestibility of forages by ruminants. J Anim Sci. (1995) 73:2774–90. doi: 10.2527/1995.7392774X,

PubMed Abstract | Crossref Full Text | Google Scholar

14. Mertens, DR, and Grant, RJ. Digestibility and intake. Forages Sci. Grassland Agric. (2020) 2:609–31. doi: 10.1002/9781119436669.CH34

Crossref Full Text | Google Scholar

15. Stuth, J, Jama, A, and Tolleson, D. Direct and indirect means of predicting forage quality through near infrared reflectance spectroscopy. Field Crop Res. (2003) 84:45–56. doi: 10.1016/S0378-4290(03)00140-0

Crossref Full Text | Google Scholar

16. Raab, C, Riesch, F, Tonn, B, Barrett, B, Meißner, M, Balkenhol, N, et al. Target-oriented habitat and wildlife management: estimating forage quantity and quality of semi-natural grasslands with Sentinel-1 and Sentinel-2 data. Remote Sens Ecol Conserv. (2020) 6:381–98. doi: 10.1002/RSE2.149

Crossref Full Text | Google Scholar

17. Zhang, X, Liang, T, Gao, J, Zhang, D, Liu, J, Feng, Q, et al. Mapping the forage nitrogen, phosphorus, and potassium contents of alpine grasslands by integrating Sentinel-2 and Tiangong-2 data. Plant Methods. (2023) 19:1–21. doi: 10.1186/S13007-023-01024-Y/FIGURES/9

Crossref Full Text | Google Scholar

18. Lugassi, R, Zaady, E, Goldshleger, N, Shoshany, M, and Chudnovsky, A. Spatial and temporal monitoring of pasture ecological quality: Sentinel-2-based estimation of crude protein and neutral detergent fiber contents. Remote Sens. (2019) 11:799. doi: 10.3390/RS11070799

Crossref Full Text | Google Scholar

19. Punalekar, SM, Thomson, A, Verhoef, A, Humphries, DJ, Reynolds, CK, Verhoef, A, et al. Assessing suitability of Sentinel-2 bands for monitoring of nutrient concentration of pastures with a range of species compositions. Agronomy. (2021) 11:1661. doi: 10.3390/AGRONOMY11081661

Crossref Full Text | Google Scholar

20. Fernandes, MHM d R, FernandesJunior, J d S, Adams, JM, Lee, M, Reis, RA, and Tedeschi, LO. Using sentinel-2 satellite images and machine learning algorithms to predict tropical pasture forage mass, crude protein, and fiber content. Sci Rep. (2024) 14:8704. doi: 10.1038/S41598-024-59160-X

Crossref Full Text | Google Scholar

21. Morse-McNabb, EM, Hasan, MF, and Karunaratne, S. A multi-variable Sentinel-2 random forest machine learning model approach to predicting perennial ryegrass biomass in commercial dairy farms in Southeast Australia. Remote Sens. (2023) 15:2915. doi: 10.3390/RS15112915

Crossref Full Text | Google Scholar

22. Orr, RJ, Murray, PJ, Eyles, CJ, Blackwell, MSA, Cardenas, LM, Collins, AL, et al. The north Wyke farm platform: effect of temperate grassland farming systems on soil moisture contents, runoff and associated water quality dynamics. Eur J Soil Sci. (2016) 67:374–85. doi: 10.1111/EJSS.12350,

PubMed Abstract | Crossref Full Text | Google Scholar

23. Segura, C, Neal, AL, Castro-Sardiňa, L, Harris, P, Rivero, MJ, Cardenas, LM, et al. Comparison of direct and indirect soil organic carbon prediction at farm field scale. J Environ Manag. (2024) 365:121573. doi: 10.1016/j.jenvman.2024.121573,

PubMed Abstract | Crossref Full Text | Google Scholar

24. Takahashi, T, Harris, P, Blackwell, MSA, Cardenas, LM, Collins, AL, Dungait, JAJ, et al. Roles of instrumented farm-scale trials in trade-off assessments of pasture-based ruminant production systems. Animal. (2018) 12:1766–76. doi: 10.1017/S1751731118000502,

PubMed Abstract | Crossref Full Text | Google Scholar

25. Gorelick, N, Hancher, M, Dixon, M, Ilyushchenko, S, Thau, D, and Moore, R. Google earth engine: planetary-scale geospatial analysis for everyone. Remote Sens Environ. (2017) 202:18–27. doi: 10.1016/J.RSE.2017.06.031

Crossref Full Text | Google Scholar

26. Segura, C, Neal, AL, Castro-Sardiňa, L, Harris, P, Rivero, MJ, Cardenas, LM, et al. Comparison of direct and indirect soil organic carbon prediction at farm field scale. J Environ Manag. (2024) 365:121573. doi: 10.1016/j.jenvman.2024.121573,

PubMed Abstract | Crossref Full Text | Google Scholar

27. Pembleton, KG, Rawnsley, RP, Turner, LR, Corkrey, R, and Donaghy, DJ. Quantifying the interactions between defoliation interval, defoliation intensity and nitrogen fertiliser application on the nutritive value of rainfed and irrigated perennial ryegrass. Crop Pasture Sci. (2017) 68:1100–11. doi: 10.1071/CP16385

Crossref Full Text | Google Scholar

28. Lee, MA. A global comparison of the nutritive values of forage plants grown in contrasting environments. J Plant Res. (2018) 131:641–54. doi: 10.1007/S10265-018-1024-Y,

PubMed Abstract | Crossref Full Text | Google Scholar

29. McAuliffe, GA, Takahashi, T, Orr, RJ, Harris, P, and Lee, MRF. Distributions of emissions intensity for individual beef cattle reared on pasture-based production systems. J Clean Prod. (2018) 171:1672–80. doi: 10.1016/j.jclepro.2017.10.113,

PubMed Abstract | Crossref Full Text | Google Scholar

30. Van Soest, PJ. Nutritional ecology of the ruminant. 2nd ed. Ithaca, NY: Cornell University Press (1994).

Google Scholar

31. Ministry of Agriculture, Fisheries and Food (MAFF). Energy allowances and feeding Systems for Ruminants. Technical bulletin 433. London: Her Majesty’s Stationery Office (1984).

Google Scholar

32. Thomas, C, and Rae, RC. Concentrate supplementation of grass silages. Grass Forage Sci. (1988) 43:223–36. doi: 10.1111/j.1365-2494.1988.tb02154.x

Crossref Full Text | Google Scholar

33. IPCC. 2019 refinement to the 2006 IPCC guidelines for National Greenhouse gas Inventories: Volume 4 agriculture, forestry and other land use. Geneva: IPCC (2019).

Google Scholar

34. Meyer, H, and Pebesma, E. Predicting into unknown space? Estimating the area of applicability of spatial prediction models. Methods Ecol Evol. (2021) 12:1796–810. doi: 10.1111/2041-210X.13650

Crossref Full Text | Google Scholar

35. Robinson, NP, Allred, BW, Jones, MO, Moreno, A, Kimball, JS, Naugle, DE, et al. Regional-scale modeling of vegetation structure and function using remote sensing and random forest. Remote Sens Environ. (2019) 228:73–91. doi: 10.1016/j.rse.2019.04.003

Crossref Full Text | Google Scholar

36. GrassCheckGB. (2025). Grass growth and quality monitoring network. Available online at: https://grasscheckgb.co.uk (accessed October 20, 2025).

Google Scholar

Keywords: sentinel-2, North Wyke farm platform, grassland, crude protein, NDF, ADF, water soluble content

Citation: Irisarri JGN, Texeira MA, Harris P and Pembleton KG (2025) Satellite-based monitoring of forage quality in grasslands of the United Kingdom using sentinel-2 data and random forest regression. Front. Vet. Sci. 12:1678123. doi: 10.3389/fvets.2025.1678123

Received: 01 August 2025; Revised: 13 November 2025; Accepted: 17 November 2025;
Published: 19 December 2025.

Edited by:

Tiago do Prado Paim, Goiano Federal Institute (IFGOIANO), Brazil

Reviewed by:

Nuh Ocak, Ondokuz Mayıs University, Türkiye
Tony Waterhouse, Scotland's Rural College, United Kingdom
Marcia Helena Machado Da Rocha Fernandes, São Paulo State University, Brazil
Karun Kaniyamattam, Texas A&M University, United States

Copyright © 2025 Irisarri, Texeira, Harris and Pembleton. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: J. G. N. Irisarri, amlyaXNhcnJAdXd5by5lZHU=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.