A scalable big data approach for remotely tracking rangeland conditions

Rangelands, covering half of the global land area, are critically degraded by unsustainable use and climate change. Despite their extensive presence, global assessments of rangeland condition and sustainability are limited. Here we introduce a novel analytical approach that combines satellite big data and statistical modeling to quantify the likelihood of changes in rangeland conditions. These probabilities are then used to assess the effectiveness of management interventions targeting rangeland sustainability. This approach holds global potential, as demonstrated in Mongolia, where the shift to a capitalist economy has led to increased livestock numbers and grazing intensity. From 1986 to 2020, heavy grazing caused a marked decline in Mongolia ’ s rangeland condition. Our evaluation of diverse management strategies, corroborated by local ground observations, further substantiates our approach. Leveraging globally available yet locally detailed satellite data, our proposed condition tracking approach provides a rapid, cost-effective tool for sustainable rangeland management.

Productivity (Aus RAPP) platform, drawing on MODIS data, offers broader global coverage 9 .However, it faces limitations due to the inherent constraints of MODIS imagery, including a reduced spatial resolution (500 m) and a shorter temporal span (available since 2000).Meanwhile, other existing methods projecting global rangeland changes rely on modeled data, presenting long-term forecasts at even coarser resolutions (>25 km) 10,11 .Although finer mapping techniques (e.g., ≤30 m pixels) exist for diverse ecosystems (e.g., forests, croplands, and wetlands [12][13][14] ), a comparable highresolution framework for global rangeland analysis remains absent.Our study aims to bridge this gap, laying the foundation for advanced rangeland monitoring and management strategies (see Table S1).
Rangeland condition tracking needs to be done in a way that can effectively support decision-making for greater sustainability.This is particularly challenging for rangelands which, compared with forests, tend to exhibit fewer natural thresholds (e.g., forest conversion) under their dominant use (grazing) and experience high inter-and intra-annual variability predominantly due to changes in rainfall.Effective decision support therefore requires a flexible but robust statistical approach that can illuminate when changes of concern to decision-makers are observed and how likely these are.
To address the above mentioned gaps in the field, this study developed a novel analytical approach to map rangeland condition changes at 30 m spatial resolution.We illustrated this approach over Mongolian rangelands, which can further extend globally in the future.Based on the new capabilities of the long-term global Landsat image archive, we incorporated dynamic statistical modeling into the proposed approach to track the rangeland condition changes as well as to quantify the likelihood of such changes.Therefore, our approach facilitates the use of accumulated satellite data to monitor multiple levels of change in the rangeland system (as defined by end users) and to identify the support for each level of change given the data.Such information can then be used to quantify management impact and define a level of certainty regarding this impact, delivering improved knowledge for global rangeland management and conservation investments.

Results
Initial status of rangeland condition Identifying rangeland condition change requires the determination of a prior or initial condition to assess this change against.We used a benchmarking method to identify the condition of rangelands prior to 1990, which quantifies local rangeland productivity relative to the potential productivity in the same Land Capability Class (LCC) of homogeneous lands with similar biogeophysical capability under similar climate conditions.A total of 4918 LCCs were defined across the continent of Asia and 1962 classes over Mongolia (Fig. 2a).Then, we identified areas within each LCC that possessed the minimum (lowest tenth percentile) or maximum (highest tenth percentile) productivity performance possible for the LCC they belonged to (red and green areas respectively in Fig. 2b).Areas in red in Fig. 2b indicate the lowest production relative to their potential, assuming the presence of rangeland vegetation.This suggests that these rangeland areas might have already experienced severe land degradation at the beginning of our study period.As a result, pixels within these regions that have 'stable condition' modeled in our subsequent analysis indicate rangelands that stay consistently degraded during the study time.In contrast, areas in green in Fig. 2b show rangelands with the highest relative production, in other words reaching nearly their full potential within each LCC.
Rangelands varied across Mongolia in their initial conditions relative to the maximum productivity possible within each LCC (Fig. 2c), which present near full productivity potential in the center of the country while reductions in parts of northern and southern Mongolia prior to 1990.In addition, transitions in land cover within Mongolia from 1992 to 2020 were investigated with land cover yearly data and observed to be few, with only a minor proportion of rangelands in the north (almost none in the south) being converted to other land cover types (Fig. 2d).These areas were https://doi.org/10.1038/s43247-024-01516-2included in our rangeland condition change analysis but excluded in the process of our management assessment.

Trajectories of rangeland condition change
The spatial distributions of overall changes in rangelands across Mongolia were derived from the 35-year window of the entire study time (1986-2020), and represented by the RGB composited belief images of the three condition models (Fig. 3a).In contrast to the spatial patterns of initial condition that indicated maximum productivity in the center of Mongolia and lower productivity elsewhere prior to 1990 (Fig. 2c), we found rangelands across the center of Mongolia with dense population and intensive grazing activities showed overall higher belief in 'decline' (red areas in Fig. 3a) compared to other areas.Specifically, most of the southern Mongolia desert region, where vegetation dynamics are being driven by climate variability, presented improved rangeland conditions for the past 35 years, with the green color indicating a high belief in 'increase' condition.
The annual time series of statistical beliefs in decline, increase, or stable rangeland status were derived from the accumulated vegetation condition changes during each previous decade (Fig. 3b), which may show details of the overall condition changes geographically presented in Fig. 3a.Such time series of beliefs in condition changes enable quantitative assessment of rangeland management impacts over both large and small areas.For example, the black box indicated in Fig. 3b is a protected area located within an overall improved condition environment.For this protected area there was a dramatic jump in the trajectory of beliefs in 'increase' (red circled point in Fig. 3b), which may indicate the implementation of a successful management intervention in 2012.

Management impacts
We evaluated the impacts of rangeland management within each protected area individually and then synthesized the findings to compare the effectiveness of various management strategies implemented across Mongolia.With the abovementioned map series of beliefs in the 'increase' condition, management impacts were quantified as the difference of beliefs in the 'increase' condition over each protected area between the prior-and postmanagement periods.Since rangelands have been protected across  Mongolia by three levels of management (National Protected Areas with four sub-management strategies, NPAs; Local Protected Areas, LPAs; and Community-Based Organisation, CBO), we focused our analysis on the country-wide effects of these different management strategies.
The spatial distributions of the belief difference in 'increase' over local protected areas before and after the management are shown in Fig. 4a, with green-colored areas indicating improved rangeland condition after the management while those in red presented no improvement.We found over 74% of the LPAs showed higher belief in the 'increase' condition after the management applied.However, such improvement over a large proportion of the protected areas could have resulted from the combined impacts of both natural climate and human management interventions.Indeed, we found there has been an increasing trend in rainfall across Mongolia since 2011 that may have positively impacted rangelands over the last decade (Fig. 4b).As a result, we need to separate the management impact from those due to inter-annual climatic variability on rangelands to better understand the benefits from rangeland management interventions.
Using all areas under the same LCC as a reference, we removed the rainfall effect by subtracting belief in the 'increase' of protected areas from that of the reference sites for the pre-and post-management periods ("Methods" section).We found after the removal of the rainfall effect, the percentage of local protected areas showing higher belief for the postmanagement period reduced to around 25% (Fig. 4c).This indicates the ecological improvement in most of the local protected areas was due to climate variability (i.e., rainfall increase) rather than management effects.More specifically, 25% of local protected areas showed condition improvement linked to management, rainfall contributed to the vegetation increase for half of the areas, and the other 25% didn't show significant improvement under either climate or management effect.

Comparison among different management strategies
In contrast to local protected areas, the national protected areas (NPA) presented 100% improvement under the context of both management and climate effects (Fig. 5a); 77% of these benefits were due to management alone after rainfall correction (Fig. 5b).Being both equally distributed across Mongolia, the management of national protected areas outperformed that of the local protected areas, with higher median belief difference in 'increase' and less variance (Fig. 5c).We also found the earlier NPA management establishment was likely to link with higher belief in 'increase' after the management (Fig. 5c).In addition, the protected areas on both levels showed land degradation around Mongolia's capital city of Ulaanbaatar, where intensive grazing activities occur such as areas in the central and eastern grassland as well as northern mountain grassland ecoregions (Figs.4c and 5b).
The effectiveness of the four different sub-management strategies under the NPA management was also analyzed in this study, including the 'strictly protected area', 'national parks', 'nature reserve', and 'natural and historical monument areas'.We found overall improvement due to all four NPA strategies, with the median belief difference in 'increase' above zero between before and after all the management applied (Fig. 5d).'Strictly protected areas' were observed to have the most effective management among all four, followed by the 'national parks' while 'nature reserve' and 'natural and historical monument areas' showed relatively lower effectiveness.The improved rangeland condition in the community-based organization areas (highlighted in red polygons in Fig. 5a, b) was observed to be solely attributed to rainfall.However, the management of these areas investigated in our study was only recently implemented (late 2018 and 2020), which likely limited the extent of any potential impact on rangeland conditions.

Discussion
The validation of our rangeland condition changes and management assessments, using a triad of methods informed by field observations and regional reports on Mongolia's rangeland health, yielded coherent results.Firstly, our statistical beliefs in 'increase', 'stable', and 'decline' conditions showed strong spatial concordance with the observed degradation levels across Mongolia (Fig. 6a).For instance, areas with high beliefs in 'decline' predominantly encompassed locations of fully and heavily degraded field observations, whereas areas with low beliefs in 'decline' largely corresponded with healthy field sites.Furthermore, the density distribution plot of validation points (Fig. 6b), which illustrates the beliefs in 'decline' across different degradation levels, unveiled a discernible pattern: the central tendency and dispersion indicate that areas with more severe degradation observed in the field are associated with higher beliefs in 'decline' in our analysis.Initially, Mongolian rangelands (1986-1990) exhibited minimal productivity reduction in central and eastern regions (Fig. 2c), yet these areas later displayed high beliefs in 'decline' from 1986 to 2016 (Fig. 6a).This shift could be attributed to Mongolia's economic transition post-1990, which saw a surge in subsistence herding and migration towards provinces near Ulaanbaatar due to job losses and a tripling of livestock numbers 15 .The confusion matrix analysis underlines our model's proficiency and challenges in classifying rangeland conditions (Table 1).The matrix indicates a User's Accuracy of 76.7% for 'Degradation' and 81.8% for 'Healthy', with corresponding Producer's Accuracies of 71.7% and 85.4%, respectively.This reflects a solid capability in identifying both degraded and healthy rangelands, supported by an Overall Accuracy of 79.9% (±2.56).However, the reduced accuracy observed in the 'Degradation' category underscores the complexities of accurately identifying areas of true degradation.This may be attributed to the limitations of remote sensing to detect early-stage degradation, the masking effects of human activity and climate change on degradation signals, and our preliminary application of basic models with a theoretical 20% threshold for change.This threshold merits adjustment by rangeland managers utilizing their practical insights.
Our results from assessing the impact of different rangeland management strategies showed good agreement with the impacts of management strategies described in the regional reports 16 .For example, the four categories of management applied to the national protected areas varied in their effectiveness with 'strictly protected areas' showing the most benefits to vegetation increase in rangelands followed by 'national parks', which outperformed 'nature reserve' and 'natural and historical monument areas' (Fig. 5d).This is because 'strictly protected areas' and 'national parks' have been managed by different protected area administrations from the other two, the former of which received more funding from the central government to have better governance 16,17 .This was further confirmed by the qualitative assessments of previous government research on the comparison among the four management strategies over different small regions via the Management Effectiveness Tracking Tool developed by WWf [18][19][20] .In addition, we found the management of national protected areas had much more positive impacts on protected rangelands than that implemented for local protected areas.This may be attributed to the limited financial capacity of local governments (Province and county levels), to provide sufficient funding for the management of local protected areas 17 .As we observed, however, the time since management establishment is an important factor in the rangeland condition trend (Fig. 5c), suggesting some of the differences between national and local protected areas captured here could be attributed to the earlier implementation of the national protected areas.Moreover, results from this study might be subject to uncertainties stemming from data limitations, methodological constraints, and local environmental factors.For example, the use of Land Capability Class, designed to assess agricultural potential based on soil properties, lacks specific linkage to ecosystem functionality or vegetation communities.As global data and methodologies evolve, further efforts are necessary to refine or replace LCC, however, our approach is flexible to such adaptations.Similar to conventional remote sensing methods, our study's focus on the greenness index mainly captures grassland cover changes, possibly missing rangeland areas with significant woody and invasive species encroachment.Future efforts could integrate additional field data on invasive species abundance, and apply ecological models to understand the spread and impact of invasive species on grasslands.In general, the ability to incorporate more diverse data sources will allow for a more comprehensive management impact analysis.
The proposed approach in this study consists of simple but robust statistical algorithms, which can be achieved by satellite observations and spatial products that are globally freely available without being dependent on the field observations.Thus, it has great future potential to extend to a global scale with even higher spatial resolution satellite imagery (e.g., sentinel-2 or planet) via cloud computing platforms such as GEE.In addition, this study enables flexible models of change to be identified.First, the linear percentage change can be replaced by a non-linear if needed.Further, the level of change of interest, e.g., ±20% investigated in this study, is flexible, and any, or multiple, values could be assigned, and more than three models could be identified and assessed within this approach.The assessment window designed in this study is also flexible and could track a single year or any length up to the entire period of the data.There is a welldocumented concern with methods that set a binary threshold (e.g., a pvalue) to assess whether an area is degraded or its condition has changed 21 .This study avoids using such thresholds and instead focuses on likelihoods associated with different extents of change.Avoiding a standardized approach to defining rangeland conditions provides flexibility, so end users can look for changes and levels of confidence in those changes, which reflect their own beliefs and choices about sustainability.

Conclusions
Despite their vast extent and the crucial resources that rangelands provide to people and nature, we currently have a limited accounting of rangeland's condition.In this study, we proposed a novel analytical approach that does not compromise on spatial resolution or time scale and avoids standardized definitions of rangeland condition change.This approach links three features critical for a significant advance in rangeland research: (1) identifying the initial state of rangeland condition before the satellite era; (2) mapping rangeland condition changes and estimating the likelihood of such changes at high spatial resolution (30 × 30 m) over long-time horizons (40 years); and (3) quantifying the impacts of different large-scale management on rangeland condition.Our initial findings show the new approach can be applied over the rangelands of Mongolia to detect and identify landscape changes especially in vegetation cover, at scales that can be used to inform the development and assessment of sustainable rangelands management programs.Details of these management programs need to be regionally specific and are beyond the scope of this work.

Methods
The proposed approach encompasses a four-stage process, as outlined below (Fig. 7), and has been developed utilizing the Google Earth Engine (GEE) platform.The data employed in this study are freely accessible on a global scale, with the majority pre-existing within GEE, while the remainder has been uploaded to the user's asset space on GEE.

Study area
We used Mongolian rangelands as a case study to test our proposed approach, which is an area that covers 80% of the country (over 120 million hectares).Rangelands in Mongolia support approximately one million nomadic herders (30% of the population) and provide refuge for a variety of rare and wide-ranging wildlife 22 .Rangelands may be degrading due to loss of traditional grazing practices, and increases in livestock numbers (tripled during the past three decades), combined with impacts of mining and infrastructure development.The degradation caused by these impacts has been accelerated by climate change and population growth 23 .In recognition

Landsat data and vegetation index selection
We used the entire Landsat archive of surface reflectance data from 1986 to 2020 (Table 2), which have been corrected for atmospheric, reflectance, topographic, and satellite sensor effects 25 .Spatially, 129 Landsat scenes were required from World Reference System path 122-144 and rows 23-31 to   26,27 .Initially vegetation index, as a biophysical measure for changes in the amount of green vegetation/photosynthesis in rangelands 28,29 , was used in this analysis to demonstrate the feasibility of the research design and approach over the spatial and temporal scales.This will provide a basis to transition to other products such as vegetation fractional cover, which with the inclusion of field observations may be more effectively linked to local management activities but are not currently available at high resolution and globally 9 .The Enhanced Vegetation Index (EVI) was developed as an alternative vegetation index of the Normalized Difference Vegetation Index (NDVI) with reduced soil and atmospheric influences and improved sensitivity in high biomass regions 30 .However, EVI requires a blue band, which can be sensitive to variations in the viewing geometry, surface albedo, sun angle, and terrain 31 .For example, here we found such erratic blue behavior often causes upward spikes in EVI over many parts of Mongolia.To address this issue, we chose to use two-band EVI2 as a surrogate of EVI to avoid data noise over complex environments such as landscapes in Mongolia.EVI2 has been proven to be a good indicator for above-ground biomass 32 , which has been widely used for Mongolian rangeland studies 33,34 .EVI2 was calculated as below 35 : Where ρ nir and ρ red represent surface reflectance in the near-infrared and red bands.Annual composites of EVI2 for the period of 1986-2020 were then computed for this analysis based on the median of the year, which is a specific data point instead of an averaged or blended value and robust against extreme values and outliers 25 .

Management data
Rangeland protected area data provided by TNC Mongolia was used for management assessment in this study, which was applied during the last three decades to balance conservation and development in Mongolia (Table 2).They are under three levels of management strategies, including National Protected Areas (NPA), Local Protected Areas (LPA), and Community-Based Organization areas (CBO).The Mongolian government aims to establish 30% of the country's land as national protected areas and has already reached 20%.TNC Mongolia has been collaborating with the TNC's global team over the past decade to conduct ecoregional assessments and landscape-level planning across Mongolia, helping central and local governments to establish new protected areas.Traditional grazing of livestock is currently allowed in all local protected areas and in parts of national protected areas except the strictly protected NPAs.In addition, the Mongolian government and TNC support local communities to actively participate in the development of rangeland management plans within local protected areas and CBO areas (in many cases these areas overlap).Communities can improve local rangeland management by implementing their traditional sustainable grazing practices, limiting the inflow of herders from other areas, and preventing the establishment of new mining concessions.

Precipitation
We obtained daily precipitation data on GEE from Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) to assess the interannual climatic variability over Mongolia (Table 2).CHIRPS is a quasiglobal rainfall dataset from 1986 to the present, spanning 50°S-50°N 36,37 .It incorporates 0.05°resolution satellite imagery with in-situ station data to create gridded rainfall time series suitable for trend analysis and seasonal drought monitoring.Annually summed precipitation time series were used in this analysis.

Biogeophysical data and land capability classification
'Land capability' is a widely used concept in agricultural science, indicating the suitability of the land for a specific use such as rangeland or rain-fed cultivation 38 .We used Land Capability Class (LCC) to form the basis for identifying the rangeland initial condition, which represents areas with similar land capability under homogeneous biogeophysical and climatic environments (Table 2).K-Means clustering, which is an unsupervised machine learning algorithm, was used to classify the LCCs based on unique and comprehensive layers of independent variables with all data available at a global scale.This data included: Land cover -Land cover map at 300 m resolution with 22 major classes from 1992 to 2020 provided by the ESA Climate Change Initiative 39 , which was derived from five different satellite missions based on the United Nations land cover classification system; Soil units -soil taxonomy at 250 m with 12 classes 40 , which is distribution of the USDA soil great groups based on machine learning predictions from global compilation of soil profiles (>350,000 training points); Terrain -landform at 250 m with 5 major classes produced using the improved Hammond Landform Classification Algorithm 41 ; Climate -Koppen-Geiger climate zones at 1 km resolution for the period of 1986-2020 42 , which was derived from an ensemble of 32 climate models based on temperature, precipitation, potential evapotranspiration, etc with 13 major classes (at 1st and 2nd levels).

Initial rangeland condition identification
Being confined by the satellite lifespan, remote sensing methods can only detect lands that have been actively affected by the ongoing processes of degradation since the 1980s 43 .This means rangelands with the initial condition of severe land degradation may show little decline during the satellite time.As such, identifying the initial condition of rangelands is the crucial first step for a better understanding of the changes in the condition of these lands during our study period.To do this, the Local Net Production Scaling (LNS) was used 44 , which as a benchmarking method works on condition that: (1) it is based on the LCC with the same rangeland productivity level, and (2) non-degraded reference sites exist in each LCC.We thus defined LCC over the whole continent of Asia instead of Mongolia to ensure there are enough reference sites.The LNS index was calculated across the averaged annual EVI2 image over the first five years of the study time (1986-1990), which was derived by scaling the actual EVI2 of each pixel relative to the higher observed EVI2 value in the frequency distribution of all pixels in the same LCC.The areas with full potential and minimum productivity performance of each LCC class were estimated by the EVI2 values at 90th and 10th percentiles.We set the LNS index of pixels with EVI2 being in the bottom and top 10th percentiles to 0 and 100 respectively to reduce the effect of very high or very low annual EVI2 values caused by atypical areas with low frequency (e.g., unmapped cultivation or bare rock).The LNS index of the other pixel percentiles was scaled linearly between these limits as follows 38,45 : If qi ≤ q0:1; LNS i ¼ 0; or If qi ≥ q0:9; LNS i ¼ 100; or If q0:1 ≤ qi ≤ q0:9; LNS i ¼ 100 ðqi À q0:1Þ=ðq0:9 À q0:1Þ; ð2Þ where q 0.1 , q 0.9 and q i are EVI2 values at 10th, 90th and in between percentiles respectively.
Then, areas with LNS < = 10 present the lowest production relative to their capability, i.e., already under severe land degradation at the start of our study period.As a result, we set the belief in 'decline' as 100% for pixels within these areas that have 'stable condition' modeled in our following analysis, indicating rangelands under consistent degradation during the study time instead of being stable.

Rangeland condition change modeling
We defined three condition models at each Landsat pixel location area (30 × 30 m) to represent the significant change (decline/increase) and stable conditions of rangelands (Fig. 2).Previous studies showed 10-20% reduction in land productivity (represented by vegetation greenness) was observed between degraded and non-degraded areas within arid and semiarid regions 46,47 .Therefore, the three models for Mongolian rangelands in this analysis were expressed as simple linear equations, i.e., 20% decline or 20% increase or similar as mean from the start to the end of the study period: Where y represents the predicted EVI2 value, y m is the mean of annual EVI2 during the study period, t is time at the yearly step, and T is the end year of the study period.The threshold z for the significant change can be modified according to the local environment by end users in the future.The three defined rangeland condition models were further applied through a moving-window method, which shifts a kernel of the 10-year period over the EVI2 time series of each pixel, sliding one year step each time to create an annual belief time series from 1995 to 2020.This allowed candidate condition models to be dynamically built at each window to model the rangeland condition changes over time.Rather than the traditional null hypothesis testing approach with an arbitrary probability threshold (e.g., P-value < 0.05) that is normally used 48 , this study used a formal model selection procedure to identify the relative support for each of the three models described in Eq. ( 3), given a 10-year time series of EVI2 values for each pixel.Specifically, Akaike Information Criterion (AIC) weights were calculated for each of the three competing models to provide a quantitative measure of belief (scaled between 0 and 1; all beliefs sum to 1) for each rangeland condition hypothesis 49 .AIC uses maximum likelihood scores as a measure of fit, enabling inference to be drawn from multiple candidate models simultaneously.Thus, with each 10-year moving window, rangeland condition change was calculated at pixel level across Mongolian rangelands.The derived time series of three weights represents the relative likelihood of the condition change type, indicating the statistical beliefs in decline, increase or stable rangeland status from the accumulated vegetation condition changes during each previous decade.

Accuracy assessment
Field observations of rangeland conditions across Mongolia for validating our results were collected from a nationwide rangeland survey program by the National Agency for Meteorology and Environmental Monitoring (NAMEM) in 2016, covering most of the vegetation types at 1040 monitoring sites across the country.The NAMEM further converted the observations to 5 different degradation levels by using an interpretation tool named Ecological Site Descriptions.These degradation levels were carefully established by comparing the vegetation health at each site against that of healthy counterparts within the same specific vegetation communities.They describe a sequence of changes from healthy rangeland conditions toward unproductive states based on key criteria such as species composition, bare soil cover, total species number, the proportion of palatable and degradation indicator species, litter accumulation, and above-ground biomass (details in Table 3).
The accuracy assessment of this study was carried out using both qualitative and quantitative methods based on field data.Qualitatively, we examined the spatial coherence between field observations and changes in rangeland conditions depicted in the map, alongside its concurrence with local environmental and management reports.Quantitatively, we employed a confusion matrix to calculate Overall Accuracy (OA), Producer's Accuracy (PA), User's Accuracy (UA), and Balanced Accuracy (BA).In constructing the confusion matrix, we generated a dominant rangeland condition change map from 1986 to 2016 using the maximum value rule across beliefs in "Increase", "Stable", and "Decline" maps.To harmonize the five degradation levels from field data with the three categories on our rangeland change map, we grouped them into two classifications: "Healthy" and "Slightly degraded" levels from field observations, along with "Increase" or "Stable" conditions from our map, were designated as the "Healthy" class; conversely, "Moderately degraded", "Heavily degraded", and "Fully degraded" levels from field data, matched with the "Decline" condition from our map, were classified as the "Degradation" class.
Separation of natural variability from rangeland management impacts Changes in rangeland condition represented as a belief in the support of the data for different models of change (i.e., decline, increase, and stable) were output as time series over a 10-year moving window and over the whole study period (1986-2020) for each pixel.These trajectories serve as a foundation for exploring condition change in rangelands under differing perspectives, including tracking historical trajectories and spatial patterns of rangeland change in response to different land management decisions.Here we utilized the time series of belief in 'increase' to quantify the effectiveness and spatial distribution of benefits associated with protected area policies across Mongolia.
To ensure the difference observed pre-and post-management periods from these belief trajectories represents a signature of the management and not inter-annual climatic variability inherent in rangeland systems, it was essential to remove the impact of this variability from the trajectories.Such a need has been recognized in previous studies 50,51 .The most common method, Residual Trends Analysis (RESTREND), compares expected plant growth based on rainfall patterns with actual satellite observations to remove the effects of climate variability.However, this method has been proven unreliable when the change observed is greater than 20% 46 .Given our study aims to be flexible to the level of change to be detected, to complement targets, decision triggers, or reporting needs of many end users in these vast rangeland systems, we had to explore another approach.We instead utilized a method based on LCC to isolate the effects of human management 38,44 .
To remove the rainfall variability from our belief calculation, we developed a continental scale reference value for each LCC across Asia for comparison with the managed areas in this study, in this case, protected area declaration and management.By subtracting the belief in an 'increase' in condition after protection with the average for all corresponding reference sites over the same time horizon, we effectively removed rainfall variability and allowed change to be more likely a signature of management implemented.The difference in the belief of an increase in protected areas was mapped to show the spatial distribution of the management impact, as well Table 3 | The basic criteria followed to classify the field observations of rangeland conditions to 5 degradation levels

Degradation levels Criteria
Healthy All native dominants are in place.
Slightly degraded Key dominants are still dominating, some grazing-sensitive forbs are in decline and grazing-resistant species are in increase.
Moderately degraded Dominants are in decline and replaced by other subdominants, number of species drops down.

Heavily degraded
Remnants of key species are thinning, and the abundance of degradation indicator species increases.

Fully degraded
Total vegetation cover is reduced or dominated by very few degradation indicator species.
as maps of belief when climate variability is not removed (Figs. 4 and 5).land cover changes, such as transitions from rangeland to other land types or vice versa, were identified and excluded from the management assessment analysis.This step was taken to prevent any misattribution of changes in rangeland conditions.

Fig. 3 |
Fig. 3 | Spatio-temporal condition changes of rangelands in Mongolia from 1986 to 2020. a RGB composite map of beliefs in 'increase', 'stable', and 'decline' conditions modeled from the 35-year window of the entire study time (n.b., the vertical stripes were due to the mosaic of satellite image scenes).b Map series stack of beliefs in 'increase' from decadal moving-window sliding one year step each time, with each year on the x-axis representing the belief in rangeland productivity increase during its previous 10-year.For example, a dramatic jump in the statistical belief (circle in red) in 2012 shows an improved rangeland condition during 2003-2012.

Fig. 2 |
Fig. 2 | Rangeland initial condition.a Land capability classes (n.b., colors are only for illustration purposes and not enough for all LCCs).b Extremes of initial rangeland productivity performance (lowest and highest tenth percentiles) across the continent of Asia.c Initial rangeland productivity across Mongolia, scaled by their maximum production.d Land cover transitions over Mongolia from 1992 to 2020.

Fig. 4 |
Fig. 4 | Belief difference in 'increase' between the prior-and post-management periods.a Spatial distributions of the combined impacts from management and climatic variability on local protected areas.b Annual rate of change in precipitation across Mongolia during 2011-2020, compared to time series of spline-smoothed yearly mean precipitation over four local protected areas (labeled as S1, S2, S3, and S4) from 1986 to 2020.c Spatial distributions of the sole management impacts on local protected areas.

Fig. 6 |
Fig. 6 | Spatial alignment between field observations and mapped rangeland condition changes across Mongolia.a Rainfall corrected belief in the 'decline' map over a 31-year period of 1986-2016 with field observations of five different degradation levels acquired in 2016.b Corresponding rainfall corrected beliefs in 'decline' for all the validation points across different degradation levels.

Fig. 5 |
Fig. 5 | Comparison of different management strategies.a Spatial distributions of the combined impacts from management and climatic variability on national protected areas.b Spatial distributions of the sole management impacts on national protected areas.Boxplots for c local (LPA) and national (NPA) protected areas colored by the time of management establishment and d protected areas of the four different management strategies under the NPA management (SPA 'strictly protected area', NP 'national parks', NR 'nature reserve', NHM 'natural and historical monument areas').Values in (a-d) are belief differences in 'increase' between the prior-and post-management periods (i.e., after the correction of the climatic effect).

Table 1 |
Accuracy statistics of the rangeland condition change map from 1986 to 2016 derived from the confusion matrix analysis 24 the importance of sustainable rangeland use, the Mongolian government has been working with The Nature Conservancy to strengthen the management of protected areas to balance conservation and development.Three levels of land protection at different spatial scales, covering most of Mongolia, have been gradually established since the 1990s including National Protected Areas (NPA), Local Protected Areas (LPA), and more recently, Community-Based Organization (CBO) areas24(Table2).The large extent of rangelands with a long-term record of different rangeland management interventions makes Mongolia an ideal region to test our analysis for global rangeland application.

Table 2 |
Summary of data used in this study

Table 2 .
mosaic and cover the entire Mongolia.This resulted in the acquisition of totally over 100,000 images from GEE's median across different Landsat sensors, including the Landsat Thematic Mapper (TM, Landsat 4-5), the Enhanced Thematic Mapper Plus (ETM+, Landsat 7) and the Operational Land Imager (OLI, Landsat 8) at 16 days and 30 m resolutions.The quality assessment flags provided in Tier 1, Collection 1 product were used to remove undesirable elements such as aerosol, cloud, cloud shadow, and inundation