Anny Cazenave
email: anny.cazenave@legos.obs-mip.fr

Global mean sea level is an integral of changes occurring in the climate system in response to unforced climate variability as well as natural and anthropogenic forcing factors. Its temporal evolution allows changes (e.g., acceleration) to be detected in one or more components. Study of the sea-level budget provides constraints on missing or poorly known contributions, such as the unsurveyed deep ocean or the still uncertain land water component. In the context of the World Climate Research Programme Grand Challenge entitled "Regional Sea Level and Coastal Impacts", an international effort involving the sea-level community worldwide has been recently initiated with the objective of assessing the various datasets used to estimate components of the sea-level budget during the altimetry era (1993 to present). These datasets are based on the combination of a broad range of space-based and in situ observations, model estimates, and algorithms. Evaluating their quality, quantifying uncertainties and identifying sources of discrepancies between component estimates is extremely useful for various applications in climate research. This effort involves several tens of scientists from about 50 research teams/institutions worldwide (www.wcrp-climate.org/grand-challenges/gc-sea-level, last access: 22 August 2018). The results presented in this paper are a synthesis of the first assessment performed during 2017-2018. We present estimates of the altimetry-based global mean sea level (average rate of 3.1 ± 0.3 mm yr -1 and acceleration of 0.1 mm yr -2 over 1993-present), as well as of the different components of the sea-level budget (http://doi.org/10.17882/54854, last access: 22 August 2018). We further examine closure of the sea-level budget, comparing the observed global mean sea level with the sum of components. Ocean thermal expansion, glaciers, Greenland and Antarctica contribute 42 %, 21 %, 15 % and 8 % to the global mean sea level over the 1993-present period. We also study the sea-level budget over 2005-present, using GRACE-based ocean mass estimates instead of the sum of individual mass components. Our results demonstrate that the global mean sea level can be closed to within 0.3 mm yr -1 (1σ ). Substantial uncertainty remains for the land water storage component, as shown when examining individual mass contributions to sea level.

Introduction

Global warming has already several visible consequences, in particular an increase in the Earth's mean surface temperature and ocean heat content [START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]IPCC, 2013), melting of sea ice, loss of mass of glaciers [START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF], and ice mass loss from the Greenland and Antarctica ice sheets (Rignot et al., 2011a;[START_REF] Shepherd | A reconciled estimate of ice-sheet mass balance[END_REF]. On average over the last 50 years, about 93 % of heat excess accumulated in the climate system because of greenhouse gas emissions has been stored in the ocean, and the remaining 7 % has been warming the atmosphere and con-tinents, and melting sea and land ice [START_REF] Von Schuckmann | Earth's energy imbalance: an imperative for monitoring[END_REF]. Because of ocean warming and land ice mass loss, sea level rises. Since the end of the last deglaciation about 3000 years ago, sea level remained nearly constant (e.g., [START_REF] Lambeck | Sea-level change from mid-Holocene to recent time: An Australian example with global implications, in: Ice Sheets, Sea Level and the Dynamic Earth[END_REF][START_REF] Lambeck | Paleoenvironmental records, geophysical modelling and reconstruction of sea level trends and variability on centennial and longer time scales[END_REF][START_REF] Kemp | Climate related sea level variations over the past two millennia[END_REF]. However, direct observations from in situ tide gauges available since the mid-to-late 19th century show that the 20th century global mean sea level has started to rise again at a rate of 1.2 to 1.9 mm yr -1 (Church and White, 2011;[START_REF] Jevrejeva | Trends and Acceleration in Global and Regional Sea Levels since 1807[END_REF]Hay et al., 2015;[START_REF] Dangendorf | Reassessment of 20th Century Global Mean Sea Level Rise[END_REF]. Since the early 1990s sea-level rise (SLR) is measured by high-precision altimeter satellites and the rate has Published by Copernicus Publications. increased to ∼ 3 mm yr -1 on average (Legeais et al., 2018;[START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF].

Accurate assessment of present-day global mean sea-level variations and its components (ocean thermal expansion, ice sheet mass loss, glaciers mass change, changes in land water storage, etc.) is important for many reasons. The global mean sea level is an integral of changes occurring in the Earth's climate system in response to unforced climate variability as well as natural and anthropogenic forcing factors, e.g., net contribution of ocean warming, land ice mass loss and changes in water storage in continental river basins. Temporal changes in the components are directly reflected in the global mean sea-level curve. If accurate enough, study of the sea-level budget provides constraints on missing or poorly known contributions, e.g., the deep ocean undersampled by current observing systems, or still uncertain changes in water storage on land due to human activities (e.g., groundwater depletion in aquifers). Global mean sea level corrected for ocean mass change in principle allows one to independently estimate temporal changes in total ocean heat content, from which the Earth's energy imbalance can be deduced (von Schuckmann et al., 2016). The sea level and/or ocean mass budget approach can also be used to constrain models of glacial isostatic adjustment (GIA). The GIA phenomenon has a significant impact on the interpretation of GRACE-based space gravimetry data over the oceans (for ocean mass change) and over Antarctica (for ice sheet mass balance). However, there is still no complete consensus on best estimates, a result of uncertainties in deglaciation models and mantle viscosity structure. Finally, observed changes in the global mean sea level and its components are fundamental for validating climate models used for projections.

In the context of the Grand Challenge entitled "Regional Sea Level and Coastal Impacts" of the World Climate Research Programme (WCRP), an international effort involving the sea-level community worldwide has been recently initiated with the objective of assessing the sea-level budget during the altimetry era (1993 to present). To estimate the different components of the sea-level budget, different datasets are used. These are based on the combination of a broad range of space-based and in situ observations. Evaluating their quality, quantifying their uncertainties and identifying the sources of discrepancies between component estimates, including the altimetry-based sea-level time series, are extremely useful for various applications in climate research.

Several previous studies have addressed the sea-level budget over different time spans and using different datasets. For example, [START_REF] Munk | Twentieth century sea level: An enigma[END_REF] found that the 20th century sea-level rise could not be closed with the data available at that time and showed that if the missing contribution were due to polar ice melt, this would be in conflict with external astronomical constraints. The enigma has been resolved in two ways. Firstly, an improved theory of rotational stability of the Earth [START_REF] Mitrovica | Reanalysis of ancient eclipse, astronomic and geodetic data: A possible route to resolving the enigma of global sea-level rise[END_REF] effectively removed the constraints proposed by [START_REF] Munk | Twentieth century sea level: An enigma[END_REF] and allows a polar ice sheet contribution to 20th century sea-level rise of as much as ∼ 1.1 mm yr -1 , with about 0.8 mm yr -1 beginning in the 20th century. In addition, more recent studies by Gregory et al. (2013) and Slangen et al. (2017), combining observations with model estimates, showed that it was possible to effectively close the 20th century sea-level budget within uncertainties, particularly over the altimetry era (e.g., [START_REF] Cazenave | Sea level budget over 2003-2008: A reevaluation from GRACE space gravimetry, satellite altimetry and Argo[END_REF][START_REF] Leuliette | Balancing the sea level budget[END_REF]Church and White, 2011;Llovel et al., 2014;Chambers et al., 2017;Dieng et al., 2017;X. Chen et al., 2017;[START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF]. Assessments of the published literature have also been performed in past IPCC (Intergovernmental Panel on Climate Change) reports (e.g., [START_REF] Church | Sea level change, in: Climate Change 2013: The Physical Science Basis[END_REF]. Building on these previous works, here we intend to provide a collective update of the global mean sea-level budget, involving the many groups worldwide interested in present-day sea-level rise and its components. We focus on observations rather than model-based estimates and consider the high-precision altimetry era starting in 1993. This era includes the period since the mid-2000s in which new observing systems, like the Argo float project (Roemmich et al., 2012) and the GRACE space gravimetry mission (Tapley et al., 2004a, b), provide improved datasets of high value for such a study. Only the global mean budget is considered here. Regional budget will be the focus of a future assessment.

Section 2 describes for each component of the sea-level budget equation the different datasets used to estimate the corresponding contribution to sea level, discusses associated errors and provides trend estimates for the two periods. Section 3 addresses the mass and sea-level budgets over the study periods. A discussion is provided in Sect. 4, followed by a conclusion.

Methods and data

In this section, we briefly present the global mean sea-level budget (Sect. 2.1) and then provide, for each term of the budget equation, an assessment of the most up-to-date published results. Multiple organizations and research groups routinely generate the basic measurements as well as the derived datasets and products used to study the sea-level budget. Sections 2.2 to 2.7 summarize the measurements and methodologies used to derive observed sea level, as well as steric and mass components. In most cases, we focus on observations but in some instances (e.g., for GIA corrections applied to the data), model-based estimates are the only available information.

Sea-level budget equation

Global mean sea level (GMSL) change as a function of time t is usually expressed by the sea-level budget equation: GMSL(t) = GMSL(t) steric + GMSL(t) ocean mass ,

(1)

where GMSL(t) steric refers to the contributions of ocean thermal expansion and salinity to sea-level change, and GMSL(t) oceanmass refers to the change in mass of the oceans. Due to water conservation in the climate system, the ocean mass term (also noted as M(t) ocean ) can further be expressed as follows:

M(t) ocean + M(t) glaciers + M(t) Greenland + M(t) Antarctica + M(t) TWS + M(t) WV + M(t) Snow

+ uncertainty = 0, (2)
where M(t) glaciers , M(t) Greenland , M(t) Antarctica , M(t) TWS , M(t) WV and M(t) Snow represent temporal changes in mass of glaciers, Greenland and Antarctica ice sheets, terrestrial water storage (TWS), atmospheric water vapor (WV), and snow mass changes. The uncertainty is a result of uncertainties in all of the estimates. For the altimetry era, many studies have investigated closure of the sea-level budget and potentially missing mass terms, for example, permafrost melting. From Eq. ( 2), we deduce the following: GMSL(t) ocean mass = -[M(t) glaciers + M(t) Greenland + M(t) Antarctica + M(t) TWS + M(t) WV + M(t) Snow

+ missing mass terms]

(3)

In the next subsections, we successively discuss the different terms of the budget (Eqs. 1 and 2) and how they are estimated from observations. We do not consider the atmospheric water vapor and snow components, assumed to be small. Two periods are considered: (1) 1993-present (i.e., the entire altimetry era) and (2) 2005-present (i.e., the period covered by both Argo and GRACE).

Altimetry-based global mean sea level over 1993-present

The launch of the TOPEX/Poseidon (T/P) altimeter satellite in 1992 led to a new paradigm for measuring sea level from space, providing for the first time precise and globally distributed sea-level measurements at 10-day intervals.

At the time of the launch of T/P, the measurements were not expected to have sufficient accuracy for measuring GMSL changes. However, as the radial orbit error decreased from ∼ 10 cm at launch to ∼ 1 cm presently, and other instrumental and geophysical corrections applied to altimetry system improved (e.g., [START_REF] Stammer | Satellite Altimetry Over Oceans and Land Surfaces[END_REF], several groups regularly provided an altimetry-based GMSL time series (e.g., [START_REF] Nerem | Estimating Mean Sea Level Change from the TOPEX and Jason Altimeter Missions[END_REF]Church et al., 2011;Ablain et al., 2015;Legeais et al., 2018). The initial T/P GMSL time series was extended with the launch of Jason-1 (2001), Jason-2 (2008) and Jason-3 (2016). By design, each of these missions has an overlap period with the previous one in order to intercompare the sea-level measurements and estimate instrument biases (e.g., [START_REF] Nerem | Estimating Mean Sea Level Change from the TOPEX and Jason Altimeter Missions[END_REF]Ablain et al., 2015). This has allowed the construction of an uninterrupted GMSL time series that is currently 25 years long.

Global mean sea-level datasets

Six groups (AVISO/CNES, SL_cci/ESA, University of Colorado, CSIRO, NASA/GSFC, NOAA) provide altimetrybased GMSL time series. All of them use 1 Hz altimetry measurements derived from T/P, Jason-1, Jason-2 and Jason-3 as reference missions. These missions provide the most accurate long-term stability at global and regional scales [START_REF] Ablain | A New Assessment of the Error Budget of Global Mean Sea Level Rate Estimated by Satellite Altimetry over 1993-2008[END_REF](Ablain et al., , 2017a)), and are all on the same historical T/P ground track. This allows computation of a long-term record of the GMSL from 1993 to present. In addition, complementary missions (ERS-1, ERS-2, Envisat, Geosat Followon, CryoSat-2, SARAL/AltiKa and Sentinel-3A) provide increased spatial resolution and coverage of high-latitude ocean areas, pole-ward of 66 • N-S latitude (e.g., the European Space Agency/ESA Climate Change Initiative/CCI sealevel dataset; Legeais et al., 2018).

The above groups adopt different approaches when processing satellite altimetry data. The most important differences concern the geophysical corrections needed to account for various physical phenomena such as atmospheric propagation delays, sea state bias, ocean tides, and the ocean response to atmospheric wind and pressure forcing. Other differences come from data editing, methods to spatially average individual measurements during orbital cycles and links between successive missions [START_REF] Masters | Comparison of Global Mean Sea Level Time Series from TOPEX/Poseidon, Jason-1, and Jason-2[END_REF]Henry et al., 2014).

Overall, the quality of the different GMSL time series is similar. Long-term trends agree well to within 6 % of the signal, approximately 0.2 mm yr -1 (see Fig. 1) within the GMSL trend uncertainty range (∼ 0.3 mm yr -1 ; see next section). The largest differences are observed at interannual timescales and during the first years (before 1999; see below). Here we use an ensemble mean GMSL based on averaging all individual GMSL time series.

Global mean sea-level uncertainties and TOPEX-A drift

Based on an assessment of all sources or uncertainties affecting satellite altimetry (Ablain et al., 2017a), the GMSL trend uncertainty (90 % confidence interval) is estimated as 0.3 to 0.4 mm yr -1 over the whole altimetry era .

The main contribution to the uncertainty is the wet tropospheric correction with a drift uncertainty in the range of 0.2-0.3 mm yr -1 (Legeais et al., 2018) over a 10-year period. To a lesser extent, the orbit error [START_REF] Couhert | Towards the 1mm/y Stability of the Radial Orbit Error at Regional Scales[END_REF] et al., 2018) and the altimeter parameters' (range, σ 0 and significant wave height -SWH) instability [START_REF] Ablain | Detection of Long-Term Instabilities on Altimeter Backscatter Coefficient Thanks to Wind Speed Data Comparisons from Altimeters and Models[END_REF] also contribute to the GMSL trend uncertainty, at the level of 0.1 mm yr -1 . Furthermore, imperfect links between successive altimetry missions lead to another trend uncertainty of about 0.15 mm yr -1 over the 1993-2017 period [START_REF] Zawadzki | Estimating a drift in TOPEX-A Global Mean Sea Level using Poseidon-1 measurements[END_REF]. Uncertainties are higher during the first decade (1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002), when T/P measurements display larger errors at climatic scales. For instance, the orbit solutions are much more uncertain due to gravity field solutions calculated without GRACE data. Furthermore, the switch from TOPEX-A to TOPEX-B in February 1999 (with no overlap between the two instrumental observations) leads to an error of ∼ 3 mm in the GMSL time series [START_REF] Escudier | Satellite radar altimetry: principle, accuracy and precision[END_REF].

However, the most significant error that affects the first 6 years (January 1993 to February 1999) of the T/P GMSL measurements is due to an instrumental drift of the TOPEX-A altimeter, not included in the formal uncertainty estimates discussed above. This effect on the GMSL time series was recently highlighted via comparisons with tide gauges [START_REF] Valladeau | Comparing Altimetry with Tide Gauges and Argo Profiling Floats for Data Quality Assessment and Mean Sea Level Studies[END_REF]Watson et al., 2015;X. Chen et al., 2017;[START_REF] Ablain | The TOPEX-A Drift and Impacts on GMSL Time Series[END_REF], via a sea-level budget approach (i.e., comparison with the sum of mass and steric components; Dieng et al., 2017) and by comparing with Poseidon-1 measurements (Lionel Zawadsky, personal communication, 2017). In a recent study, Beckley et al. (2017) asserted that the corresponding error on the 1993-1998 GMSL resulted from incorrect onboard calibration parameters.

All approaches conclude that during the period January 1993 to February 1999, the altimetry-based GMSL was overestimated. TOPEX-A drift correction was estimated to be close to 1.5 mm yr -1 (in terms of sea-level trend) with an uncertainty of ±0.5 to ±1.0 mm yr -1 (Watson et al., 2015;X. Chen et al., 2017;Dieng et al., 2017). Beckley et al. (2017) proposed to not apply the suspect onboard calibration correction on TOPEX-A measurements. The impact of this approach is similar to the TOPEX-A drift correction estimated by Dieng et al. (2017) and [START_REF] Ablain | The TOPEX-A Drift and Impacts on GMSL Time Series[END_REF]. In the latter study, accurate comparison between TOPEX-A-based GMSL and tide gauge measurements leads to a drift correction of about -1.0 mm yr -1 between January 1993 and July 1995, and +3.0 mm yr -1 between August 1995 and February 1999, with an uncertainty of 1.0 mm yr -1 (with a 68 % confidence level, see Table 1).

Global mean sea-level variations

The ensemble mean GMSL rate after correcting for the TOPEX-A drift (for all of the proposed corrections) amounts to 3.1 mm yr -1 over 1993-2017 (Fig. 2). This corresponds to a mean sea-level rise of about 7.5 cm over the whole altimetry period. More importantly, the GMSL curve shows a net acceleration, estimated to be at 0.08 mm yr -2 (X. Chen et al., 2017;Dieng et al., 2017) and0 Removing the trend from the GMSL time series highlights interannual variations (not shown). Their magnitudes depend on the period (+3 mm in 1998-1999, -5 mm in 2011-2012 and +10 mm in 2015-2016) and are well correlated in time with El Niño and La Niña events [START_REF] Nerem | Estimating Mean Sea Level Change from the TOPEX and Jason Altimeter Missions[END_REF][START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF]Cazenave et al., 2014). However, substantial differences (of 1-3 mm) exist between the six detrended GMSL time series. This issue needs further investigation.

For the sea-level budget assessment (Sect. 3), we will use the ensemble mean GMSL time series corrected for the TOPEX-A drift using the [START_REF] Ablain | The TOPEX-A Drift and Impacts on GMSL Time Series[END_REF] correction.

Comparison with tide gauges

Prior to 1992, global sea-level rise estimates relied on the tide gauge measurements, and it is worth mentioning past attempts to produce global sea-level reconstructions utilizing these measurements (e.g., [START_REF] Gornitz | Global sea level trend in the past century[END_REF][START_REF] Bartnett | The estimation of "global" sea level change: A problem of uniqueness[END_REF][START_REF] Douglas | Global sea level rise[END_REF][START_REF] Douglas | Global sea rise: a redetermination[END_REF]Douglas, , 2001). Here we focus on global sea-level reconstructions that overlap with satellite altimetry data over a substantial common time span. Some of these reconstructions rely on tide gauge data only [START_REF] Jevrejeva | Nonlinear trends and multi-year cycle in sea level records[END_REF][START_REF] Jevrejeva | Trends and Acceleration in Global and Regional Sea Levels since 1807[END_REF][START_REF] Merrifield | An anomalous recent acceleration of global sea level rise[END_REF][START_REF] Wenzel | Reconstruction of regional mean sea level anomalies from tide gauges using neural networks[END_REF][START_REF] Ray | Experiments in reconstructing twentieth-century sea levels[END_REF]Hamlington et al., 2011;[START_REF] Spada | New estimates of secular sea level rise from tide gauge data and GIA modelling[END_REF][START_REF] Thompson | A unique asymmetry in the pattern of recent sea level change[END_REF][START_REF] Dangendorf | Reassessment of 20th Century Global Mean Sea Level Rise[END_REF][START_REF] Frederikse | A consistent sea-level reconstruction and its budget on basin and global scales over 1958-2014[END_REF]. In addition, there are reconstructions that jointly use satellite altimetry, tide gauge records (Church andWhite, 2006, 2011) and reconstructions, which combine tide gauge records with ocean models (Meyssignac et al., 2011) or physics-based and model-derived geometries of the contributing processes (Hay et al., 2015).

For the period since 1993, with most of the world coastlines densely sampled, the rates of sea-level rise from all tide-gauge-based reconstructions and estimates from satellite altimetry agree within their specific uncertainties, e.g., rates of 3.0 ± 0.7 mm yr -1 (Hay et al. 2015), 2.8 ± 0.5 mm yr -1 (Church and White, 2011;[START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF], 3.1±0.6 mm yr -1 [START_REF] Jevrejeva | Trends and Acceleration in Global and Regional Sea Levels since 1807[END_REF], 3.1±1.4 mm yr -1 [START_REF] Dangendorf | Reassessment of 20th Century Global Mean Sea Level Rise[END_REF] and the estimate from satellite altimetry 3.2 ± 0.4 mm yr -1 [START_REF] Nerem | Estimating Mean Sea Level Change from the TOPEX and Jason Altimeter Missions[END_REF][START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]. However, classical tide-gauge-based reconstructions still tend to overestimate the interannual to decadal variability of global mean sea level (e.g., [START_REF] Calafat | On the ability of global sea level reconstructions to determine trends and variability[END_REF]Dangendorf et al., 2015;[START_REF] Natarov | Regional influences on reconstructed global mean sea level[END_REF] [START_REF] Frederikse | A consistent sea-level reconstruction and its budget on basin and global scales over 1958-2014[END_REF]. Uncertainties involved in specifying vertical land motion corrections at tide gauges also impact tide gauge reconstructions [START_REF] Jevrejeva | Trends and Acceleration in Global and Regional Sea Levels since 1807[END_REF]Wöppelmann and Marcos, 2016;Hamlington et al., 2016). [START_REF] Frederikse | A consistent sea-level reconstruction and its budget on basin and global scales over 1958-2014[END_REF] also recently demonstrated that both global mean sea level reconstructed from tide gauges and the sum of steric and mass contributors show a good agreement with altimetry estimates for the overlapping period 1993-2014.

Steric sea level

Steric sea-level variations result from temperature-(T ) and salinity-(S) related density changes in sea water associated with volume expansion and contraction. These are referred to as thermosteric and halosteric components. Despite clear detection of regional salinity changes and the dominance of the salinity effect on density changes at high latitudes [START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF], the halosteric contribution to presentday global mean steric sea-level rise is negligible, as the ocean's total salt content is essentially constant over multidecadal timescales [START_REF] Gregory | Predictions of global and regional sea-level rise using AOGCMs with and without flux adjustment[END_REF]. Hence, in this study, we essentially consider the thermosteric sea-level component.

Averaged over the 20th century, ocean thermal expansion associated with ocean warming has been the largest contribution to global mean sea-level rise [START_REF] Church | Sea level change, in: Climate Change 2013: The Physical Science Basis[END_REF]. This remains true for the altimetry period starting in the year 1993 (e.g., X. Chen et al., 2017;Dieng et al., 2017;[START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF]. But total land ice mass loss (from glaciers, Greenland and Antarctica) during this period now dominates the sea-level budget (see Sect. 3).

Until the mid-2000s, the majority of ocean temperature data were retrieved from shipboard measurements. These include vertical temperature profiles along research cruise tracks from the surface sometimes all the way down to the bottom layer (e.g., [START_REF] Purkey | Warming of global abyssal and deep southern ocean waters between the 1990s and 2000s: Contributions to global heat and sea level rise budget[END_REF] and upperocean broad-scale measurements from ships of opportunity (Abraham et al., 2013). These upper-ocean in situ temperature measurements, however, are limited to the upper 700 m depth due to common use of expandable bathythermographs (XBTs). Although the coverage has been improved through time, large regions characterized by difficult meteorological conditions remained under-sampled, in particular the southern hemispheric oceans and the Arctic area.

Thermosteric datasets

Over the altimetry era, several research groups have produced gridded time series of temperature data for different depth levels, based on XBTs (with additional data from mechanical bathythermographs -MBTs -and conductivitytemperature-depth -CTD -devices and moorings) and Argo float measurements. The temperature data have further been used to provide thermosteric sea-level products. These differ because of different strategies adopted for data editing, temporal and spatial data gaps filling, mapping methods, baseline climatology, and instrument bias corrections (in particular the time-to-depth correction for XBT data, [START_REF] Boyer | Sensitivity of global ocean heat content estimates to mapping methods, XBT bias cor-rections, and baseline climatology[END_REF].

The global ocean in situ observing system has been dramatically improved through the implementation of the international Argo program of autonomous floats, delivering a unique insight into the interior ocean from the surface down to 2000 m depth of the ice-free global ocean (Roemmich et al., 2012;[START_REF] Riser | Fifteen years of ocean observations with the global Argo array[END_REF]. More than 80 % of initially planned full deployment of Argo float program was achieved during the year 2005, with quasi global coverage of the icefree ocean by the start of 2006. At present, more than 3800 floats provide systematic T and S data, with quasi (60 • S-60 • N latitude) global coverage down to 2000 m depth. A full overview on in situ ocean temperature measurements is given for example in Abraham et al. (2013).

In this section, we consider a set of 11 direct (in situ) estimates, publicly available over the entire altimetry era, to review global mean thermosteric sea-level rise and, ultimately, to construct an ensemble mean time series. These datasets are as follows:

1 Their characteristics are presented in Table 2.

Individual estimates

All in situ estimates compiled in this study show a steady rise in global mean thermosteric sea level, independent of depth integration and decadal or multidecadal periods (Figs. 4 and Interannual to decadal variability during the altimeter era (since 1993) is similar for both 0-700 and 700-2000 m, with larger amplitude in the upper ocean (Figs. 4 and5, right panels). For the 0-700 m, there is an apparent change in amplitude before and after the Argo era (since 2005), mostly due to a maximum (2-4 mm) around 2001-2004, except for one estimate. Higher amplitude and larger spread in variability between estimates before the Argo era is a symptom of the much sparser in situ coverage of the global ocean. Interannual variability over the Argo era (Figs. 4 and5, right panels) is mainly modulated by El Niño-Southern Oscillation (ENSO) phases in the upper 500 m of the ocean, particularly for the Pacific, the largest ocean basin (Roemmich and Gilson, 2011;[START_REF] Roemmich | Multidecadal change of the South Pacific gyre circulation[END_REF]Johnson and Birnbaum, 2017).

In terms of depth contribution, on average, the upper 300 m explains the same percentage (almost 70 %) of the 0-700 m linear rate over both altimetry and Argo eras, but the contribution from the 0-700 to 0-2000 m varies: about 75 % for 1993-2016 and 65 % for 2005-2016. Thus, the 700-2000 m contribution increases by 10 % during the Argo decade, when the number of observations within 700-2000 m has significantly increased.

Ensemble mean thermosteric sea level

Given that the global mean thermosteric sea-level anomaly estimates compiled for this study are not necessarily referenced to the same baseline climatology, they cannot be directly averaged together to create an ensemble mean. To circumvent this limitation, we created an ensemble mean in three steps, as explained below.

Firstly, we detrended the individual time series by removing a linear trend for 1993-2016 and averaged together to obtain an "ensemble mean variability time series". Sec- ondly, we averaged together the corresponding linear trends of the individual estimates to obtain an "ensemble mean linear rate". Thirdly, we combined this "ensemble mean linear rate" with the "ensemble mean variability time series" to obtain the final ensemble mean time series. We applied the same steps for the Argo era (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016).

To maximize the number of individual estimates used in the final full-depth ensemble mean time series, the three steps above were actually divided into depth integrations and then summed. For the Argo era, we summed 0-2000 m (nine estimates) and ≥ 2000 m (one estimate). For the altimetry era, we summed 0-700 m (six estimates), 700-2000 (four esti-mates) and ≥ 2000 m (one estimate), although there is no statistical difference if the calculation was only based on the sum of 0-2000 m (4 estimates) and ≥ 2000 m (1 estimate). There is also no statistical difference between the full-depth ensemble mean time series created for the Altimeter and Argo eras during their overlapping years (since 2005).

Figure 6 shows the full-depth ensemble mean time series over 1993-2015 and 2005-2015. It reveals a global mean thermosteric sea-level rise of about 30 mm over 1993-2016 (24 years) or about 18 mm over 2005-2016 (12 years), with a record high in 2015. These thermosteric changes are equiv- alent to a linear rate of 1.32 ± 0.4 and 1.31 ± 0.4 mm yr -1 respectively.

Figure 7 shows thermosteric sea-level trends for each of the datasets used over

Glaciers

Glaciers have strongly contributed to sea-level rise during the 20th century -around 40 % -and will continue to be an important part of the projected sea-level change during the 21st century -around 30 % (Kaser et al., 2006;Church et al., 2013;[START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF][START_REF] Marzeion | Attribution of global glacier mass loss to anthropogenic and natural causes[END_REF][START_REF] Zemp | Historically unprecedented global glacier decline in the early 21st century[END_REF][START_REF] Huss | A new model for global glacier change and sea-level rise[END_REF]. Because glaciers are time-integrated dynamic systems, a response lag of at least 10 years to a few hundred years is observed between changes in climate forcing and glacier shape, mainly depending on glacier length and slope [START_REF] Johannesson | Time-Scale for Adjustment of Glaciers to Changes in Mass Balance[END_REF][START_REF] Bahr | Response time of glaciers as a function of size and mass balance: 1. Theory[END_REF]. Today, glaciers are globally (a notable exception is the Karakoram-Kunlun Shan region, e.g., Brun et al., 2017) in a strong disequilibrium with the current climate and are losing mass, due essentially to the global warming in the second half of the 20th century (Marzeion et al., 2018). Global glacier mass changes are derived from in situ measurements of glacier mass changes or glacier length changes. Remote sensing methods measure elevation changes over entire glaciers based on differencing digital elevation models (DEMs) from satellite imagery between two epochs (or at points from repeat altimetry), surface flow velocities for determination of mass fluxes and glacier mass changes from space-based gravimetry. Mass balance modeling driven by climate observations is also used [START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], provide a review of these different methods).

Glacier contribution to sea level is primarily the result of their surface mass balance and dynamic adjustment, plus iceberg discharge and frontal ablation (below sea level) in the case of marine-terminating glaciers. The sum of worldwide glacier mass balances does not correspond to the total glacier contribution to sea-level change for the following reasons:

-Glacier ice below sea level does not contribute to sealevel change, apart from a small lowering when replacing ice with seawater of a higher density. Total volume of glacier ice below sea level is estimated to be 10-60 mm sea-level equivalent (SLE, Huss and Farinotti, 2012;Haeberli and Linsbauer, 2013;[START_REF] Huss | A new model for global glacier change and sea-level rise[END_REF].

-There is incomplete transfer of melting ice from glaciers to the ocean: meltwater stored in lakes or wetlands, meltwater intercepted by natural processes and human activities (e.g., drainage to lakes and aquifers in endorheic basins, impoundment in reservoirs, agriculture use of freshwater, [START_REF] Loriaux | Evolution of glacial lakes from the Northern Patagonia Icefield and terrestrial water storage in a sealevel rise context, Global Planet[END_REF]Käab et al., 2015).

Despite considerable progress in observing methods and spatial coverage [START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], estimating glacier con-tribution to sea-level change remains challenging due to the following reasons:

-The number of regularly observed glaciers (in the field) remains very low (0.25 % of the 200 000 glaciers of the world have at least one observation and only 37 glaciers have multidecade-long observations, [START_REF] Zemp | Historically unprecedented global glacier decline in the early 21st century[END_REF].

-Uncertainty of the total glacier ice mass remains high (Fig. 8, Grinsted, 2013;[START_REF] Pfeffer | The Randolph Glacier Inventory: a globally complete inventory of glaciers[END_REF]Farinotti et al., 2017;Frey et al. 2014).

-Uncertainties in glacier inventories and DEMs are not negligible. Sources of uncertainties include debriscovered glaciers, disappearance of small glaciers, positional uncertainties, wrongly mapped seasonal snow, rock glaciers, voids and artifacts in DEMs [START_REF] Paul | Combining satellite multispectral image data and a digital elevation model for mapping of debriscovered glaciers[END_REF]Bahr and Radić, 2012).

-Uncertainties of satellite retrieval algorithms from space-based gravimetry and regional DEM differencing are still high, especially for global estimates [START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF][START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF]Chambers et al., 2017).

-Uncertainties of global glacier modeling (e.g., initial conditions, model assumptions and simplifications, local climate conditions; Marzeion et al., 2012).

-Knowledge about some processes governing mass balance (e.g., wind redistribution and metamorphism, sublimation, refreezing, basal melting) and dynamic processes (e.g., basal hydrology, fracking, surging) remains limited (Farinotti et al., 2017).

An annual assessment of glacier contribution to sea-level change is difficult to perform from ground-based or spacebased observations apart from space-based gravimetry, due to the sparse and irregular observation of glaciers, and the difficulty of accurately assessing the annual mass balance variability. Global annual averages are highly uncertain because of the sparse coverage, but successive annual balances are uncorrelated and therefore averages over several years are known with greater confidence.

Glacier datasets

The following datasets are considered, with a focus on the trends of annual mass changes:

1. update of [START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF] [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF], from satellite gravimetry and altimetry, and glaciological records, called G16;

2. update of Marzeion et al. (2012) [START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], from global glacier modeling and mass balance observations, called M17;

3. update of [START_REF] Cogley | Geodetic and direct mass-balance measurements: comparison and joint analysis[END_REF] [START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], from geodetic and direct mass-balance measurements, called C17;

4. update of [START_REF] Leclercq | Estimating the glacier contribution to sea-level rise for the period 1800-2005[END_REF][START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], from glacier length changes, called L17;

5. average of GRACE-based estimates of [START_REF] Marzeion | Observation-Based Estimates of Global Glacier Mass Change and Its Contribution to Sea-Level Change[END_REF], from spatial gravimetry measurements, called M17-G.

In general it is not possible to align measurements of glacier mass balance with the calendar. Most in situ measurements are for glaciological years that extend between successive annual minima of the glacier mass at the end of the summer melt season. Geodetic measurements have start and end dates several years apart and are distributed irregularly through the calendar year; some are corrected to align with annual mass minima but most are not. Consequently, measurements discussed here for 1993-2016 (the altimetry era) and 2005-2016 (the GRACE and Argo era) are offset by up to a few months from the nominal calendar years.

Peripheral glaciers around the Greenland and Antarctic ice sheets are not treated in detail in this section (see Sects. 2.5 and 2.6 for mass-change estimates that combine the peripheral glaciers with the Greenland ice sheet and Antarctic ice sheet respectively). This is primarily because of the lack of observations (especially ground-based measurements) and also because of the high spatial variability of mass balance in those regions, and the slightly different climate (e.g., precipitation regime) and processes (e.g., refreezing). In the past, these regions have often been neglected. However, Radić and Hock (2010) estimated the total ice mass of peripheral glaciers around Greenland and Antarctica as 191 ± 70 mm SLE, with an actual contribution to sea-level rise of around 0.23 ± 0.04 mm yr -1 (Radić and Hock, 2011). [START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF] found a contribution from Greenland and Antarctic peripheral glaciers equal to 0.12 ± 0.05 mm yr -1 .

Note that some new or updated datasets for peripheral glaciers surrounding polar ice sheets are under development and will hopefully be available in coming years in order to incorporate Greenland and Antarctic peripheral glaciers in the estimates of global glacier mass changes.

Methods

No globally complete observational dataset exists for glacier mass changes (except GRACE estimates; see below). Any calculation of the global glacier contribution to sea-level change has to rely on spatial interpolation or extrapolation or both, or to consider limited knowledge of responses to climate change (due to the heterogeneous spatial distribution of glaciers around the world). Consequently, most observational methods to derive glacier sea-level contribution must extend local observations (in situ or satellite) to a larger region. Thanks to the recent global glacier outline inventory (Randolph Glacier Inventory -RGI -first release in 2012) as well as global climate observations, glacier modeling can now also be used to estimate the contribution of glaciers to sea level (Marzeion et al., 2012;[START_REF] Huss | A new model for global glacier change and sea-level rise[END_REF][START_REF] Maussion | The Open Global Glacier Model (OGGM) v1.0[END_REF]. Still, those global modeling methods need to globalize local observations and glacier processes which require fundamental assumptions and simplifications. Only GRACE-based gravimetric estimates are global but they suffer from large uncertainties in retrieval algorithms (signal leakage from hydrology, GIA correction) and coarse spatial resolution, not resolving smaller glaciated mountain ranges or those peripheral to the Greenland ice sheet.

The DEM differencing method is not yet global, but regional, and can hopefully in the near future be applied globally. This method needs also to convert elevation changes to mass changes (using assumptions on snow and ice densities). In contrast, very detailed glacier surface mass balance and glacier dynamic models are today far from being applicable globally, mainly due to the lack of crucial observations (e.g., meteorological data, glacier surface velocity and thickness) and of computational power for the more demanding theoretical models. However, somewhat simplified approaches are currently being developed to make the best use of the steadily increasing datasets. Modeling-based estimates suffer also from the large spread in estimates of the actual global glacier ice mass (Fig. 8). The mean value is 469±146 mm SLE, with recent studies converging towards a range of values between 400 and 500 mm SLE global glacier ice mass. But as mentioned above, a part of this ice mass will not contribute to sea level.

Results (trends)

Table 3 presents most recent estimates of trends in global glacier mass balances.

The ensemble mean contribution of glaciers to sea-level rise for the time period 1993-2016 is 0.65 ± 0.051 mm yr -1 SLE and 0.74 ± 0.18 mm yr -1 for the time period 2005-2016 (uncertainties are averaged). Different studies refer to dif- The main source of uncertainty is that the vast majority of glaciers are unmeasured, which makes interpolation or extrapolation necessary, whether for in situ or satellite measurements, as well as for glacier modeling. Other main contributions to uncertainty in the ensemble mean stem from methodological differences, such as the downscaling of atmospheric forcing required for glacier modeling, the separation of glacier mass change to other mass change in the spatial gravimetry signal and the derivation of observational estimates of mass change from different raw measurements (e.g., length and volume changes, mass balance measurements, and geodetic methods), all with their specific uncertainties.

Greenland

Ice sheets are the largest potential source of future sea-level rise and represent the largest uncertainty in projections of future sea level. Almost all land ice (∼ 99.5 %) is locked in the ice sheets, with a volume in sea-level equivalent (SLE) terms of 7.4 m for Greenland and 58.3 m for Antarctica. It has been estimated that approximately 25 % to 30 % of the total land ice contribution to sea-level rise over the last decade came from the Greenland ice sheet (e.g., Dieng et al., 2017;[START_REF] Box | Sea level rise contribution from Arctic land ice: 1850-2100, Snow, Water, Ice and Permafrost in the Arctic (SWIPA) 2017[END_REF].

There are three main methods that can be used to estimate the mass balance of the Greenland ice sheet: (1) measurement of changes in elevation of the ice surface over time (dh / dt) either from imagery or altimetry; (2) the mass budget or input-output method (IOM), which involves estimating the difference between the surface mass balance and ice discharge; and (3) consideration of the redistribution of mass via gravity anomaly measurements, which only became viable with the launch of GRACE in 2002. Uncertainties due to the GIA correction are small in Greenland compared Barletta et al. (2013) 2003-2016 GRACE [START_REF] Groh | The method of tailored sensitivity kernels for GRACE mass change estimates[END_REF] 2003-2015 GRACE Update from [START_REF] Luthcke | Antarctica, Greenland and Gulf of Alaska landice evolution from an iterated GRACE global mascon solution[END_REF] 2003-2015 GRACE Update from [START_REF] Sasgen | Timing and origin of recent regional ice-mass loss in Greenland[END_REF] 2003-2016 GRACE Update from [START_REF] Schrama | A mascon approach to assess ice sheet and glacier mass balances and their uncertainties from GRACE data[END_REF] 2003-2016 GRACE Update from van den Broeke et al. (2016) 1993-2016 Input-output method (IOM) Wiese et al. (2016a, b) 2003-2016 GRACE Update from [START_REF] Wouters | GRACE observes small-scale mass loss in Greenland[END_REF] 2003-2016 GRACE to Antarctica: on the order of ±20 Gt yr -1 mass equivalent (Khan et al., 2016). Prior to 2003, mass trends are reliant on IOM and altimetry. Both techniques have limited sampling in time and/or space for parts of the satellite era (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002) and errors for this earlier period are, therefore, higher (van den Broeke et al., 2016;Hurkmans et al., 2014). The consistency between the three methods mentioned above was demonstrated for Greenland by [START_REF] Sasgen | Timing and origin of recent regional ice-mass loss in Greenland[END_REF] for the period 2003-2009. Ice-sheet-wide estimates showed excellent agreement although there was less consistency at a basin scale. We have, therefore, high confidence and relatively low uncertainties in the mass rates for the Greenland ice sheet in the satellite era (see also [START_REF] Bamber | The land ice contribution to sea level during the satellite era[END_REF].

Datasets considered for the assessment

This assessment of sea-level budget contribution from the Greenland ice sheet considers the datasets shown in Table 4.

Methods and analyses

All but one of these datasets are based on GRACE data and therefore provide annual time series from ∼ 2002 onwards. The one exception uses IOM (van Den Broeke et al., 2016) to give an annual mass time series for a longer time period (1993 onwards).

Notwithstanding this, each group has chosen their own approach to estimate mass balance from GRACE observations. As the aim of this global sea-level budget assessment is to compile existing results (rather than undertake new analyses), we have not imposed a specific methodology. Instead, we asked for the contributed datasets to reflect each group's 'best estimate' of annual trends for Greenland using the method(s) they have published.

Greenland contains glaciers and ice caps (GIC) around the margins of the main ice sheet, often referred to as peripheral GIC (PGIC), which are a significant proportion of the total mass imbalance (circa 15-20 %) [START_REF] Bolch | Mass loss of Greenland's glaciers and ice caps 2003-2008 revealed from ICESat laser altimetry data[END_REF]. Some studies consider the mass balance of the ice sheets and the PGIC separately but there has been, in general, no consistency in the treatment of PGIC and many studies do not

Results

There is generally a good level of agreement between the datasets (Fig. 9), and taken together they provide an average estimate of 171 Gt yr -1 of ice mass loss (or sea-level budget contribution) from Greenland for the period 1993 to 2016, increasing to 272 Gt yr -1 for the period 2005 to 2016 (Table 5).

All the datasets illustrate the previously documented accelerating mass loss up to 2012 (Rignot et al., 2011a;[START_REF] Velicogna | Increasing rates of ice mass loss from the Greenland and Antarctic ice sheets revealed by GRACE[END_REF] . In 2012, the ice sheet experienced exceptional surface melting reaching as far as the summit (Nghiem et al., 2012) and a record mass loss, since at least 1958, of over 400 Gt (van Den Broeke et al., 2016). The following years, however, show a reduced loss (not more than 270 Gt in any year). Inclusion of the years since 2012 in the 2005-2016 trend estimate reduces the overall rate of mass loss acceleration and its statistical significance. There is greater divergence in the GRACE time series for 2016. We associate this with the degradation of the satellites as they came to-Earth Syst. Sci. Data, 10, 1551Data, 10, -1590Data, 10, , 2018 www.earth-syst-sci-data.net/10/1551/2018/ wards the end of their mission. For 2005-2012, it might be inferred that there is a secular trend towards greater mass loss and from 2010 to 2012 the value is relatively constant. Interannual variability in mass balance of the ice sheet is driven, primarily, by the surface mass balance (i.e., atmospheric weather) and it is apparent that the magnitude of this year-to-year variability can be large: it exceeded 360 Gt (or 1 mm sea-level equivalent) between 2012 and 2013. Caution is required, therefore, in extrapolating trends from a short record such as this.

Antarctica

The annual turnover of mass of Antarctica is about 2200 Gt yr -1 (over 6 mm yr -1 of SLE), 5 times larger than in Greenland (Wessem et al., 2017). In contrast to Greenland, ice and snow melt have a negligible influence on Antarctica's mass balance, which is therefore completely controlled by the balance between snowfall accumulation in the drainage basins and ice discharge along the periphery. The continent is also 7 times larger than Greenland, which makes satellite techniques absolutely essential to survey the continent. Interannual variations in accumulation are large in Antarctica, showing decadal to multidecadal variability, so that many years of data are required to extract trends, and missions limited to only a few years may produce misleading results (e.g., Rignot et al., 2011a, b).

As in Greenland, the estimation of the mass balance has employed a variety of techniques, including (1) the gravity method with GRACE since April 2002 until the end of the mission in late 2016; (2) the IOM method using a series of Landsat and synthetic-aperture radar (SAR) satellites for measuring ice motion along the periphery (Rignot et al., 2011a, b), ice thickness from airborne depth radar sounders such as Operation IceBridge (Leuschen, 2014a), and reconstructions of surface mass balance using regional atmospheric climate models constrained by re-analysis data (RACMO, MAR and others); and (3) a radar or laser altimetry method which mixes various satellite altimeters and correct ice elevation changes with density changes from firm models. The largest uncertainty in the GRACE estimate in Antarctica is the GIA, which is larger than in Greenland, and a large fraction of the observed signal. The IOM method compares two large numbers with large uncertainties to estimate the mass balance as the difference. In order to detect an imbalance at the 10 % level, surface mass balance and ice discharge need to be estimated with a precision typically of 5 to 7 %. The altimetry method is limited to areas of shallow slope; hence, it is difficult to use in the Antarctic Peninsula and in the deep interior of the Antarctic continent due to unknown variations of the penetration depth of the signal in snow and firn. The only method that expresses the partitioning of the mass balance between surface processes and dynamic processes is the IOM method (e.g., Rignot et al., 2011a). The gravity method is an integrand method which does not suffer from the limitations of surface mass balance models but is limited in spatial resolution (e.g., [START_REF] Velicogna | Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time-variable gravity data[END_REF]. The altimetry method provides independent evidence of changes in ice dynamics, e.g., by revealing rapid ice thinning along the ice streams and glaciers revealed by ice motion maps, as opposed to large-scale variations reflecting a variability in surface mass balance [START_REF] Mcmillan | Increased ice losses from Antarctica detected by Cryosat-2[END_REF].

All these techniques have improved in quality over time and have accumulated a decade to several decades of observations, so that we are now able to assess the mass balance of the Antarctic continent using methods with reasonably low uncertainties and multiple lines of evidence as the methods are largely independent, which increases confidence in the results (see recent publication by the IMBIE Team, 2018). There is broad agreement in the mass loss from the Antarctic Peninsula and West Antarctica; most residual uncertainties are associated with East Antarctica as the signal is relatively small compared to the uncertainties, although most estimates tend to indicate a low contribution to sea level (e.g., [START_REF] Shepherd | A reconciled estimate of ice-sheet mass balance[END_REF].

Datasets considered for the assessment

This assessment considers the datasets shown in Table 6.

In Table 6, the negative trend estimate by Zwally et al. ( 2016) is not added. It is worth noting that including it would only slightly reduce the ensemble mean trend.

Methods and analyses

The datasets used in this assessment are Antarctica mass balance time series generated using different approaches. Two estimates are a joint inversion of GRACE, altimetry and GPS data [START_REF] Martín-Español | Spatial and temporal Antarctic Ice Sheet mass trends, glacio-isostatic adjustment, and surface processes from a joint inversion of satellite altimeter, gravity, and GPS data[END_REF] as well as GRACE and CryoSat data (Forsberg et al., 2017). Two methods are mascon solutions obtained from the GRACE intersatellite rangerate measurements over equal-area spherical caps covering the Earth's surface [START_REF] Luthcke | Antarctica, Greenland and Gulf of Alaska landice evolution from an iterated GRACE global mascon solution[END_REF][START_REF] Wiese | JPL GRACE Mascon Ocean, Ice, and Hydrology Equivalent Water Height RL05M. 1 CRI Filtered[END_REF], three estimates use the GRACE spherical harmonics solutions [START_REF] Velicogna | Regional acceleration in ice mass loss from Greenland and Antarctica using GRACE time-variable gravity data[END_REF][START_REF] Wiese | JPL GRACE Mascon Ocean, Ice, and Hydrology Equivalent Water Height RL05M. 1 CRI Filtered[END_REF][START_REF] Wouters | Limits in detecting acceleration of ice sheet mass loss due to climate variability[END_REF] and one uses gridded GRACE products [START_REF] Wouters | Limits in detecting acceleration of ice sheet mass loss due to climate variability[END_REF].

All GRACE time series were provided as monthly time series (except for the one using the Martín-Español et al., 2016, method, which was provided as annual estimates). In addition, different groups use different GIA corrections, therefore the spread of the trend solutions also represents the error associated with the GIA correction which, in Antarctica, is the largest source of uncertainty. [START_REF] Wouters | Limits in detecting acceleration of ice sheet mass loss due to climate variability[END_REF] used their own GIA solution [START_REF] Sasgen | Joint inversion estimate of regional glacial isostatic adjustment in Antarctica considering a lateral varying Earth structure (ESA STSE Project REGINA)[END_REF] In addition, [START_REF] Groh | The method of tailored sensitivity kernels for GRACE mass change estimates[END_REF] did not include the peripheral glaciers and ice caps, while all other estimates do.

Table 6 shows the Antarctic contribution to sea level during 2005-2015 from the different GRACE solutions, and for the input and output method. There is a single IOM-based dataset that provides trends for the period 1993-2015 (update of Rignot et al., 2011a). For the period 2005-2015, we calculated the annual sea-level contribution from Antarctica using GRACE and IOM estimates (Table 7).

As we are interested in evaluating the long-term trend and interannual variability of the Antarctic contribution to sea level, for each GRACE dataset available in monthly time se- ries, we first removed the annual and subannual components of the signal by applying a 13-month averaging filter and we then used the smoothed time series to calculate annual mass change. Figure 10 shows the annual sea-level contribution from Antarctica calculated from the GRACE-derived estimates and for the input-output method. The GRACE mean annual estimates are calculated as the mean of the annual contributions from the different groups, and the associated error calculated as the sum of the spread of the annual estimates and the mean annual error.

Results

There is generally broad agreement between the GRACE datasets (Fig. 10), as most of the differences between GRACE estimates are caused by differences in the GIA correction. We find a reasonable agreement between GRACE and the IOM estimates although the IOM estimates indicate higher losses. Taken together, these estimates yield an average of 0.42 mm yr -1 sea-level budget contribution from Antarctica for the period 2005 to 2015 (Table 7) and 0.25 mm yr -1 sea level for the time period 1993-2005, where the latter value is based on IOM only.

All the datasets illustrate the previously documented accelerating mass loss of Antarctica (Rignot et al., 2011a, b;[START_REF] Velicogna | Increasing rates of ice mass loss from the Greenland and Antarctic ice sheets revealed by GRACE[END_REF]. In 2005-2010, the ice sheet experienced ice mass loss driven by an increase in mass loss in the Amundsen Sea sector of West Antarctica [START_REF] Mouginot | Sustained increase in ice discharge from the Amundsen Sea Embayment, West Antarctica, from 1973 to 2013[END_REF]. The following years showed a reduced increase in mass loss, as colder ocean conditions prevailed in the Amundsen Sea embayment sector of West Antarctica in 2012-2013 which reduced the melting of the ice shelves in Earth Syst. Sci. Data, 10, 1551Data, 10, -1590Data, 10, , 2018 www.earth-syst-sci-data.net/10/1551/2018/ front of the glaciers [START_REF] Dutrieux | Strong sensitivity of Pine Island ice shelf melting to climatic variability[END_REF]. Divergence in the GRACE time series is observed after 2015 due to the degradation of the satellites towards the end of the mission.

The large interannual variability in mass balance in 2005-2015, characteristic of Antarctica, nearly masks out the trend in mass loss, which is more apparent in the longer time series than in short time series. The longer record highlights the pronounced decadal variability in ice sheet mass balance in Antarctica, demonstrating the need for multidecadal time series in Antarctica, which have been obtained only by IOM and altimetry. The interannual variability in mass balance is driven almost entirely by surface mass balance processes. The mass loss of Antarctica, about 200 Gt yr -1 in recent years, is only about 10 % of its annual turnover of mass (2200 Gt yr -1 ), in contrast with Greenland where the mass loss has been growing rapidly to nearly 100 % of the annual turnover of mass. This comparison illustrates the challenge of detecting mass balance changes in Antarctica, but at the same time, that satellite techniques and their interpretation have made tremendous progress over the last 10 years, producing realistic and consistent estimates of the mass using a number of independent methods [START_REF] Bamber | The land ice contribution to sea level during the satellite era[END_REF]the IMBIE Team, 2018).

Terrestrial water storage

Human transformations of the Earth's surface have impacted the terrestrial water balance, including continental patterns of river flow and water exchange between land, atmosphere and ocean, ultimately affecting global sea level. For instance, massive impoundment of water in man-made reservoirs has reduced the direct outflow of water to the sea through rivers, while groundwater abstractions, wetland and lake storage losses, deforestation, and other land use changes have caused changes to the terrestrial water balance, including changing evapotranspiration over land, leading to net changes in landocean exchanges [START_REF] Chao | Impact of artificial reservoir water impoundment on global sea level[END_REF]Wada et al., 2012a, b;Konikow, 2011;[START_REF] Church | Sea level change, in: Climate Change 2013: The Physical Science Basis[END_REF]Döll et al., 2014a, b). Overall, the combined effects of direct anthropogenic processes have reduced land water storage, increasing the rate of sea-level rise by 0.3-0.5 mm yr -1 during recent decades [START_REF] Church | Sea level change, in: Climate Change 2013: The Physical Science Basis[END_REF][START_REF] Gregory | Twentieth-Century Global-Mean Sea Level Rise: Is the Whole Greater than the Sum of the Parts?[END_REF]Wada et al., 2016). Additionally, recent work has shown that climatedriven changes in water stores can perturb the rate of sealevel change over interannual to decadal timescales, making global land mass budget closure sensitive to varying observational periods (Cazenave et al., 2014;Dieng et al., 2015a;[START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF][START_REF] Rietbroek | Revisiting the contemporary sea-level budget on global and regional scales[END_REF]). Here we discuss each of the major component contributions from land, with a summary in Table 8, and estimate the net terrestrial water storage contribution to sea level.

Direct anthropogenic changes in terrestrial water storage

Water impoundment behind dams

Global groundwater depletion

Groundwater currently represents the largest secular trend component to the land water storage budget. The rate of groundwater depletion (GWD) and its contribution to sea level has been subject to debate [START_REF] Gregory | Twentieth-Century Global-Mean Sea Level Rise: Is the Whole Greater than the Sum of the Parts?[END_REF]Taylor et al., 2013). In the IPCC AR4 (Solomon et al., 2007), the contribution of nonfrozen terrestrial waters (including GWD) to sea-level variation was not considered due to its perceived uncertainty (Wada et al., 2016). Observations from GRACE opened a path to monitor total water storage changes, including groundwater in data-scarce regions [START_REF] Strassberg | Comparison of seasonal terrestrial water storage variations from GRACE with groundwater-level measurements from the High Plains Aquifer (USA)[END_REF][START_REF] Rodell | Satellite-based estimates of groundwater depletion in India[END_REF][START_REF] Tiwari | Dwindling groundwater resources in northern India, from satellite gravity observations[END_REF][START_REF] Shepherd | A reconciled estimate of ice-sheet mass balance[END_REF][START_REF] Shamsudduha | Monitoring groundwater storage changes in the highly seasonal humid tropics: Validation of GRACE measurements in the Bengal Basin[END_REF]Voss et al., 2013). Some studies have also applied global hydrological models in combination with the GRACE data (see Wada et al., 2016, for a review).

Earlier estimates of GWD contribution to sea level range from 0.075 to 0.30 mm yr -1 [START_REF] Sahagian | Direct anthropogenic contributions to sea level rise in the twentieth century[END_REF][START_REF] Gornitz | Sea-level rise: A review of recent past and near-future trends[END_REF][START_REF] Gornitz | Sea Level Rise: History and Consequences[END_REF]Foster and Loucks, 2006). More recently, Wada et al. (2012b), using hydrological modeling, estimated that the contribution of GWD to global sea level increased from 0.035 (±0.009) to 0.57 (±0.09) mm yr -1 during the 20th century and projected that it would further increase to 0.82 (±0.13) mm yr -1 by 2050. Döll et al. (2014b) used hydrological modeling, well observations and GRACE satellite gravity anomalies to estimate a 2000-2009 global GWD of 113 km 3 yr -1 (0.314 mm yr -1 SLE). This value represents the impact of human groundwater withdrawals only and does not consider the effect of climate variability on groundwater storage. A study by Konikow (2011) estimated global GWD to be 145 (±39) km 3 yr -1 (0.41 ± 0.1 mm yr -1 SLE) during 1991-2008 based on measurements of changes in groundwater storage from in situ observations, calibrated groundwater modeling, GRACE satellite data and extrapolation to unobserved aquifers.

An assumption of most existing global estimates of GWD impacts on sea-level change is that nearly 100 % of the GWD ends up in the ocean. However, groundwater pumping can also perturb regional climate due to land-atmosphere interactions (Lo and Famiglietti, 2013). A recent study by Wada et al. (2016) used a coupled land-atmosphere model simulation to track the fate of water pumped from underground and found it more likely that ∼ 80 % of the GWD ends up in the ocean over the long term, while 20 % re-infiltrates and remains in land storage. They estimated an updated contribution of GWD to global sea-level rise ranging from 0.02 (±0.004) mm yr -1 in 1900 to 0.27 (±0.04) mm yr -1 in 2000 (Fig. 11). This indicates that previous studies had likely overestimated the cumulative contribution of GWD to global SLR during the 20th century and early 21st century by 5-10 mm.

Land cover and land-use change

Humans have altered a large part of the land surface, replacing about 40 % of natural vegetation by anthropogenic land cover such as crop fields or pasture. Such land cover change can affect terrestrial hydrology by changing the infiltrationto-runoff ratio and can impact subsurface water dynamics by modifying recharge and increasing groundwater storage [START_REF] Strassberg | Comparison of seasonal terrestrial water storage variations from GRACE with groundwater-level measurements from the High Plains Aquifer (USA)[END_REF]. The combined effects of anthropogenic land cover changes on land water storage can be quite complex. Using a combined hydrological and water resource model, Bosmans et al. (2017) estimated that land cover change between 1850 and 2000 has contributed to a discharge increase of 1058 km 3 yr -1 , on the same order of magnitude as the effect of human water use. These recent model results suggest that land-use change is an important topic for further investigation in the future. So far, this contribution remains highly uncertain.

Deforestation and afforestation

At present, large losses in tropical forests and moderate gains in temperate-boreal forests result in a net reduction of global forest cover (FAO, 2015;[START_REF] Keenan | Dynamics of global forest area: Results from the FAO Global Forest Resources Assessment 2015[END_REF][START_REF] Macdicken | Global Forest Resources Assessment, What, why and how?[END_REF][START_REF] Sloan | Forest Resources Assessment of 2015 shows positive global trends but forest loss and degradation persist in poor tropical countries[END_REF]. Net deforestation releases carbon and water stored in both biotic tissues and Earth Syst. Sci. Data, 10, 1551Data, 10, -1590Data, 10, , 2018 www.earth-syst-sci-data.net/10/1551/2018/ soil, which leads to sea-level rise through three primary processes: deforestation-induced runoff increases [START_REF] Gornitz | Effects of anthropogenic intervention in the land hydrologic cycle on global sea level rise, Global Planet[END_REF], carbon loss-related decay and plant storage loss, and complex climate feedbacks [START_REF] Butt | Evidence that deforestation affects the onset of the rainy season in Rondonia, Brazil[END_REF][START_REF] Chagnon | Contemporary climate change in the Amazon[END_REF][START_REF] Nobre | Amazon Deforestation and Climate Change in a Coupled Model Simulation[END_REF][START_REF] Shukla | Amazon Deforestation and Climate Change[END_REF]Spracklen et al., 2012). Due to these three causes, and if uncertainties from the land-atmospheric coupling are excluded, a summary by Wada et al. (2016) suggests that the current net global deforestation leads to an upper-bound contribution of ∼ 0.035 mm yr -1 SLE.

Wetland degradation

Wetland degradation contributes to sea level primarily through (i) direct water drainage or removal from standing inundation, soil moisture and plant storage, and (ii) water release from vegetation decay and peat combustion. Wada et al. (2016) consider a recent wetland loss rate of 0.565 % yr -1 since 1990 (Davidson, 2014) and a present global wetland area of 371 mha averaged from three databases: Matthews natural wetlands [START_REF] Matthews | Methane emission from natural wetlands: Global distribution, area, and environmental characteristics of sources[END_REF], ISLSCP (Darras, 1999), and DISCover [START_REF] Belward | The IGBP-DIS global 1-km land-cover data set DISCover: A project overview, Photogramm[END_REF][START_REF] Lovel | The IGBP-DIS global 1 km land cover data set, DISCover: first results[END_REF]. They assume a uniform 1 m depth of water in wetlands [START_REF] Milly | Terrestrial water-storage contributions to sealevel rise and variability[END_REF], to estimate a contribution of recent global wetland drainage to sea level of 0.067 mm yr -1 . Wada et al. ( 2016) apply a wetland area and loss rate as used for assessing wetland water drainage, to determine that the annual reduction of wetland carbon stock since 1990, if completely emitted, releases water equivalent to 0.003-0.007 mm yr -1 SLE. Integrating the impacts of wetland drainage, oxidation and peat combustion, Wada et al. (2016) suggest that the recent global wetland degradation results in an upper bound of 0.074 mm yr -1 SLE.

Lake storage changes

Lakes store the greatest mass of liquid water on the terrestrial surface [START_REF] Oki | Global hydrological cycles and world water resources[END_REF]), yet, because of their "dynamic" nature (Sheng et al., 2016;[START_REF] Wang | Drained thaw lake basin recovery on the western Arctic Coastal Plain of Alaska using high-resolution digital elevation models and remote sensing imagery[END_REF], their overall contribution to sea level remains uncertain. In the past century, perhaps the greatest contributor in global lake storage was the Caspian Sea [START_REF] Milly | Terrestrial water-storage contributions to sealevel rise and variability[END_REF], where the water level exhibits substantial oscillations attributed to meteorological, geological, and anthropogenic factors [START_REF] Ozyavas | A possible connection of Caspian Sea level fluctuations with meteorological factors and seismicity[END_REF]Chen et al., 2017a). Assuming the lake level variation kept pace with groundwater changes [START_REF] Sahagian | Direct anthropogenic contributions to sea level rise in the twentieth century[END_REF], the overall contribution of the Caspian Sea, including both surface and groundwater storage variations through 2014, has been about 0.03 mm yr -1 SLE since 1900, 0.075 (±0.002) mm yr -1 since 1995 and 0.109 (±0.004) mm yr -1 since 2002. Additionally, between 1960 and 1990, the water storage in the Aral Sea Basin declined at a striking rate of 64 km 3 yr -1 , equivalent to 0.18 mm yr -1 SLE [START_REF] Sahagian | Global physical effects of anthropogenic hydrological alterations: sea level and water redistribution[END_REF][START_REF] Sahagian | Direct anthropogenic contributions to sea level rise in the twentieth century[END_REF][START_REF] Vörösmarty | Anthropogenic disturbance of the terrestrial water cycle[END_REF] due mostly to upstream water diversion for irrigation (Perera, (Schwatke et al., 2015). Assuming that groundwater drainage has kept pace with lake level reduction [START_REF] Sahagian | Direct anthropogenic contributions to sea level rise in the twentieth century[END_REF], the Aral Sea has contributed 0.0358 (±0.0003) mm yr -1 to the recent sealevel rise.

Water cycle variability

Natural changes in the interannual to decadal cycling of water can have a large effect on the apparent rate of sea-level change over decadal and shorter time periods [START_REF] Milly | Contribution of climate-driven change in continental water storage to recent sealevel rise[END_REF][START_REF] Lettenmaier | Land waters and sea level[END_REF][START_REF] Llovel | Global land water storage change from GRACE over 2002-2009; Inference on sea level[END_REF]. For instance, ENSO-driven modulations of the global water cycle can be important in decadal-scale sea-level budgets and can mask underlying secular trends in sea level (Fasullo et al., 2013;Cazenave et al., 2014;[START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF]. Sea-level variability due to climate-driven hydrology represents a super-imposed variability on the secular rates of global mean sea-level rise. While this term can be large and is important in the interpretation of the sea-level record, it is arguably the most difficult term in the land water budget to quantify.

Net terrestrial water storage

GRACE-based estimates

Measurements of non-ice-sheet continental land mass from GRACE satellite gravity have been presented in several recent studies (Jensen et al., 2013[START_REF] Rietbroek | Revisiting the contemporary sea-level budget on global and regional scales[END_REF][START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF][START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] and can be used to 2016), GRACE shows a total gain in land water storage during the 2002-2014 period, corresponding to a sea-level trend of -0.33±0.16 mm yr -1 SLE (modified from [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF]. These trends include all human-driven and climate-driven processes in Table 1 and can be used to close the land water budget over the study period. constrain a global land mass budget. Note that these "topdown" estimates contain both climate-driven and direct anthropogenic driven effects, which makes them most useful in assessing the total impact of land water storage changes and closing the budget of all contributing terms. GRACE observations, when averaged over the whole land domain following [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF], indicate a total TWS change (including glaciers) over the 2002-2014 study period of approximately +0.32 ± 0.13 mm yr -1 SLE (i.e., ocean gaining mass). Global mountain glaciers have been estimated to lose mass at a rate of 0.65 ± 0.09 mm yr -1 (e.g., [START_REF] Gardner | A reconciled estimate of glacier contributions to sea level rise: 2003 to 2009[END_REF][START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF] during that period, such that a mass balance indicates that global glacier-free land gained water at a rate of -0.33 ± 0.16 mm yr -1 SLE (i.e., ocean losing mass; Fig. 12). A roughly similar estimate was found from GRACE using glacier-free river basins globally (-0.21 ± 0.09 mm yr -1 ) [START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF]. Thus, the GRACE-based net TWS estimates suggest a negative sealevel contribution from land over the GRACE period (Table 8). However, mass change estimates from GRACE incorporate uncertainty from all potential error sources that arise in processing and postprocessing of the data, including from the GIA model, and from the geocenter and mean pole corrections.

Estimates based on global hydrological models

Global land water storage can also be estimated from global hydrological models (GHMs) and global land surface models. These compute water, or water and energy balances, at the Earth's surface, yielding time variations of water storage in response to prescribed atmospheric data (temperature, humidity and wind) and the incident water and energy fluxes from the atmosphere (precipitation and radiation). Meteorological forcing is usually based on atmospheric model reanalysis. Model uncertainties result from several factors. Recent work has underlined the large differences among different state-of-the-art precipitation datasets (Beck et al., 2017), with large impacts on model results at seasonal (Schellekens et al., 2017) and longer timescales [START_REF] Felfelani | Natural and human-induced terrestrial water storage change: A global analysis using hydrological models and GRACE[END_REF]. Another source of uncertainty is the treatment of subsurface storage in soils and aquifers, as well as dynamic changes in storage capacity due to representation of frozen soils and permafrost, the complex effects of dynamic vegetation, atmospheric vapor pressure deficit estimation and an insufficiently deep soil column. A recent study by [START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] compared water storage trends from five global land surface models and two global hydrological models to GRACE storage trends and found that models estimated the opposite trend in net land water storage to GRACE over the 2002-2014 period. These authors attributed this discrepancy to model deficiencies, in particular soil depth limitations. These combined error sources are responsible for a range of storage trends across models of approximately 0.5±0.2 mm yr -1 SLE. In terms of global land average, model differences can cause up to ∼ 0.4 mm yr -1 SLE uncertainty.

Synthesis

Based on the different approaches to estimate the net land water storage contribution, we estimate that the corresponding sea-level rate ranges from -0.33 to 0.23 mm yr -1 during the period of 2002-2014/15 due to water storage changes (Table 8). According to GRACE, the net TWS change (i.e., not including glaciers) over the period 2002-2014 shows a negative contribution to sea level of -0.33 and -0.21 mm yr -1 by [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF] and [START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] respectively. Such a negative signal is not currently reproduced by hydrological models which estimate slightly positive trends over the same period (see Table 8). It is to be noted, however, that looking at trends only over periods on the order of a decade may not be appropriate due the strong interannual variability of TWS at basin and global scales. For example, Fig. 5 from [START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] (see also Fig. S9 from their Supplement), which compares GRACE TWS and model estimates over large river basins over 2002-2014, clearly shows that the discrepancies between GRACE and models occur at the end of the record for the majority of basins. This is particularly striking for the Amazon basin (the largest contributor to TWS), for which GRACE and models agree reasonably well until 2011, and then depart significantly, with GRACE TWS showing a strongly positive trend since then, unlike the models. Such a divergence at the end of the record is also noticed for several other large basins (see Scanlon et al., 2018, Fig. S9). No clear explanation can be provided yet, even though one may question the quality of the meteorological forcing used by hydrological models for the recent years. But this calls for some caution when comparing GRACE and other models on the basis of trends only because of the dominant interannual variability of the TWS component. Much more work is needed to understand differences among models, and between models and GRACE. Of all components entering in the sea-level budget, the TWS contribution currently appears to be the most uncertain one.

Glacial isostatic adjustment

The Earth's dynamic response to the waxing and waning of the late-Pleistocene ice sheets is still causing isostatic disequilibrium in various regions of the world. The accompanying slow process of GIA is responsible for regional and global fluctuations in relative and absolute sea level, 3-D crustal deformations, and changes in the Earth's gravity field (for a review, see [START_REF] Spada | Glacial isostatic adjustment and contemporary sea level rise: An overview[END_REF]. To isolate the contribution of current climate change, geodetic observations must be corrected for the effects of GIA (King et al., 2010). These are obtained by solving the "sea-level equation" [START_REF] Farrell | On postglacial sea level[END_REF][START_REF] Mitrovica | On post-glacial sea level: I. General theory[END_REF]. The sea level can be expressed as S = N -U , where S is the rate of change in sea level relative to the solid Earth, N is the geocentric rate of sea-level change, and U is the vertical rate of displacement of the solid Earth. The sea-level equation accounts for solid Earth deformational, gravitational and rotational effects on sea level, which are sensitive to the Earth's mechanical properties and to the melting chronology of continental ice. Forward GIA modeling, based on the solution of the sealevel equation, provides predictions of unique spatial patterns (or fingerprints; see [START_REF] Plag | Inversion of global tide gauge data for present-day ice load changes[END_REF] of relative and geocentric sea-level change (e.g., [START_REF] Milne | Identifying the causes of sea-level change[END_REF]Kopp et al., 2015). During recent decades, the two fundamental components of GIA modeling have been progressively constrained from the observed history of relative sea level during the Holocene (see, e.g., [START_REF] Lambeck | Sea Level Change Through the Last Glacial Cycle[END_REF][START_REF] Peltier | Global glacial isostasy and the surface of the ice-age Earth: the ICE-5G (VM2) model and GRACE[END_REF]. In the context of climate change, the importance of GIA has been recognized since the mid-1980s, when the awareness of global sea-level rise stimulated the evaluation of the isostatic contribution to tide gauge observations (see Table 1 in [START_REF] Spada | New estimates of secular sea level rise from tide gauge data and GIA modelling[END_REF]. Subsequently, GIA models have been applied to the study of the pattern of sealevel change from satellite altimetry [START_REF] Tamisiea | Ongoing glacial isostatic contributions to observations of sea level change[END_REF], and since 2002 to the study of the gravity field variations from GRACE. Our primary goal here is to analyze GIA model outputs that have been used to infer global mean sea-level change and ice sheet volume change from geodetic datasets during the altimetry era. These outputs are the sea-level variations detected by satellite altimetry across oceanic regions (n), the ocean mass change (w) and the modern ice sheets mass balance from GRACE. We also discuss the GIA correction that needs to be applied to GRACE-based land water storage changes. The GIA correction applied to tide-gaugebased sea-level observations at the coastlines is not discussed here. Since GIA evolves on timescales of millennia (e.g., [START_REF] Turcotte | Geodynamics[END_REF], the rate of change of all the isostatic signals can be considered constant on the timescale of interest.

GIA correction to altimetry-based sea level

Unlike tide gauges, altimeters directly sample the sea surface in a geocentric reference frame. Nevertheless, GIA contributes significantly to the rates of absolute sea-level change observed over the "altimetry era", which require a correction N gia that is obtained by solving the SLE (e.g., [START_REF] Spada | Glacial isostatic adjustment and contemporary sea level rise: An overview[END_REF]. As discussed in detail by [START_REF] Tamisiea | Ongoing glacial isostatic contributions to observations of sea level change[END_REF], N gia is sensitive to the assumed rheological profile of the Earth and to the history of continental glacial ice sheets. The variance of N gia over the surface of the oceans is much reduced, being primarily determined by the change in the Earth's gravity potential, apart from a spatially uniform shift. As discussed by [START_REF] Spada | Spectral analysis of sea level during the altimetry era, and evidence for GIA and glacial melting fingerprints[END_REF], the GIA contribution N gia is strongly affected by variations in the centrifugal potential associated with Earth's rotation, whose fingerprint is www.earth-syst-sci-data.net/10/1551/2018/ Earth Syst. Sci. Data, 10, 1551-1590, 2018 GRACE-based estimates of total land water storage (not including glaciers) -0.20 to -0.33 [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF][START_REF] Rietbroek | Revisiting the contemporary sea-level budget on global and regional scales[END_REF][START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] (±0.09-0.16) dominated by a spherical harmonic contribution of degree l = 2 and order m = ±1. Since N gia has a smooth spatial pattern, the global GIA correction to altimetry data can be obtained by simply subtracting its average n =< N gia > over the ocean sampled by the altimetry missions. The computation of the GIA contribution N gia has been the subject of various investigations, based on different GIA models. The estimate by [START_REF] Peltier | Global glacial isostatic adjustment and modern instrumental records of relative sea level history[END_REF] of n = -0.30 mm yr -1 is based on the ICE-4G (VM2) GIA model. Such a value has been adopted in the majority of studies estimating the GMSL rise from altimetry. Since n appears to be small compared to the global mean sea-level rise from altimetry (∼ 3 mm yr -1 ), a more precise evaluation has not been of concern until recently. However, it is important to notice that n is of comparable magnitude as the GMSL trend uncertainty, currently estimated to be ∼ 0.3 mm yr -1 (see Sect. 2.2). In Table 9a, we summarize the values of n according to works in the literature where various GIA model models and averaging methods have been employed. Based on values in Table 9a for which a standard deviation is available, the average of n (weighted by the inverse of associated errors), assumed to represent the best estimate, is n = (-0.29 ± 0.02) mm yr -1 , where the uncertainty corresponds to 2σ .

GIA correction to GRACE-based ocean mass

GRACE observations of present-day gravity variations are sensitive to GIA, due to the sheer amount of rock material that is transported by GIA throughout the mantle and the resulting changes in surface topography, especially over the formerly glaciated areas. The continuous change in the gravity field results in a nearly linear signal in GRACE observations. Since the gravity field is determined by global mass redistribution, GIA models used to correct GRACE data need to be global as well, especially when the region of interest is represented by all ocean areas. To date, the only global ice reconstruction publicly available is provided by the University of Toronto. Their latest product, named ICE-6G, has been published and distributed in 2015 [START_REF] Peltier | Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model[END_REF]; note that the ice history has been simultaneously constrained with a specific Earth model, named VM5a. During the early period of the GRACE mission, the available Toronto model was ICE-5G (VM2) [START_REF] Peltier | Global glacial isostasy and the surface of the ice-age Earth: the ICE-5G (VM2) model and GRACE[END_REF] Over most of the oceans, the GIA signature is much smaller than over the continents. However, once integrated over the global ocean, the signal w due to GIA is about -1 mm yr -1 of equivalent sea-level change [START_REF] Chambers | Ocean mass from GRACE and glacial isostatic adjustment[END_REF], which is of the same order of magnitude as the total ocean mass Earth Syst. Sci. Data, 10, 1551-1590, 2018 www.earth-syst-sci-data.net/10/1551/2018/ change induced by increased ice melt [START_REF] Leuliette | Balancing the sea level budget[END_REF]. The main uncertainty in the GIA contribution to ocean mass change estimates, apart from the general uncertainty in ice history and Earth mechanical properties, originates from the importance of changes in the orientation of the Earth's rotation axis [START_REF] Chambers | Ocean mass from GRACE and glacial isostatic adjustment[END_REF][START_REF] Tamisiea | Ongoing glacial isostatic contributions to observations of sea level change[END_REF]. Different choices in implementing the so-called "rotational feedback" lead to significant changes in the resulting GIA contribution to GRACE estimates. The issue of properly accounting for rotational effects has not been settled yet [START_REF] Mitrovica | The rotational stability of an ice-age earth[END_REF][START_REF] Peltier | On the origins of Earth rotation anomalies: New insights on the basis of both "paleogeodetic" data and Gravity Recovery and Climate Experiment (GRACE) data[END_REF][START_REF] Mitrovica | Ice age Earth rotation[END_REF][START_REF] Martinec | The rotational feedback linearmomentum balance in glacial isostatic adjustment[END_REF]. Table 9b summarizes the values of the mass-rate GIA contribution w according to the literature, where various models and averaging methods are employed. The weighted average of the values in Table 9b, for which an assessment of the standard deviation is available, is w = (-1.44 ± 0.36) mm yr -1 (the uncertainty is 2σ ), which we assume represents the preferred estimate.

GIA correction to GRACE-based terrestrial water storage

As discussed in the previous section, the GIA correction to apply to GRACE over land is significant, especially in regions formerly covered by the ice sheets (Canada and Scandinavia). Over Canada, GIA models significantly differ. This is illustrated in Fig. 13, which shows difference between two models of GIA correction to GRACE over land, the A et al. ( 2013) and [START_REF] Peltier | On the origins of Earth rotation anomalies: New insights on the basis of both "paleogeodetic" data and Gravity Recovery and Climate Experiment (GRACE) data[END_REF] models. We see that over the majority of the land areas, differences are small, except over northern Canada, in particular around the Hudson Bay, where differences larger than ± 20 mm yr -1 SLE are noticed. This may affect GRACE-based TWS estimates over Canadian river basins.

When averaged over the whole land surface as done in some studies to estimate the combined effect of land water storage and glacier melting from GRACE (e.g., [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF] see Sect. 2.7), the GIA correction ranges from ∼ 0.5 to 0.7 mm yr -1 (in mm yr -1 SLE). Values for different GIA models are given in Table 9c.

GIA correction to GRACE-based ice sheet mass balance

The GRACE gravity field observations allow the determination of mass balances of ice sheets and large glacier systems with inaccuracy similar or superior to the input-output method or satellite laser and radar altimetry [START_REF] Shepherd | A reconciled estimate of ice-sheet mass balance[END_REF]. However, GRACE ice-mass balances rely on successfully separating and removing the apparent mass change related to GIA. While the GIA correction is small compared to the mass balance for the Greenland ice sheet (ca. < 10 %), its magnitude and uncertainty in Antarctica is on the order of the ice-mass balance itself (e.g., [START_REF] Martín-Español | Spatial and temporal Antarctic Ice Sheet mass trends, glacio-isostatic adjustment, and surface processes from a joint inversion of satellite altimeter, gravity, and GPS data[END_REF].

Particularly for today's glaciated areas, GIA remains poorly resolved due to the sparse data constraining the models, leading to large uncertainties in the climate history, the geometry and retreat chronology of the ice sheet, as well as the Earth structure. The consequences are ambiguous GIA predictions, despite fitting the same observational data. There are two principal approaches towards resolving GIA underneath the ice sheets. Empirical estimates can be derived that make use of the different sensitivities of satellite observations to ice-mass changes and GIA (e.g., [START_REF] Riva | Glacial isostatic adjustment over Antarctica from combined ICESat and GRACE satellite data[END_REF]Riva et al., , 2010;;[START_REF] Wu | Simultaneous estimation of global present-day water transport and glacial isostatic adjustment[END_REF]. Alternatively, GIA can be modeled numerically by forcing an Earth model with a fixed ice retreat scenario (e.g., [START_REF] Peltier | Closure of the budget of global sea level rise over the GRACE era: the importance and magnitudes of the required corrections for global glacial isostatic adjustment[END_REF][START_REF] Whitehouse | A new glacial isostatic adjustment model for Antarctica: calibrating the deglacial model using observations of relative sea-level and present-day uplift rates[END_REF] or with output from a thermodynamic ice sheet model [START_REF] Gomez | A 3-D coupled ice sheet-sea level model applied to Antarctica through the last 40 ky[END_REF][START_REF] Konrad | Potential of the solid-Earth response for limiting long-term West Antarctic Ice Sheet retreat in a warming climate[END_REF]. Values of GIA-induced apparent mass change for Greenland and Antarctica as listed in the literature should be applied with caution (Table 9d) when applying them to GRACE mass balances. Each of these estimates may rely on a different GRACE postprocessing strategy and may differ in the approach used for solving the gravimetric inverse problem (mascon analysis, forward-modeling, averaging kernels). Of particular concern is the modeling and filtering of the pole tide correction caused by the rotational variations related to GIA, affecting coefficients of harmonic degree l = 2 and order m = ±1. As mentioned above, agreement on the modeling of the rotational feedback has not been reached within the GIA community. Furthermore, the pole tide correction applied during the determination gravity-field solutions differs between the GRACE processing centers and may not be consistent with the GIA correction listed. This inconsistency may introduce a significant bias in the ice-mass balance estimates (e.g., Sasgen et al., 2013, Supplement). [START_REF] Wahr | The pole tide and its effect on GRACE time-variable gravity measurements: Implications for estimates of surface mass variations[END_REF] presented recommendations on how to treat the pole tides in GRACE analysis. However, a systematic intercomparison of the GIA predictions in terms of their low-degree coefficients and their consistency with the GRACE processing standards still need to be done. The GRACE-based ocean mass, Antarctica mass and terrestrial water storage changes are very model dependent. As these GIA corrections cannot be assessed from independent information, they represent a large source of uncertainties to the sea-level budget components based on GRACE.

Ocean mass change from GRACE

Since 2002, GRACE satellite gravimetry has provided a revolutionary means for measuring global mass change and redistribution at monthly intervals with unprecedented accuracy, and offered the opportunity to directly estimate ocean mass change due to water exchange between the ocean and other components of the Earth (e.g., ice sheets, mountain glaciers, terrestrial water). GRACE time-variable gravity data have been successfully applied in a series of studies of ice mass balance of polar ice sheets (e.g., [START_REF] Velicogna | Measurements of Time-Variable Gravity Show Mass Loss in Antarctica[END_REF][START_REF] Luthcke | Recent Greenland Ice Mass Loss by Drainage System from Satellite Gravity Observations[END_REF] and mountain glaciers (e.g., Tamisiea et al., 2005;[START_REF] Chen | Patagonia Icefield Melting Observed by GRACE[END_REF] and their contributions to global sea-level change. GRACE data can also be used to directly study long-term oceanic mass change or nonsteric sea-level change (e.g., [START_REF] Willis | Assessing the globally averaged sea level budget on seasonal to interannual time scales[END_REF][START_REF] Leuliette | Closing the sea level rise budget with altimetry, Argo, and GRACE[END_REF][START_REF] Cazenave | Sea level budget over 2003-2008: A reevaluation from GRACE space gravimetry, satellite altimetry and Argo[END_REF], and provide a unique opportunity to study interannual or long-term TWS change and its potential impacts on sea-level change (Richey et al., 2015;[START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF].

GRACE time-variable gravity data can be used to quantify ocean mass change from three different main approaches. One is through measuring ice mass balance of polar ice sheets and mountain glaciers and variations of TWS, and their contributions to the GMSL (e.g., [START_REF] Velicogna | Measurements of Time-Variable Gravity Show Mass Loss in Antarctica[END_REF][START_REF] Schrama | A mascon approach to assess ice sheet and glacier mass balances and their uncertainties from GRACE data[END_REF]. The second approach is to directly quantify ocean mass change using ocean basin mask (kernel) (e.g., [START_REF] Luthcke | Calculating trends from GRACE in the presence of large changes in continental ice storage and ocean mass[END_REF][START_REF] Llovel | Global land water storage change from GRACE over 2002-2009; Inference on sea level[END_REF][START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]. In the ocean basin kernel approach, coastal ocean areas within certain distance (e.g., 300 or 500 km) from the coast are excluded, in order to minimize contaminations from mass change signal over the land (e.g., glacial mass loss and TWS change). The third approach solves mass changes on land and over ocean at the same time via forward modeling (e.g., Chen et al., 2013;Yi et al., 2015). The forward modeling is a global inversion to reconstruct the "true" mass change magnitudes over land and ocean with geographical constraint of locations of the mass change signals, and can help effectively reduce leakage between land and ocean (Chen et al., 2013).

Estimates of ocean mass changes from GRACE are subject to a number of major error sources. These include (1) leakage errors from the larger signals over ice sheets and land hydrology due to GRACE's low spatial resolution (of at least a few to several hundred kilometers) and the need for coastal masking, (2) spatial filtering of GRACE data to reduce spatial noise, (3) errors and biases in geophysical model correc-tions (e.g., GIA, atmospheric mass) that need to be removed from GRACE observations to isolate oceanic mass change and/or polar ice sheets and mountain glaciers mass balance, and (4) residual measurement errors in GRACE gravity measurements, especially those associated with GRACE lowdegree gravity changes. In addition, how to deal with the absent degree-1 terms, i.e., geocenter motion in GRACE gravity fields, is expected to affect estimates of GRACE-based oceanic mass rates and ice mass balances.

With a different treatment of the GRACE land-ocean signal leakage effect through global forward modeling, Chen et al. ( 2013) estimated ocean mass rates using GRACE RL05 time-variable gravity solutions over the period 2005-2011. They demonstrated that the ocean mass change contributes up to 1.80 ± 0.47 mm yr -1 (over the same period), which is significantly larger than previous estimates over about the same period. Yi et al. (2015) further confirmed that correct calibration of GRACE data and appropriate treatment of GRACE leakage bias are critical to improve the accuracy of GRACE-estimated ocean mass rates. Table 10 summarizes different estimates of GRACE ocean mass rates. The uncertainty estimates of the listed studies (Table 10) are computed from different methods, with different considerations of error sources into the error budget, and represent different confidence levels.

As demonstrated in Chen et al. (2013), different treatments of just the degree-2 spherical harmonics of the GRACE gravity solution alone can lead to substantial differences in GRACE-estimated ocean mass rates (ranging from 1.71 to 2.17 mm yr -1 ). Similar estimates from GRACE gravity solutions from different data processing centers can also be different. In the meantime, long-term degree-1 spherical harmonics variation, representing long-term geocenter motion and neglected in some of the previous studies (due to the lack of accurate observations) are also expected to have a non-negligible effect on GRACE-derived ocean mass rates (Chen et al., 2013). Different methods for computing ocean mass change using GRACE data may also lead to different estimates (Chen et al., 2013;[START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]Jensen et al., 2013).

To help better understand the potential and uncertainty of GRACE satellite gravimetry in quantification of the ocean mass rate, Table 11 provides a comparison of GRACEestimated ocean mass rates over the period January 2005 to December 2016 based on different GRACE data products and different data processing methods, including the CSR, GFZ and JPL GRACE RL05 spherical harmonic solutions (i.e., the so-called GSM solutions), as well as CSR, JPL and GSFC mascon solutions (the available GSFC mascons only cover the period up to July 2016). The three GRACE GSM results (CSR, GFZ and JPL) are updates from Johnson and [START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF], with degree-2 zonal term replaced by satellite laser ranging results [START_REF] Cheng | Monthly estimates of C20 from 5 SLR satellites based on GRACE RL05 models[END_REF], geocenter motion from [START_REF] Swenson | Estimating geocenter variations from a combination of GRACE and ocean model output[END_REF]), GIA model from A et al. (2013), an averaging kernel with a land mask Table 9. Estimated contributions of GIA to the rate of absolute sea-level change observed by altimetry (a), to the rate of mass change observed by GRACE over the global oceans (b), to the rate of mass change observed by GRACE over land (c), and to Greenland and Antarctic ice sheets (c), during the altimetry era. The GIA corrections are expressed in millimeters per year SLE except over Greenland and Antarctica where values are given in gigatons per year (ice mass equivalent). Most of the GIA contributions are expressed as a value ± 1 standard deviation; a few others are given in terms of a plausible range, and for some the uncertainties are not specified. Tamisiea (2011) (Figs. 3 and4) -0.5 to -1.9 -0.9 to -1.5

Ocean average of a range of estimates varying the Earth model, and based on a restricted set, respectively. that extends out 300 km, and no destriping or smoothing, as described in [START_REF] Rhein | Observations: Ocean, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[END_REF]. An update of GRACE ocean mass rate from Chen et al. (2013) is also included for comparisons, which is based on the CSR GSM solutions using forward modeling (a global inversion approach), with similar treatments of the degree-2 zonal term, geocenter motion and GIA effects.

The JPL mascon ocean mass rate is computed from all mascon grids over the ocean, and the GSFC mascon ocean mass rate is computed from all ocean mascons, with the Mediterranean, Black and Red seas excluded. A coastline resolution improvement (CRI) filter is already applied in the JPL mascons to reduce leakage [START_REF] Wiese | JPL GRACE Mascon Ocean, Ice, and Hydrology Equivalent Water Height RL05M. 1 CRI Filtered[END_REF], and in both the GSFC and JPL mascon solutions, the ocean and land are separately defined [START_REF] Luthcke | Antarctica, Greenland and Gulf of Alaska landice evolution from an iterated GRACE global mascon solution[END_REF]Watkins et al., 2015). For the CSR mascon results, an averaging kernel with a land mask that extends out 200 km is applied to reduced leakage (Chen et al., 2017b). Similar treatments or corrections of degree-2 zonal term, geocenter motion and GIA effects are also applied in the three mascon solutions. When solving GRACE mascon solutions, the GRACE GAD fields (representing ocean bottom pressure changes, or combined atmospheric and oceanic mass changes) have been added back to the mascon solutions. To correctly quantify ocean mass change using GRACE mascon solutions, the means of the GAD fields over the oceans, which represents mean atmospheric mass changes over the ocean (as ocean mass is conserved in the GAD fields) need to be removed from GRACE mascon solutions. The removal of GAD average over the ocean in GRACE mascon solutions has very minor or negligible effect (of ∼ 0.02 mm yr -1 ) on ocean mass rate estimates, but is important for studying GMSL change at seasonal timescales.

Over the 12-year period (2005-2016), the three GRACE GSM solutions show pretty consistent estimates of ocean mass rate, in the range of 2.3 to 2.5 mm yr -1 . Greater differences are noticed for the mascon solutions. The GSFC mascons show the largest rate of 2.61 mm yr -1 . The CSR and JPL mascon solutions show relatively smaller ocean mass rates of 1.76 and 2.02 mm yr -1 , respectively, over the studied period. Based on the same CSR GSM solutions, the forward modeling and basin kernel estimates agree reasonably well (2.52 vs. 2.44 mm yr -1 ). In addition to the degree-2 zonal term, geocenter motion, and GIA correction, the degree-2, order-1 spherical harmonics of the current GRACE RL05 solutions are affected by the definition of the reference mean pole in GRACE pole tide correction [START_REF] Wahr | The pole tide and its effect on GRACE time-variable gravity measurements: Implications for estimates of surface mass variations[END_REF]. This mean pole correction, excluded in all estimates listed in Table 11 (for fair comparison), is estimated to contribute ∼ -0.11 mm yr -1 to GMSL. How to reduce errors from the different sources plays a critical role in estimating ocean mass change from GRACE time-variable gravity data.

GRACE satellite gravimetry has brought a completely new era for studying global ocean mass change. Owing to the extended record of GRACE gravity measurements (now over 15 years), improved understanding of GRACE gravity data and methods for addressing GRACE limitations (e.g., leakage and low-degree spherical harmonics), and improved knowledge of background geophysical signals (e.g., GIA), GRACE-derived ocean mass rates from different studies in recent years show clearly increased consistency (Table 11). Most of the results agree well with independent observations from satellite altimeter and Argo floats, although the uncertainty ranges are still large. The GRACE Follow-On (FO) mission was launched in May 2018. The GRACE and GRACE-FO together are expected to provide at least over 2 (or even 3) decades of time-variable gravity measurements. Continuous improvements of GRACE data quality (in future releases) and background geophysical models are also expected, which will help improve the accuracy GRACE observed ocean mass change.

For the sea-level budget assessment over the GRACE period, we use the ensemble mean.

Sea-level budget results

In Sect. 2, we have presented the different terms of the sealevel budget equation, mostly based on published estimates (and in some cases, from their updates). We now use them to examine the closure of the sea-level budget. For all terms, we only consider ensemble mean values.

3.1 Entire altimetry era (1993-present)

Trend estimates over 1993-present

Because it is now clear that the GMSL and some components are accelerating (e.g., [START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF], we propose to characterize the long-term variations of the time series by both a trend and an acceleration. We start by looking at trends. Table 12 gathers the trends estimated in Sect. 2. The end year is not always the same for all components (see Sect. 2). Thus the word "present" means either 2015 or 2016 depending on the component. As no trend estimate is available for the entire altimetry era for the terrestrial water storage contribution, we do not consider this component. The residual trend (GMSL minus sum of components trend) may then provide some constraint on the TWS contribution.

Results presented in Table 12 are discussed in detail in Sect. 4.

Acceleration

The GMSL acceleration estimated in Sect. 2.2 using Ablain et al.'s (2017b) TOPEX-A drift correction amounts to 0.10 mm yr -2 for the 1993-2017 time span. This value is in good agreement with the Nerem et al. ( 2018) estimate (of 0.084 ± 0.025 mm yr -2 ) over nearly the same period, after removal of the interannual variability of the GMSL. In [START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF], acceleration of individual components are also estimated as well as acceleration of the sum of components. The latter agrees well with the GMSL acceleration. Here we do not estimate the acceleration of the component ensemble means because time series are not always available. We leave this for a future assessment.

GRACE and Argo period (2005-present)

Sea-level budget using GRACE-based ocean mass

If we consider the ensemble mean trends for the GMSL, thermosteric and ocean mass components given in Sects. 2.2, 2.3 and 2.9 over 2005-present, we find agreement (within error bars) between the observed GMSL (3.5 ± 0.2 mm yr -1 ) and the sum of Argo-based thermosteric plus GRACE-based ocean mass (3.6 ± 0.4 mm yr -1 ) (see Table 13). The residual (GMSL minus sum of components) trend amounts to -0.1 mm yr -1 . Thus in terms of trends, the sea-level budget appears closed over this time span within quoted uncertainties.

Trend estimates over 2005-present from estimates of individual contributions

Table 13 gathers trends of individual components of the sealevel budget over 2005-present, as well as the trend of the sum of components and residuals (GMSL minus sum of components). As for the longer period, ensemble mean values are considered for each component.

As for Table 12, the results presented in Table 13 are discussed in detail in Sect. 4.

Year-to-year budget over 2005-present using GRACE-based ocean mass

We now examine the year-to-year sea-level and mass budgets. Table 14 provides annual mean values for the ensemble mean GMSL, GRACE-based ocean mass and Argo-based thermosteric component. The components are expressed as anomalies and their reference is arbitrary. So to compare with the GMSL, a constant offset for all years was applied to the thermosteric and ocean mass annual means. The reference year (where all values are set to zero) is 2003.

Figure 14 shows the sea-level budget over 2005-2015 in terms of an annual bar chart using values given in Table 14. It compares for years 2005 to 2016 the annual mean GMSL (blue bars) and annual mean sum of thermosteric and GRACE-based ocean mass (red bars). Annual residuals are also shown (green bars). These are either positive or negative depending on the years. The trend of these annual residuals is estimated to be 0.135 mm yr -1 .

In Fig. 15 is also shown the annual sea-level budget over 2005-2015 but now using the individual components for the mass terms. As we have no annual estimates for TWS, we ignore it, so that the total mass includes only glaciers, Greenland and Antarctica. The annual residuals thus include the TWS component in addition to the missing contributions (e.g., deep ocean warming). For years 2006 to 2011, the residuals are negative, an indication of a negative TWS to sea level as suggested by GRACE results [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF][START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF]. But as of 2012, the residuals become

Discussion

The results presented in Sect. 2 for the components of the sea-level budget are based on syntheses of the recently pub- components, allowing us to also look at the mass budget. For TWS, as discussed in Sect. 2.7, GRACE provides a negative trend contribution to sea level over the last decade (i.e., increase on water storage on land) attributed to internal natural variability [START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF], unlike hydrological models that lead to a small (possibly not significantly different from zero) positive contribution to sea level over the same period. Assuming that GRACE observations are perfect, such discrepancies could be attributed to the inability of models to correctly account for uncertainties in meteorological forcing and inadequate modeling of soil storage capacity (see discussion in Sect. 2.7). However, when looking at the sea-level budget over the GRACE time span and using the GRACEbased TWS, we find a rather large positive residual trend (> 0.5 mm yr -1 ) that needs to be explained. Since GRACEbased ocean mass is supposed to represent all mass terms, one may want to attribute this residual trend to an additional contribution of the deep ocean to the abyssal contribution already taken into account here, but possibly underestimated because of incomplete monitoring by current observing systems. If such a large positive contribution from the deep ocean (meaning ocean warming) is real (which is unlikely, given the high implied heat storage), this has to be confirmed by independent approaches, e.g., using ocean reanalysis, and eventually model-based and top-of-the-atmosphere estimates of the Earth energy imbalance.

In addition to mean trends over the period, we also looked at the annual budget for all years, starting in 2005. For most components, annual mean values are provided during the Argo-GRACE era, except for the terrestrial water storage component. However, the sea-level budget based on GRACE ocean mass (plus ocean thermal expansion; Fig. 14) includes the TWS contribution. As shown in Fig. 14, yearly residuals are small, suggesting near closure of the sea-level budget. The residual trend amounts to 0.13 mm yr -1 . It could be interpreted as an additional deep ocean contribution not accounted by the SIO estimate (see Sect. 2.3). However, when looking at Fig. 14, we note that yearly residuals are either positive or negative, an indication of interannual variability that can hardly be explained by a deep ocean contribution. The residual trend derived from the difference (GMSL minus sum of components) (Table 13) amounts -0.1±0.3 mm yr -1 , suggesting a sea-level budget closed within 0.3 mm yr -1 over 2005-present, with no substantial deep ocean contribution.

Figure 16 compares GRACE ocean mass to the sum of mass components (excluding TWS, for the reasons mentioned above). In principle, this mass budget may provide a constraint on the TWS contribution. The corresponding residual trend amounts to 0.14 mm yr -1 over the GRACE period, a value that disagrees with the above-quoted GRACEbased TWS estimates. However, it is worth noting that the GRACE-based TWS trend is very dependent on the considered time span because of the strong interannual variability; a recent study by [START_REF] Palanisamy | Land water storage changes over world river basins from GRACE and global hydrological models[END_REF], based on 347 land river basins, found GRACE-based TWS trend of zero over 2005-2015. Given the remaining data uncertainties, any robust conclusion can hardly be reached so far. That being said, more work is needed to clarify the sign discrepancy between GRACE-based and model-based TWS estimates.

Concluding remarks

As mentioned in the introduction, the global mean sea-level budget has been the subject of numerous previous studies, including successive IPCC assessments of the published literature. What is new in the effort presented here is that it involves the international community currently studying present-day sea level and its components. Moreover, it relies on a large variety of datasets derived from different spacebased and in situ observing systems. The near closure of the sea-level budget, as reported here over the GRACE and Argo era, suggests that no large systematic errors affect these independent observing systems, including the satellite altimetry system. Study of the sea-level budget allows improved understanding of the different processes causing sea-level rise, such as ocean warming and land ice melt. When accuracy in-creases, it will offer an integrative view of the response of the Earth system to natural and anthropogenic forcing and internal variability, and provide an independent constraint on the current Earth energy imbalance. Validation of climate models against observations is another important application of this kind of assessment (e.g., Slangen et al., 2017).

However, important uncertainties still remain, which affect several terms of the budget; for example the GIA correction applied to GRACE data over Antarctica or the net land water storage contribution to sea level. The latter results from a variety of factors but is dominated by groundwater pumping and natural climate variability. Both terms are still uncertain and accurately quantifying them remains a challenge.

Several ongoing international projects related to sea level should provide, in the near future, improved estimates of the components of the sea-level budget. This is the case, for example, of the ice sheet mass balance intercomparison exercise (IMBIE, second assessment), a community effort supported by NASA (National Aeronautics and Space Administration) and ESA, dedicated to reconciling satellite measurements of ice sheet mass balance (The IMBIE Team, 2018). This is also the case for the ongoing ESA Sea Level Budget Closure project [START_REF] Horwath | ESA Climate Change Initiative (CCI) Sea Level Budget Closure (SLBC_cci) Sea Level Budget Closure Assessment Report D3.1[END_REF] that uses a number of space-based essential climate variables (ECVs) reprocessed during the last few years in the context of the ESA Climate Change Initiative project. The recently launched GRACE follow-on mission will lengthen the current mass component time series, with hopefully increased precision and resolution. Finally, the deep Argo project, still in an experimental phase, will provide important information on the deep ocean heat content in the coming years. Availability of this new dataset will provide new insights into the total thermosteric component of the sea-level budget, allowing other missing or poorly known contributions to be constrained from the evaluation of the budget.

The sea-level budget assessment discussed here essentially relies on trend estimates. But annual budget estimates have been proposed for the first time over the GRACE-Argo era. It is planned to provide updates of the global sea-level budget every year, as done for more than a decade for the global carbon budget (Le Queré et al., 2018). In the next assessments, updates of all components will be considered, accounting for improved evaluation of the raw data, improved processing and corrections, use of ocean reanalysis, etc. The need for additional information where gaps exist should also be considered. As a closing remark, study of the sea-level budget in terms of time series and not just trends, as done here, will be required.

level budget (i.e., altimetry-based sea level, tide gauges, thermal expansion, glaciers, Greenland, Antarctica, terrestrial water storage, glacial isostatic adjustment, ocean mass from GRACE). Thanks to the team leaders (in alphabetic order) MA, JB, NC, JC, CD, SJ, JTR, KvS, GS, IV and RvdW, who interacted with their team members, collected all needed information, provided a synthesized assessment of the literature and when needed, updated the published results. The coordinators AC and BM collected those materials and prepared a first draft of the paper, but all authors contributed to its refinement and to the discussion of the results. Special thanks are addressed to JB, EB, GC, JC, GJ (PMEL Contribution Number 4776), BM, FP, RP and ES for improving the successive versions of the paper, and to HP for providing all figures presented in Sect. 3.

The views, opinions and findings contained in this paper are those of the authors and should not be construed as an official position, policy or decision of the NOAA, US Government or other institutions.

Figure 1 .

 1 Figure 1. Evolution of GMSL time series from six different groups' (AVISO/CNES, SL_cci/ESA, University of Colorado, CSIRO, NASA/GSFC, NOAA) products. Annual signals are removed and 6-month smoothing applied. All GMSL time series are centered in 1993 with zero mean. A GIA correction of -0.3 mm yr -1 has been subtracted from each dataset.

  .084 ± 0.025 mm yr -2[START_REF] Nerem | Climate Change Driven Accelerated Sea Level Rise Detected In The Altimeter Era[END_REF] (noteWatson et al., 2015 found a smaller acceleration after correcting for the instrumental bias over a shorter period up to the end of 2014.). GMSL trends calculated over 10-year moving windows illustrate this acceleration (Fig.3). GMSL trends are close to 2.5 mm yr -1 over 1993-2002 and 3.0 mm yr -1 over 1996-2005. After a slightly smaller trend over 2002-2011, the 2008-2017 trend reaches 4.2 mm yr -1 . Uncertainties (90 % confidence interval) associated with these 10-year trends regularly decrease through time from 1.3 mm yr -1 over 1993-2002 (corresponding to T/P data) to 0.65 mm yr -1 for 2008-2017 (corresponding to Jason-2 and Jason-3 data).

Figure 2 .

 2 Figure 2. Evolution of ensemble mean GMSL time series (average of the six GMSL products from AVISO/CNES, SL_cci/ESA, University of Colorado, CSIRO, NASA/GSFC and NOAA). On the black, red and green curves, the TOPEX-A drift correction is applied respectively based on Ablain et al. (2017b), Watson et al. (2015) and Dieng et al. (2017), and Beckley et al. (2017). Annual signal removed and 6-month smoothing applied; GIA correction also applied. Uncertainties (90 % confidence interval) of correlated errors over a 1-year period are superimposed for each individual measurement (shaded area).

Figure 3 .

 3 Figure 3. Ensemble mean GMSL trends calculated over 10-year moving windows. On the black, red and green curves, the TOPEX-A drift correction is applied respectively based on Ablain et al. (2017b), Watson et al. (2015) and Dieng et al. (2017), and Beckley et al. (2017). Uncorrected GMSL trends are shown by the blue curve. The shaded area represents trend uncertainty over 10-year periods (90 % confidence interval).

Figure 4 .

 4 Figure 4. Left panels: annual mean global mean thermosteric anomaly time series since 1970, from various research groups (color) and for three depth integrations: 0-700 m (top), 700-2000 m (middle) and below 2000 m (bottom). Vertical dashed lines are plotted along 1993 and 2005. For comparison, all time series were offset arbitrarily. Right panels: respective linearly detrended time series for 1993-2015. Black bold dashed line is the ensemble mean and gray shadow bar the ensemble spread (1 standard deviation). Units are millimeters.

Figure 5 .

 5 Figure 5. Left panel: annual mean global mean thermosteric anomaly time series since 2004, from various research groups (color) in the upper 2000 m. A vertical dashed line is plotted along 2005. For comparison, all time series were offset arbitrarily. Right panel: respective linearly detrended time series for 2005-2015. Black bold dashed line is the ensemble mean and gray shadow bar the ensemble spread (1 standard deviation). Units are millimeters.

Figure 6 .

 6 Figure 6. Ensemble mean time series for global mean thermosteric anomaly, for three depth integrations (a) and for 0-2000 m and full depth (b). In the bottom panel, dashed lines are for the 1993-2015 period whereas solid lines are for 2005-2015. Error bars represent the ensemble spread (standard deviation). Units are millimeters.

  the 1993-2015 (a) and 2005-2015 (b) time spans and different depth ranges (including full depth), as well as associated ensemble mean trends. The full depth ensemble mean trend amounts to 1.3 ± 0.4 mm yr -1 over 2005-2015. It is similar to the 1993-2015 ensemble mean trend, suggesting negligible acceleration of the thermosteric component over the altimetry era.

Figure 7 .

 7 Figure 7. Linear rates of global mean thermosteric sea level for depth integrations (x axis), individual estimates and ensemble means, over 1993-2015 (a) and 2005-2015 (b). Ensemble mean rates with a black circle were used in the estimation of the time series described in Sect. 2.3.4. Error bars are standard deviation due to spread of the estimates except for ≥ 2000 m. Units are millimeters per year.

Figure 8 .

 8 Figure 8. Evolution of global glacier ice mass estimates from different studies published over the past 2 decades, based on different observations and methods. The red marks correspond to IPCC reports. We clearly see the most recent publications lead to less scattered results. Note that Antarctica and Greenland peripheral glaciers are taken into account in this figure.

Figure 9 .

 9 Figure 9. Greenland annual mass change from 1993 to 2016. The medium blue region shows the range of estimates from the datasets listed in Table1. The lighter blue region shows the range of estimates when stated errors are included, to provide upper and lower bounds. The dark blue line shows the mean mass trend.

  , as did Martín-Español et al. (2016); Luthcke et al. (2013), Velicogna et al. (2014), and Groh and Horwath (2016) used IJ05-R2 (Ivins et al., 2013). Wouter et al. (2013) used Whitehouse et al. (2012), and Wiese et al. (2016b) used A et al. (2013).

Figure 10 .

 10 Figure 10. Antarctic annual sea-level contribution during 2005 to 2015. The black squares are the mean annual sea level calculated using the GRACE datasets listed in Table6. The darker blue band shows the range of estimates from the datasets. The light blue band accounts for the error in the different GRACE estimates. The brown squares are the annual sea-level contribution calculated using the input-output method (updated fromRignot et al., 2011a); the light brown band is the associated error.

Figure 11 .

 11 Figure 11. Time series of the estimated annual contribution of terrestrial water storage change to global sea level over the period 1900-2014 (rates in millimeters per year of SLE) (modified from Wada et al., 2016).

Figure 12 .

 12 Figure 12. An example of trends in land water storage from GRACE observations, April 2002 to November 2014. Glaciers and ice sheets are excluded. Shown are the global map (gigatons per year), zonal trends and full time series of land water storage (in mm yr -1 SLE). Following methods detailed in Reager et al. (2016), GRACE shows a total gain in land water storage during the 2002-2014 period, corresponding to a sea-level trend of -0.33±0.16 mm yr -1 SLE (modified from[START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF]. These trends include all human-driven and climate-driven processes in Table1and can be used to close the land water budget over the study period.

  . However, different groups have independently computed GIA model solutions based on the Toronto ice history reconstruction, by using different implementations of GIA codes and somehow different Earth models. The most widely used model is the one by Paulson et al. (2007), later updated by A et al. (2013). Both studies use a deglaciation history based on ICE-5G, but differ for the viscosity profile of the mantle: A et al. (2013) use a 3-D compressible Earth with VM2 viscosity profile and a PREM-based elastic structure used by Peltier (2004), whereas Paulson et al. (2007) use an incompressible Earth with self-gravitation, and a Maxwell 1-D multilayer mantle.

Figure 13 .

 13 Figure 13. Difference map between two models of GIA correction to GRACE over land: A et al. (2013) versus Peltier et al. (2015). Units in millimeters per year of SLE.

  positive and on average over 2005-2015, the residual trend amounts to +0.28 mm yr -1 , a value larger than when using GRACE ocean mass.Finally, Fig.16presents the mass budget. It compares annual GRACE-based ocean mass to the sum of the mass components, without TWS as in Fig.15. The residual trend over the 2005-2015 time span is 0.14 mm yr -1 . It may dominantly represent the TWS contribution. From one year to another residuals can be either positive or negative, suggesting important interannual variability in the TWS or even in the deep ocean.

Figure 14 .

 14 Figure 14. Annual sea level (blue bars) and sum of thermal expansion (full depth) and GRACE ocean mass component (red bars). Black vertical bars are associated uncertainties. Annual residuals (green bars) are also shown.

Figure 15 .

 15 Figure 15. Annual global mean sea level (blue bars) and sum components without TWS (full depth thermal expansion + glaciers + Greenland + Antarctica) (red bars). Black vertical bars are associated uncertainties. Annual residuals (green bars) are also shown.

Figure 16 .

 16 Figure 16. Annual GRACE-based ocean mass (red bars) and sum components without TWS (full depth thermal expansion + glaciers + Greenland + Antarctica) (blue bars). Annual residuals (green bars) are also shown.

Table 1 .

 1 TOPEX-A GMSL drift corrections proposed by different studies.

  compared to global mean sea level from satellite altimetry, due to limited and uneven spatial sampling of the global ocean afforded by the tide gauge network. Sea-level rise being non uniform, spatial variability of sea-level measured at tide gauges is evidenced by 2-D reconstruction methods. The most widely used approach is the use of empirical orthogonal functions (EOFs) calibrated with the satellite altimetry data (e.g.,[START_REF] Church | A 20th century acceleration in global sea-level rise[END_REF]. Alternatively,[START_REF] Choblet | Probabilistic surface reconstruction of coastal sea level rise during the twentieth century[END_REF] implemented a Bayesian inference method based on a Voronoi tessellation of the Earth's surface to reconstruct sea level during the 20th century. Considerable uncertainties remain, however, in long-term assessments due to poorly sampled ocean basins such as the South Atlantic, or regions which are significantly influenced by open-ocean circulation (e.g., subtropical North Atlantic)

www.earth-syst-sci-data.net/10/1551/2018/ Earth Syst. Sci. Data, 10, 1551-1590, 2018

Table 2 .

 2 Compilation of available in situ datasets from different originators and/or contributors.The table indicates the time span covered by the data, the depth of integration, as well as the temporal resolution and latitude coverage.

	Product/institution	Period	Depth integration (m)			Temporal resolu-	Reference
								tion/latitudinal
								range
				0-700	700-	0-2000	≥ 2000
					2000
	1	CORA	1993-2016	Y	Y	Y	-	Monthly 60 • S-	http://marine.
								60 • N	copernicus.eu/
									services-portfolio/
									access-to-products/
	2	CSIRO (RSOI) 2004-2017	Y/E (0-300)	Y/E	Y/E	-	Monthly 65 • S-	Roemmich et
								65 • N	al. (2015), Wijffels
									et al. (2016)
	3	CSIRO/	1970-2017	Y/E (0-300)	-	-	-	Yearly (3-year	Domingues et al.
		ACECRC/						running mean)	(2008), Church et
		IMAS-UTAS						65 • S-65 • N	al. (2011)
	4	ICCES	1970-2016	Y/E (0-300)	Y/E	Y/E	-	Yearly 89 • S-	Cheng et al. (2017)
								89 • N
	5	ICDC	1993-2016	Y (1993)	-	Y (2005)	-	Monthly	Gouretski and
									Koltermann (2007)
	6	IPRC	2005-2016	-	-	Y	-	Monthly	http://apdrc.
									soest.hawaii.edu/
									projects/argo (last
									access: 22 August
									2018)
	7	JAMSTEC	2005-2016	-	-	Y	-	Monthly	Hosoda et al.
									(2008)
	8	MRI/JMA	1970-2016	Y/E (0-300)	Y/E	Y/E	-	Yearly 89 • S-	Ishii et al.
			(rel. to 1961-					89 • N	(2009, 2017)
			1990 averages)
	9	NCEI/NOAA	1970-2016	Y/E	Y/E	Y/E	-	Yearly 89 • S-	Levitus et al.
								89 • N	(2012)
	10 SIO	2005-2016	-	-	Y	-	Monthly	Roemmich and
									Gilson (2009)
	11 SIO	1990-2010 (as	-	-	-	Y/E	Linear trend 89 • S-	Purkey and
		(Deep-	of Jan 2018)					89 • N, as an aggre-	Johnson (2010)
		abyssal)						gation of 32 deep
								ocean basins

Table 3 .

 3 Glacier contribution to sea level; all data are in millimeters per year of SLE.

		1993-2016	2005-2016
		mm yr -1 SLE mm yr -1 SLE
	G16		0.70 ± 0.070 a
	M17	0.68 ± 0.032	0.80 ± 0.048
	C17	0.63 ± 0.070 0.75 ± 0.070 b
	L17		0.84 ± 0.640 c
	M17-G		0.61 ± 0.070 d

a The time period of G16 is 2002-2014. b The time period of C17 is 2003-2009. c The time period of L17 is 2003-2009. d The time period of M17-G is 2002/2005-2013/2015 because this value is an average of different estimates.

ferent time periods. However, because of the probable low variability of global annual glacier changes, compared to other components of the sea-level budget, averaging trends for slightly different time periods is appropriate.

Table 4 .

 4 Datasets considered in the Greenland mass balance assessment, as well as covered time span and type of observations.

	Reference	Time period Method
	Update from

Table 5 .

 5 Annual time series of Greenland mass change (GT yr -1 , negative values mean decreasing mass). The mass is calculated as the weighted mean based on the stated error value for each year. The error for each year is calculated as the mean of all stated 1σ errors divided by sqrt(N ) where N is the number of datasets available for that year, assuming that the errors are uncorrelated.

	The stan-

Table 6 .

 6 Datasets considered in this assessment of the Antarctica mass change, and associated trends for the 2005-2015 and 1993-2015 expressed in millimeters per year of SLE. Positive values mean positive contribution to sea level (i.e., sea-level rise).

	2005-2015	1993-2015
	trend (mm yr -1 ) trend (mm yr -1 )

Table 7 .

 7 Annual sea-level contribution from Antarctica during 2005-2015 from GRACE and input-output method (IOM) calculated as described above and expressed in millimeters per year of SLE. Also shown is the mean of the estimate from the two methods; associated errors are the mean of the two estimated errors. Positive values mean positive contribution to sea level (i.e., sea-level rise).

		GRACE	IOM	Mean
		(mm yr -1 )	(mm yr -1 )	(mm yr -1 )
	Year	SLE	SLE	SLE
	2005	-0.34 ± 0.47 -0.51 ± 0.16 -0.42 ± 0.31
	2006	0.04 ± 0.36	0.23 ± 0.16	0.14 ± 0.26
	2007	0.58 ± 0.42	0.68 ± 0.16	0.63 ± 0.29
	2008	0.22 ± 0.29	0.35 ± 0.16	0.29 ± 0.22
	2009	0.09 ± 0.26	0.42 ± 0.16	0.26 ± 0.21
	2010	0.74 ± 0.30	0.59 ± 0.16	0.67 ± 0.23
	2011	0.15 ± 0.39	0.30 ± 0.16	0.23 ± 0.27
	2012	0.25 ± 0.30	0.64 ± 0.16	0.44 ± 0.23
	2013	0.63 ± 0.38	0.67 ± 0.16	0.65 ± 0.27
	2014	0.78 ± 0.46	0.69 ± 0.16	0.73 ± 0.31
	2015	0.09 ± 0.77	0.50 ± 0.16	0.29 ± 0.46
	Average estimate
	2005-2015	0.38 ± 0.06	0.46 ± 0.05	0.42 ± 0.06

Table 8 .

 8 Estimates of TWS components due to human intervention and net TWS based on hydrological models and GRACE.

	2002-2014/15
	(mm yr -1 ) SLE
	(positive values mean

Table 3

 3

		-0.30 ± 0.02	Average of three groups of four values obtained by variants of the analysis
	)	-0.29 ± 0.03	procedure, using ICE-5G(VM2), over a global ocean, in the range of latitudes
		-0.28 ± 0.02	66 • S to 66 • N and 60 • S to 60 • N, respectively.
	Tamisiea (2011)	-0.15 to -0.45	Simple average over the oceans for a range of estimates obtained varying the
	(Fig. 2)	-0.20 to -0.50	Earth model parameters, over a global ocean and between latitudes 66 • S and
			66 • N.
	Huang (2013)	-0.26 ± 0.07	Average from an ensemble of 14 GIA models over a global ocean and between
	(Table 3.6)	-0.27 ± 0.08	latitude from 66 • S to 66 • N.
	Spada (2017)	-0.32 ± 0.08	Based on four runs of the sea-level equation solver SELEN (Spada and Stoc-
	(Table 1)		chi, 2007) using model ICE-5G(VM2), with different assumptions in solving
			the SLE.
	(b) GIA contribution to GRACE mass rate of change over the oceans
	Reference	GIA	Notes
		(mm yr -1 SLE)
	Peltier (2009)	-1.60 ± 0.30	Average of values from 12 corrections for variants of the analysis procedure,
	(Table 3)		using ICE-5G (VM2).
	Chambers et	-1.45 ± 0.35	Average over the oceans for a range of estimates produced by varying the Earth
	al. (2010)		models.
	(Table 1)

Table 9 .

 9 Continued. (d) GIA contribution to GRACE mass rate of the ice sheets

	Reference	Greenland GIA	Notes
		(Gt yr -1 )
	Simpson et	-3 ± 12 m	Thermodynamic sheet/solid Earth model, 1-D (uncoupled); constrained by ge-
	al. (2009) r		omorphology; inversion results in Sutterley et al. (2014).
	Peltier (2009)	-4 d	Ice load reconstruction/solid Earth model, 1-D (ICE-5G/similar to VM2);
	(ICE-5G) b		Greenland component of ICE-5G (13 Gt yr -1 ) + Laurentide component of
			ICE-5G (-17 Gt yr -1 ); inversion results in Khan et al. (2016), Discussion.
	Khan et al. (2016)	15 ± 10 f	Ice load reconstruction/solid Earth model, 1-D (uncoupled); constrained with
	(GGG-1D) a		geomorphology and GPS; Greenland component (+32 Gt yr -1 ) + Laurentide
			component of ICE-5G (-17 Gt yr -1 ); inversion results in Khan et al. (2016),
			Discussion.
	Fleming and	3 d	Ice load reconstruction/solid Earth model, 1-D (uncoupled); constrained with
	Lambeck (2004) a		geomorphology; Greenland component (+20 Gt yr -1 ) + Laurentide compo-
	(Green1)		nent of ICE-5G (-17 Gt yr -1 ); inversion in Sasgen et al. (2012, Supplement).
	Wu et al. (2010) b	-69 ± 19 m	Joint inversion estimate based on GPS, satellite laser ranging, and very long
			baseline interferometry, and bottom pressure from ocean model output; inver-
			sion results in Sutterley et al. (2014).
	Reference	Antarctica GIA	Notes
		(Gt yr -1 )
	Whitehouse et al.	60 f	Thermodynamic sheet/solid Earth model, 1-D (uncoupled); constrained by
	(2012)		geomorphology; inversion results in Shepherd et al. (2012), Supplement
	(W12a) a		(Fig. S8).
	Ivins et al. (2013)	40-65 f	Ice load reconstruction/solid Earth model, 1-D; constrained by geomorphol-
	(IJ05_R2) a		ogy and GPS uplift rates; Ivins et al. (2013); inversion results in Shepherd et
			al. (2012), Supplement (Fig. S8).
	Peltier (2009)	140-180 f	Ice load reconstruction/solid Earth model ICE-5G(VM2); constrained by geo-
	(ICE-5G) b		morphology; inversion results in Shepherd et al. (2012), Supplement (Fig. S8).
	Argus et al. (2014)	107 f	Ice load reconstruction/solid Earth model ICE-6G(VM5a); constrained by ge-
	(ICE-6G) b		omorphology and GPS; theory recently corrected by Purcell et al. (2016); in-
			version results in Argus et al. (2014), conclusion 7.8.
	Sasgen et al. (2017)	55 ± 22 f	Joint inversion estimate based on GRACE, altimetry, GPS and viscoelastic
	(REGINA) a		response functions; lateral heterogeneous Earth model parameters; inversion
			results in Sasgen et al. (2017), Table 1.
	Gunter et al. (2014)	ca. 64 ± 40 a	Joint inversion estimate based on GRACE, altimetry, GPS and regional cli-
	(G14) a	(multimodel	mate model output; conversion of uplift to mass using average rock density;
		uncert.)	inversion results in, Gunter et al. (2014) Table 1.
	Martín-Español et al.	55 ± 8	Joint inversion estimate based on GRACE, altimetry, GPS and regional cli-
	(2016) (RATES) a	45 ± 7*	mate model output; inversion results in Sasgen et al. (2017), * is improved for
			GIA of smaller spatial scales; inversion results in Martin-Español et al. (2016),
			Fig. 6.

a Regional model. b Global model. c Mascon inversion. d Forward modeling inversion. e Averaging kernel inversion. f Inversion method not specified.

Table 10 .

 10 Recently published (since 2013) estimates of GRACE-based ocean mass rates (GIA corrected). Most of the listed studies use either the A13(A et al., 2013) or Paulson07[START_REF] Paulson | Inference of mantle viscosity from GRACE and relative sea level data[END_REF] GIA model.

	Ocean mass

Table 11 .

 11 Ocean mass trends (in mm yr -1 ) estimated from GRACE for the period January 2005-December 2016 (the GSFC mascon solutions cover up to July 2016). The uncertainty is based on 2 times the sigma of least-squares fitting.

	Ocean mass

Table 12 .

 12 Trend estimates for individual components of the sealevel budget, sum of components and GMSL minus sum of components over 1993-present. Uncertainties of the sum of components and residuals represent rooted mean squares of components errors, assuming that errors are independent.

		Trends (mm yr -1 )
	Component	1993-present
	1.GMSL (TOPEX-A drift corrected)	3.07 ± 0.37
	2. Thermosteric sea level (full depth)	1.3 ± 0.4
	3. Glaciers	0.65 ± 0.15
	4. Greenland	0.48 ± 0.10
	5. Antarctica	0.25 ± 0.10
	6. TWS	/
	7. Sum of components	2.7 ± 0.23
	(without TWS →2.+3.+4.+5.)
	8. GMSL minus sum of components	0.37 ± 0.3
	(without TWS)

Table 13 .

 13 Trend estimates for individual components of the sea-level budget, sum of components and GMSL minus sum of components over 2005-present. TWS from GRACE (mean of[START_REF] Reager | A decade of sea level rise slowed by climate-driven hydrology[END_REF][START_REF] Scanlon | Global models under-estimate large decadal declining and rising water storage trends relative to GRACE satellite data[END_REF] -0.27 ± 0.15 7. Sum of components (2.+3.+4.+5.+6.)2.95 ± 0.21 8. Sum of components (thermosteric full depth + GRACE-based ocean mass)3.6 ± 0.4 9. GMSL minus sum of components (including GRACE-based TWS → 2.+3.+4.+5.+6.) 0.55 ± 0.3 10.GMSL minus sum of components (without GRACE-based TWS → 2.+3.+4.+5.) 0.28 ± 0.2 11. GMSL minus sum of components (thermosteric full depth + GRACE-based ocean mass) -0.1 ± 0.3

	Trend (mm yr -1 )

Table 14 .

 14 Annual mean values for the ensemble mean GMSL and sum of components (GRACE-based ocean mass and Argo-based thermosteric, full depth). Constant offset applied to the sum of components. The reference year (where all values are set to zero) is 2003.

		Ensemble	Sum of	GMSL
		mean GMSL components	minus sum of
	Year	(mm)	(mm) components (mm)
	2005	7.00	8.78	-0.78
	2006	10.25	10.78	-0.53
	2007	10.51	11.35	-0.85
	2008	15.33	15.07	0.25
	2009	18.78	18.88	-0.10
	2010	20.64	20.53	0.11
	2011	20.91	21.38	-0.48
	2012	31.10	29.33	1.77
	2013	33.40	33.87	-0.47
	2014	36.65	36.22	0.43
	2015	46.34	45.69	0.65

www.earth-syst-sci-data.net/10/1551/2018/ Earth Syst. Sci. Data, 10, 1551-1590, 2018

Earth Syst. Sci.Data, 10, 1551Data, 10, -1590Data, 10, , 2018 www.earth-syst-sci-data.net/10/1551/2018/

Data availabilityThe data sets used in this study are freely available at https://doi.org/10.17882/54854. We provide annual mean Earth Syst. Sci.Data, 10, 1551Data, 10, -1590Data, 10, , 2018 www.earth-syst-sci-data.net/10/1551/2018/

Author contributions. This community assessment was initiated by AC and BM as a contribution to the Grand Challenge "Regional Sea Level and Coastal Impacts" of the World Climate Research Programme (WCRP). The results presented in this paper were prepared by nine different teams dedicated to the various terms of the seawww.earth-syst-sci-data.net/10/1551/2018/ Earth Syst. Sci.Data, 10, 1551Data, 10, -1590Data, 10, , 2018

Acknowledgements. We are grateful to the anonymous reviewer for his/her thorough comments that helped to improve the paper.

Edited by: Giuseppe M. R. Manzella Reviewed by: one anonymous referee

Competing interests. The authors declare that they have no conflict of interest.