Search

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

Changes in ocean-driven basal melting have a key influence on the stability of ice shelves, the mass loss from the ice sheet, ocean circulation, and global sea level rise. Coupled ice sheet–ocean models play a critical role in understanding future ice sheet evolution and examining the processes governing ice sheet responses to basal melting. However, as a new approach, coupled ice sheet–ocean systems come with new challenges, and the impacts of solutions implemented to date have not been investigated. An emergent feature in several contributing coupled models to the 1st Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1) was a time-varying oscillation in basal melt rates. Here, we use a recently developed coupling framework, FISOC (v1.1), to connect the modified ocean model ROMSIceShelf (v1.0) and ice sheet model Elmer/Ice (v9.0), to investigate the origin and implications of the feature and, more generally, the impact of coupled modeling strategies on the simulated basal melt in an idealized ice shelf cavity based on the MISOMIP setup. We found the spatial-averaged basal melt rates (3.56 m yr−1) oscillated with an amplitude ∼0.7 m yr−1 and approximate period of ∼6 years between year 30 and 100 depending on the experimental design. The melt oscillations emerged in the coupled system and the standalone ocean model using a prescribed change of cavity geometry. We found that the oscillation feature is closely related to the discretized ungrounding of the ice sheet, exposing new ocean, and is likely strengthened by a combination of positive buoyancy–melt feedback and/or melt–geometry feedback near the grounding line, and the frequent coupling of ice geometry and ocean evolution. Sensitivity tests demonstrate that the oscillation feature is always present, regardless of the choice of coupling interval, vertical resolution in the ocean model, tracer properties of cells ungrounded by the retreating ice sheet, or the dependency of friction velocities to the vertical resolution. However, the amplitude, phase, and sub-cycle variability of the oscillation varied significantly across the different configurations. We were unable to ultimately determine whether the feature arises purely due to numerical issues (related to discretization) or a compounding of multiple physical processes amplifying a numerical artifact. We suggest a pathway and choices of physical parameters to help other efforts understand the coupled ice sheet–ocean system using numerical models.

Ice shelves play a key role in the dynamics of marine ice sheets by buttressing grounded ice and limiting rates of ice flux to the oceans. In response to recent climatic and oceanic change, ice shelves fringing the West Antarctic Ice Sheet (WAIS) have begun to fragment and retreat, with major implications for ice-sheet stability. Here, we focus on the Thwaites Eastern Ice Shelf (TEIS), the remaining pinned floating extension of Thwaites Glacier. We show that TEIS has undergone a process of fragmentation in the last 5 years, including brittle failure along a major shear zone, formation of tensile cracks on the main body of the shelf, and a release of tabular bergs on both the eastern and western flanks. Simulations with the Helsinki Discrete Element Model (HiDEM) show that this pattern of failure is associated with high backstress from a submarine pinning point at the distal edge of the shelf. We show that a significant zone of shear, upstream of the main pinning point, developed in response to the rapid acceleration of the shelf between 2002 and 2006, seeding damage on the shelf. Subsequently, basal melting and positive feedback between damage and strain rates weakened TEIS, allowing damage to accumulate. Thus, although backstress on TEIS has likely diminished over time as the pinning point shrunk, accumulation of damage has ensured that the ice in the shear zone remained the weakest link in the system. Experiments with the BISICLES ice-sheet model indicate that additional damage to or unpinning of TEIS is unlikely to trigger significantly increased ice loss from WAIS, but the calving response to the loss of TEIS remains highly uncertain. It is widely recognised that ice-shelf fragmentation and collapse can be triggered by hydrofracturing and/or unpinning from ice-shelf margins or grounding points. Our results indicate a third mechanism, backstress triggered failure, that can occur if and when an ice shelf is no longer able to withstand stress imposed by pinning points. In most circumstances, pinning points are essential for ice-shelf stability, but as ice shelves thin and weaken, the concentration of backstress in damaged ice upstream of a pinning point may provide the seeds of their demise.

The Amery Ice Shelf (AIS), East Antarctica, has a layered structure, due to the presence of both meteoric and marine ice. In this study, the thermal structure of the AIS and its spatial pattern are evaluated and analysed through borehole observations and numerical simulations with Elmer/Ice, a full-Stokes ice sheet/shelf model. In the area with marine ice, a near-isothermal basal layer up to 120 m thick is observed, which closely conforms to the pressure-dependent freezing temperature of seawater. In the area experiencing basal melting, large temperature gradients, up to −0.36 ∘C m−1, are observed at the base. Three-dimensional (3-D) steady-state temperature simulations with four different basal mass balance (BMB) datasets for the AIS reveal a high sensitivity of ice shelf thermal structure to the distribution of BMB. We also construct a one-dimensional (1-D) transient temperature column model to simulate the process of an ice column moving along a flowline with corresponding boundary conditions, which achieves slightly better agreement with borehole observations than the 3-D simulations. Our simulations reveal internal cold ice advected from higher elevations by the AIS's main inlet glaciers, warming downstream along the ice flow, and we suggest the thermal structures dominated by these cold cores may commonly exist among Antarctic ice shelves. For the marine ice, the porous structure of its lower layer and interactions with ocean below determine the local thermal regime and give rise to the near-isothermal phenomenon. The limitations in our simulations identify the need for ice shelf–ocean coupled models with improved thermodynamics and more comprehensive boundary conditions. Given the temperature dependence of ice rheology, the depth-averaged ice stiffness factor B(T′)‾ derived from the most realistic simulated temperature field is presented to quantify the influence of the temperature distribution on ice shelf dynamics. The full 3-D temperature field provides a useful input to future modelling studies.

A number of important questions concern processes at the margins of ice sheets where multiple components of the Earth system, most crucially ice sheets and oceans, interact. Such processes include thermodynamic interaction at the ice–ocean interface, the impact of meltwater on ice shelf cavity circulation, the impact of basal melting of ice shelves on grounded ice dynamics and ocean controls on iceberg calving. These include fundamentally coupled processes in which feedback mechanisms between ice and ocean play an important role. Some of these mechanisms have major implications for humanity, most notably the impact of retreating marine ice sheets on the global sea level. In order to better quantify these mechanisms using computer models, feedbacks need to be incorporated into the modelling system. To achieve this, ocean and ice dynamic models must be coupled, allowing runtime information sharing between components. We have developed a flexible coupling framework based on existing Earth system coupling technologies. The open-source Framework for Ice Sheet–Ocean Coupling (FISOC) provides a modular approach to coupling, facilitating switching between different ice dynamic and ocean components. FISOC allows fully synchronous coupling, in which both ice and ocean run on the same time step, or semi-synchronous coupling in which the ice dynamic model uses a longer time step. Multiple regridding options are available, and there are multiple methods for coupling the sub-ice-shelf cavity geometry. Thermodynamic coupling may also be activated. We present idealized simulations using FISOC with a Stokes flow ice dynamic model coupled to a regional ocean model. We demonstrate the modularity of FISOC by switching between two different regional ocean models and presenting outputs for both. We demonstrate conservation of mass and other verification steps during evolution of an idealized coupled ice–ocean system, both with and without grounding line movement.

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

The marine-terminating outlet in Basin 3, Austfonna ice cap, has been accelerating since the mid-1990s. Stepwise multi-annual acceleration associated with seasonal summer speed-up events was observed before the outlet entered the basin-wide surge in autumn 2012. We used multiple numerical models to explore hydrologic activation mechanisms for the surge behaviour. A continuum ice dynamic model was used to invert basal friction coefficient distributions using the control method and observed surface velocity data between April 2012 and July 2014. This has provided input to a discrete element model capable of simulating individual crevasses, with the aim of finding locations where meltwater entered the glacier during the summer and reached the bed. The possible flow paths of surface meltwater reaching the glacier bed as well as those of meltwater produced at the bed were calculated according to the gradient of the hydraulic potential.The inverted friction coefficients show the unplugging of the stagnant ice front and expansion of low-friction regions before the surge reached its peak velocity in January 2013. Crevasse distribution reflects the basal friction pattern to a high degree. The meltwater reaches the bed through the crevasses located above the margins of the subglacial valley and the basal melt that is generated mainly by frictional heating flows either to the fast-flowing units or potentially accumulates in an overdeepened region. Based on these results, the mechanisms facilitated by basal meltwater production, crevasse opening and the routing of meltwater to the bed are discussed for the surge in Basin 3.

Floating ice shelves exert a stabilizing force onto the inland ice sheet. However, this buttressing effect is diminished by the fracture process, which on large scales effectively softens the ice, accelerating its flow, increasing calving, and potentially leading to ice shelf breakup. We add a continuum damage model (CDM) to the BISICLES ice sheet model, which is intended to model the localized opening of crevasses under stress, the transport of those crevasses through the ice sheet, and the coupling between crevasse depth and the ice flow field and to carry out idealized numerical experiments examining the broad impact on large-scale ice sheet and shelf dynamics. In each case we see a complex pattern of damage evolve over time, with an eventual loss of buttressing approximately equivalent to halving the thickness of the ice shelf. We find that it is possible to achieve a similar ice flow pattern using a simple rule of thumb: introducing an enhancement factor ∼ 10 everywhere in the model domain. However, spatially varying damage (or equivalently, enhancement factor) fields set at the start of prognostic calculations to match velocity observations, as is widely done in ice sheet simulations, ought to evolve in time, or grounding line retreat can be slowed by an order of magnitude.

The association between exacerbation aetiology and exacerbation frequency is poorly understood. We analysed 2-year follow-up data from a prospective observational study of patients with chronic obstructive pulmonary disease (COPD) (www.clinicaltrials.gov identifier number NCT01360398) to evaluate year-to-year variation in exacerbation frequency and related aetiology. A total of 127 patients underwent blood and sputum sampling monthly and at exacerbation to detect respiratory infections and eosinophilic inflammation; 103 continued into year 2 and 88 completed both years. The most common bacterial species at stable state and exacerbation was Haemophilus influenzae. Among infrequent exacerbators (one exacerbation per year), the incidence of viral infection at exacerbation was high (60.0% (95% CI 35.1–81.7%) in year 1 and 78.6% (53.4–94.2%) in year 2). Those with more frequent exacerbations tended to have higher relative incidence of bacterial than viral infection. Patients with at least two additional exacerbations in year 2 versus year 1 had a higher risk of H. influenzae colonisation at stable state than those with at least two fewer exacerbations, as detected by culture (OR 1.43 (95% CI 0.71–2.91) versus 0.63 (0.40–1.01), p=0.06) and PCR (1.76 (95% CI 0.88–3.51) versus 0.56 (0.37–0.86), p

Understanding the response of fast flowing ice streams or outlet glaciers to changing climate is crucial in order to make reliable projections of sea level change over the coming decades. Motion of fast outlet glaciers occurs largely through basal motion governed by physical processes at the glacier bed, which are not yet fully understood. Various subglacial mechanisms have been suggested for fast flow but common to most of the suggested processes is the requirement of presence of liquid water, and thus temperate conditions. We use a combination of modelling, field, and remote observations in order to study links between different heat sources, the thermal regime and basal sliding in fast flowing areas on Vestfonna ice cap. A special emphasis lies on Franklinbreen, a fast flowing outlet glacier which has been observed to accelerate recently. We use the ice flow model Elmer/Ice including a Weertman type sliding law and a Robin inverse method to infer basal friction parameters from observed surface velocities. Firn heating, i.e. latent heat release through percolation of melt water, is included in our model; its parameterisation is calibrated with the temperature record of a deep borehole. We found that strain heating is negligible, whereas friction heating is identified as one possible trigger for the onset of fast flow. Firn heating is a significant heat source in the central thick and slow flowing area of the ice cap and the essential driver behind the ongoing fast flow in all outlets. Our findings suggest a possible scenario of the onset and maintenance of fast flow on the Vestfonna ice cap based on thermal processes and emphasise the role of latent heat released through refreezing of percolating melt water for fast flow. However, these processes cannot yet be captured in a temporally evolving sliding law. In order to simulate correctly fast flowing outlet glaciers, ice flow models not only need to account fully for all heat sources, but also need to incorporate a sliding law that is not solely based on the basal temperature, but also on hydrology and/or sediment physics.

The outlet glacier of Basin 3 (B3) of Austfonna ice cap, Svalbard, is one of the fastest outlet glaciers in Svalbard, and shows dramatic changes since 1995. In addition to previously observed seasonal summer speed-up associated with the melt season, the winter speed of B3 has accelerated approximately fivefold since 1995. We use the Elmer/Ice full-Stokes model for ice dynamics to infer spatial distributions of basal drag for the winter seasons of 1995, 2008 and 2011. This "inverse" method is based on minimising discrepancy between modelled and observed surface velocities, using satellite remotely sensed velocity fields. We generate steady-state temperature distributions for 1995 and 2011. Frictional heating caused by basal sliding contributes significantly to basal temperatures of the B3 outlet glacier, with heat advection (a longer-timescale process than frictional heating) also being important in the steady state. We present a sensitivity experiment consisting of transient simulations under present-day forcing to demonstrate that using a temporally fixed basal drag field obtained through inversion can lead to thickness change errors of the order of 2 m year−1. Hence it is essential to incorporate the evolution of basal processes in future projections of the evolution of B3. Informed by a combination of our inverse method results and previous studies, we hypothesise a system of processes and feedbacks involving till deformation and basal hydrology to explain both the seasonal accelerations (short residence time pooling of meltwater at the ice–till interface) and the ongoing interannual speed-up (gradual penetration of water into the till, reducing till strength).

Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.