Greening drylands despite warming consistent with carbon dioxide fertilization effect.

The rising atmospheric CO2 concentration leads to a CO2 fertilization effect on plants-that is, increased photosynthetic uptake of CO2 by leaves and enhanced water-use efficiency (WUE). Yet, the resulting net impact of CO2 fertilization on plant growth and soil moisture (SM) savings at large scale is poorly understood. Drylands provide a natural experimental setting to detect the CO2 fertilization effect on plant growth since foliage amount, plant water-use and photosynthesis are all tightly coupled in water-limited ecosystems. A long-term change in the response of leaf area index (LAI, a measure of foliage amount) to changes in SM is likely to stem from changing water demand of primary productivity in water-limited ecosystems and is a proxy for changes in WUE. Using 34-year satellite observations of LAI and SM over tropical and subtropical drylands, we identify that a 1% increment in SM leads to 0.15% (±0.008, 95% confidence interval) and 0.51% (±0.01, 95% confidence interval) increments in LAI during 1982-1998 and 1999-2015, respectively. The increasing response of LAI to SM has contributed 7.2% (±3.0%, 95% confidence interval) to total dryland greening during 1999-2015 compared to 1982-1998. The increasing response of LAI to SM is consistent with the CO2 fertilization effect on WUE in water-limited ecosystems, indicating that a given amount of SM has sustained greater amounts of photosynthetic foliage over time. The LAI responses to changes in SM from seven dynamic global vegetation models are not always consistent with observations, highlighting the need for improved process knowledge of terrestrial ecosystem responses to rising atmospheric CO2 concentration.

Comparison of forest above-ground biomass from dynamic global vegetation models with spatially explicit remotely sensed observation-based estimates.

Gaps in our current understanding and quantification of biomass carbon stocks, particularly in tropics, lead to large uncertainty in future projections of the terrestrial carbon balance. We use the recently published GlobBiomass data set of forest above-ground biomass (AGB) density for the year 2010, obtained from multiple remote sensing and in situ observations at 100 m spatial resolution to evaluate AGB estimated by nine dynamic global vegetation models (DGVMs). The global total forest AGB of the nine DGVMs is 365 ± 66 Pg C, the spread corresponding to the standard deviation between models, compared to 275 Pg C with an uncertainty of ~13.5% from GlobBiomass. Model-data discrepancy in total forest AGB can be attributed to their discrepancies in the AGB density and/or forest area. While DGVMs represent the global spatial gradients of AGB density reasonably well, they only have modest ability to reproduce the regional spatial gradients of AGB density at scales below 1000 km. The 95th percentile of AGB density (AGB95 ) in tropics can be considered as the potential maximum of AGB density which can be reached for a given annual precipitation. GlobBiomass data show local deficits of AGB density compared to the AGB95 , particularly in transitional and/or wet regions in tropics. We hypothesize that local human disturbances cause more AGB density deficits from GlobBiomass than from DGVMs, which rarely represent human disturbances. We then analyse empirical relationships between AGB density deficits and forest cover changes, population density, burned areas and livestock density. Regression analysis indicated that more than 40% of the spatial variance of AGB density deficits in South America and Africa can be explained; in Southeast Asia, these factors explain only ~25%. This result suggests TRENDY v6 DGVMs tend to underestimate biomass loss from diverse and widespread anthropogenic disturbances, and as a result overestimate turnover time in AGB.

State of the science in reconciling top-down and bottom-up approaches for terrestrial CO2 budget.

Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: "top-down" atmospheric inversions and "bottom-up" biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere-biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink-source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top-down and bottom-up approaches and yield more robust knowledge of regional CO2 budgets.

Simulating the recent impacts of multiple biotic disturbances on forest carbon cycling across the United States.

Biotic disturbances (BDs, for example, insects, pathogens, and wildlife herbivory) substantially affect boreal and temperate forest ecosystems globally. However, accurate impact assessments comprising larger spatial scales are lacking to date although these are critically needed given the expected disturbance intensification under a warming climate. Hence, our quantitative knowledge on current and future BD impacts, for example, on forest carbon (C) cycling, is strongly limited. We extended a dynamic global vegetation model to simulate ecosystem response to prescribed tree mortality and defoliation due to multiple biotic agents across United States forests during the period 1997-2015, and quantified the BD-induced vegetation C loss, that is, C fluxes from live vegetation to dead organic matter pools. Annual disturbance fractions separated by BD type (tree mortality and defoliation) and agent (bark beetles, defoliator insects, other insects, pathogens, and other biotic agents) were calculated at 0.5° resolution from aerial-surveyed data and applied within the model. Simulated BD-induced C fluxes totaled 251.6 Mt C (annual mean: 13.2 Mt C year-1 , SD ±7.3 Mt C year-1 between years) across the study domain, to which tree mortality contributed 95% and defoliation 5%. Among BD agents, bark beetles caused most C fluxes (61%), and total insect-induced C fluxes were about five times larger compared to non-insect agents, for example, pathogens and wildlife. Our findings further demonstrate that BD-induced C cycle impacts (i) displayed high spatio-temporal variability, (ii) were dominated by different agents across BD types and regions, and (iii) were comparable in magnitude to fire-induced impacts. This study provides the first ecosystem model-based assessment of BD-induced impacts on forest C cycling at the continental scale and going beyond single agent-host systems, thus allowing for comparisons across regions, BD types, and agents. Ultimately, a perspective on the potential and limitations of a more process-based incorporation of multiple BDs in ecosystem models is offered.

Adaptation of global land use and management intensity to changes in climate and atmospheric carbon dioxide.

Land use contributes to environmental change, but is also influenced by such changes. Climate and atmospheric carbon dioxide (CO2 ) levels' changes alter agricultural crop productivity, plant water requirements and irrigation water availability. The global food system needs to respond and adapt to these changes, for example, by altering agricultural practices, including the crop types or intensity of management, or shifting cultivated areas within and between countries. As impacts and associated adaptation responses are spatially specific, understanding the land use adaptation to environmental changes requires crop productivity representations that capture spatial variations. The impact of variation in management practices, including fertiliser and irrigation rates, also needs to be considered. To date, models of global land use have selected agricultural expansion or intensification levels using relatively aggregate spatial representations, typically at a regional level, that are not able to characterise the details of these spatially differentiated responses. Here, we show results from a novel global modelling approach using more detailed biophysically derived yield responses to inputs with greater spatial specificity than previously possible. The approach couples a dynamic global vegetative model (LPJ-GUESS) with a new land use and food system model (PLUMv2), with results benchmarked against historical land use change from 1970. Land use outcomes to 2100 were explored, suggesting that increased intensity of climate forcing reduces the inputs required for food production, due to the fertilisation and enhanced water use efficiency effects of elevated atmospheric CO2 concentrations, but requiring substantial shifts in the global and local patterns of production. The results suggest that adaptation in the global agriculture and food system has substantial capacity to diminish the negative impacts and gain greater benefits from positive outcomes of climate change. Consequently, agricultural expansion and intensification may be lower than found in previous studies where spatial details and processes consideration were more constrained.

Large uncertainty in carbon uptake potential of land-based climate-change mitigation efforts.

Most climate mitigation scenarios involve negative emissions, especially those that aim to limit global temperature increase to 2°C or less. However, the carbon uptake potential in land-based climate change mitigation efforts is highly uncertain. Here, we address this uncertainty by using two land-based mitigation scenarios from two land-use models (IMAGE and MAgPIE) as input to four dynamic global vegetation models (DGVMs; LPJ-GUESS, ORCHIDEE, JULES, LPJmL). Each of the four combinations of land-use models and mitigation scenarios aimed for a cumulative carbon uptake of ~130 GtC by the end of the century, achieved either via the cultivation of bioenergy crops combined with carbon capture and storage (BECCS) or avoided deforestation and afforestation (ADAFF). Results suggest large uncertainty in simulated future land demand and carbon uptake rates, depending on the assumptions related to land use and land management in the models. Total cumulative carbon uptake in the DGVMs is highly variable across mitigation scenarios, ranging between 19 and 130 GtC by year 2099. Only one out of the 16 combinations of mitigation scenarios and DGVMs achieves an equivalent or higher carbon uptake than achieved in the land-use models. The large differences in carbon uptake between the DGVMs and their discrepancy against the carbon uptake in IMAGE and MAgPIE are mainly due to different model assumptions regarding bioenergy crop yields and due to the simulation of soil carbon response to land-use change. Differences between land-use models and DGVMs regarding forest biomass and the rate of forest regrowth also have an impact, albeit smaller, on the results. Given the low confidence in simulated carbon uptake for a given land-based mitigation scenario, and that negative emissions simulated by the DGVMs are typically lower than assumed in scenarios consistent with the 2°C target, relying on negative emissions to mitigate climate change is a highly uncertain strategy.

Process-oriented analysis of dominant sources of uncertainty in the land carbon sink.

The observed global net land carbon sink is captured by current land models. All models agree that atmospheric CO2 and nitrogen deposition driven gains in carbon stocks are partially offset by climate and land-use and land-cover change (LULCC) losses. However, there is a lack of consensus in the partitioning of the sink between vegetation and soil, where models do not even agree on the direction of change in carbon stocks over the past 60 years. This uncertainty is driven by plant productivity, allocation, and turnover response to atmospheric CO2 (and to a smaller extent to LULCC), and the response of soil to LULCC (and to a lesser extent climate). Overall, differences in turnover explain ~70% of model spread in both vegetation and soil carbon changes. Further analysis of internal plant and soil (individual pools) cycling is needed to reduce uncertainty in the controlling processes behind the global land carbon sink.

Evaluating the Interplay Between Biophysical Processes and Leaf Area Changes in Land Surface Models.

Land Surface Models (LSMs) are essential to reproduce biophysical processes modulated by vegetation and to predict the future evolution of the land-climate system. To assess the performance of an ensemble of LSMs (JSBACH, JULES, ORCHIDEE, CLM, and LPJ-GUESS) a consistent set of land surface energy fluxes and leaf area index (LAI) has been generated. Relationships of interannual variations of modeled surface fluxes and LAI changes have been analyzed at global scale across climatological gradients and compared with those obtained from satellite-based products. Model-specific strengths and deficiencies were diagnosed for tree and grass biomes. Results show that the responses of grasses are generally well represented in models with respect to the observed interplay between turbulent fluxes and LAI, increasing the confidence on how the LAI-dependent partition of net radiation into latent and sensible heat are simulated. On the contrary, modeled forest responses are characterized by systematic bias in the relation between the year-to-year variability in LAI and net radiation in cold and temperate climates, ultimately affecting the amount of absorbed radiation due to LAI-related effects on surface albedo. In addition, for tree biomes, the relationships between LAI and turbulent fluxes appear to contradict the experimental evidences. The dominance of the transpiration-driven over the observed albedo-driven effects might suggest that LSMs have the incorrect balance of these two processes. Such mismatches shed light on the limitations of our current understanding and process representation of the vegetation control on the surface energy balance and help to identify critical areas for model improvement.