ALBERT

Description: During the four most recent glacial maxima, atmospheric CO2 has been lowered by about 90–100 ppm with respect to interglacial concentrations. It is likely that most of the atmospheric CO2 deficit was stored in the ocean. Changes in the biological pump, which are related to the efficiency of the biological carbon uptake in the surface ocean and/or of the export of organic carbon to the deep ocean, have been proposed as a key mechanism for the increased glacial oceanic CO2 storage. The biological pump is strongly constrained by the amount of available surface nutrients. In models, it is generally assumed that the ratio between elemental nutrients, such as phosphorus, and carbon (C∕P ratio) in organic material is fixed according to the classical Redfield ratio. The constant Redfield ratio appears to approximately hold when averaged over basin scales, but observations document highly variable C∕P ratios on regional scales and between species. If the C∕P ratio increases when phosphate availability is scarce, as observations suggest, this has the potential to further increase glacial oceanic CO2 storage in response to changes in surface nutrient distributions. In the present study, we perform a sensitivity study to test how a phosphate-concentration-dependent C∕P ratio influences the oceanic CO2 storage in an Earth system model of intermediate complexity (cGENIE). We carry out simulations of glacial-like changes in albedo, radiative forcing, wind-forced circulation, remineralization depth of organic matter, and mineral dust deposition. Specifically, we compare model versions with the classical constant Redfield ratio and an observationally motivated variable C∕P ratio, in which the carbon uptake increases with decreasing phosphate concentration. While a flexible C∕P ratio does not impact the model's ability to simulate benthic δ13C patterns seen in observational data, our results indicate that, in production of organic matter, flexible C∕P can further increase the oceanic storage of CO2 in glacial model simulations. Past and future changes in the C∕P ratio thus have implications for correctly projecting changes in oceanic carbon storage in glacial-to-interglacial transitions as well as in the present context of increasing atmospheric CO2 concentrations.

Description: The terrestrial biosphere is thought to be a key component in the climatic variability seen in the palaeo-record. It has a direct impact on surface temperature through changes in surface albedo and evapotranspiration (so-called biogeophysical effects) and, in addition, has an important indirect effect through changes in vegetation and soil carbon storage (biogeochemical effects) and hence modulates the concentrations of greenhouse gases in the atmosphere. The biogeochemical and biogeophysical effects generally have opposite signs, meaning that the terrestrial biosphere could potentially have played only a very minor role in the dynamics of the glacial–interglacial cycles of the late Quaternary. Here we use a fully coupled dynamic atmosphere–ocean–vegetation general circulation model (GCM) to generate a set of 62 equilibrium simulations spanning the last 120 kyr. The analysis of these simulations elucidates the relative importance of the biogeophysical versus biogeochemical terrestrial biosphere interactions with climate. We find that the biogeophysical effects of vegetation account for up to an additional −0.91 °C global mean cooling, with regional cooling as large as −5 °C, but with considerable variability across the glacial–interglacial cycle. By comparison, while opposite in sign, our model estimates of the biogeochemical impacts are substantially smaller in magnitude. Offline simulations show a maximum of +0.33 °C warming due to an increase of 25 ppm above our (pre-industrial) baseline atmospheric CO2 mixing ratio. In contrast to shorter (century) timescale projections of future terrestrial biosphere response where direct and indirect responses may at times cancel out, we find that the biogeophysical effects consistently and strongly dominate the biogeochemical effect over the inter-glacial cycle. On average across the period, the terrestrial biosphere has a −0.26 °C effect on temperature, with −0.58 °C at the Last Glacial Maximum. Depending on assumptions made about the destination of terrestrial carbon under ice sheets and where sea level has changed, the average terrestrial biosphere contribution over the last 120 kyr could be as much as −50 °C and −0.83 °C at the Last Glacial Maximum.

Description: Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high-resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled general circulation models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel long-term continuous projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations, which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted, which would not be possible with other modelling techniques on the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling surface air temperature – SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 kyr (thousand years) (by modelling SAT and precipitation). A range of CO2 scenarios are prescribed for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy sea surface temperature (SST) data from different sites and emulator results, finds that relatively similar CO2 concentrations are estimated based on sites at lower latitudes, whereas higher-latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ∼ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost. However, there are uncertainties associated with the approach used, including the inability of the emulator to capture deviations from a quasi-stationary response to the forcing, such as transient adjustments of the deep-ocean temperature and circulation, in addition to its limited range of fixed ice sheet configurations and its requirement for prescribed atmospheric CO2 concentrations.

Description: During the four most recent glacial cycles, atmospheric CO2 during glacial maxima has been lowered by about 90–100 ppm with respect to interglacials. There is widespread consensus that most of this carbon was partitioned in the ocean. It is however still debated which processes were dominant in achieving this increased carbon storage. In this paper, we use an Earth system model of intermediate complexity to constrain the range in ocean carbon storage for an ensemble of ocean circulation equilibrium states. We do a set of simulations where we run the model to pre-industrial equilibrium, but where we achieve different ocean circulation by changing forcing parameters such as wind stress, ocean diffusivity and atmospheric heat diffusivity. As a consequence, the ensemble members also have different ocean carbon reservoirs, global ocean average temperatures, biological pump efficiencies and conditions for air-sea CO2 disequilibrium. We analyse changes in total ocean carbon storage and separate it into contributions by the solubility pump, the biological pump and the CO2 disequilibrium component. We also relate these contributions to differences in strength of ocean overturning circulation. In cases with weaker circulation, we see that the ocean's capacity for carbon storage is larger. Depending on which ocean forcing parameter that is tuned, the origin of the change in carbon storage is different. When wind stress or ocean vertical diffusivity is changed, the response of the biological pump gives the most important effect on ocean carbon storage, whereas when atmospheric heat diffusivity or ocean horizontal diffusivity is changed, the solubility pump and the disequilibrium component are also important and sometimes dominant. Finally, we do a drawdown experiment, where we investigate the capacity for increased carbon storage by maximising the efficiency of the biological pump in our ensemble members. We conclude that different initial states for an ocean model result in different capacities for ocean carbon storage, due to differences in the ocean circulation state. This could explain why it is difficult to achieve comparable responses of the ocean carbon pumps in model intercomparison studies, where the initial states vary between models. The drawdown experiment highlights the importance of the strength of the biological pump in the control state for model studies of increased biological efficiency.

Description: The terrestrial biosphere is thought to be a key component in the climatic variability seen in the paleo record. It has a direct impact on surface temperature through changes in surface albedo and evapotranspiration (so called biogeophysical effects) and in addition, has an important indirect effect through changes in vegetation and soil carbon storage (biogeochemical effects) and hence modulates the concentrations of greenhouse gases in the atmosphere. The biogeochemical and biogeophysical effects generally have opposite signs meaning that the terrestrial biosphere could potentially have played only a very minor role in the dynamics of the glacial-interglacial cycles of the late Quaternary. Here we use a fully-coupled dynamic atmosphere-ocean-vegetation General Circulation Model (GCM) to generate a set of 62 simulations spanning the last 120 ka. The analysis of these simulations elucidates the relative importance of the biogeophysical versus biogeochemical terrestrial biosphere interactions with climate. We find that the biogeophysical effects of vegetation account for up to an additional −0.84 °C global mean cooling, with regional cooling as large as −5 °C, but with considerable variability across the glacial-interglacial cycle. By comparison, while opposite in sign, our model estimates of the biogeochemical impacts are substantially smaller in magnitude. Offline simulations show a maximum of +0.33 °C warming due an increase of 25 ppm above our (pre-industrial) baseline atmospheric CO2 mixing ratio. In contrast to shorter (century) time-scale projections of future terrestrial biosphere response where direct and indirect responses may at times, cancel out, we find that the biogeophysical effects consistently and strongly dominate the biogeochemical effect over the inter-glacial cycle. In addition, depending on the assumptions about soil carbon under ice-sheets and sea level rise, we find a range in terrestrial carbon storage change from a reduction in LGM carbon storage of −440 PgC, to a gain of +37 PgC. We suggest that prevailing uncertainties allow for only a small net transfer of carbon between terrestrial biosphere and ocean atmosphere implying that explaining the observed CO2 ice core record could be rather simpler than previously thought.

Description: Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled General Circulation Models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel long-term “continuous” projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted which would not be possible with other modelling techniques at the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 thousand years (kyr) (by modelling SAT and precipitation). A range of CO2 scenarios are modelled for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy data and emulator results and using proxy CO2 records, finds that relatively similar concentrations are produced at lower latitudes, although higher latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ~ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost.

Description: The response of the Earth System to greenhouse-gas driven warming is of critical importance for the future trajectory of our planetary environment. Hypethermal events – past climate transients with significant global-scale warming – can provide insights into the nature and magnitude of these responses. The largest hyperthermal of the Cenozoic was the Palaeocene-Eocene Thermal Maximum (PETM ~ 56 Ma). Here we present a new high-resolution cyclostratigraphy for the classic PETM section at Zumaia, Spain. With this new age model we are able to demonstrate that detrital sediment accumulation rates within this continental margin section increased more than four-fold during the PETM, representing a radical change in regional hydrology that drove dramatic increases in terrestrial to marine sediment flux. During the body of the PETM, orbital-scale variations in bulk sediment Si/Fe ratios are evidence for the continued orbital pacing of sediment erosion and transport processes, most likely linked to precession controls on sub-tropical hydroclimates. Most remarkable is that detrital accumulation rates remain high throughout the body of the PETM, and even reach peak values during the recovery phase of the characteristic PETM carbon isotope excursion (CIE). Using a series of Earth System Model inversions, we demonstrate that the silicate weathering feedback alone is insufficient to recover the PETM CIE, and that active organic carbon burial is required to match the observed dynamics of the CIE. Further, that the period of maximum organic carbon sequestration coincides with the peak in detrital accumulation rates observed at Zumaia. Based on these results, we hypothesize that precession controls on tropical and sub-tropical hydroclimates, and the sediment dynamics associated with this variation, play a significant role in the timing of the rapid climate and CIE recovery from peak-PETM conditions.

Description: The average depth in the ocean at which the majority of sinking organic matter particles remineralise is a fundamental parameter in the oceans role in regulating atmospheric CO2. Observed spatial patterns in sinking fluxes and relationships between the fluxes of different particles in the modern ocean have widely been used to invoke controlling mechanisms with important implications for CO2 regulation. However, such analyses are limited by the sparse spatial sampling of the available sediment trap data. Here we explore whether model ocean circulation rates, in the form of a transport matrix, can be used to derive remineralisation rates and sinking particle flux curves from the much more highly resolved observations of dissolved nutrient concentrations. Initially we use the Earth system model GENIE to generate a synthetic tracer dataset to explore the methods and its sensitivity to key sources of uncertainty arising from errors in the tracer observations and in the model circulation. We use a perturbed physics ensemble to generate 54 different estimates of circulation to explore errors associated with model transport rates. We find that reconstructed remineralisation rates are highly sensitive to both errors in observations and our ensemble uncertainty in model circulation rates such that a simple inversion does not provide a robust estimate of particulate flux profiles. Inferred remineralisation rates are particularly sensitive to differences between the "observed" and modelled transport fluxes because remineralisation rates are 3–4 magnitudes smaller than circulation rates. We also find that when inferring particle flux curves from remineralisation rates the cycling of dissolved organic matter also creates biases that have a similar magnitude and spatial variability to flux curves diagnosed using sediment trap data. We end with a discussion on the potential future directions and pitfalls of estimating remineralisation rates using model circulation schemes.

Description: A new model of terrestrial rock weathering – the Rock Geochemical Model (RokGeM) – was developed for incorporation into the GENIE Earth System modelling framework. In this paper we describe the model. We consider a range of previously devised parameterizations, ranging from simple dependencies on global mean temperature following Berner et al. (1983), to spatially explicit dependencies on run-off and temperature (GKWM, Bluth and Kump, 1994; GEM-CO2, Amiotte-Suchet et al., 2003) – fields provided by the energy-moisture balance atmosphere model component in GENIE. Using long-term carbon cycle perturbation experiments, we test the effects of a wide range of model parameters, including whether or not the atmosphere was "short-circuited" in the carbon cycle; the sensitivity and feedback strength of temperature and run-off on carbonate and silicate weathering; different river-routing schemes; 0-D (global average) vs. 2-D (spatially explicit) weathering schemes; and the lithology dependence of weathering. Included are details of how to run the model and visualize the results.

Description: Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.