Open Access
Refine
Language
- English (28)
Is part of the Bibliography
- yes (28)
Keywords
- West Antarctica (3)
- Greenland (2)
- climate (2)
- climate change (2)
- climate-change (2)
- collapse (2)
- model (2)
- ocean (2)
- sea-level rise (2)
- shelf (2)
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.
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.
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.
Here we report on a cyclic, physical ice-discharge instability in the Parallel Ice Sheet Model, simulating the flow of a three-dimensional, inherently buttressed ice-sheet-shelf system which periodically surges on a millennial timescale. The thermomechanically coupled model on 1 km horizontal resolution includes an enthalpy-based formulation of the thermodynamics, a nonlinear stress-balance-based sliding law and a very simple subglacial hydrology. The simulated unforced surging is characterized by rapid ice streaming through a bed trough, resulting in abrupt discharge of ice across the grounding line which is eventually calved into the ocean. We visualize the central feedbacks that dominate the subsequent phases of ice buildup, surge and stabilization which emerge from the interaction between ice dynamics, thermodynamics and the subglacial till layer. Results from the variation of surface mass balance and basal roughness suggest that ice sheets of medium thickness may be more susceptible to surging than relatively thin or thick ones for which the surge feedback loop is damped. We also investigate the influence of different basal sliding laws (ranging from purely plastic to nonlinear to linear) on possible surging. The presented mechanisms underlying our simulations of self-maintained, periodic ice growth and destabilization may play a role in large-scale ice-sheet surging, such as the surging of the Laurentide Ice Sheet, which is associated with Heinrich events, and ice-stream shutdown and reactivation, such as observed in the Siple Coast region of West Antarctica.
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.
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 - (CMIP)-M-4). 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 degrees 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 degrees C by 2300.
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.
In order to provide probabilistic projections of the future evolution of the Atlantic Meridional Overturning Circulation (AMOC), we calibrated a simple Stommel-type box model to emulate the output of fully coupled three-dimensional 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 degrees 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 degrees 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.
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 degrees 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 mmyr(-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.