Concurrent and lagged effects of spring greening on seasonal carbon gain and water loss across the Northern Hemisphere.

Spring greening has been widely observed across the Northern Hemisphere (NH) using a remotely sensed vegetation index (e.g., the normalized difference vegetation index, NDVI). However, there is still a debate on the ecological effects of spring greening on seasonal carbon and water budgets. This study jointly investigated the concurrent and lagged effects of spring greening on carbon gain (gross primary productivity, GPP) and water loss (evapotranspiration, ET) in the summer-active ecosystems at mid and high latitudes of NH using remote sensing and multimodel ensemble data during 1982-2013. The results showed that the collective promotion of spring greening to concurrent GPP and ET is widespread despite variations in magnitude and significance. Both beneficial and adverse lagged effects of spring greening on summer GPP commonly appear with an obvious spatial heterogeneity and difference among climate-plant types. However, the expected significant suppression of spring greening to summer GPP was rarely observed even in the areas where spring ET was significantly promoted by spring greening. Nevertheless, when drought was taken into account, the response patterns of spring water use to spring greening varied to some extent, and the adverse lagged effect of spring greening to summer GPP appeared or strengthened in some regions, especially during the years with dry summer. Given the predicted warming of the climate and more frequent climatic extremes, the adverse effect of spring greening should be given more attention.

Nonlinear dynamics of river runoff elucidated by horizontal visibility graphs.

Horizontal Visibility Graphs (HVGs) are a recently developed method to construct networks from time series. The values of the time series are considered as the nodes of the network and are linked to each other if there is no larger value between them, such as they can "see" each other. The network properties reflect the nonlinear dynamics of the time series. For some classes of stochastic processes and for periodic time series, analytical results can be obtained for network-derived quantities such as the degree distribution, the local clustering coefficient distribution, the mean path length, and others. HVGs have the potential to discern between deterministic-chaotic and correlated-stochastic time series. Here, we investigate the sensitivity of the HVG methodology to properties and pre-processing of real-world data, i.e., time series length, the presence of ties, and deseasonalization, using a set of around 150 runoff time series from managed rivers at daily resolution from Brazil with an average length of 65 years. We show that an application of HVGs on real-world time series requires a careful consideration of data pre-processing steps and analysis methodology before robust results and interpretations can be obtained. For example, one recent analysis of the degree distribution of runoff records reported pronounced sub-exponential "long-tailed" behavior of North American rivers, whereas another study of South American rivers showed hyper-exponential "short-tailed" behavior resembling correlated noise. We demonstrate, using the dataset of Brazilian rivers, that these apparently contradictory results can be reconciled by minor differences in data-preprocessing (here: small differences in subtracting the seasonal cycle). Hence, data-preprocessing that is conventional in hydrology ("deseasonalization") changes long-term correlations and the overall runoff dynamics substantially, and we present empirical consequences and extensive simulations to investigate these issues from a HVG methodological perspective. After carefully accounting for these methodological aspects, the HVG analysis reveals that the river runoff dataset shows indeed complex behavior that appears to stem from a superposition of short-term correlated noise and "long-tailed behaviour," i.e., highly connected nodes. Moreover, the construction of a dam along a river tends to increase short-term correlations in runoff series. In summary, the present study illustrates the (often substantial) effects of methodological and data-preprocessing choices for the interpretation of river runoff dynamics in the HVG framework and its general applicability for real-world time series.

Increasing extreme precipitation variability plays a key role in future record-shattering event probability.

Climate events that break records by large margins are a threat to society and ecosystems. Climate change is expected to increase the probability of such events, but quantifying these probabilities is challenging due to natural variability and limited data availability, especially for observations and very rare extremes. Here we estimate the probability of precipitation events that shatter records by a margin of at least one pre-industrial standard deviation. Using large ensemble climate simulations and extreme value theory, we determine empirical and analytical record shattering probabilities and find they are in high agreement. We show that, particularly in high emission scenarios, models project much higher record-shattering precipitation probabilities in a changing relative to a stationary climate by the end of the century for almost all the global land, with the strongest increases in vulnerable regions in the tropics. We demonstrate that increasing variability is an essential driver of near-term increases in record-shattering precipitation probability, and present a framework that quantifies the influence of combined trends in mean and variability on record-shattering behaviour in extreme precipitation. Probability estimates of record-shattering precipitation events in a warming world are crucial to inform risk assessment and adaptation policies.

Enhanced global carbon cycle sensitivity to tropical temperature linked to internal climate variability.

The sensitivity of atmospheric CO2 growth rate to tropical temperature (γT) has almost doubled between 1959 and 2011, a trend that has been linked to increasing drought in the tropics. However, γT has declined since then. Understanding whether these variations in γT reflect forced changes or internal climate variability in the carbon cycle is crucial for future climate projections. We show that doubling sensitivity events can arise in simulations by Earth system models with perturbed initial conditions but are likely explained by internal climate variability. We show that the doubling sensitivity event is associated with the occurrence of a few, but very strong, El Niño events, such as 1982/83 and 1997/98. Such extreme events result in concurrent carbon release by tropical and extratropical ecosystems, increasing the variance of the global land carbon sink and its apparent sensitivity to tropical temperature. Our results imply that the doubling sensitivity does not necessarily indicate a change in carbon cycle response to climate change.

Heat extremes in Western Europe increasing faster than simulated due to atmospheric circulation trends.

Over the last 70 years, extreme heat has been increasing at a disproportionate rate in Western Europe, compared to climate model simulations. This mismatch is not well understood. Here, we show that a substantial fraction (0.8 °C [0.2°-1.4 °C] of 3.4 °C per global warming degree) of the heat extremes trend is induced by atmospheric circulation changes, through more frequent southerly flows over Western Europe. In the 170 available simulations from 32 different models that we analyzed, including 3 large model ensembles, none have a circulation-induced heat trend as large as observed. This can be due to underestimated circulation response to external forcing, or to a systematic underestimation of low-frequency variability, or both. The former implies that future projections are too conservative, the latter that we are left with deep uncertainty regarding the pace of future summer heat in Europe. This calls for caution when interpreting climate projections of heat extremes over Western Europe, in view of adaptation to heat waves.

Climatic and soil factors explain the two-dimensional spectrum of global plant trait variation.

Plant functional traits can predict community assembly and ecosystem functioning and are thus widely used in global models of vegetation dynamics and land-climate feedbacks. Still, we lack a global understanding of how land and climate affect plant traits. A previous global analysis of six traits observed two main axes of variation: (1) size variation at the organ and plant level and (2) leaf economics balancing leaf persistence against plant growth potential. The orthogonality of these two axes suggests they are differently influenced by environmental drivers. We find that these axes persist in a global dataset of 17 traits across more than 20,000 species. We find a dominant joint effect of climate and soil on trait variation. Additional independent climate effects are also observed across most traits, whereas independent soil effects are almost exclusively observed for economics traits. Variation in size traits correlates well with a latitudinal gradient related to water or energy limitation. In contrast, variation in economics traits is better explained by interactions of climate with soil fertility. These findings have the potential to improve our understanding of biodiversity patterns and our predictions of climate change impacts on biogeochemical cycles.

Robust detection of forced warming in the presence of potentially large climate variability.

Climate warming is unequivocal and exceeds internal climate variability. However, estimates of the magnitude of decadal-scale variability from models and observations are uncertain, limiting determination of the fraction of warming attributable to external forcing. Here, we use statistical learning to extract a fingerprint of climate change that is robust to different model representations and magnitudes of internal variability. We find a best estimate forced warming trend of 0.8°C over the past 40 years, slightly larger than observed. It is extremely likely that at least 85% is attributable to external forcing based on the median variability across climate models. Detection remains robust even when evaluated against models with high variability and if decadal-scale variability were doubled. This work addresses a long-standing limitation in attributing warming to external forcing and opens up opportunities even in the case of large model differences in decadal-scale variability, model structural uncertainty, and limited observational records.

Past warming trend constrains future warming in CMIP6 models.

Future global warming estimates have been similar across past assessments, but several climate models of the latest Sixth Coupled Model Intercomparison Project (CMIP6) simulate much stronger warming, apparently inconsistent with past assessments. Here, we show that projected future warming is correlated with the simulated warming trend during recent decades across CMIP5 and CMIP6 models, enabling us to constrain future warming based on consistency with the observed warming. These findings carry important policy-relevant implications: The observationally constrained CMIP6 median warming in high emissions and ambitious mitigation scenarios is over 16 and 14% lower by 2050 compared to the raw CMIP6 median, respectively, and over 14 and 8% lower by 2090, relative to 1995-2014. Observationally constrained CMIP6 warming is consistent with previous assessments based on CMIP5 models, and in an ambitious mitigation scenario, the likely range is consistent with reaching the Paris Agreement target.

Half a degree and rapid socioeconomic development matter for heatwave risk.

While every society can be exposed to heatwaves, some people suffer far less harm and recover more quickly than others from their occurrence. Here we project indicators of global heatwave risk associated with global warming of 1.5 and 2 °C, specified by the Paris agreement, for two future pathways of societal development representing low and high vulnerability conditions. Results suggest that at the 1.5 °C warming level, heatwave exposure in 2075 estimated for the population living in low development countries is expected to be greater than exposure at the warming level of 2 °C for the population living in very high development countries. A similar result holds for an illustrative heatwave risk index. However, the projected difference in heatwave exposure and the illustrative risk index for the low and very high development countries will be significantly reduced if global warming is stabilized below 1.5 °C, and in the presence of rapid social development.

Extreme heat-related mortality avoided under Paris Agreement goals.

In key European cities, stabilizing climate warming at 1.5 °C would decrease extreme heat-related mortality by 15-22% per summer compared with stabilization at 2 °C.

Diagnosing the Dynamics of Observed and Simulated Ecosystem Gross Primary Productivity with Time Causal Information Theory Quantifiers.

Data analysis and model-data comparisons in the environmental sciences require diagnostic measures that quantify time series dynamics and structure, and are robust to noise in observational data. This paper investigates the temporal dynamics of environmental time series using measures quantifying their information content and complexity. The measures are used to classify natural processes on one hand, and to compare models with observations on the other. The present analysis focuses on the global carbon cycle as an area of research in which model-data integration and comparisons are key to improving our understanding of natural phenomena. We investigate the dynamics of observed and simulated time series of Gross Primary Productivity (GPP), a key variable in terrestrial ecosystems that quantifies ecosystem carbon uptake. However, the dynamics, patterns and magnitudes of GPP time series, both observed and simulated, vary substantially on different temporal and spatial scales. We demonstrate here that information content and complexity, or Information Theory Quantifiers (ITQ) for short, serve as robust and efficient data-analytical and model benchmarking tools for evaluating the temporal structure and dynamical properties of simulated or observed time series at various spatial scales. At continental scale, we compare GPP time series simulated with two models and an observations-based product. This analysis reveals qualitative differences between model evaluation based on ITQ compared to traditional model performance metrics, indicating that good model performance in terms of absolute or relative error does not imply that the dynamics of the observations is captured well. Furthermore, we show, using an ensemble of site-scale measurements obtained from the FLUXNET archive in the Mediterranean, that model-data or model-model mismatches as indicated by ITQ can be attributed to and interpreted as differences in the temporal structure of the respective ecological time series. At global scale, our understanding of C fluxes relies on the use of consistently applied land models. Here, we use ITQ to evaluate model structure: The measures are largely insensitive to climatic scenarios, land use and atmospheric gas concentrations used to drive them, but clearly separate the structure of 13 different land models taken from the CMIP5 archive and an observations-based product. In conclusion, diagnostic measures of this kind provide data-analytical tools that distinguish different types of natural processes based solely on their dynamics, and are thus highly suitable for environmental science applications such as model structural diagnostics.