Inflow of Warm Circumpolar Deep Water in the Central Amundsen Shelf*

JUNE 2010 W Å HLIN ET AL. 1427 Inflow of Warm Circumpolar Deep Water in the Central Amundsen Shelf* A. K. WÅHLIN Department of Earth Sciences, University of Gothenburg, Göteborg, Sweden X. YUAN Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York G. BJÖRK AND C. NOHR Department of Earth Sciences, University of Gothenburg, Göteborg, Sweden (Manuscript received 14 January 2010, in final form 25 February 2010) ABSTRACT The thinning and acceleration of the West Antarctic Ice Sheet has been attributed to basal melting induced by intrusions of relatively warm salty water across the continental shelf. A hydrographic section including lowered acoustic Doppler current profiler measurements showing such an inflow in the channel leading to the Getz and Dotson Ice Shelves is presented here. The flow rate was 0.3 0.4 Sv (1 Sv [ 10 6 m 3 s 21 ), and the subsurface heat loss was estimated to be 1.2 1.6 TW. Assuming that the inflow persists throughout the year, it corresponds to an ice melt of 110 130 km 3 yr 21, which exceeds recent estimates of the net ice glacier ice volume loss in the Amundsen Sea. The results also show a 100 150-m-thick intermediate water mass consisting of Circumpolar Deep Water that has been modified (cooled and freshened) by subsurface melting of ice shelves and/or icebergs. This water mass has not previously been reported in the region, possibly because of the paucity of historical data. 1. Introduction The glaciers that drain the West Antarctic Ice Sheet into the Amundsen Sea are accelerating, presumably because of melting of the oceanic ice shelves (e.g., Pritchard et al. 2009; Shepherd et al. 2004; Thomas 2004; Jacobs et al. 1996; Rignot and Thomas 2002; Walker et al. 2007). The suspect heat source for such glacial melting (Walker et al. 2007; Thoma et al. 2008) is the relatively warm Circumpolar Deep Water (CDW) found in the adjacent oceanic domain, with temperature above 1.88C near the southern boundary of the Antarctic Circumpolar Current (Orsi et al. 1995). Farther to the south, the Lower CDW available to the Amundsen Sea has neutral density (Jackett and McDougall 1997) between 28.03 and * Lamont-Doherty Earth Observatory Publication Number 7337. Corresponding author address: Dr. A. K. Wåhlin, Department of Earth Sciences, University of Gothenburg, PB 460, SWE-405 30 Göteborg, Sweden. E-mail: awahlin@gu.se 28.27 kg m 23 (Whitworth et al. 1998). Even after mixing with colder surface water as it enters the continental shelf, Modified CDW (MCDW) is still warmer than 18C in the Amundsen Sea. MCDW with temperatures approximately 38C above freezing has been observed to penetrate across the east Amundsen shelf into the ice shelves through a deep channel there (Walker et al. 2007). The central Amundsen shelf is separated from the east Amundsen shelf by a 400 500-m-shallow bank (Fig. 1), which extends all the way to Bear Peninsula, where it also divides the deeper ice shelf basins into a western part and an eastern part (Nitsche et al. 2007). The circulation of MCDW in the central Amundsen shelf is not known. The only existing estimate of MCDW inflow to the Amundsen shelf was done in the eastern channel (Fig. 1) leading to the east Amundsen shelf (Walker et al. 2007). Geostrophic calculations using a hydrographic transect across the channel showed an inflow of 0.25 Sv (1 Sv [ 10 6 m 3 s 21 ). Here, new hydrographic data, complemented with lowered acoustic Doppler current profiler (LADCP) velocity measurements, across the western channel leading into the central Amundsen shelf are presented for the first time. DOI: 10.1175/2010JPO4431.1 Ó 2010 American Meteorological Society

1428 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 40 FIG. 1. Map of the region with the CTD stations. Purple dots show the cross section from 2008 expedition, red dots show the drifting yo-yo stations, and pink and green dots show the historical data. The yellow dashed line shows the ice edge on 20 Dec 2008. The arrow shows the on-shelf direction (408 deviation from the compass north). 2. Methods During December 2008, a survey was conducted in the central Amundsen shelf, including one section across the western channel and three yo-yo stations (Fig. 1). Hydrographic data were sampled with the conductivity temperature depth (CTD) system from Sea-Bird Electronics, Inc., (SBE) the SBE 9111 CTD with double sensor pack: SBE 31 temperature, SBE 4C conductivity, and SBE 43 oxygen sensors. The sensors were calibrated before and after the cruise, and the data were compensated for the drift in salinity, temperature, and oxygen. Current shear was measured using a Teledyne RD Instruments, Inc., 300-kHz Workhorse Sentinel LADCP attached to the rosette sampler. The LADCP had bottom track everywhere where MCDW was detected. The velocity data were detided using the barotropic tide model (Circum-Antarctic Tidal Simulation Model CATS2008b) described in Padman et al. (2002), using all 10 available tidal components (further details of the model can be found online at http://www.esr.org/;padman/antarctic_ Tides/cats0201.html). Historical data from the World Ocean Database (Boyer et al. 2006) have also been included. All stations located on the shelf between 718 and 758S and between 1308 and 1008W that were available in the database were used. The data collected by Walker et al. 2007 were not present in the database.

JUNE 2010 W Å HLIN ET AL. 1429 FIG. 2. Sections across the western channel looking to the northwest (see Fig. 1 for station positions): (a) temperature (8C) (b) salinity (psu) (c) dissolved oxygen (ml L 21 ), and (d) along-channel velocity component (m s 21 ). Negative velocities indicate flow toward the ice shelves (see Fig. 1 for exact orientation). Thin white lines are constant neutral density (kg m 23 ), and thick white line shows the boundary of the MCDW (i.e., neutral density 5 28.03 kg m 23 ). 3. Results Figure 2 shows across-channel sections (note that the easternmost part, stations 33 and 34, were occupied 3 days before the rest of the section) of temperature T, salinity S, dissolved oxygen O 2, and the detided along-channel velocity component. Thin white lines show contours of neutral density g N (Jackett and McDougall 1997), and the thick white line shows g N 5 28.03kgm 23,which marks the boundary between MCDW and overlying water masses (Whitworth et al. 1998). Three main water masses can be seen in the acrosschannel sections. Extending from the surface down to 200 400-m depth is cold Antarctic Surface Water [oftenalsoreferredtoaswinterwater(ww);jenkins 1999], which is defined as having g N, 8.03 kg m 23 and potential temperature u,21.68c (Bindoff et al. 2000). It is characterized by oxygen values higher than 5mLL 21 and has blue color in the temperature section and green to blue colors in the salinity section. At the bottom, in the northeastern part of the channel, there is a layer of comparatively warm and salty MCDW with g N. 28.03 kg m 23, T 18C, S 34.7 psu, and O 2 4.3 ml L 21. Below the WW and above the MCDW (where the MCDW layer is present), there is an approximately 150-m-thick layer of water with 27.75, g N, 28.03 kg m 23, T 08C, S 34.3 psu, and oxygen values below 4.8 ml L 21. The warmest and saltiest water is concentrated along the eastern flank of the channel and shoaling toward the northeast end, giving rise to a significant upward slope of the neutral density surfaces at depths below 250 m. The slope is consistent with a southward geostrophic flow along the channel as evident in the LADCP data. The along-channel velocities from LADCP measurements (Fig. 2d) show generally currents in the southerly direction from the bottom to 100 150 m above bottom. The southward current is strongest, with speeds up to 10 cm s 21 at stations 56 and 57, at the deepest part of the channel, and at station 63, farther up the slope. Five

1430 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 40 obtained from the (detided) LADCP measurements in the surface. Figure 3 shows the sum of the barotropic and the baroclinic part. Both LADCP and geostrophic calculations indicate that approximately 0.3 0.4 Sv of MCDW entered onto the central Amundsen shelf through the western channel, which can be compared to the geostrophic transport of 0.25 Sv measured by Walker et al. (2007) for the same water mass in the eastern channel. Figure 4 shows T S diagrams as well as oxygen concentration of the data from the present cruise and historical data from the World Ocean Database (Boyer et al. 2006). The historical data have been divided into those taken in the vicinity of ice shelves (green dots in Fig. 1) and others (pink dots in Fig. 1). Also shown in Fig. 4 is the surface freezing temperature (red line). The solid black line is the Gade line discussed below. 4. Meltwater mixture FIG. 3. (top) Vertical mean of the along-channel velocity component (cm s 21 ) measured at the yo-yo stations (see Fig. 1 for positions) along with the barotropic tide from the tidal model (Padman et al. 2002). (bottom) Cumulative along-channel transport (Sv) of MCDW (i.e., water with neutral density exceeding 28.03 kg m 23 ) integrated across the channel. Negative values indicate flow toward the ice shelves. The dashed line is detided LADCP data, and the solid line is the geostrophic transport. Error bars show the cumulative integrated standard error from the LADCP. of the stations, 58 62, were very closely spaced, and only the temporal average of the stations is plotted in Fig. 2. During the cruise, three time series yo-yo stations were occupied (Fig. 1), during which the ship drifted with the ice and CTD casts were performed every 2 h. A comparison between the velocity measured by the LADCP (the vertical average in the water column) at the yo-yo stations and the on-shelf tidal current component, calculated by the tide model in Padman et al. (2002), is shown in Fig. 3. The tidal currents are generally less than 3 cm s 21, and the observed velocities are consistent with the tide model results. The cumulative along-channel transport of MCDW (i.e., water with g N. 28.03 kg m 23 bounded by the thick white line in Figs. 2a c) is shown in the bottom panel of Fig. 3. The dashed line is the (detided) LADCP data, and the solid line is the geostrophic transport. The geostrophic transport was referenced to the surface; that is, the baroclinic pressure gradient obtained from the CTD section was used to calculate the baroclinic part of the geostrophic velocity, and the barotropic part was When relatively warm ocean water comes into subsurface contact with ice, the ice will warm up until it begins to melt. Unless the ice is freshly calved and colder than approximately 2308C, the energy for warming the ice to freezing temperature (i.e., the sensible heat transfer from ocean to ice) is small relative to the energy for the melting of the ice (i.e., the latent heat transfer from ocean to ice). Assuming that the sensible heat transfer is negligible relative to the latent heat transfer and that the volume of meltwater is small relative to the volume of ocean water (and hence the energy for warming the meltwater is also negligible relative to the latent heat transfer), then the ocean meltwater mixture that results from such a process will have temperature T P and salinity S P given by (e.g., Gade 1979; Jenkins 1999) T P (S P ) 5 T OCEAN 1 L F 1 S OCEAN, (1) C P S P where T OCEAN and S OCEAN are, respectively, temperature and salinity of the ocean water prior to the melting; L F 5 334 kj kg 21 is the latent heat of fusion for ice; and C P 5 3.97 kj kg 21 K 21 is the specific heat of water with salinity of 34.7 psu, temperature of 18C, and pressure of 400 dbar. Equation (1) is hereafter referred to as the Gade line. It has been plotted in Fig. 4 (black solid line) using the characteristic T and S for MCDW: that is, T OCEAN 18C and S OCEAN 34.7 psu. Water that is a mixture between the MCDW and the ice shelf meltwater, referred to here as MCDW meltwater mixture, will fall close to the Gade line. In Jenkins (1999), the hydrographic properties of water in the east Amundsen shelf and the effect of subsurface glacier melting (from the Pine Island Glacier) on the water

JUNE 2010 W Å HLIN ET AL. 1431 FIG. 4. Salinity (psu) temperature (8C) plot of the CTD data with colors indicating dissolved oxygen (ml L 21 ) according to color bar. The red line is the surface freezing temperature, and the black line shows the Gade line [see Eq. (1)]. Thin black lines show g N (kg m 23 ; Jackett and McDougall 1997) according to labels, and the two green lines show g N 5 28.03 kg m 23 and g N 5 28.27 kg m 23, which bound the MCDW. (a) Data from the 2008 expedition (color) and historical data (gray). Light gray indicates stations away from the ice shelves marked with pink dots in Fig. 1, and dark gray indicates stations close to ice shelves marked with green dots in Fig. 1. (b) Data from the 2008 expedition only. (c) Historical data from stations away from the ice shelves marked with pink dots in Fig. 1. (d) Historical data from stations close to the ice shelves marked with green dots in Fig. 1. Gray markers in (c) and (d) indicate that no oxygen data were available. masses were examined. Based on stations in the east Amundsen shelf, it was shown that the water on the shelf that is not affected by subsurface ice melt is a mixture of MCDW (with T 18C, S 34.7 psu, O 2 4.3 ml L 21 ) and WW (with T 21.78C, S 34 psu, O 2 7mLL 21 ) that falls on a straight line in the S O 2 plane and in the T O 2 plane. Subsurface melting of the Pine Island Glacier lowers the salinity and the temperature but does not increase the oxygen content as much as mixing with WW does, and hence the MCDW meltwater mixture deviates from the mixing line between MCDW and WW. In Fig. 5, the mixing line between MCDW and WW (Jenkins 1999; Jenkins et al. 1997) has been plotted in the S O 2 plane and in the T O 2 plane. Also shown are the present dataset (circles) and the historical data (plus signs). The color denotes R, the deviation from the Gade line in degrees Celsius, given by qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi R 5 [T i T P (S i )] 2, (2) where T i and S i are, respectively, the measured temperature and salinity and T P (S i ) is the temperature calculated from Eq. (1). The present data, unlike the historical data, show a strong signal of MCDW meltwater mixture in the T S, T O 2, and S O 2 planes. In the T S plane, the strong signal is indicated by the fact that the data points follow

1432 JOURNAL OF PHYSICAL OCEANOGRAPHY FIG. 5. Historical data (plus signs) and the 2008 data (circles) plotted in the (a) T (8C) O2 (ml L21) and (b) S (psu) O2 (ml L21) planes. The black dashed lines show the mixing line between WW and MCDW. (c) Crosschannel distribution of R (8C). In (a) (c), colors indicate deviation (8C) from the Gade line according to the color bar. VOLUME 40

JUNE 2010 W Å HLIN ET AL. 1433 the Gade line and have low oxygen values. In the T O 2 and S O 2 planes, it is indicated by the fact that the data have lower oxygen values than the mixing line between MCDW and WW and that they have low R values. The only historical stations that show the presence of the MCDW meltwater mixture are three stations very close to Pine Island Glacier (Jenkins 1999; Jacobs et al. 1996). Data from these three stations are marked in Fig. 5. To determine the spatial distribution of MCDW meltwater mixture, an across-channel section plot of R [cf. Eq. (2)] has been constructed (Fig. 5, bottom). The MCDW meltwater mixture is clearly identifiable and occupies a 100 150-m-thick layer between the incoming MCDW (bounded by the thick white line) and the WW (yellow and red colors). The MCDW meltwater mixture is relatively warm (T 08C), salty (S 34.3 psu), and oxygen poor (O 2, 4.8 ml L 21 ) relative to the WW but colder and fresher than the MCDW. Previous observations of MCDW meltwater mixture (e.g., in the east Amundsen shelf; Jacobs et al. 1996; Jenkins 1999) show thin intrusions of MCDW meltwater mixture interleaving into the ambient water. A distinct water layer consisting of MCDW meltwater mixture has not been previously reported in the region. The MCDW meltwater mixture is approximately 18C cooler than the MCDW. The heat transferred to the ice shelves or icebergs through subsurface melting is hence estimated to be 1.2 1.6 TW (using the measured transport of 0.3 0.4 Sv), equivalent to melting 110 140 km 3 ice annually if the melting persists throughout the year. This value can be compared, for example, to the recent estimate of the dynamic thinning of the glaciers feeding the Dotson and Crosson Ice Shelves (Pritchard et al. 2009). Using a satellite-borne laser altimeter, the thinning in 2003 07 in the most rapidly moving parts was estimated to exceed 9 m yr 21, which gives a net volume loss there of approximately 50 km 3 yr 21. Regardless of whether the subsurface ice melt occurs underneath the ice shelves or by melting of icebergs, we expect the mean flow of the MCDW meltwater mixture to be directed out from the shelf region. Because the cross-sectional area occupied by this water mass is larger than the one occupied by MCDW, the average velocity in the layer will be smaller. Using a cross-sectional area of approximately 25 km 2 and the measured inflow of 0.3 0.4 Sv gives an average outward velocity in the MCDW meltwater mixture layer of less than 1.5 cm s 21 probably too slow to detect. 5. Discussion The present results show the existence of a large body of water produced by MCDW that has been cooled and freshened in a subsurface melting process of ice shelves and/or icebergs. This MCDW meltwater mixture occupies a 100 150-m-thick intermediate layer between the MCDW and the WW. Previous observations of MCDW meltwater mixture on the shelf (e.g., Jenkins 1999; Jacobs et al. 1996) has been in close proximity to the ice shelves (e.g., Jenkins and Jacobs 2008; Potter et al. 1984, 1988) where it has formed thin intrusions into the ambient water. However, only approximately 6 10 of the historical stations from the Amundsen shelf were located in the deep channels, and most of these stations were in the eastern channel (Fig. 1). Hence, it is possible that the MCDW meltwater mixture has been present on the shelf previously but not detected. In March 2003 (Walker et al. 2007), the entrance of the eastern channel was transected with closely spaced CTD stations that show the presence and inflow of MCDW onto the shelf. However, the layer of WW extended all the way down to the MCDW; if there was any MCDW meltwater mixture present, it occupied only a very thin layer. The MCDW meltwater mixture may exit the shelf at another location: for example, flowing westward along the coast first, steered by the earth s rotation. Planned future field studies will attempt to better resolve the spatial distribution and shorter-term time changes of the meltwaterenriched deep flow. Acknowledgments. The cruise Oden Southern Ocean 2008/09 was carried out on IB Oden, and we are grateful for all of the help and logistic support from the crew. Bengt Liljebladh, Kevin Pedigo, and Jeremy Lucke provided expert technical assistance with preparations as well as during the cruise. The cruise was supported by National Science Foundation (NSF), Swedish Polar Research Secretariat, and the Swedish Research Council (SRC). The data from the World Ocean Database were compiled by the National Oceanographic Data Center of NOAA. Wåhlin, Nohr, and Björk were funded by SRC, and Yuan was supported by NSF Grant ANT08-36137. REFERENCES Bindoff, N. L., M. A. Rosenberg, and M. J. Warner, 2000: On the circulation and water masses over the Antarctic continental slope and rise between 808 and 1508 E. Deep-Sea Res. II, 47, 2299 2326. Boyer, T. P., and Coauthors, 2006: World Ocean Database 2005. NOAA Atlas NESDIS 60, 190 pp. Gade, H., 1979: Melting of ice in sea water: A primitive model with application to the Antarctic ice shelf and icebergs. J. Phys. Oceanogr., 9, 189 198. Jackett, D. R., and T. R. McDougall, 1997: A neutral density variable for the world s oceans. J. Phys. Oceanogr., 27, 237 263. Jacobs, S., H. Hellmer, and A. Jenkins, 1996: Antarctic ice sheet melting in the southeast Pacific. Geophys. Res. Lett., 23, 957 960.

1434 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 40 Jenkins, A., 1999: The impact of melting ice on ocean waters. J. Phys. Oceanogr., 29, 2370 2381., and S. Jacobs, 2008: Circulation and melting beneath George VI Ice Shelf, Antarctica. J. Geophys. Res., 113, C04013, doi:10.1029/ 2007JC004449., D. G. Vaughan, S. S. Jacobs, H. H. Hellmer, and J. R. Keys, 1997: Glaciological and oceanographic evidence of high melt rates beneath Pine Island Glacier, West Antarctica. J. Glaciol., 43, 114 121. Nitsche, F. O., S. S. Jacobs, R. D. Larter, and K. Gohl, 2007: Bathymetry of the Amundsen Sea continental shelf: Implications for geology, oceanography, and glaciology. Geochem. Geophys. Geosyst., 8, Q10009, doi:10.1029/2007gc001694. Orsi, A., T. Whitworth, and W. Nowlin, 1995: On the meridional extent and fronts of the Antarctic Circumpolar Current. Deep- Sea Res. I, 42, 641 673. Padman, L., H. A. Fricker, R. Coleman, S. Howard, and S. Erofeeva, 2002: A new tidal model for the Antarctic ice shelves and seas. Ann. Glaciol., 34, 247 254. Potter, J. R., J. G. Paren, and J. Loynes, 1984: Glaciological and oceanographic calculations of the mass balance and oxygen isotope ratio of a melting ice shelf. J. Glaciol., 30, 161 170., M. H. Talbot, and J. G. Paren, 1988: Oceanic regimes at the ice fronts of George VI Sound, Antarctic Peninsula. Cont. Shelf Res., 8, 347 362, doi:10.1016/0278-4343(88)90008-8. Pritchard, H. D., R. J. Arthern, D. G. Vaughan, and L. A. Edwards, 2009: Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets. Nature, 461, 971 975, doi:10.1038/nature08471. Rignot, E., and R. Thomas, 2002: Mass balance of polar ice sheets. Science, 297, 1502 1506. Shepherd, A., D. J. Wingham, J. A. D. Mansley, and H. F. J. Corr, 2004: Inland thinning of Pine Island Glacier, West Antarctica. Science, 291, 862 864. Thoma, M., A. Jenkins, D. Holland, and S. Jacobs, 2008: Modelling Circumpolar Deep Water intrusions on the Amundsen Sea continental shelf, Antarctica. Geophys. Res. Lett., 35, L18602, doi:10.1029/2008gl034939. Thomas, R. H., 2004: Force-perturbation analysis of recent thinning and acceleration of Jakobshavn Isbræ, Greenland. J. Glaciol., 50, 57 66. Walker, D. P., M. A. Brandon, A. Jenkins, J. T. Allen, J. A. Dowdeswell, and J. Evans, 2007: Oceanic heat transport onto the Amundsen Sea shelf through a submarine glacial trough. Geophys. Res. Lett., 34, L02602, doi:10.1029/2006gl028154. Whitworth, T., III, A. H. Orsi, S.-J. Kim, and W. D. Nowlin Jr., 1998: Water masses and mixing near the Antarctic Slope Front. Ocean, Ice, and Atmosphere: Interactions at the Antarctic Continental Margin, S. S. Jacobs and R. F. Weiss, Eds., Antarctic Research Series, Vol. 75, Amer. Geophys. Union, 1 27.