Open Access
Refine
Has Fulltext
- no (21) (remove)
Document Type
- Article (21)
Language
- English (21)
Is part of the Bibliography
- yes (21)
Keywords
- climate change (2)
- baroclinicity (1)
- chlimate impacts (1)
- coupled climate model (1)
- cultural heritage (1)
- monsoon failure (1)
- non-linear dynamics (1)
- sea-level rise (1)
- stochastic model (1)
- storm tracks (1)
Institute
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.
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).
Both climate-change damages and climate-change mitigation will incur economic costs. While the risk of severe damages increases with the level of global warming (Dell et al., 2014; IPCC, 2014b, 2018; Lenton et al., 2008), mitigating costs increase steeply with more stringent warming limits (IPCC, 2014a; Luderer et al., 2013; Rogelj et al., 2015). Here, we show that the global warming limit that minimizes this century's total economic costs of climate change lies between 1.9 and 2 ∘C, if temperature changes continue to impact national economic growth rates as observed in the past and if instantaneous growth effects are neither compensated nor amplified by additional growth effects in the following years. The result is robust across a wide range of normative assumptions on the valuation of future welfare and inequality aversion. We combine estimates of climate-change impacts on economic growth for 186 countries (applying an empirical damage function from Burke et al., 2015) with mitigation costs derived from a state-of-the-art energy–economy–climate model with a wide range of highly resolved mitigation options (Kriegler et al., 2017; Luderer et al., 2013, 2015). Our purely economic assessment, even though it omits non-market damages, provides support for the international Paris Agreement on climate change. The political goal of limiting global warming to “well below 2 degrees” is thus also an economically optimal goal given above assumptions on adaptation and damage persistence.
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.
A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.
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.
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.
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.