# Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands

Jinfeng Chang, Philippe Ciais, Thomas Gasser, Pete Smith, Mario Herrero, Petr Havlík, Michael Obersteiner, Bertrand Guenet, Daniel Goll, Wei Li, et al.

► To cite this version:

Jinfeng Chang, Philippe Ciais, Thomas Gasser, Pete Smith, Mario Herrero, et al. Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands. Nature Communications, 2021, 12 (1), 10.1038/s41467-020-20406-7. hal-03136901

HAL Id: hal-03136901

https://hal.science/hal-03136901

Submitted on 17 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
ARTICLE

https://doi.org/10.1038/s41467-020-20406-7

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands

Jinfeng Chang ID 1,2,3✉, Philippe Ciais ID 1, Thomas Gasser 2, Pete Smith ID 4, Mario Herrero ID 5, Petr Havlík ID 2, Michael Obersteiner ID 2, Bertrand Guenet ID 1, Daniel S. Goll ID 1, Wei Li 1, Victoria Naipal ID 6, Shushi Peng ID 7, Chunjing Qiu 1, Hanqin Tian ID 8, Nicolas Viovy 1, Chao Yue 9 & Dan Zhu ID 1

Grasslands absorb and release carbon dioxide (CO2), emit methane (CH4) from grazing livestock, and emit nitrous oxide (N2O) from soils. Little is known about how the fluxes of these three greenhouse gases, from managed and natural grasslands worldwide, have contributed to past climate change, or the roles of managed pastures versus natural grasslands. Here, global trends and regional patterns of the full greenhouse gas balance of grasslands are estimated for the period 1750 to 2012. A new spatially explicit land surface model is applied, to separate the direct effects of human activities from land management and the indirect effects from climate change, increasing CO2 and regional changes in nitrogen deposition. Direct human management activities are simulated to have caused grasslands to switch from a sink to a source of greenhouse gas, because of increased livestock numbers and accelerated conversion of natural lands to pasture. However, climate change drivers contributed a net carbon sink in soil organic matter, mainly from the increased productivity of grasslands due to increased CO2 and nitrogen deposition. The net radiative forcing of all grasslands is currently close to neutral, but has been increasing since the 1960s. Here, we show that the net global climate warming caused by managed grassland cancels the net climate cooling from carbon sinks in sparsely grazed and natural grasslands. In the face of future climate change and increased demand for livestock products, these findings highlight the need to use sustainable management to preserve and enhance soil carbon storage in grasslands and to reduce greenhouse gas emissions from managed grasslands.

<footnote n="1">¹Laboratoire des Sciences du Climat et de l'Environnement, LSCE/IPSL, CEA-CNRS-UVSQ, Université Paris-Saclay, 91191 Gif-sur-Yvette, France. ²International Institute for Applied Systems Analysis, A-2361 Laxenburg, Austria. ³College of Environmental and Resource Sciences, Zhejiang University, 310058 Hangzhou, China. ⁴Institute of Biological & Environmental Sciences, University of Aberdeen, 23 St Machar Drive, Aberdeen AB24 3UU, UK. ⁵Commonwealth Scientific and Industrial Research Organization, St Lucia, QLD 4067, Australia. ⁶Ludwig-Maximilian University, Munich, Germany. ⁷Sino-French Institute for Earth System Science, College of Urban and Environmental Sciences, Peking University, 100871 Beijing, China. ⁸International Center for Climate and Global Change Research and School of Forestry and Wildlife Sciences, Auburn University, Auburn, AL, USA. ⁹State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Northwest A&F University, 712100 Yangling, Shaanxi, China. ✉email: changjf@zju.edu.cn</footnote>
Grasslands are managed worldwide to support livestock production, with widely contrasting practices and intensity. Today, grass forage accounts for nearly half of the global intake of livestock¹, although the proportion of grass-fed animals is decreasing as more use is made of crop-based animal feed². Over the last century, grassland management has intensified across the world, with the number of domestic ruminants increasing dramatically—from 1.4 billion head to 3.4 billion³. This historical increase of domestic livestock was, however, preceded in the late nineteenth century by a massive extirpation of wild grazers that were hunted, killed by diseases, or confined by expanding agricultural lands⁴. The biomass of wild mammals today is an order of magnitude lower than that of livestock⁵. Changes in management intensity and the reduction of wild animal numbers alter the greenhouse gas (GHG) emissions and sinks of the world's grasslands, and, together, determine the net radiative forcing (RF) of grasslands on climate change. The intensification of grassland management increases CH₄ emissions from livestock and decreases those from wild grazers. As management becomes more intensive, N₂O emissions, from the acceleration of nitrogen turnover and the fertilization of pasturelands, also increase. However, soil organic carbon stocks can increase or decrease alongside more intensive management, through complex interactions between grazing, soil carbon inputs and decomposition processes. Meta-analyses and literature reviews indicate that land-use⁶ and land-management changes⁷⁻¹⁰ have significant impacts on grassland soil carbon that differ between regions and climate zones⁶,⁹. Site-level data suggest that soil carbon stocks can increase, as long as grazing remains light to moderate⁸,¹¹⁻¹³, but decrease under overgrazing, causing degradation¹⁴,¹⁵.

Assessing the historical GHG budgets of the different grassland biomes and their drivers is of key importance for understanding the trade-offs between grassland services on food security and climate mitigation and how current grassland management could be improved to meet low-warming climate targets. Global scale GHG emissions and sinks have been assessed for all the terrestrial ecosystems¹⁶, the agricultural sector¹⁷ and croplands¹⁸, but grasslands have been ignored, mainly because dealing with pasture management differences, and how they change with time, was deemed too complicated a problem. A few regional-scale

estimates of the GHG budget of specific grassland types have been published¹⁹⁻²¹, but an assessment of the full GHG balance of the world's grasslands, separating the direct human effect of intensified management and the indirect human effects of climate change, is still lacking.

Here, we address this research gap by quantifying the changes in carbon storage and GHG fluxes in natural and managed grasslands over the years 1750–2012. We consider the following processes: soil organic matter and plant productivity changes driven by historical shifts in livestock and wild grazers; fire and its interaction with grazing; soil carbon losses by water erosion; land-use change emissions related to grassland (pasture created by deforestation and the conversion of grassland to cropland); CH₄ emissions from animals and excreta; N₂O emissions from animal excreta, manure and mineral fertilizer applications; and atmospheric nitrogen deposition (Fig. 1). We use the spatially explicit, process-based ecosystem model ORCHIDEE-GM, which includes parameterizations of all these processes (Supplementary Methods 1). Although the model has previously been evaluated at regional and global scales²²,²³, here, we provide additional evaluation for the effects of overgrazing, leading to degradation and soil organic carbon loss, and the depth distribution of soil carbon incorporation. The model is run over the whole globe for the historical period 1860 to 2012 at a spatial resolution of 0.5 × 0.5°. It is forced by variable CO₂, changing climate, nitrogen deposition, variable areas of pasture, and changes of animal density and per capita grass intake²² (see Supplementary Table 1 for the simulation protocol and the time span of the input data). In the following, GHG fluxes are expressed in CO₂ equivalent (CO₂e) based on the global warming potentials of each gas on a 100-year time horizon²⁴ to better combine the effects of emissions of different gases. The climate effect (i.e., RF) due to historical emissions and sinks of the three GHGs and to albedo changes related to grassland is calculated using a compact Earth system model OSCAR v3.1 (ref. ²⁵) that accounts for the different lifetimes of the GHGs (see ‘Methods’).

## Results and discussion

**The GHG budget of grasslands.** Estimates of the grassland fluxes of CO₂, CH₄ and N₂O for the last five decades from ORCHIDEE-GM were evaluated against site observations,

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Illustration of the processes and the greenhouse gas fluxes that are accounted for in this study. The box with an arrow in the top left of the figure presents the drivers that used as input to the ORCHIDEE-GM v3.2. The processes and the greenhouse gas fluxes that are related to grassland are shown as red and blue arrows, respectively. The red-to-blue layout at the bottom of the figure indicates management intensity from high (red; intensively managed grassland) to low (blue; sparsely grazed and natural grassland).</figcaption></figure>
and regional GHG inventories (Supplementary Discussion 1 and Supplementary Table 2). We also compared our simulated GHG budgets with estimates from regional grassland models. The non-CO₂ emissions estimated here were found to be consistent with all five inventories we compiled¹,²,⁶⁻²⁹. We compare simulated carbon budgets, that is, soil carbon changes, with published estimates at national/regional scale and with data from long-term ecological sites known as LTER (Supplementary Table 2). Observed soil carbon change and CO₂ budgets are uncertain and spatially variable. The ORCHIDEE-GM model captures the sign of the carbon budget in 32 out of 42 studies (76%) and its results are within the 1-sigma uncertainty of 22 studies (52%). We found that the model results are within the 1-sigma uncertainty of literature estimates for a number of regions covering 42% of the global grassland area and 86% of the global mean grassland carbon budget during the period of 1980–2012, as simulated by the ORCHIDEE-GM. Model results are entirely outside of the 1-sigma uncertainty of literature estimates in regions covering only 0.3% of the global grassland area and 3% of the grassland carbon budget. This means our model tends to capture the grassland carbon budget in the important regions for a global assessment. The results from this comparison give us some confidence in using this model for long-term historical simulations. We also compared the vertical profile of the simulated grassland soil carbon change against those from a meta-analysis based on stable carbon isotope signature observations of soil profiles³⁰ (see Supplementary Discussion 2). The results indicate that the modelled soil carbon age distribution and the recent soil carbon incorporation (since 1965) between the topsoil (0–30 cm) and subsoil (30–100 cm) are comparable with those from isotope observations of soil profiles³⁰. This gives us some confidence in the values of residence time of soil organic carbon and carbon sequestration in the model. Further, a detailed analysis of how the model simulates regional losses of soil carbon and productivity in response to overgrazing was performed by comparing with the observed loss of productivity and soil carbon in site studies¹¹,³¹,³² and by comparing the simulated overgrazing patterns with maps of grassland degradation derived from satellite data³³ and regional surveys³⁴ (see Supplementary Discussion 3). At the global scale, we compared, for information, the net carbon balance of grass-dominated grid cells in the ORCHIDEE-GM with the results of generic dynamic global vegetation models (TRENDY project) used in annual global carbon budget assessments. Note that the other global models do not incorporate detailed descriptions of changes in livestock, management practices and their feedbacks on ecosystem biogeochemistry (Fig. 1)³⁵. We further compared the ORCHIDEE-GM net grassland-atmosphere CO₂ flux against a top-down atmospheric inversion carbon sink³⁵ after subtraction of forest³⁶ and cropland CO₂ budgets¹⁸. The results show that our estimate fell within the range of this observation-based large-scale constraint. Finally, net primary productivity (NPP), the trend of NPP and whole ecosystem carbon turnover in the model were compared against satellite-based (GIMMS NDVI³³) and observation-driven data sets respectively (MODIS NPP³⁷,³⁸, GIMMS NPP³⁹ and three soil carbon stock data sets at 1 m depth⁴⁰). The results show that it is unlikely that our model overestimates the carbon sink of grasslands, given that: (i) it does not systematically overestimate NPP (i.e., not too high a carbon input) and (ii) simulates faster whole ecosystem carbon turnover than observation-driven estimates (i.e., not too low a carbon output; Supplementary Discussion 1).

Regarding the model's capability of simulating overgrazed grasslands, and the response of productivity and soil organic carbon to different grazing intensities (Supplementary Discussion 3), we found that ORCHIDEE-GM captured the decreasing leaf area index (LAI) trend due to overgrazing in the western United States, southern Brazil, Argentina and Australia³³, but underestimated overgrazing effects in central Asia and Sub-Saharan Africa.

The historical model simulation also captures the areas classified as overgrazed in surveys from the 1980s in the Global Assessment of Soil Degradation (GLASOD) database³⁴. The historical simulation showed an increase in overgrazed areas since the 1860s with the largest expansion of overgrazed grassland occurring in Africa. The modelled responses of soil organic carbon to different grazing intensities are comparable with those from meta-analysis based on in situ observations¹¹,³¹,³². The above analysis gives us some confidence in using this model to simulate the global responses of grassland soil carbon and productivity to management.

The evolution of GHG fluxes during the last two and a half centuries (Fig. 2) shows that emissions of CH₄ and N₂O from grasslands have increased by a factor of 2.5 since 1750. This increase is mainly due to the intensification of management, with increased livestock numbers, greater turnover of manure and, in some regions, mineral nitrogen fertilizer addition. The largest historical emissions from deforestation to pasture occurred after the 1930s, following the large expansion of grazing land in Latin America⁴¹. Yet emissions from deforestation to pasture seem to show a decreasing trend after the 1970s in Latin America and East and Southeast Asia (Supplementary Fig. 1). In North America, Europe, Russia and South Asia, historical emissions from the conversion of grassland to cropland dominate land-use change emissions attributed to grassland. In contrast, grasslands globally have persistently absorbed CO₂ from the atmosphere, resulting in an increase in soil carbon storage (Fig. 2).

Decadal fluctuations have occurred in the GHG balance of grasslands, with values ranging from a net GHG source of 0.6 ± 1.3 Gt CO₂e yr⁻¹ in the 1970s [all results are presented as mean ± standard deviation unless otherwise specified] to a source of 1.8 ± 0.7 Gt CO₂e yr⁻¹ during the 2000s (Fig. 2). Decadal variability of climate and of the area of deforestation are the two most important drivers that result in switching the sign of the decadal GHG balance of grasslands. In contrast, increasing atmospheric CO₂ concentration, nitrogen deposition and management practices²² have slow trends that result in multi-decadal changes of GHG budgets of the same sign.

<figure><img src="image_4.png" /></figure>

**Fig. 2** The decadal greenhouse gas (GHG) fluxes of global grassland during the period 1750–2009. Light and dark blue bars represent CO₂ fluxes from managed and sparsely grazed grassland, respectively; orange and yellow represent CH₄ fluxes from managed (domestic livestock) and sparsely grazed (wild grazers) grassland, respectively; light green and dark green represent N₂O fluxes from managed and sparsely grazed grassland, respectively; and pink and purple represents land-use change emissions related to grassland from deforestation to pasture and from conversion of grassland to cropland, respectively. Black dots and their error bars indicate net total GHG balance and its 1-sigma uncertainty. Red squares and their error bars indicate the anthropogenic GHG balance after subtracting pre-industrial GHG fluxes. Negative values indicate GHG sinks and positive values indicate GHG sources.
The increase of CH₄ and N₂O emissions from grasslands, shown in Fig. 2, is the sum of increased emissions from livestock²⁶ (Supplementary Fig. 2) and decreased emissions from the reduced numbers of wild grazers (Supplementary Methods 2 and Supplementary Fig. 3). In our simulation we found that hunting, disease and reduced habitats due to expanding agriculture caused a decrease in CH₄ and N₂O emissions from wild grazers, from 0.7 ± 0.1 Gt CO₂e yr⁻¹ in 1750 to 0.3 ± 0.07 Gt CO₂e yr⁻¹ in 1900. Today, CH₄ and N₂O emissions are largely dominated by grass-fed livestock, and they amount to 1.6 ± 0.3 and 0.8 ± 0.3 Gt CO₂e yr⁻¹, respectively in 2005. These fluxes comprise 53 ± 9% and 42 ± 14%, respectively, of the emissions from the full livestock sector, that is, when including both grass-fed animals and those fed on crop-products⁴².

The net carbon sink in grasslands worldwide intensified over the last century (Fig. 2), mainly driven by North America, Europe and Russia (Supplementary Figs. 1 and 4). These increasing soil carbon sinks were due to the interaction between indirect human activities, like rising CO₂ concentration, climate change (e.g., warming at high latitudes leading to higher LAI⁴³ and grassland productivity⁴⁴), atmospheric nitrogen deposition, and direct human activities like recent decreases of livestock numbers and pasture abandonment in Europe and Russia. Sparsely grazed and natural grasslands account for 80% of the total cumulative carbon sink of the world's grasslands, and explain most of the current global sink (Fig. 2). Carbon emissions from deforestation to pasture are mostly contributed by South America, and East and Southeast Asia (70% and 21%, respectively), while carbon emissions from conversion of grassland to cropland are dominated by North America (39%), Europe (36%) and South Asia (21%). Over the last decade, managed grasslands are found to be a net GHG source of 2.0 ± 0.4 Gt CO₂e yr⁻¹ (land-use change emissions were excluded here to enable comparisons with the cropland estimates in ref. ¹⁸), that is, in mean value, similar to croplands which are a large GHG source of 2.0 ± 2.2 Gt CO₂e yr⁻¹ (ref. ¹⁸).

To disentangle the effects of management and climate change drivers on the carbon balance of grasslands, we conducted a series of simulations for the period 1860–2012 with one factor fixed at a time to enable us to attribute the contribution of different drivers (Supplementary Discussion 4). The increase of CO₂ by 105 ppm since 1860 caused a net sink of CO₂ (accounting for 85% of the net sink from all drivers including land-use change averaged for the period 1860–2012; Supplementary Fig. 5). The contribution of climate change is a carbon sink in some regions, but a carbon source in others, although, globally, it is a net sink (22%). The contribution of N deposition is comparable, being a global net sink (24% of the total net sink from all drivers). Human-caused increases of livestock densities, with larger sized animals over time, and extirpation of wild grazers, account for 39% of the global net sink including all drivers. In contrast, land-use change-related direct human activities associated with the conversion of forest to pasture and of grassland to cropland is responsible for a large net loss of carbon to the atmosphere, which offsets 141% of the net sink. A residual sink of 69% is due to nonlinear interactions of all these drivers. Net soil carbon removals induced by water erosion over grasslands resulted in a small source of 0.03 [0.03–0.06] Gt CO₂e yr⁻¹ in the 1860s, growing to a significant source of 0.23 [0.22–0.28] Gt CO₂e yr⁻¹ in the 1990s (Supplementary Discussion 5). These net removals comprise water erosion impacts on carbon removal and soil compensatory sinks from continuing litter input concurrent with less soil respiration due to topsoil carbon removal. Regarding fire regime changes resulting from intensified grassland management, the results of our simulation show that, altogether, management effects contributed 35% and 60% of the overall simulated downward trend of burned area during the period 2000–2012 in northern Africa and South America, respectively, through decreased fuel load as animal

density increases (Supplementary Discussion 6). In northern Africa, this effect is larger than the contribution from expanding cropland from which fire is assumed to be excluded⁴⁵.

Though various direct human activities and indirect effects are accounted for in our estimates, some other processes that are not considered may cause further uncertainties: (1) we account for fluxes from mineral soils only: carbon losses from peatland drainage for pasture are not considered. Peatlands drained for agriculture and forestry are substantial GHG sources, especially in the past few decades⁴⁶. However, to our knowledge, historical emissions from pasture expansion on organic soils are not available; (2) though fire processes are simulated specifically, the formation of pyrogenic carbon is not considered in the model. Pyrogenic carbon, as the recalcitrant by-products of fire, can be stored in terrestrial ecosystem for centuries to millennia. However, it should be noted that the amount of pyrogenic carbon created from non-woody biomass burning (e.g., grasses burning) is small⁴⁷ due to the relatively complete burning of non-woody components; (3) ozone's impacts on the grassland carbon cycle are not considered. Elevated ozone (O₃) concentrations may cause leaf damage, reduce biomass growth, and could contribute to changes in species composition⁴⁸. However, a global data set and functions that quantify its effect on grassland productivity and soil carbon dynamics is not available.

Figure 3 shows the spatial distribution of the GHG balance of grasslands and its trend over the past three decades. Net GHG sinks are located in temperate North America and Eurasian grasslands, especially in the tundra and dry grassland regions that have sparse livestock populations, and small CH₄ and N₂O emissions but significant carbon uptake (Supplementary Fig. 6). Conversely, net GHG sources are found in regions of intensive management with high livestock densities. Regional trends also have contrasts. On the one hand, grasslands in the central United States, eastern Europe and Russia show a decreasing trend in GHG emissions (Fig. 3b) coincident with decreasing livestock populations (Supplementary Fig. 2). In these regions, decreased livestock numbers have resulted in smaller CH₄ and N₂O emissions and enhanced soil carbon storage (Fig. 3c). On the other hand, China, Mongolia, eastern Africa and southern Brazil show positive trends of GHG emissions due to their rapidly increasing ruminant populations and, in the case of Brazil, the expansion of pasture over carbon-rich ecosystems (Fig. 3c).

**The net climate effect of grasslands.** The GHG balance of grasslands presented in the previous section contains pre-industrial fluxes and perturbations by human activities including both direct effects from land use and management changes and indirect effects from climate change, increased CO₂ and anthropogenic atmospheric nitrogen deposition. After subtracting pre-industrial emissions (see 'Methods'), the anthropogenic climate change effect of the grassland biome is found to be neutral (Fig. 2), that is, an RF of 12 ± 105 mW m⁻² as calculated by the OSCAR Earth system model (Fig. 4a) (Supplementary Methods 2). In spite of this, the historical increase in livestock alone caused a substantial warming (147 ± 27 mW m⁻² by CH₄ and N₂O emissions) partly balanced by a cooling from the reduced number of wild grazers –47 ± 11 mW m⁻² (CH₄ and N₂O; Supplementary Table 3). The importance of the negative RF component from lost wild grazers has been overlooked in previous studies. Deforestation to pasture and the conversion of grassland to cropland caused net warming of 81 ± 26 and 27 ± 9 mW m⁻², respectively. But the most intriguing result is that the cooling effect of carbon sinks (–194 ± 99 mW m⁻²) nearly offsets the warming effect of land-use change emissions and CH₄ and N₂O sources (209 ± 39 mW m⁻²), when their RF is compared (Supplementary Table 3). We show below that anthropogenic carbon sinks are mainly located in sparsely grazed and natural grasslands, whereas CO₂ and non-CO₂ sources prevail in managed grasslands.
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Global distribution of grassland greenhouse gas (GHG) balance during the period 1981-2010. **a** Mean value, **b** its trend, and **c** the trends of CO₂ fluxes versus non-CO₂ fluxes. Negative values in **a**, the total GHG fluxes, indicate GHG sinks and positive values indicate GHG sources. Negative values in **b**, the trend of total GHG fluxes, indicate enhanced GHG sinks or decreased GHG emissions, while positive values indicate increased GHG emissions or attenuated GHG sinks.</figcaption></figure>

Regionally, grasslands in North America, Russia and Oceania have a small negative RF (Fig. 5) owing to increased carbon storage in soils (Supplementary Fig. 1). The grasslands of Sub-Saharan Africa also have a small net negative RF from the reduced number of wild grazers (reduced CH₄ and N₂O emissions in sparsely grazed and natural grasslands). Grasslands in South Asia have a positive RF due to a decreased albedo from pasture reduction, grassland conversion to cropland and livestock increase. The largest positive RF was ascribed to grasslands in Latin America, due to increases in livestock CH₄ emissions, N₂O from soils and deforestation-induced CO₂ losses (Supplementary Fig. 1). Grassland in other regions has small positive or negative RF resulting from carbon sinks nearly compensating non-CO₂ emissions.

**The cancelling climate effects of managed and sparsely grazed grasslands.** Globally, the positive RF from managed grasslands, including emissions from land-use change related to grasslands, equals 205 ± 48 mW m⁻² (Fig. 4a). This management-induced RF counter-balances the negative RF from sparsely grazed grasslands (–193 ± 80 mW m⁻²) (see ‘Methods’). This result illustrates the important climate cooling service provided by sparsely grazed areas, and the key role played by quickly increasing livestock numbers, associated with more CH₄ and N₂O emissions, in determining the contemporary net RF of the grassland biome.

Managed grasslands in all regions, except Russia, have a net warming effect on climate (Fig. 5). The direct human effect of grassland management on climate was analysed by comparison
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Radiative forcing induced by direct and indirect human activities on grasslands worldwide a in the year 2012 and b for the period 1900–2012.</figcaption></figure>

The 'Total' column/line is the total contribution of grassland to the global radiative forcing (RF), accounting for the land-use change induced albedo change related to grassland, land-use change emission related to grassland (CO₂ E_LUC) and the three grassland greenhouse gas (GHG) fluxes (CO₂ soil sink, CH₄ and N₂O). In a, the RF of each component is further divided for managed and sparsely grazed grasslands, and is shown as an empty bar outlined in the same colour with (managed grasslands) and without (sparsely grazed grasslands) hatching, respectively. The stacked bar for CO₂ E_LUC includes, from bottom to top, the effect from deforestation to pasture and conversion of grassland to cropland. The stacked bar for CH₄ includes, from bottom to top, the direct effect from CH₄, and its associated effects through tropospheric O₃ and lifetime change induced by ozone precursors, and stratospheric H₂O. The stacked bar for N₂O includes, from bottom to top, the direct effect from N₂O, and its associated effects through stratospheric O₃. The uncertainties shown in a are 1-sigma standard deviation. The coloured shaded areas in b are the 1-sigma confidence interval derived from our uncertainty assessment (see 'Methods').

with a baseline simulation in which all grasslands would remain natural without livestock changes or fertilization (see 'Methods' and Supplementary Table 1). The results indicate that management intensification caused 9% less soil carbon storage since the pre-industrial period (set to before 1750), because grazing and mowing reduced the carbon input to soils (i.e., overall carbon outputs from grassland ecosystem increased) more than they increased input via the stimulation of plant productivity through the generation of new leaves (Supplementary Figs. 7 and 8). In particular, regions with overgrazing have a decreased productivity and either smaller soil carbon gains or carbon losses (see Supplementary Discussion 3).

In contrast, regions with sparsely grazed grasslands provided an important climate cooling service, as in North America, Russia and Oceania (Fig. 5), through indirect climate change effects that enhanced carbon sinks. This result makes it clear that countries should assess not only the GHG budgets of their managed pastures (such as in the current national GHG reporting rules to the UNFCCC) but also the sinks/sources of sparsely grazed rangelands, steppes, tundra and natural grassland systems. Full GHG reporting for the entire territory of each country could facilitate the assessment of progress towards the goals of the Paris Agreement⁴⁹ and better link national GHG budgets to the observed growth rates of GHGs in the atmosphere. Moving to a complete estimate of all GHG fluxes from grasslands calls for global models to incorporate a realistic representation of management practices, including livestock grazing, degradation, erosion and fires.

Today, on a global basis, managed grasslands provide half of the feed needed for livestock, with the other half coming from crop-based feed and residues. To fully assess the future climate impact of the livestock sector, or more generally, of the entire agricultural sector, the crop-feed-related GHG balance or emissions since 1750 would need to be included: a requirement that could be built upon data already compiled for the last decade¹⁸.

Despite increased livestock and management intensification, intriguingly, grasslands worldwide are found to exert no warming effect on climate. This conclusion is thanks to the presence of intensified carbon sinks, especially over sparsely grazed grasslands, which mainly result from the increased productivity of grasslands exposed to increased CO₂ and nitrogen deposition. Maintaining this

cooling service will be important for meeting the Paris Agreement targets. However, the recent trend (Fig. 4b) reveals that global grasslands are transitioning from a net cooling towards a net warming effect on climate. This trend can be attributed to the recent grassland management intensification for livestock production and the conversion of tropical forest to pasture. In the context of low-warming climate targets, the mitigating or amplifying role of grasslands will depend on: (1) future changes in grass-fed livestock numbers; (2) the stability of accumulated soil carbon in grasslands; and (3) whether carbon storage can be further increased over time, or if it will saturate, as observed in long-term experiments⁵⁰,⁵¹. For example, a projected continual growth in bovine and ovine meat production and consumption in the near future, especially an accelerated growth of milk demand and production⁵², would result in a continual increase in grass biomass demand and thus increases in CH₄ and N₂O emissions. The future grassland net CO₂ flux will depend on changes in environmental drivers such as warming, increased CO₂ concentration, management intensity⁵³, nitrogen addition, including fertilizer and atmospheric nutrient deposition, future land-use⁶,⁵⁴, and fire regimes. To prevent further warming contributions from managed grasslands, sustainable management through optimized grazing, pasture improvement and restoration of degraded pasture⁷,⁵⁵,⁵⁶ are critical to enhance grassland soil carbon sequestration. Halting deforestation to pasture is also a priority in countries like Brazil where moderate intensification could avoid extensive cattle production and clearing of forests. Technologies reducing livestock CH₄ emissions, such as improving feeding practices, specific agents and dietary additives, and longer-term management change and animal breeding⁵⁷,⁵⁸, and practices increasing nitrogen use efficiency and thus reducing N₂O emissions from fertilization and manure (e.g., as collected in ref. ⁵⁹), are also key for climate mitigation and should be explored in integrated assessments.

## Methods

**Methods overview.** This study assesses CO₂, CH₄ and N₂O fluxes in grassland ecosystems (including ecosystem-level fluxes from grasslands, and farm-level CH₄ and N₂O fluxes from management activities, i.e., grass forage consumption and manure management). The process-based land surface model ORCHIDEE-GM v3.2 (refs. ²²,⁶⁰) was used to simulate all carbon cycling flux and carbon pool components within ecosystems, between farm and ecosystem (i.e., grass forage from mown
<figure><img src="image_7.png" /><figcaption>**Fig. 5** The contribution of different regions to the global radiative forcing of grasslands. The ‘Total’ column is the total contributions of grassland to the global radiative forcing (RF), accounting for the land-use change induced albedo change related to grassland, land-use change emission ($CO_2 E_{LUC}$) from deforestation to pasture and conversion of grassland to cropland, and the three grassland GHG fluxes ($CO_2$ soil sink, $CH_4$ and $N_2O$). The RF of each component is further divided for managed and sparsely grazed grasslands, and is shown as an empty bar outlined in the same colour with (managed grasslands) and without (sparsely grazed grasslands) hatching, respectively. The uncertainties are 1-sigma standard deviation. The vertical-scale of the regional contribution has been adjusted so that the effects can be easily seen.</figcaption></figure>

grasslands and managed manure applied to pasture), and $CH_4$ and $N_2O$ emissions produced at the ecosystem and farm level. ORCHIDEE-GM v3.2 simulated the soil carbon changes of 11 soil layers. Litter input (structural and metabolic, above and below ground, is discretized in each model layer and decomposition is simulated from the fraction of three pools with different turnovers in each layer, modified by the temperature and soil moisture content of each layer which are simulated every 30 min by the land surface model. The model distinguishes explicitly between sparsely grazed and managed grassland, and calculates emissions from domestic livestock²² and wild grazers⁴. Populations of wild grazers since 1750 were reconstructed on a 0.5° global grid (Supplementary Methods 1). Pre-industrial $CH_4$ and $N_2O$ fluxes are subtracted from historical values to quantify post-industrial anthropogenic fluxes. To calculate the RF induced by each separate anthropogenic flux (nine different regions, three gases with a separation between land-use change emissions and soil carbon storage for $CO_2$, managed pastures and sparsely grazed grasslands), the compact Earth system model OSCAR v3.1 (refs. ²⁵,⁶¹,⁶²) was used. This model is specifically designed to attribute regional and global climate change to multiple forcings and feedbacks⁶³⁻⁶⁵. Here, OSCAR is driven by GHG flux distributions from ORCHIDEE-GM. These drivers are integrated to disentangle the effect of each climate forcer for managed grasslands (mown, grazed by livestock, fertilized) versus sparsely grazed and natural grasslands for nine regions of the globe. To calculate the contribution of grassland to the land-use change induced albedo change and associated RF, we used a diagnostic method that combines present-day satellite data of surface albedo with historical land-cover maps⁶⁶ (see also ref. ⁶¹).

**The estimation of GHG fluxes and their uncertainty.** ORCHIDEE-GM v3.2 (refs. ²²,⁶⁰), which includes specific parameterizations of grassland management⁶⁷,

was applied over the entire globe, at a spatial resolution of 0.5° × 0.5°, with variable $CO_2$, climate, atmospheric nitrogen deposition, land-use change and management intensity (including manure and mineral fertilizer application, grazing density and managed area)²². Historical changes of grass biomass consumption (biomass intake during grazing and harvested grass forage for confined livestock, including landless production systems¹) were reconstructed from historical data¹,²², and prescribed as input data to the model²². A detailed description of the model, inputs and simulation setups can be found in Supplementary Methods 1. We defined grassland areas as a dynamic variable including all lands that have been grasslands (natural or human managed) since 1750. For the period of 1860–2012, grassland areas were defined by prescribed historical land-cover maps⁶⁸ (Supplementary Table 4), and managed (i.e., mown, grazed and fertilized) and sparsely grazed grassland areas were reconstructed using the method detailed in ref. ²² (Supplementary Fig. 2 and Supplementary Methods 1). We assumed no changes in grassland areas between 1750 and 1860 due to lack of data. Grassland area changed annually across history in the historical land-cover maps⁶⁸ from 51.7 million km² in 1860 to 49.7 million km² in 2012. Managed grasslands are defined as all grasslands used to feed livestock (i.e. grazed and/or mown for grass forage, fertilized in some regions), and sparsely grazed grasslands are natural grasslands not affected by livestock but may be sparsely grazed by wild grazers. A detailed discussion on the grassland and managed grassland areas, and the distributions used here, can be found in Supplementary Discussion 7. Fire⁶⁹,⁷⁰ and interactions between grazing and fire are explicitly represented in ORCHIDEE-GM (Supplementary Discussion 6). Though the model does not simulate erosion explicitly, we account for grassland carbon fluxes due to water erosion by using estimates from a model framework combining the carbon cycle of ORCHIDEE-MICT⁷¹ (a similar version that ORCHIDEE-GM v3.2 is built on) and a modified version of the Revised Universal Soil Loss Equation (RUSLE) model that is applicable at a large spatial
scale<sup>72</sup> (Supplementary Discussion 5). It should be noted that only mineral soils were considered in this study, since the model does not simulate peatlands (organic soils). Grassland GHG fluxes due to land-use change from sparsely grazed grasslands to managed grasslands are explicitly simulated by the model, given reconstructed historical management intensity changes. For each year, we estimated the grassland GHG balance, including new grasslands converted from other land-cover types (like forests and croplands), and land-use change emissions (CO₂ E<sub>LUC</sub>) related to grassland. CH₄ and N₂O emissions following the conversion of natural ecosystems to pasture and the increased livestock density, caused by livestock digestion and excreta deposition, were accounted for through the changes in grass biomass consumption.

In ORCHIDEE-GM, carbon stocks and associated carbon fluxes for land-cover transitions involving grasslands are attributed to the original source land cover type if it is a grassland type. For instance, the carbon balance of a cropland created from a grassland is reported under grassland. The carbon balance of a forest cleared to pasture is also reported under grassland. However, in the latter case, legacy emissions from harvested wood products during deforestation to pasture were not included in the land use-induced component of the grassland GHG balance.

In our estimates of the net grassland GHG balance induced by land use change, with the above definitions, carbon losses due to deforestation to pasture, which is dominated by biomass loss<sup>73</sup> rather than changes in soil organic carbon<sup>6</sup>, and due to conversion of grassland to cropland were separated with factorial simulations (CO₂ E<sub>LUC</sub>) as in the DGVMs that were run for the global annual carbon budget assessments from the Global Carbon Project<sup>35</sup>. Specifically, three sets of simulations were performed with ORCHIDEE-GM. The first set was forced with a time-invariant pre-industrial land-cover distribution. The second simulation was forced with historical land-use transitions from forest to grassland (with invariant cropland areas as in 1860, in order not to include cropland land-use change GHG emissions), climate, atmospheric CO₂ concentration, and nitrogen deposition. The third simulation was forced with historical land-use transitions from grassland to cropland (with invariant forest areas as in 1860), climate, atmospheric CO₂ concentration, and nitrogen deposition. Land-use change emissions from deforestation to pasture were diagnosed by the difference between the first and the second simulations. Land-use change emissions from conversion of grassland to cropland were diagnosed by the difference between the first and the third simulations. We note that this approach to quantify land cover change emissions involving grasslands includes a foregone sink named the lost atmospheric sink capacity<sup>74,75</sup> and thus may give higher values in the recent period than emissions calculated by land use bookkeeping models (e.g. ref. 76).

The grassland GHG balance at ecosystem and farm scale includes CO₂, CH₄ and N₂O fluxes of grassland ecosystems and off-site CO₂, CH₄ and N₂O emissions from the digestion of harvested forage by livestock and manure management. The CO₂ and CH₄ fluxes of the grassland ecosystem were simulated by ORCHIDEE-GM. Tier 1 or Tier 2 methods of the 2006 IPCC guidelines for national greenhouse gas inventories<sup>77</sup> were used for calculating N₂O emissions from managed soils, CH₄ emissions from digestion of harvested forage, and N₂O and CH₄ emissions from manure management. The activity data required for estimating N₂O and CH₄ emissions through IPCC Tier 1 or Tier 2 methods were either simulated by ORCHIDEE-GM v3.2 (like grazed biomass, excreta nitrogen deposited as urine and dung by grazing animals on grassland, harvested grass forage, and excreta nitrogen from housed livestock consuming harvested grass forage), or model inputs (like manure and fertilizer application, and atmospheric nitrogen deposition). It should be noted that potential CH₄ and N₂O uptake by grassland soils is not taken into account here. The grassland GHG balance (unit: Gt CO₂e yr⁻¹) was calculated as

$$GHG = F_{CO_2-C} \times \frac{44}{12} + E_{LUC-grass} \times \frac{44}{12} + F_{CH_4-C} \times \frac{16}{12} + GWP_{CH_4} + F_{N_2-O-N} \times \frac{44}{28} + GWP_{N_2-O} \quad (1)$$

where the molecular weight conversion factors 44/12, 16/12 and 44/28 are used to convert the mass of CO₂-C, CH₄-C and N₂O-N into CO₂, CH₄ and N₂O, respectively. GWP<sub>CH₄</sub> (Gt CO₂e per Gt CH₄) and GWP<sub>N₂O</sub> (Gt CO₂e per Gt N₂O) are global warming potential constants used to give the integrated RF of CH₄ and N₂O in terms of their CO₂ equivalent. Here, we adopt the GWP for a time horizon of 100 years, with inclusion of climate-carbon feedbacks (i.e., GWP100 AR5 + OSCAR + IRF and REs updates in Table 2 of ref. 24) of 34 and 267 for CH₄ and N₂O, respectively. F<sub>CO₂-C</sub>, F<sub>CH₄-C</sub> and F<sub>N₂O-N</sub> are annual fluxes of CO₂, CH₄ and N₂O, respectively (unit: Pg C yr⁻¹ or Pg N yr⁻¹), between the grassland ecosystem/farm system and the atmosphere. E<sub>LUC-grass</sub> (unit: Pg C yr⁻¹) is the land-use change emissions from deforestation to pasture. All the GHG flux calculations in this study were made on a grid-cell basis. A detailed description of the estimation of the GHG fluxes can be found in Supplementary Methods 1.

The uncertainty in the regional and global GHG fluxes and their components was assessed using a Monte Carlo approach (n = 10,000). Given the scarcity of evaluation data for changes in grassland soil carbon, we estimated the uncertainty of the grassland carbon budget by comparing the simulated carbon balance with a compilation of observations at site and regional level obtained from an extensive literature survey (Supplementary Discussion 1). The relative uncertainty was estimated as the relative standard error of the mean distance between the simulated carbon budgets and the observations collated from the literature, with model output sampled across the grid-points or regions with available data. In total, 22 observation-based studies (Supplementary Table 2) are used for estimating the relative uncertainty, while regional estimates from ecosystem models, or estimates giving only a range of observations (i.e., without mean values), are excluded. We derived a relative uncertainty of 46%. It is noteworthy that the relative uncertainty derived from model

structural uncertainty from 16 TRENDYv6 models (42%; Supplementary Table 5) is as large as that estimated by model-data comparison. Here, the relative uncertainty for land-use change emissions related to grassland (31%) is derived directly from the estimates of the TRENDYv6 models<sup>35</sup> (standard deviation divided by multi-model ensemble mean; Supplementary Table 5). For each element of the Monte Carlo approach, we randomly selected regional/global CO₂ fluxes (e.g., grassland CO₂ fluxes and land-use change emissions) from their distributions. For CH₄ and N₂O, we assessed regional and global CH₄ and N₂O emission uncertainties due to emission-related factors (Supplementary Methods 1) through a Monte Carlo approach (n = 10,000). For each element of the Monte Carlo approach, we randomly selected emission-related factors from their distributions. Following a suggestion by ref. 77 (Vol 4, Chapter 3), a normal distribution was assumed for factors with symmetric uncertainty, and a lognormal distribution was assumed for factors with asymmetric ranges of uncertainty.

**The estimation of anthropogenic GHG fluxes.** Following the IPCC<sup>78</sup>, we defined the year 2012 as the present day and the year 1750 to be the end of pre-industrial times. Given the existence of domestic livestock and wild grazers in the pre-industrial era, CH₄ and N₂O emissions from both domestic livestock and wild grazers before 1750 do not contribute to the global RF. They were estimated in this study to emit 13.8 ± 2.7 Tg C yr⁻¹ of CH₄ and 0.9 ± 0.3 Tg N yr⁻¹ of N₂O, for a total of 1.0 ± 0.2 Gt CO₂e yr⁻¹ (converted using GWP100 with climate-carbon feedback<sup>24</sup>), while global grasslands were assumed to be CO₂ neutral during pre-industrial times. Therefore, for the purpose of attributing the contribution of grassland's GHG fluxes to RF, we defined the human-induced (anthropogenic) GHG fluxes of grassland by subtracting pre-industrial CH₄ and N₂O emissions from the historical/contemporary GHG fluxes. The net CO₂ fluxes from sparsely grazed grasslands were also considered to be human-induced, because they are driven by indirect human drivers.

The anthropogenic GHG fluxes from grassland, for the period 1750–2012, were used by OSCAR v3.1 to calculate the grassland contribution to RF. The total (i.e. anthropogenic + pre-industrial) GHG fluxes for the period 2060–2010 were either directly simulated by ORCHIDEE-GM, or estimated through IPCC Tier 1 or Tier 2 methods using active data simulated by ORCHIDEE-GM or as model inputs. They were extrapolated over the period 1750–1859 following a set of assumptions. First, grassland is assumed to be CO₂ neutral with net CO₂ fluxes of 0 Gt CO₂ yr⁻¹. Second, CH₄ and N₂O emissions from managed grasslands (i.e., mainly from domestic ruminants, and to a lesser extent N₂O emissions from atmospheric nitrogen deposition) were varied following the human population variation in each region, with the emissions of 1860 used as the reference value. Third, CH₄ and N₂O emissions from sparsely grazed grasslands (i.e., mainly from wild grazers, and to a small extent from N₂O emissions from atmospheric nitrogen deposition) were varied, following the reconstructed regional wild grazer population, over the period 1800–1859 (see Supplementary Methods 1 for details), but kept constant, at the 1800 value, over the period 1750–1799. The human population density since 1750 was retrieved from the History Database of the Global Environment (HYDE) compiled by the Netherlands Environmental Assessment Agency<sup>41</sup>. The HYDE gridded data are available for the beginning of each decade from 1750 until 1860. Annual data were linearly interpolated within each decade.

**Attribution of RF induced by anthropogenic GHG fluxes from grasslands.** In all the simulations by OSCAR v3.1 (ref. 25), climate change was forced and based on observations. The model was fed by reconstructions of past climate change following HadCRUT4 (ref. 79 and updates in 2018) for global temperature, HadISST1 (ref. 80 and updates in 2018) for sea surface temperature, and CRU<sup>81</sup> for regional land temperature and precipitation. The oldest 30-year data period was assumed to correspond to the pre-industrial climate. The anthropogenic GHG fluxes from grassland since 1750 were used as input for OSCAR v3.1. The uncertainty in those simulations was assessed via a Monte Carlo approach (n = 10,000) to the OSCAR simulations, in which biophysical parameters of the model were randomly drawn (see ref. 61) and associated to every single set of anthropogenic grassland GHG fluxes of the first Monte Carlo ensemble. Because OSCAR naturally produces a range of uncertainty larger than we considered to be reasonable for the three GHGs simulated in this study, we weighted our Monte Carlo ensemble using existing emission inventory data sets as constraints (see Supplementary Methods 2 for details).

The attribution is made to nine regions and eight sectors. Sectors are defined as emissions of CO₂, CH₄ and N₂O from managed versus sparsely grazed grasslands, and of CO₂ from land-use change emissions of deforestation to pasture and of conversion of grassland to cropland. The attribution protocol requires that, in addition to a control experiment, 73 simulations are made, one for each combination of region and sector, plus one for the rest (i.e. all the remaining emitting sectors unrelated to grasslands). In each of these simulations, the same emissions as in the control experiment are prescribed, except that a small fraction of 1% of the region's sectoral emissions are removed. The difference between the control simulation and each of those factorial simulations, once normalized by the sum of the 73 differences, provides the contribution of each region and sector. This approach to attribution is called the normalized marginal method, and it is described further by ref. 65 and references therein.

The RF induced by albedo changes in grassland-related land cover (RF<sup>LCC</sup><sub>grass</sub>) was assessed using the same method as used by OSCAR v3.1 (ref. 25) derived from ref. 66. The RF<sup>LCC</sup><sub>grass</sub> was calculated, using a time step of one year, following annual
changes of grassland area across history in the historical land-cover maps⁶⁸. Here, grassland-related land-cover change was defined as both grasslands converted from (grassland increase) and to (grassland decrease) other land-cover groups (i.e., other biomes). In this case, any RF between two biomes should be equally attributed to the biome with decreased area and the biome with increased area. To ensure that the sum of biome-specific RF (i.e., $\sum_b RF_b^{LCC}$) was equal to the RF from albedo change induced by all land-cover changes, and therefore to avoid double counting, we divided the final $RF_b^{LCC}$ grass by 2, as in Supplementary Methods 2, Equation (45). The uncertainty of $RF_b^{LCC}$ grass was assessed as the standard deviation of the 18 estimates of $RF_b^{LCC}$ grass derived from the combinations of different data sets: three land-cover climatologies, three radiation flux climatologies, and two albedo climatologies (see Supplementary Methods 2 for detail). The total current global RF induced by land-use change induced albedo estimated in this study (-0.09 ± 0.02 W m⁻²; Supplementary Fig. 9) is within the uncertainty range of the IPCC estimate of -0.15 ± 0.10 W m⁻² (ref. ⁷⁸). Values of $RF_b^{LCC}$ grass for managed and sparsely grazed grasslands (see above for the definitions) were also separated, and the same was done for the anthropogenic GHG fluxes and their RF values.

## Data availability

The data that support the findings of this study are available at the public Data Repository of the International Institute of Applied Systems Analysis (IIASA DARE; https://dare.iiasa.ac.at/110; DOI: 10.22022/esm/11-2020.110).

## Code availability

The source code for ORCHIDEE-GM v3.2 is available at https://doi.org/10.14768/20190319001.1. The source code of OSCAR is available at https://github.com/tgasser/OSCAR (https://doi.org/10.5281/zenodo.4282113). The code used to generate all the results of this study and the specific scripts used to execute the attribution with OSCAR are available from the corresponding author upon request.

Received: 4 June 2020; Accepted: 1 December 2020;
Published online: 05 January 2021

## References

1. Herrero, M. et al. Biomass use, production, feed efficiencies, and greenhouse gas emissions from global livestock systems. *Proc. Natl Acad. Sci. USA* **110**, 20888–20893 (2013).

2. Hasha, G. *Livestock Feeding and Feed Imports in the European Union: A Decade of Change* (US Department of Agriculture, Economic Research Service, 2002).

3. HYDE database, https://themasites.pbl.nl/tridion/en/themasites/hyde/landusedata/livestock/index-2.html (2019).

4. Smith, F. A. et al. Exploring the influence of ancient and historic megaherbivore extirpations on the global methane budget. *Proc. Natl Acad. Sci. USA* **113**, 874–879 (2016).

5. Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on Earth. *Proc. Natl Acad. Sci. USA* **115**, 6506–6511 (2018).

6. Li, W. et al. Temporal response of soil organic carbon after grassland-related land-use change. *Glob. Change Biol.* https://doi.org/10.1111/gcb.14328 (2018).

7. Conant, R. T., Cerri, C. E., Osborne, B. B. & Paustian, K. Grassland management impacts on soil carbon stocks: a new synthesis. *Ecol. Appl.* **27**, 662–668 (2017).

8. Conant, R. T., Paustian, K. & Elliott, E. T. Grassland management and conversion into grassland: effects on soil carbon. *Ecol. Appl.* **11**, 343–355 (2001).

9. McSherry, M. E. & Ritchie, M. E. Effects of grazing on grassland soil carbon: a global review. *Glob. Change Biol.* **19**, 1347–1357 (2013).

10. Wang, X. et al. Grazing improves C and N cycling in the Northern Great Plains: a meta-analysis. *Sci. Rep.* **6**, 33190 (2016).

11. Abdalla, M. et al. Critical review of the impacts of grazing intensity on soil organic carbon storage and other soil quality indicators in extensively managed grasslands. *Agricult. Ecosyst. Environ.* **253**, 62–81 (2018).

12. Hewins, D. B. et al. Grazing and climate effects on soil organic carbon concentration and particle-size association in northern grasslands. *Sci. Rep.* **8**, 1336 (2018).

13. Soussana, J. F., Tallec, T. & Blanfort, V. Mitigating the greenhouse gas balance of ruminant production systems through carbon sequestration in grasslands. *Animal* **4**, 334–350 (2010).

14. Ganjegunte, G. K. et al. Soil organic carbon composition in a northern mixed-grass prairie. *Soil Sci. Soc. Am. J.* **69**, 1746–1756 (2005).

15. Khresat, S., Rawajfih, Z. & Mohammad, M. Land degradation in north-western Jordan: causes and processes. *J. Arid Environ.* **39**, 623–629 (1998).

16. Tian, H. Q. et al. The terrestrial biosphere as a net source of greenhouse gases to the atmosphere. *Nature* **531**, 225 (2016).

17. Tubiello, F. N. et al. The FAOSTAT database of greenhouse gas emissions from agriculture. *Environ. Res. Lett.* **8**, 015009 (2013).

18. Carlson, K. M. et al. Greenhouse gas emissions intensity of global croplands. *Nat. Clim. Change* **7**, 63–68 (2017).

19. Chang, J. et al. The greenhouse gas balance of European grasslands. *Glob. Change Biol.* **21**, 3748–3761 (2015).

20. Schulze, E. D. et al. Importance of methane and nitrous oxide for Europe's terrestrial greenhouse-gas balance. *Nat. Geosci.* **2**, 842–850 (2009).

21. Soussana, J. F. et al. Full accounting of the greenhouse gas (CO₂, N₂O, CH₄) budget of nine European grassland sites. *Agricult. Ecosyst. Environ.* **121**, 121–134 (2007).

22. Chang, J. et al. Combining livestock production information in a process-based vegetation model to reconstruct the history of grassland management. *Biogeosciences* **13**, 3757–3776 (2016).

23. Chang, J. et al. Modeled changes in potential grassland productivity and in grass-fed ruminant livestock density in Europe over 1961–2010. *PLoS ONE* **10**, e0127554 (2015).

24. Gasser, T. et al. Accounting for the climate-carbon feedback in emission metrics. *Earth Syst. Dyn.* **8**, 235 (2017).

25. Gasser, T. et al. Historical CO₂ emissions from land use and land cover change and their uncertainty. *Biogeosciences* **17**, 4075–4101 (2020).

26. FAOSTAT. *Online Statistical Service* (Food and Agriculture Organization (FAO), 2017).

27. Dangal, S. R. S. et al. Global nitrous oxide emissions from pasturelands and rangelands: magnitude, spatio-temporal patterns and attribution. *Glob. Biogeochem. Cycles* **33**, https://doi.org/10.1029/2018GB006091 (2019).

28. EDGARv4.3.2. http://edgar.jrc.ec.europa.eu/, last accessed 2018 (2018).

29. Oenema, O., Velthof, G., Yamulki, S. & Jarvis, S. Nitrous oxide emissions from grazed grassland. *Soil Use Manag.* **13**, 288–295 (1997).

30. Balesdent, J. et al. Atmosphere-soil carbon transfer as a function of soil depth. *Nature* **559**, 599–602 (2018).

31. Conant, R. T. & Paustian, K. Potential soil carbon sequestration in overgrazed grassland ecosystems. *Glob. Biogeochem. Cycles* **16**, https://doi.org/10.1029/2001gb001661 (2002).

32. Dlamini, P., Chivenge, P. & Chaplot, V. Overgrazing decreases soil organic carbon stocks the most under dry climates and low soil pH: a meta-analysis shows. *Agricult. Ecosyst. Environ.* **221**, 258–269 (2016).

33. Tucker, C. J. et al. An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. *Int. J. Remote Sens.* **26**, 4485–4498 (2005).

34. Oldeman, L., Hakkeling, R., Sombroek, W. & Batjes, N. Global assessment of human-induced soil degradation (GLASOD). *World Map of the Status of Human-Induced Soil Degradation* (Winand Staring Centre-ISSSFAO-ITC, Wageningen, Netherlands, 1991).

35. Le Quéré, C. et al. Global Carbon Budget 2017. *Earth Syst. Sci. Data* **10**, 405–448 (2018).

36. Pan, Y. et al. A large and persistent carbon sink in the world's forests. *Science* **333**, 988–993 (2011).

37. Zhao, M. & Running, S. W. Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. *Science* **329**, 940–943 (2010).

38. Zhao, M. S., Heinsch, F. A., Nemani, R. R. & Running, S. W. Improvements of the MODIS terrestrial gross and net primary production global data set. *Remote Sens. Environ.* **95**, 164–176 (2005).

39. Smith, W. K. et al. Large divergence of satellite and Earth system model estimates of global terrestrial CO₂ fertilization. *Nat. Clim. Change* **6**, 306–310 (2016).

40. Tifafi, M., Guenet, B. & Hatté, C. Large differences in global and regional total soil carbon stock estimates based on SoilGrids, HWSD, and NCSCD: intercomparison and evaluation based on field data from USA, England, Wales, and France. *Glob. Biogeochem. Cycles* **32**, 42–56 (2018).

41. Klein Goldewijk, K., Beusen, A., Doelman, J. & Stehfest, E. Anthropogenic land use estimates for the Holocene-HYDE 3.2. *Earth Syst. Sci. Data* **9**, 927 (2017).

42. Gerber, P. J. et al. *Tackling Climate Change Through Livestock: A Global Assessment of Emissions and Mitigation Opportunities* (Food and Agriculture Organization of the United Nations (FAO), 2013).

43. Zhu, Z. et al. Greening of the Earth and its drivers. *Nat. Clim. Change* **6**, 791–795 (2016).

44. Zhang, K. et al. Satellite-based model detection of recent climate-driven changes in northern high-latitude vegetation productivity. *J. Geophys. Res. Biogeosci.* **113**, https://doi.org/10.1029/2007JG000621 (2008).

45. Andela, N. & van der Werf, G. R. Recent trends in African fires driven by cropland expansion and El Niño to La Niña transition. *Nat. Clim. Change* **4**, 791–795 (2014).

46. Leifeld, J., Wüst-Galley, C. & Page, S. Intact and managed peatland soils as a source and sink of GHGs from 1850 to 2100. *Nat. Clim. Change* **9**, 945–947 (2019).
47. Jones, M. W., Santín, C., van der Werf, G. R. & Doerr, S. H. Global fire emissions buffered by the production of pyrogenic carbon. *Nat. Geosci.* **12**, 742–747 (2019).

48. ASHMORE, M. R. Assessing the future global impacts of ozone on vegetation. *Plant, Cell Environ.* **28**, 949–964 (2005).

49. Grassi, G. et al. Reconciling global-model estimates and country reporting of anthropogenic forest CO2 sinks. *Nat. Clim. Change* **8**, 914–920 (2018).

50. Johnston, A. E., Poulton, P. R. & Coleman, K. Soil organic matter: its importance in sustainable agriculture and carbon dioxide fluxes. *Adv. Agron.* **101**, 1–57 (2009).

51. Smith, P. Do grasslands act as a perpetual sink for carbon? *Glob. Change Biol.* **20**, 2708–2711 (2014).

52. Alexandratos, N. & Bruinsma, J. *World Agriculture Towards 2030/2050: The 2012 Revision*. ESA working paper (FAO, Rome, 2012).

53. Chang, J. F. et al. Effect of climate change, CO₂ trends, nitrogen addition, and land-cover and management intensity changes on the carbon balance of European grasslands. *Glob. Change Biol.* **22**, 338–350 (2016).

54. Pugh, T. A. M. et al. Simulated carbon emissions from land-use change are substantially enhanced by accounting for agricultural management. *Environ. Res. Lett.* **10**, 124008 (2015).

55. Herrero, M. et al. Greenhouse gas mitigation potentials in the livestock sector. *Nat. Clim. Change* **6**, 452–461 (2016).

56. Smith, P. et al. Greenhouse gas mitigation in agriculture. *Philos. Trans. R. Soc. B Biol. Sci.* **363**, 789–813 (2008).

57. Monteny, G.-J., Bannink, A. & Chadwick, D. Greenhouse gas abatement strategies for animal husbandry. *Agricult. Ecosyst. Environ.* **112**, 163–170 (2006).

58. Soliva, C. R., Takahashi, J. & Kreuzer, M. (eds) Greenhouse gases and animal agriculture: An update. International Congress Series No. 1293, pp. 1–378. (Elsevier, The Netherlands, 2006).

59. Smith, P. et al. Agriculture. In Climate Change 2007: Mitigation. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Metz, B. et al. eds), (Cambridge University Press, Cambridge, United Kingdom and New York, NY, 2007).

60. Chang, J. et al. Future productivity and phenology changes in European grasslands for different warming levels: implications for grassland management and carbon balance. *Carbon Balance Manag.* **12**, 11 (2017).

61. Gasser, T. et al. The compact Earth system model OSCAR v2.2: description and first results. *Geosci. Model Dev.* **10**, 271–319 (2017).

62. Gasser, T. et al. Path-dependent reductions in CO2 emission budgets caused by permafrost carbon release. *Nat. Geosci.* **11**, 830–835 (2018).

63. Ciais, P. et al. Attributing the increase in atmospheric CO₂ to emitters and absorbers. *Nat. Clim. Change* **3**, 926 (2013).

64. Gasser, T. *Attribution régionalisée des causes anthropiques du changement climatique* (Université Pierre et Marie Curie-Paris VI, 2014).

65. Li, B. G. et al. The contribution of China's emissions to global climate forcing. *Nature* **531**, 357–361 (2016).

66. Bright, R. M. & Kvalevag, M. M. Technical Note: Evaluating a simple parameterization of radiative shortwave forcing from surface albedo change. *Atmos. Chem. Phys.* **13**, 11169–11174 (2013).

67. Chang, J. F. et al. Incorporating grassland management in ORCHIDEE: model description and evaluation at 11 eddy-covariance sites in Europe. *Geosci. Model Dev.* **6**, 2165–2181 (2013).

68. Peng, S. et al. Sensitivity of land use change emission estimates to historical land use and land cover mapping. *Glob. Biogeochem. Cycles* **31**, 626–643 (2017).

69. Yue, C., Ciais, P., Cadule, P., Thonicke, K. & van Leeuwen, T. Modelling the role of fires in the terrestrial carbon balance by incorporating SPITFIRE into the global vegetation model ORCHIDEE-Part 2: carbon emissions and the role of fires in the global carbon balance. *Geosci. Model Dev.* **8**, 1321–1338 (2015).

70. Yue, C. et al. Modelling the role of fires in the terrestrial carbon balance by incorporating SPITFIRE into the global vegetation model ORCHIDEE—Part 1: simulating historical global burned area and fire regimes. *Geoscientific Model Dev.* **7**, 2747–2767 (2014).

71. Guimberteau, M. et al. ORCHIDEE-MICT (v8.4.1), a land surface model for the high latitudes: model description and validation. *Geosci. Model Dev.* **11**, 121–163 (2018).

72. Naipal, V. et al. Global soil organic carbon removal by water erosion under climate change and land use change during AD 1850–2005. *Biogeosciences* **15**, 4459–4480 (2018).

73. Houghton, R. & Nassikas, A. A. Global and regional fluxes of carbon from land use and land cover change 1850–2015. *Glob. Biogeochem. Cycles* **31**, 456–472 (2017).

74. Pongratz, J., Reick, C. H., Raddatz, T. & Claussen, M. Effects of anthropogenic land cover change on the carbon cycle of the last millennium. *Glob. Biogeochem. Cycles* **23**, https://doi.org/10.1029/2009gb003488 (2009).

75. Stocker, B. D. & Joos, F. Quantifying differences in land use emission estimates implied by definition discrepancies. *Earth Syst. Dyn.* **6**, 731–744 (2015).

76. Houghton, R. A. et al. Carbon emissions from land use and land-cover change. *Biogeosciences* **9**, 5125–5142 (2012).

77. IPCC. 2006 IPCC Guidelines for National Greenhouse Gas Inventories (Institute for Global Environmental Strategies Hayama, Japan, 2006).

78. Myhre, G. et al. Anthropogenic and natural radiative forcing. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Stocker, T. F. et al. eds) (Cambridge University Press, Cambridge, 2013).

79. Morice, C. P., Kennedy, J. J., Rayner, N. A. & Jones, P. D. Quantifying uncertainties in global and regional temperature change using an ensemble of observational estimates: The HadCRUT4 data set. *J. Geophys. Res. Atmospheres* **117**, D08101 (2012).

80. Rayner, N. et al. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. *J. Geophys. Res. Atmos.* **108**, 4407.

81. Harris, I., Jones, P. D., Osborn, T. J. & Lister, D. H. Updated high-resolution grids of monthly climatic observations—the CRU TS3.10 Dataset. *Int. J. Climatol.* **34**, 623–642 (2014).

## Acknowledgements

J.C., D.S.G., T.G., and P.C. acknowledge support from the ERC Synergy grant ERC-2013-SyG-610028 IMBALANCE-P. The input of P.S. contributes to the Belmont Forum/FACCE-JPI funded DEVIL project (NE/M021327/1), the UK NERC funded projects U-GRASS (NE/M016900/1) and Soils-RGRREAT (NE/P019455/1) and the EU funded H2020-EU.3.5.1. project VERIFY (No 776810). B.G. received support from the French Agence Nationale de la Recherche (ANR) Convergence Lab Changement climatique et usage des terres (CLAND). This work was performed using HPC resources from GENCI-TGCC (Grant 2017-A0030106328). We thank Dr. William Kolby Smith for sharing NPP data. We thank all people of the TRENDY modelling group who provided TRENDYv6 data used in this study. We thank John Gash for his effort on English language editing.

## Author contributions

J.C. and P.C. designed the study. J.C. and T.G. carried out modelling and analysis. J.C., P.C., and T.G. wrote most of the manuscript. S.P. provided reconstructed land cover change data. P.S., M.H., P.H., M.O., B.G., D.S.G., W.L., V.N., S.P., C.Q., H.T., N.V., C.Y., and D.Z. contributed significantly to the final revisions of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-20406-7.

Correspondence and requests for materials should be addressed to J.C.

Peer review information *Nature Communications* thanks Carol Barford, Longlong Xia, and other, anonymous, reviewers for their contributions to the peer review of this work.

Reprints and permission information is available at http://www.nature.com/reprints

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021