Researcher profile

Frank Lunkeit

Frank Lunkeit contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.ao-ph cond-mat.stat-mech physics.flu-dyn physics.geo-ph physics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging

7works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Predicting Climate Change using Response Theory: Global Averages and Spatial Patterns

The provision of accurate methods for predicting the climate response to anthropogenic and natural forcings is a key contemporary scientific challenge. Using a simplified and efficient open-source general circulation model of the atmosphere featuring O($10^5$) degrees of freedom, we show how it is possible to approach such a problem using nonequilibrium statistical mechanics. Response theory allows one to practically compute the time-dependent measure supported on the pullback attractor of the climate system, whose dynamics is non-autonomous as a result of time-dependent forcings. We propose a simple yet efficient method for predicting - at any lead time and in an ensemble sense - the change in climate properties resulting from increase in the concentration of CO$_2$ using test perturbation model runs. We assess strengths and limitations of the response theory in predicting the changes in the globally averaged values of surface temperature and of the yearly total precipitation, as well as in their spatial patterns. The quality of the predictions obtained for the surface temperature fields is rather good, while in the case of precipitation a good skill is observed only for the global average. We also show how it is possible to define accurately concepts like the the inertia of the climate system or to predict when climate change is detectable given a scenario of forcing. Our analysis can be extended for dealing with more complex portfolios of forcings and can be adapted to treat, in principle, any climate observable. Our conclusion is that climate change is indeed a problem that can be effectively seen through a statistical mechanical lens, and that there is great potential for optimizing the current coordinated modelling exercises run for the preparation of the subsequent reports of the Intergovernmental Panel for Climate Change.

preprint2016arXiv

The role of atmospheric greenhouse gases, orbital parameters, and southern ocean gateways: an idealized model study

A set of idealized experiments are performed to analyze the competing effects of declining atmospheric CO2 concentrations, the opening of an ocean gateway, and varying orbital parameters. These forcing mechanisms, which influence the global mean climate state, may have played a role for triggering climate transitions of the past (for example during the Eocene-Oligocene climate transition and the build-up of the Antarctic Ice Sheet). Sensitivity simulations with a coupled atmosphere-ocean general circulation model are conducted to test these three forcings and their roles for the global climate. The simulations are carried out under idealized conditions to focus on the essentials. The combination of all three forcings triggers a climate transition which resembles the onset of the Antarctic glaciation. In particular, the temperatures in the southern high latitudes decrease and snow accumulates constantly. Moreover, the relative importance of each possible forcing is explored. All three of the mechanisms (atmospheric CO2 decrease, opening of the ocean gateways, and changing orbital parameters) cool the climate of the southern polar continent. The change induced by the orbital parameters depends strongly on the magnitude of the change in obliquity.

preprint2015arXiv

Global instability in the Ghil-Sellers model

The Ghil-Sellers energy balance model of Earth's climate, features -- for a considerable range of the solar intensity -- two stable climate states (a warm and a cold snowball Earth), where the bistability results from the celebrated ice-albedo feedback. The unstable solution is obtained and characterized in this paper. We find such unstable states by applying for the first time in a geophysical context the so-called edge tracking method that has been used for studying multiple coexisting states in shear flows. We examine robustness, efficiency, and accuracy properties of the edge tracking algorithm. We find that the procedure is the most efficient when taking a single bisection per cycle. Due to the strong diffusivity of the system trajectories of transient dynamics, initialized between the stable states with respect to the mean temperature, are confined to the heteroclininc trajectory, one which connects the fixed unstable and stable states, after relatively short transient times. This constraint dictates a functional relationship between observables. We characterize such a relationship between the global average temperature and a descriptor of nonequilibrium thermodynamics, the large scale temperature gradient between low and high latitudes. We find that a maximum of the temperature gradient is realized at the same value of the average temperature, about 270 K, largely independent of the strength of incoming solar radiation. We also examine the structural properties of the system defined by the equilibria depending on the solar strength. We construct new bifurcation diagrams in terms of quantities relevant for describing the thermodynamic disequilibrium, such as the temperature gradient and the material entropy production due to heat transport. We compare our results for the EBM to results for the intermediate complexity GCM PlaSim and find an interesting qualitative agreement.

preprint2014arXiv

A new framework for climate sensitivity and prediction: a modelling perspective

The sensitivity of climate models to increasing CO2 concentration and the climate response at decadal time scales are still major factors of uncertainty for the assessment of the long and short term effects of anthropogenic climate change. While the relative slow progress on these issues is partly due to the inherent inaccuracies of numerical climate models, this also hints at the need for stronger theoretical foundations to the problem of studying climate sensitivity and performing climate change predictions with numerical models. Here we demonstrate that it is possible to use Ruelle's response theory to predict the impact of an arbitrary CO2 forcing scenario on the global surface temperature of a general circulation model. Response theory puts the concept of climate sensitivity on firm theoretical grounds, and addresses rigorously the problem of predictability at different time scales. Conceptually, our results show that performing climate change experiments with general circulation models is a well defined problem from a physical and mathematical point of view. Practically, our results show that considering one single CO2 forcing scenario is enough to construct operators able to predict the response of climatic observables to any other CO2 forcing scenario, without the need to perform additional numerical simulations. We also introduce a general relationship between climate sensitivity and climate response at different time scales, thus providing an explicit definition of the inertia of the system at different time scales. While what we report here refers to the linear response, the general theory allows for treating nonlinear effects as well. Our results pave the way for redesigning and interpreting climate change experiments from a radically new perspective.

preprint2014arXiv

The Impact of Oceanic Heat Transport on the Atmospheric Circulation: a Thermodynamic Perspective

The present study investigates how global thermodynamic properties of the climate system are affected by the changes in the intensity of the imposed oceanic heat transport in an atmospheric general circulation model in aqua-planet configuration. Increasing the poleward oceanic heat transport results in an overall increase in the surface temperature and a decrease in the equator-to-pole surface temperature difference as a result of the ice-albedo feedback. Following the classical ansatz by Stone, the atmospheric heat transport changes in such a way that the total poleward heat transport remains almost unchanged. We also find that the efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport which suggests that the climate system becomes less efficient and turns into a state of reduced entropy production, as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperature, thus reducing the availability in the climate system and bringing it closer to a state of thermal equilibrium.

preprint2013arXiv

Towards a Fluctuation Theorem in an Atmospheric Circulation Model

An investigation of the distribution of finite time trajectory divergence is performed on an Atmospheric Global Circulation Model. The distribution of the largest local Lyapunov exponent shows a significant probability for negative values over time spans up to 10 days. This effect is present for resolutions up to wave numbers l=42 (~ 250km). The probability for a negative local largest Lyapunov exponent decreases over time, similarly to the predictions of the Fluctuation Theorem for entropy production. The model used is hydrostatic with variable numbers of vertical levels and different horizontal resolutions.

preprint2009arXiv

Thermodynamic Analysis of Snowball Earth Hysteresis Experiment: Efficiency, Entropy Production, and Irreversibility

We present an extensive thermodynamic analysis of a hysteresis experiment performed on a simplified yet Earth-like climate model. We slowly vary the solar constant by 20% around the present value and detect a substantial bistability: for a large range of values the realization of snowball (SB) or of warm (W) climate conditions depend on the history of the system. Using recent results on the global climate thermodynamics, we show that that the two regimes feature radically different properties. The efficiency of the climate machine increases with decreasing solar constant in W climate conditions, whereas the opposite takes place in the SB regime. Instead, entropy production is increasing with the solar constant in both branches of climate conditions, and its value is about 4 times as large in the W branch than in the corresponding C state. Finally, the degree of irreversibility of the system is much higher in the W conditions, with an explosive growth in the upper range of the considered values of solar constants. Whereas in the SB regime a dominating role is played by changes in the meridional albedo contrast, in the W climate regime changes in the intensity of latent heat fluxes are crucial for determining the observed properties. This substantiates the importance of addressing correctly the variations of the hydrological cycle in a changing climate. An interpretation of the transitions at the boundary of the bistable region based upon macro-scale thermodynamic properties is also proposed. Our results support the adoption of a new generation of diagnostic tools based on the 2nd law of thermodynamics for auditing climate model and outline a set of parameterizations to be used in conceptual and intermediate complexity models or for the reconstruction of the past climate conditions.