Search Results

Now showing 1 - 10 of 14
  • Item
    On the effect of orbital forcing on mid-Pliocene climate, vegetation and ice sheets
    (München : European Geopyhsical Union, 2013) Willeit, M.; Ganopolski, A.; Feulner, G.
    We present results from modelling of the mid-Pliocene warm period (3.3–3 million years ago) using the Earth system model of intermediate complexity CLIMBER-2 analysing the effect of changes in boundary conditions as well as of orbital forcing on climate. First we performed equilibrium experiments following the PlioMIP (Pliocene Model Intercomparison Project) protocol with a CO2 concentration of 405 ppm, reconstructed mid-Pliocene orography and vegetation and a present-day orbital configuration. Simulated global Pliocene warming is about 2.5 °C, fully consistent with results of atmosphere–ocean general circulation model simulations performed for the same modelling setup. A factor separation analysis attributes 1.5 °C warming to CO2, 0.3 °C to orography, 0.2 °C to ice sheets and 0.4 °C to vegetation. Transient simulations for the entire mid-Pliocene warm period with time-dependent orbital forcing as well as interactive ice sheets and vegetation give a global warming varying within the range 1.9–2.8 °C. Ice sheet and vegetation feedbacks in synergy act as amplifiers of the orbital forcing, transforming seasonal insolation variations into an annual mean temperature signal. The effect of orbital forcing is more significant at high latitudes, especially during boreal summer, when the warming over land varies in the wide range from 0 to 10 °C. The modelled ice-sheet extent and vegetation distribution also show significant temporal variations. Modelled and reconstructed data for Northern Hemisphere sea-surface temperatures and vegetation distribution show the best agreement if the reconstructions are assumed to be representative for the warmest periods during the orbital cycles. This suggests that low-resolution Pliocene palaeoclimate reconstructions can reflect not only the impact of increased CO2 concentrations and topography changes but also the effect of orbital forcing. Therefore, the climate (Earth system) sensitivity estimates from Pliocene reconstructions which do not account for the effect of orbital forcing can be biased toward high values.
  • Item
    The importance of snow albedo for ice sheet evolution over the last glacial cycle
    (München : European Geopyhsical Union, 2018) Willeit, Matteo; Ganopolski, Andrey
    The surface energy and mass balance of ice sheets strongly depends on the amount of solar radiation absorbed at the surface, which is mainly controlled by the albedo of snow and ice. Here, using an Earth system model of intermediate complexity, we explore the role played by surface albedo for the simulation of glacial cycles. We show that the evolution of the Northern Hemisphere ice sheets over the last glacial cycle is very sensitive to the representation of snow albedo in the model. It is well known that the albedo of snow depends strongly on snow grain size and the content of light-absorbing impurities. Excluding either the snow aging effect or the dust darkening effect on snow albedo leads to an excessive ice build-up during glacial times and consequently to a failure in simulating deglaciation. While the effect of snow grain growth on snow albedo is well constrained, the albedo reduction due to the presence of dust in snow is much more uncertain because the light-absorbing properties of dust vary widely as a function of dust mineral composition. We also show that assuming slightly different optical properties of dust leads to very different ice sheet and climate evolutions in the model. Conversely, ice sheet evolution is less sensitive to the choice of ice albedo in the model. We conclude that a proper representation of snow albedo is a fundamental prerequisite for a successful simulation of glacial cycles.
  • Item
    A simple stress-based cliff-calving law
    (Göttingen : Copernicus GmbH, 2019) Schlemm, T.; Levermann, A.
    Over large coastal regions in Greenland and Antarctica the ice sheet calves directly into the ocean. In contrast to ice-shelf calving, an increase in calving from grounded glaciers contributes directly to sea-level rise. Ice cliffs with a glacier freeboard larger than ≈100 m are currently not observed, but it has been shown that such ice cliffs are increasingly unstable with increasing ice thickness. This cliff calving can constitute a self-amplifying ice loss mechanism that may significantly alter sea-level projections both of Greenland and Antarctica. Here we seek to derive a minimalist stress-based parametrization for cliff calving from grounded glaciers whose freeboards exceed the 100 m stability limit derived in previous studies. This will be an extension of existing calving laws for tidewater glaciers to higher ice cliffs.

    To this end we compute the stress field for a glacier with a simplified two-dimensional geometry from the two-dimensional Stokes equation. First we assume a constant yield stress to derive the failure region at the glacier front from the stress field within the glacier. Secondly, we assume a constant response time of ice failure due to exceedance of the yield stress. With this strongly constraining but very simple set of assumptions we propose a cliff-calving law where the calving rate follows a power-law dependence on the freeboard of the ice with exponents between 2 and 3, depending on the relative water depth at the calving front. The critical freeboard below which the ice front is stable decreases with increasing relative water depth of the calving front. For a dry water front it is, for example, 75 m. The purpose of this study is not to provide a comprehensive calving law but to derive a particularly simple equation with a transparent and minimalist set of assumptions.

  • Item
    Scaling of instability timescales of Antarctic outlet glaciers based on one-dimensional similitude analysis
    (Göttingen : Copernicus GmbH, 2019) Levermann, A.; Feldmann, J.
    Recent observations and ice-dynamic modeling suggest that a marine ice-sheet instability (MISI) might have been triggered in West Antarctica. The corresponding outlet glaciers, Pine Island Glacier (PIG) and Thwaites Glacier (TG), showed significant retreat during at least the last 2 decades. While other regions in Antarctica have the topographic predisposition for the same kind of instability, it is so far unclear how fast these instabilities would unfold if they were initiated. Here we employ the concept of similitude to estimate the characteristic timescales of several potentially MISI-prone outlet glaciers around the Antarctic coast. Our results suggest that TG and PIG have the fastest response time of all investigated outlets, with TG responding about 1.25 to 2 times as fast as PIG, while other outlets around Antarctica would be up to 10 times slower if destabilized. These results have to be viewed in light of the strong assumptions made in their derivation. These include the absence of ice-shelf buttressing, the one-dimensionality of the approach and the uncertainty of the available data. We argue however that the current topographic situation and the physical conditions of the MISI-prone outlet glaciers carry the information of their respective timescale and that this information can be partially extracted through a similitude analysis.
  • Item
    Modeling Antarctic tides in response to ice shelf thinning and retreat
    (Hoboken, NJ : Blackwell Publishing Ltd, 2014) Rosier, S.H.R.; Green, J.A.M.; Scourse, J.D.; Winkelmann, R.
    Tides play an important role in ice sheet dynamics by modulating ice stream velocity, fracturing, and moving ice shelves and mixing water beneath them. Any changes in ice shelf extent or thickness will alter the tidal dynamics through modification of water column thickness and coastal topography but these will in turn feed back onto the overall ice shelf stability. Here, we show that removal or reduction in extent and/or thickness of the Ross and Ronne-Filchner ice shelves would have a significant impact on the tides around Antarctica. The Ronne-Filchner appears particularly vulnerable, with an increase in M2 amplitude of over 0.5 m beneath much of the ice shelf potentially leading to tidally induced feedbacks on ice shelf/sheet dynamics. These results highlight the importance of understanding tidal feedbacks on ice shelves/streams due to their influence on ice sheet dynamics.
  • Item
    Committed sea-level rise under the Paris Agreement and the legacy of delayed mitigation action
    (London : Nature Publishing Group, 2018) Mengel, M.; Nauels, A.; Rogelj, J.; Schleussner, C.-F.
    Sea-level rise is a major consequence of climate change that will continue long after emissions of greenhouse gases have stopped. The 2015 Paris Agreement aims at reducing climate-related risks by reducing greenhouse gas emissions to net zero and limiting global-mean temperature increase. Here we quantify the effect of these constraints on global sea-level rise until 2300, including Antarctic ice-sheet instabilities. We estimate median sea-level rise between 0.7 and 1.2 m, if net-zero greenhouse gas emissions are sustained until 2300, varying with the pathway of emissions during this century. Temperature stabilization below 2 °C is insufficient to hold median sea-level rise until 2300 below 1.5 m. We find that each 5-year delay in near-term peaking of CO2 emissions increases median year 2300 sea-level rise estimates by ca. 0.2 m, and extreme sea-level rise estimates at the 95th percentile by up to 1 m. Our results underline the importance of near-term mitigation action for limiting long-term sea-level rise risks.
  • Item
    MIS-11 duration key to disappearance of the Greenland ice sheet
    (London : Nature Publishing Group, 2017) Robinson, A.; Alvarez-Solas, J.; Calov, R.; Ganopolski, A.; Montoya, M.
    Palaeo data suggest that Greenland must have been largely ice free during Marine Isotope Stage 11 (MIS-11). However, regional summer insolation anomalies were modest during this time compared to MIS-5e, when the Greenland ice sheet likely lost less volume. Thus it remains unclear how such conditions led to an almost complete disappearance of the ice sheet. Here we use transient climate-ice sheet simulations to simultaneously constrain estimates of regional temperature anomalies and Greenland's contribution to the MIS-11 sea-level highstand. We find that Greenland contributed 6.1 m (3.9-7.0 m, 95% credible interval) to sea level, ∼7 kyr after the peak in regional summer temperature anomalies of 2.8 °C (2.1-3.4 °C). The moderate warming produced a mean rate of mass loss in sea-level equivalent of only around 0.4 m per kyr, which means the long duration of MIS-11 interglacial conditions around Greenland was a necessary condition for the ice sheet to disappear almost completely.
  • Item
    Synthesizing long-term sea level rise projections – the MAGICC sea level model v2.0
    (München : European Geopyhsical Union, 2017) Nauels, Alexander; Meinshausen, Malte; Mengel, Matthias; Lorbacher, Katja; Wigley, Tom M.L.
    Sea level rise (SLR) is one of the major impacts of global warming; it will threaten coastal populations, infrastructure, and ecosystems around the globe in coming centuries. Well-constrained sea level projections are needed to estimate future losses from SLR and benefits of climate protection and adaptation. Process-based models that are designed to resolve the underlying physics of individual sea level drivers form the basis for state-of-the-art sea level projections. However, associated computational costs allow for only a small number of simulations based on selected scenarios that often vary for different sea level components. This approach does not sufficiently support sea level impact science and climate policy analysis, which require a sea level projection methodology that is flexible with regard to the climate scenario yet comprehensive and bound by the physical constraints provided by process-based models. To fill this gap, we present a sea level model that emulates global-mean long-term process-based model projections for all major sea level components. Thermal expansion estimates are calculated with the hemispheric upwelling-diffusion ocean component of the simple carbon-cycle climate model MAGICC, which has been updated and calibrated against CMIP5 ocean temperature profiles and thermal expansion data. Global glacier contributions are estimated based on a parameterization constrained by transient and equilibrium process-based projections. Sea level contribution estimates for Greenland and Antarctic ice sheets are derived from surface mass balance and solid ice discharge parameterizations reproducing current output from ice-sheet models. The land water storage component replicates recent hydrological modeling results. For 2100, we project 0.35 to 0.56m (66% range) total SLR based on the RCP2.6 scenario, 0.45 to 0.67m for RCP4.5, 0.46 to 0.71m for RCP6.0, and 0.65 to 0.97m for RCP8.5. These projections lie within the range of the latest IPCC SLR estimates. SLR projections for 2300 yield median responses of 1.02m for RCP2.6, 1.76m for RCP4.5, 2.38m for RCP6.0, and 4.73m for RCP8.5. The MAGICC sea level model provides a flexible and efficient platform for the analysis of major scenario, model, and climate uncertainties underlying long-term SLR projections. It can be used as a tool to directly investigate the SLR implications of different mitigation pathways and may also serve as input for regional SLR assessments via component-wise sea level pattern scaling.
  • Item
    Greenland ice sheet model parameters constrained using simulations of the Eemian Interglacial
    (München : European Geopyhsical Union, 2011) Robinson, A.; Calov, R.; Ganopolski, A.
    Using a new approach to force an ice sheet model, we performed an ensemble of simulations of the Greenland Ice Sheet evolution during the last two glacial cycles, with emphasis on the Eemian Interglacial. This ensemble was generated by perturbing four key parameters in the coupled regional climate-ice sheet model and by introducing additional uncertainty in the prescribed "background" climate change. The sensitivity of the surface melt model to climate change was determined to be the dominant driver of ice sheet instability, as reflected by simulated ice sheet loss during the Eemian Interglacial period. To eliminate unrealistic parameter combinations, constraints from present-day and paleo information were applied. The constraints include (i) the diagnosed present-day surface mass balance partition between surface melting and ice discharge at the margin, (ii) the modeled present-day elevation at GRIP; and (iii) the modeled elevation reduction at GRIP during the Eemian. Using these three constraints, a total of 360 simulations with 90 different model realizations were filtered down to 46 simulations and 20 model realizations considered valid. The paleo constraint eliminated more sensitive melt parameter values, in agreement with the surface mass balance partition assumption. The constrained simulations resulted in a range of Eemian ice loss of 0.4–4.4 m sea level equivalent, with a more likely range of about 3.7–4.4 m sea level if the GRIP δ18O isotope record can be considered an accurate proxy for the precipitation-weighted annual mean temperatures.
  • Item
    Heinrich event 1: An example of dynamical ice-sheet reaction to oceanic changes
    (München : European Geopyhsical Union, 2011) Álvarez-Solas, J.; Montoya, M.; Ritz, C.; Ramstein, G.; Charbit, S.; Dumas, C.; Nisancioglu, K.; Dokken, T.; Ganopolski, A.
    Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.