Mechanisms driving variability in the ocean forcing of Pine Island Glacier

Pine Island Glacier (PIG) terminates in a rapidly melting ice shelf, and ocean circulation and temperature are implicated in the retreat and growing contribution to sea level rise of PIG and nearby glaciers. However, the variability of the ocean forcing of PIG has been poorly constrained due to a lack of multi-year observations. Here we show, using a unique record close to the Pine Island Ice Shelf (PIIS), that there is considerable oceanic variability at seasonal and interannual timescales, including a pronounced cold period from October 2011 to May 2013. This variability can be largely explained by two processes: cumulative ocean surface heat fluxes and sea ice formation close to PIIS; and interannual reversals in ocean currents and associated heat transport within Pine Island Bay, driven by a combination of local and remote forcing. Local atmospheric forcing therefore plays an important role in driving oceanic variability close to PIIS

Biophysical Factors Affecting the Distribution of Demersal Fish around the Head of a Submarine Canyon Off the Bonney Coast, South Australia

We sampled the demersal fish community of the Bonney Canyon, South Australia at depths (100–1,500 m) and locations that are poorly known. Seventy-eight species of demersal fish were obtained from 12 depth-stratified trawls along, and to either side, of the central canyon axis. Distributional patterns in species richness and biomass were highly correlated. Three fish assemblage groupings, characterised by small suites of species with narrow depth distributions, were identified on the shelf, upper slope and mid slope. The assemblage groupings were largely explained by depth (ρw = 0.78). Compared to the depth gradient, canyon-related effects are weak or occur at spatial or temporal scales not sampled in this study. A conceptual physical model displayed features consistent with the depth zonational patterns in fish, and also indicated that canyon upwelling can occur. The depth zonation of the fish assemblage was associated with the depth distribution of water masses in the area. Notably, the mid-slope community (1,000 m) coincided with a layer of Antarctic Intermediate Water, the upper slope community (500 m) resided within the core of the Flinders Current, and the shelf community was located in a well-mixed layer of surface water (<450 m depth)

Ice shelf basal melt rates around Antarctica from simulations and observations

We introduce an explicit representation of Antarctic ice shelf cavities in the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) ocean retrospective analysis; and compare resulting basal melt rates and patterns to independent estimates from satellite observations. Two simulations are carried out: The first is based on the original ECCO2 vertical discretization; the second has higher vertical resolution particularly at the depth range of ice shelf cavities. The original ECCO2 vertical discretization produces higher than observed melt rates and leads to a misrepresentation of Southern Ocean water mass properties and transports. In general, thicker levels at the base of the ice shelves lead to increased melting because of their larger heat capacity. This strengthens horizontal gradients and circulation within and outside the cavities and, in turn, warm water transports from the shelf break to the ice shelves. The simulation with more vertical levels produces basal melt rates (1735 ± 164 Gt/a) and patterns that are in better agreement with observations. Thinner levels in the sub-ice-shelf cavities improve the representation of a fresh/cold layer at the ice shelf base and of warm/salty water near the bottom, leading to a sharper pycnocline and reduced vertical mixing underneath the ice shelf. Improved water column properties lead to more accurate melt rates and patterns, especially for melt/freeze patterns under large cold-water ice shelves. At the 18 km grid spacing of the ECCO2 model configuration, the smaller, warm-water ice shelves cannot be properly represented, with higher than observed melt rates in both simulations

Ice shelf basal melt rates around Antarctica from simulations and observations

We introduce an explicit representation of Antarctic ice shelf cavities in the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) ocean retrospective analysis; and compare resulting basal melt rates and patterns to independent estimates from satellite observations. Two simulations are carried out: The first is based on the original ECCO2 vertical discretization; the second has higher vertical resolution particularly at the depth range of ice shelf cavities. The original ECCO2 vertical discretization produces higher than observed melt rates and leads to a misrepresentation of Southern Ocean water mass properties and transports. In general, thicker levels at the base of the ice shelves lead to increased melting because of their larger heat capacity. This strengthens horizontal gradients and circulation within and outside the cavities and, in turn, warm water transports from the shelf break to the ice shelves. The simulation with more vertical levels produces basal melt rates (1735 ± 164 Gt/a) and patterns that are in better agreement with observations. Thinner levels in the sub-ice-shelf cavities improve the representation of a fresh/cold layer at the ice shelf base and of warm/salty water near the bottom, leading to a sharper pycnocline and reduced vertical mixing underneath the ice shelf. Improved water column properties lead to more accurate melt rates and patterns, especially for melt/freeze patterns under large cold-water ice shelves. At the 18 km grid spacing of the ECCO2 model configuration, the smaller, warm-water ice shelves cannot be properly represented, with higher than observed melt rates in both simulations

A damage mechanics assessment of the Larsen B ice shelf prior to collapse: Toward a physically-based calving law

Calving is a primary process of mass ablation for glaciers and ice sheets, though it still eludes a general physical law. Here, we propose a calving framework based on continuum damage mechanics coupled with the equations of viscous deformation of glacier ice. We introduce a scalar damage variable that quantifies the loss of load-bearing surface area due to fractures and that feeds back with ice viscosity to represent fracture-induced softening. The calving law is a standard failure criterion for viscous damaging materials and represents a macroscopic brittle instability quantified by a critical or threshold damage. We constrain this threshold using the Ice Sheet System Model (ISSM) by inverting for damage on the Larsen B ice shelf prior to its 2002 collapse. By analyzing the damage distribution in areas that subsequently calved, we conclude that calving occurs after fractures have reduced the load-bearing capacity of the ice by 60 ± 10%. © 2012. American Geophysical Union. All Rights Reserved

On the influence of adequate Weddell Sea characteristics in a large-scale global ocean circulation model

Global ocean circulation models usually lack an adequate consideration of high-latitude processes due to a limited model domain or insufficient resolution. Without the processes in key areas of the global thermohaline circulation, the characteristics and flow of deep and bottom waters cannot be modeled realistically. In this study, a high-resolution (~20 km) ocean model focused on the Weddell Sea sector of the Southern Ocean is combined with a low-resolution (2°× 2°) global ocean model applying the state estimation technique. Temperature, salinity, and velocity data on two Weddell Sea sections from the regional model are used as constraints for the large-scale model in addition to satellite altimetry and sea-surface temperatures. The differences between the model with additional constraints and without document that the Weddell Sea circulation exerts significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. Furthermore, a warming trend in the period 19932001 was found in the Weddell Sea and adjacent basins in agreement with float measurements in the upper Southern Ocean. Teleconnections to the North Atlantic are suggested but need further studies to demonstrate their statistical significance

Observed thinning of Totten Glacier is linked to coastal polynya variability.

Analysis of ICESat-1 data (2003-2008) shows significant surface lowering of Totten Glacier, the glacier discharging the largest volume of ice in East Antarctica, and less change on nearby Moscow University Glacier. After accounting for firn compaction anomalies, the thinning appears to coincide with fast-flowing ice indicating a dynamical origin. Here, to elucidate these observations, we apply high-resolution ice-ocean modelling. Totten Ice Shelf is simulated to have higher, more variable basal melting rates. We link this variability to the volume of cold water, originating in polynyas upon sea ice formation, reaching the sub-ice-shelf cavity. Hence, we propose that the observed increased thinning of Totten Glacier is due to enhanced basal melting caused by a decrease in cold polynya water reaching its cavity. We support this hypothesis with passive microwave data of polynya extent variability. Considering the widespread changes in sea ice conditions, this mechanism could be contributing extensively to ice-shelf instability

Exploration of Antarctic Ice Sheet 100-year contribution to sea level rise and associated model uncertainties using the ISSM framework

Estimating the future evolution of the Antarctic Ice Sheet (AIS) is critical for improving future sea level rise (SLR) projections. Numerical ice sheet models are invaluable tools for bounding Antarctic vulnerability; yet, few continental-scale projections of century-scale AIS SLR contribution exist, and those that do vary by up to an order of magnitude. This is partly because model projections of future sea level are inherently uncertain and depend largely on the model's boundary conditions and climate forcing, which themselves are unknown due to the uncertainty in the projections of future anthropogenic emissions and subsequent climate response. Here, we aim to improve the understanding of how uncertainties in model forcing and boundary conditions affect ice sheet model simulations. With use of sampling techniques embedded within the Ice Sheet System Model (ISSM) framework, we assess how uncertainties in snow accumulation, ocean-induced melting, ice viscosity, basal friction, bedrock elevation, and the presence of ice shelves impact continental-scale 100-year model simulations of AIS future sea level contribution. Overall, we find that AIS sea level contribution is strongly affected by grounding line retreat, which is driven by the magnitude of ice shelf basal melt rates and by variations in bedrock topography. In addition, we find that over 1.2 m of AIS global mean sea level contribution over the next century is achievable, but not likely, as it is tenable only in response to unrealistically large melt rates and continental ice shelf collapse. Regionally, we find that under our most extreme 100-year warming experiment generalized for the entire ice sheet, the Amundsen Sea sector is the most significant source of model uncertainty (1032 mm 6σ spread) and the region with the largest potential for future sea level contribution (297 mm). In contrast, under a more plausible forcing informed regionally by literature and model sensitivity studies, the Ronne basin has a greater potential for local increases in ice shelf basal melt rates. As a result, under this more likely realization, where warm waters reach the continental shelf under the Ronne ice shelf, it is the Ronne basin, particularly the Evans and Rutford ice streams, that are the greatest contributors to potential SLR (161 mm) and to simulation uncertainty (420 mm 6spread)