Implications of climate variability for the detection of multiple equilibria and for rapid transitions in the atmosphere-vegetation system

Paleoclimatic records indicate a decline of vegetation cover in the Western Sahara at the end of the African Humid Period (about 5,500 years before present). Modelling studies have shown that this phenomenon may be interpreted as a critical transition that results from a bifurcation in the atmosphere-vegetation system. However, the stability properties of this system are closely linked to climate variability and depend on the climate model and the methods of analysis. By coupling the Planet Simulator (PlaSim), an atmosphere model of intermediate complexity, with the simple dynamic vegetation model VECODE, we assess previous methods for the detection of multiple equilibria, and demonstrate their limitations. In particular, a stability diagram can yield misleading results because of spatial interactions, and the system's steady state and its dependency on initial conditions are affected by atmospheric variability and nonlinearities. In addition, we analyse the implications of climate variability for the abruptness of a vegetation decline. We find that a vegetation collapse can happen at different locations at different times. These collapses are possible despite large and uncorrelated climate variability. Because of the nonlinear relation between vegetation dynamics and precipitation the green state is initially stabilised by the high variability. When precipitation falls below a critical threshold, the desert state is stabilised as variability is then also decreased. © 2011 The Author(s)

Early warning signals of tipping points in periodically forced systems

This is the final version of the article. Available from the European Geosciences Union via the DOI in this record.The prospect of finding generic early warning signals of an approaching tipping point in a complex system has generated much interest recently. Existing methods are predicated on a separation of timescales between the system studied and its forcing. However, many systems, including several candidate tipping elements in the climate system, are forced periodically at a timescale comparable to their internal dynamics. Here we use alternative early warning signals of tipping points due to local bifurcations in systems subjected to periodic forcing whose timescale is similar to the period of the forcing. These systems are not in, or close to, a fixed point. Instead their steady state is described by a periodic attractor. For these systems, phase lag and amplification of the system response can provide early warning signals, based on a linear dynamics approximation. Furthermore, the Fourier spectrum of the system's time series reveals harmonics of the forcing period in the system response whose amplitude is related to how nonlinear the system's response is becoming with nonlinear effects becoming more prominent closer to a bifurcation. We apply these indicators as well as a return map analysis to a simple conceptual system and satellite observations of Arctic sea ice area, the latter conjectured to have a bifurcation type tipping point. We find no detectable signal of the Arctic sea ice approaching a local bifurcation.The research leading to these results has received funding from the European Union Seventh Framework Programme FP7/2007-2013 under grant agreement no. 603864 (HELIX). We are grateful to Peter Ashwin, Peter Cox, Michel Crucifix, Vasilis Dakos, Henk Dijkstra, Jan Sieber, Marten Scheffer and Appy Sluijs for the fruitful discussions over beers and balls

Abrupt Climate Change in an Oscillating World.

This is the final version of the article. Available from Nature Publishing Group via the DOI in this record.The notion that small changes can have large consequences in the climate or ecosystems has become popular as the concept of tipping points. Typically, tipping points are thought to arise from a loss of stability of an equilibrium when external conditions are slowly varied. However, this appealingly simple view puts us on the wrong foot for understanding a range of abrupt transitions in the climate or ecosystems because complex environmental systems are never in equilibrium. In particular, they are forced by diurnal variations, the seasons, Milankovitch cycles and internal climate oscillations. Here we show how abrupt and sometimes even irreversible change may be evoked by even small shifts in the amplitude or time scale of such environmental oscillations. By using model simulations and reconciling evidence from previous studies we illustrate how these phenomena can be relevant for ecosystems and elements of the climate system including terrestrial ecosystems, Arctic sea ice and monsoons. Although the systems we address are very different and span a broad range of time scales, the phenomena can be understood in a common framework that can help clarify and unify the interpretation of abrupt shifts in the Earth system.This work was carried out under the program of the Netherlands Earth System Science Centre (NESSC), financially supported by the Ministry of Education, Culture and Science (OCW). We are grateful to Chris Huntingford for his constructive comments that helped us to improve the manuscript. We would also like to acknowledge Michel Crucifix, Henk Dijkstra, and Peter Cox for their helpful comments. S.B. is eternally grateful to Nina Engelhardt and the University of Edinburgh for the inspiring working conditions

Statistical indicators of Arctic sea-ice stability-prospects and limitations

This is the final version of the article. Available from the European Geosciences Union via the DOI in this record.We examine the relationship between the mean and the variability of Arctic sea-ice coverage and volume in a large range of climates from globally ice-covered to globally ice-free conditions. Using a hierarchy of two column models and several comprehensive Earth system models, we consolidate the results of earlier studies and show that mechanisms found in simple models also dominate the interannual variability of Arctic sea ice in complex models. In contrast to predictions based on very idealised dynamical systems, we find a consistent and robust decrease of variance and autocorrelation of sea-ice volume before summer sea ice is lost. We attribute this to the fact that thinner ice can adjust more quickly to perturbations. Thereafter, the autocorrelation increases, mainly because it becomes dominated by the ocean water's large heat capacity when the ice-free season becomes longer. We show that these changes are robust to the nature and origin of climate variability in the models and do not depend on whether Arctic sea-ice loss occurs abruptly or irreversibly. We also show that our climate is changing too rapidly to detect reliable changes in autocorrelation of annual time series. Based on these results, the prospects of detecting statistical early warning signals before an abrupt sea-ice loss at a "tipping point" seem very limited. However, the robust relation between state and variability can be useful to build simple stochastic climate models and to make inferences about past and future sea-ice variability from only short observations or reconstructions.This work was carried out under the programme of the Netherlands Earth System Science Centre (NESSC), financially supported by the Ministry of Education, Culture and Science (OCW). We also acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modelling groups for producing and making available their model output. We thank Vasilis Dakos for helping to apply his early warnings R package and Chao Li for making available the MPI-ESM model output. S. B. gratefully acknowledges Arie Staal for his fruitful and revealing approaches to savour scientific achievements. We are also indebted to Till Wagner and Ian Eisenman for their valuable comments and their very amiable and cooperative spirit. Finally, we acknowledge two anonymous reviewers who helped us to improve the manuscript

Differentiable Programming for Earth System Modeling

Earth System Models (ESMs) are the primary tools for investigating future Earth system states at time scales from decades to centuries, especially in response to anthropogenic greenhouse gas release. State-of-the-art ESMs can reproduce the observational global mean temperature anomalies of the last 150 years. Nevertheless, ESMs need further improvements, most importantly regarding (i) the large spread in their estimates of climate sensitivity, i.e., the temperature response to increases in atmospheric greenhouse gases, (ii) the modeled spatial patterns of key variables such as temperature and precipitation, (iii) their representation of extreme weather events, and (iv) their representation of multistable Earth system components and their ability to predict associated abrupt transitions. Here, we argue that making ESMs automatically differentiable has huge potential to advance ESMs, especially with respect to these key shortcomings. First, automatic differentiability would allow objective calibration of ESMs, i.e., the selection of optimal values with respect to a cost function for a large number of free parameters, which are currently tuned mostly manually. Second, recent advances in Machine Learning (ML) and in the amount, accuracy, and resolution of observational data promise to be helpful with at least some of the above aspects because ML may be used to incorporate additional information from observations into ESMs. Automatic differentiability is an essential ingredient in the construction of such hybrid models, combining process-based ESMs with ML components. We document recent work showcasing the potential of automatic differentiation for a new generation of substantially improved, data-informed ESMs.Comment: 17 pages, 2 figure

Climatic impacts of land-use change due to crop yield increases and a universal carbon tax from a scenario model*

Future land cover will have a significant impact on climate and is strongly influenced by the extent of agricultural land use. Differing assumptions of crop yield increase and carbon pricing mitigation strategies affect projected expansion of agricultural land in future scenarios. In the representative concentration pathway 4.5 (RCP4.5) from phase 5 of the Coupled Model Intercomparison Project (CMIP5), the carbon effects of these land cover changes are included, although the biogeophysical effects are not. The afforestation in RCP4.5 has important biogeophysical impacts on climate, in addition to the land carbon changes, which are directly related to the assumption of crop yield increase and the universal carbon tax. To investigate the biogeophysical climatic impact of combinations of agricultural crop yield increases and carbon pricing mitigation, five scenarios of land-use change based on RCP4.5 are used as inputs to an earth system model [Hadley Centre Global Environment Model, version 2-Earth System (HadGEM2-ES)]. In the scenario with the greatest increase in agricultural land (as a result of no increase in crop yield and no climate mitigation) there is a significant -0.49 K worldwide cooling by 2100 compared to a control scenario with no land-use change. Regional cooling is up to -2.2 K annually in northeastern Asia. Including carbon feedbacks from the land-use change gives a small global cooling of -0.067 K. This work shows that there are significant impacts from biogeophysical land-use changes caused by assumptions of crop yield and carbon mitigation, which mean that land carbon is not the whole story. It also elucidates the potential conflict between cooling from biogeophysical climate effects of land-use change and wider environmental aims

Conditional diffusion models for downscaling & bias correction of Earth system model precipitation

Climate change exacerbates extreme weather events like heavy rainfall and flooding. As these events cause severe losses of property and lives, accurate high-resolution simulation of precipitation is imperative. However, existing Earth System Models (ESMs) struggle with resolving small-scale dynamics and suffer from biases, especially for extreme events. Traditional statistical bias correction and downscaling methods fall short in improving spatial structure, while recent deep learning methods lack controllability over the output and suffer from unstable training. Here, we propose a novel machine learning framework for simultaneous bias correction and downscaling. We train a generative diffusion model in a supervised way purely on observational data. We map observational and ESM data to a shared embedding space, where both are unbiased towards each other and train a conditional diffusion model to reverse the mapping. Our method can be used to correct any ESM field, as the training is independent of the ESM. Our approach ensures statistical fidelity, preserves large-scale spatial patterns and outperforms existing methods especially regarding extreme events and small-scale spatial features that are crucial for impact assessments