Strong Sensitivity of Pine Island Ice-Shelf Melting to Climatic Variability

See allHide authors and affiliations

Science  10 Jan 2014:
Vol. 343, Issue 6167, pp. 174-178
DOI: 10.1126/science.1244341

Cold Glacier Growth

Pine Island Glacier in Antarctica has thinned significantly during the last two decades and has provided a measurable contribution to sea-level rise as a result. Both glacier dynamics and climate are thought to be responsible for thinning, but exactly how they influence the glacier are incompletely known. Dutrieux et al. (p. 174, published online 2 January) provide another layer of detail to our understanding of the process through observations of ocean temperatures in the surrounding waters. The thermocline adjacent in the sea adjacent to the glacier calving front (where ice is discharged) lowered by 250 meters in the austral summer of 2012. This change exposed the bottom of the ice shelf to colder surface waters rather than to the warmer, deeper layer, thereby reducing heat transfer from the ocean to the overlying ice and decreasing basal melting of the ice by more than 50% compared to 2010. Those 2012 ocean conditions were partly caused by a strong La Niña event, thus illustrating how important atmospheric variability is for regulating how the Antarctic Ice Sheet responds to climate change.


Pine Island Glacier has thinned and accelerated over recent decades, significantly contributing to global sea-level rise. Increased oceanic melting of its ice shelf is thought to have triggered those changes. Observations and numerical modeling reveal large fluctuations in the ocean heat available in the adjacent bay and enhanced sensitivity of ice-shelf melting to water temperatures at intermediate depth, as a seabed ridge blocks the deepest and warmest waters from reaching the thickest ice. Oceanic melting decreased by 50% between January 2010 and 2012, with ocean conditions in 2012 partly attributable to atmospheric forcing associated with a strong La Niña event. Both atmospheric variability and local ice shelf and seabed geometry play fundamental roles in determining the response of the Antarctic Ice Sheet to climate.

Austral summer observations in the Amundsen Sea, West Antarctica, show that lightly modified, warm (0.5° to 1.2°C) and saline (>34.6) Circumpolar Deep Water (CDW), 2° to 4°C above the in situ freezing point, pervades a network of glacially scoured seabed troughs (1) (Fig. 1A). The CDW reaches nearby Antarctic glaciers and delivers heat to the base of their 200- to 1000-m-thick ice shelves (24). It is overlain by a 200- to 300-m-thick layer of cold (–1.5°C) and fresh (salinity <34.4) Winter Water (WW, Fig. 2A) that is seasonally replenished by interaction with the atmosphere and sea ice.

Fig. 1 Temperature maximum, seabed elevation, and selected observation locations in the Amundsen Sea.

(A) Maximum subsurface temperature deduced from a climatology of observations taken between 1994 and 2012 (supplementary materials). Black, gray, and white contours indicate seabed elevations (33) of –1000, –500, and –400 m, respectively. The dashed black box shows the area in (B). (B) Detailed bathymetry on the eastern Amundsen Sea continental shelf. The white line shows –750-m seabed elevation. Color-coded dots (squares) show ocean station positions used in Fig. 2 (Fig. 3, A and B). The dashed black box shows the area in (C). The PIG ice shelf has been removed to show the seabed elevation beneath (supplementary materials). (C) The 2009 coastline and ice-shelf grounding line (9) are indicated by the black line. The yellow line shows the position of the grounding line in 1996 (7). Color-coded squares indicate ocean stations used to compute PIG ice shelf melt budgets. Simulated freshwater flux into the ocean (color) for 2009 boundary conditions.

Fig. 2 Hydrographic properties in PIB.

(A) Cruise-average, 25-dbar averaged temperature (continuous lines; shading indicates 1 SD for selected years) and salinity (dashed lines) profiles from January to March 1994, 2000, 2007, 2009, 2010, and 2012 in the wider Pine Island Bay (individual profile locations are shown in Fig. 1B; number of profiles per year is indicated in parentheses). Averages are performed in density-space and the resulting profiles regridded into pressure-space using the average stratification for each year. Dashed gray lines roughly indicate the 2009 ice-front draft and the seabed ridge crest, respectively. Inset: the same data in potential temperature-salinity space, with the surface freezing line indicated in dashed blue, the 2009 CDW-glacial-ice melt line in dashed red, and the 2009 CDW-WW mixing line in dashed black. Black-dotted contours of potential density anomaly are also shown, with the isopycnals corresponding to the 2009 ice-front draft and seabed ridge crest indicated in dotted gray. (B) Mean 2012 offset from other years in potential temperature (continuous) and meltwater concentration (dashed) for all other observed summers [color-coded as in (A)]. The difference is computed in density-space and then plotted in pressure-space using the mean 2009 stratification. (C) Same as (A), but for the near-shelf edge area to the east (see Fig. 1B). (D) Same as (A), but for the near-shelf edge area to the west (see Fig. 1B).

Pine Island Glacier (PIG), a major outlet glacier feeding one such ice shelf, has shown apparently continuous thinning (5, 6) and intermittent acceleration (79) from 1973 to 2009. During this period, its ice shelf has also thinned (6, 1012), and the reduction in buttressing driven by oceanic melting is believed to be responsible for the changes inland. Earlier analysis indicated that a higher CDW volume and temperature in Pine Island Bay (PIB) in January 2009 caused an increase in ice-shelf melting and in the associated meltwater-driven circulation, relative to 1994 (2). The lack of subannual variability in CDW temperature during 1-year-long measurement in PIB (1) and the long-term correlation between the oceanic melting and the mass loss required to sustain thinning of the ice shelf gave the impression that the ice-ocean system had shown progressive change over the last two decades. This is consistent with a positive geometrical feedback, with oceanic melt enlarging the cavity under the ice shelf, allowing stronger circulation and further melting.

However, such ice-ocean systems are likely to be more complex. The glacier’s rapid change over the last few decades was probably triggered by its ungrounding from the top of a seabed ridge transverse to the ice flow at some time before the 1970s (4). Subsequent migration of the glacier’s grounding line (13) down the seabed slope upstream from the ridge crest was probably an inevitable response (14), which had a major impact on cavity geometry, but projections of future evolution remain unclear (1517). In this context, climate-driven variability of the delivery of heat to the ice-sheet margin may be important (18). To understand the relative roles of glacier dynamics and climate variability in driving past, present, and future glacier evolution requires knowledge of the spatial (19, 20) (Fig. 1C) and temporal variability of ocean-driven melting.

Ocean simulations with a varying degree of horizontal resolution and using poorly constrained seabed and ice geometry produce a broad spectrum of variability in CDW volume in PIB (21) and oceanic melt of PIG (22). Some studies suggest that basal melt is sensitive to the strength of the circulation in the Amundsen Sea (22). Others associate its variability with changing wind-driven inflow of warm CDW at the continental shelf edge and temperature variability at the calving front of the glacier (21, 23) influenced by far-field forcing of the atmospheric circulation (18). All these ocean-atmosphere sensitivities are linked dynamically, but a clear consensus on the most important driver of melt variability is lacking.

Hydrographic observations made in January 2012 in the eastern Amundsen Sea complement previous surveys from 1994, 2000, 2007, 2009 (2), and 2010 (24) (Fig. 1B) and offer a new perspective. Indeed, though near-bottom 2012 CDW temperatures of 1.2°C are very similar to the warmest records (2007, 2009, and 2010), the available heat reaching the PIG calving front is significantly reduced. The top of the thermocline (the depth at which temperature first increases sharply with increasing depth) is about 250 m deeper compared with that of any other year for which measurements exist (Figs. 2A and 3, A and B). Furthermore, the 2012 thermocline is sharpened such that the temperature increases from –0.8° to 0.8°C over only 180 m, compared with 250 to 350 m for other summers (Fig. 2A).

Fig. 3 Observed and simulated hydrography and circulation in 2009 and 2012.

(A) Section of observed and simulated 2009 potential temperatures (color) and salinity (black contours) along the eastern Amundsen Sea trough and underneath the PIG ice shelf. White lines show the surface-referenced 27.47 and 27.75 isopycnals. The panel shows observations outside the PIG cavity, and simulation results within it. Observations are linearly interpolated from profiles (black triangles) indicated in Fig. 1B. (B) Same as (A) but for the 2012 observations and simulation. (C) Modeled potential temperature (color) and velocity (black vectors; every fifth vector is shown) averaged within 50 m of the seabed for the 2009 simulation. White vectors show the corresponding velocity observed by Autosub (binned on the model grid, see also fig. S2A). The cyan line indicates the position of the section used in (A) and (B). The white line indicates 750-m seabed depth. (D) Same as (C), but for the difference between the 2012 and the 2009 simulations.

The temperature change occurring between 2010 and 2012 is largely associated with a salinity change, implying either that anomalies in the regional surface buoyancy or wind forcing cooled the upper thermocline through diapycnal exchanges (25) or that remote modifications in the atmosphere-ocean system led to decreased isopycnal advection of CDW from the continental shelf edge (21), or both. Thermocline depth changes in PIB are generally expected to occur at intraseasonal to interannual time scales (21, 22), but the few synoptic summer observations available do not allow assessment of the origin, magnitude, frequency, or duration of such variability. Hydrographic profiles taken at the continental shelf edge in the troughs leading to PIB (Fig. 2, C and D) show that the western trough has a deeper and more variable thermocline than the eastern trough, consistent with the warmest inflows being found in the east (Fig. 1A), but heavy sea ice prevented sampling in the eastern trough in 2012.

Another observation is that a temperature decrease of 0.4° to 1°C occurred in 2012 at the location of previously identified meltwater outflows (2) (around 400-m depth, Fig. 2B), and this must be associated with changes in the ice shelf–ocean interaction. Derived ice-shelf meltwater concentrations (see supplementary materials) were 30 to 50% lower at these locations than during previous warm summers, and a weaker secondary meltwater maximum developed at around 600-m depth (Fig. 2B and fig. S3). The reduced meltwater concentration is identifiable in temperature-salinity space (Fig. 2A, inset), where previous warm summer observations depart far less from a two-component mixture of CDW and meltwater. Geostrophic budgets of dissolved oxygen, temperature, and salinity at the ice front (26) show that meltwater production dropped from ~80 km3 year−1 in 2009 (2) and 2010 (24) to 37 km3 year−1 in 2012 (Table 1), assuming that synoptic sections are representative of a short-term balanced state. This amounts to a 53% decrease in meltwater production, making the basal melting of PIG in summer 2012 the lowest on record, even below that estimated in 1994, when the warmest, deepest waters were ~0.2°C cooler at the glacier front.

Table 1 Balanced-transport estimates.

For each set of observations, transports are estimated by balancing geostrophic fluxes of temperature, salinity, and dissolved oxygen budgets at the calving front of the ice shelf (see supplementary materials), allowing for a barotropic adjustment of the geostrophic velocity profiles (last column). Seawater transport into the cavity (“In”) and net transport out of the cavity (“Diff”, positive outward) are given in mSv (1 mSv = 103 m3 s−1) and in equivalent of ice (Iceq), using a water density of 1030 kg m−3 and an ice density of 920 kg m−3. Meltwater transports are also shown, using the calculated meltwater fraction and balanced geostrophic velocity. Both transports do not depart significantly from each other, providing a consistency check. The 2007 observations were taken farther away from the ice front and thus are reflective of the entire ice shelf (27), whereas other years probably represent only the fast-flowing southern part of the ice shelf (see supplementary materials). Transport estimates are typically subject to a methodological uncertainty of about 10% (2).

View this table:

To clarify links between ocean and basal melt variability, we use a state-of-the-art regional ocean–ice shelf model at 400-m resolution with the latest ice and seabed geometry (see supplementary materials). When 2009 hydrographic conditions are imposed at its lateral boundaries (Fig. 3A), the model reproduces the main oceanographic features observed under the ice shelf by an autonomous submarine (4). The ocean circulation and water properties are divided in two by the ridge (Fig. 3, A and C). Offshore of the ridge, warm CDW circulates cyclonically after entering the cavity from the north, and then meanders back toward the ridge. Only the upper part of the CDW can flow over the ridge and reach the ice-shelf grounding line. This CDW melts the ice, creating a colder, fresher, buoyant meltwater plume that rises along the ice-shelf base (Fig. 3A). The vorticity created by this process and the cavity geometry impose a vigorous cyclonic circulation on the inshore side of the ridge, with near-seabed velocity reaching 0.2 m s−1 (Fig. 3C). The boundary between the inner and outer cavities is therefore marked by a dynamical front, further emphasizing (2, 4) the crucial role played by the ridge in controlling the ocean circulation and its interaction with the ice shelf. For 2009, the model simulates reasonable basal melt patterns (Fig. 1C) compared with observations (19, 20), producing 105 km3 year−1 of meltwater over the entire ice shelf (27, 28) and 86 km3 year−1 over its more commonly quoted fast-flowing southern part (2, 19, 24) (see also supplementary materials).

In sensitivity experiments, 1994 and 2012 hydrographic conditions were imposed at the lateral boundaries while all other features (including ice geometry) were kept the same. For 2012, the simulated cooling of the sub-ice cavity and diminished basal melting are consistent with the observed decrease in meltwater production. The lowered thermocline in 2012 reduces the amount of heat flowing over the ridge and cools the water that reaches the grounding line by 0.3°C (Fig. 3B). In turn, melting is less pronounced and the buoyancy-driven circulation is concurrently decreased, particularly in the inner cavity (Fig. 3D). Overall, the model estimates a 31 to 38% decrease in meltwater flux from the ice shelf in response to the cooling at its boundaries (table S1). For 1994, when deepest waters are cooler than in 2012 but mid-depth waters are warmer, the model estimates only an 8 to 10% decrease in melting. These simulations neglect surface forcing and are limited by the imperfect seabed and constant ice geometry, which may explain the smaller reduction in melting than that observed. However, the model confirms the large sensitivity of the present ice-ocean system to the depth of the thermocline relative to the ridge that is implied by the observations. The inability of the simulation to proportionally reproduce the 1994 melting level also suggests the importance of unresolved changes in cavity geometry between 1994 and 2009 (13).

If prolonged, the ocean conditions observed in January 2012 would have profound implications for the PIG ice shelf. Continuation of a deep thermocline would reverse the current ice-shelf thinning (as ice advection overcompensates the weakened basal melt), potentially permitting a readvance of the grounding line. However, conditions in the months leading up to the 2012 observations were unusual (see also supplementary materials). Most of 2011 was marked by strong easterly wind anomalies over the Amundsen Sea (Fig. 4, A and B), weakening the typically cyclonic wind stress curl over the continental shelf (Fig. 4B) and curtailing the westerlies north of the shelf edge that are thought to enhance onshore CDW transport (21). Integrated over the preceding year, the zonal wind north of the continental shelf even reversed to easterly in 2011 (red line in Fig. 4B), an occurrence that is unique in the reanalysis wind record dating back to 1979. Such conditions potentially decreased the flux of CDW onto the continental shelf and enhanced coastal downwelling, creating the 2012 thermohaline configuration with a relatively thin CDW layer in PIB.

Fig. 4 Exceptional forcing before the 2012 anomaly.

(A) Monthly-mean 10-m-height zonal wind velocity (34) (black) and its running integral over the preceding 12 months (red), averaging over the area indicated by the black-dotted box in (B). Shaded areas indicate 1 SD around the mean for each. The green line shows the sea-surface temperature monthly anomaly in the Niño 3.4 area. Dates of ocean observations are indicated by colored vertical lines. (B) 2011 annual mean of anomalies in surface wind (vector) and wind-stress curl (color) with respect to the 1979–2011 monthly climatology of each. (C) 2011 annual mean of anomalies in atmospheric geopotential height at 200 hPa (color) and sea surface zonal wind (contour) with respect to the 1979–2011 monthly climatology. White (black) contour are negative (positive), with 1 m s−1 interval. The white and magenta boxes indicate the area depicted in (B) and the Niño 3.4 area, respectively.

The wind variability over the Amundsen Sea has both local and remote origins. In particular, convective anomalies in the equatorial Pacific troposphere are known to radiate an atmospheric wave train toward the area (29), affecting winds at the continental shelf edge. The strong easterly anomaly in January 2012 can be attributed to exceptionally weak convection in the western tropical Pacific associated with a major La Niña event developing in late 2011 (Fig. 4C). The relevance of remote forcing from the tropical Pacific is illustrated by the statistically significant correlation (0.48) between the Niño 3.4 sea-surface temperature (SST) anomaly and the zonal wind over the continental shelf edge (Fig. 4A, see also supplementary materials). Although other factors affect the thermohaline structure in PIB, such as eddy variability at the shelf edge (30), coastal downwelling (1), and buoyancy forcing (25, 26, 31), the conjunction of the January 2012 ocean conditions and the atmospheric anomalies in the preceding months confirms the important role of tropical forcing in determining the variability of melt rates under the PIG ice shelf (18).

Neighboring ice shelves closer to the continental shelf edge and exposed to greater ocean variability are naturally prone to a larger variability in melting (32). However, the presence of the ridge in front of the grounding line of PIG enhances its sensitivity to changes in oceanic and climatic forcing. Specifically, the blocking effect of the ridge makes the vertical distribution of heat a key element of the melting response. This study therefore stresses the importance of both local geometry and climate variability in determining ice-shelf melting. With the added realization that this melting varies appreciably over kilometer scales (19, 20) and that oceanic variability is expected on intraseasonal to interannual time scales, it is clear that progress in the understanding and prediction of ice-sheet contributions to sea-level rise requires observations and models that capture a wide range of spatial and temporal scales.

Supplementary Materials

Materials and Methods

Figs. S1 to S10

Tables S1

References (3560)

References and Notes

  1. Acknowledgments: P.D. and J.D.R. were supported by Natural Environment Research Council grants NE/G001367/1, NE/H02333X/1, and NE/J005770/1. H.K.H. and S.H.L. were supported by Korea Polar Research Institute grants PP12010 and PP13020. The European Centre for Medium-Range Weather Forecasts is acknowledged for serving the ERA-Interim reanalysis data set, and Center for Remote Sensing of Ice Sheets (with support from NSF grant ANT-0424589 and NASA grant NNX10AT68G) is acknowledged for generating and serving radar observations.
View Abstract

Navigate This Article