Browsing by Author "Levermann, A."
Now showing 1 - 20 of 23
Results Per Page
Sort Options
- ItemAdaptation required to preserve future high-end river flood risk at present levels(Washington : American Association for the Advancement of Science (A A A S), 2018) Willner, S.N.; Levermann, A.; Zhao, F.; Frieler, K.Earth’s surface temperature will continue to rise for another 20 to 30 years even with the strongest carbon emission reduction currently considered. The associated changes in rainfall patterns can result in an increased flood risk worldwide. We compute the required increase in flood protection to keep high-end fluvial flood risk at present levels. The analysis is carried out worldwide for subnational administrative units. Most of the United States, Central Europe, and Northeast and West Africa, as well as large parts of India and Indonesia, require the strongest adaptation effort. More than half of the United States needs to at least double their protection within the next two decades. Thus, the need for adaptation to increased river flood is a global problem affecting industrialized regions as much as developing countries.
- ItemClimate change under a scenario near 1.5° C of global warming: Monsoon intensification, ocean warming and steric sea level rise(München : European Geopyhsical Union, 2011) Schewe, J.; Levermann, A.; Meinshausen, M.We present climatic consequences of the Representative Concentration Pathways (RCPs) using the coupled climate model CLIMBER-3α, which contains a statistical-dynamical atmosphere and a three-dimensional ocean model. We compare those with emulations of 19 state-of-the-art atmosphere-ocean general circulation models (AOGCM) using MAGICC6. The RCPs are designed as standard scenarios for the forthcoming IPCC Fifth Assessment Report to span the full range of future greenhouse gas (GHG) concentrations pathways currently discussed. The lowest of the RCP scenarios, RCP3-PD, is projected in CLIMBER-3α to imply a maximal warming by the middle of the 21st century slightly above 1.5 °C and a slow decline of temperatures thereafter, approaching today's level by 2500. We identify two mechanisms that slow down global cooling after GHG concentrations peak: The known inertia induced by mixing-related oceanic heat uptake; and a change in oceanic convection that enhances ocean heat loss in high latitudes, reducing the surface cooling rate by almost 50%. Steric sea level rise under the RCP3-PD scenario continues for 200 years after the peak in surface air temperatures, stabilizing around 2250 at 30 cm. This contrasts with around 1.3 m of steric sea level rise by 2250, and 2 m by 2500, under the highest scenario, RCP8.5. Maximum oceanic warming at intermediate depth (300–800 m) is found to exceed that of the sea surface by the second half of the 21st century under RCP3-PD. This intermediate-depth warming persists for centuries even after surface temperatures have returned to present-day values, with potential consequences for marine ecosystems, oceanic methane hydrates, and ice-shelf stability. Due to an enhanced land-ocean temperature contrast, all scenarios yield an intensification of monsoon rainfall under global warming.
- ItemConsistent increase in Indian monsoon rainfall and its variability across CMIP-5 models(München : European Geopyhsical Union, 2013) Menon, A.; Levermann, A.; Schewe, J.; Lehmann, J.; Frieler, K.The possibility of an impact of global warming on the Indian monsoon is of critical importance for the large population of this region. Future projections within the Coupled Model Intercomparison Project Phase 3 (CMIP-3) showed a wide range of trends with varying magnitude and sign across models. Here the Indian summer monsoon rainfall is evaluated in 20 CMIP-5 models for the period 1850 to 2100. In the new generation of climate models, a consistent increase in seasonal mean rainfall during the summer monsoon periods arises. All models simulate stronger seasonal mean rainfall in the future compared to the historic period under the strongest warming scenario RCP-8.5. Increase in seasonal mean rainfall is the largest for the RCP-8.5 scenario compared to other RCPs. Most of the models show a northward shift in monsoon circulation by the end of the 21st century compared to the historic period under the RCP-8.5 scenario. The interannual variability of the Indian monsoon rainfall also shows a consistent positive trend under unabated global warming. Since both the long-term increase in monsoon rainfall as well as the increase in interannual variability in the future is robust across a wide range of models, some confidence can be attributed to these projected trends.
- ItemA critical humidity threshold for monsoon transitions(München : European Geopyhsical Union, 2012) Schewe, J.; Levermann, A.; Cheng, H.Monsoon systems around the world are governed by the so-called moisture-advection feedback. Here we show that, in a minimal conceptual model, this feedback implies a critical threshold with respect to the atmospheric specific humidity qo over the ocean adjacent to the monsoon region. If qo falls short of this critical value qoc, monsoon rainfall over land cannot be sustained. Such a case could occur if evaporation from the ocean was reduced, e.g. due to low sea surface temperatures. Within the restrictions of the conceptual model, we estimate qoc from present-day reanalysis data for four major monsoon systems, and demonstrate how this concept can help understand abrupt variations in monsoon strength on orbital timescales as found in proxy records.
- ItemDelaying future sea-level rise by storing water in Antarctica(München : European Geopyhsical Union, 2016) Frieler, K.; Mengel, M.; Levermann, A.Even if greenhouse gas emissions were stopped today, sea level would continue to rise for centuries, with the long-term sea-level commitment of a 2 °C warmer world significantly exceeding 2 m. In view of the potential implications for coastal populations and ecosystems worldwide, we investigate, from an ice-dynamic perspective, the possibility of delaying sea-level rise by pumping ocean water onto the surface of the Antarctic ice sheet. We find that due to wave propagation ice is discharged much faster back into the ocean than would be expected from a pure advection with surface velocities. The delay time depends strongly on the distance from the coastline at which the additional mass is placed and less strongly on the rate of sea-level rise that is mitigated. A millennium-scale storage of at least 80 % of the additional ice requires placing it at a distance of at least 700 km from the coastline. The pumping energy required to elevate the potential energy of ocean water to mitigate the currently observed 3 mm yr−1 will exceed 7 % of the current global primary energy supply. At the same time, the approach offers a comprehensive protection for entire coastlines particularly including regions that cannot be protected by dikes.
- ItemEmulating Atlantic overturning strength for low emission scenarios: Consequences for sea-level rise along the North American east coast(München : European Geopyhsical Union, 2011) Schleussner, C.F.; Frieler, K.; Meinshausen, M.; Yin, J.; Levermann, A.In order to provide probabilistic projections of the future evolution of the Atlantic Meridional Overturning Circulation (AMOC), we calibrated a simple Stommeltype box model to emulate the output of fully coupled threedimensional atmosphere-ocean general circulation models (AOGCMs) of the Coupled Model Intercomparison Project (CMIP). Based on this calibration to idealised global warming scenarios with and without interactive atmosphere-ocean fluxes and freshwater perturbation simulations, we project the future evolution of the AMOC mean strength within the covered calibration range for the lower two Representative Concentration Pathways (RCPs) until 2100 obtained from the reduced complexity carbon cycle-climate model MAGICC 6. For RCP3-PD with a global mean temperature median below 1.0 C warming relative to the year 2000, we project an ensemble median weakening of up to 11% compared to 22% under RCP4.5 with a warming median up to 1.9 C over the 21st century. Additional Greenland meltwater of 10 and 20 cm of global sea-level rise equivalent further weakens the AMOC by about 4.5 and 10 %, respectively. By combining our outcome with a multi-model sea-level rise study we project a dynamic sea-level rise along the New York City coastline of 4 cm for the RCP3-PD and of 8 cm for the RCP4.5 scenario over the 21st century. We estimate the total steric and dynamic sea-level rise for New York City to be about 24 cm until 2100 for the RCP3-PD scenario, which can hold as a lower bound for sea-level rise projections in this region, as it does not include ice sheet and mountain glacier contributions.
- ItemEnhanced Atlantic subpolar gyre variability through baroclinic threshold in a coarse resolution model(München : European Geopyhsical Union, 2012) Mengel, M.; Levermann, A.; Schleussner, C.-F.; Born, A.Direct observations, satellite measurements and paleo records reveal strong variability in the Atlantic subpolar gyre on various time scales. Here we show that variations of comparable amplitude can only be simulated in a coupled climate model in the proximity of a dynamical threshold. The threshold and the associated dynamic response is due to a positive feedback involving increased salt transport in the subpolar gyre and enhanced deep convection in its centre. A series of sensitivity experiments is performed with a coarse resolution ocean general circulation model coupled to a statistical-dynamical atmosphere model which in itself does not produce atmospheric variability. To simulate the impact of atmospheric variability, the model system is perturbed with freshwater forcing of varying, but small amplitude and multi-decadal to centennial periodicities and observational variations in wind stress. While both freshwater and wind-stress-forcing have a small direct effect on the strength of the subpolar gyre, the magnitude of the gyre's response is strongly increased in the vicinity of the threshold. Our results indicate that baroclinic self-amplification in the North Atlantic ocean can play an important role in presently observed SPG variability and thereby North Atlantic climate variability on multi-decadal scales.
- ItemEstimating the near-surface permafrost-carbon feedback on global warming(München : European Geopyhsical Union, 2012) Schneider von Deimling, T.; Meinshausen, M.; Levermann, A.; Huber, V.; Frieler, K.; Lawrence, D.M.; Brovkin, V.Thawing of permafrost and the associated release of carbon constitutes a positive feedback in the climate system, elevating the effect of anthropogenic GHG emissions on global-mean temperatures. Multiple factors have hindered the quantification of this feedback, which was not included in climate carbon-cycle models which participated in recent model intercomparisons (such as the Coupled Carbon Cycle Climate Model Intercomparison Project – C4MIP) . There are considerable uncertainties in the rate and extent of permafrost thaw, the hydrological and vegetation response to permafrost thaw, the decomposition timescales of freshly thawed organic material, the proportion of soil carbon that might be emitted as carbon dioxide via aerobic decomposition or as methane via anaerobic decomposition, and in the magnitude of the high latitude amplification of global warming that will drive permafrost degradation. Additionally, there are extensive and poorly characterized regional heterogeneities in soil properties, carbon content, and hydrology. Here, we couple a new permafrost module to a reduced complexity carbon-cycle climate model, which allows us to perform a large ensemble of simulations. The ensemble is designed to span the uncertainties listed above and thereby the results provide an estimate of the potential strength of the feedback from newly thawed permafrost carbon. For the high CO2 concentration scenario (RCP8.5), 33–114 GtC (giga tons of Carbon) are released by 2100 (68 % uncertainty range). This leads to an additional warming of 0.04–0.23 °C. Though projected 21st century permafrost carbon emissions are relatively modest, ongoing permafrost thaw and slow but steady soil carbon decomposition means that, by 2300, about half of the potentially vulnerable permafrost carbon stock in the upper 3 m of soil layer (600–1000 GtC) could be released as CO2, with an extra 1–4 % being released as methane. Our results also suggest that mitigation action in line with the lower scenario RCP3-PD could contain Arctic temperature increase sufficiently that thawing of the permafrost area is limited to 9–23 % and the permafrost-carbon induced temperature increase does not exceed 0.04–0.16 °C by 2300.
- ItemFracture-induced softening for large-scale ice dynamics(München : European Geopyhsical Union, 2014) Albrecht, T.; Levermann, A.Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e.g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.
- ItemA framework for the cross-sectoral integration of multi-model impact projections: Land use decisions under climate impacts uncertainties(München : European Geopyhsical Union, 2015) Frieler, K.; Levermann, A.; Elliott, J.; Heinke, J.; Arneth, A.; Bierkens, M.F.P.; Ciais, P.; Clark, D.B.; Deryng, D.; Döll, P.; Falloon, P.; Fekete, B.; Folberth, C.; Friend, A.D.; Gellhorn, C.; Gosling, S.N.; Haddeland, I.; Khabarov, N.; Lomas, M.; Masaki, Y.; Nishina, K.; Neumann, K.; Oki, T.; Pavlick, R.; Ruane, A.C.; Schmid, E.; Schmitz, C.; Stacke, T.; Stehfest, E.; Tang, Q.; Wisser, D.; Huber, V.; Piontek, F.; Warszawski, L.; Schewe, J.; Lotze-Campen, H.; Schellnhuber, H.J.Climate change and its impacts already pose considerable challenges for societies that will further increase with global warming (IPCC, 2014a, b). Uncertainties of the climatic response to greenhouse gas emissions include the potential passing of large-scale tipping points (e.g. Lenton et al., 2008; Levermann et al., 2012; Schellnhuber, 2010) and changes in extreme meteorological events (Field et al., 2012) with complex impacts on societies (Hallegatte et al., 2013). Thus climate change mitigation is considered a necessary societal response for avoiding uncontrollable impacts (Conference of the Parties, 2010). On the other hand, large-scale climate change mitigation itself implies fundamental changes in, for example, the global energy system. The associated challenges come on top of others that derive from equally important ethical imperatives like the fulfilment of increasing food demand that may draw on the same resources. For example, ensuring food security for a growing population may require an expansion of cropland, thereby reducing natural carbon sinks or the area available for bio-energy production. So far, available studies addressing this problem have relied on individual impact models, ignoring uncertainty in crop model and biome model projections. Here, we propose a probabilistic decision framework that allows for an evaluation of agricultural management and mitigation options in a multi-impact-model setting. Based on simulations generated within the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP), we outline how cross-sectorally consistent multi-model impact simulations could be used to generate the information required for robust decision making. Using an illustrative future land use pattern, we discuss the trade-off between potential gains in crop production and associated losses in natural carbon sinks in the new multiple crop- and biome-model setting. In addition, crop and water model simulations are combined to explore irrigation increases as one possible measure of agricultural intensification that could limit the expansion of cropland required in response to climate change and growing food demand. This example shows that current impact model uncertainties pose an important challenge to long-term mitigation planning and must not be ignored in long-term strategic decision making.
- ItemFuture changes in extratropical storm tracks and baroclinicity under climate change(Bristol : IOP, 2014) Lehmann, J.; Coumou, D.; Frieler, K.; Eliseev, A.V.; Levermann, A.The weather in Eurasia, Australia, and North and South America is largely controlled by the strength and position of extratropical storm tracks. Future climate change will likely affect these storm tracks and the associated transport of energy, momentum, and water vapour. Many recent studies have analyzed how storm tracks will change under climate change, and how these changes are related to atmospheric dynamics. However, there are still discrepancies between different studies on how storm tracks will change under future climate scenarios. Here, we show that under global warming the CMIP5 ensemble of coupled climate models projects only little relative changes in vertically averaged mid-latitude mean storm track activity during the northern winter, but agree in projecting a substantial decrease during summer. Seasonal changes in the Southern Hemisphere show the opposite behaviour, with an intensification in winter and no change during summer. These distinct seasonal changes in northern summer and southern winter storm tracks lead to an amplified seasonal cycle in a future climate. Similar changes are seen in the mid-latitude mean Eady growth rate maximum, a measure that combines changes in vertical shear and static stability based on baroclinic instability theory. Regression analysis between changes in the storm tracks and changes in the maximum Eady growth rate reveal that most models agree in a positive association between the two quantities over mid-latitude regions.
- ItemA High-End Estimate of Sea Level Rise for Practitioners(Hoboken, NJ : Wiley-Blackwell, 2022) van de Wal, R.S.W.; Nicholls, R J.; Behar, D.; McInnes, K.; Stammer, D.; Lowe, J.A.; Church, J.A.; DeConto, R.; Fettweis, X.; Goelzer, H.; Haasnoot, M.; Haigh, I.D.; Hinkel, J.; Horton, B.P.; James, T.S.; Jenkins, A.; LeCozannet, G.; Levermann, A.; Lipscomb, W.H.; Marzeion, B.; Pattyn, F.; Payne, A.J.; Pfeffer, W.T.; Price, S.F.; Seroussi, H.; Sun, S.; Veatch, W.; White, K.Sea level rise (SLR) is a long-lasting consequence of climate change because global anthropogenic warming takes centuries to millennia to equilibrate for the deep ocean and ice sheets. SLR projections based on climate models support policy analysis, risk assessment and adaptation planning today, despite their large uncertainties. The central range of the SLR distribution is estimated by process-based models. However, risk-averse practitioners often require information about plausible future conditions that lie in the tails of the SLR distribution, which are poorly defined by existing models. Here, a community effort combining scientists and practitioners builds on a framework of discussing physical evidence to quantify high-end global SLR for practitioners. The approach is complementary to the IPCC AR6 report and provides further physically plausible high-end scenarios. High-end estimates for the different SLR components are developed for two climate scenarios at two timescales. For global warming of +2°C in 2100 (RCP2.6/SSP1-2.6) relative to pre-industrial values our high-end global SLR estimates are up to 0.9 m in 2100 and 2.5 m in 2300. Similarly, for a (RCP8.5/SSP5-8.5), we estimate up to 1.6 m in 2100 and up to 10.4 m in 2300. The large and growing differences between the scenarios beyond 2100 emphasize the long-term benefits of mitigation. However, even a modest 2°C warming may cause multi-meter SLR on centennial time scales with profound consequences for coastal areas. Earlier high-end assessments focused on instability mechanisms in Antarctica, while here we emphasize the importance of the timing of ice shelf collapse around Antarctica. This is highly uncertain due to low understanding of the driving processes. Hence both process understanding and emission scenario control high-end SLR.
- ItemInteraction of marine ice-sheet instabilities in two drainage basins: Simple scaling of geometry and transition time(München : European Geopyhsical Union, 2015) Feldmann, J.; Levermann, A.The initiation of a marine ice-sheet instability (MISI) is generally discussed from the ocean side of the ice sheet. It has been shown that the reduction in ice-shelf buttressing and softening of the coastal ice can destabilize a marine ice sheet if the bedrock is sloping upward towards the ocean. Using a conceptional flow-line geometry, we investigate the possibility of whether a MISI can be triggered from the direction of the ice divide as opposed to coastal forcing and explore the interaction between connected basins. We find that the initiation of a MISI in one basin can induce a destabilization in the other. The underlying mechanism of basin interaction is based on dynamic thinning and a consecutive motion of the ice divide which induces a thinning in the adjacent basin and a successive initiation of the instability. Our simplified and symmetric topographic setup allows scaling both the geometry and the transition time between both instabilities. We find that the ice profile follows a universal shape that is scaled with the horizontal extent of the ice sheet and that the same exponent of 1/2 applies for the scaling relation between central surface elevation and horizontal extent as in the pure shallow ice approximation (Vialov profile). Altering the central bed elevation, we find that the extent of grounding-line retreat in one basin determines the degree of interaction with the other. Different scenarios of basin interaction are discussed based on our modeling results as well as on a conceptual flux-balance analysis. We conclude that for the three-dimensional case, the possibility of drainage basin interaction on timescales on the order of 1 kyr or larger cannot be excluded and hence needs further investigation.
- ItemKinematic first-order calving law implies potential for abrupt ice-shelf retreat(München : European Geopyhsical Union, 2012) Levermann, A.; Albrecht, T.; Winkelmann, R.; Martin, M.A.; Haseloff, M.; Joughin, I.Recently observed large-scale disintegration of Antarctic ice shelves has moved their fronts closer towards grounded ice. In response, ice-sheet discharge into the ocean has accelerated, contributing to global sea-level rise and emphasizing the importance of calving-front dynamics. The position of the ice front strongly influences the stress field within the entire sheet-shelf-system and thereby the mass flow across the grounding line. While theories for an advance of the ice-front are readily available, no general rule exists for its retreat, making it difficult to incorporate the retreat in predictive models. Here we extract the first-order large-scale kinematic contribution to calving which is consistent with large-scale observation. We emphasize that the proposed equation does not constitute a comprehensive calving law but represents the first-order kinematic contribution which can and should be complemented by higher order contributions as well as the influence of potentially heterogeneous material properties of the ice. When applied as a calving law, the equation naturally incorporates the stabilizing effect of pinning points and inhibits ice shelf growth outside of embayments. It depends only on local ice properties which are, however, determined by the full topography of the ice shelf. In numerical simulations the parameterization reproduces multiple stable fronts as observed for the Larsen A and B Ice Shelves including abrupt transitions between them which may be caused by localized ice weaknesses. We also find multiple stable states of the Ross Ice Shelf at the gateway of the West Antarctic Ice Sheet with back stresses onto the sheet reduced by up to 90 % compared to the present state.
- ItemMechanism for potential strengthening of Atlantic overturning prior to collapse(München : European Geopyhsical Union, 2014) Ehlert, D.; Levermann, A.The Atlantic meridional overturning circulation (AMOC) carries large amounts of heat into the North Atlantic influencing climate regionally as well as globally. Palaeo-records and simulations with comprehensive climate models suggest that the positive salt-advection feedback may yield a threshold behaviour of the system. That is to say that beyond a certain amount of freshwater flux into the North Atlantic, no meridional overturning circulation can be sustained. Concepts of monitoring the AMOC and identifying its vicinity to the threshold rely on the fact that the volume flux defining the AMOC will be reduced when approaching the threshold. Here we advance conceptual models that have been used in a paradigmatic way to understand the AMOC, by introducing a density-dependent parameterization for the Southern Ocean eddies. This additional degree of freedom uncovers a mechanism by which the AMOC can increase with additional freshwater flux into the North Atlantic, before it reaches the threshold and collapses: an AMOC that is mainly wind-driven will have a constant upwelling as long as the Southern Ocean winds do not change significantly. The downward transport of tracers occurs either in the northern sinking regions or through Southern Ocean eddies. If freshwater is transported, either atmospherically or via horizontal gyres, from the low to high latitudes, this would reduce the eddy transport and by continuity increase the northern sinking which defines the AMOC until a threshold is reached at which the AMOC cannot be sustained. If dominant in the real ocean this mechanism would have significant consequences for monitoring the AMOC.
- ItemParameterization for subgrid-scale motion of ice-shelf calving fronts(München : European Geopyhsical Union, 2011) Albrecht, T.; Martin, M.; Haseloff, M.; Winkelmann, R.; Levermann, A.In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS.
- ItemThe Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 1: Model description(München : European Geopyhsical Union, 2011) Winkelmann, R.; Martin, M.A.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.We present the Potsdam Parallel Ice Sheet Model (PISM-PIK), developed at the Potsdam Institute for Climate Impact Research to be used for simulations of large-scale ice sheet-shelf systems. It is derived from the Parallel Ice Sheet Model (Bueler and Brown, 2009). Velocities are calculated by superposition of two shallow stress balance approximations within the entire ice covered region: the shallow ice approximation (SIA) is dominant in grounded regions and accounts for shear deformation parallel to the geoid. The plug-flow type shallow shelf approximation (SSA) dominates the velocity field in ice shelf regions and serves as a basal sliding velocity in grounded regions. Ice streams can be identified diagnostically as regions with a significant contribution of membrane stresses to the local momentum balance. All lateral boundaries in PISM-PIK are free to evolve, including the grounding line and ice fronts. Ice shelf margins in particular are modeled using Neumann boundary conditions for the SSA equations, reflecting a hydrostatic stress imbalance along the vertical calving face. The ice front position is modeled using a subgrid-scale representation of calving front motion (Albrecht et al., 2011) and a physically-motivated calving law based on horizontal spreading rates. The model is tested in experiments from the Marine Ice Sheet Model Intercomparison Project (MISMIP). A dynamic equilibrium simulation of Antarctica under present-day conditions is presented in Martin et al. (2011).
- ItemThe Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet(München : European Geopyhsical Union, 2011) Martin, M.A.; Winkelmann, R.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model.
- ItemProjecting Antarctic ice discharge using response functions from SeaRISE ice-sheet models(München : European Geopyhsical Union, 2014) Levermann, A.; Winkelmann, R.; Nowicki, S.; Fastook, J.L.; Frieler, K.; Greve, R.; Hellmer, H.H.; Martin, M.A.; Meinshausen, M.; Mengel, M.; Payne, A.J.; Pollard, D.; Sato, T.; Timmermann, R.; Wang, W.L.; Bindschadler, R.A.The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02–0.14 m; 90% range: 0.0–0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04–0.21 m; 90% range: 0.01–0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66% range: 0.04–0.17 m; 90% range: 0.02–0.25 m) for RCP-2.6 and 0.15 m (66% range: 0.07–0.28 m; 90% range: 0.04–0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.
- ItemThe role of the North Atlantic overturning and deep ocean for multi-decadal global-mean-temperature variability(München : European Geopyhsical Union, 2014) Schleussner, C.F.; Runge, J.; Lehmann, J.; Levermann, A.Earth's climate exhibits internal modes of variability on various timescales. Here we investigate multi-decadal variability of the Atlantic meridional overturning circulation (AMOC), Northern Hemisphere sea-ice extent and global mean temperature (GMT) in an ensemble of CMIP5 models under control conditions. We report an inter-annual GMT variability of about ±0.1° C originating solely from natural variability in the model ensemble. By decomposing the GMT variance into contributions of the AMOC and Northern Hemisphere sea-ice extent using a graph-theoretical statistical approach, we find the AMOC to contribute 8% to GMT variability in the ensemble mean. Our results highlight the importance of AMOC sea-ice feedbacks that explain 5% of the GMT variance, while the contribution solely related to the AMOC is found to be about 3%. As a consequence of multi-decadal AMOC variability, we report substantial variations in North Atlantic deep-ocean heat content with trends of up to 0.7 × 1022 J decade−1 that are of the order of observed changes over the last decade and consistent with the reduced GMT warming trend over this period. Although these temperature anomalies are largely density-compensated by salinity changes, we find a robust negative correlation between the AMOC and North Atlantic deep-ocean density with density lagging the AMOC by 5 to 11 yr in most models. While this would in principle allow for a self-sustained oscillatory behavior of the coupled AMOC–deep-ocean system, our results are inconclusive about the role of this feedback in the model ensemble.