Global cotton production under climate change – Implications for yield and water consumption

Being an extensively produced natural fiber on earth, cotton is of importance for economies. Although the plant is broadly adapted to varying environments, the growth of and irrigation water demand on cotton may be challenged by future climate change. To study the impacts of climate change on cotton productivity in different regions across the world and the irrigation water requirements related to it, we use the process-based, spatially detailed biosphere and hydrology model LPJmL (Lund Potsdam Jena managed land). We find our modeled cotton yield levels in good agreement with reported values and simulated water consumption of cotton production similar to published estimates. Following the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) protocol, we employ an ensemble of five general circulation models under four representative concentration pathways (RCPs) for the 2011 2099 period to simulate future cotton yields. We find that irrigated cotton production does not suffer from climate change if CO2 effects are considered, whereas rainfed production is more sensitive to varying climate conditions. Considering the overall effect of a changing climate and CO2 fertilization, cotton production on current cropland steadily increases for most of the RCPs. Starting from _ 65 million tonnes in 2010, cotton production for RCP4.5 and RCP6.0 equates to 83 and 92 million tonnes at the end of the century, respectively. Under RCP8.5, simulated global cotton production rises by more than 50% by 2099. Taking only climate change into account, projected cotton production considerably shrinks in most scenarios, by up to one-Third or 43 million tonnes under RCP8.5. The simulation of future virtual water content (VWC) of cotton grown under elevated CO2 results for all scenarios in less VWC compared to ambient CO2 conditions. Under RCP6.0 and RCP8.5, VWC is notably decreased by more than 2000m3 t1 in areas where cotton is produced under purely rainfed conditions. By 2040, the average global VWC for cotton declines in all scenarios from currently 3300 to 3000m3 t1, and reduction continues by up to 30% in 2100 under RCP8.5. While the VWC decreases by the CO2 effect, elevated temperature acts in the opposite direction. Ignoring beneficial CO2 effects, global VWC of cotton would increase for all RCPs except RCP2.6, reaching more than 5000m3 t1 by the end of the simulation period under RCP8.5. Given the economic relevance of cotton production, climate change poses an additional stress and deserves special attention. Changes in VWC and water demands for cotton production are of special importance, as cotton production is known for its intense water consumption. The implications of climate impacts on cotton production on the one hand and the impact of cotton production on water resources on the other hand illustrate the need to assess how future climate change may affect cotton production and its resource requirements. Our results should be regarded as optimistic, because of high uncertainty with respect to CO2 fertilization and the lack of implementing processes of boll abscission under heat stress. Still, the inclusion of cotton in LPJmL allows for various large-scale studies to assess impacts of climate change on hydrological factors and the implications for agricultural production and carbon sequestration. © 2021 BMJ Publishing Group. All rights reserved

Future tree survival in European forests depends on understorey tree diversity

Climate change heavily threatens forest ecosystems worldwide and there is urgent need to understand what controls tree survival and forests stability. There is evidence that biodiversity can enhance ecosystem stability (Loreau and de Mazancourt in Ecol Lett 16:106–115, 2013; McCann in Nature 405:228–233, 2000), however it remains largely unclear whether this also holds for climate change and what aspects of biodiversity might be most important. Here we apply machine learning to outputs of a flexible-trait Dynamic Global Vegetation Model to unravel the effects of enhanced functional tree trait diversity and its sub-components on climate-change resistance of temperate forests (http://www.pik-potsdam.de/~billing/video/Forest_Resistance_LPJmLFIT.mp4). We find that functional tree trait diversity enhances forest resistance. We explain this with 1. stronger complementarity effects (~ 25% importance) especially improving the survival of trees in the understorey of up to + 16.8% (± 1.6%) and 2. environmental and competitive filtering of trees better adapted to future climate (40–87% importance). We conclude that forests containing functionally diverse trees better resist and adapt to future conditions. In this context, we especially highlight the role of functionally diverse understorey trees as they provide the fundament for better survival of young trees and filtering of resistant tree individuals in the future

Constraining modelled global vegetation dynamics and carbon turnover using multiple satellite observations

The response of land ecosystems to future climate change is among the largest unknowns in the global climate-carbon cycle feedback. This uncertainty originates from how dynamic global vegetation models (DGVMs) simulate climate impacts on changes in vegetation distribution, productivity, biomass allocation, and carbon turnover. The present-day availability of a multitude of satellite observations can potentially help to constrain DGVM simulations within model-data integration frameworks. Here, we use satellite-derived datasets of the fraction of absorbed photosynthetic active radiation (FAPAR), sun-induced fluorescence (SIF), above-ground biomass of trees (AGB), land cover, and burned area to constrain parameters for phenology, productivity, and vegetation dynamics in the LPJmL4 DGVM. Both the prior and the optimized model accurately reproduce present-day estimates of the land carbon cycle and of temporal dynamics in FAPAR, SIF and gross primary production. However, the optimized model reproduces better the observed spatial patterns of biomass, tree cover, and regional forest carbon turnover. Using a machine learning approach, we found that remaining errors in simulated forest carbon turnover can be explained with bioclimatic variables. This demonstrates the need to improve model formulations for climate effects on vegetation turnover and mortality despite the apparent successful constraint of simulated vegetation dynamics with multiple satellite observations

Improving the LPJmL4-SPITFIRE vegetation–fire model for South America using satellite data

Vegetation fires influence global vegetation distribution, ecosystem functioning, and global carbon cycling. Specifically in South America, changes in fire occurrence together with land-use change accelerate ecosystem fragmentation and increase the vulnerability of tropical forests and savannas to climate change. Dynamic global vegetation models (DGVMs) are valuable tools to estimate the effects of fire on ecosystem functioning and carbon cycling under future climate changes. However, most fire-enabled DGVMs have problems in capturing the magnitude, spatial patterns, and temporal dynamics of burned area as observed by satellites. As fire is controlled by the interplay of weather conditions, vegetation properties, and human activities, fire modules in DGVMs can be improved in various aspects. In this study we focus on improving the controls of climate and hence fuel moisture content on fire danger in the LPJmL4-SPITFIRE DGVM in South America, especially for the Brazilian fire-prone biomes of Caatinga and Cerrado. We therefore test two alternative model formulations (standard Nesterov Index and a newly implemented water vapor pressure deficit) for climate effects on fire danger within a formal model–data integration setup where we estimate model parameters against satellite datasets of burned area (GFED4) and aboveground biomass of trees. Our results show that the optimized model improves the representation of spatial patterns and the seasonal to interannual dynamics of burned area especially in the Cerrado and Caatinga regions. In addition, the model improves the simulation of aboveground biomass and the spatial distribution of plant functional types (PFTs). We obtained the best results by using the water vapor pressure deficit (VPD) for the calculation of fire danger. The VPD includes, in comparison to the Nesterov Index, a representation of the air humidity and the vegetation density. This work shows the successful application of a systematic model–data integration setup, as well as the integration of a new fire danger formulation, in order to optimize a process-based fire-enabled DGVM. It further highlights the potential of this approach to achieve a new level of accuracy in comprehensive global fire modeling and prediction

Simulating the effect of tillage practices with the global ecosystem model LPJmL (version 5.0-tillage)

The effects of tillage on soil properties, crop productivity, and global greenhouse gas emissions have been discussed in the last decades. Global ecosystem models have limited capacity to simulate the various effects of tillage. With respect to the decomposition of soil organic matter, they either assume a constant increase due to tillage or they ignore the effects of tillage. Hence, they do not allow for analysing the effects of tillage and cannot evaluate, for example, reduced tillage or no tillage (referred to here as “no-till”) practises as mitigation practices for climate change. In this paper, we describe the implementation of tillage-related practices in the global ecosystem model LPJmL. The extended model is evaluated against reported differences between tillage and no-till management on several soil properties. To this end, simulation results are compared with published meta-analyses on tillage effects. In general, the model is able to reproduce observed tillage effects on global, as well as regional, patterns of carbon and water fluxes. However, modelled N fluxes deviate from the literature values and need further study. The addition of the tillage module to LPJmL5 opens up opportunities to assess the impact of agricultural soil management practices under different scenarios with implications for agricultural productivity, carbon sequestration, greenhouse gas emissions, and other environmental indicators

Climate-induced hysteresis of the tropical forest in a fire-enabled Earth system model

Tropical rainforests are recognized as one of the terrestrial tipping elements which could have profound impacts on the global climate, once their vegetation has transitioned into savanna or grassland states. While several studies investigated the savannization of, e.g., the Amazon rainforest, few studies considered the influence of fire. Fire is expected to potentially shift the savanna-forest boundary and hence impact the dynamical equilibrium between these two possible vegetation states under changing climate. To investigate the climate-induced hysteresis in pan-tropical forests and the impact of fire under future climate conditions, we employed the Earth system model CM2Mc, which is biophysically coupled to the fire-enabled state-of-the-art dynamic global vegetation model LPJmL. We conducted several simulation experiments where atmospheric CO2 concentrations increased (impact phase) and decreased from the new state (recovery phase), each with and without enabling wildfires. We find a hysteresis of the biomass and vegetation cover in tropical forest systems, with a strong regional heterogeneity. After biomass loss along increasing atmospheric CO2 concentrations and accompanied mean surface temperature increase of about 4 ∘C (impact phase), the system does not recover completely into its original state on its return path, even though atmospheric CO2 concentrations return to their original state. While not detecting large-scale tipping points, our results show a climate-induced hysteresis in tropical forest and lagged responses in forest recovery after the climate has returned to its original state. Wildfires slightly widen the climate-induced hysteresis in tropical forests and lead to a lagged response in forest recovery by ca. 30 years

CM2Mc-LPJmL v1.0: biophysical coupling of a process-based dynamic vegetation model with managed land to a general circulation model

The terrestrial biosphere is exposed to land-use and climate change, which not only affects vegetation dynamics but also changes land–atmosphere feedbacks. Specifically, changes in land cover affect biophysical feedbacks of water and energy, thereby contributing to climate change. In this study, we couple the well-established and comprehensively validated dynamic global vegetation model LPJmL5 (Lund–Potsdam–Jena managed Land) to the coupled climate model CM2Mc, the latter of which is based on the atmosphere model AM2 and the ocean model MOM5 (Modular Ocean Model 5), and name it CM2Mc-LPJmL. In CM2Mc, we replace the simple land-surface model LaD (Land Dynamics; where vegetation is static and prescribed) with LPJmL5, and we fully couple the water and energy cycles using the Geophysical Fluid Dynamics Laboratory (GFDL) Flexible Modeling System (FMS). Several improvements to LPJmL5 were implemented to allow a fully functional biophysical coupling. These include a sub-daily cycle for calculating energy and water fluxes, conductance of the soil evaporation and plant interception, canopy-layer humidity, and the surface energy balance in order to calculate the surface and canopy-layer temperature within LPJmL5. Exchanging LaD with LPJmL5 and, therefore, switching from a static and prescribed vegetation to a dynamic vegetation allows us to model important biospheric processes, including fire, mortality, permafrost, hydrological cycling and the impacts of managed land (crop growth and irrigation). Our results show that CM2Mc-LPJmL has similar temperature and precipitation biases to the original CM2Mc model with LaD. The performance of LPJmL5 in the coupled system compared to Earth observation data and to LPJmL offline simulation results is within acceptable error margins. The historical global mean temperature evolution of our model setup is within the range of CMIP5 (Coupled Model Intercomparison Project Phase 5) models. The comparison of model runs with and without land-use change shows a partially warmer and drier climate state across the global land surface. CM2Mc-LPJmL opens new opportunities to investigate important biophysical vegetation–climate feedbacks with a state-of-the-art and process-based dynamic vegetation model

Drivers and patterns of land biosphere carbon balance reversal

The carbon balance of the land biosphere is the result of complex interactions between land, atmosphere and oceans, including climatic change, carbon dioxide fertilization and land-use change. While the land biosphere currently absorbs carbon dioxide from the atmosphere, this carbon balance might be reversed under climate and land-use change (‘carbon balance reversal’). A carbon balance reversal would render climate mitigation much more difficult, as net negative emissions would be needed to even stabilize atmospheric carbon dioxide concentrations. We investigate the robustness of the land biosphere carbon sink under different socio-economic pathways by systematically varying climate sensitivity, spatial patterns of climate change and resulting land-use changes. For this, we employ a modelling framework designed to account for all relevant feedback mechanisms by coupling the integrated assessment model IMAGE with the process-based dynamic vegetation, hydrology and crop growth model LPJmL. We find that carbon balance reversal can occur under a broad range of forcings and is connected to changes in tree cover and soil carbon mainly in northern latitudes. These changes are largely a consequence of vegetation responses to varying climate and only partially of land-use change and the rate of climate change. Spatial patterns of climate change as deduced from different climate models, substantially determine how much pressure in terms of global warming and land-use change the land biosphere will tolerate before the carbon balance is reversed. A reversal of the land biosphere carbon balance can occur as early as 2030, although at very low probability, and should be considered in the design of so-called peak-and-decline strategies.Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347Peer Reviewe

Variable tree rooting strategies are key for modelling the distribution, productivity and evapotranspiration of tropical evergreen forests

A variety of modelling studies have suggested tree rooting depth as a key variable to explain evapotranspiration rates, productivity and the geographical distribution of evergreen forests in tropical South America. However, none of those studies have acknowledged resource investment, timing and physical constraints of tree rooting depth within a competitive environment, undermining the ecological realism of their results. Here, we present an approach of implementing variable rooting strategies and dynamic root growth into the LPJmL4.0 (Lund-Potsdam-Jena managed Land) dynamic global vegetation model (DGVM) and apply it to tropical and sub-tropical South America under contemporary climate conditions. We show how competing rooting strategies which underlie the trade-off between above- and below-ground carbon investment lead to more realistic simulation of intra-annual productivity and evapotranspiration and consequently of forest cover and spatial biomass distribution. We find that climate and soil depth determine a spatially heterogeneous pattern of mean rooting depth and below-ground biomass across the study region. Our findings support the hypothesis that the ability of evergreen trees to adjust their rooting systems to seasonally dry climates is crucial to explaining the current dominance, productivity and evapotranspiration of evergreen forests in tropical South America.Fil: Sakschewski, Boris. Potsdam Institute for Climate Impact Research; AlemaniaFil: Von Bloh, Werner. Humboldt-Universität zu Berlin; AlemaniaFil: Drüke, Markus. Humboldt-Universität zu Berlin; AlemaniaFil: Sörensson, Anna. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Ruscica, Romina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Langerwisch, Fanny. Universitat Potsdam; AlemaniaFil: Billing, Maik. Universidade Federal de Santa Catarina; BrasilFil: Bereswill, Sarah. Universidade Estadual de Campinas; BrasilFil: Hirota, Marina. Potsdam Institute for Climate Impact Research; AlemaniaFil: Oliveira, Rafael Silva. Potsdam Institute for Climate Impact Research; AlemaniaFil: Heinke, Jens. Potsdam Institute for Climate Impact Research; AlemaniaFil: Thonicke, Kirsten. Potsdam Institute for Climate Impact Research; Alemani

Variable tree rooting strategies are key for modelling the distribution, productivity and evapotranspiration of tropical evergreen forests

A variety of modelling studies have suggested tree rooting depth as a key variable to explain evapotranspiration rates, productivity and the geographical distribution of evergreen forests in tropical South America. However, none of those studies have acknowledged resource investment, timing and physical constraints of tree rooting depth within a competitive environment, undermining the ecological realism of their results. Here, we present an approach of implementing variable rooting strategies and dynamic root growth into the LPJmL4.0 (Lund-Potsdam-Jena managed Land) dynamic global vegetation model (DGVM) and apply it to tropical and sub-tropical South America under contemporary climate conditions. We show how competing rooting strategies which underlie the trade-off between above- and below-ground carbon investment lead to more realistic simulation of intra-annual productivity and evapotranspiration and consequently of forest cover and spatial biomass distribution. We find that climate and soil depth determine a spatially heterogeneous pattern of mean rooting depth and below-ground biomass across the study region. Our findings support the hypothesis that the ability of evergreen trees to adjust their rooting systems to seasonally dry climates is crucial to explaining the current dominance, productivity and evapotranspiration of evergreen forests in tropical South America