Michel Crucifix

Thermohaline circulation hysteresis: a model intercomparison

We present results from an intercomparison of 11 different climate models of intermediate complexity, in which the North Atlantic Ocean was subjected to slowly varying changes in freshwater input. All models show a characteristic hysteresis response of the thermohaline circulation to the freshwater forcing; which can be explained by Stommels salt advection feedback. The width of the hysteresis curves varies between 0.2 and 0.5 Sv in the models. Major differences are found in the location of present-day climate on the hysteresis diagram. In seven of the models, present-day climate for standard parameter choices is found in the bi-stable regime, in four models this climate is in the mono-stable regime. The proximity of the present-day climate to the Stommel bifurcation point, beyond which North Atlantic Deep Water formation cannot be sustained, varies from less than 0.1 Sv to over 0.5 Sv.

Last glacial maximum ocean thermohaline circulation: PMIP2 model intercomparisons and data constraints

The ocean thermohaline circulation is important for transports of heat and the carbon cycle. We present results from PMIP2 coupled atmosphere-ocean simulations with four climate models that are also being used for future assessments. These models give very different glacial thermohaline circulations even with comparable circulations for present. An integrated approach using results from these simulations for Last Glacial Maximum (LGM) with proxies of the state of the glacial surface and deep Atlantic supports the interpretation from nutrient tracers that the boundary between North Atlantic Deep Water and Antarctic Bottom Water was much shallower during this period. There is less constraint from this integrated reconstruction regarding the strength of the LGM North Atlantic overturning circulation, although together they suggest that it was neither appreciably stronger nor weaker than modern. Two model simulations identify a role for sea ice in both hemispheres in driving the ocean response to glacial forcing.

Estimating shortwave radiative forcing and response in climate models

Abstract Feedback analysis in climate models commonly involves decomposing any change in the system’s energy balance into radiative forcing terms due to prescribed changes, and response terms due to the radiative effects of changes in model variables such as temperature, water vapor, clouds, sea ice, and snow. The established “partial radiative perturbation” (PRP) method allows an accurate separation of these terms, but requires processing large volumes of model output with an offline version of the model’s radiation code. Here, we propose an “approximate PRP” (APRP) method for the shortwave that provides an accurate estimate of the radiative perturbation, but derived from a quite modest amount of monthly mean model output. The APRP method is based on a simplified shortwave radiative model of the atmosphere, where surface absorption and atmospheric scattering and absorption are represented by means of three parameters that are diagnosed for overcast and clear-sky portions of each model grid cell. The accu...

Climate of the last millennium: a sensitivity study

Seventy-one sensitivity experiments have been performed using a two-dimensional sectoraveragedglobal climate model to assess the potential impact of six different factors on the lastmillennium climate and in particular on the surface air temperature evolution. Both natural(i.e, solar and volcanism) and anthropogenically-induced (i.e. deforestation, additional greenhousegases, and tropospheric aerosol burden) climate forcings have been considered. Comparisons of climate reconstructions with model results indicate that all the investigated forcings are needed to simulate the surface air temperature evolution. Due to uncertainties inhistorical climate forcings and temperature reconstructions, the relative importance of a particularforcing in the explanation of the recorded temperature variance is largely function of theforcing time series used. Nevertheless, our results indicate that whatever the historical solar andvolcanic reconstructions may be, these externally driven natural climate forcings are unable togive climate responses comparable in magnitude and time to the late—20th-century temperaturewarming while for earlier periods combination of solar and volcanic forcings can explain theLittle Ice Age and the Medieval Warm Period. Only the greenhouse gas forcing allows themodel to simulate an accelerated warming rate during the last three decades. The best guesssimulation (largest similarity with the reconstruction) for the period starting 1850 AD requireshowever to include anthropogenic sulphate forcing as well as the impact of deforestation toconstrain the magnitude of the greenhouse gas twentieth century warming to better fit theobservation. On the contrary, prior to 1850 AD mid-latitude land clearance tends to reinforcethe Little Ice age in our simulations.

STABILITY ANALYSIS OF THE CLIMATE-VEGETATION SYSTEM IN THE NORTHERN HIGH LATITUDES

The stability of the climate-vegetation system in the northern high latitudesis analysed with three climate system models of different complexity: A comprehensive 3-dimensional model of the climate system, GENESIS-IBIS, and two Earth system models of intermediate complexity (EMICs), CLIMBER-2 andMoBidiC. The biogeophysical feedback in the latitudinal belt 60–70° N, although positive, is not strong enough to support multiple steady states: A unique equilibriumin the climate-vegetation system is simulated by all the models on a zonal scale for present-day climate and doubled CO2 climate.EMIC simulations with decreased insolation also reveal a unique steady state. However, the climate sensitivity to tree cover,Δ TF, exhibits non-linear behaviour within the models. For GENESIS-IBIS and CLIMBER-2, Δ TF islower for doubled CO2 climate than for present-day climate due to a shorter snow season and increased relative significance ofthe hydrological effect of forest cover. For the EMICs, Δ TF is higher for low tree fraction than for high treefraction, mainly due to a time shift in spring snow melt in response to changes in tree cover. The climate sensitivity to tree coveris reduced when thermohaline circulation feedbacks are accounted for in the EMIC simulations. Simpler parameterizations of oceanic processes have opposite effects onΔ TF: Δ TF is lower in simulations with fixed SSTs and higher in simulations with mixed layer oceans. Experiments with transient CO2 forcing show climate and vegetation not in equilibrium in the northern high latitudes at the end of the 20thcentury. The delayed response of vegetation and accelerated global warming lead to rather abrupt changes in northern vegetation cover in the first halfof the 21st century, when vegetation cover changes at double the present day rate.

Interglacials of the last 800,000 years

Interglacials, including the present (Holocene) period, are warm, low land ice extent (high sea level), end-members of glacial cycles. Based on a sea level definition, we identify eleven interglacials in the last 800,000 years, a result that is robust to alternative definitions. Data compilations suggest that despite spatial heterogeneity, Marine Isotope Stages (MIS) 5e (last interglacial) and 11c (~400 ka ago) were globally strong (warm), while MIS 13a (~500 ka ago) was cool at many locations. A step change in strength of interglacials at 450 ka is apparent only in atmospheric CO2 and in Antarctic and deep ocean temperature. The onset of an interglacial (glacial termination) seems to require a reducing precession parameter (increasing Northern Hemisphere summer insolation), but this condition alone is insufficient. Terminations involve rapid, nonlinear, reactions of ice volume, CO2, and temperature to external astronomical forcing. The precise timing of events may be modulated by millennial-scale climate change that can lead to a contrasting timing of maximum interglacial intensity in each hemisphere. A variety of temporal trends is observed, such that maxima in the main records are observed either early or late in different interglacials. The end of an interglacial (glacial inception) is a slower process involving a global sequence of changes. Interglacials have been typically 10–30 ka long. The combination of minimal reduction in northern summer insolation over the next few orbital cycles, owing to low eccentricity, and high atmospheric greenhouse gas concentrations implies that the next glacial inception is many tens of millennia in the future.

Enhanced ice sheet growth in Eurasia owing to adjacent ice-dammed lakes

Large proglacial lakes cool regional summer climate because of their large heat capacity, and have been shown to modify precipitation through mesoscale atmospheric feedbacks, as in the case of Lake Agassiz. Several large ice-dammed lakes, with a combined area twice that of the Caspian Sea, were formed in northern Eurasia about 90,000 years ago, during the last glacial period when an ice sheet centred over the Barents and Kara seas blocked the large northbound Russian rivers. Here we present high-resolution simulations with an atmospheric general circulation model that explicitly simulates the surface mass balance of the ice sheet. We show that the main influence of the Eurasian proglacial lakes was a significant reduction of ice sheet melting at the southern margin of the Barents–Kara ice sheet through strong regional summer cooling over large parts of Russia. In our simulations, the summer melt reduction clearly outweighs lake-induced decreases in moisture and hence snowfall, such as has been reported earlier for Lake Agassiz. We conclude that the summer cooling mechanism from proglacial lakes accelerated ice sheet growth and delayed ice sheet decay in Eurasia and probably also in North America.

Using the past to constrain the future: how the palaeorecord can improve estimates of global warming

Climate sensitivity is defined as the change in global mean equilibrium temperature after a doubling of atmospheric CO2 concentration and provides a simple measure of global warming. An early estimate of climate sensitivity, 1.5—4.5°C, has changed little subsequently, including the latest assessment by the Intergovernmental Panel on Climate Change. The persistence of such large uncertainties in this simple measure casts doubt on our understanding of the mechanisms of climate change and our ability to predict the response of the climate system to future perturbations. This has motivated continued attempts to constrain the range with climate data, alone or in conjunction with models. The majority of studies use data from the instrumental period (post-1850), but recent work has made use of information about the large climate changes experienced in the geological past. In this review, we first outline approaches that estimate climate sensitivity using instrumental climate observations and then summarize attempts to use the record of climate change on geological timescales. We examine the limitations of these studies and suggest ways in which the power of the palaeoclimate record could be better used to reduce uncertainties in our predictions of climate sensitivity.

Is the astronomical forcing a reliable and unique pacemaker for climate? A conceptual model study

There is evidence that ice age cycles are paced by astronomical forcing, suggesting some kind of synchronisation phenomenon. Here, we identify the type of such synchronisation and explore systematically its uniqueness and robustness using a simple paleoclimate model akin to the van der Pol relaxation oscillator and dynamical system theory. As the insolation is quite a complex quasiperiodic signal involving different frequencies, the traditional concepts used to define synchronisation to periodic forcing are no longer applicable. Instead, we explore a different concept of generalised synchronisation in terms of (coexisting) synchronised solutions for the forced system, their basins of attraction and instabilities. We propose a clustering technique to compute the number of synchronised solutions, each of which corresponds to a different paleoclimate history. In this way, we uncover multistable synchronisation (reminiscent of phase- or frequency-locking to individual periodic components of astronomical forcing) at low forcing strength, and monostable or unique synchronisation at stronger forcing. In the multistable regime, different initial conditions may lead to different paleoclimate histories. To study their robustness, we analyse Lyapunov exponents that quantify the rate of convergence towards each synchronised solution (local stability), and basins of attraction that indicate critical levels of external perturbations (global stability). We find that even though synchronised solutions are stable on a long term, there exist short episodes of desynchronisation where nearby climate trajectories diverge temporarily (for about 50 kyr). As the attracting trajectory can sometimes lie close to the boundary of its basin of attraction, a small perturbation could quite easily make climate to jump between different histories, reducing the predictability. Our study brings new insight into paleoclimate dynamics and reveals a possibility for the climate system to wander throughout different climatic histories related to preferential synchronisation regimes on obliquity, precession or combinations of both, all over the history of the Pleistocene.

TESTING A PARTICLE FILTER TO RECONSTRUCT CLIMATE CHANGES OVER THE PAST CENTURIES

We implement a data-assimilation method based on a particle filter in the coupled climate model LEVECLIM focusing on decadal to centennial time scales. Several tests are performed with particle filtering using pseudo-observations obtained from a twin experiment with the model, as well as using real-data observations over the last century. These tests demonstrate that it is possible to obtain a model output well correlated with the observations at the large scale at a reasonable cost