Bracco, 1049095 Wang, 1049114 Zaliapin, 1049092
The project team is made up of climate dynamicists and of applied mathematicians. The investigators (a) formulate a mathematical theory of climate sensitivity and (b) devise a set of optimization algorithms for general circulation models and Earth System Models. The team brings together strengths in dynamical systems, partial differential equations, and numerical methods, with depth and broad coverage in the study of atmospheric, oceanic, and climate dynamics. The project's three main objectives are to: (i) continue developing powerful new methods for the fundamental understanding of climate sensitivity and predictability; (ii) extend earlier work of the investigators on modes of low-frequency variability associated with the El Nino-Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO), interannual as well as decadal; and (iii) combine items (i) and (ii) in analyzing the sensitivity and predictability of these modes when subjected to climate change. All three objectives are pursued across a full hierarchy of models, from conceptual "toy" models through intermediate climate models and on to Earth system Models of Intermediate Complexity.
Ghil and his associates have recently worked on extending the theory of random dynamical systems and applying it to the climate system. This theory allows one to (1) investigate the effect of random perturbations ("weather") on nonlinear dynamical systems ("climate variability"); (2) evaluate the robustness and sensitivity of a random dynamical system to changes in either the system or its forcing, whether deterministic (e.g., slow, anthropogenic changes in greenhouse gas or aerosol concentrations) or stochastic (e.g., volcanic eruptions); and (3) obtain sharper results on the system's predictability by accounting for the effect of the random perturbations. Methods developed for the systematic study of parameter dependence in a streamlined global circulation model, the ICTP-AGCM, have promising parallels to results published by the PIs and co-workers on idealized models. The team obtains rigorous results on the latter kinds of models, as well as on random dynamical system bifurcations, sensitivity, and predictability, while extending the ICTP-AGCM results to models of intermediate complexity like SPEEDO, and eventually to full Earth System Models like the Community Climate System Model (CCSM). This work leads to a deeper understanding of the causes and mechanisms of climate sensitivity. It also provides efficient ways to evaluate and improve the ability of global circulation models and Earth System Models to simulate past and present climate, and to predict our environment's future evolution. It helps strengthen the basis for robust climate projections on decade-to-century time scales, and it provides a systematic way to evaluate and improve both deterministic and stochastic parameterizations in such models. The results of this work have implications for other areas in which complex deterministic dynamics interacts with external forcing, deterministic as well as random. This situation characterizes the life and socio-economic sciences, as well as climate science and the geosciences. Strong interactions across disciplinary boundaries -- among team members themselves and with colleagues in other areas -- help accelerate the transfer of new methods and results to other disciplines.
The climate system comprises the atmosphere, oceans, ice masses and the biosphere. Its intrinsic complexity leads to chaotic and random behavior on many scales of time and space. In addition, it is an open system that exchanges heat and momentum with the Sun and the interplanetary medium, as well as being affected by human activities, like greenhouse gas and aerosol generation, that change in time. The complex dynamics of this open system requires novel tools for its analysis. In this project, we have developed such novel tools and brought them to bear on the interactions between the climate system’s intrinsic dynamics and its time-? dependent forcing, both natural and anthropogenic. Random dynamical systems and their attractors. The proper tools for studying this interaction arise from the theory of nonlinear, deterministically and stochastically forced dynamical systems, called random dynamical systems (RDS) for short. We have shown that this interaction can take novel, spectacular forms even in highly simplified systems, such as the stochastically forced Lorenz system. The original form of this system does not include any time-dependent forcing and it still exhibits the by now well known "butterfly effect," i.e., sensitive dependence on its initial state. In the presence of random forcing, the original system’s strange attractor, with its two butterfly wings, becomes a time-dependent object. Four snapshots of the evolution in time of this object are shown in Figure 1. The figure illustrates the nature of the interactions between the nonlinear deterministic dynamics, with the two wings still visible, and of the random forcing, with jitter around the wings. Smooth and rough parameter dependence in climate models. Atmospheric, oceanic and coupled climate models are based on the numerically discretized version of the governing equations of large-?scale planetary flows. Many details, though, like the evolution of clouds and their interaction with the radiation and the large-scale flow, cannot be derived from first principles. The net effects of such subgrid-scale phenomena are modeled on a semi-empirical basis and the so-called parameterizations employed in this modeling process involve many uncertainties in parameter values. Hence it is of paramount importance to evaluate the sensitivity of model results to these uncertainties. The dependence of certain model statistics on parameter values can be either fairly smooth or, to the contrary, quite sudden or rough. We studied this behavior in an El Niño–Southern Oscillation (ENSO) model of intermediate complexity and showed that the sensitive parameter dependence of the statistics of model observables is related to the dominant ENSO modes of variability. A crucial result is that the model statistics are more sensitive for regimes in which the climate system’s low-frequency variability (LFV) is more pronounced, and autocorrelations decay more slowly; LFV here refers to the presence of periodic or near-periodic processes within the system, such as the seasonal cycle or ENSO’s quasi-biennial mode; see Figure 2. Systematic parameter optimization in climate models. Global climate models (GCMs) are an increasingly important source of information for policymakers. GCMs are the "upper echelons" of the climate modeling hierarchy, and the systematic and automated optimization of the semi-?empirical parameters mentioned above is of paramount importance for their reliability. Given the tenths of such parameters in an IPCC-class GCM, and the computational cost of running such GCMs for any given parameter setting, it is necessary to develop meta-models — also called emulators or response surfaces — that replicate the results of the GCMs at much lower cost. Meta-modeling had previously been applied by members of our team to the sensitivity of the atmospheric GCM of the International Centre for Theoretical Physics (ICTP) in Trieste, Italy. Their meta-model accurately estimated global spatial averages of common fields of climatic interest, including precipitation, low- and high-level winds, temperature at various levels, sea level pressure and geopotential height, while providing a computationally cheap strategy to explore the influence of a fairly wide range of parameter settings. In the recently completed project we are reporting on, these results have been further expanded to high-dimensional spatial fields, in order to examine the considerable regional-scale biases to which current GCMs are still subject. These biases may vary substantially depending on the climate variable considered, as well as on the performance metric adopted. Common dilemmas in parameter optimization for GCMs are associated with model revisions that yield improvement in one field or regional pattern or season, but degradation in another, or improvement in the model climatology but degradation in its ability to reproduce interannual variability. Challenges are posed by the high dimensionality of a GCM’s output fields and by the large number of adjustable parameters. We have shown that the use of the meta- model in the optimization strategy helps visualize trade-offs at a regional level, e.g., how mismatches between sensitivity and spatial error fields yield regional errors under minimization of global objective functions; see Figure 3.