Sentinel-3 Delay-Doppler altimetry over Antarctica
- 1UK Centre for Polar Observation and Modelling, Centre of Excellence in Environmental Data Science, Lancaster University, Lancaster, LA1 4YW, UK
- 2Department of Space and Climate Physics, Mullard Space Science Laboratory, University College London, Gower Street, London, WC1E 6BT, UK
- 3Centre for Polar Observation and Modelling, University of Leeds, Leeds, LS2 9JT, UK
- 4IsardSAT Ltd, Surrey Space Incubator, 40 Occam Road, The Surrey Research Park, Guildford, Surrey, GU2 7YG, UK
- 5CLS, 11 Rue Hermes, Parc Technologique du Canal, 31520 Ramonville Saint-Agne, France
- 6SERCO, c/o ESA ESRIN, Largo Galileo Galilei, 1, 00044 Frascati RM, Italy
- 7DEIMOS, c/o ESA ESRIN, Largo Galileo Galilei, 1, 00044 Frascati RM, Italy
- 8ESA ESRIN, Largo Galileo Galilei, 1, 00044 Frascati RM, Italy
Correspondence: Malcolm McMillan (firstname.lastname@example.org)
The launch of Sentinel-3A in February 2016 represented the beginning of a new long-term series of operational satellite radar altimeters, which will provide Delay-Doppler altimetry measurements over ice sheets for decades to come. Given the potential benefits that these satellites can offer to a range of glaciological applications, it is important to establish their capacity to monitor ice sheet elevation and elevation change. Here, we present the first analysis of Sentinel-3 Delay-Doppler altimetry over the Antarctic ice sheet, and assess the accuracy and precision of retrievals of ice sheet elevation across a range of topographic regimes. Over the low-slope regions of the ice sheet interior, we find that the instrument achieves both an accuracy and a precision of the order of 10 cm, with ∼98 % of the data validated being within 50 cm of co-located airborne measurements. Across the steeper and more complex topography of the ice sheet margin, the accuracy decreases, although analysis at two coastal sites with densely surveyed airborne campaigns shows that ∼60 %–85 % of validated data are still within 1 m of co-located airborne elevation measurements. We then explore the utility of the Sentinel-3A Delay-Doppler altimeter for mapping ice sheet elevation change. We show that with only 2 years of available data, it is possible to resolve known signals of ice dynamic imbalance and to detect evidence of subglacial lake drainage activity. Our analysis demonstrates a new, long-term source of measurements of ice sheet elevation and elevation change, and the early potential of this operational system for monitoring ice sheet imbalance for decades to come.
Accurate knowledge of ice sheet topography and regional changes in ice volume is essential for developing a process-based understanding of ice sheet evolution and for monitoring the response of Greenland and Antarctica to climate change (Davis et al., 2005; Price et al., 2011; Shepherd et al., 2004). For the past quarter of a century, satellite radar altimeters have provided near-continuous coverage of Earth's polar regions, yielding detailed topographic information of ice sheets (Bamber et al., 2009; Bamber and Bindschadler, 1997; Helm et al., 2014; Remy et al., 1989; Slater et al., 2018), together with estimates of changes in ice sheet volume (Davis and Ferguson, 2004; Helm et al., 2014; Johannessen et al., 2005) and mass (McMillan et al., 2014, 2016; Shepherd et al., 2012; Wingham et al., 2006b; Zwally et al., 2011). By resolving changes at the scale of individual glacier basins, these satellites have been able to identify emerging signals of imbalance (Flament and Rémy, 2012; Wingham et al., 2009), loci of rapid ice loss (Hurkmans et al., 2014; McMillan et al., 2016; Sørensen et al., 2015; Zwally et al., 2005), and the regional contribution of ice sheets to global sea level rise (Shepherd et al., 2012; The IMBIE team, 2018; Zwally et al., 2011).
During the earlier part of this 25-year record, missions carried conventional low-resolution, or pulse-limited, instruments, including those flown on board the ERS-1, ERS-2, Envisat, and SARAL satellites. These systems, which were originally developed to measure the ocean geoid, flew to a latitude of ∼81∘ and provided a ground footprint of approximately 2 km2 (corresponding to a Ku-band pulse-limited footprint over a flat, orthogonal surface, a 320 MHz measured bandwidth, and a satellite altitude of 700–800 km). The size of this footprint, together with the large area illuminated by the radar antenna beam (∼200 km2), meant that correctly locating the origin of the surface reflection in regions of complex terrain could be challenging. In 2010, the first dedicated ice radar altimetry mission, CryoSat-2, was launched, with two improvements in system design that were specifically aimed at enhancing altimeter performance in areas of steep and complex ice margin terrain. Synthetic aperture radar (SAR), or Delay-Doppler, processing delivered a four-fold improvement in along-track resolution, to approximately 400 m, and interferometric techniques were used to locate the origin of the surface reflection in the across-track plane (Raney, 1998; Wingham et al., 2006a). These developments, in conjunction with the unique long-period, high-inclination orbit, have delivered improved coverage of Earth's ice sheets (McMillan et al., 2014, 2017) and yielded greater confidence in determining their ongoing evolution.
With CryoSat-2 now operating far beyond its original design lifetime of 3.5 years, and the recent initiation of the Copernicus operational programme of satellites, there is presently a need to establish the utility of the new class of Sentinel-3 Delay-Doppler radar altimeters (Donlon et al., 2014) for the purpose of monitoring ice sheet change. The first of four satellites, Sentinel-3A, was launched on 16 February 2016 and was followed by Sentinel-3B on 25 April 2018. Each satellite provides coverage up to a latitude of 81.35∘, with 385 orbital revolutions per cycle, yielding an on-the-ground revisit time of 27 days (Table 1), together with a 4-day sub-cycle. The main altimetry payload, SRAL, is a Ku-band SAR altimeter, which provides elevation measurements with a resolution of ∼300 m along track by ∼1.6–3 km across track, depending upon the surface roughness (Chelton et al., 1989). To date, the focus of Sentinel-3A exploitation has been upon retrievals over ocean and inland water surfaces, with early studies demonstrating its capability to retrieve fine-scale (∼20 km) oceanographic features (Heslop et al., 2017), to increase the quality of river level and discharge estimates in central Africa (Bogning et al., 2018), and to resolve near-coastal sea surface heights (Bonnefond et al., 2018).
Over ice sheet surfaces, Sentinel-3 is unique among altimeters because it operates in Delay-Doppler mode across all regions. This mode of operation contrasts with CryoSat-2, which operates in low-resolution mode over the interior of each ice sheet and SAR interferometric mode at coastal locations. As a consequence, although no interferometric information is available to aid Sentinel-3 retrievals around the ice sheet margins, high-resolution measurements are for the first time routinely acquired throughout the ice sheet interior. Given these different operating modes across both inland and coastal ice sheet regions, together with the future longevity of the EU Copernicus programme of operational satellites, it is imperative that early assessments of the accuracy and precision of the instrument are made over ice sheet surfaces, to establish the basis for glaciological applications of these data. Here we provide a first evaluation of Sentinel-3 Delay-Doppler measurements over Antarctica, to determine its utility for monitoring ice sheet surfaces.
To evaluate the performance of the Sentinel-3 Delay-Doppler altimeter across a range of topographic regimes, we selected four study sites across Antarctica for detailed analysis (Fig. 1). Two sites – the areas surrounding Lake Vostok and Dome C (Fig. 1b and c) – are located within the interior of the East Antarctic ice sheet and are characterized by relatively simple topography. These sites allowed us to evaluate the performance of Delay-Doppler altimetry in regions representative of a large part of the Antarctic interior. Both sites have low and relatively uniform topographic slopes, with an average and standard deviation of 0.09 and 0.05∘ (Lake Vostok) and 0.04 and 0.03∘ (Dome C), respectively, based upon a 1 km digital elevation model (DEM) (Slater et al., 2018). Furthermore, the flat ice surface directly above Lake Vostok represents an established validation site for new altimetry missions (Richter et al., 2014; Schröder et al., 2017; Shuman et al., 2006), although it is important to note that our study site does extend beyond the region floating in hydrostatic equilibrium, so as to incorporate areas of grounded ice with topography that is more representative of the ice sheet interior. To assess performance in regions of steeper and more complex topography, we then selected two coastal sites covering parts of Dronning Maud Land and Wilkes Land (Fig. 1d and e). These locations were chosen because of the availability of airborne campaigns that could be used as independent validation. Both sites have an order of magnitude steeper and less uniform topography than the inland sites, with the mean and standard deviation of the surface slope being 0.50 and 0.94∘ (Dronning Maud Land) and 0.40 and 0.51∘ (Wilkes Land), respectively. Finally, beyond these focused, site-specific studies, we also conducted several continent-wide analyses, in order to better understand the performance of the Sentinel-3A altimeter across a broader range of topographic regimes.
To evaluate the accuracy of Sentinel-3A elevation measurements, we began by analysing 14 cycles of Sentinel-3A SRAL data acquired between December 2016 and December 2017. Our processing followed a standard chain, beginning with the 20 Hz waveform data provided by the European Space Agency (ESA) within their freely distributed “enhanced” data file and generated using their Processing Baseline 2.27. Firstly, an estimate of the waveform noise was made from the mean power of the lowest six waveform samples, and waveforms where this value exceeded 0.3 of the maximum recorded power were rejected and did not pass to the subsequent processing. Next, each remaining waveform was oversampled by a factor of 100 using a spline interpolation, and then the leading edge of each waveform was identified based upon the first set of waveform samples that satisfied the following criteria: (1) a normalized power that exceeded the noise floor (defined as being 0.05 above the mean normalized power of the lowest 6 samples); (2) a change in normalized power from the noise floor to the next waveform peak that was greater than 0.2; and (3) an increase in power that ended with a defined waveform peak (such that there was a decrease in power at delay times beyond the peak location). Each waveform that had a leading edge satisfying these criteria was then retracked using several empirical retrackers, namely a threshold on the offset centre of gravity amplitude (TCOG) (Wingham et al., 1986), a threshold first maximum retracker (TFMRA) (Helm et al., 2014), and a maximum gradient of the first leading edge retracker (Gray et al., 2015). For the first two solutions, a threshold of 50 % of the leading edge power (defined as the leading edge maximum power minus the noise floor) was used, with the aim of providing a stable retracking point across both low-slope and more complex topographic surfaces. More specifically, we chose this mid-power threshold as a balance between minimizing the sensitivity to noise at the start of the waveform leading edge and reducing the impact of radar speckle (which is more apparent near to the waveform peak due to its multiplicative nature). For the majority of this study, we have focused on reporting results produced by the TCOG retracking because of the continuity it provides with the ground segments of past European Space Agency (ESA) missions and the broadly consistent results between all three of the retrackers tested. However, for completeness we do also report statistics from all retrackers within our independent validation exercise.
After retracking, Level-2 instrument and geophysical corrections were applied to each range measurement to account for the distance between the antenna and satellite centre of mass, dry and wet troposphere delays, ionosphere delays, solid Earth tide, ocean loading tide and polar tide, plus ocean tide and the inverse barometer effect over floating ice. These corrections are all included within the enhanced data product, and further details can be found within the product specification (Sentinel-3 Core PDGS Instrument Processing Facility (IPF), 2015). These geophysical corrections are provided at 1 Hz sampling, and so we used linear interpolation to resample these fields to the native 20 Hz rate of the altimeter measurements. The echoing point was then relocated to the point of closest approach (Roemer et al., 2007) within the SAR beam footprint using a DEM derived from 7 years of CryoSat-2 data (Slater et al., 2018), with echoing points that were relocated by more than ∼8 km, and therefore at the edge of antenna beamwidth, removed (∼3.7 % of data, based on statistics from a single cycle), together with relocated elevations that deviated by more than 100 m from the DEM (∼3.0 % of data, based on statistics from a single cycle). This dataset then formed the basis of the subsequent validation activities.
To assess the utility of the Sentinel-3A altimeter over ice surfaces, we firstly investigated the precision of the SRAL measurements by assessing their repeatability in space and time. For this purpose we performed two sets of analysis: (1) an evaluation of repeated profiles that crossed subglacial Lake Vostok, a site that provides a stable and low-slope surface that is well established for validation studies (Richter et al., 2014; Schröder et al., 2017; Shuman et al., 2006) and (2) a continent-wide single-cycle cross-over analysis, to evaluate the repeatability of measurements at locations where ascending and descending satellite passes intersect (Wingham et al., 1998; Zwally et al., 1989).
4.1 Shot-to-shot precision at Lake Vostok
We analysed repeated altimeter profiles that crossed the ice surface above the Lake Vostok site in East Antarctica, in order to assess the SRAL instrument precision. The smooth, flat surface above the lake minimizes the influence of topography and allowed us to focus primarily on the performance of the SRAL instrument itself, and specifically to understand the impact of radar speckle, small-scale variations in the firn backscattering properties and the influence of retracker imprecision on the SAR altimeter measurements. Between December 2016 and December 2017, the satellite made 14 passes over the lake, and so, focusing on two ground tracks that crossed the flat (< 0.01∘) central portion of the lake (Fig. 2), we assessed the repeatability of these measurements in space and time. For each ground track, we used the 14 repeated elevation profiles to compute (1) the mean elevation profile, (2) the residual elevations from the mean profile, and (3) the standard deviations of all elevation measurements within 400 m intervals along track (Fig. 2). Together, these provide an assessment of the instrument shot-to-shot precision over ice sheet surfaces under the influence of minimal topography and also the repeatability of measurements through time. We find that, over the first year of routine operations, the SAR altimeter has operated with sub-decimetre precision. On average, the 14-cycle 1σ standard deviation along both tracks was 7 cm and rarely fell outside of the range of 5–10 cm along the entirety of the track segments analysed (Fig. 2).
4.2 Single-cycle cross-over analysis
Next, we used a single-cycle cross-over analysis to assess the repeatability of measurements at all locations where ascending and descending satellite passes crossed. Elevation differences at orbital cross-overs are commonly used as a metric for measurement precision and integrate a number of factors, including spatially uncorrelated orbit errors, retracker imprecision, the impact of radar speckle, echo relocation errors, and any sensitivity to anisotropic scattering within the near-surface snowpack (Armitage et al., 2014; Remy et al., 2012; Wingham et al., 1998). We selected two 27-day cycles of Sentinel-3A data (cycles 12 and 24) and for each cycle computed the elevation difference at all locations where ascending and descending passes crossed. Specifically, we (1) removed outlier measurements that deviated by more than 50 m from our DEM (Slater et al., 2018), (2) identified crossing points in the remaining dataset as the intersection between two consecutive measurements of an ascending pass and two consecutive measurements of a descending pass, and (3) computed the elevation difference by interpolating the bracketing ascending and descending records to the cross-over location. Finally, for each cycle we binned all the cross-over differences within 0.2∘ intervals of the surface slope to investigate the relationship between the magnitude of the surface slope and the cross-over elevation precision. The results of this analysis are shown in Fig. 3 and Table 2.
In total, approximately 90 000 cross-overs were identified during each cycle. At the continental scale, the distribution of these elevation differences is non-normal (0.01 % significance level). Cross-overs from both cycles have a median difference of < 1 cm in magnitude, and a higher than normal proportion of the differences are clustered around this central value, reflecting the good repeatability of measurements across the low-slope interior of the ice sheet. For example, across the entirety of the ice sheet, the median absolute deviation of cross-overs from the median value is ∼0.3 m, and 80 % of all cross-over differences are less than 1.8 m. At the extremes of these distributions, however, there are a number of outliers, with around 5 % of cross-overs having an elevation difference that exceeds 10 m. Generally, the magnitude of the cross-over differences increases with surface slope (Fig. 3d), with the largest differences occurring in regions with steep and complex coastal topography. In these regions, the processes of locating the echoing point within the beam footprint and of retracking complex multi-peaked waveforms become more challenging. These remain active topics of research, which are likely to deliver further improvements to ice sheet SAR altimetry in the future. In the meantime, we note that for many glaciological applications it may be beneficial to remove these outliers and, although not the focus of this study, we touch upon possible filtering strategies within Sect. 5.3.
5.1 Reference dataset and methods
To conduct an independent evaluation of the accuracy of our Sentinel-3A ice sheet measurements, we used elevation data acquired by the Airborne Topographic Mapper (ATM) and Riegl Laser Altimeter (RLA) instruments carried on Operation IceBridge campaigns flown between 2009 and 2016. At our two inland study sites we used the Level-2 ATM product, which provides surface elevation measurements with an along-track sampling every 0.25 s (equating to a ∼30 m interval for a typical aircraft velocity) and an 80 m across-track platelet at nadir. The Level-2 product was chosen because the processing includes a smoothing of the Level-1b data, which reduces the impact of uncorrelated shot-to-shot noise on our validation dataset and also brings the measurement cell closer to the resolution of the SAR footprint. The ATM measurements have been estimated to have a vertical accuracy and precision of 7 and 3 cm, respectively (Martin et al., 2012). At our two coastal sites, where ATM measurements have not been acquired, we instead used RLA acquisitions. This instrument has a smaller ground footprint of 25 m along track by 1 m across track and a slightly larger reported accuracy of 12 cm (Blankenship et al., 2012).
1 Vostok IceBridge measurements were acquired during the months of 10/2013 and 11/2013. 2 Dome C IceBridge measurements were acquired during the months of 10/2013 and 11/2013. 3 Dronning Maud Land IceBridge measurements were acquired during the month of 2/2011. 4 Wilkes Land IceBridge measurements were acquired during the months of 1/2009, 2/2009, 12/2009, 1/2010, 12/2010, 1/2011, 12/2011, 11/2012 and 12/2012.
To compute elevation differences between our Sentinel-3A and IceBridge datasets, we identified IceBridge records within a 200 m search radius of each satellite measurement. Where multiple IceBridge records existed within the search radius, we selected the closest measurement. Alternative methods for selecting IceBridge measurements were also tested, such as bilinear interpolation of multiple surrounding measurements, but this approach produced a less comprehensive set of comparison points from which to generate our assessment statistics. As part of this process, we identified and removed anomalous IceBridge elevation records that deviated by more than 100 m from an independent DEM (Slater et al., 2018). This step removed 7.3 % of the total Antarctic IceBridge dataset. We then corrected for elevation differences arising from the spatial and temporal separation of the satellite and airborne measurements. In the case of the former, we constructed a bicubic interpolation of the surrounding 4 pixel × 4 pixel area of the DEM surface and used this to estimate the difference in elevation between the satellite and airborne measurement locations. For the latter, we used an estimate of the local rate of elevation change (McMillan et al., 2014). The magnitude of the elevation change correction is small (< 1 cm yr−1) at the Vostok, Dome C, and Dronning Maud Land sites. At our Wilkes Land site, the magnitude of the correction is larger (8 cm yr−1), and it is therefore possible that inaccuracies in the correction could contribute, in part, to the differences between the airborne and satellite measurements at this site. For example, a 10 % error in the correction over a 5-year period would equate to a 4 cm error in the corrected IceBridge elevation. Finally, we computed the Sentinel-3 minus IceBridge elevation difference for each measurement pair, and so generated a set of statistics for each study site (Table 3). Because the differences, particularly at coastal sites, are not normally distributed and exhibit higher clustering around the central value, together with a greater proportion of outliers, we principally use the median and median absolute deviation (MAD) from the median as measures of the bias and dispersion, respectively. We choose to use the MAD because for non-normal distributions this statistic provides a more representative measure of the midpoint of the dispersion. We do, nonetheless, report both the MAD and the standard deviation within Table 3.
5.2 Evaluation at inland sites
At the inland sites of Lake Vostok and Dome C, we find very good agreement between the Sentinel-3A and airborne datasets (Fig. 4 and Table 3). At Lake Vostok, the median bias between our TCOG solution and the airborne data is 1 cm, and the MAD dispersion of the differences is 13 cm. At Dome C, the bias is larger (20 cm), but the dispersion of the differences is smaller (6 cm). The differing bias between the two sites is investigated in more detail in Sect. 5.4. Comparing the results from the different retrackers, we find variations of approximately 10–30 cm in the median bias, which reflects differences in the algorithms used to select the retracking point on the waveform leading edge. There is nonetheless relatively little difference between retrackers in the MAD of the elevation differences, which is typically of the order of 10 cm in magnitude (Table 3). This is consistent with our previous analysis of the instrument precision above Lake Vostok and suggests that at these relatively low-slope inland sites, uncorrelated sources of error, for example due to imprecision of the retracker, radar speckle, the process of measurement relocation, or small-scale variations in snowpack characteristics, have not significantly affected the SAR altimeter elevation measurements. In total, we find that more than 97 % of validation points (TCOG retracker) have an elevation difference of less than 50 cm (Table 3), and that 70 % (Vostok) and 49 % (Dome C) have a difference of less than 20 cm.
5.3 Evaluation at coastal sites
At the coastal sites of Dronning Maud Land and Wilkes Land the differences between the Sentinel-3A and airborne datasets are, as expected, more widely dispersed than at our inland study locations (Fig. 4 and Table 3). At these sites, the more rugged coastal topography can produce complex waveforms, as energy is often returned from several distinct surfaces within the illuminated beam footprint. These factors represent well-established challenges for radar altimetry, both for retracking algorithms and for the procedure of correctly locating the on-the-ground origin of the derived elevation measurement. SAR altimetry, due to its smaller ground footprint, has the potential to be less affected by these topographic influences, and indeed we find that the overall median biases relative to IceBridge remain small (Table 3), namely 0.03 and 0.12 m at Dronning Maud Land and Wilkes Land, respectively (TCOG retracking). The magnitude of these biases is comparable to those found at our inland sites, suggesting that for a metric that is robust to outliers, no systematic bias is introduced as large-scale topographic complexity increases.
For our coastal sites, the dispersion of the elevation differences relative to IceBridge is larger, as indicated by the MAD values of 0.30 and 0.74 m (TCOG retracking) at Dronning Maud Land and Wilkes Land, respectively (Table 3). Nonetheless, these first results demonstrate that even in these more challenging regions, the MAD precision of SAR elevation measurements is well below 1 m. At these sites, we find that ∼60 %–85 % of the validated satellite elevation measurements (TCOG retracking) are within 1 m of their airborne counterpart, and 92 %–98 % are within 10 m (Table 3). As is evident from these statistics, and also the standard deviation of the differences (Table 3), there are a small number of outlying measurements that exhibit larger deviations from the airborne validation data. Given that the primary focus of this study is on assessing data quality, we have chosen not to remove these outliers, although we note that for many future applications it may be beneficial to implement filtering procedures to do so. For example, when we consider an Antarctic-wide evaluation (see Sect. 5.4), we find that removing Sentinel-3A points that deviate by more than 10 m from our pre-existing DEM achieves a 70 % reduction in the standard deviation of the Sentinel-3 minus IceBridge differences. Although a small proportion of outliers still remain even after this filtering step (∼1 % of data deviate by more than 20 m), either because the DEM fails to identify them as outliers or because the IceBridge data themselves are inaccurate, it is evident that a simple post-processing strategy such as this may be beneficial for many glaciological applications.
5.4 Influence of surface topography on Sentinel-3 measurements
To better understand the variations between the Sentinel-3A and IceBridge measurements at our different sites, we considered the influence of both large-scale (wavelength much greater than the beam footprint, i.e. surface slope) and small-scale (wavelength equal to, or shorter than, the beam footprint) topography. To investigate the former, we expanded our site-specific analysis to calculate elevation differences between one complete cycle of Sentinel-3A data (cycle 12, elevations derived using the TCOG retracker) and all ATM and RLA measurements acquired across the Antarctic ice sheet between 2009 and 2016, in order to assess the sensitivity of the elevation differences to surface slope (Fig. 5). At the continent scale we find a Sentinel-3A minus IceBridge median elevation difference of 0.06 m and a MAD of the elevation differences of 1.06 m. It is important to note that these statistics do not represent an unbiased sample of the total ice sheet distribution. The IceBridge surveys have more frequently been flown across steep and more complex ice margin topography (Fig. 5b), resulting in a median slope at the locations where the Sentinel-3A data are validated by IceBridge that is 37 % higher than the median slope of the ice sheet as a whole.
Nonetheless, the IceBridge dataset covers a range of surface slopes and allows us to investigate the relationship between the magnitude of the surface slope and the Sentinel-3A–IceBridge elevation differences. Specifically, we grouped the validation data within 0.2∘ slope intervals (using the slope model presented in Slater et al., 2018), and computed the median absolute elevation difference across each band (Fig. 5). Unsurprisingly, the magnitude of the median elevation difference increases with the magnitude of the surface slope, mirroring the trend in precision that was apparent in our cross-over data and reflecting the challenges of retracking and reliably locating the echoing point in areas of steeply sloping terrain. For surface slopes below 0.4∘, which constitutes 74 % of the ice sheet area, the median absolute difference between Sentinel-3 and IceBridge is less than 1 m, whereas for highly sloping terrain greater than 1∘ the difference increases to ∼10 m. We note that these statistics are specifically related to the current processing baseline and that expected improvements to the Sentinel-3 processing chain are likely to improve performance in highly sloping areas in the future (see further discussion below).
Considering the validation statistics across all four study sites, the pattern of increasing dispersion of elevation differences at coastal locations is consistent with our understanding that measurement precision degrades with increasing topographic complexity. In contrast, it was less expected that, for any given retracker, we find a difference of ∼20 cm in the bias recorded at the Lake Vostok and Dome C sites, especially because the large-scale topographic characteristics (surface slope < 0.1∘) and climatological setting (cold, dry ice sheet interior) are similar. Although an inaccurate elevation rate correction could introduce a bias into our measurements (due to the difference in date between the Sentinel-3A and IceBridge acquisitions), the median elevation rates at these inland sites are small (< 1 cm yr−1 at both Dome C and Vostok), and so we believe it is unlikely that this is the source of the observed difference. To investigate other possible physical explanations for this difference, we therefore used the IceBridge data to assess the finer-scale topography at both of these inland sites. Specifically, we estimated the long wavelength signal along each elevation profile by fitting a quadratic curve to the data, and then plotted the elevation residuals having removed the modelled long-wavelength topography (Fig. 6). For this analysis, we focused specifically on the inland sites because (1) they exhibited a differing, and unexplained, bias relative to IceBridge and (2) at long wavelengths they are relatively flat, and so any impact of small-scale roughness is likely to be more evident than at coastal sites with more complex long-wavelength topography. The airborne flight lines (Fig. 6) show that Dome C presents a much rougher surface at 100–500 m length scales, with amplitudes typically ranging from ∼5 to 30 cm. In comparison, the amplitude of oscillations at Lake Vostok is much smaller, typically 1–5 cm. Along these profiles the standard deviation of the residuals is 1.7 and 6.9 cm for Lake Vostok and Dome C, respectively, indicating that by this metric of surface roughness, Dome C is ∼4 times rougher than Lake Vostok. It is therefore possible that part of the larger bias at Dome C can be explained by the rougher surface and the tendency of the satellite altimeter, given its larger footprint, to be more influenced by the local topographic peaks than the airborne instrument. Finally, it should be noted that any of the retracking algorithms can likely be tuned to reduce the bias at a particular site, for example by selecting a higher retracking point on the waveform leading edge that is closer to the theoretical mean return from a surface with these roughness characteristics. However, we reiterate that our philosophy here is to use a conservative retracking threshold that is likely to deliver robust and stable results across all types of topographic regimes, and one that is therefore well suited to delivering reliable continent-wide estimates of surface elevation change through time.
Our analysis has provided the first comprehensive assessment of ice sheet elevation measurements that have been derived using Sentinel-3A Delay-Doppler (SAR mode) altimetry and an initial demonstration of their accuracy and precision across a range of topographic regimes. Within a wider geophysical context, one of the principle uses of altimetry data is to determine changes in ice sheet elevation over time (Flament and Rémy, 2012; Shepherd et al., 2012; Shepherd and Wingham, 2007; Zwally et al., 2005). Although the available time span of Sentinel-3A acquisitions is short for detailed glaciological interpretation of any signals, it is nonetheless of interest to establish (1) the extent to which the system has provided a stable measurement platform in time and (2) whether the precision and accuracy of the SAR mode of operation is sufficient to be able to resolve known signals and modes of glaciological change. As a preliminary investigation of these questions, we therefore applied a modified model-fit method (McMillan et al., 2014, 2016) to all Sentinel-3A data acquired up to and including cycle 32, in order to explore the potential of these data for mapping elevation changes of the Antarctic ice sheet. In summary, we firstly removed an a priori estimate of elevation from each measurement using an auxiliary DEM (Slater et al., 2018), rejecting records that deviated by more than 50 m from the same DEM. We then used the resulting elevation residuals to simultaneously solve for linear spatial and temporal rates of elevation change on a 5 km × 5 km grid. We rejected grid cells for which the model produced a poor or geophysically unrealistic fit to the data, defined to be where the root mean square of the observed-minus-modelled residuals exceeded 2 m, the absolute rate of elevation change exceeded 10 m yr−1, the spatial gradient computed from the elevation residuals exceeded 5∘, or where fewer than 20 measurements constrained the model fit.
Using this method, we determined an estimate of the rate of ice sheet surface elevation change across a total area of 5 061 700 km2, constituting 42.3 % of the ice sheet (Fig. 7). Across large parts of the slow-flowing ice sheet interior, the derived rates of elevation change are low. This agrees with numerous recent studies (Flament and Rémy, 2012; Helm et al., 2014; McMillan et al., 2014) and provides an early indication that the Sentinel-3 instrument and orbital configuration is suitable for mapping changes across the low relief ice sheet interior. Although we believe that the Sentinel-3 record is still too short to perform a detailed, ice sheet-wide, quantitative inter-comparison relative to previously published altimeter datasets, we do find evidence that Delay-Doppler altimetry is able to map the higher, dynamically driven, rates of elevation change that are occurring across coastal regions of the ice sheet (Flament and Rémy, 2012; Helm et al., 2014; McMillan et al., 2014). In particular, we observe widespread elevation change across the fast-flowing ice streams draining into the Amundsen Sea Sector of West Antarctica, with rates of surface lowering of 2–3 m yr−1 close to the grounding line of Pine Island Glacier and 4–5 m yr−1 upstream of the grounding lines of the Thwaites Glacier and Smith Glacier. In comparison to inland regions, the coverage is generally poorer across the steeper ice margin regions and the more mountainous terrain of the Antarctic Peninsula. This is expected given the lower precision of elevation measurements in these locations, the wider track spacing, and the short time period over which the trends are being computed. Additionally, it is important to note that the current approach to performing the waveform windowing and SAR multi-looking within the Level-1B processing chain of the ground segment is not fully optimized for these challenging ice regions. Further refinements to this processing step, namely to adjust the windowing during the Doppler beam stacking to account for large variations in the satellite tracker range, are currently being implemented and are expected to deliver future improvements in data retrieval in these regions. Based upon this preliminary assessment, however, there is good reason to expect that Sentinel-3 Delay-Doppler altimetry will prove to be an effective tool for mapping ice sheet elevation change.
Finally, we investigated the capability of SAR altimetry to make precise measurements of surface elevation change within the inland regions of the ice sheet, where previously only low-resolution-mode observations have been available. Specifically, we focused on a small region (of the order of tens of kilometres) of anomalously high-elevation changes at a location within the interior of the East Antarctic ice sheet (location shown in Fig. 7a). We analysed 28 cycles of Sentinel-3A data passing over this region, grouping data from all cycles within 340 m intervals along track and again employing a model-fit approach to isolate the temporal evolution of the ice surface with a 27-day repeat frequency (Moholdt et al., 2010; Smith et al., 2009). On either side of this region of high-elevation change, we find a high level of repeatability of the SAR elevation measurements, giving us confidence in the precision of the instrument and our method for isolating along-track elevation changes through time (Fig. 7b). Over the feature itself, we resolve a progressive, spatially coherent movement of the ice sheet surface, with a total lowering of ∼1.7 m over a period of 16 months. Transient changes in elevation at this spatial scale are widely understood to be caused by subglacial lake drainage (Smith et al., 2009), and as such, our observations provide a first indication of the capability of Sentinel-3 SAR altimetry to systematically monitor such events.
We have undertaken a first assessment of the utility of Sentinel-3 Delay-Doppler (SAR mode) altimetry for measuring ice sheet elevation and elevation change using the standard ESA Level-1b product and our own Level-2 processing chain. Analysis of repeated acquisitions over the Lake Vostok validation site indicates that, over the first year of routine operations, the instrument has operated with sub-decimetre precision. Through validation with airborne campaigns, we find small median biases in elevation, typically of the order 1–10 cm, at both inland and coastal sites. The dispersion of elevation residuals, measured with respect to the validation data, is of the order of 10 cm at inland sites, increasing to ∼1 m at coastal sites with more complex topography. This reflects the main challenges associated with processing radar altimetry data in complex ice margin regions, namely (1) reliably retracking multi-peaked waveforms that arise when multiple distinct surface reflections are captured within the receive window and (2) accurately establishing the location of the echoing point within the SAR beam footprint. These represent principle avenues of future research within the field of ice sheet Delay-Doppler altimetry. Nonetheless, the accuracy achieved in even these complex ice margin regions is encouraging and expected to improve further as refinements are made to the operational ground segment processing. Finally, we have shown the capability of Sentinel-3, albeit with the relatively short record of data currently available, to resolve the known signals of elevation change that currently dominate Antarctica's contribution to sea level rise and to monitor subglacial lake activity. Together, our analysis demonstrates the early promise of Sentinel-3 SAR altimetry as a platform for the long-term, operational monitoring of Earth's ice sheets.
The Sentinel-3A altimetry data used in this study are freely available through the Copernicus Open Access Hub (https://scihub.copernicus.eu/dhus/\#/home, Copernicus Open Access Hub, 2018). The IceBridge airborne altimetry data used in this study are freely available from the US National Snow and Ice Data Center (https://nsidc.org/, Blankenship et al., 2012; Studinger, 2014). The CryoSat-2 DEM used in this study is freely distributed by the UK NERC Centre for Polar Observation and Modelling (http://www.cpom.ucl.ac.uk/csopr/icesheets2, Slater et al., 2018).
MM designed the experiments. MM and AM processed and analysed the data. MM prepared the manuscript with contributions from AS, MR, AA, and JB, and all authors reviewed the manuscript.
The authors declare that they have no conflict of interest.
This work was supported by the UK NERC Centre for Polar Observation and
Modelling, the European Space Agency contract “SEOM – Sentinel-3
Performance Improvements for ICE sheets” (contract number
4000115201/15/I-BG), and the Sentinel-3 Mission Performance Centre.
We thank the editor Ginny Catannia and two anonymous reviewers for their
comments, which have substantially improved the manuscript.
Edited by: Ginny Catania
Reviewed by: two anonymous referees
Armitage, T. W. K., Wingham, D. J., and Ridout, A. L.: Meteorological Origin of the Static Crossover Pattern Present in Low-Resolution-Mode CryoSat-2 Data Over Central Antarctica, IEEE Geosci. Remote S., 11, 1295–1299, https://doi.org/10.1109/LGRS.2013.2292821, 2014.
Bamber, J. L. and Bindschadler, R. A.: An improved elevation dataset for climate and ice-sheet modelling: validation with satellite imagery, Ann. Glaciol., 25, 438–444, 1997.
Bamber, J. L., Gomez-Dans, J. L., and Griggs, J. A.: A new 1 km digital elevation model of the Antarctic derived from combined satellite radar and laser data – Part 1: Data and methods, The Cryosphere, 3, 101–111, https://doi.org/10.5194/tc-3-101-2009, 2009.
Blankenship, D. D., Kempf, S. D., Young, D. A., Roberts, J. L., van Ommen, T., Forsberg, R., Siegert, M. J., Palmer, S. J., and Dowdeswell, J. A.: IceBridge Riegl Laser Altimeter L2 Geolocated Surface Elevation Triplets, Version 1, Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center, https://doi.org/10.5067/JV9DENETK13E, 2012 (updated 2013).
Bogning, S., Frappart, F., Blarel, F., Niño, F., Mahé, G., Bricquet, J.-P., Seyler, F., Onguéné, R., Etamé, J., Paiz, M.-C., and Braun, J.-J.: Monitoring Water Levels and Discharges Using Radar Altimetry in an Ungauged River Basin: The Case of the Ogooué, Remote Sens., 10, 350, https://doi.org/10.3390/rs10020350, 2018.
Bonnefond, P., Laurain, O., Exertier, P., Boy, F., Guinle, T., Picot, N., Labroue, S., Raynal, M., Donlon, C., Féménias, P., Parrinello, T., and Dinardo, S.: Calibrating the SAR SSH of Sentinel-3A and CryoSat-2 over the Corsica Facilities, Remote Sens., 10, 92, https://doi.org/10.3390/rs10010092, 2018.
Chelton, D. B., Walsh, E. J., and MacArthur, J. L.: Pulse Compression and Sea Level Tracking in Satellite Altimetry, J. Atmos. Ocean. Tech., 6, 407–438, https://doi.org/10.1175/1520-0426(1989)006<0407:PCASLT>2.0.CO;2, 1989.
Copernicus Open Access Hub: Copernicus Open Access Data Hub System 0.14.7-2, created by a Serco and GAEL Systems consortium under a contract with the European Space Agency, funded by the EU and ESA, Dataset, available at: https://scihub.copernicus.eu/dhus/#/home, last access: Septeber 2018.
Davis, C. and Ferguson, A.: Elevation change of the Antarctic ice sheet, 1995–2000, from ERS-2 satellite radar altimetry, IEEE T. Geosci. Remote, 42, 2437–2445, https://doi.org/10.1109/TGRS.2004.836789, 2004.
Davis, C. H., Li, Y. H., McConnell, J. R., Frey, M. M., and Hanna, E.: Snowfall-driven growth in East Antarctic ice sheet mitigates recent sea-level rise, Science, 308, 1898–1901, https://doi.org/10.1126/science.1110662, 2005.
Donlon, C., Berruti, B., Buongiorno, A., Ferreira, M.-H., Féménias, P., Frerick, J., Goryl, P., Klein, U., Laur, H., Mavrocordatos, C., Nieke, J., Rebhan, H., Seitz, B., Stroede, J., and Sciarra, R.: The Global Monitoring for Environment and Security (GMES) Sentinel-3 mission, Remote Sens. Environ., 120, 37–57, https://doi.org/10.1016/j.rse.2011.07.024, 2014.
Flament, T. and Rémy, F.: Dynamic thinning of Antarctic glaciers from along-track repeat radar altimetry, J. Glaciol., 58, 830–840, https://doi.org/10.3189/2012JoG11J118, 2012.
Gray, L., Burgess, D., Copland, L., Demuth, M. N., Dunse, T., Langley, K., and Schuler, T. V.: CryoSat-2 delivers monthly and inter-annual surface elevation change for Arctic ice caps, The Cryosphere, 9, 1895–1913, https://doi.org/10.5194/tc-9-1895-2015, 2015.
Haran, T., Bohlander, J., Scambos, T., Painter, T., and Fahnestock, M.: MODIS mosaic of Antarctica (MOA) image map, National Snow and Ice Data Center, Boulder, Colorado, USA 2006.
Helm, V., Humbert, A., and Miller, H.: Elevation and elevation change of Greenland and Antarctica derived from CryoSat-2, The Cryosphere, 8, 1539–1559, https://doi.org/10.5194/tc-8-1539-2014, 2014.
Heslop, E. E., Sánchez-Román, A., Pascual, A., Rodríguez, D., Reeve, K. A., Faugère, Y., and Raynal, M.: Sentinel-3A Views Ocean Variability More Accurately at Finer Resolution, Geophys. Res. Lett., 44, 12367–12374, https://doi.org/10.1002/2017GL076244, 2017.
Hurkmans, R. T. W. L., Bamber, J. L., Davis, C. H., Joughin, I. R., Khvorostovsky, K. S., Smith, B. S., and Schoen, N.: Time-evolving mass loss of the Greenland Ice Sheet from satellite altimetry, The Cryosphere, 8, 1725–1740, https://doi.org/10.5194/tc-8-1725-2014, 2014.
Johannessen, O. M., Khvorostovsky, K., Miles, M. W., and Bobylev, L. P.: Recent ice-sheet growth in the interior of Greenland, Science, 310, 1013–1016, https://doi.org/10.1126/science.1115356, 2005.
Martin, C. F., Krabill, W. B., Manizade, S. S., Russell, R. L., Sonntag, J. G., Swift, R. N., and Yungel, J. K.: Airborne Topographic Mapper Calibration Procedures and Accuracy Assessment, NASA Technical Report NASA/TM/u20132012-215891, Goddard Space Flight Center, Greenbelt, Maryland 20771, available at: https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/20120008479.pdf (last access: 9 November 2015), 2012.
McMillan, M., Shepherd, A., Sundal, A., Briggs, K., Muir, A., Ridout, A., Hogg, A., and Wingham, D.: Increased ice losses from Antarctica detected by CryoSat-2, Geophys. Res. Lett., 41, 1–7, https://doi.org/10.1002/2014GL060111, 2014.
McMillan, M., Leeson, A., Shepherd A., Briggs, K., Armitage, T. W. K., Hogg, A., Kuipers Munneke, P., van den Broeke, M., Noël, B., van de Berg, W. J., Ligtenberg, S., Horwath, M., Groh, A., Muir, A., and Gilbert, L.: A high-resolution record of Greenland mass balance, Geophys. Res. Lett., 43, 7002–7010, https://doi.org/10.1002/2016GL069666, 2016.
McMillan, M., Shepherd, A., Muir, A., Gaudelli, J., Hogg, A. E., and Cullen, R.: Assessment of CryoSat-2 interferometric and non-interferometric SAR altimetry over ice sheets, Adv. Space Res., 62, 1281–1291, https://doi.org/10.1016/j.asr.2017.11.036, 2017.
Moholdt, G., Nuth, C., Hagen, J. O., and Kohler, J.: Recent elevation changes of Svalbard glaciers derived from ICESat laser altimetry, Remote Sens. Environ., 114, 2756–2767, https://doi.org/10.1016/j.rse.2010.06.008, 2010.
Price, S. F., Payne, A. J., Howat, I. M., and Smith, B. E.: Committed sea-level rise for the next century from Greenland ice sheet dynamics during the past decade, P. Natl. Acad. Sci. USA, 108, 8978–83, https://doi.org/10.1073/pnas.1017313108, 2011.
Raney, K.: The delay/doppler radar altimeter, IEEE T. Geosci. Remote, 36, 1578–1588, https://doi.org/10.1109/36.718861, 1998.
Remy, F., Mazzega, P., Houry, S., Brossier, C., and Minster, J. F.: Mapping of the Topography of Continental Ice by Inversion of Satellite-Altimeter Data, J. Glaciol., 35, 98–107, https://doi.org/10.3189/002214389793701419, 1989.
Remy, F., Flament, T., Blarel, F., and Benveniste, J.: Radar altimetry measurements over antarctic ice sheet: A focus on antenna polarization and change in backscatter problems, Adv. Space Res., 50, 998–1006, https://doi.org/10.1016/j.asr.2012.04.003, 2012.
Richter, A., Popov, S. V., Fritsche, M., Lukin, V. V., Matveev, A. Y., Ekaykin, A. A., Lipenkov, V. Y., Fedorov, D. V., Eberlein, L., Schröder, L., Ewert, H., Horwath, M., and Dietrich, R.: Height changes over subglacial Lake Vostok, East Antarctica: Insights from GNSS observations, J. Geophys. Res., 119, 2460–2480, https://doi.org/10.1002/2014JF003228.Received, 2014.
Roemer, S., Legrésy, B., Horwath, M., and Dietrich, R.: Refined analysis of radar altimetry data applied to the region of the subglacial Lake Vostok/Antarctica, Remote Sens. Environ., 106, 269–284, https://doi.org/10.1016/j.rse.2006.02.026, 2007.
Schröder, L., Richter, A., Fedorov, D. V., Eberlein, L., Brovkov, E. V., Popov, S. V., Knöfel, C., Horwath, M., Dietrich, R., Matveev, A. Y., Scheinert, M., and Lukin, V. V.: Validation of satellite altimetry by kinematic GNSS in central East Antarctica, The Cryosphere, 11, 1111–1130, https://doi.org/10.5194/tc-11-1111-2017, 2017.
Sentinel-3 Core PDGS Instrument Processing Facility (IPF): Implementation Product Data Format Specification – SRAL/MWR Level 1 and 2 Instrument Products, available at: https://sentinel.esa.int/documents/247904/2753172/Sentinel-3-Product-Data-Format-Specification-Level-2-Land (last access: 24 September 2018), 2015.
Shepherd, A. and Wingham, D. Recent sea-level contributions of the Antarctic and Greenland ice sheets, Science, 315, 1529–1532, 2007.
Shepherd, A., Wingham, D., and Rignot, E.: Warm ocean is eroding West Antarctic Ice Sheet, Geophys. Res. Lett., 31, L23402, https://doi.org/10.1029/2004gl021106, 2004.
Shepherd, A., Ivins, E. R., Geruo, A., Barletta, V. R., Bentley, M. J., Bettadpur, S., Briggs, K. H., Bromwich, D. H., Forsberg, R., Galin, N., Horwath, M., Jacobs, S., Joughin, I., King, M. A., Lenaerts, J. T., Li, J., Ligtenberg, S. R., Luckman, A., Luthcke, S. B., McMillan, M., Meister, R., Milne, G., Mouginot, J., Muir, A., Nicolas, J.P., Paden, J., Payne, A. J., Pritchard, H., Rignot, E., Rott, H., Sørensen, L. S., Scambos, T. A., Scheuchl, B., Schrama, E. J., Smith, B., Sundal, A. V., Van Angelen, J. H., Van De Berg, W. J., Van Den Broeke, M. R., Vaughan, D. G., Velicogna, I., Wahr, J., Whitehouse, P. L., Wingham, D. J., Yi, D., Young, D., and Zwally, H. J.: A reconciled estimate of ice-sheet mass balance, Science, 338, 1183–1189, https://doi.org/10.1126/science.1228102, 2012.
Shuman, C. A., Zwally, H. J., Schutz, B. E., Brenner, A. C., DiMarzio, J. P., Suchdeo, V. P., and Fricker, H. A.: ICESat Antarctic elevation data: Preliminary precision and accuracy assessment, Geophys. Res. Lett., 33, L07501, https://doi.org/10.1029/2005GL025227, 2006.
Slater, T., Shepherd, A., McMillan, M., Muir, A., Gilbert, L., Hogg, A. E., Konrad, H., and Parrinello, T.: A new digital elevation model of Antarctica derived from CryoSat-2 altimetry, The Cryosphere, 12, 1551–1562, https://doi.org/10.5194/tc-12-1551-2018, 2018 (data available at: http://www.cpom.ucl.ac.uk/csopr/icesheets2/dems.html, last access: 15 February 2019).
Smith, B., Fricker, H. A., Joughin, I., and Tulaczyk, S.: An inventory of active subglacial lakes in Antarctica detected by ICESat (2003–2008), J. Glaciol., 55, 573–595, 2009.
Sørensen, L. S., Simonsen, S. B., Meister, R., Forsberg, R., Levinsen, J. F., and Flament, T.: Envisat-derived elevation changes of the Greenland ice sheet, and a comparison with ICESat results in the accumulation area, Remote Sens. Environ., 160, 56–62, https://doi.org/10.1016/j.rse.2014.12.022, 2015.
Studinger, M.: IceBridge ATM L2 Icessn Elevation, Slope, and Roughness, Version 2. Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center, https://doi.org/10.5067/CPRXXK3F39RV, 2014 (updated 2018).
The IMBIE team: Mass balance of the Antarctic Ice Sheet from 1992 to 2017, Nature, 558, 219–222, 2018.
Wingham, D., Ridout, A., Scharroo, R., Arthern, R., and Shum, C.: Antarctic elevation change from 1992 to 1996, Science, 282, 456–458, 1998.
Wingham, D., Wallis, D., and Shepherd, A.: Spatial and temporal evolution of Pine Island Glacier thinning, 1995–2006, Geophys. Res. Lett., 36, L17501, https://doi.org/10.1029/2009gl039126, 2009.
Wingham, D. J., Rapley, C. G., and Griffiths, H. D.: New techniques in satellite altimeter tracking systems, Proceedings of the IGARSS Symposium, Zurich, 8–11 September 1986, SP-254, 1339–1344, 1986.
Wingham, D. J., Francis, C. R., Baker, S., Bouzinac, C., Brockley, D., Cullen, R., de Chateau-Thierry, P., Laxon, S. W., Mallow, U., Mavrocordatos, C., Phalippou, L., Ratier, G., Rey, L., Rostan, F., Viau, P., and Wallis, D. W.: CryoSat: A mission to determine the fluctuations in Earth's land and marine ice fields, in: Natural Hazards and Oceanographic Processes from Satellite Data, edited by: Singh, M. and Shea, R. P., Elsevier Science Ltd., 37, 841–871, https://doi.org/10.1016/j.asr.2005.07.027, 2006a.
Wingham, D. J., Shepherd, A., Muir, A., and Marshall, G. J.: Mass balance of the Antarctic ice sheet, Philos. T. Roy. Soc. A, 364, 1627–1635, https://doi.org/10.1098/rsta.2006.1792, 2006b.
Zwally, H. J., Bindschadler, R. A., Brenner, A. C., Major, J. A., and Marsh, J. G.: Growth of greenland ice sheet: measurement, Science, 246, 1587–1589, https://doi.org/10.1126/science.246.4937.1587, 1989.
Zwally, H. J., Giovinetto, M. B., Li, J., Cornejo, H. G., Beckley, M. A., Brenner, A. C., Saba, J. L., and Yi, D.: Mass changes of the Greenland and Antarctic ice sheets and shelves and contributions to sea-level rise?: 1992–2002, J. Glaciol., 51, 1992–2002, 2005.
Zwally, H. J., Li, J., Brenner, A. C., Beckley, M., Cornejo, H. G., DiMarzio, J., Giovinetto, M. B., Neumann, T. A., Robbins, J., Saba, J. L., Yi, D., and Wang, W.: Greenland ice sheet mass balance: distribution of increased mass loss with climate warming; 2003–07 versus 1992–2002, J. Glaciol., 57, 88–102, https://doi.org/10.3189/002214311795306682, 2011.