Introduction

Global climate change has a great impact on ice shelves, because they are sensitive to changes in air and ocean temperature or circulation near Antarctica. The break-up of ice shelves in the Antarctic Peninsula is regarded as a signal of regional climate change (Reference Vaughan and DoakeVaughan and Doake, 1996). The area of the ice shelves around the Antarctic Peninsula has shrunk by >12 500 km² over the past three decades (Reference Scambos, Hulbe, Fahnestock, Domack, Burnett, Leventer, Conley, Kirby and BindschadlerScambos and others, 2003; Reference DomackDomack and others, 2005), resulting in the acceleration and thinning of tributary glaciers (Reference De Angelis and SkvarcaDe Angelis and Skvarca, 2003; Reference Rignot, Casassa, Gogineni, Krabill, Rivera and ThomasRignot and others, 2004; Reference Scambos, Bohlander, Shuman and SkvarcaScambos and others, 2004). Intrusion of warmer ocean water beneath the ice shelves in the Amundsen Sea sector of West Antarctica increased basal melting, resulting in ice-shelf thinning (Reference Shepherd, Wingham and RignotShepherd and others, 2004) and the accelerated discharge of inland ice mass towards the ocean (e.g. Reference RignotRignot, 1998, Reference Rignot2001; Reference Payne, Vieli, Shepherd, Wingham and RignotPayne and others, 2004; Reference Shepherd, Wingham and RignotShepherd and others, 2004; Reference ThomasThomas and others, 2004). This may lead to the instability of the West Antarctic ice sheet, which is grounded on bedrock below sea level (Reference OppenheimerOppenheimer, 1998; Reference Payne, Vieli, Shepherd, Wingham and RignotPayne and others, 2004). Understanding the dynamics of ice shelves and their interaction with ocean is therefore fundamentally important for predicting the contribution of the Antarctic ice sheet to future global sea-level rise (Reference Pritchard, Arthern, Vaughan and EdwardsPritchard and others, 2009). In this study, we examine the basal melt/freeze pattern of the Amery Ice Shelf which is fed by one of the largest glacier systems in East Antarctica.

The Amery Ice Shelf is the largest ice shelf in East Antarctica and drains the grounded ice of the Lambert Glacier/Amery Ice Shelf system. Scientific investigation of the Amery Ice Shelf began in the mid-1950s (Reference Mellor and McKinnonMellor and McKinnon, 1960). A series of surveys on the ice shelf and the ocean cavity beneath the shelf and in Prydz Bay have subsequently investigated ice-shelf dynamics and mass budget and the ice/ocean interaction, including patterns of melting and freezing at the base of the ice shelf (e.g. Reference MorganMorgan, 1972; Reference Budd, Corry and JackaBudd and others, 1982; Reference Fricker, Popov, Allison and YoungFricker and others, 2001, Reference Fricker2002; Reference AllisonAllison, 2003; Reference King, Coleman, Morgan and HurdKing and others, 2007, Reference King2009; Reference Galton-Fenzi, Maraldi, Coleman and HunterGalton-Fenzi and others, 2008; Reference Craven, Allison, Fricker and WarnerCraven and others, 2009).

Melting and freezing beneath the Amery Ice Shelf have also been studied using numerical ocean models and hydrographic observations. Reference Hellmer and JacobsHellmer and Jacobs (1992) reported a mean melting rate as high as 0.65 m a⁻¹ by modeling the sub-Amery ocean thermohaline circulation. This melt rate is equivalent to the removal of 23 G ta⁻¹ of basal ice. Using a three-dimensional (3-D) numerical ocean model, Reference Williams, Grosfeld, Warner, Gerdes and DetermannWilliams and others (2001) simulated the ocean cavity beneath the Amery Ice Shelf and obtained net melting rates of 5.8 and 18.0 Gt a⁻¹ for two different boundary conditions. Both model simulations also showed regions of basal freezing under the ice shelf, with marine ice accretion of several Gt a⁻¹. Based on hydrographic observations collected near the front of the Amery Ice Shelf, Reference Wong, Bindoff, Forbes, Jacobs and WeissWong and others (1998) reported that the net amount of ice loss from the ice-shelf bottom ranged from 10.7 to 21.9 G ta⁻¹. Using a coupled ice/ocean model, Reference HellmerHellmer (2004) estimated the spatial average basal melting rate (0.35 m a⁻¹)and total basal mass loss (17.65 G ta⁻¹) for the Amery Ice Shelf.

In the upper reaches of the ice shelf, the southern grounding line (Fig. 1) was redefined from a hydrostatic equilibrium model (Reference FrickerFricker and others, 2002) and interferometric synthetic aperture radar (InSAR) analysis (Reference RignotRignot, 2002). This grounding line extends the Amery Ice Shelf by ∼240 km upstream, compared with the grounding line position reported by Reference Budd, Corry and JackaBudd and others (1982). Using InSAR techniques, Reference Rignot and JacobsRignot and Jacobs (2002) estimated a melt rate of 31 ± 5 m ice a⁻¹ near the southern grounding line. Using remote-sensing datasets and in situ measurements, Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others (2007) reported the change in pattern of basal mass flux from the southern grounding line to the ice-shelf front along the flowbands of Lambert, Mellor and Fisher glaciers. From in situ borehole temperature and salinity measurements at AM02 (Fig. 1), Reference CravenCraven and others (2004) estimated a minimum melt rate of 0.5 m a⁻¹ at the ice-shelf base at that location. However, a complete basal melting- and freezing-rate distribution over the entire ice shelf has not been reported in these previous studies.

Fig. 1. Location of the Amery Ice Shelf, Antarctica (black box in inset). The green curves show the outline of the ice shelf and four ice rises, and the front of the ice shelf as defined by the RAMP synthetic aperture radar (SAR) image mosaic acquired in 1997 (background image: Reference JezekJezek, 1999). The dark red curve shows the southern grounding line. The light red curves show the gates where ice flux across the ice front was calculated. The orange dots are GPS stations occupied during the 1988–99 GPS campaigns (Reference KingKing, 2002, table 5-12), and the red squares are locations mentioned in the text. The color shading (scale lower left) shows the GLAS/ICESat 500 m digital elevation model (J. DiMarzio and others, http://nsidc.org/data/nsidc-0304.html).

In this paper, various remote-sensing datasets and in situ measurements have been assimilated in a Geographic Information System environment to map and analyze the basal melting/freezing patterns over the entire Amery Ice Shelf. The calving front of the Amery Ice Shelf was defined from the RADARSAT-1 Antarctic Mapping Project (RAMP) image mosaic derived in 1997 (Reference JezekJezek, 1999; Fig. 1). The grounding line of the ice shelf (Fig. 1) is based on that of Reference FrickerFricker and others (2002), and further improved by incorporating the southern grounding line position mapped by the InSAR technique (Reference RignotRignot, 2002), and the grounding line extracted by identifying the surface slope breaks between the grounded ice sheet and the floating ice shelf in Moderate Resolution Imaging Spectroradiometer (MODIS) images (personal communication from N. Young, 2009). The slope breaks are on the landward side of the grounding zone, and define the southern boundary of the ice shelf. The original polar stereographic projection was transformed to the Lambert azimuthal equal-area projection to calculate the area of the Amery Ice Shelf. The area is estimated to be 60 785 km², which is similar to the area estimate of Reference Galton-Fenzi, Maraldi, Coleman and HunterGalton-Fenzi and others (2008).

Glaciological volume conservation methods (e.g. Reference Jenkins and DoakeJenkins and Doake, 1991; Reference Joughin and PadmanJoughin and Padman, 2003) were employed and are discussed below. Based on the recent remote-sensing data, we find that net melting accounts for about 50% of the total mass loss of the ice shelf, which is far greater than a typical ratio (about 30%) of the basal melting to the total mass loss for ice shelves (Reference Jacobs, Hellmer and JenkinsJacobs and others, 1996). Basal melting and freezing under the Amery Ice Shelf play a critical role in the ice-shelf mass balance, and the interaction at the ice/ocean interface beneath the ice shelf may be more intensive than was thought previously.

Methods

The basal melting and freezing calculation is based on the assumption that the ice shelf is in a steady balance. A recent study by Reference KingKing and others (2009) suggests that the surface elevation (and hence thickness) has not changed significantly over four decades since 1968. Although short-term fluctuations of an elevation magnitude less than 1.0 m may exist, the ice shelf has been almost in balance over the longer term.

The ice mass conservation can be expressed as

where Z is the thickness of the ice shelf, V is the horizontal velocity, ∇ is the two-dimensional gradient operator, is the surface net accumulation rate and is the basal melting/freezing rate (negative for basal freezing). When the ice shelf is in steady state (∂Z/∂t = 0), the horizontal divergence of the volume flux is exactly balanced by the sum of surface accumulation and basal melting/freezing rates. The equilibrium basal melting/freezing rate can then be expressed as

For a floating ice shelf, the ice-flow speed and direction at any depth can be assumed to be the same as at the surface (Reference BuddBudd, 1966). For an ice column in steady state, the total ice inflow into the column is equal to the total outflow from mass conservation (e.g. Reference Jenkins and DoakeJenkins and Doake, 1991). Thus, for an ice column with surface area S, x-direction ice fluxes Φ _x and Φ _{x +Δx} , and y-direction ice fluxes Φ _y and Φ _{y+ Δy} , Equation (2) can be expressed as

Φ _x and Φ _{x +Δx} are calculated from the thickness Z and the x-direction component of ice velocity V. Positive values of indicate basal melting beneath the ice shelf, and negative values indicate freezing.

Three flux gates normal to the velocity vector were defined at the ice-shelf front (Fig. 1) to calculate the ice discharge and overall mass balance for the ice shelf. The ice flux Φ across each gate can be calculated as

where i (1, 2, …, n) is the index for gridcells along the flux gate, ΔX is the width of gridcells, and V_i , Z_i and ρ_i are respectively ice velocity, ice thickness and column-averaged ice density for gridcell i.

The ice thickness is derived by applying the hydrostatic equilibrium equation,

where H_i is the surface elevation for cell i, ρ _w is the sea-water density and other terms are the same as in Equation (4).

Data

The datasets used in our basal melting/freezing analysis include (1) InSAR velocity data, (2) surface elevation data, (3) column-averaged ice density, (4) ice-thickness data and (5) annual snow accumulation data.

Data Processing and Error Analysis

We first attempted to calculate the basal melting and freezing rates at a 400 m cell-size grid. But because of random errors in the velocity and ice-thickness fields, some melting/freezing rates were unrealistically variable and high (up to hundreds of metres per year). Although differences between the GPS and InSAR velocities are mostly <10 m a⁻¹, a few differences are >20 m a⁻¹. Errors in ice thickness are due to errors in the DEM and the column-averaged ice density, and while most ICESat-DEM elevations are generally within ±5 m of the GPS measurements (Fig. 4), a few elevation outliers are up to 10 m different. Therefore, we calculated the rates over larger areas, rather than the original 400 m gridcells. We defined 56 squares on the Amery Ice Shelf, each with an area of 25 km × 25 km (Fig. 5). The use of 25 km squares as the basic spatial unit for our numerical analysis increases the signal-to-noise ratio through averaging and making each square estimate more reliable. Around the edge of the ice shelf, an additional 20 smaller rectangles were defined so that the basal melting and freezing data have a relatively homogeneous distribution over the entire ice shelf (Fig. 5). In order to test the sensitivity to square size, we also calculated the basal melting rates on 20 km × 20 km squares as discussed below.

Fig. 4. Histograms and normal distribution of number of differences between the GPS and InSAR velocities (a), and between the ICESat-DEM elevations (b).

Fig. 5. Numbered areas for the calculation of the basal melting and freezing rates under the Amery Ice Shelf (see Table 3).

For each square we determined the fluxes in both the x and y directions (Equation (3)) from the averaged ice velocity and thickness values along each side of the square and converted the melting/freezing rates to water equivalent. The four ice rises were excluded from our calculation (Fig. 1).

The basal melting and freezing rates, , for the 76 defined squares and rectangles on the Amery Ice Shelf are calculated using Equation (3). The value shown in Table 3 is assigned to the geometrical center of each square or rectangle. The spatial distributions of the melting/freezing rates beneath the Amery Ice Shelf shown in Figure 6 were obtained by interpolating the 76 points (Table 3) together with 18 points along flowbands originating from Lambert, Mellor and Fisher glaciers from a previous study of Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others (2007). The additional 18 points from Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others (2007) were calculated using the AIS-DEM and the same velocity data and density model. These values have a similar uncertainty to the present study. They are used to supplement melting rates near the ice front and the southern grounding line areas. The interpolation was made using an ordinary kriging method in ArcGIS, and the interpolation result is a smooth melting/freezing rate grid with a 400 m cell size. The empirical variogram was modeled as an exponential function that was applied when spatial autocorrelation decreases exponentially with increasing distance. This function may provide the best fit through the points as the melting rates decrease rapidly downstream from the southern grounding line. The fit parameters were: number of lags =12, range = 35 000 m and nugget = 0.1 m a⁻¹.A quadrant neighborhood was defined to enclose the points used to predict values at unmeasured locations, and in each quadrant/sector of the search neighborhood a maximum of five and a minimum of two data points were used to contribute to an interpolated value.

Fig. 6. Spatial distribution of the basal melting/freezing rates (m a⁻¹) beneath the Amery Ice Shelf. Red/brown colors show the basal freezing area, and blue colors show basal melting areas. Two longitudinal bands of the maximum freezing rates are marked with A and B.

We assess the errors and their propagation in our calculation of basal melting and freezing rates. The main error sources involved in calculating the ice thickness Z are errors in the column-averaged ice density ρ_i and the ICESat-DEM elevation H. The relative error, ΔZ/Z, in ice thickness is approximated by

where elevations are typically ∼100 m, and accurate to better than 3 m with the improved geoid model. The average density of the ice column is typically 908 kg m⁻³ and has an uncertainty of ∼5 kg m⁻³ (Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others, 2007). The uncertainty of sea-water density is probably <2 kg m⁻³. Using these values, ΔZ/Z is 5.2%. The relative ice-thickness differences between the hydrostatic estimates and borehole measurements are 2.7% at AM01 (with a borehole measured thickness of 479 ± 2 m), 1.1% at AM02 (373 ± 2 m), 2.5% at AM03 (722 ± 2 m) and 6.3% at AM04 (603 ± 2 m) (Reference CravenCraven and others, 2004, Reference Craven, Allison, Fricker and Warner2009; personal communication from M. Craven, 2009). This comparison with the in situ borehole measurements suggests that the estimated ice-thickness uncertainty of 5.2% is reasonable.

The errors involved in calculating the basal melting and freezing rates include the errors in ice thickness, ice velocity components, column-averaged ice density and annual net accumulation. The error of ice velocity components is estimated to be about 10 m a⁻¹. Uncertainties of ice fluxes (Φ _x , Φ _y ) were estimated to be about 10%. Over the northern portion of the Amery Ice Shelf, the uncertainty of the basal melting/freezing rate is possibly up to 0.5 m a⁻¹, and the error of annual net accumulation could make an important contribution to the errors of basal melting and freezing rates where they are smaller than ±0.5 m a⁻¹. Over the southern portion, the uncertainty of annual net accumulation is insignificant compared to the basal melting rate, and the uncertainty of the basal melting/freezing rate is possibly up to 20%, mainly due to the errors in ice-thickness and ice-velocity components. An overall uncertainty of 15% is assumed for the total basal melting and freezing.

Although the trend of elevation change over the Amery Ice Shelf between 1968 and 2007 is not significantly different from zero, there have been two shorter periods of average elevation change exceeding ±0.1 m a⁻¹ (Reference KingKing and others, 2009). The velocity may also have slowed slightly (∼2.2 m a⁻¹ or ∼0.6%) between 1968 and 1999 (Reference King, Coleman, Morgan and HurdKing and others, 2007). Over the 4 year time-span (2000–04) between the InSAR velocity data acquisition and ICESat-DEM data acquisition, elevation and velocity changes may have respectively reached about 0.4 and 8.8 m a⁻¹ if the trends were the same as before 2000. These changes have an insignificant effect on the estimate of the basal melting and freezing, since they are systematic changes and relatively small compared with the elevation (typically 100 m) and ice velocity (typically >300 m a⁻¹)over the ice shelf.

Results and Discussion

The total area of net melting under the Amery Ice Shelf is about 34 700 km² (Fig. 6). The melting rate has a maximum value of 25.0 ± 4.0 m a⁻¹ averaged over an area of 180 km² near the southern grounding line and decreases rapidly downstream. Basal freezing primarily occurs beneath the northwest part of the ice shelf. The total area of net freezing is about 26 100 km², and the area-averaged freezing rate can be as high as 2.4 ± 0.4 m a⁻¹.The total basal melting under the entire ice shelf is 62.5 ± 9.3 Gt a⁻¹, and the total basal freezing is 16.1 ± 2.4 Gt a⁻¹, giving a total net basal mass loss of 46.4 ± 6.9 Gt a⁻¹.

When the melting/freezing rates were calculated on 20 km × 20 km squares, the total area of net melting was 33 700 km², with total basal melting of 62.2 Gt a⁻¹.The total area of net freezing was 27 100 km², with total basal freezing of 17.6 Gt a⁻¹. The total net mass loss was 44.6 Gt a⁻¹. The distribution of melting and freezing is very similar to that in Figure 6, with some small differences mostly around the edges of the ice shelf. This indicates that our results are relatively insensitive to the size of the squares over which the melting/freezing rate calculations are performed.

For the whole Lambert Glacier/Amery Ice Shelf system, which includes Antarctic drainage basins 9, 10 and 11, as defined by Reference Giovinetto and ZwallyGiovinetto and Zwally (2000), the ice flux across the grounding line onto the ice shelf is 88.9 ± 8.9 Gt a⁻¹ (Reference WenWen and others, 2008). From the modified Giovinetto compilation, the mean annual accumulation rate over the whole Amery Ice Shelf surface is 162 kg m⁻² a⁻¹,giving a total accumulation of 9.8 ± 1.0 Gt a⁻¹. The ice fluxes, Φ, were calculated for three flux gates at the ice-shelf front according to Equation (4). The total ice flux across the gates at the ice-shelf front is 47.2 ± 3.6 Gt a⁻¹. If the Amery Ice Shelf is in steady state, then total net basal melting of 51.5 ± 9.6 G ta⁻¹ is required to balance the sum of the inflow and the integrated accumulation less the mass discharge. This is consistent with the net basal melting of 46.4 ± 6.9 Gt a⁻¹ estimated above.

The freezing beneath the Amery Ice Shelf is distributed within the northwest of the ice shelf (Fig. 6), which is consistent with the marine ice distribution mapped previously by Reference Fricker, Popov, Allison and YoungFricker and others (2001) and Reference Wang, Wen, Liu, Jezek and CsathoWang and others (2006). The freezing area mapped from this research is 26 100 km², and the area of marine ice from Reference Wang, Wen, Liu, Jezek and CsathoWang and others (2006) is 26 400 km². The thickest marine ice forms two longitudinal bands, oriented along the ice-flow direction. The distribution of the maximum freezing rates (marked with A and B in Fig. 6) corresponds to the locations of the thickest marine ice, where the freezing rates range between 1.5 ± 0.2 and 2.4 ± 0.4 m a⁻¹. The freezing rates between the two longitudinal bands are relatively low, and slight melting is detected in one patch, though the melting rates are <0.5 m a⁻¹.This patch of melting is not consistent with the marine ice distribution, and may be caused by the combined errors in the ice thickness, the column-averaged ice density, annual net accumulation and velocity data, which may result in an uncertainty of the basal melting/freezing rate of up to 0.5 m a⁻¹. Based on the temperature and salinity gradients of the ocean/ice interface measured at AM02, Reference CravenCraven and others (2004) estimated a minimum melting rate of 0.5 m a⁻¹, compared to our result of 0.3 ± 0.2 m a⁻¹. Reference Craven, Allison, Fricker and WarnerCraven and others (2009) determined an average accretion rate of marine ice of 1.1 ± 0.2 m a⁻¹ between AM01 and AM04, and inferred a similar average rate upstream of 1.3 ± 0.2 m a⁻¹ between AM04 and JP (Fig. 1), with a reference density of 920 kg m⁻³. Our estimates at the corresponding locations are 0.86 ± 0.2 and 0.71 ± 0.2 m a⁻¹,respectively.

Melting of 1.0–1.6 m a⁻¹ near the ice front (Fig. 6) is caused by tidal pumping and the seasonally warmer waters of the coastal current (Reference Jacobs, Hellmer, Doake, Jenkins and FrolichJacobs and others, 1992). It varies, however, along the ice front from east to west. At the front of the thick bands of accreted marine ice, the extent and magnitude of melting is smaller or zero. The sub-freezing ice-shelf water beneath the bands possibly mitigates the melting associated with the tidal pumping (Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others, 2007).

Our estimates for the total melting, freezing and net melting beneath the Amery Ice Shelf are considerably larger than the previous results from numerical simulations (Reference Hellmer and JacobsHellmer and Jacobs, 1992; Reference Williams, Grosfeld, Warner, Gerdes and DetermannWilliams and others, 2001; Reference HellmerHellmer, 2004). These model estimates, however, all used a smaller and shallower sub-shelf cavity, based on the grounding line position of Reference Budd, Corry and JackaBudd and others (1982). The grounding line of the Amery Ice Shelf has been redefined and shifted about 240 km upstream (Reference FrickerFricker and others, 2002; Reference RignotRignot, 2002). Modeling studies of simplified ice-shelf cavities by Reference Holland, Jenkins and HollandHolland and others (2008) show the strong effect of cavity shape on basal melting rates. Reference WilliamsWilliams (1999) found that the melting and freezing pattern from a sub-iceshelf model varies significantly when the cavity is extended to that defined by Reference FrickerFricker and others (2002). The different dimensions of the sub-ice cavity as defined by the previous position of the grounding line are responsible for these smaller simulated values (Reference Galton-Fenzi, Maraldi, Coleman and HunterGalton-Fenzi and others, 2008).

Melting processes beneath ice shelves are governed by the transport of ocean heat and by the sea-water freezing temperature (Reference DoakeDoake, 1976). The ice draft at the grounding line location defined by Reference Budd, Corry and JackaBudd and others (1982) is only about 700 m, whereas at the redefined southern grounding line it is >2000 m. The in situ sea-water freezing point is lower at greater depth, and the difference between this and the ocean temperature is greater, resulting in a higher melting rate. From our present study the area-averaged melting rate is up to 25 ± 4.0 m a⁻¹ at the redefined grounding line, which is similar to the 31 ± 5 m ice a⁻¹ reported by Reference Rignot and JacobsRignot and Jacobs (2002), whereas it is only 0.5–1.0 m a⁻¹ at the location given by Reference Budd, Corry and JackaBudd and others (1982). More than 80% of the continental ice from Lambert, Mellor and Fisher glaciers is lost through basal melting beneath the southernmost 240 km of the Amery Ice Shelf (Reference Wen, Jezek, Csatho, Herzfeld, Farness and HuybrechtsWen and others, 2007; Reference Craven, Allison, Fricker and WarnerCraven and others, 2009).

Our estimate of the net melting beneath the shelf is also larger than that estimated from oceanographic measurements in Prydz Bay (Reference Wong, Bindoff, Forbes, Jacobs and WeissWong and others, 1998). A possible reason for this difference is the small number of hydro-graphic measurements across the ice-shelf front available to Reference Wong, Bindoff, Forbes, Jacobs and WeissWong and others (1998), and seasonal and interannual variability in the melting rate.

Conclusions

We have assessed the distribution of basal melting and freezing rates under the Amery Ice Shelf by integrating a variety of datasets. Our results show that the areal extents of the basal melting and freezing under the ice shelf are 34 700 and 26 100 km², respectively. The average melting rate over the melting region is ∼1.8 ± 0.3 m a⁻¹,and the maximum melt rate is ∼25.0 ± 4.0 m a⁻¹ (averaged over an area of 180 km²) under the deepest ice near the southern grounding line. The total amount of melt is 62.5 ± 9.3 Gt a⁻¹. Basal freezing is located in the northwest region of the ice shelf, and the highest freezing rate is 2.4 ± 0.4 m a⁻¹. The total basal accretion of marine ice is estimated to be 16.2 ± 2.4 Gt a⁻¹. The total net loss of mass beneath the ice shelf is 46.4 ± 6.9 Gt a⁻¹, which is consistent with the difference (51.5 ± 10.1 Gt a⁻¹) between mass input to the ice shelf (the sum of inland ice inflow plus the total surface accumulation on the ice shelf) and mass output through the ice-shelf front.

A map of the melting and freezing rates was generated using a kriging interpolation method. This map shows that the basal freezing and corresponding basal marine ice are characterized by the two longitudinal bands along the ice-flow direction. This pattern is in agreement with the marine ice distribution mapped by Reference Fricker, Popov, Allison and YoungFricker and others (2001) and Reference Wang, Wen, Liu, Jezek and CsathoWang and others (2006). The net melting accounts for about 50% of the total mass loss from the ice shelf. This is far greater than a typical ratio (about 30%) of the basal melting to the total mass loss for ice shelves (Reference Jacobs, Hellmer and JenkinsJacobs and others, 1996). Basal melting and freezing under the Amery Ice Shelf play a critical role in the ice-shelf mass balance, and the interaction at the ice/ocean interface beneath the ice shelf may be more intensive than was thought previously.