Search

The planetary boundaries framework defines a safe operating space for humanity. To date, these boundaries have mostly been investigated separately, and it is unclear whether breaching one boundary can lead to the transgression of another. By employing a dynamic global vegetation model, we systematically simulate the strength and direction of the effects of different transgression levels of the climate change boundary (using climate output from ten phase 6 of the Coupled Model Intercomparison Project models for CO _2 levels ranging from 350 ppm to 1000 ppm). We focus on climate change-induced shifts of Earth’s major forest biomes, the control variable for the land-system change boundary, both by the end of this century and, to account for the long-term legacy effect, by the end of the millennium. Our simulations show that while staying within the 350 ppm climate change boundary co-stabilizes the land-system change boundary, breaching it (>450 ppm) leads to critical transgression of the latter, with greater severity the higher the ppm level rises and the more time passes. Specifically, this involves a poleward treeline shift, boreal forest dieback (nearly completely within its current area under extreme climate scenarios), competitive expansion of temperate forest into today’s boreal zone, and a slight tropical forest extension. These interacting changes also affect other planetary boundaries (freshwater change and biosphere integrity) and provide feedback to the climate change boundary itself. Our quantitative process-based study highlights the need for interactions to be studied for a systemic operationalization of the planetary boundaries framework.

Future flood and drought risks have been predicted to transition from moderate to high levels at global warmings of 1.5 °C and 2.0 °C above pre-industrial levels, respectively. However, these results were obtained by approximating the equilibrium climate using transient simulations with steadily warming. This approach was recently criticised due to the warmer global land temperature and higher mean precipitation intensities of the transient climate in comparison with the equilibrium climate. Therefore, it is unclear whether floods and droughts projected under a transient climate can be systematically substituted for those occurring in an equilibrated climate. Here, by employing a large ensemble of global hydrological models (HMs) forced by global climate models, we assess the validity of estimating flood and drought characteristics under equilibrium climates from transient simulations. Differences in flood characteristics under transient and equilibrium climates could be largely ascribed to natural variability, indicating that the floods derived from a transient climate reasonably approximate the floods expected in an equally warm, equilibrated climate. By contrast, significant differences in drought intensity between transient and equilibrium climates were detected over a larger global land area than expected from natural variability. Despite the large differences among HMs in representing the low streamflow regime, we found that the drought intensities occurring under a transient climate may not validly represent the intensities in an equally warm equilibrated climate for approximately 6.7% of the global land area.

Limiting mean global warming to well below 2 °C will probably require substantial negative emissions (NEs) within the 21st century. To achieve these, bioenergy plantations with subsequent carbon capture and storage (BECCS) may have to be implemented at a large scale. Irrigation of these plantations might be necessary to increase the yield, which is likely to put further pressure on already stressed freshwater systems. Conversely, the potential of bioenergy plantations (BPs) dedicated to achieving NEs through CO _2 assimilation may be limited in regions with low freshwater availability. This paper provides a first-order quantification of the biophysical potentials of BECCS as a negative emission technology contribution to reaching the 1.5 °C warming target, as constrained by associated water availabilities and requirements. Using a global biosphere model, we analyze the availability of freshwater for irrigation of BPs designed to meet the projected NEs to fulfill the 1.5 °C target, spatially explicitly on areas not reserved for ecosystem conservation or agriculture. We take account of the simultaneous water demands for agriculture, industries, and households and also account for environmental flow requirements (EFRs) needed to safeguard aquatic ecosystems. Furthermore, we assess to what extent different forms of improved water management on the suggested BPs and on cropland may help to reduce the freshwater abstractions. Results indicate that global water withdrawals for irrigation of BPs range between ∼400 and ∼3000 km ^3 yr ^−1 , depending on the scenario and the conversion efficiency of the carbon capture and storage process. Consideration of EFRs reduces the NE potential significantly, but can partly be compensated for by improved on-field water management.

To support implementation of the Paris Agreement, the new HAPPI ensemble of 20 bias-corrected simulations of four climate models was used to drive two global hydrological models, WaterGAP and LPJmL, for assessing freshwater-related hazards and risks in worlds approximately 1.5 °C and 2 °C warmer than pre-industrial. Quasi-stationary HAPPI simulations are better suited than transient CMIP-like simulations for assessing hazards at the two targeted long-term global warming (GW) levels. We analyzed seven hydrological hazard indicators that characterize freshwater-related hazards for humans, freshwater biota and vegetation. Using a strict definition for significant differences, we identified for all but one indicator that areas with either significantly wetter or drier conditions (calculated as percent changes from 2006–2015) are smaller in the 1.5 °C world. For example, 7 day high flow is projected to increase significantly on 11% and 21% of the global land area at 1.5 °C and 2 °C, respectively. However, differences between hydrological hazards at the two GW levels are significant on less than 12% of the area. GW affects a larger area and more people by increases—rather than by decreases—of mean annual and 1-in-10 dry year streamflow, 7 day high flow, and groundwater recharge. The opposite is true for 7 day low flow, maximum snow storage, and soil moisture in the driest month of the growing period. Mean annual streamflow shows the lowest projected percent changes of all indicators. Among country groups, low income countries and lower middle income countries are most affected by decreased low flows and increased high flows, respectively, while high income countries are least affected by such changes. The incremental impact between 1.5 °C and 2 °C on high flows would be felt most by low income and lower middle income countries, the effect on soil moisture and low flows most by high income countries.

Global-scale hydrological models are routinely used to assess water scarcity, flood hazards and droughts worldwide. Recent efforts to incorporate anthropogenic activities in these models have enabled more realistic comparisons with observations. Here we evaluate simulations from an ensemble of six models participating in the second phase of the Inter-Sectoral Impact Model Inter-comparison Project (ISIMIP2a). We simulate monthly runoff in 40 catchments, spatially distributed across eight global hydrobelts. The performance of each model and the ensemble mean is examined with respect to their ability to replicate observed mean and extreme runoff under human-influenced conditions. Application of a novel integrated evaluation metric to quantify the models’ ability to simulate timeseries of monthly runoff suggests that the models generally perform better in the wetter equatorial and northern hydrobelts than in drier southern hydrobelts. When model outputs are temporally aggregated to assess mean annual and extreme runoff, the models perform better. Nevertheless, we find a general trend in the majority of models towards the overestimation of mean annual runoff and all indicators of upper and lower extreme runoff. The models struggle to capture the timing of the seasonal cycle, particularly in northern hydrobelts, while in southern hydrobelts the models struggle to reproduce the magnitude of the seasonal cycle. It is noteworthy that over all hydrological indicators, the ensemble mean fails to perform better than any individual model—a finding that challenges the commonly held perception that model ensemble estimates deliver superior performance over individual models. The study highlights the need for continued model development and improvement. It also suggests that caution should be taken when summarising the simulations from a model ensemble based upon its mean output.

Human impacts increasingly affect the global hydrological cycle and indeed dominate hydrological changes in some regions. Hydrologists have sought to identify the human-impact-induced hydrological variations via parameterizing anthropogenic water uses in global hydrological models (GHMs). The consequently increased model complexity is likely to introduce additional uncertainty among GHMs. Here, using four GHMs, between-model uncertainties are quantified in terms of the ratio of signal to noise (SNR) for average river flow during 1971–2000 simulated in two experiments, with representation of human impacts (VARSOC) and without (NOSOC). It is the first quantitative investigation of between-model uncertainty resulted from the inclusion of human impact parameterizations. Results show that the between-model uncertainties in terms of SNRs in the VARSOC annual flow are larger (about 2% for global and varied magnitude for different basins) than those in the NOSOC, which are particularly significant in most areas of Asia and northern areas to the Mediterranean Sea. The SNR differences are mostly negative (−20% to 5%, indicating higher uncertainty) for basin-averaged annual flow. The VARSOC high flow shows slightly lower uncertainties than NOSOC simulations, with SNR differences mostly ranging from −20% to 20%. The uncertainty differences between the two experiments are significantly related to the fraction of irrigation areas of basins. The large additional uncertainties in VARSOC simulations introduced by the inclusion of parameterizations of human impacts raise the urgent need of GHMs development regarding a better understanding of human impacts. Differences in the parameterizations of irrigation, reservoir regulation and water withdrawals are discussed towards potential directions of improvements for future GHM development. We also discuss the advantages of statistical approaches to reduce the between-model uncertainties, and the importance of calibration of GHMs for not only better performances of historical simulations but also more robust and confidential future projections of hydrological changes under a changing environment.

Global hydrological models (GHMs) have been applied to assess global flood hazards, but their capacity to capture the timing and amplitude of peak river discharge—which is crucial in flood simulations—has traditionally not been the focus of examination. Here we evaluate to what degree the choice of river routing scheme affects simulations of peak discharge and may help to provide better agreement with observations. To this end we use runoff and discharge simulations of nine GHMs forced by observational climate data (1971–2010) within the ISIMIP2a project. The runoff simulations were used as input for the global river routing model CaMa-Flood. The simulated daily discharge was compared to the discharge generated by each GHM using its native river routing scheme. For each GHM both versions of simulated discharge were compared to monthly and daily discharge observations from 1701 GRDC stations as a benchmark. CaMa-Flood routing shows a general reduction of peak river discharge and a delay of about two to three weeks in its occurrence, likely induced by the buffering capacity of floodplain reservoirs. For a majority of river basins, discharge produced by CaMa-Flood resulted in a better agreement with observations. In particular, maximum daily discharge was adjusted, with a multi-model averaged reduction in bias over about 2/3 of the analysed basin area. The increase in agreement was obtained in both managed and near-natural basins. Overall, this study demonstrates the importance of routing scheme choice in peak discharge simulation, where CaMa-Flood routing accounts for floodplain storage and backwater effects that are not represented in most GHMs. Our study provides important hints that an explicit parameterisation of these processes may be essential in future impact studies.

The insufficiency of water resources to meet the needs of food production is a pressing issue that is likely to increase in importance in the future. Improved understanding of historical developments can provide a basis for addressing future challenges. In this study we analyse how hydroclimatic variation, cropland expansion and evolving agricultural practices have influenced the potential for food self-sufficiency within the last century. We consider a food production unit (FPU) to have experienced green–blue water (GBW) scarcity if local renewable green (in soils) and blue water resources (in rivers, lakes, reservoirs, aquifers) were not sufficient for producing a reference food supply of 3000 kcal with 20% animal products for all inhabitants. The number of people living in FPUs affected by GBW scarcity has gone up from 360 million in 1905 (21% of world population at the time) to 2.2 billion (34%) in 2005. During this time, GBW scarcity has spread to large areas and become more frequent in regions where it occurs. Meanwhile, cropland expansion has increased green water availability for agriculture around the world, and advancements in agronomic practices have decreased water requirements of producing food. These efforts have improved food production potential and thus eased GBW scarcity considerably but also made possible the rapid population growth of the last century. The influence of modern agronomic practices is particularly striking: if agronomic practices of the early 1900s were applied today, it would roughly double the population under GBW scarcity worldwide.

Human land use and anthropogenic climate change (CC) are placing mounting pressure on natural ecosystems worldwide, with impacts on biodiversity, water resources, nutrient and carbon cycles. Here, we present a quantitative macro-scale comparative analysis of the separate and joint dual impacts of land use and land cover change (LULCC) and CC on the terrestrial biosphere during the last ca. 300 years, based on simulations with a dynamic global vegetation model and an aggregated metric of simultaneous biogeochemical, hydrological and vegetation-structural shifts. We find that by the beginning of the 21st century LULCC and CC have jointly caused major shifts on more than 90% of all areas now cultivated, corresponding to 26% of the land area. CC has exposed another 26% of natural ecosystems to moderate or major shifts. Within three centuries, the impact of LULCC on landscapes has increased 13-fold. Within just one century, CC effects have caught up with LULCC effects.

Climate warming affects permafrost soil carbon pools in two opposing ways: enhanced vegetation growth leads to higher carbon inputs to the soil, whereas permafrost melting accelerates decomposition and hence carbon release. Here, we study the spatial and temporal dynamics of these two processes under scenarios of climate change and evaluate their influence on the carbon balance of the permafrost zone. We use the dynamic global vegetation model LPJmL, which simulates plant physiological and ecological processes and includes a newly developed discrete layer energy balance permafrost module and a vertical carbon distribution within the soil layer. The model is able to reproduce the interactions between vegetation and soil carbon dynamics as well as to simulate dynamic permafrost changes resulting from changes in the climate. We find that vegetation responds more rapidly to warming of the permafrost zone than soil carbon pools due to long time lags in permafrost thawing, and that the initial simulated net uptake of carbon may continue for some decades of warming. However, once the turning point is reached, if carbon release exceeds uptake, carbon is lost irreversibly from the system and cannot be compensated for by increasing vegetation carbon input. Our analysis highlights the importance of including dynamic vegetation and long-term responses into analyses of permafrost zone carbon budgets.

In our globalizing world, the geographical locations of food production and consumption are becoming increasingly disconnected, which increases reliance on external resources and their trade. We quantified to what extent water and land constraints limit countries’ capacities, at present and by 2050, to produce on their own territory the crop products that they currently import from other countries. Scenarios of increased crop productivity and water use, cropland expansion (excluding areas prioritized for other uses) and population change are accounted for. We found that currently 16% of the world population use the opportunities of international trade to cover their demand for agricultural products. Population change may strongly increase the number of people depending on ex situ land and water resources up to about 5.2 billion (51% of world population) in the SRES A2r scenario. International trade will thus have to intensify if population growth is not accompanied by dietary change towards less resource-intensive products, by cropland expansion, or by productivity improvements, mainly in Africa and the Middle East. Up to 1.3 billion people may be at risk of food insecurity in 2050 in present low-income economies (mainly in Africa), if their economic development does not allow them to afford productivity increases, cropland expansion and/or imports from other countries.

This modelling study demonstrates at what level of global mean temperature rise (Δ T _g ) regions will be exposed to significant decreases of freshwater availability and changes to terrestrial ecosystems. Projections are based on a new, consistent set of 152 climate scenarios (eight Δ T _g trajectories reaching 1.5–5 ° C above pre-industrial levels by 2100, each scaled with spatial patterns from 19 general circulation models). The results suggest that already at a Δ T _g of 2 ° C and mainly in the subtropics, higher water scarcity would occur in >50% out of the 19 climate scenarios. Substantial biogeochemical and vegetation structural changes would also occur at 2 ° C, but mainly in subpolar and semiarid ecosystems. Other regions would be affected at higher Δ T _g levels, with lower intensity or with lower confidence. In total, mean global warming levels of 2 ° C, 3.5 ° C and 5 ° C are simulated to expose an additional 8%, 11% and 13% of the world population to new or aggravated water scarcity, respectively, with >50% confidence (while ∼1.3 billion people already live in water-scarce regions). Concurrently, substantial habitat transformations would occur in biogeographic regions that contain 1% (in zones affected at 2 ° C), 10% (3.5 ° C) and 74% (5 ° C) of present endemism-weighted vascular plant species, respectively. The results suggest nonlinear growth of impacts along with Δ T _g and highlight regional disparities in impact magnitudes and critical Δ T _g levels.