
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39300072
52508
10.1038/s41467-024-52508-x
Article
Temperature overshoot responses to ambitious forestation in an Earth System Model
http://orcid.org/0000-0002-4448-2967
Moustakis Yiannis yiannis.moustakis@geographie.uni-muenchen.de

1
http://orcid.org/0000-0002-3603-4283
Nützel Tobias 1
Wey Hao-Wei 2
Bao Wenkai 1
Pongratz Julia 13
1 https://ror.org/05591te55 grid.5252.0 0000 0004 1936 973X Ludwig-Maximilians-Universität in Munich, Munich, Germany
2 https://ror.org/02h2x0161 grid.15649.3f 0000 0000 9056 9663 GEOMAR Helmholtz Centre for Ocean Research Kiel, Kiel, Germany
3 https://ror.org/05esem239 grid.450268.d 0000 0001 0721 4552 Max Planck Institute for Meteorology, Hamburg, Germany
19 9 2024
19 9 2024
2024
15 82357 12 2023
7 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
Despite the increasing relevance of temperature overshoot and the rather ambitious country pledges on Afforestation/Reforestation globally, the mitigation potential and the Earth system responses to large-scale non-idealized Afforestation/Reforestation patterns under a high overshoot scenario remain elusive. Here, we develop an ambitious Afforestation/Reforestation scenario by harnessing 1259 Integrated Assessment Model scenarios, restoration potential maps, and biodiversity constraints, reaching 595 Mha by 2060 and 935 Mha by 2100. We then force the Max Planck Institute’s Earth System Model with this scenario which yields a reduction of peak temperature by 0.08 oC, end-of-century temperature by 0.2 oC, and overshoot duration by 13 years. Afforestation/Reforestation in the range of country pledges globally could thus constitute a useful mitigation tool in overshoot scenarios in addition to fossil fuel emission reductions, but socio-ecological implications need to be scrutinized to avoid severe side effects.

Forestation can be a valuable tool for mitigating climate change. Using an Earth System Model, Moustakis et al. show that ambitious deployment in the range of country pledges can mitigate global temperature under an overshoot pathway by 0.2 oC, but highlight the associated socioeconomic tradeoffs.

Subject terms

Climate-change mitigation
Climate and Earth system modelling
Carbon cycle
https://doi.org/10.13039/501100002347 Bundesministerium für Bildung und Forschung (Federal Ministry of Education and Research) 01LS2101A Moustakis Yiannis issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Reaching net-zero greenhouse gas emissions as required by the Paris Agreement’s goals will require the deployment of Carbon Dioxide Removal (CDR) methods, compensating at least for hard-to-abate emissions1–3. Pursuing a more ambitious 1.5 oC (compared to pre-industrial levels) goal implies a more limited remaining emission budget, and thus rapid decarbonization, major societal transformations, and large-scale CDR application1,4. Thus, CDR deployment has been at the spotlight of scientific research recently5,6. Future scenarios generated by Integrated Assessment Models (IAMs) that keep to 1.5 oC by 2100 typically exhibit excess positive emissions early on which are then sharply reduced, leading to a temperature overshoot (i.e., period during which mean global temperature exceeds intended warming level, here 1.5 oC, before returning to it)7,8. This delay of early action in IAMs emerges due to the high upfront investments and near-term mitigation costs required to limit the overshoot9,10, but caution is needed since the normalization of the overshoot idea could possibly facilitate a political flexibility associated with a lack of action11. Overshoot is becoming increasingly more relevant due to incredible country pledges12,13, and concerns of political feasibility and delay of the required climate action14,15. Despite that, uncertainty still remains with regards to overshoot dynamics such as the peak temperature and overshoot duration and the associated Earth system responses and adaptation needs, which should thus be further studied9,16.

Afforestation/Reforestation (AR), which has been practiced for decades, can be a useful mitigation tool17 and combined with forest management practices constitutes virtually all (99.9%) of currently applied CDR18. In the literature a wide range of possible AR sequestration potentials has been reported, varying depending on the model used, land availability, CO2 levels, and assumptions regarding the capacity of different biomes in sequestering CO2 under a changing climate7,19,20. Apart from the potential of durably sequestering carbon (C) and mitigating warming, uncertainty remains regarding the permanence of C storage21, possible impacts of large-scale AR on local climates22, ecosystem biodiversity19, food security23, and the associated societal risks7. Still, recently the Land-Gap Report24 has estimated that long-term and net-zero country-level mitigation pledges amount to quite ambitious land-based CDR commitments, which would require 633 Mha of AR globally by 2060, in addition to 551 Mha for the restoration of degraded ecosystems. This estimation has recently been updated reaching likely up to 490 Mha of AR and up to 570 Mha of restoration25. Even though such levels of AR might seem rather high and their feasibility and associated socioeconomic risks have been questioned24, it is only through decisive action and ambitious policies that we can diverge from the high-risk climate path that we are currently following13. Hence, the mitigation potential and impacts of large-scale AR in the range of country pledges should be further assessed.

The complexity of assessing ambitious AR deployment emerges since both the changes in atmospheric CO2 following CDR application and the biogeophysical effects via altered moisture and surface energy fluxes can strongly affect the Earth system26. To capture these effects and feedbacks and thus get more reliable estimates of AR potentials and impacts, fully coupled Earth System Models (ESMs) need to be employed, which solve for the exchange of mass and energy between the land surface, the atmosphere and the ocean27. The assessment of biogeophysical effects has mostly been done for land use change, land management, and related land-cover changes (short “land use” in the following) in general. Studies investigated local or remote temperature26,28, precipitation29, and circulation changes30, caused by changes in surface roughness and albedo31, moisture recycling and net radiation22, as well as their interplay with the biogeochemical effect (cooling) of the associated decrease in atmospheric CO232.

Even though significant progress has been made in isolating and investigating such effects, studies in many cases employ fully idealized land use scenarios including massive (de)forestation globally, and cannot be thus used to confidently assess the effects of more realistic patterns33–35. This is especially the case with the multitude of results reported by deforestation studies, since it has recently been demonstrated that the responses of local climate to forest loss and gain are asymmetric36. Modeling efforts with a focus on assessing the impacts and mitigation potentials of more realistic AR scenarios have been rather limited. Some have been based on the available land use change scenarios associated with the Representative Concentration Pathway 4.5 (RCP4.5) from the Coupled Model Intercomparison Project 5 (CMIP5)37–39. In RCP4.5, however, the total ~800Mha forest area increase by 2100 also includes a forest expansion that the models simulate in response to climate change, and typically, by experimental design (comparison is only possible between various land use scenarios), AR can be assessed only together with land transitions other than AR, including avoided deforestation, and thus it is not possible to isolate the effect of CDR39. Recently, Matthews et al.40 examined the mitigation potential of forest regrowth to 1920 levels by 2056, followed by a release of sequestered carbon during the second half of the century, and demonstrated that high AR levels, even when of temporary nature, can mitigate global temperature increase under a low overshoot scenario. However, Matthews et al.40 employed an intermediate complexity climate model, which does not fully represent all the relevant processes, and an idealized spatiotemporal AR pattern. In the ongoing Land Use Model Intercomparison Project (LUMIP)41 the land use pattern of the Shared Socioeconomic Pathway (SSP) SSP1-2.6 is considered. Even though this was intended to include substantial AR based on the marker IAM (IMAGE) scenario42, translation to ESMs has failed to reproduce that signal yielding moderate AR levels43, thus falling behind the ambitious amounts suggested by country pledges globally.

Overall, uncertainty still remains regarding the Earth system responses under temperature overshoot pathways16. Reported results show a significant variability of responses among ESMs, whose statistical significance, as well as the signal emergence from internal natural variability have not been robustly assessed, due to the lack of multiple (within and across model) ensemble members16,44. Yet, in most AR studies high future emission scenarios are used34,39. Importantly, such approaches suffer from the paradox of applying ambitious AR -implying strong coordinated climate action-, while following high-end emission scenarios -implying delayed and insufficient climate action. The normalization of the idea of coupling largely unmitigated warming scenarios with ambitious AR could connotate the establishment of AR as an alternative to emission reductions, and needs to be treated with caution45. From a biophysical perspective, high-end scenarios induce an enhanced CO2 fertilization effect, thus strengthening vegetation carbon uptake, land surface fluxes, mitigation potentials and biogeochemical cooling34,38,39,41. Typically, they also yield intensified disturbances (e.g., wildfires) which could compromise the permanence of carbon stored in vegetation46. Evidently, the role of AR under scenarios with stronger emission reductions yielding a weaker vegetation fertilization, and hence less biogeochemical cooling, as well as lower disturbance rates, thus suffering less from CDR permanence overestimation, remains unclear.

To fill this gap, we investigate to what extent an ambitious constrained AR scenario in the range of total country pledges globally can mitigate a temperature overshoot emission scenario and alter its dynamics (duration, peak and end-of-century temperature), and whether biogeochemically-induced cooling can emerge from internal variability and compensate for any possible AR-induced biogeophysical warming. We first develop an AR scenario in the range of country pledges globally reaching 595Mha of AR in 2060 and 935 Mha by 2100, and then run fully coupled simulations with the Max Planck Institute’s Earth System Model (MPI-ESM). Simulations are performed under a high-overshoot SSP5-3.4os emissions scenario which follows SSP5-8.5 until 2040, followed by rapid decarbonization, and reaching net-zero around 2070 and net-negative thereafter. To take into account the technoeconomic considerations, biogeophysical potentials, and biodiversity concerns, we harvest the recently published AR6 Scenarios Database (AR6-SD)47, restoration potential maps, and maps of human influence. Beyond these considerations we also discuss in detail the possible socioeconomic barriers and tradeoffs of such ambitious AR implementation.

Our results show that ambitious AR in the range of country pledges globally can mitigate a temperature overshoot scenario by reducing the peak and end-of-century temperature by 0.08 oC and 0.2 oC respectively, as well as the duration of the overshoot by 13 years. Overall, biogeochemical cooling due to the reduction of CO2 dominates, and at least compensates for any biogeophysically-induced warming.

Results

Afforestation/Reforestation pattern

Following AR6-SD estimates (Fig. 1, “Methods”), AR is applied at the expense of grazing land globally, while cropland areas and natural grassy and shrubby biomes remain unchanged. By considering restoration potential maps and biodiversity proxies (Fig. 2, “Methods”), the resulting AR scenario (Fig. 3) exhibits a strong increase in forest area over the 21st century, with slightly higher dynamics in the first half, reaching 595 (935) Mha of forest area increase by 2060 (2100) over current (2015) levels. This occurs at the expense of heavily managed pastures, which decrease by 472 (575) Mha and more lightly managed rangelands, which decrease by 123 (360) Mha. Latin America (256 Mha in 2100) and Africa (215 Mha) emerge as hotspots of forest area increase, followed by the Organization for Economic Cooperation and Development (OECD90) developed countries (201 Mha), Asia (162 Mha), and Eastern European countries and the reforming economies of the former Soviet Union (101 Mha). At the country level, the United States (142 Mha), Brazil (113 Mha), China (110 Mha), Russian Federation (64 Mha), and Argentina (24 Mha) carry out in total 48% of AR by 2100 (Supplementary Fig. 1). The scenario includes reforesting 514 Mha in 2100 of previously deforested areas globally, corresponding to 60% of total historical deforestation (851 Mha during 1850–2014 in MPI-ESM following the Land Use Harmonization 248 historical pattern) (Supplementary Fig. 2). The remaining 421 Mha are afforestation, which is employed by priority over the less biodiverse rangelands (Fig. 2 & Supplementary Fig. 2).Fig. 1 Percentiles of the AR6 scenario database.

The figures present the different percentiles of global (a) forest cover, (b) grazing land, and (c) cropland change in Mha compared to 2015. Each line does not correspond to a specific scenario, but rather represents a different percentile of area change for each year, estimated by pooling all AR6 Scenario Database scenarios together. Highlighted is the 90th percentile for forest area change and its complementary threshold, the 10th percentile for cropland and grazing land change. The probability density plots for 2040, 2060, 2080, and 2100 indicative of the spread of the area change distribution are also shown in light gray. The country pledges global estimate presented in the Land-Gap report24, and its 2023 update25 are shown with the green bar.

Fig. 2 Overview of the AR scenario development.

Given that total rangelands (a) can have different levels of management and grazing intensity and hence biodiversity richness, we use the Very Low Human Influence (VLHI) map based on Riggio et al.127 (b) to filter out grazing land that can be considered closer to a pristine state. Combined with pasture, the remaining available grazing land (c) constitutes the total grazing land that is considered for AR (d) in our framework. Available rangeland is further categorized into 4 biodiversity groups (e), based on the Low Human Influence (LHI) map based on Riggio et al.127 (f). The Griscom et al.19 restoration potential map (GRS, (g)) and Atlas of Forest and Landscape Restoration Opportunities126 (ATL, (h)) are further used to guide the spatial pattern of AR, given an annual AR target for each world region. Allocation of the AR target across the gridcells is performed using the iterative process presented in the “Methods” section. The resulting AR pattern by 2100 for the target set in this study is presented in (i).

Fig. 3 Results of the AR scenario development (see “Methods” and Fig. 2).

The panels in (a) show the timeseries of forest area increase (in Mha) and how AR is allocated between pasture and rangelands i) globally and at the regional level for ii) Organization for Economic Cooperation and Development (OECD90) developed countries, (iii) Eastern Europe and reforming economies of the former Soviet Union (EEu & RefSU), (iv) Latin America (LAC), (v) Middle East and Africa (MAF), and (vi) Asia (the 5 world economic regions typically used in IAMs). The country pledges global estimate presented in the Land-Gap report24, and its 2023 update25 are also shown with the green bar in (i). In the right-column maps the resulting change in cover fraction of (b) forest cover (corresponding to Fig. 2i), (c) heavily managed grazing land (Pasture), and (d) lightly managed grazing land (Rangeland) from 2015 to 2100 is shown (other vegetation/land use including croplands remain unaltered).

Mitigation of global temperature

The average global 2 m air temperature (compared to the 1850–1900 MPI-ESM average) and the global carbon fluxes are shown in Fig. 4. Under the reference (REF) simulation (constant land use at 2015 levels, see “Methods”), temperature overshoot lasts for ~65 ± 2 years (mean ± standard deviation estimated from smoothing and bootstrapping—see “Methods”), from 2035 up to 2099. The probability that end-of-century temperature will be within 1.5 °C over pre-industrial levels is 46%, with an average temperature of 1.51 °C (minimum—maximum of 5-year period from ensemble members: 1.27–1.87 °C—see “Methods”). In year 2058 ± 2, a peak temperature of 2.06 °C (1.89–2.29 °C) is reached, which follows the year during which atmospheric carbon reaches its maximum value (~ + 980 GtCO2, 126 ppm in 2057 compared to 2015). Under the AR simulation, temperature overshoot lasts for ~52 ± 2 years, from 2036 up to 2087. The probability that the 1.5 °C target is returned to by the end of the century is 90%, with an average temperature of 1.31 °C (1.02–1.59 °C). In year 2058 ± 2, a peak temperature of 1.98 °C (1.8–2.23 °C) is reached, staying below the Paris Agreement 2 °C target with a probability of 48%. Peak temperature lags ~3 years after the maximum in atmospheric carbon. Comparing the REF and AR simulations suggests that the isolated effect of AR on average global temperature yields peak temperature reduced by ~0.08 oC, overshoot duration by ~13 years, and end-of-century temperature by ~0.2 oC. The impact of AR on global temperature emerges already in ~2052 ± 2 years, when the temperature difference with REF starts becoming statistically significant (Fig. 4a). At that point, AR has reached 495 Mha, atmospheric carbon is lower by ~110 GtCO2 (30 PgC, 14 ppm), and oceanic carbon by ~22 GtCO2 (6 PgC), due to a total increase of ~132 GtCO2 (36 PgC) in the land carbon sink (Fig. 4b). The temperature trajectory under AR is smoother, rather reaching a plateau before it starts leveling off. By examining different percentiles of global average daily temperature a similar behavior is demonstrated (Supplementary Fig. 3).Fig. 4 Earth system responses.

a Average annual global 2 m air temperature (difference compared to pre-industrial era, expressed here as the 1850–1900 average) is plotted for each ensemble member (light lines) for both AR (green color) and REF (blue color) scenarios. The thick lines represent the ensemble mean (10 members) for each scenario. The colored arrows depict the overshoot duration for each scenario. The bright green shaded region indicates that there is no statistically significant difference between the two scenarios, while the purple shading suggests that a statistically significant difference exists (significance estimated as described in “Methods”). b The changes in global ensemble mean total atmospheric (Catmo), land (Cland), and oceanic (Cocean) carbon (GtCO2) compared to 2015 are shown for each scenario. All ensemble members are plotted in light color, but overall variability is weak compared to the temperature variability. c Map showing the ensemble mean difference in carbon (GtCO2 in each gridcell) stored on land between the AR and the REF scenario in 2100. Blue color indicates an increase in land carbon under AR.

CO2 removal

The land C stock divergence of AR from REF is notable within few decades and continues until the end of the century (Fig. 4b). Over gridcells where AR is applied the land sequesters in total ~395 GtCO2 (108 PgC) more compared to REF by 2100, with less carbon being sequestered over the forested rangelands, which are typically more arid (Fig. 4c). This is the net result of ~420 GtCO2 (115 PgC) sink enhancement over 912 Mha of our AR area, and a weakening of the land C sink by ~25 GtCO2 (7 PgC) over the remaining 23 Mha, mostly in the tropics (despite the forest area increase). Overall, ~249 GtCO2 (68 PgC) correspond to the partial reversal of historical deforestation. Over gridcells where AR is not applied, a total weakening of CO2 uptake by ~13 GtCO2 (4 PgC) emerges, most prominently over the Amazon. Therefore, the net land sink increase compared to REF is ~382 GtCO2 (104 PgC) by 2100 in total, reducing the atmospheric CO2 level by ~281 GtCO2 (77 PgC, 36 ppm), and the ocean sink by ~101 GtCO2 (27 PgC).

The rate of change in land C under both REF and AR mostly follows the signal of change in atmospheric C, which is driven by fossil fuel emissions (Supplementary Fig. 4). Under REF, the land C stock transits from being a sink to a source in ~2090, whereas under AR the land C remains a sink, even though with weak and ever decreasing sequestration rates, becoming a source only in 2100. This suggests that under AR increased forest cover coupled with decreased CO2 atmospheric levels compared to REF offer a ~10-year buffer before the land C stock transits to a source behavior, allowing for more C to be sequestered during this period. However, the reverse image is observed when the ocean C sink is examined (Fig. 4 & Supplementary Fig. 4). Here, not only is ocean C uptake consistently weaker under AR compared to REF, due to the lower CO2 concentration as a consequence of the CDR application, but in fact the ocean sink transits to a source behavior under AR in ~2085. Under REF the ocean C sink does not transit to a source behavior, even though by 2100 the rate of increase is almost zero. This demonstrates that the ocean tends to partly offset the efficiency of applied CDR. On average, for every additional 100 GtCO2 (27 PgC) sequestered in land under AR compared to REF 26 GtCO2 (7 PgC) less carbon is taken up by the ocean, leading to a net removal of 74 GtCO2 (20 PgC) from the atmosphere, which suggests a 74% efficiency.

Overall, the applied AR manages to sequester ~4.49 GtCO2/year more over land compared to REF on average, but eventually reduces atmospheric CO2 only by ~3.32 GtCO2/year on average, due to compensation by the ocean CO2 fluxes. By 2100, this corresponds to a ~41 GtCO2 (11 PgC) gain over land compared to REF and ~0.02 °C temperature reduction for every 100 Mha of forest area increase, which also suggests a ~0.05 °C temperature reduction for every 100 GtCO2 (27 PgC) additionally sequestered in land compared to REF. CDR efficiency expressed as temperature reduction per 100 GtCO2 removed from the atmosphere compared to REF increases across time, reaching ~0.07 °C in 2100 (Supplementary Fig. 5). Efficiency expressed in terms of C removal per 100 Mha increases until ~2080 and then plateaus until the end of the century, reaching on average ~30 GtCO2 (8 PgC) (Supplementary Fig. 5).

AR sequestration rate defined as the rate of land C increase when AR and REF are compared shows a sharp increase to 3.9 GtCO2/year in 2030 (27 GtCO2 cumulatively), rising further to 6.3 GtCO2/year in 2050 (121 GtCO2) (Fig. 5). These high rates are sustained and strengthened until ~2065 (7.1 GtCO2/year, 222 GtCO2), before they start levelling off, reaching 1.7 GtCO2/year (382 GtCO2) in 2100 even though at this point the individual land C stocks under REF and AR have both transited to a source behavior. When it comes to the rate of decrease in atmospheric CO2 when AR and REF are compared, a removal rate of 3.39 GtCO2/year is achieved by 2030 (24 GtCO2), 5.3 GtCO2/year by 2050 (103 GtCO2), and 5.6 GtCO2/year by 2065 (185 GtCO2). By 2100, atmospheric CO2 rates reach −0.21 GtCO2/year (281 GtCO2 in total), since additively the ocean and land C sources under AR are stronger than the land C source under REF.Fig. 5 Sequestration potentials of Afforestation/Reforestation.

Sequestration potentials of AR based on this study, previous literature, and AR6-SD are shown, expressed as (a) cumulative carbon sequestration through AR over time (GtCO2) and (b) annual sequestration rate through AR (GtCO2/year). The green line represents the ensemble mean (10 members) additional sequestration by land compared to REF achieved in this study and the dark red line the corresponding removal (positive) from the atmosphere, both smoothed with a Savitzky-Golay122 filter. The 5–95% percentiles (purple shading), as well as the mean (dark purple line) of AR sequestration estimates of AR6-SD are also presented, estimated based on a subset of 300 scenarios for which explicit information on C sequestration from AR is available. The black dots (and the associated minimum-maximum ranges with vertical thin black lines—if applicable) represent (mean) AR sequestration estimates based on the compiled list of 40 previous studies (see “Methods”, and Supplementary Data 1). The most widely cited studies among them19,20,67 and the most recent estimate of Mo et al.75 which are explicitly mentioned in the “Discussion” section, are shown with different color. The technical and economic mitigation potentials suggested by IPCC WGIII Ch.7174 are also shown for comparison. To aid visualization, the studies have been randomly shifted in time by ±2.5 years.

Changes in regional hydroclimate

Gridcell-level results (Fig. 6, Supplementary Fig. 6) show patches of cooling emerging as early as 2030, when the mean temperature pattern under AR is already statistically significantly different compared to REF (see “Methods”, and Supplementary Fig. 7). Regions with significant cooling during 2030–2050 emerge not only for mean, but also for more extreme temperature, mostly over the northern high-latitudes, and over Latin America, and central Africa (Fig. 6, top row). Such cooling becomes progressively more robust and widespread over time, and is already dominant during the period around the peak of global temperature (2050–2070) (Fig. 6, mid row), especially over Africa, Latin America, and the northern high-latitudes. The pattern of mitigation is consistent across the mean, and the high and low temperature extremes, but differs in magnitude over the northern-latitudes where low temperature extremes are reduced more strongly. Patches of warming at the sites of forestation in Eurasia and North America emerge, which are however statistically insignificant. By the end of the century, cooling is dominant globally for mean temperature (Fig. 6, bottom row). Cooling is stronger over Africa, Latin America south of the Amazon rainforest, northern Asia, and North America. A similar pattern is obtained for the extreme high temperature, albeit with a weaker temperature mitigation over the northern high-latitudes. On the contrary, mitigation for the extreme low percentile is stronger over the northern high-latitudes, and a small patch of statistically significant warming emerges in eastern Asia over sites of forestation.Fig. 6 Spatiotemporal pattern of temperature change.

The differences in 10th percentile (low temperature extremes, left), mean (mid), and 90th percentile (high temperature extremes, right) of 2 m air temperature between AR and REF simulations during 2030–2050 (top), the period around peak warming (2050–2070) (mid), and end-of-century (2090–2100) are shown. A negative difference (blue color) indicates that temperature is lower in the AR scenario. Dots indicate regions where the difference is statistically insignificant at the 5% level, estimated with a two-tailed Student’s t-test after correcting for lag-1 temporal autocorrelation123 (see ”Methods”).

Our results further suggest a wetter hydroclimate over forested regions emerging during the mid-century (Fig. 7). This is most clearly evident with increases in relative humidity and cloudiness over sites of forestation, which start emerging already from 2030 onwards, and are especially pronounced and widespread over Latin and North America, Africa, and east Asia. This is accompanied by statistically significant but less widespread patterns of precipitation and evapotranspiration increases evident over Latin and North America, Africa, and east Asia, where the hotspots of changes in cloudiness and relative humidity are reported. Over sites of forestation reduced albedo leads to an increase in surface net radiation, thus resulting in an increase in latent heat fluxes commensurate with the reported increase in evapotranspiration and in a reduction of sensible heat flux (Supplementary Fig. 8). An increase in albedo is reported over the Arctic already from the mid-century, which can be explained by the increase in sea-ice, which is however noisy at the global scale (Supplementary Fig. 9). Surface ocean pH reaches a 0.03 increase under AR compared to REF (Supplementary Fig. 9).Fig. 7 Spatiotemporal pattern of changes in hydroclimate.

From left to right, the differences in mean precipitation (mm/day), evapotranspiration (mm/day), cloud cover fraction, and relative humidity (%) between AR and REF simulations during 2030–2050 (top), the period around peak warming (2050–2070) (mid), and end-of-century (2090–2100) (bottom) are shown. A negative value indicates a reduction in the AR scenario. Dots indicate regions where the difference is statistically insignificant at the 5% level, estimated with a two-tailed Student’s t-test after correcting for lag-1 temporal autocorrelation123 (see “Methods”).

Discussion

Mitigation potential of ambitious Afforestation/Reforestation

Here we go beyond simply using marker IAM scenarios and rather utilize an ensemble of 1259 land use scenarios, which to our knowledge is unprecedented. This allows us to navigate through the more optimistic spectrum of the wide ensemble of opportunity that IAMs can offer49 and design an AR scenario aligned in spirit with the ambitious climate policies that are necessary to overcome climate uncertainty and yield robust results in terms of achieving our climate goals13,50. For the first time, the mitigation potential of ambitious AR in the range of country pledges24 is studied under an overshoot scenario with a fully coupled high-complexity ESM. This is particularly important since delayed climate action15, insufficient Nationally Determined Contributions12, the policy credibility gap13, political feasibility concerns14, and the possible overstatement of CDR potentials in IAMs51, suggest that reaching 1.5 oC with no or limited overshoot is less likely. This increases the possibility and policy-relevance of high temperature overshoot pathways11,12 unless immediate drastic action is taken. Our multi-ensemble member approach allows for a robust statistical treatment of the results, adding useful insight to the uncertainty surrounding the responses of the Earth system to temperature overshoot pathways16.

Even though direct comparisons with previous studies employing different spatiotemporal AR patterns, models, and/or emission pathways cannot be made easily, our findings are in rough agreement with previous studies. Matthews et al.40 reported a 0.02 oC reduction of peak warming and a lower efficiency of 0.03 oC for every 100 GtCO2 sequestered in land, although under the SSP1-1.9 emission pathway and an idealized AR pattern. Sonntag et al.39 reported similar results with MPI-ESM yielding a 0.27 oC end-of-century temperature reduction, as a consequence of an 800 Mha forest area increase concurrent with 100 Mha of avoided deforestation compared to the reference scenario under RCP8.5 emissions. Using a reduced complexity climate model Dooley et al.52 demonstrated an upper limit of 0.25 oC temperature reduction in 2100, but a negligible impact on peak temperature by ecosystem restoration efforts under a low overshoot scenario. However, their scenario yields peak temperatures in 2030–2050, which is too early for the signal to emerge, since time is needed not only for the CDR to be scaled up, but also for carbon to accumulate. This agrees with our results where the signal emerges only in 2052, despite ambitious deployment early on, and only after 132 GtCO2 have been accumulated in land. Under a high-emission scenario Arora and Montenegro53 reported a global temperature reduction of 0.25 oC estimated with an ESM, after increasing forest cover following an idealized pattern by 1010 Mha globally due to an increase of the land C sink by 440 GtCO2, roughly agreeing in magnitude with our estimate. Recently, based on the LUMIP simulation employing the land use pattern SSP1-2.6 under emissions following SSP5-8.5, Loughran et al.43 reported insignificant changes in global temperature across 6 ESMs following a land C sink increase with an intermodel range of 37–220 GtCO2. However, part of this apparent uncertainty may be related to the lack of ensemble members, and the significant differences in forest area increase across models.

The extent to which our findings are model-specific remains to be assessed, by utilizing a multitude of available ESMs, each with different process representation and climate sensitivity under an emission-driven configuration. Both peak and end-of-century temperature under REF are at the lower range of the CMIP6 SSP5-3.4os multi-model ranges (5 models not including MPI-ESM) of 2–4.35 oC and 1.39–3.47 oC, respectively, as recently reported by Asaadi et al.54, suggesting that MPI-ESM is among the models with the stronger cooling behavior as a response to negative emissions. However, direct comparisons cannot be made with confidence, since in our study land use under REF is constant, while it follows the trajectory of SSP5-3.4os in Asaadi et al.54. It should also be noted that the ranges reported therein are not based on multiple ensemble members. Results from the Zero Emissions Commitment Model Intercomparison Porject (ZECMIP)55 have demonstrated that MPI-ESM shows a stronger cooling following the cessation of emissions compared to other models56. Despite the sensitivity of MPI-ESM to negative emissions not having been directly compared against other models so far16,54 and uncertainty remaining, it is the additional sequestration under AR and the consequent feedbacks that determine temperature mitigation compared to REF. Therefore, it is rather the feedback strength within MPI-ESM that would play the key role. In particular, both land carbon gain due to increased atmospheric CO2 concentration as well as land carbon loss due to elevated temperature in MPI-ESM are below the CMIP6-model average57. As a result, MPI-ESM is placed close to the average of CMIP6-models in terms of carbon sequestration when the total effect of elevated CO2 is considered. However, MPI-ESM has also a slightly below-average Transient Climate Response compared to other ESMs58, indicating that the cooling reported here is possibly below anticipated average ESM responses, for given levels of carbon sequestration. All this adds confidence to the robustness of our results, demonstrating the effectiveness of large-scale AR at reducing temperature even under strongly reduced emissions scenarios.

Changes in rainfall can affect terrestrial ecosystem water and carbon fluxes in various ways, depending on the ecosystem type59. De Hertog et al.60 reported on the local and non-local biogeophysical effects of idealized forestation on the water cycle with the MPI-ESM and the Community Earth System Model (CESM) ESMs, with both models reporting increased evapotranspirative fluxes and precipitation over the tropics and subtropics agreeing with our results, albeit with CESM not yielding increased water fluxes over the higher latitudes. A wetter hydroclimate with increases in rainfall due to AR has also been recently reported based on a kilometer-scale model61. However, based on an idealized global forestation scenario with CESM Portmann et al.30 reported that pronounced precipitation and cloud cover changes did not coincide with hotspots of forest area increase, but were rather driven by large-scale circulation adjustments.

When it comes to the obtained temperature pattern at the gridcell level, direct comparisons with previous findings cannot be made, since they significantly differ in the AR spatiotemporal pattern and other land use, and the experimental setups, whereas the number of relevant ESM-based studies is limited. For example, under a SSP5-3.4os scenario where AR is applied instead of bioenergy crops after 2040, Melnikova et al.44 reported patches of strong warming over the mid- and high- latitudes compared to a fixed land use scenario, even though locally at many forested gridcells cooling emerged. However, this was a concentration-driven run where the land C uptake does not result in a reduction of atmospheric CO2 and thus the results reported therein are indicative of the biogeophysical effects only. De Hertog et al.34 suggested that the total biogeophysical effect of AR in ESMs (including MPI-ESM) under idealized AR patterns tends to yield a temperature increase over the mid- and high-latitudes and a cooling over the (sub)tropics. According to our results following an emission-driven approach, this signal is dominated or at least compensated for locally by global biogeochemical cooling caused by the reduction of atmospheric CO2 levels, as the land takes up more CO2. This agrees qualitatively with Sonntag et al.39 who also used MPI-ESM and is noteworthy because here we employ a scenario with stronger emission reductions and hence significantly lower CO2 fertilization rates. However, the extent to which this can be consistent across models remains to be assessed.

An important emerging feature in our simulations is the early rise and continuous buildup of the land C stock, with high productivity reached already within decades, followed by a continuous forest C uptake, as the woody and soil C pools saturate only slowly62. However, despite the continuous increase of forest cover, the signal of AR sequestration is dominated by the emission trajectory that reaches net-negative thus drastically reducing atmospheric CO2 concentrations, and consequently the strength of CO2 fertilization. Nevertheless, AR offers a notable buffer delaying the transition of the land C sink from a sink to a source by ~10 years. This demonstrates that a rapid deployment of AR early on can provide a continuous C sink acting on both short and long timescales. However, atmospheric C removal is less than total land C uptake due to the emerging carbon-cycle feedbacks63. In particular, oceanic C uptake is dependent on atmospheric CO2 partial pressure and hence part of the increase in the land C sink is compensated for by less uptake by the ocean as atmospheric CO2 concentrations decrease, or even transits to being a source64–66, as has also been reported here.

Across the literature a wide range of potentials has been reported (Fig. 5, see “Methods”, and Supplementary Data 1), with our estimates of cumulative C sequestration under AR compared to REF being towards the upper range of AR6-SD (300 scenarios for which information is available) throughout the 21st century. This potential is within the upper range of estimates found in previous literature during early- and mid-century but falls within the mid-range by 2100. When it comes to yearly sequestration rates, our results lie within the upper range of AR6-SD and previous literature estimates during the early- and mid-century, before levelling off and reaching the lower-range of estimates by 2100. Notably, our estimate reaches ~50% of the estimate of Bastin et al.67, who consider an AR area of similar magnitude (~900 Mha) yielding 752 GtCO2, and received wide attention and criticism68–74. Similarly, our estimate reaches ~55% of the most recent comprehensive estimate of Mo et al.75 who employed several ground-sourced and satellite-derived approaches and suggested that 694 GtCO2 could be sequestered through reforestation, in agreement with Bastin et al.67. However, the AR pattern employed here also includes afforestation over land that is inevitably less productive, which could partially explain this difference, while it should be noted that the estimates of Mo et al.75 are based on current climatic conditions. Our estimate remains within the lower boundary of the 2.7–17.9 GtCO2/year range reached by 2030 reported by Griscom et al.19, who however consider a potential reforestation area up to 1779 Mha in their calculations, which could partly explain the divergence. Roe et al.20 have reported a potential from 0.5 GtCO2/year, which corresponds to a low ESM-based estimate, up to 10.1 GtCO2/year, which includes global reforestation constrained by food security and biodiversity based on Griscom et al.19, and on average aligns with the results reported here.

It should be noted that literature estimates are the outcome of different combinations of available land, temporal dynamics of AR deployment, emission trajectories, and assumptions with regards to the capacity of forests and soils to sequester C across different biomes, and hence such a wide range can be expected7,19,20,67,76,77. Most importantly, only 2 of the studies compiled in Fig. 5 employ a fully coupled ESM in an emission-driven setup that accounts for climate and carbon cycle feedbacks39,53, whereas the vast majority includes either offline estimates of sequestration potentials and thus is not easily relatable to changes in atmospheric CO2 concentration and climate19,67, or IAM-based estimates (the AR6-SD scenarios, and e.g., ref. 76,78). IAMs have different configurations and represent the climate and carbon-cycle feedbacks with different levels of complexity79, but they lack the detail of process representation available in ESMs, which is necessary to robustly capture the efficiency of AR, especially when it comes to carbon fluxes under overshoot scenarios80.

Even though pooling multiple AR6-SD scenarios together eliminates the bias of selecting a single IAM scenario, and utilizing regional data allows for the regional-level dynamics to be preserved, it should still be acknowledged that with the approach employed here, the internal consistency found in single IAM scenarios is inevitably lacking. Disaggregating to the gridcell level based on biodiversity and restoration potential maps, not implementing land use transitions other than AR, and independetly selecting an emission trajectory as is done here can also break internal IAM logic. For example, while opting for the SSP5-3.4os emission trajectory allows for a robust assessment of temperature overshoot dynamics, achieving net-negative emissions in SSP5-3.4os is heavily based on carbon capture and storage (CCS)16. This suggests that additional land pressure would be exerted by the need for large-scale bioenergy CCS application, thus increasing competition with AR over land. However, it should be noted that inconsistencies are introduced even when a single IAM scenario is run with an ESM, due to differences in spatial resolution, land use patterns, and discrepancies in the representation of the biosphere and the carbon cycle81,82.

Contextualizing ambitious Afforestation/Reforestation

The present study does not propose a spatiotemporal AR pattern to be followed, but rather aims at overcoming the limitations of highly idealized patterns and of scenarios employing only moderate levels of AR that fall short of the range of total country pledges, by harnessing a multitude of available IAM-generated scenarios coupled with restoration potential maps and biodiversity proxies34,38–41. Using the AR6-SD ensemble allows us to consider constraints of land competition, and technoeconomic, environmental, and societal feasibility, to the extent that those are accounted for in the IAMs, and the scenario here is hence spatiotemporally plausible only to the extent that IAM scenarios can be treated as such. Concerns regarding the uncertainties and weaknesses of IAMs have been reported83, such as their lack of social and institutional considerations, and their inability to conceive more radical societal reorganization and policy challenging strategies, thus not being able to capture the full socioeconomic possibility space84–86. However, it should be noted that the IAMs still remain the main tool at hand when it comes to contextualizing the technical, social, and economic developments in the world, which is necessary for developing a constrained AR pattern that goes beyond fully idealized setups.

At the same time, it should be acknowledged that future AR will inevitably include reversing historical deforestation, not only because such regions can naturally support growth, but also if afforestation is to be avoided87. This increases the relevance of patterns aimed at prioritizing reforestation, as is the case here. Nevertheless, the feasibility of setting such an ambitious global AR target in the range of country pledges can be questioned24,25,52,88. Overall, the pledges by 9 countries amount to 90% of total AR pledges globally (Supplementary Fig. 1), including the United States (25%), Canada (4%) and China (4%), whose net-zero targets are considered of lower confidence rating, and Saudi Arabia (42%), India (4%), Ethiopia (4%) and Australia (4%), whose net-zero targets are of much lower confidence rating based on Rogelj et al.13. The scenario employed here shows less concentration of AR with 39 countries amounting to 90% of the global target (Supplementary Fig. 1). However, this also relies on countries whose net-zero targets are considered of lower (United States, China, Russian Federation, Colombia) and much lower (Brazil, Argentina) confidence rating13. This credibility issue can be even more alarming when one considers that total pledges would likely grow in the future, as more countries set their long-term strategies25.

In our scenario grazing land is significantly reduced and no deforestation occurs, possibly implying a strong dietary shift and reduction of meat consumption89, coupled with a sustainable intensification of remaining grazing lands and regulations to prevent deforestation elsewhere as a compensating mechanism for lost agricultural land90. This fits with a land-sparing paradigm that can indeed offer high mitigation potentials91 and benefit biodiversity92,93; however, the extent to which this can be considered realistic from a socioeconomic perspective is not investigated here. It should also be noted that given fair burden sharing considerations, our scenario, which largely includes forest expansion over the global South, would require financing AR deployment over these regions by the more developed countries94, as is also suggested by the pledge of Saudi Arabia. Until this day, lack of public and private finance has been one of the key barriers to meeting global restoration needs, and scaling-up funding would be a big challenge95.

Socioeconomic factors could also possibly pose barriers to AR implementation or a threat to the permanence of a newly planted forest, while the possible societal consequences and associated tradeoffs should also be considered45,96. In particular, the rates of sequestration achieved here can carry a high risk of exceeding sustainability thresholds97, while the strong dependence of IAMs (and consequently of our scenario) on land-based mitigation in the Global South to reach ambitious climate targets might pose high risks to food security and raise equity concerns88. Comparing the AR pattern employed here with socioeconomic factors such as poverty, population density, governance, land tenure, and indigenous and community land is thus crucial (see “Methods”, and Supplementary Fig. 10). In particular, 6.5% (61 Mha) of the AR in our scenario is applied over gridcells with a significant extent of indigenous and community land (>30%). Planting forests in these regions can often violate the will of indigenous people who can be strongly tied to their land spiritually, financially, and/or are nutritionally dependent on local food production24, and lead to forced physical or economic displacement98. At the same time, 17.6% (165 Mha) of AR is deployed over regions where more than 25% of the population live below the international poverty threshold. Implementing AR over poverty-stricken regions can carry the risk of depriving people of their livelihoods and exacerbating poverty, even though positive outcomes on livelihoods have also been reported99–101. Additionally, 2% (19 Mha) of AR is applied over regions with a population density higher than 500 people/km2 where rates of human disturbance and deforestation can be higher102,103, thus likely jeopardizing permanence.

At the same time, 18% (169 Mha) of AR is applied over regions with poorer governance (governance indicator <0.3), and 16% (151 Mha) over regions with insecure land (>50% of survey respondents perceiving their land or property as being insecure). Poor governance and weak rule of law are considered to be significant barriers to successful implementation and permanence of AR, facilitating forest loss and degradation78,104. Land tenure insecurity can significantly threaten the establishment and permanence of a newly planted forest105, since clear and secure land tenure can determine farmers’ decision to adopt reforestation and maintain a forest over the long-term101,106.

Outlook

Societies have been (re)planting trees for centuries, and AR has reasonably emerged as an efficient way to naturally store C6, while at the same time restoring ecosystems that offer a wide variety of services to our communities and support life107. Our results based on MPI-ESM clearly demonstrate that ambitious AR in the range of country pledges can mitigate a temperature overshoot by lowering global peak temperature, overshoot duration, and end-of-century temperature. Temperature is mostly reduced locally at the sites of forestation, suggesting that the overall biogeochemical cooling effect of the increased land C sink can dominate, or at least compensate for any biogeophysically-induced warming that can arise. Mitigation emerges both for average and for extreme temperature, yet to infer upon the appropriateness of large-scale AR, one should also consider potential AR-induced changes in climate hazard exposure, such as droughts108, rainfall variability and extremes59,109, heat stress110, and compound events111, that could increase adaptation needs locally or regionally112. However, such a multi-dimensional investigation is beyond the scope of the current study and could be the focus of future work, potentially also employing state-of-the-art high-resolution convection-permitting models that can take into account the effects of fine-scale topography, and the local features of atmospheric convection on climatic extremes61,113. Future work could also focus on investigating to what extent the model behavior reported herein can be model- or pattern-specific by comparing multiple ESMs and AR patterns under an emission-driven configuration.

The mitigation potential demonstrated here constitutes ambitious AR as a useful complementary short- and long-term mitigation tool for climate action even under a scenario with strongly reduced emissions, where fertilization of vegetation by CO2 is weaker. However, the scale of mitigation achievable by such high-ambition AR scenarios—with 0.08 oC and 0.2 oC decrease of global mean peak and end-of-century temperature respectively—clearly shows that AR does not alleviate the need for high ambitions in emission reduction52. Importantly, even though a normative judgment on the desirability of ambitious AR is not made here, the results demonstrate the possible socioeconomic tradeoffs associated with it, as well as the significant barriers to implementation and possible threats to permanence that exist.

Methods

Model and experimental setup

We run fully coupled ESM simulations using the CMIP6 configuration of MPI-ESM (MPI-ESM-1-2.01p7-LR)114. MPI-ESM and its land surface component JSBACH have been widely applied, evaluated, and compared against observations and other models, and have been generally found to perform well for various key land surface variables115–117, and for both biogeophysical and biogeochemical effects35,118. The model is run in an emission-driven setup, which means that atmospheric CO2 concentration is not prescribed, but is rather dynamically updated. The fully coupled emission-driven setup offers a dynamic interaction of water, carbon, and energy fluxes between the land, the atmosphere, and the ocean, thus allowing for a full representation of the complex biogeophysical/biogeochemical effects and feedbacks of AR within the entire Earth system.

Our experimental setup includes an AR and a reference (REF) scenario spanning from 2015 to 2100, following SSP5-3.4os, except for land use. SSP5-3.4os reaches a total radiative forcing of 3.4 W/m2, and is an overshoot scenario that follows SSP5-8.5 emissions until 2040, followed by rapid decarbonization, thus reaching net-zero around 2070 and then net-negative (~−3.8 GtCO2/year)16,119. To fully isolate the total effects of AR on the Earth system, which is estimated as the mere difference between the AR and REF scenarios, we keep land use constant at the 2015 state under REF. Hence, we do not include any avoided deforestation, or other land use transitions, which would further complicate the accurate estimation of CDR. Towards that direction, even though increased forest area and wood biomass availability could impact the demand for wood products and hence harvest rates via market mechanisms120, we keep wood harvest amount constant at 2015 levels for both scenarios. At the same time, the amount of vegetated surface is kept constant for both scenarios, and competition among different Plant Functional Types (PFTs) over vegetated land is switched off, thus not allowing forest to expand as a response to a warming climate, in contrast to the CMIP6 model configuration. This allows for prescribing forest area in each timestep, thus being able to fully isolate the signal of AR, which is the main aim of this study. First, accounting naturally expanding vegetated land could raise a concern as to whether CDR can be claimed over land that becomes available due to our sole failure to decarbonize our economy. Second, the permanence of stored carbon over such land is in fact questionable since, assuming further climate mitigation, forest would contract again, suggesting carbon fluxes towards the atmosphere. We run in total 20 model runs; 10 ensemble members for each scenario, to allow for a robust probabilistic treatment of our results.

Probabilistic treatment

To robustly assess the global dynamics of the overshoot and its characteristics, a probabilistic framework is employed. To estimate the time of signal emergence, the overlapping coefficient (OVL)121 is employed, which is a measure of similarity between the probability density functions (pdfs) of two populations f1(x) and f2(x), and is calculated as follows:1 OVL=∫minf1x,f2(x)dx

Over a temporally moving window, global mean 2 m air temperature data from all ensemble members are pooled together, creating two samples, whose OVL is estimated with a gaussian kernel density function. This process is repeated by bootstrapping 1000 times, and a mean OVL is estimated. To estimate an OVL threshold below which dissimilarity between AR and REF distributions can be assumed, an autoregressive model (AR (1)) fitted to REF global temperature is used. At every temporal window, first 1000 pairs of stochastically generated temperature timeseries of trend and sample size equal to REF are created, for each of whom OVL is then estimated. The lower 5th percentile of the 1000 OVL values estimated from the synthetic data (coming from the same known AR(1) process) is determined as the OVL threshold, below which two data samples can be considered to come from dissimilar distributions (at the 5% significance level). The year of signal emergence is then calculated as the year when the mean OVL between REF and AR data starts becoming consistently lower than the significance OVL threshold. This process is repeated with a window length ranging from 5 to 10 years and a mean year of emergence is estimated, to make sure that the choice of window length does not bias the results.

To robustly characterize end-of-century temperature while taking into account interannual variability, the full end-of-century 5-year period (2096–2100) is considered. Similarly, to get robust estimates on peak temperature and avoid biases by likely extremely hot or cold years later or early on, the full 5-year consecutive period yielding the highest average temperature for each ensemble member is considered, and the temperature data are then pooled together. The choice of the period length does not bias the results, as shown in Supplementary Fig. 11.

To robustly capture overshoot duration, the signal needs to be isolated from the noise of internal variability. Hence, a smoothing Savitzky-Golay122 filter is applied. Overshoot duration is then estimated with a varying smoothing window length, ranging from 5 to 10 years, to alleviate any bias induced by the choice of window length. This process is repeated 1000 times, where each time a random combination of 10 ensemble members with replacement are chosen. The mean overshoot duration is then estimated as the mean of the values obtained by bootstrapping. To get a clearer picture of the responses of the full spectrum of temperature variability, different percentiles of global average daily temperature are treated with the same approach.

At the grid level, the statistical significance of changes in 2 m temperature (Fig. 6) or other hydroclimatic variables (Fig. 7, Supplementary Fig. 8) during different time periods is inferred at the 5% level. To do this, for each scenario yearly values from all 10 ensemble members during that period are pooled together, and a two-tailed Student’s t-test adjusted to account for lag-1 temporal autocorrelation123 is applied. To assess field significance (i.e., whether the annual average temperature field under AR is statistically significantly different from REF) the probability of at least one false positive among the multitude of gridcell-level tests is accounted for124. If a single gridcell is found to be statistically significantly different after accounting for this probability, then the global null hypothesis that a statistically significant difference between the AR and REF temperature field does not exist can be rejected125. To increase the robustness of the approach and test its sensitivity, a) temperature data at the gridcell level are pooled over a moving window with length ranging from 1 to 10 years, and b) field significance is declared when the null hypothesis is consistently rejected for a consecutive number of years ranging from 1 to 5, defining as year of emergence the starting year of that period. The results are presented in Supplementary Fig. 7.

AR scenario development

To design an ambitious constrained AR scenario following technoeconomic, environmental, and societal constraints, we leverage (a) the AR6 Scenarios Database (AR6-SD)47, (b) the restoration potential map published by Griscom et al.19 (hereafter GRS), (c) the Atlas of Forest and Landscape Restoration Opportunities (hereafter ATL)126, and (d) the Very Low and Low Human Influence maps published by Riggio et al.127 (hereafter VLHI and LHI, respectively). All maps are regridded to the MPI-ESM resolution with a conservative remapping algorithm.

The AR6-SD includes a multitude of IAM-generated scenarios that consider the pressure exerted on land by rising population and food demand, carbon pricing strategies, energy policies and economic costs47. Here, we use the 1259 scenarios for which explicit information on global forest area is available to capture the feasible range of global AR79. The change in forest area, when positive, can represent not only active AR deployment, but also abandonment of agricultural land, and natural regrowth and succession dynamics, depending on the IAM20. For each year, we estimate the forest area change compared to the previous year for each scenario, and pool all the values together. To obtain an ambitious AR target at the global level which is in the range of estimates of total country pledges24,25 we use the 90th percentile of the pooled global yearly forest area change estimates, which cumulatively reaches 595 Mha by 2060, and 935 Mha by 2100 (Fig. 1). This forest area increase would roughly correspond to its complementary threshold, the 10th percentile of cropland and pasture area change, which suggests that AR is mostly applied at the expense of grazing lands across AR6-SD, while croplands remain relatively stable (Fig. 1). To mimic that behavior, AR only replaces grazing land in our scenario. It is important to note that we do not select individual scenarios, but rather a high percentile across all values each year, and that we refrain from making any judgment upon the likelihood of the individual AR6-SD scenarios, which are treated as different equiprobable narratives across the feasibility space79, nor do we consider the high percentile utilized here as being less probable than the mean49.

Given the lack of pledges for many countries globally, and that the significantly high pledge of Saudi Arabia (42% of global) can only be met through international offsets in addition to domestic land-based CDR25, the spatial disaggregation of global AR pledges is not made based on the individual country pledges themselves, but rather regionally based on AR6-SD estimates. Unavoidably, this creates a spatial mismatch between the scenario employed here and pledges at the country level (Supplementary Fig. 1). However, it should be noted that employing the AR6-SD allows for an explicit treatment of the temporal evolution of AR across the century based on technoeconomic considerations, instead of simply interpolating pledges up to 2060, and arbitrarily extrapolating thereafter.

Here, we use information available in AR6-SD to disaggregate the global yearly AR target across the 5 world economic regions shown in Fig. 3; a regionalization typically used in IAMs. Even though different regionalizations are also available in AR6-SD, the one chosen here offered the most available scenarios (1124). In a similar fashion, for every region we pool available scenarios together and select the 90th percentile of the pooled regional yearly forest area change. Given that a different number of scenarios is available at the regional compared to the global level, and that a percentile for every region is chosen rather than a specific scenario consistent across regions, we rescale the regional estimates so that their sum matches the yearly global target (preserving their relative magnitude with respect to their sum). Interpolation of the 5-year AR6-SD data to yearly both at the global and regional level is performed with a piecewise cubic hermite interpolating polynomial algorithm.

The regional-scale yearly AR target is further distributed spatially, from regional level to the MPI-ESM grid by following the GRS and ATL maps. We consider the ATL restoration potential that includes reforestation and tree integration into mixed-use landscapes and excludes croplands126. GRS is a subset of ATL further constrained by food security and biodiversity concerns and excludes a range of areas from AR: croplands, regions where forest is not the native cover type (afforestation), regions with dense rural population, regions with intensive management, and boreal regions that could potentially yield strong local warming19. Both GRS and ATL potentials per gridcell are expressed as fractions of the gridcell area. Using these maps to guide our scenarios takes implicitly into account the considerations therein. For a given year and region, the AR target is first distributed across gridcells prioritized by GRS. The gridcells considered for AR are the ones where (a) grazing land is available, (b) the remaining GRS potential is greater than zero, (c) forest area has not already increased more than 10% in that year, to avoid excessive AR rates that could be institutionally challenging or incompatible with sustainable local socioeconomic changes, and (d) at least 20% forest cover has historically been or could potentially be sustained, to not use climatically unsuitable regions that would largely require management practices. The target is uniformly applied to the gridcells considered; however, the amount of AR applied over each is limited by (a–d). AR applied within each step reduces the remaining GRS potential by the same amount. This is an iterative process until either the regional target is met, or there are no remaining gridcells to be considered because (a–d) are no longer satisfied anywhere in the given region. In the latter case, the remaining AR target is filled by gridcells prioritized by ATL with a similar iterative process, and if still needed, the rest of the gridcells (satisfying (a–d)) can also be used to meet the AR yearly target. As a result, the AR pattern is constrained first by GRS (362 Mha in the resulting pattern), and then by ATL (375 Mha), but meeting the regional AR target each year can also mean increasing forest cover beyond total restoration potential (GRS and ATL combined) and reaching other gridcells (not prioritized by GRS or ATL), if needed (198 Mha). Reforestation of historically deforested sites is thus strongly prioritized in our approach, but reaching an ambitious target while preserving croplands globally constitutes afforestation (i.e., increasing tree cover beyond 1850 levels) and utilization of land not prioritized by restoration potential maps inevitable.

To account for biodiversity protection, heavily managed pastures are replaced first within every iterative step of the distributive process, and then the more lightly managed rangelands are replaced progressively, starting from the less to the more biodiverse. To protect old-growth grassy biomes which can be of huge ecological importance and exhibit high biodiversity richness87, we characterize rangeland biodiversity by using the VLHI map to determine which rangelands can be considered as being closer to a pristine, and thus as more biodiverse state, and exclude them. Natural grassy biomes which are not characterized as grazing land are also excluded. Then, we use the LHI map to split remaining rangelands into 4 biodiversity groups, based on the level of agreement among the datasets used by Riggio et al.127. By using these maps as a proxy for biodiversity, we assume that the higher the dataset agreement is, the more likely a gridcell experiences low human influence, and hence it is more likely less heavily managed, with less infrastructure, human stress, and cattle densities, weaker land fragmentation, and thus, probably closer to its pristine state, and hence likely more biodiverse127. Even though anthropogenic disturbance is known to be a major driver of biodiversity loss128, biodiversity richness is not necessarily associated with the level of human influence129. At the same time, although abandoning pastures can promote biodiversity, in some cases this does not translate to ecosystems returning to their initial biodiverse state130. Even though we acknowledge these shortcomings, a more detailed treatment of biodiversity concerns is beyond the scope of this study, and the usage of largely available metrics of human influence is preferred here127.

Evidently, in our study no biophysical or climatic limit is applied when considering gridcells eligible for AR, other than whether forest cover has or could potentially (based on MPI-ESM estimates) be at least 20%. Stricter limitations would reduce the amount of grazing land used for AR, since grazing land is often found in semi-arid, transitional regimes (Fig. 2). Hence, increasing forest cover over such areas possibly implies the need of some form of forest management in the real world to support forest growth. Still, carbon sequestration is climate-dependent in MPI-ESM, and more arid regions typically yield lower productivity and sequestration rates (Fig. 4c), and hence afforesting more arid rangelands would not contribute to the mitigation potential in the same way as the more productive areas included in our scenario.

It should also be noted that our simulations start in 2015 and do not include historical land use transitions having occurred later than 2014. According to FAO data131, forest area has decreased by ~30 Mha from 2015 to 2021. By that time in our model AR amounts to ~5 Mha, which corresponds to ~0.5% of total AR application, thus yielding a net difference of ~35 Mha of forest area. By 2021 land C stocks are increased by ~3 GtCO2 (0.81 PgC) in the AR compared to the REF scenario, constituting 0.78% of total additional sequestration (382 GtCO2). Clearly, the bulk of forest area increase and additional C sequestration under AR occur after 2025 in our simulations, which suggests that our results remain relevant with respect to overshoot dynamics and mitigation, regardless of this discrepancy with recent historical data.

An overview of the methodology used to develop our scenario is shown in Fig. 2 and the detailed resulting pattern is presented in Fig. 3 and Supplementary Figs. 1 & 2.

To test the sensitivity of the AR scenario development to the different constraints employed, different configurations are considered. In particular, we develop scenarios where: a) only the restoration map published by Griscom et al.19 (GRS) and the Atlas of Forest and Landscape Restoration Opportunities126 (ATL) map (i.e., no explicit biodiversity consideration), b) only GRS, and biodiversity maps127, and c) no restoration or biodiversity maps are used to constrain the scenario development.

Results are shown in Supplementary Fig. 12. Even though differences arise between the configurations, the patterns share some similarity with the AR scenario employed in this study. It should be noted that pastures are an integral part of both GRS and ATL (the first being a subset of the latter). As a result, prioritizing pastures can partly compensate for not using restoration potential maps to guide AR, due to the inherent consistency between the two. At the same time, consideration of biodiversity does not affect the partitioning between rangelands and pastures that are given up, since biodiversity is only considered for rangelands. Therefore, since pastures are generally given up first in our algorithm regardless of the configuration tested, the rough AR pattern is not heavily sensitive to the different constraints employed. However, considering biodiversity affects the spatial pattern of AR, mostly since it determines the regions where rangelands can be considered closer to a pristine state, and thus excluded from AR. Even though the AR pattern by 2100 unavoidably converges to the available grazing land, the differences across time can be more pronounced as the partitioning between pastures and rangelands, and the specifics of each configuration change.

Literature estimates of sequestration potential

In Fig. 5, the AR sequestration estimates demonstrated here are compared with potentials that have been reported in previous literature. The list of 40 studies presented includes (a) a subset of the comprehensive list of studies provided by “CO2 removal.org”5,7,132 for which explicit quantitative information on AR sequestration potentials at the global scale is provided (19 studies in total)19,39,76,133–148, (b) studies more recent than5,7,132, which are not included in the comprehensive list (16 studies)20,52,67,75,77,78,149–158, and (c) studies published prior to5,7,132 that we identified as missing from the comprehensive list (5 studies)53,159–162. The complete list of studies including the estimates of AR sequestration potentials are listed in Supplementary Data 1.

Across the identified studies we have included only the ones for which cumulative sequestration AR potentials over a defined period, and/or instant sequestration fluxes at a defined year or period can be confidently estimated, based on the information available in each study. Studies for which sequestration potentials cannot be attributed to an increase in forest cover and mostly include other nature-based climate solutions or contain only regional estimates have been omitted. Importantly, modeling studies that are employing idealized AR patterns34,163–165 are not considered.

For studies reporting multiple sequestration potentials, we define the range of estimates as [minimum-maximum] and calculate the mean potential across all values. For studies providing a range of values, we also estimate their mean. It should be noted, however, that not all estimates are entirely independent from one another, but quite a few studies rather build upon prior knowledge provided in the literature. As a result, the pool of studies is not treated as a statistical sample of independent values, and hence descriptive statistics (e.g., mean, standard deviation) are not employed.

The AR6 Scenarios Database (AR6-SD) estimates of AR sequestration potentials is based on a subset of 300 scenarios for which: a) explicit information on AR sequestration is available, b) the three exclusion criteria defined by Prütz et al.166 are met. Therefore, this pool of available scenarios includes only a subset of the 1259 scenarios used to guide the spatiotemporal design of the AR pattern of this study.

Socioeconomic indicators

To aid in the discussion of the possible socioeconomic risks and trade-offs associated with the implementation of the AR scenario employed in this study, several spatially explicit socioeconomic indicators are used. In particular, we use:Governance indicator: The Worldwide Governance Indicators (WGI) provided by the World Bank are used167. The WGI includes country-level indicators on voice and accountability, political stability and absence of violence/terrorism, government effectiveness, regulatory quality, rule of law, and control of corruption over the period 1996–2022. Missing values are filled by taking the arithmetic mean of the specific world subregion that a country belongs to, based on the 22 geographical subregions defined in the United Nations Geoscheme by the United Nations Statistics Division. Following Andrijevic et al.168 a composite governance indicator is estimated by taking the arithmetic mean of the six WGI components for 2022. The composite indicator is found to correlate well (>0.8 correlation) with all of the underlying components.

Land tenure insecurity indicator: The Prindex 2020 global dataset169,170 of land tenure insecurity is used. This indicator describes land tenure and property insecurity as perceived by individuals over 140 countries, expressed as a percentage (%) of the survey’s respondents. Missing values are filled as described above.

Poverty indicator: The poverty headcount ratio (%) at $2.15 (2017 purchasing power parity, lineup est. of 2019) from the Global Subnational Atlas of Poverty (GSAP)171 is used. The data are regridded to the resolution of MPI-ESM with conservative remapping. Missing values are filled as described above.

Indigenous and community land: Indigenous and community land data from Landmark172 are used. The data are regridded to the resolution of MPI-ESM with conservative remapping. In particular, the land that is considered here consists of:i. Indigenous land—Acknowledged by government: Documented and not documented.

ii. Indigenous land— Not acknowledged by government: Held or used with formal land claim submitted and held or used under customary tenure.

iii. Community land—Acknowledged by government: Documented and not documented.

iv. Community land—Not acknowledged by government: Held or used with formal land claim submitted and held or used under customary tenure.

Population density: The WorldPop population count173 is used. The data are regridded from the ~1 km resolution to the resolution of MPI-ESM by summing, and are then divided by gridcell area to obtain population density.

An overview of the socioeconomic indicators overlaid with the AR pattern, and some indicative estimates of how AR is distributed across regions with various levels of governance, land tenure insecurity, poverty, perscentage of indigenous and community land, and population density are presented in Supplementary Fig. 10.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52508-x.

Acknowledgements

Yiannis Moustakis would like to acknowledge the contributions of Raphael Ganzenmüller and Wolfgang Obermeier, whose comments have helped improve this work. This work was funded by the German Federal Ministry of Education and Research (BMBF), projects CDRSynTra (No. 01LS2101A) and STEPSEC (01LS2102A). This work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under project ID bm1241.

Author contributions

Yiannis Moustakis: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Data Curation, Writing—Original Draft, Review & Editing, Visualization Tobias Nützel: Methodology, Software, Validation, Data Curation, Writing—Review & Editing Hao-Wei Wey: Methodology, Writing—Review & Editing Wenkai Bao: Data Curation, Writing—Review & Editing Julia Pongratz: Methodology, Resources, Writing—Review & Editing, Supervision, Project administration, Funding acquisition.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

The AR6 Scenarios Database (AR6-SD)47 is freely available at AR6 Scenario Explorer and Database hosted by IIASA. The Atlas of Forest and Landscape Restoration Opportunities126 is freely available at https://www.wri.org/data/atlas-forest-and-landscape-restoration-opportunities. The (Very) Low Human Influence maps of Riggio et al.127, and the restoration potential map of Griscom et al.19 are freely available with the respective publications. The Worldwide Governance Indicators (WGI) are freely available at https://www.worldbank.org/en/publication/worldwide-governance-indicators167. The Prindex 2020 global dataset is freely available at https://www.prindex.net/data/170. The poverty headcount ratio is freely available at https://pipmaps.worldbank.org/en/data/datatopics/poverty-portal/home171. The Indigenous and community land data are freely available at https://www.landmarkmap.org/172. The WorldPop population count is freely available at https://www.worldpop.org/173. A repository with data supporting this publication is published in Zenodo at: 10.5281/zenodo.12533125.

Code availability

The Max Planck Insitute’s Earth System Model (MPI-ESM-1-2.01p7-LR) is made available under a version of the MPI-M software license agreement (the license and information on how to access the code can be found here: https://code.mpimet.mpg.de/projects/mpi-esm-license). Python 3.11.2 has been used for all data analysis.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Babiker, M. et al. 2022: Cross-sectoral perspectives. in Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds P. R. Shukla, J. Skea, R. Slade, A. Al Khourdajie, R. van Diemen, D. McCollum, M. Pathak, S. Some, P. Vyas, R. Fradera, M. Belkacemi, A. Hasija, G. Lisboa, S. Luz, J. Malley) (Cambridge University Press, Cambridge, UK and New York, 2022).
2. Luderer G Residual fossil CO2 emissions in 1.5–2 °C pathways Nat. Clim. Change 2018 8 626 633 10.1038/s41558-018-0198-6
Luderer, G. et al. Residual fossil CO2 emissions in 1.5–2 °C pathways. Nat. Clim. Change 8, 626–633 (2018).
3. Peters GP The challenge to keep global warming below 2 °C Nat. Clim. Change 2013 3 4 6 10.1038/nclimate1783
Peters, G. P. et al. The challenge to keep global warming below 2 °C. Nat. Clim. Change 3, 4–6 (2013).
4. Rogelj J Scenarios towards limiting global mean temperature increase below 1.5 °C Nat. Clim. Change 2018 8 325 332 10.1038/s41558-018-0091-3
Rogelj, J. et al. Scenarios towards limiting global mean temperature increase below 1.5 °C. Nat. Clim. Change 8, 325–332 (2018).
5. Minx JC Negative emissions—Part 1: research landscape and synthesis Environ. Res. Lett. 2018 13 063001 10.1088/1748-9326/aabf9b
Minx, J. C. et al. Negative emissions—Part 1: research landscape and synthesis. Environ. Res. Lett. 13, 063001 (2018).
6. Smith, S. M. et al. The State of Carbon Dioxide Removal. 1st edn, 10.17605/OSF.IO/W3B4Z (2023).
7. Fuss S Negative emissions—Part 2: costs, potentials and side effects Environ. Res. Lett. 2018 13 063002 10.1088/1748-9326/aabf9f
Fuss, S. et al. Negative emissions—Part 2: costs, potentials and side effects. Environ. Res. Lett. 13, 063002 (2018).
8. Schleussner C-F Science and policy characteristics of the Paris Agreement temperature goal Nat. Clim. Change 2016 6 827 835 10.1038/nclimate3096
Schleussner, C.-F. et al. Science and policy characteristics of the Paris Agreement temperature goal. Nat. Clim. Change 6, 827–835 (2016).
9. Bauer N Exploring risks and benefits of overshooting a 1.5 °C carbon budget over space and time Environ. Res. Lett. 2023 18 054015 10.1088/1748-9326/accd83
Bauer, N. et al. Exploring risks and benefits of overshooting a 1.5 °C carbon budget over space and time. Environ. Res. Lett. 18, 054015 (2023).
10. Riahi K Cost and attainability of meeting stringent climate targets without overshoot Nat. Clim. Chang. 2021 11 1063 1069 10.1038/s41558-021-01215-2
Riahi, K. et al. Cost and attainability of meeting stringent climate targets without overshoot. Nat. Clim. Chang. 11, 1063–1069 (2021).
11. Geden O Löschel A Define limits for temperature overshoot targets Nat. Geosci. 2017 10 881 882 10.1038/s41561-017-0026-z
Geden, O. & Löschel, A. Define limits for temperature overshoot targets. Nat. Geosci. 10, 881–882 (2017).
12. Riahi, K. et al. Mitigation pathways compatible with long-term goals. in Climate change 2022: Mitigation of climate change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, UK and New York, 2022).
13. Rogelj J Credibility gap in net-zero climate targets leaves world at high risk Science 2023 380 1014 1016 10.1126/science.adg6248 37289874
Rogelj, J. et al. Credibility gap in net-zero climate targets leaves world at high risk. Science 380, 1014–1016 (2023).37289874
14. Jewell J Cherp A On the political feasibility of climate change mitigation pathways: is it too late to keep warming below 1.5 °C? WIREs Clim. Change 2020 11 e621 10.1002/wcc.621
Jewell, J. & Cherp, A. On the political feasibility of climate change mitigation pathways: is it too late to keep warming below 1.5 °C? WIREs Clim. Change 11, e621 (2020).
15. Johansson DJA Azar C Lehtveer M Peters GP The role of negative carbon emissions in reaching the Paris climate targets: The impact of target formulation in integrated assessment models Environ. Res. Lett. 2020 15 124024 10.1088/1748-9326/abc3f0
Johansson, D. J. A., Azar, C., Lehtveer, M. & Peters, G. P. The role of negative carbon emissions in reaching the Paris climate targets: The impact of target formulation in integrated assessment models. Environ. Res. Lett. 15, 124024 (2020).
16. Melnikova I Carbon cycle response to temperature overshoot beyond 2 °C: an analysis of CMIP6 models Earth’s Future 2021 9 e2020EF001967 10.1029/2020EF001967
Melnikova, I. et al. Carbon cycle response to temperature overshoot beyond 2 °C: an analysis of CMIP6 models. Earth’s Future 9, e2020EF001967 (2021).
17. Smith P Land-management options for greenhouse gas removal and their impacts on ecosystem services and the sustainable development goals Annu. Rev. Environ. Resour. 2019 44 255 286 10.1146/annurev-environ-101718-033129
Smith, P. et al. Land-management options for greenhouse gas removal and their impacts on ecosystem services and the sustainable development goals. Annu. Rev. Environ. Resour. 44, 255–286 (2019).
18. Pongratz, J. et al. Chapter 7: Current Levels of CDR. in The State of Carbon Dioxide Removal 2024. 2nd edn (eds. Smith, S. M. et al.). https://www.stateofcdr.org/ (2024) 10.17605/OSF.IO/ZXSKB (2024).
19. Griscom BW Natural climate solutions Proc. Natl Acad. Sci. 2017 114 11645 11650 10.1073/pnas.1710465114 29078344
Griscom, B. W. et al. Natural climate solutions. Proc. Natl Acad. Sci. 114, 11645–11650 (2017).29078344
20. Roe S Contribution of the land sector to a 1.5 °C world Nat. Clim. Chang. 2019 9 817 828 10.1038/s41558-019-0591-9
Roe, S. et al. Contribution of the land sector to a 1.5 °C world. Nat. Clim. Chang. 9, 817–828 (2019).
21. Schwartz, N. B., Aide, T. M., Graesser, J., Grau, H. R. & Uriarte, M. Reversals of reforestation across Latin America limit climate mitigation potential of tropical forests. Front. For. Glob. Change 3, 85 (2020).
22. Pongratz J Land use effects on climate: current state, recent progress, and emerging topics Curr. Clim. Change Rep. 2021 7 99 120 10.1007/s40641-021-00178-y
Pongratz, J. et al. Land use effects on climate: current state, recent progress, and emerging topics. Curr. Clim. Change Rep. 7, 99–120 (2021).
23. Boysen LR The limits to global-warming mitigation by terrestrial carbon removal Earth’s Future 2017 5 463 474 10.1002/2016EF000469
Boysen, L. R. et al. The limits to global-warming mitigation by terrestrial carbon removal. Earth’s Future 5, 463–474 (2017).
24. Dooley, K. et al. The Land Gap Report 2022. (2022).
25. Self, A., Burdon, R., Lewis, J., Riggs, P. & Dooley, K. The Land Gap Report: 2023 Update. (2023).
26. Winckler J Reick CH Pongratz J Robust identification of local biogeophysical effects of land-cover change in a global climate model J. Clim. 2017 30 1159 1176 10.1175/JCLI-D-16-0067.1
Winckler, J., Reick, C. H. & Pongratz, J. Robust identification of local biogeophysical effects of land-cover change in a global climate model. J. Clim. 30, 1159–1176 (2017).
27. Bonan GB Doney SC Climate, ecosystems, and planetary futures: the challenge to predict life in Earth system models Science 2018 359 eaam8328 10.1126/science.aam8328 29420265
Bonan, G. B. & Doney, S. C. Climate, ecosystems, and planetary futures: the challenge to predict life in Earth system models. Science 359, eaam8328 (2018).29420265
28. Jones AD Greenhouse gas policy influences climate via direct effects of land-use change J. Clim. 2013 26 3657 3670 10.1175/JCLI-D-12-00377.1
Jones, A. D. et al. Greenhouse gas policy influences climate via direct effects of land-use change. J. Clim. 26, 3657–3670 (2013).
29. Meier R Empirical estimate of forestation-induced precipitation changes in Europe Nat. Geosci. 2021 14 473 478 10.1038/s41561-021-00773-6
Meier, R. et al. Empirical estimate of forestation-induced precipitation changes in Europe. Nat. Geosci. 14, 473–478 (2021).
30. Portmann R Global forestation and deforestation affect remote climate via adjusted atmosphere and ocean circulation Nat. Commun. 2022 13 5569 10.1038/s41467-022-33279-9 36195588
Portmann, R. et al. Global forestation and deforestation affect remote climate via adjusted atmosphere and ocean circulation. Nat. Commun. 13, 5569 (2022).36195588
31. Winckler J Reick CH Bright RM Pongratz J Importance of Surface Roughness for the Local Biogeophysical Effects of Deforestation J. Geophys. Res. Atmos. 2019 124 8605 8618 10.1029/2018JD030127
Winckler, J., Reick, C. H., Bright, R. M. & Pongratz, J. Importance of Surface Roughness for the Local Biogeophysical Effects of Deforestation. J. Geophys. Res. Atmos. 124, 8605–8618 (2019).
32. Pongratz, J., Reick, C. H., Raddatz, T. & Claussen, M. Biogeophysical versus biogeochemical climate response to historical anthropogenic land cover change. Geophys. Res. Lett. 37, L08702 (2010).
33. Boysen LR Global climate response to idealized deforestation in CMIP6 models Biogeosciences 2020 17 5615 5638 10.5194/bg-17-5615-2020
Boysen, L. R. et al. Global climate response to idealized deforestation in CMIP6 models. Biogeosciences 17, 5615–5638 (2020).
34. De Hertog SJ The biogeophysical effects of idealized land cover and land management changes in Earth system models Earth Syst. Dyn. 2023 14 629 667 10.5194/esd-14-629-2023
De Hertog, S. J. et al. The biogeophysical effects of idealized land cover and land management changes in Earth system models. Earth Syst. Dyn. 14, 629–667 (2023).
35. Winckler J Lejeune Q Reick CH Pongratz J Nonlocal effects dominate the global mean surface temperature response to the biogeophysical effects of deforestation Geophys. Res. Lett. 2019 46 745 755 10.1029/2018GL080211
Winckler, J., Lejeune, Q., Reick, C. H. & Pongratz, J. Nonlocal effects dominate the global mean surface temperature response to the biogeophysical effects of deforestation. Geophys. Res. Lett. 46, 745–755 (2019).
36. Su Y Asymmetric influence of forest cover gain and loss on land surface temperature Nat. Clim. Chang. 2023 13 823 831 10.1038/s41558-023-01757-7
Su, Y. et al. Asymmetric influence of forest cover gain and loss on land surface temperature. Nat. Clim. Chang. 13, 823–831 (2023).
37. Kracher D Nitrogen-related constraints of carbon uptake by large-scale forest expansion: simulation study for climate change and management scenarios Earth’s Future 2017 5 1102 1118 10.1002/2017EF000622
Kracher, D. Nitrogen-related constraints of carbon uptake by large-scale forest expansion: simulation study for climate change and management scenarios. Earth’s Future 5, 1102–1118 (2017).
38. Sonntag S Quantifying and comparing effects of climate engineering methods on the Earth system Earth’s Future 2018 6 149 168 10.1002/2017EF000620
Sonntag, S. et al. Quantifying and comparing effects of climate engineering methods on the Earth system. Earth’s Future 6, 149–168 (2018).
39. Sonntag S Pongratz J Reick CH Schmidt H Reforestation in a high-CO2 world—higher mitigation potential than expected, lower adaptation potential than hoped for Geophys. Res. Lett. 2016 43 6546 6553 10.1002/2016GL068824
Sonntag, S., Pongratz, J., Reick, C. H. & Schmidt, H. Reforestation in a high-CO2 world—higher mitigation potential than expected, lower adaptation potential than hoped for. Geophys. Res. Lett. 43, 6546–6553 (2016).
40. Matthews HD Temporary nature-based carbon removal can lower peak warming in a well-below 2 °C scenario Commun. Earth Environ. 2022 3 1 8 10.1038/s43247-022-00391-z
Matthews, H. D. et al. Temporary nature-based carbon removal can lower peak warming in a well-below 2 °C scenario. Commun. Earth Environ. 3, 1–8 (2022).
41. Lawrence DM The Land Use Model Intercomparison Project (LUMIP) contribution to CMIP6: rationale and experimental design Geosci. Model Dev. 2016 9 2973 2998 10.5194/gmd-9-2973-2016
Lawrence, D. M. et al. The Land Use Model Intercomparison Project (LUMIP) contribution to CMIP6: rationale and experimental design. Geosci. Model Dev. 9, 2973–2998 (2016).
42. van Vuuren DP Energy, land-use and greenhouse gas emissions trajectories under a green growth paradigm Glob. Environ. Change 2017 42 237 250 10.1016/j.gloenvcha.2016.05.008
van Vuuren, D. P. et al. Energy, land-use and greenhouse gas emissions trajectories under a green growth paradigm. Glob. Environ. Change 42, 237–250 (2017).
43. Loughran TF Limited mitigation potential of forestation under a high emissions scenario: results from multi-model and single model ensembles J. Geophys. Res. 2023 128 e2023JG007605 10.1029/2023JG007605
Loughran, T. F. et al. Limited mitigation potential of forestation under a high emissions scenario: results from multi-model and single model ensembles. J. Geophys. Res. 128, e2023JG007605 (2023).
44. Melnikova I Ciais P Tanaka K Vuichard N Boucher O Relative benefits of allocating land to bioenergy crops and forests vary by region Commun. Earth Environ. 2023 4 1 12 10.1038/s43247-023-00866-7 37325084
Melnikova, I., Ciais, P., Tanaka, K., Vuichard, N. & Boucher, O. Relative benefits of allocating land to bioenergy crops and forests vary by region. Commun. Earth Environ. 4, 1–12 (2023).37325084
45. Holl KD Brancalion PHS Tree planting is not a simple solution Science 2020 368 580 581 10.1126/science.aba8232 32381704
Holl, K. D. & Brancalion, P. H. S. Tree planting is not a simple solution. Science 368, 580–581 (2020).32381704
46. Parry IM Ritchie PDL Cox PM Evidence of localised Amazon rainforest dieback in CMIP6 models Earth Syst. Dyn. 2022 13 1667 1675 10.5194/esd-13-1667-2022
Parry, I. M., Ritchie, P. D. L. & Cox, P. M. Evidence of localised Amazon rainforest dieback in CMIP6 models. Earth Syst. Dyn. 13, 1667–1675 (2022).
47. Byers, E. et al. AR6 Scenarios Database. Intergovernmental Panel on Climate Change10.5281/zenodo.7197970 (2022).
48. Hurtt GC Harmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6 Geosci. Model Dev. 2020 13 5425 5464 10.5194/gmd-13-5425-2020
Hurtt, G. C. et al. Harmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6. Geosci. Model Dev. 13, 5425–5464 (2020).
49. Huppmann D Rogelj J Kriegler E Krey V Riahi K A new scenario resource for integrated 1.5 °C research Nat. Clim. Change 2018 8 1027 1030 10.1038/s41558-018-0317-4
Huppmann, D., Rogelj, J., Kriegler, E., Krey, V. & Riahi, K. A new scenario resource for integrated 1.5 °C research. Nat. Clim. Change 8, 1027–1030 (2018).
50. Sognnaes I A multi-model analysis of long-term emissions and warming implications of current mitigation efforts Nat. Clim. Chang. 2021 11 1055 1062 10.1038/s41558-021-01206-3
Sognnaes, I. et al. A multi-model analysis of long-term emissions and warming implications of current mitigation efforts. Nat. Clim. Chang. 11, 1055–1062 (2021).
51. Warszawski L All options, not silver bullets, needed to limit global warming to 1.5 °C: a scenario appraisal Environ. Res. Lett. 2021 16 064037 10.1088/1748-9326/abfeec
Warszawski, L. et al. All options, not silver bullets, needed to limit global warming to 1.5 °C: a scenario appraisal. Environ. Res. Lett. 16, 064037 (2021).
52. Dooley K Nicholls Z Meinshausen M Carbon removals from nature restoration are no substitute for steep emission reductions One Earth 2022 5 812 824 10.1016/j.oneear.2022.06.002
Dooley, K., Nicholls, Z. & Meinshausen, M. Carbon removals from nature restoration are no substitute for steep emission reductions. One Earth 5, 812–824 (2022).
53. Arora VK Montenegro A Small temperature benefits provided by realistic afforestation efforts Nat. Geosci. 2011 4 514 518 10.1038/ngeo1182
Arora, V. K. & Montenegro, A. Small temperature benefits provided by realistic afforestation efforts. Nat. Geosci. 4, 514–518 (2011).
54. Asaadi A Carbon cycle feedbacks in an idealized simulation and a scenario simulation of negative emissions in CMIP6 Earth system models Biogeosciences 2024 21 411 435 10.5194/bg-21-411-2024
Asaadi, A. et al. Carbon cycle feedbacks in an idealized simulation and a scenario simulation of negative emissions in CMIP6 Earth system models. Biogeosciences 21, 411–435 (2024).
55. Jones CD The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) contribution to C4MIP: quantifying committed climate changes following zero carbon emissions Geosci. Model Dev. 2019 12 4375 4385 10.5194/gmd-12-4375-2019
Jones, C. D. et al. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) contribution to C4MIP: quantifying committed climate changes following zero carbon emissions. Geosci. Model Dev. 12, 4375–4385 (2019).
56. MacDougall AH Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO2 Biogeosciences 2020 17 2987 3016 10.5194/bg-17-2987-2020
MacDougall, A. H. et al. Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO2. Biogeosciences 17, 2987–3016 (2020).
57. Arora VK Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models Biogeosciences 2020 17 4173 4222 10.5194/bg-17-4173-2020
Arora, V. K. et al. Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models. Biogeosciences 17, 4173–4222 (2020).
58. Meehl GA Context for interpreting equilibrium climate sensitivity and transient climate response from the CMIP6 Earth system models Sci. Adv. 2020 6 eaba1981 10.1126/sciadv.aba1981 32637602
Meehl, G. A. et al. Context for interpreting equilibrium climate sensitivity and transient climate response from the CMIP6 Earth system models. Sci. Adv. 6, eaba1981 (2020).32637602
59. Moustakis Y Fatichi S Onof C Paschalis A Insensitivity of ecosystem productivity to predicted changes in fine-scale rainfall variability J. Geophys. Res. 2022 127 e2021JG006735 10.1029/2021JG006735
Moustakis, Y., Fatichi, S., Onof, C. & Paschalis, A. Insensitivity of ecosystem productivity to predicted changes in fine-scale rainfall variability. J. Geophys. Res. 127, e2021JG006735 (2022).
60. De Hertog SJ Effects of idealized land cover and land management changes on the atmospheric water cycle Earth Syst. Dyn. 2024 15 265 291 10.5194/esd-15-265-2024
De Hertog, S. J. et al. Effects of idealized land cover and land management changes on the atmospheric water cycle. Earth Syst. Dyn. 15, 265–291 (2024).
61. Buechel, M. et al. Hydrometeorological response to afforestation in the UK: findings from a kilometer-scale climate model. Environ. Res. Lett. 19, 6 (2024).
62. Luyssaert S Old-growth forests as global carbon sinks Nature 2008 455 213 215 10.1038/nature07276 18784722
Luyssaert, S. et al. Old-growth forests as global carbon sinks. Nature 455, 213–215 (2008).18784722
63. Jones CD Simulating the Earth system response to negative emissions Environ. Res. Lett. 2016 11 095012 10.1088/1748-9326/11/9/095012
Jones, C. D. et al. Simulating the Earth system response to negative emissions. Environ. Res. Lett. 11, 095012 (2016).
64. Chimuka VR Nzotungicimpaye C-M Zickfeld K Quantifying land carbon cycle feedbacks under negative CO2 emissions Biogeosciences 2023 20 2283 2299 10.5194/bg-20-2283-2023
Chimuka, V. R., Nzotungicimpaye, C.-M. & Zickfeld, K. Quantifying land carbon cycle feedbacks under negative CO2 emissions. Biogeosciences 20, 2283–2299 (2023).
65. Vichi M Navarra A Fogli PG Adjustment of the natural ocean carbon cycle to negative emission rates Clim. Change 2013 118 105 118 10.1007/s10584-012-0677-0
Vichi, M., Navarra, A. & Fogli, P. G. Adjustment of the natural ocean carbon cycle to negative emission rates. Clim. Change 118, 105–118 (2013).
66. Zickfeld K Azevedo D Mathesius S Matthews HD Asymmetry in the climate–carbon cycle response to positive and negative CO2 emissions Nat. Clim. Chang. 2021 11 613 617 10.1038/s41558-021-01061-2
Zickfeld, K., Azevedo, D., Mathesius, S. & Matthews, H. D. Asymmetry in the climate–carbon cycle response to positive and negative CO2 emissions. Nat. Clim. Chang. 11, 613–617 (2021).
67. Bastin J-F The global tree restoration potential Science 2019 365 76 79 10.1126/science.aax0848 31273120
Bastin, J.-F. et al. The global tree restoration potential. Science 365, 76–79 (2019).31273120
68. Bastin J-F Forest restoration: transformative trees—response Science 2019 366 317 317 10.1126/science.aaz2148 31624205
Bastin, J.-F. et al. Forest restoration: transformative trees—response. Science 366, 317–317 (2019).31624205
69. Delzeit R Forest restoration: expanding agriculture Science 2019 366 316 317 10.1126/science.aaz0705 31624203
Delzeit, R. et al. Forest restoration: expanding agriculture. Science 366, 316–317 (2019).31624203
70. Lewis SL Mitchard ETA Prentice C Maslin M Poulter B Comment on “The global tree restoration potential” Science 2019 366 eaaz0388 10.1126/science.aaz0388 31624179
Lewis, S. L., Mitchard, E. T. A., Prentice, C., Maslin, M. & Poulter, B. Comment on “The global tree restoration potential”. Science 366, eaaz0388 (2019).31624179
71. Luedeling E Forest restoration: overlooked constraints Science 2019 366 315 315 10.1126/science.aay7988 31624202
Luedeling, E. et al. Forest restoration: overlooked constraints. Science 366, 315–315 (2019).31624202
72. Sheil D Forest restoration: transformative trees Science 2019 366 316 317 10.1126/science.aay7309 31624204
Sheil, D. et al. Forest restoration: transformative trees. Science 366, 316–317 (2019).31624204
73. Skidmore AK Wang T de Bie K Pilesjö P Comment on “The global tree restoration potential” Science 2019 366 eaaz0111 10.1126/science.aaz0111 31780528
Skidmore, A. K., Wang, T., de Bie, K. & Pilesjö, P. Comment on “The global tree restoration potential”. Science 366, eaaz0111 (2019).31780528
74. Veldman JW Comment on “The global tree restoration potential” Science 2019 366 eaay7976 10.1126/science.aay7976 31624182
Veldman, J. W. et al. Comment on “The global tree restoration potential”. Science 366, eaay7976 (2019).31624182
75. Mo L Integrated global assessment of the natural forest carbon potential Nature 2023 624 92 101 10.1038/s41586-023-06723-z 37957399
Mo, L. et al. Integrated global assessment of the natural forest carbon potential. Nature 624, 92–101 (2023).37957399
76. Humpenöder F Investigating afforestation and bioenergy CCS as climate change mitigation strategies Environ. Res. Lett. 2014 9 064029 10.1088/1748-9326/9/6/064029
Humpenöder, F. et al. Investigating afforestation and bioenergy CCS as climate change mitigation strategies. Environ. Res. Lett. 9, 064029 (2014).
77. Roe S Land-based measures to mitigate climate change: potential and feasibility by country Glob. Change Biol. 2021 27 6025 6058 10.1111/gcb.15873
Roe, S. et al. Land-based measures to mitigate climate change: potential and feasibility by country. Glob. Change Biol. 27, 6025–6058 (2021).
78. Doelman JC Afforestation for climate change mitigation: potentials, risks and trade-offs Glob. Change Biol. 2020 26 1576 1591 10.1111/gcb.14887
Doelman, J. C. et al. Afforestation for climate change mitigation: potentials, risks and trade-offs. Glob. Change Biol. 26, 1576–1591 (2020).
79. Guivarch, C. et al. IPCC, 2022: Annex III: Scenarios and modelling methods. in Proc. IPCC, 2022: Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, Cambridge, UK, 2022).
80. Melnikova I Ciais P Boucher O Tanaka K Assessing carbon cycle projections from complex and simple models under SSP scenarios Clim. Change 2023 176 168 10.1007/s10584-023-03639-5
Melnikova, I., Ciais, P., Boucher, O. & Tanaka, K. Assessing carbon cycle projections from complex and simple models under SSP scenarios. Clim. Change 176, 168 (2023).
81. Cheng Y A bioenergy-focused versus a reforestation-focused mitigation pathway yields disparate carbon storage and climate responses Proc. Natl Acad. Sci. 2024 121 e2306775121 10.1073/pnas.2306775121 38315850
Cheng, Y. et al. A bioenergy-focused versus a reforestation-focused mitigation pathway yields disparate carbon storage and climate responses. Proc. Natl Acad. Sci. 121, e2306775121 (2024).38315850
82. Melnikova I Impact of bioenergy crop expansion on climate–carbon cycle feedbacks in overshoot scenarios Earth Syst. Dyn. 2022 13 779 794 10.5194/esd-13-779-2022
Melnikova, I. et al. Impact of bioenergy crop expansion on climate–carbon cycle feedbacks in overshoot scenarios. Earth Syst. Dyn. 13, 779–794 (2022).
83. Weyant J Some contributions of integrated assessment models of global climate change Rev. Environ. Econ. Policy 2017 11 115 137 10.1093/reep/rew018
Weyant, J. Some contributions of integrated assessment models of global climate change. Rev. Environ. Econ. Policy 11, 115–137 (2017).
84. Proctor, J. C. Expanding the possible: exploring the role for heterodox economics in integrated climate-economy modeling. Rev. Evol. Polit. Econ. 4, 537–557 (2023).
85. van Beek L Oomen J Hajer M Pelzer P van Vuuren D Navigating the political: an analysis of political calibration of integrated assessment modelling in light of the 1.5 °C goal Environ. Sci. Policy 2022 133 193 202 10.1016/j.envsci.2022.03.024
van Beek, L., Oomen, J., Hajer, M., Pelzer, P. & van Vuuren, D. Navigating the political: an analysis of political calibration of integrated assessment modelling in light of the 1.5 °C goal. Environ. Sci. Policy 133, 193–202 (2022).
86. van Beek L Hajer M Pelzer P van Vuuren D Cassen C Anticipating futures through models: the rise of Integrated Assessment Modelling in the climate science-policy interface since 1970 Glob. Environ. Change 2020 65 102191 10.1016/j.gloenvcha.2020.102191
van Beek, L., Hajer, M., Pelzer, P., van Vuuren, D. & Cassen, C. Anticipating futures through models: the rise of Integrated Assessment Modelling in the climate science-policy interface since 1970. Glob. Environ. Change 65, 102191 (2020).
87. Bond WJ Ancient grasslands at risk Science 2016 351 120 122 10.1126/science.aad5132 26744392
Bond, W. J. Ancient grasslands at risk. Science 351, 120–122 (2016).26744392
88. Jaiswal S Nagarajan A Mythri A Projecting a food insecure world: equity implications of land-based mitigation in IPCC mitigation pathways Environ. Sci. Policy 2024 155 103724 10.1016/j.envsci.2024.103724
Jaiswal, S., Nagarajan, A. & Mythri, A. Projecting a food insecure world: equity implications of land-based mitigation in IPCC mitigation pathways. Environ. Sci. Policy 155, 103724 (2024).
89. Eshel G Stainier P Shepon A Swaminathan A Environmentally optimal, nutritionally sound, protein and energy conserving plant based alternatives to U.S Meat. Sci. Rep. 2019 9 10345 10.1038/s41598-019-46590-1 31395893
Eshel, G., Stainier, P., Shepon, A. & Swaminathan, A. Environmentally optimal, nutritionally sound, protein and energy conserving plant based alternatives to U.S. Meat. Sci. Rep. 9, 10345 (2019).31395893
90. Fuchs R Brown C Rounsevell M Europe’s Green Deal offshores environmental damage to other nations Nature 2020 586 671 673 10.1038/d41586-020-02991-1 33106645
Fuchs, R., Brown, C. & Rounsevell, M. Europe’s Green Deal offshores environmental damage to other nations. Nature 586, 671–673 (2020).33106645
91. Lamb A The potential for land sparing to offset greenhouse gas emissions from agriculture Nat. Clim. Change 2016 6 488 492 10.1038/nclimate2910
Lamb, A. et al. The potential for land sparing to offset greenhouse gas emissions from agriculture. Nat. Clim. Change 6, 488–492 (2016).
92. Ekroos, J. et al. Sparing land for biodiversity at multiple spatial scales. Front. Ecol. Evolut. 3, 145 (2016).
93. Kremen C Reframing the land-sparing/land-sharing debate for biodiversity conservation Ann. N. Y. Acad. Sci. 2015 1355 52 76 10.1111/nyas.12845 26213864
Kremen, C. Reframing the land-sparing/land-sharing debate for biodiversity conservation. Ann. N. Y. Acad. Sci. 1355, 52–76 (2015).26213864
94. Fyson CL Baur S Gidden M Schleussner C-F Fair-share carbon dioxide removal increases major emitter responsibility Nat. Clim. Chang. 2020 10 836 841 10.1038/s41558-020-0857-2
Fyson, C. L., Baur, S., Gidden, M. & Schleussner, C.-F. Fair-share carbon dioxide removal increases major emitter responsibility. Nat. Clim. Chang. 10, 836–841 (2020).
95. Löfqvist S Garrett RD Ghazoul J Incentives and barriers to private finance for forest and landscape restoration Nat. Ecol. Evol. 2023 7 707 715 10.1038/s41559-023-02037-5 37165107
Löfqvist, S., Garrett, R. D. & Ghazoul, J. Incentives and barriers to private finance for forest and landscape restoration. Nat. Ecol. Evol. 7, 707–715 (2023).37165107
96. Schultz B Recognizing the equity implications of restoration priority maps Environ. Res. Lett. 2022 17 114019 10.1088/1748-9326/ac9918
Schultz, B. et al. Recognizing the equity implications of restoration priority maps. Environ. Res. Lett. 17, 114019 (2022).
97. Deprez A Sustainability limits needed for CO2 removal Science 2024 383 484 486 10.1126/science.adj6171 38301011
Deprez, A. et al. Sustainability limits needed for CO2 removal. Science 383, 484–486 (2024).38301011
98. Erbaugh JT Global forest restoration and the importance of prioritizing local communities Nat. Ecol. Evol. 2020 4 1472 1476 10.1038/s41559-020-01282-2 32839542
Erbaugh, J. T. et al. Global forest restoration and the importance of prioritizing local communities. Nat. Ecol. Evol. 4, 1472–1476 (2020).32839542
99. Adams C Rodrigues ST Calmon M Kumar C Impacts of large-scale forest restoration on socioeconomic status and local livelihoods: what we know and do not know Biotropica 2016 48 731 744 10.1111/btp.12385
Adams, C., Rodrigues, S. T., Calmon, M. & Kumar, C. Impacts of large-scale forest restoration on socioeconomic status and local livelihoods: what we know and do not know. Biotropica 48, 731–744 (2016).
100. Le HD Smith C Herbohn J Harrison S More than just trees: assessing reforestation success in tropical developing countries J. Rural Stud. 2012 28 5 19 10.1016/j.jrurstud.2011.07.006
Le, H. D., Smith, C., Herbohn, J. & Harrison, S. More than just trees: assessing reforestation success in tropical developing countries. J. Rural Stud. 28, 5–19 (2012).
101. Ota L Herbohn J Gregorio N Harrison S Reforestation and smallholder livelihoods in the humid tropics Land Use Policy 2020 92 104455 10.1016/j.landusepol.2019.104455
Ota, L., Herbohn, J., Gregorio, N. & Harrison, S. Reforestation and smallholder livelihoods in the humid tropics. Land Use Policy 92, 104455 (2020).
102. Borda-Niño M Meli P Brancalion PHS Drivers of tropical forest cover increase: a systematic review Land Degrad. Dev. 2020 31 1366 1379 10.1002/ldr.3534
Borda-Niño, M., Meli, P. & Brancalion, P. H. S. Drivers of tropical forest cover increase: a systematic review. Land Degrad. Dev. 31, 1366–1379 (2020).
103. Busch J Ferretti-Gallon K What drives and stops deforestation, reforestation, and forest degradation? An updated meta-analysis Rev. Environ. Econ. Policy 2023 17 217 250 10.1086/725051
Busch, J. & Ferretti-Gallon, K. What drives and stops deforestation, reforestation, and forest degradation? An updated meta-analysis. Rev. Environ. Econ. Policy 17, 217–250 (2023).
104. Tacconi L Williams DA Corruption and anti-corruption in environmental and resource management Annu. Rev. Environ. Resour. 2020 45 305 329 10.1146/annurev-environ-012320-083949
Tacconi, L. & Williams, D. A. Corruption and anti-corruption in environmental and resource management. Annu. Rev. Environ. Resour. 45, 305–329 (2020).
105. Rakotonarivo OS Resolving land tenure security is essential to deliver forest restoration Commun. Earth Environ. 2023 4 1 8 10.1038/s43247-023-00847-w 37325084
Rakotonarivo, O. S. et al. Resolving land tenure security is essential to deliver forest restoration. Commun. Earth Environ. 4, 1–8 (2023).37325084
106. Arvola A What drives smallholder tree growing? Enabling conditions in a changing policy environment For. Policy Econ. 2020 116 102173 10.1016/j.forpol.2020.102173
Arvola, A. et al. What drives smallholder tree growing? Enabling conditions in a changing policy environment. For. Policy Econ. 116, 102173 (2020).
107. Brockerhoff EG Forest biodiversity, ecosystem functioning and the provision of ecosystem services Biodivers. Conserv 2017 26 3005 3035 10.1007/s10531-017-1453-2
Brockerhoff, E. G. et al. Forest biodiversity, ecosystem functioning and the provision of ecosystem services. Biodivers. Conserv 26, 3005–3035 (2017).
108. Cook BI Mankin JS Anchukaitis KJ Climate change and drought: from past to future Curr. Clim. Change Rep. 2018 4 164 179 10.1007/s40641-018-0093-2
Cook, B. I., Mankin, J. S. & Anchukaitis, K. J. Climate change and drought: from past to future. Curr. Clim. Change Rep. 4, 164–179 (2018).
109. Moustakis Y Papalexiou SM Onof CJ Paschalis A Seasonality, intensity, and duration of rainfall extremes change in a warmer climate Earth’s Future 2021 9 e2020EF001824 10.1029/2020EF001824
Moustakis, Y., Papalexiou, S. M., Onof, C. J. & Paschalis, A. Seasonality, intensity, and duration of rainfall extremes change in a warmer climate. Earth’s Future 9, e2020EF001824 (2021).
110. Vargas Zeppetello LR Raftery AE Battisti DS Probabilistic projections of increased heat stress driven by climate change Commun. Earth Environ. 2022 3 1 7 10.1038/s43247-022-00524-4
Vargas Zeppetello, L. R., Raftery, A. E. & Battisti, D. S. Probabilistic projections of increased heat stress driven by climate change. Commun. Earth Environ. 3, 1–7 (2022).
111. Zscheischler J A typology of compound weather and climate events Nat. Rev. Earth Environ. 2020 1 333 347 10.1038/s43017-020-0060-z
Zscheischler, J. et al. A typology of compound weather and climate events. Nat. Rev. Earth Environ. 1, 333–347 (2020).
112. Berrang-Ford L A systematic global stocktake of evidence on human adaptation to climate change Nat. Clim. Chang. 2021 11 989 1000 10.1038/s41558-021-01170-y
Berrang-Ford, L. et al. A systematic global stocktake of evidence on human adaptation to climate change. Nat. Clim. Chang. 11, 989–1000 (2021).
113. Moustakis Y Onof CJ Paschalis A Atmospheric convection, dynamics and topography shape the scaling pattern of hourly rainfall extremes with temperature globally Commun. Earth Environ. 2020 1 1 9 10.1038/s43247-020-0003-0
Moustakis, Y., Onof, C. J. & Paschalis, A. Atmospheric convection, dynamics and topography shape the scaling pattern of hourly rainfall extremes with temperature globally. Commun. Earth Environ. 1, 1–9 (2020).
114. Mauritsen T Developments in the MPI-M Earth System Model version 1.2 (MPI-ESM1.2) and its response to increasing CO2 J. Adv. Modeling Earth Syst. 2019 11 998 1038 10.1029/2018MS001400
Mauritsen, T. et al. Developments in the MPI-M Earth System Model version 1.2 (MPI-ESM1.2) and its response to increasing CO2. J. Adv. Modeling Earth Syst. 11, 998–1038 (2019).
115. Brovkin V Evaluation of vegetation cover and land-surface albedo in MPI-ESM CMIP5 simulations J. Adv. Modeling Earth Syst. 2013 5 48 57 10.1029/2012MS000169
Brovkin, V. et al. Evaluation of vegetation cover and land-surface albedo in MPI-ESM CMIP5 simulations. J. Adv. Modeling Earth Syst. 5, 48–57 (2013).
116. Eyring V Taking climate model evaluation to the next level Nat. Clim. Change 2019 9 102 110 10.1038/s41558-018-0355-y
Eyring, V. et al. Taking climate model evaluation to the next level. Nat. Clim. Change 9, 102–110 (2019).
117. Hagemann S Loew A Andersson A Combined evaluation of MPI-ESM land surface water and energy fluxes J. Adv. Modeling Earth Syst. 2013 5 259 286 10.1029/2012MS000173
Hagemann, S., Loew, A. & Andersson, A. Combined evaluation of MPI-ESM land surface water and energy fluxes. J. Adv. Modeling Earth Syst. 5, 259–286 (2013).
118. Duveiller G Biophysics and vegetation cover change: a process-based evaluation framework for confronting land surface models with satellite observations Earth Syst. Sci. Data 2018 10 1265 1279 10.5194/essd-10-1265-2018
Duveiller, G. et al. Biophysics and vegetation cover change: a process-based evaluation framework for confronting land surface models with satellite observations. Earth Syst. Sci. Data 10, 1265–1279 (2018).
119. O’Neill BC The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6 Geosci. Model Dev. 2016 9 3461 3482 10.5194/gmd-9-3461-2016
O’Neill, B. C. et al. The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6. Geosci. Model Dev. 9, 3461–3482 (2016).
120. Favero A Daigneault A Sohngen B Forests: carbon sequestration, biomass energy, or both? Sci. Adv. 2020 6 eaay6792 10.1126/sciadv.aay6792 32232153
Favero, A., Daigneault, A. & Sohngen, B. Forests: carbon sequestration, biomass energy, or both? Sci. Adv. 6, eaay6792 (2020).32232153
121. Inman HF Bradley EL The overlapping coefficient as a measure of agreement between probability distributions and point estimation of the overlap of two normal densities Commun. Stat. Theory Methods 1989 18 3851 3874 10.1080/03610928908830127
Inman, H. F. & Bradley, E. L. The overlapping coefficient as a measure of agreement between probability distributions and point estimation of the overlap of two normal densities. Commun. Stat. Theory Methods 18, 3851–3874 (1989).
122. Savitzky, A. & Golay, M. J. Smoothing and differentiation of data by simplified least squares procedures. Anal. Chem. 36, 1627–1639 (1964).
123. Zwiers FW von Storch H Taking serial correlation into account in tests of the mean J. Clim. 1995 8 336 351 10.1175/1520-0442(1995)008<0336:TSCIAI>2.0.CO;2
Zwiers, F. W. & von Storch, H. Taking serial correlation into account in tests of the mean. J. Clim. 8, 336–351 (1995).
124. Bonferroni C Teoria statistica delle classi e calcolo delle probabilita Pubblicazioni del. R. Istituto Super. di Sci. Economiche e Commericiali di Firenze 1936 8 3 62
Bonferroni, C. Teoria statistica delle classi e calcolo delle probabilita. Pubblicazioni del. R. Istituto Super. di Sci. Economiche e Commericiali di Firenze 8, 3–62 (1936).
125. Cortés J Mahecha M Reichstein M Brenning A Accounting for multiple testing in the analysis of spatio-temporal environmental data Environ. Ecol. Stat. 2020 27 293 318 10.1007/s10651-020-00446-4
Cortés, J., Mahecha, M., Reichstein, M. & Brenning, A. Accounting for multiple testing in the analysis of spatio-temporal environmental data. Environ. Ecol. Stat. 27, 293–318 (2020).
126. Minnemeyer, S., Laestadius, L., Sizer, N., Saint-Laurent, C. & Potapov, P. Atlas of Forest and Landscape Restoration Opportunities. https://www.wri.org/data/atlas-forest-and-landscape-restoration-opportunities#project-tabs (2014).
127. Riggio J Global human influence maps reveal clear opportunities in conserving Earth’s remaining intact terrestrial ecosystems Glob. Change Biol. 2020 26 4344 4356 10.1111/gcb.15109
Riggio, J. et al. Global human influence maps reveal clear opportunities in conserving Earth’s remaining intact terrestrial ecosystems. Glob. Change Biol. 26, 4344–4356 (2020).
128. Swift TL Hannon SJ Critical thresholds associated with habitat loss: a review of the concepts, evidence, and applications Biol. Rev. 2010 85 35 53 10.1111/j.1469-185X.2009.00093.x 19930172
Swift, T. L. & Hannon, S. J. Critical thresholds associated with habitat loss: a review of the concepts, evidence, and applications. Biol. Rev. 85, 35–53 (2010).19930172
129. Pimm SL Jenkins CN Li BV How to protect half of Earth to ensure it protects sufficient biodiversity Sci. Adv. 2018 4 eaat2616 10.1126/sciadv.aat2616 30167461
Pimm, S. L., Jenkins, C. N. & Li, B. V. How to protect half of Earth to ensure it protects sufficient biodiversity. Sci. Adv. 4, eaat2616 (2018).30167461
130. Filazzola A The effects of livestock grazing on biodiversity are multi-trophic: a meta-analysis Ecol. Lett. 2020 23 1298 1309 10.1111/ele.13527 32369874
Filazzola, A. et al. The effects of livestock grazing on biodiversity are multi-trophic: a meta-analysis. Ecol. Lett. 23, 1298–1309 (2020).32369874
131. FAO. F. A. O. FAOSTAT Land Use. License: CC BY-NC-SA 3.0 IGO. Extracted from: https://www.fao.org/faostat/en/#data/RL/visualize (2024).
132. Nemet GF Negative emissions—Part 3: Innovation and upscaling Environ. Res. Lett. 2018 13 063003 10.1088/1748-9326/aabff4
Nemet, G. F. et al. Negative emissions—Part 3: Innovation and upscaling. Environ. Res. Lett. 13, 063003 (2018).
133. Brown, S. Management of Forests for Mitigation of Greenhouse Gas Emissions. (1996).
134. Canadell JG Raupach MR Managing Forests for Climate Change Mitigation Science 2008 320 1456 1457 10.1126/science.1155458 18556550
Canadell, J. G. & Raupach, M. R. Managing Forests for Climate Change Mitigation. Science 320, 1456–1457 (2008).18556550
135. Cannell MGR Carbon sequestration and biomass energy offset: theoretical, potential and achievable capacities globally, in Europe and the UK Biomass. Bioenergy 2003 24 97 116 10.1016/S0961-9534(02)00103-4
Cannell, M. G. R. Carbon sequestration and biomass energy offset: theoretical, potential and achievable capacities globally, in Europe and the UK. Biomass. Bioenergy 24, 97–116 (2003).
136. Dixon RK Winjum JK Andrasko KJ Lee JJ Schroeder PE Integrated land-use systems: Assessment of promising agroforest and alternative land-use practices to enhance carbon conservation and sequestration Clim. Change 1994 27 71 92 10.1007/BF01098474
Dixon, R. K., Winjum, J. K., Andrasko, K. J., Lee, J. J. & Schroeder, P. E. Integrated land-use systems: Assessment of promising agroforest and alternative land-use practices to enhance carbon conservation and sequestration. Clim. Change 27, 71–92 (1994).
137. Houghton RA Byers B Nassikas AA A role for tropical forests in stabilizing atmospheric CO2 Nat. Clim. Change 2015 5 1022 1023 10.1038/nclimate2869
Houghton, R. A., Byers, B. & Nassikas, A. A. A role for tropical forests in stabilizing atmospheric CO2. Nat. Clim. Change 5, 1022–1023 (2015).
138. Humpenöder F Land-use and carbon cycle responses to moderate climate change: implications for land-based mitigation? Environ. Sci. Technol. 2015 49 6731 6739 10.1021/es506201r 25939014
Humpenöder, F. et al. Land-use and carbon cycle responses to moderate climate change: implications for land-based mitigation? Environ. Sci. Technol. 49, 6731–6739 (2015).25939014
139. Kreidenweis U Afforestation to mitigate climate change: impacts on food prices under consideration of albedo effects Environ. Res. Lett. 2016 11 085001 10.1088/1748-9326/11/8/085001
Kreidenweis, U. et al. Afforestation to mitigate climate change: impacts on food prices under consideration of albedo effects. Environ. Res. Lett. 11, 085001 (2016).
140. Lenton TM The potential for land-based biological CO2 removal to lower future atmospheric CO2 concentration Carbon Manag. 2010 1 145 160 10.4155/cmt.10.12
Lenton, T. M. The potential for land-based biological CO2 removal to lower future atmospheric CO2 concentration. Carbon Manag. 1, 145–160 (2010).
141. Moore JC Jevrejeva S Grinsted A Efficacy of geoengineering to limit 21st century sea-level rise Proc. Natl Acad. Sci. 2010 107 15699 15703 10.1073/pnas.1008153107 20798055
Moore, J. C., Jevrejeva, S. & Grinsted, A. Efficacy of geoengineering to limit 21st century sea-level rise. Proc. Natl Acad. Sci. 107, 15699–15703 (2010).20798055
142. Nilsson S Schopfhauser W The carbon-sequestration potential of a global afforestation program Clim. Change 1995 30 267 293 10.1007/BF01091928
Nilsson, S. & Schopfhauser, W. The carbon-sequestration potential of a global afforestation program. Clim. Change 30, 267–293 (1995).
143. Richards KR Stokes C A review of forest carbon sequestration cost studies: a dozen years of research Clim. Change 2004 63 1 48 10.1023/B:CLIM.0000018503.10080.89
Richards, K. R. & Stokes, C. A review of forest carbon sequestration cost studies: a dozen years of research. Clim. Change 63, 1–48 (2004).
144. Sohngen B Mendelsohn R An optimal control model of forest carbon sequestration Am. J. Agric. Econ. 2003 85 448 457 10.1111/1467-8276.00133
Sohngen, B. & Mendelsohn, R. An optimal control model of forest carbon sequestration. Am. J. Agric. Econ. 85, 448–457 (2003).
145. Strengers BJ Van Minnen JG Eickhout B The role of carbon plantations in mitigating climate change: potentials and costs Clim. Change 2008 88 343 366 10.1007/s10584-007-9334-4
Strengers, B. J., Van Minnen, J. G. & Eickhout, B. The role of carbon plantations in mitigating climate change: potentials and costs. Clim. Change 88, 343–366 (2008).
146. Thomson AM César Izaurralde R Smith SJ Clarke LE Integrated estimates of global terrestrial carbon sequestration Glob. Environ. Change 2008 18 192 203 10.1016/j.gloenvcha.2007.10.002
Thomson, A. M., César Izaurralde, R., Smith, S. J. & Clarke, L. E. Integrated estimates of global terrestrial carbon sequestration. Glob. Environ. Change 18, 192–203 (2008).
147. van Minnen JG Strengers BJ Eickhout B Swart RJ Leemans R Quantifying the effectiveness of climate change mitigation through forest plantations and carbon sequestration with an integrated land-use model Carbon Balance Manag. 2008 3 3 10.1186/1750-0680-3-3 18412946
van Minnen, J. G., Strengers, B. J., Eickhout, B., Swart, R. J. & Leemans, R. Quantifying the effectiveness of climate change mitigation through forest plantations and carbon sequestration with an integrated land-use model. Carbon Balance Manag. 3, 3 (2008).18412946
148. Winjum JK Dixon RK Schroeder PE Estimating the global potential of forest and agroforest management practices to sequester carbon Water Air Soil Pollut. 1992 64 213 227 10.1007/BF00477103
Winjum, J. K., Dixon, R. K. & Schroeder, P. E. Estimating the global potential of forest and agroforest management practices to sequester carbon. Water Air Soil Pollut. 64, 213–227 (1992).
149. Austin KG The economic costs of planting, preserving, and managing the world’s forests to mitigate climate change Nat. Commun. 2020 11 5946 10.1038/s41467-020-19578-z 33262324
Austin, K. G. et al. The economic costs of planting, preserving, and managing the world’s forests to mitigate climate change. Nat. Commun. 11, 5946 (2020).33262324
150. Dooley K Kartha S Land-based negative emissions: risks for climate mitigation and impacts on sustainable development Int Environ. Agreem. 2018 18 79 98 10.1007/s10784-017-9382-9
Dooley, K. & Kartha, S. Land-based negative emissions: risks for climate mitigation and impacts on sustainable development. Int Environ. Agreem. 18, 79–98 (2018).
151. Hasegawa T Fujimori S Ito A Takahashi K Careful selection of forest types in afforestation can increase carbon sequestration by 25% without compromising sustainability Commun. Earth Environ. 2024 5 1 10 10.1038/s43247-024-01336-4
Hasegawa, T., Fujimori, S., Ito, A. & Takahashi, K. Careful selection of forest types in afforestation can increase carbon sequestration by 25% without compromising sustainability. Commun. Earth Environ. 5, 1–10 (2024).
152. Lal R The carbon sequestration potential of terrestrial ecosystems J. Soil Water Conserv. 2018 73 145A 152A 10.2489/jswc.73.6.145A
Lal, R. et al. The carbon sequestration potential of terrestrial ecosystems. J. Soil Water Conserv. 73, 145A–152A (2018).
153. Lewis SL Wheeler CE Mitchard ETA Koch A Restoring natural forests is the best way to remove atmospheric carbon Nature 2019 568 25 28 10.1038/d41586-019-01026-8 30940972
Lewis, S. L., Wheeler, C. E., Mitchard, E. T. A. & Koch, A. Restoring natural forests is the best way to remove atmospheric carbon. Nature 568, 25–28 (2019).30940972
154. Littleton EW Dynamic modelling shows substantial contribution of ecosystem restoration to climate change mitigation Environ. Res. Lett. 2021 16 124061 10.1088/1748-9326/ac3c6c
Littleton, E. W. et al. Dynamic modelling shows substantial contribution of ecosystem restoration to climate change mitigation. Environ. Res. Lett. 16, 124061 (2021).
155. Nolan CJ Field CB Mach KJ Constraints and enablers for increasing carbon storage in the terrestrial biosphere Nat. Rev. Earth Environ. 2021 2 436 446 10.1038/s43017-021-00166-8
Nolan, C. J., Field, C. B. & Mach, K. J. Constraints and enablers for increasing carbon storage in the terrestrial biosphere. Nat. Rev. Earth Environ. 2, 436–446 (2021).
156. Walker WS The global potential for increased storage of carbon on land Proc. Natl Acad. Sci. 2022 119 e2111312119 10.1073/pnas.2111312119 35639697
Walker, W. S. et al. The global potential for increased storage of carbon on land. Proc. Natl Acad. Sci. 119, e2111312119 (2022).35639697
157. Weber J Chemistry-albedo feedbacks offset up to a third of forestation’s CO2 removal benefits Science 2024 383 860 864 10.1126/science.adg6196 38386743
Weber, J. et al. Chemistry-albedo feedbacks offset up to a third of forestation’s CO2 removal benefits. Science 383, 860–864 (2024).38386743
158. Wilkinson, K., Chissel, C. & Foley, J. The Drawdown Review: Climate Solutions for a New Decade. (Project Drawdown, 2020).
159. Caldecott, B., Lomax, G. & Workman, M. Stranded carbon assets and negative emissions technologies. Stranded Assets Programme (2015).
160. McLaren D A comparative global assessment of potential negative emissions technologies Process Saf. Environ. Prot. 2012 90 489 500 10.1016/j.psep.2012.10.005
McLaren, D. A comparative global assessment of potential negative emissions technologies. Process Saf. Environ. Prot. 90, 489–500 (2012).
161. Sathaye, J., Makundi, W., Dale, L., Chan, P. & Andrasko, K. GHG Mitigation potential, costs and benefits in global forests: a dynamic partial equilibrium approach. Energy J., 3, 127–162 (2006).
162. Smith P Biophysical and economic limits to negative CO2 emissions Nat. Clim. Change 2016 6 42 50 10.1038/nclimate2870
Smith, P. et al. Biophysical and economic limits to negative CO2 emissions. Nat. Clim. Change 6, 42–50 (2016).
163. Bathiany S Claussen M Brovkin V Raddatz T Gayler V Combined biogeophysical and biogeochemical effects of large-scale forest cover changes in the MPI earth system model Biogeosciences 2010 7 1383 1399 10.5194/bg-7-1383-2010
Bathiany, S., Claussen, M., Brovkin, V., Raddatz, T. & Gayler, V. Combined biogeophysical and biogeochemical effects of large-scale forest cover changes in the MPI earth system model. Biogeosciences 7, 1383–1399 (2010).
164. Jayakrishnan KU Bala G A comparison of the climate and carbon cycle effects of carbon removal by afforestation and an equivalent reduction in fossil fuel emissions Biogeosciences 2023 20 1863 1877 10.5194/bg-20-1863-2023
Jayakrishnan, K. U. & Bala, G. A comparison of the climate and carbon cycle effects of carbon removal by afforestation and an equivalent reduction in fossil fuel emissions. Biogeosciences 20, 1863–1877 (2023).
165. Swann ALS Fung IY Chiang JCH Mid-latitude afforestation shifts general circulation and tropical precipitation Proc. Natl Acad. Sci. 2012 109 712 716 10.1073/pnas.1116706108 22190490
Swann, A. L. S., Fung, I. Y. & Chiang, J. C. H. Mid-latitude afforestation shifts general circulation and tropical precipitation. Proc. Natl Acad. Sci. 109, 712–716 (2012).22190490
166. Prütz R Strefler J Rogelj J Fuss S Understanding the carbon dioxide removal range in 1.5 °C compatible and high overshoot pathways Environ. Res. Commun. 2023 5 041005 10.1088/2515-7620/accdba
Prütz, R., Strefler, J., Rogelj, J. & Fuss, S. Understanding the carbon dioxide removal range in 1.5 °C compatible and high overshoot pathways. Environ. Res. Commun. 5, 041005 (2023).
167. Kaufmann, D. & Kraay, A. Worldwide Governance Indicators, 2023 Update (www.govindicators.org) (2023).
168. Andrijevic M Crespo Cuaresma J Muttarak R Schleussner C-F Governance in socioeconomic pathways and its role for future adaptive capacity Nat. Sustain 2020 3 35 41 10.1038/s41893-019-0405-0
Andrijevic, M., Crespo Cuaresma, J., Muttarak, R. & Schleussner, C.-F. Governance in socioeconomic pathways and its role for future adaptive capacity. Nat. Sustain 3, 35–41 (2020).
169. Feyertag, J. et al. Prindex Comparative Report: A Global Assessment of Perceived Tenure Security from 140 Countries. https://www.prindex.net/reports/prindex-comparative-report-july-2020/ (2020).
170. Prindex. Prindex 2020 global dataset. (2020).
171. World Bank. Geospatial Poverty Portal. World Bank Group. (2023).
172. Dubertret, F. & Alden Wily, L. Percent of Indigenous and Community Lands. Data file from LandMark: The Global Platform of Indigenous and Community Lands. www.landmarkmap.org (2015).
173. WorldPop. WorldPop (www.worldpop.org—School of Geography and Environmental Science, University of Southampton; Department of Geography and Geosciences, University of Louisville; Departement de Geographie, Universite de Namur) and Center for International Earth Science Information Network (CIESIN), Columbia University). (2018).
174. Nabuurs, G.-J. et al. Agriculture, Forestry and Other Land Uses (AFOLU). in Climate change 2022: Mitigation of climate change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge University Press, Cambridge, UK and New York, 2022).
