1. Introduction
A recent, combined reanalysis of satellite-based regional altimetry, interferometry, and gravimetry datasets (Shepherd et al. 2012, hereafter S2012) has estimated the mass balance of ice sheets covering Greenland (GIS), the Antarctic Peninsula (APIS), West Antarctica (WAIS), and East Antarctica (EAIS) over the period 1992–2011 as well as during various epochs within this two-decade time window (Table 1). The uncertainties ascribed to these estimates account for the limited and unique time span associated with each of the three space geodetic techniques; observation errors; and, through modeling, uncertainties in corrections for ice sheet surface mass balance and glacial isostatic adjustment (GIA). Over the period 2000–11, for example, the total estimated mass flux from these ice sheets, −298 ± 58 Gt yr−1 (1-σ uncertainty), equates to a eustatic sea level (ESL) rise of 0.83 ± 0.16 mm yr−1. [We adopt the term “eustatic” as the geographically uniform shift in sea level over the oceans that would yield a volume equal to the meltwater addition. Gregory et al. (2013) have suggested the term “barystatic” for this quantity.] In this paper, we demonstrate that this estimate may be independently tested by invoking the constraint provided by the time rate of change of the zonal harmonic of Earth's gravitational potential at spherical harmonic degree 4
Estimates of secular trends in the low-degree zonal harmonics of the geopotential based on satellite laser ranging (SLR) have been available for over a quarter of a century, initially at spherical harmonic degree 2 (Yoder et al. 1983; Rubincam 1984) and later at higher degrees (Cheng et al. 1989). The first analyses of the
A significant change in the secular trend of the
The spherical harmonic basis function at degree 2 and order zero has the same sign in both high northern and southern latitudes. As such, the
In this regard, the
The
2. Analysis and results
a. GIA predictions
The results in Fig. 1 show the predicted GIA-induced signal in

Components of present-day GIA contributions to
Citation: Journal of Climate 26, 17; 10.1175/JCLI-D-13-00078.1
In addition to the total GIA-induced perturbation in each harmonic (solid lines in Fig. 1), we decompose the perturbation into contributions from the GIS, AIS, and all other ice sheets included in the ICE-5G inventory (e.g., Laurentide, Fennoscandian, etc.). The results indicate that, although ice sheets other than the AIS and GIS dominate the GIA-induced signal in the

Dependence of present-day GIA contributions to
Citation: Journal of Climate 26, 17; 10.1175/JCLI-D-13-00078.1
The decomposition in Fig. 1 highlights another important advantage of using the observed
We begin by focusing on the period 2000–11 and compute a GIA correction based on a combination of our GIA modeling and the W12a Antarctic GIA model (Whitehouse et al. 2012). The results in Fig. 2 indicate that the GIA contribution to

SLR-estimated
Citation: Journal of Climate 26, 17; 10.1175/JCLI-D-13-00078.1
We note that the Antarctic ice volume change in the W12a model (Whitehouse et al. 2012) is significantly smaller than in the ICE-5G history. The difference amounts to approximately 10 m of equivalent eustatic sea level rise, or ~10% of the total non-Antarctic ice volume in ICE-5G. This reduction would have to be compensated by an increase in the excess volume of Northern Hemisphere ice relative to ICE-5G in order to maintain a fit to far-field relative sea level records (Austermann et al. 2013). This suggests that the amplitude of our computed non-Antarctic GIA contribution (Fig. 2) will be biased low by ~10%. Accounting for this effect yields GIA-corrected trends in
b. Polar mass change contribution to 

We next turn to the estimate of polar ice sheet mass balance based on the comprehensive S2012 analysis of satellite-based measurements (Table 1). What signal in
The

Alternatively, in inferring polar ice sheet mass balance, we can avoid entirely the contaminating effect of GIA by considering the change in the trend of the
3. Conclusions
We have demonstrated that estimates of net polar ice sheet mass balance may be independently tested by invoking the long-neglected constraint associated with the rate of change of the zonal harmonic degree-4 of Earth's potential. The
We thank Felix Landerer and two anonymous reviewers for helpful comments and suggestions and Minkang Cheng for providing the time series of SLR zonal harmonics used in the present study. The study was funded by Harvard University, Princeton University, the Canadian Institute for Advanced Research, the Natural Sciences and Engineering Research Council of Canada, and NSF Grant EAR-1014606.
REFERENCES
Austermann, J., , J. X. Mitrovica, , K. Latychev, , and G. A. Milne, 2013: Barbados-based estimate of ice volume at Last Glacial Maximum affected by subducted plate. Nat. Geosci., 6, 553–557, doi:10.1038/ngeo1859.
Benjamin, D., , J. Wahr, , R. D. Ray, , G. D. Egbert, , and S. D. Desai, 2006: Constraints on mantle anelasticity from geodetic observations, and implications for the J2 anomaly. Geophys. J. Int., 165, 3–16.
Cheng, M. K., , and B. D. Tapley, 2004: Variations in the Earth's oblateness during the past 28 years. J. Geophys. Res., 109, B09402, doi:10.1029/2004JB003028.
Cheng, M. K., , R. J. Eanes, , C. K. Shum, , B. E. Schutz, , and B. D. Tapley, 1989: Temporal variations in low degree zonal harmonics from Starlette orbit analysis. Geophys. Res. Lett., 16, 393–396.
Cox, C. M., , and B. F. Chao, 2002: Detection of a large-scale mass redistribution in the terrestrial system since 1998. Science, 297, 831–833.
Dickey, J. O., , S. L. Marcus, , O. de Viron, , and I. Fukumori, 2002: Recent Earth oblateness variations: Unraveling climate and postglacial rebound effects. Science, 298, 1975–1977.
Dziewonski, A. M., , and D. L. Anderson, 1981: Preliminary Reference Earth Model (PREM). Phys. Earth Planet. Inter., 25, 297–356.
Fang, M., , B. H. Hager, , and T. A. Herring, 1996: Surface deformation caused by pressure changes in the fluid core. Geophys. Res. Lett., 23, 1493–1496.
Fleming, K., , and K. Lambeck, 2004: Constraints on the Greenland ice sheet since the Last Glacial Maximum from sea-level observations and glacial-rebound models. Quat. Sci. Rev., 23, 1053–1077.
Gregory, J. M., and Coauthors, 2013: Twentieth-century global-mean sea level rise: Is the whole greater than the sum of the parts? J. Climate, 26, 4476–4499.
Hofmann-Wellenhof, B., , and H. Moritz, 2006: Physical Geodesy. 2nd ed. Springer-Verlag Wien, 403 pp.
Ivins, E. R., , C. G. Sammis, , and C. F. Yoder, 1993: Deep mantle viscous structure with prior estimate and satellite constraint. J. Geophys. Res., 98, 4579–4609.
Jacob, T., , J. Wahr, , W. T. Pfeffer, , and S. Swenson, 2012: Recent contributions of glaciers and ice caps to sea level rise. Nature, 482, 514–518.
James, T. S., , and E. R. Ivins, 1997: Global geodetic signatures of the Antarctic ice sheet. J. Geophys. Res., 102, 605–633.
Kaser, G., , J. G. Cogley, , M. B. Dyurgerov, , M. F. Meier, , and A. Ohmura, 2006: Mass balance of glaciers and ice caps: Consensus estimates for 1961–2004. Geophys. Res. Lett., 33, L19501, doi:10.1029/2006GL027511.
Kendall, R. A., , J. X. Mitrovica, , and G. A. Milne, 2005: On post-glacial sea level: II. Numerical formulation and comparative results on spherically symmetric models. Geophys. J. Int., 161, 679–706.
Lambeck, K., , C. Smither, , and P. Johnston, 1998: Sea-level change, glacial rebound and mantle viscosity for northern Europe. Geophys. J. Int., 134, 102–144.
Marcus, S. L., , J. O. Dickey, , J. K. Willis, , and F. Seitz, 2009: Earth oblateness changes reveal land ice contributions to interannual sea level variability. Geophys. Res. Lett., 36, L23608, doi:10.1029/2009GL041130.
Mitrovica, J. X., , and W. R. Peltier, 1993: Present-day secular variations in the zonal harmonics of the Earth's geopotential. J. Geophys. Res., 98, 4509–4526.
Mitrovica, J. X., , and A. M. Forte, 2004: A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data. Earth Planet. Sci. Lett., 225, 177–189.
Nakada, M., , and K. Lambeck, 1989: Late Pleistocene and Holocene sea-level change in the Australian region and mantle rheology. Geophys. J. Int., 96, 497–517.
Nerem, R. S., , and J. Wahr, 2011: Recent changes in the Earth's oblateness driven by Greenland and Antarctic ice mass loss. Geophys. Res. Lett., 38, L13501, doi:10.1029/2011GL047879.
Peltier, W. R., 1974: The impulse response of a Maxwell Earth. Rev. Geophys., 12, 649–669.
Peltier, W. R., 1983: Constraint on deep mantle viscosity from LAGEOS acceleration data. Nature, 304, 434–436.
Peltier, W. R., 2004: Global glacial isostasy and the surface of the ice-age Earth: The ICE-5G (VM2) model and GRACE. Annu. Rev. Earth Planet. Sci., 32, 111–149.
Rubincam, D. P., 1984: Postglacial rebound observed by LAGEOS and the effective viscosity of the lower mantle. J. Geophys. Res., 89, 1077–1087.
Sabadini, R., , D. A. Yuen, , and P. Gasperini, 1988: Mantle rheology and satellite signatures from present-day glacial forcings. J. Geophys. Res., 93, 437–447.
Shepherd, A., and Coauthors, 2012: A reconciled estimate of ice-sheet mass balance. Science, 338, 1183–1189.
Trupin, A. S., , and R. P. Panfili, 1997: Mass balance of polar ice from long wavelength features of the Earth's gravitational field. Surv. Geophys., 18, 313–326.
Vermeersen, L. L. A., , R. Sabadini, , R. Devoti, , V. Luceri, , P. Rutigliano, , C. Sciarretta, , and G. Bianco, 1998: Mantle viscosity inferences from joint inversions of Pleistocene deglaciation-induced changes in geopotential with a new SLR analysis and polar wander. Geophys. Res. Lett., 25, 4261–4264.
Whitehouse, P. L., , M. J. Bentley, , G. A. Milne, , M. A. King, , and I. D. Thomas, 2012: A new glacial isostatic adjustment model for Antarctica: Calibrated and tested using observations of relative sea-level change and present-day uplift rates. Geophys. J. Int., 190, 1464–1482.
Wu, P., , and W. R. Peltier, 1984: Pleistocene deglaciation and the Earth's rotation: A new analysis. Geophys. J. Roy. Astron. Soc., 76, 753–791.
Yoder, C. F., , J. G. Williams, , J. O. Dickey, , B. E. Schutz, , R. J. Eanes, , and B. D. Tapley, 1983: Secular variation of Earth's gravitational harmonic J2 coefficient from LAGEOS and non-tidal acceleration of earth rotation. Nature, 303, 757–762.
Yuen, D. A., , and R. Sabadini, 1985: Viscosity stratification of the lower mantle as inferred from the J2 observation. Ann. Geophys., 3, 647–654.