ALBERT

Description: A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P − E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP∕dT) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Description: Antarctic ice volume has varied substantially during the late Quaternary, with reconstructions suggesting a glacial ice sheet extending to the continental shelf break and interglacial sea level highstands of several meters. Throughout this period, changes in the Antarctic Ice Sheet were driven by changes in atmospheric and oceanic conditions and global sea level; yet, so far modeling studies have not addressed which of these environmental forcings dominate and how they interact in the dynamical ice sheet response. Here, we force an Antarctic Ice Sheet model with global sea level reconstructions and transient, spatially explicit boundary conditions from a 408 ka climate model simulation, not only in concert with each other but, for the first time, also separately. We find that together these forcings drive glacial–interglacial ice volume changes of 12–14 ms.l.e., in line with reconstructions and previous modeling studies. None of the individual drivers – atmospheric temperature and precipitation, ocean temperatures, or sea level – single-handedly explains the full ice sheet response. In fact, the sum of the individual ice volume changes amounts to less than half of the full ice volume response, indicating the existence of strong nonlinearities and forcing synergy. Both sea level and atmospheric forcing are necessary to create full glacial ice sheet growth, whereas the contribution of ocean melt changes is found to be more a function of ice sheet geometry than climatic change. Our results highlight the importance of accurately representing the relative timing of forcings of past ice sheet simulations and underscore the need for developing coupled climate–ice sheet modeling frameworks that properly capture key feedbacks.

Description: Despite intense efforts, the mechanisms that drive glacial–interglacial changes in atmospheric pCO2 are not fully understood. Here, we aim at quantifying the potential contribution of aeolian dust deposition changes to the atmospheric pCO2 drawdown during the Last Glacial Maximum (LGM). To this end, we use the Max Planck Institute Ocean Model (MPIOM) and the embedded Hamburg Ocean Carbon Cycle model (HAMOCC), including a new parameterization of particle ballasting that accounts for the acceleration of sinking organic soft tissue in the ocean by higher-density biogenic calcite and opal particles, as well as mineral dust. Sensitivity experiments with reconstructed LGM dust deposition rates indicate that the acceleration of detritus by mineral dust played a small role in atmospheric pCO2 variations during glacial–interglacial cycles – on the order of 5 ppmv, compared to the reconstructed ∼80 ppmv rise in atmospheric pCO2 during the last deglaciation. The additional effect of the LGM dust deposition, namely the enhanced fertilization by the iron that is associated with the glacial dust, likely played a more important role; although the full iron fertilization effect can not be estimated in the particular model version used here due to underestimated present-day non-diazotroph iron limitation, fertilization of diazotrophs in the tropical Pacific already leads to an atmospheric pCO2 drawdown of around 10 ppmv.

Description: Despite intense efforts, the mechanisms that drive glacial–interglacial changes in atmospheric pCO2 are not fully understood. Here, we aim at quantifying the potential contribution of aeolian dust deposition changes to the atmospheric pCO2 drawdown during the Last Glacial Maximum (LGM). To this end, we use the ocean circulation and carbon cycle model MPIOM/HAMOCC, including a new parameterisation of particle ballasting that accounts for the acceleration of sinking organic soft tissue in the ocean by higher density biogenic calcite and opal particles, as well as mineral dust. Sensitivity experiments with reconstructed LGM dust deposition rates indicate that the acceleration of detritus by mineral dust likely played a small role for atmospheric pCO2 variations during glacial–interglacial cycles – on the order of 5ppmv, compared to the reconstructed ∼80ppmv-rise of atmospheric pCO2 during the last deglaciation. The additional effect of the LGM dust deposition, namely the enhanced fertilisation by the iron that is associated with the glacial dust, played a more important role – leading to a pCO2-drawdown by more than 8ppmv in our LGM sensitivity experiments despite an underestimated iron-limitation in the used model setup.

Description: Antarctic ice volume has varied substantially during the Quaternary, with reconstructions suggesting a glacial ice sheet extending to the continental shelf break, and interglacial sea level highstands of several meters. Throughout this period, changes in the Antarctic ice sheet were driven by changes in atmospheric and oceanic conditions and global sea level, yet so far, modeling studies have not addressed which of these environmental forcings dominate, and how they interact in the dynamical ice sheet response. Here we force an Antarctic ice sheet model with global sea level reconstructions and transient, spatially explicit boundary conditions from a 408 ka climate model simulation, not only in concert with each other but, for the first time, also separately. We find that together, these forcings drive glacial-interglacial ice volume changes of 12–14 m SLE, in line with reconstructions and previous modeling studies. None of the individual drivers – atmospheric temperature and precipitation, ocean temperatures, sea level – single-handedly explains the full ice sheet response. In fact, the sum of the individual ice volume changes amounts to less than half of the full ice volume response, indicating the existence of strong nonlinearities and forcing synergy. Both sea level and atmospheric forcing are necessary to create full glacial ice sheet growth, whereas the contribution of ocean melt changes is found to be more a function of ice sheet geometry than climatic change. Our results highlight the importance of accurately representing the relative timing of forcings of past ice sheet simulations, and underscore the need for developing coupled climate-ice sheet modeling frameworks that properly capture key feedbacks.

Description: It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice sheets have not been fully understood. Testing theories of ice ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240 to 170 ka (ka: thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate the transient evolution of Northern Hemisphere and Southern Hemisphere ice sheets during the MIS 7–6 period in response to orbital and greenhouse gas forcing. For a range of model parameters, the simulations capture the evolution of global ice volume well within the range of reconstructions. Over the MIS 7–6 period, it is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary-wave–ice-sheet topography feedback can trigger an unabated and unrealistic ice sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate–ice sheet modeling, with such coupled paleo-simulations providing opportunities to constrain such parameters.

Description: A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Description: Reduction of anthropogenic CO2 emissions alone will not sufficiently restrict global warming and enable the 1.5°C goal of the Paris agreement to be met. To effectively counteract climate change, measures to actively remove carbon dioxide from the atmosphere are required. Artificial upwelling has been proposed as one such carbon dioxide removal technique. By fueling primary productivity in the surface ocean with nutrient-rich deep water, it could potentially enhance downward fluxes of particulate organic carbon (POC) and carbon sequestration. In this study we investigated the effect of different intensities of artificial upwelling combined with two upwelling modes (recurring additions vs. one singular addition) on POC export, sinking matter stoichiometry and remineralization depth. We carried out a 39 day-long mesocosm experiment in the subtropical North Atlantic, where we fertilized oligotrophic surface waters with different amounts of deep water. The total nutrient inputs ranged from 1.6 to 11.0 μmol NO3– L–1. We found that on the one hand POC export under artificial upwelling more than doubled, and the molar C:N ratios of sinking organic matter increased from values around Redfield (6.6) to ∼8–13, which is beneficial for potential carbon dioxide removal. On the other hand, sinking matter was remineralized at faster rates and showed lower sinking velocities, which led to shallower remineralization depths. Particle properties were more favorable for deep carbon export in the recurring upwelling mode, while in the singular mode the C:N increase of sinking matter was more pronounced. In both upwelling modes roughly half of the produced organic carbon was retained in the water column until the end of the experiment. This suggests that the plankton communities were still in the process of adjustment, possibly due to the different response times of producers and consumers. There is thus a need for studies with longer experimental durations to quantify the responses of fully adjusted communities. Finally, our results revealed that artificial upwelling affects a variety of sinking particle properties, and that the intensity and mode with which it is applied control the strength of the effects.

Description: It is widely accepted that orbital variations are responsible for the generation of glacial cycles during the late Pleistocene. However, the relative contributions of the orbital forcing compared to CO2 variations and other feedback mechanisms causing the waxing and waning of ice sheets have not been fully understood. Testing theories of ice ages beyond statistical inferences, requires numerical modeling experiments that capture key features of glacial transitions. Here, we focus on the glacial buildup from Marine Isotope Stage (MIS) 7 to 6 covering the period from 240 to 170 ka (ka: thousand years before present). This transition from interglacial to glacial conditions includes one of the fastest Pleistocene glaciation–deglaciation events, which occurred during MIS 7e–7d–7c (236–218 ka). Using a newly developed three-dimensional coupled atmosphere–ocean–vegetation–ice sheet model (LOVECLIP), we simulate the transient evolution of Northern Hemisphere and Southern Hemisphere ice sheets during the MIS 7–6 period in response to orbital and greenhouse gas forcing. For a range of model parameters, the simulations capture the evolution of global ice volume well within the range of reconstructions. Over the MIS 7–6 period, it is demonstrated that glacial inceptions are more sensitive to orbital variations, whereas terminations from deep glacial conditions need both orbital and greenhouse gas forcings to work in unison. For some parameter values, the coupled model also exhibits a critical North American ice sheet configuration, beyond which a stationary-wave–ice-sheet topography feedback can trigger an unabated and unrealistic ice sheet growth. The strong parameter sensitivity found in this study originates from the fact that delicate mass imbalances, as well as errors, are integrated during a transient simulation for thousands of years. This poses a general challenge for transient coupled climate–ice sheet modeling, with such coupled paleo-simulations providing opportunities to constrain such parameters.