Sea level change is one of the major consequences of climate change and is projected to affect coastal communities around the world. Here, global mean sea level (GMSL) change estimated by 12 climate models from phase 5 of the World Climate Research Programme’s Climate Model Intercomparison Project (CMIP5) is compared to observational estimates for the period 1900–2015. Observed and simulated individual contributions to GMSL change (thermal expansion, glacier mass change, ice sheet mass change, landwater storage change) are analyzed and compared to observed GMSL change over the period 1900–2007 using tide gauge reconstructions, and over the period 1993–2015 using satellite altimetry estimates. The model-simulated contributions explain 50% ± 30% (uncertainties 1.65σ unless indicated otherwise) of the mean observed change from 1901–20 to 1988–2007. Based on attributable biases between observations and models, a number of corrections are proposed, which result in an improved explanation of 75% ± 38% of the observed change. For the satellite era (from 1993–97 to 2011–15) an improved budget closure of 102% ± 33% is found (105% ± 35% when including the proposed bias corrections). Simulated decadal trends increase over the twentieth century, both in the thermal expansion and the combined mass contributions (glaciers, ice sheets, and landwater storage). The mass components explain the majority of sea level rise over the twentieth century, but the thermal expansion has increasingly contributed to sea level rise, starting from 1910 onward and in 2015 accounting for 46% of the total simulated sea level change.
Sea level change is one of the most well-known consequences of climate change, affecting coastal communities and ecosystems worldwide. Changes in sea level are the result of changes in different components of the climate system: the ocean, the land, the atmosphere, and the cryosphere. Therefore, the study of sea level change is like a jigsaw puzzle, requiring a complete and integrative view of the climate system on a range of spatial and temporal scales (Fig. 1; Church et al. 2013a; Cazenave et al. 2017).
On a global scale, there are several processes that contribute to long-term sea level change. A major process is the thermal expansion and contraction of the ocean water, caused by density changes due to temperature changes. Ocean salinity changes, which also cause density variations, are negligible on a global scale (Lowe and Gregory 2006). Second, changes in land ice mass of glaciers and ice sheets (surface mass balance and ice dynamical processes) contribute to sea level change. In addition, changes in landwater storage due to human activity should be taken into account, most importantly reservoir building and groundwater extraction (Fig. 1, left). Natural changes in landwater storage, such as those resulting from snow cover, surface water, or soil moisture, are large on seasonal time scales and are found to be important up to decadal time scales (Reager et al. 2016) but are assumed to be small on century time scales (Wada et al. 2016; Church et al. 2013a).
This paper focuses on the global mean contributions to twentieth-century sea level change. However, it is important to note that sea level changes vary spatially, and additional processes come into play on smaller spatial and temporal scales (Fig. 1; e.g., Slangen et al. 2014a). At a regional scale this includes, for instance, changes in the ocean circulation (e.g., due to wind stress) or local ocean water density variations (both temperature and salinity driven). Other causes are changes in Earth’s gravitational field and vertical land motion due to redistribution of mass between land and the ocean—both present-day ice mass change and glacial isostatic adjustment in response to ice mass loss after the Last Glacial Maximum (Spada 2017). See Part II of this paper by Meyssignac et al. (2017, hereafter M17) for a detailed discussion of these regional sea level change patterns. Even smaller scales are required when studying coastal sea level change and potential changes in sea level extremes (Fig. 1, right; e.g., Cannaby et al. 2016).
The main goal of this paper is to evaluate model simulations of the physical processes that contribute to global mean sea level (GMSL) change and to establish how well they can explain the observed GMSL change in the twentieth century. A comparison of the observed and simulated regional distribution of sea level change at individual tide gauge locations is presented in the accompanying paper (M17).
In recent decades, one of the major questions in sea level research has been the closure of the twentieth-century (global mean) sea level budget (Munk 2002; Church et al. 2011; Gregory et al. 2013a; Jevrejeva et al. 2017). There are several parts to this question: first, the sum of the observations for the individual components (Fig. 1, left) tends to underestimate the total observed change inferred from tide gauges, raising the issue of different types of uncertainties in the observations. Second, the sum of model-based sea level contributions also tends to underestimate the total observed change for the twentieth century, probably because of uncertainties in both the models and the observations. Past observations of sea level change relying on tide gauges are spatially and temporally sparse, difficult to quality-control, and biased to the Northern Hemisphere. They may therefore not be fully representative of GMSL (Thompson et al. 2016) until satellite data started to become available in the early 1990s. The model-based contributions, on the other hand, may not fully account for all climate variability, such as multidecadal variations in the ocean, the delayed response of glaciers and ice sheets to externally driven climate change, and imperfections in the applied forcings and/or the model responses.
Recently, several large community efforts (Church et al. 2011; Gregory et al. 2013a) have explored a wide range of observational estimates of all contributions to sea level change, and were able to close the twentieth-century observational budget within uncertainties. In this paper, we focus on the second part of the problem, namely reconciling the simulated model estimates with observed GMSL change. This is an important topic, as a better understanding of and ability to model past sea level change increases confidence in the models’ ability to project future sea level changes.
We build on and extend previous work that compared model estimates to observations of GMSL rise (e.g., Church et al. 2013b; Slangen et al. 2016). With respect to previous work, the time period has been extended to 2015, allowing for the longest possible time series of satellite observations to be included. The model estimates of the individual contributions to sea level rise have been updated and now use more recent information. For instance, the glacier model now uses a more recent and reliable version of the Randolph Glacier Inventory than in Church et al. (2013b) and the estimate of groundwater extraction is updated following recent publications (Döll et al. 2014; Wada et al. 2016).
We present a consistent model dataset that uses the same set of climate models to estimate the glacier contribution, the surface mass balance of the Greenland and Antarctic ice sheets, and the ocean thermal expansion. This is a different approach from that in the work of Church et al. (2013b), where thermal expansion and the glacier contribution were simulated with different climate model ensembles. In addition, we provide separate estimates for contributions of surface mass balance (based on CMIP5 climate model output) and ice dynamics of the Greenland and Antarctic ice sheets, which in Church et al. (2013b) were estimated with a constant of 0.1 mm yr−1. Here, the use of a consistent set of models for an increased number of contributions allows us to evaluate the individual model simulations rather than the model ensemble alone.
The simulated sea level changes are completed by extended observation-based estimates of ice sheet dynamical processes and landwater storage changes. These estimates are compared to tide gauge reconstructions for the period 1900–2007 and to satellite altimetry for the period 1993–2015. This allows us to quantify and study the acceleration in sea level rise from the early 1990s onward.
We also discuss potential additional contributions that should be included in the comparison between model-based estimates and observations. These contributions are all based on observational evidence and are currently absent or poorly represented in climate models. In addition, observational estimates are particularly sparse due to limited accessibility and the large surface area of the ice sheets. While it has been generally assumed that the historical contributions of the ice sheets are small compared to thermal expansion and glacier mass change, a recent study suggests a more substantial contribution of 25 ± 9 mm (1σ) from the Greenland ice sheet over the twentieth century arising from internal climate variability and/or the response to increased radiative forcing (Kjeldsen et al. 2015).
We first present the observational estimates of GMSL and the climate models and experiments that we use in our comparisons (section 2). The individual contributions to GMSL change, both simulated and observed, are discussed in section 3. In section 4 we combine the model-simulated estimates and compare them to tide gauge reconstructions (1900–2007) and satellite measurements (1993–2015).
a. Sea level observations
For most of the twentieth century, tide gauges are the only source providing sea level measurements. This is not ideal, as tide gauges (particularly the ones going further back in time) are distributed unevenly around the world and are confined to coastal locations. It is therefore not correct to take a simple average, as this would result in a biased global mean sea level record (Thompson et al. 2016). Instead, several methods have been developed to reconstruct the sea level field and obtain an estimate of GMSL change back to 1900 using tide gauges. Here we use four different reconstructions. The first one is from Church and White (2011), an update of their previous work (Church and White 2006) in which they reconstructed historical sea levels by deriving empirical orthogonal functions (EOFs) from satellite data and used the EOFs to find the best fit to the tide gauge records. These data were downloaded from www.cmar.csiro.au/sealevel/ and provide yearly values from 1880 to 2013. Ray and Douglas (2011) provide the second reconstruction time series (downloaded from www.psmsl.org) with yearly values for 1900–2007. They use a method similar to that of Church et al. but with differences in details of the methodology. With respect to Church et al. (2013b) we add two more recent tide gauge reconstructions, constructed using different methods. The reconstruction from Jevrejeva et al. (2014) (downloaded from www.psmsl.org) provides monthly values for 1807–2010. This reconstruction uses the “virtual station” method, where 1277 tide gauge records are divided into 14 regions to overcome a geographical bias in the data, and these regions are then averaged to obtain the global mean change. Finally, the reconstruction from Hay et al. (2015) (downloaded from http://www.nature.com/nature/journal/v517/n7535/source_data/nature14093-f2.xls; their Fig. 2 source data) provides yearly values for 1900–2010. These authors have used Kalman smoothing to process the spatially and temporally sparse tide gauge data and combine this with spatial fingerprints from the different sea level contributions.
Since 1993, satellite altimetry has provided more globally complete observations of sea level change (covering at least 66°S–66°N and up to 82°N/S for Envisat). Here we use three different global mean time series, which all include measurements from the consecutive TOPEX/Poseidon, Jason-1, and Jason-2 satellite missions, but with different choices in instrumental and geophysical corrections and in the algorithms used to compute the GMSL time series (Masters et al. 2012; Henry et al. 2014). We use a time series from the European Space Agency (ESA) sea level climate change initiative (CCI) (Ablain et al. 2015; http://www.esa-sealevel-cci.org), which in addition includes data from the ERS-1 and -2 and Envisat satellite missions. The data have been corrected for the seasonal signal (annual and semiannual), inverse barometer (IB), and glacial isostatic adjustment (GIA; −0.3 mm yr−1; Tamisiea 2011). The second time series, from CSIRO (Church and White 2011; www.cmar.csiro.au/sealevel/), is also corrected for the seasonal signal, IB, and GIA. The third time series is from Watson et al. (2015), who suggested a bias correction for an inferred drift in the early part of the altimeter record, based on the difference between the CSIRO time series and tide gauge observations (where the tide gauges were corrected for vertical land movement using GPS data), which mainly influences the TOPEX part of the satellite record. This bias correction is a different approach than in the other two satellite time series (CSIRO and ESA-CCI), both of which are completely independent from tide gauge measurements.
b. Model data
To model the twentieth-century contributions to sea level change, we mainly use data from phase 5 of the Climate Model Intercomparison Project (CMIP5; Taylor et al. 2012). These data are used to evaluate ocean thermal expansion and to compute the contributions from glaciers and ice sheet surface mass balance using offline models. To cover the period 1900–2015 (as all four reconstructions provide data from 1900 onward), the historical model simulations (1850–2005) were extended using the representative concentration pathway (RCP) projections from the CMIP5 database (Moss et al. 2010), where we used the RCP8.5 scenario out of the four scenarios available (RCP2.6, 4.5, 6.0, and 8.5). This RCP projects a radiative forcing increase of >8.5 W m−2 in 2100 relative to preindustrial conditions. The choice for the RCP8.5 scenario was based on availability, as most models are available for RCP4.5 and 8.5, and less for RCP2.6 and 6.0. However, this is not critical to our results, as the different scenarios only start to diverge significantly after the year 2030 (Church et al. 2013a).
The following CMIP5 variables were used: “zostoga” (global mean ocean thermal expansion), “to” (ocean temperature), “so” (ocean salinity), “tas” and “ta” (air temperature at surface and at 600 hPa, respectively), “pr” (precipitation), “prsn” (snowfall), and “evspsbl” (evaporation). The ocean variables (zostoga, tos, so) were dedrifted by computing linear fits on the full control run (which is forced by nonevolving preindustrial conditions) and then subtracting the time-corresponding part of the fit of each models’ preindustrial control simulation from the historical simulation. The use of a linear fit ensures that the drift is removed but no physical trends are, which is particularly relevant for the regional change in the companion paper (M17). In addition, Hobbs et al. (2016) showed that the use of a linear fit is an adequate dedrifting method particularly for global quantities, and there is no additional gain from using higher-order fits. Dedrifting of ocean variables is necessary to account for spurious trends that result from the (deep) ocean not being in equilibrium with the forcing conditions (e.g., Sen Gupta et al. 2013) and/or to correct for imperfect representation of the global energy budget (Hobbs et al. 2016). We present annual time series unless stated otherwise. Details on how the CMIP5 variables were used to compute each of the sea level contributions will be discussed in section 3.
To obtain a consistent set of models for all contributions, the selection of atmosphere–ocean general circulation models (AOGCMs; Table 1) is based on the availability of all of the abovementioned required variables, resulting in a set of 12 models. From each model, the first available realization (r1i1p1) was used. Comparisons for the thermal expansion contribution show that this set of models and realizations is representative of the larger CMIP5 ensemble as used in Slangen et al. (2014a, 2015). The same set of 12 models is used in the M17 companion paper, which focuses on regional sea level patterns rather than the global mean.
3. Contributions to twentieth-century global mean sea level change
In this section, we discuss the processes that influence GMSL change. We will explain how each contribution is evaluated, and compare the simulated contributions to observations where available.
a. Thermal expansion
Thermal expansion is one of the major contributors to twentieth-century GMSL (Church et al. 2013a) and is the only contribution that is simulated directly in the CMIP5 models. Averaged over the global oceans the density change due to salinity variations largely cancels out, and steric GMSL change is driven by temperature variations alone (Lowe and Gregory 2006).
We use two different ways to estimate model-based thermal expansion. The first uses the CMIP5 variable zostoga as it is provided by the CMIP5 modeling groups, which represents the thermal expansion over the entire model ocean basin and full depth of the ocean. However, most of the observational estimates of thermal expansion are only provided for the top 700 m of the ocean. Therefore, we also use a thermal expansion estimate derived from annual three-dimensional model fields of ocean potential temperature (to) in combination with a salinity (so) climatology, using the 1980 United Nations Educational, Scientific and Cultural Organization (UNESCO) international equation of state (IES80), following Melet and Meyssignac (2015). We use each models’ native grid for the computation, for two different depths: 0–700 m (TS700) and full depth (TSfull). Semienclosed basins (e.g., the Mediterranean Sea, Black Sea, Red Sea, Caspian Sea, Baltic Sea, Hudson Bay, and Great Lakes) are excluded from the TS computation, as ocean temperatures in some of the CMIP5 models are less reliable in these regions (Melet and Meyssignac 2015).
Five of the 12 climate models have no volcanic forcing included in their preindustrial control run (Table 1). As a result of the absence of volcanic forcing in the control run, the introduction of volcanic forcing in the historical simulations results in a cooling of the ocean, which leads to an underestimate of the thermal expansion for the twentieth century (Gregory 2010; Gregory et al. 2013b). This effect is visible in historical natural-only simulations (Slangen et al. 2015), which are not expected to show a long-term thermal expansion trend. To estimate the magnitude of the response to the introduction of volcanic forcing, we derive the global mean trend over 1850–2005 of the historical natural-only forced simulations for these models where possible (three models for zostoga and four for TS; Table 2), and otherwise use a constant value of 0.1 mm yr−1 as suggested in Gregory et al. (2013b) (two models for zostoga and one for TS). The derived values (Table 2) are then used to correct the historical simulations. This is a refinement in the methodology compared to Church et al. (2013b), where a correction of 0.1 mm yr−1 was added to all the model results.
The model results are compared to three observational reconstructions of thermosteric change. These reconstructions are all derived from in situ temperature measurements in the ocean, with the largest coverage and smallest uncertainties at the present day, and sparser data with larger uncertainties going back in time. For the upper 700 m of the ocean, we use annual time series from the three different reconstructions, which all use ocean temperature profiles as input, but the data are selected and processed in different ways to arrive at a global mean estimate:
Variations among the reconstructions can arise due to differences in the 1) input data and quality control procedures, 2) application and approach to correct interplatform biases, 3) choice of reference climatology, and 4) choice of mapping methods (Palmer et al. 2010; Abraham et al. 2013; Boyer et al. 2016).
To estimate the full depth thermal expansion, we use the 0–2000-m time series of NODC to compute the 700–2000-m contribution (1957–2013) and add a constant of 0.18 mm yr−1 for that same period, based on Purkey et al. (2014), to account for change in the deep ocean below 2000 m. All time series were downloaded on 9 June 2016. The common time period between the different observational time series is 1957–2012 (Table 3).
For the full depth of the ocean, the differences between TSfull and zostoga are small for most models (Fig. 2a). Note that TSfull (corrected) has an ensemble mean change from 1901–20 to 1996–2015 of 46 ± 33 mm (the uncertainty indicating the 1.65σ CMIP5 model spread), which is 3 mm larger than the zostoga ensemble mean, and each volcanically corrected ensemble is 4 mm larger than the respective uncorrected ensemble (Table 3). While the uncorrected CMIP5 ensemble mean is at the lower end of the observational range (Table 3), the corrected ensemble mean provides quite a good fit to the observations of Domingues et al. (2008) and NODC, although the difference is probably not significant given the uncertainty ranges of both models and observations. However, the spread in the models is larger than the observational uncertainties, which is consistent with studies where the observed ocean heat content was compared to climate model simulations (Cheng et al. 2016; Gleckler et al. 2016).
There are some individual models where the differences are larger between zostoga and TSfull, such as GFDL-CM3, MIROC5, MIROC-ESM, and GISS-E2-R. The main reason for this is probably the difference in land–ocean mask (i.e., inclusion/exclusion of marginal seas): some of these models have no connection between marginal seas and the open ocean, causing overestimated thermal expansion in the marginal seas (Melet and Meyssignac 2015).
Focusing on the upper 700 m of the ocean, the volcanic corrected TS700 model ensemble mean change (Fig. 2b, Table 3) agrees best with the total change in the NODC observational time series (21 ± 5 mm from 1957–61 to 2008–12), while the Ishii observational time series shows a smaller change (15 ± 5 mm) which is reproduced by some members of the climate model ensemble. When comparing observations to climate models, it is important to note that while historical external radiative forcing is used to drive the climate models, the models are not expected to match the historical internal climate variability (this includes, for instance, the timing of El Niños). This explains part of the discrepancies in the temporal variability of the observations compared to the models on annual to decadal time scales. However, there are also some consistent features. The response to volcanic forcing is present in the simulated thermal expansion (as this is part of the prescribed external forcing) and also in all three observation time series. For instance, there is a decrease between 1963 and 1970 that can be attributed to the Mount Agung volcanic eruption in 1963, and similar responses occur after the El Chichón and Mount Pinatubo eruptions (Fig. 2). The response between the models and volcanic eruptions varies because of different model responses to aerosol forcing (Slangen et al. 2015) and also because some volcanic forcing was longer sustained than others (e.g., Mount Agung; von Schuckmann et al. 2016). Note that, because of observational uncertainties, the statistical significance of the observed response is limited by the quality and coverage of the available pre-Argo observations (Cheng et al. 2017).
The mass loss of glaciers and ice caps around the world represents another major contribution to GMSL (Church et al. 2013a; Marzeion et al. 2017). The glacier contribution is simulated with a global glacier mass balance model (Marzeion et al. 2012, 2015). The model is driven by CMIP5 temperature (tas) and precipitation (pr), using the same CMIP5 models and simulations as for the thermal expansion contribution. The present-day glacier area used in the model is the Randolph Glacier Inventory, version 4.0 (Arendt et al. 2014), which is the state-of-the-art global glacier inventory. Mainly as a result of the updated glacier inventory and improvements in the digital elevation model (Marzeion et al. 2015; Slangen et al. 2017), these new glacier model results find a reduced glacier contribution compared to the results presented in Marzeion et al. (2012) and Church et al. (2013b). Antarctic peripheral glaciers are excluded from the model (see section 5). The glacier model first iteratively determines each glacier’s area and vertical extent in 1900, and is then forced by CMIP5 data to model the glacier evolution over the twentieth century and the resulting contributions to GMSL.
We use three different observation-based estimates of the global glacier mass contribution to GMSL. The first estimate is derived from glacier length records (Leclercq et al. 2011), which were updated with additional length records in Greenland and extended to 2010 as described in Marzeion et al. (2015). The second estimate is based on geodetic and direct mass change observations (Cogley 2009), using release 1301 as in Marzeion et al. (2015). Third, we use a modeled mass change estimate from the updated global glacier mass balance model of Marzeion et al. (2012, 2015) driven by gridded climate observations [Climatic Research Unit Time Series (CRU TS 3.22); Harris et al. 2014]. As this is based on observed temperature and precipitation data, the resulting modeled glacier mass balance changes follow the historical climate variability.
The CMIP5 ensemble mean contribution of the glaciers over the twentieth century (from 1901–20 to 1996–2015) is 55 ± 13 mm (Fig. 3; 1.65σ uncertainty based on CMIP5 model spread). The observational estimate by Leclercq et al. (2011) for 1900–2010 (86 ± 36 mm, gray shading) is larger than the modeled estimate for that same period (64 ± 17 mm for 1900–2010). The glacier-model estimate forced by temperature and precipitation from the CRU gridded data is only slightly larger for 1900–2010 (67 ± 13 mm; red solid line in Fig. 3) but larger than any of the individual models for the first half of the twentieth century. The differences are reduced in the second half of the twentieth century, when the Cogley (2009) observations are also available. The Cogley (2009) data suggest a contribution of 42 mm for 1950–2015 (blue shading, uncertainty estimated at 10% of the total value), compared to 42 ± 10 mm for the CMIP5 ensemble for that same period (uncertainty based on the CMIP5 ensemble spread).
In the beginning of the twentieth century, both the Leclercq et al. (2011) observations and the CRU-driven model results show a different behavior than the CMIP5-driven model results (Fig. 3). Despite larger uncertainties in earlier observational estimates, the combination of glacier observations and CRU-driven model results suggests that there is some common variability in the observations (both in the glacier length records from Leclercq and in the CRU data) that is not reproduced in the CMIP5 models. A closer look at the regions (Marzeion et al. 2012) shows that most of the difference in the global mean is caused by a discrepancy in mass change in the glaciated regions around Greenland, corresponding to a period of warming and strong glacier melt that is present in the CRU-forced model and also found in observations (Bjørk et al. 2012), but not replicated by the CMIP5 models. There may be several reasons for this. As was mentioned before, the climate models are not forced to follow the “real world” internal climate variability; they are only forced by changes in external radiative forcing (e.g., through changing greenhouse gas concentrations and aerosol concentrations). Another reason may be a bias that was found in the atmospheric circulation of the climate models, which do not reproduce warm airflow over the south of Greenland, resulting in too little melt in this region (Fettweis et al. 2013). Third, the spatial patterns of the historical (natural and anthropogenic) aerosol forcing that drive the CMIP5 models are more uncertain going back in time (which may cause discrepancies in regional response of the model), and the (regional) response of climate models to aerosols in general is relatively uncertain compared to the response to greenhouse gas forcing. We therefore propose a bias correction to account for the reduced response in the climate models. We follow the approach of Slangen et al. (2016), and compute the ensemble mean difference between the CMIP5 models and the observation-based estimates [averaging the available twentieth-century Leclercq et al. (2011) and CRU-based observational estimates] to account for additional mass change in glaciers around Greenland in the early twentieth century. The proposed correction amounts to 15 ± 10 mm over the period from 1901–20 to 1996–2015 (Fig. 3, green dashed line; 1.65σ uncertainty from CMIP5 model spread). This correction is consistent with a 20-mm correction proposed in Church et al. (2013b) to account for internal variability in Greenland temperatures, and both corrections have the largest rates between 1930 and 1960. As the difference is very small after 1960, the correction is set to zero for the period 1960–2015.
c. Ice sheets
1) Surface mass balance
(i) Greenland surface mass balance
The modeled Greenland surface mass balance (SMB) contribution to twentieth-century GMSL is estimated using a regional statistical downscaling technique, which accounts for the nonuniform distribution of SMB change over Greenland (M17). Quadratic relations between the Greenland SMB change and CMIP5 variables (annual snowfall, CMIP5 variable prsn, and atmospheric summer temperature at 600 hPa, CMIP5 variable ta) were derived using CMIP5-forced MARv3.5 regional climate model simulations (Fettweis et al. 2013). MARv3.5 has only been forced with three CMIP5 models (MIROC5, NorESM1-M, and CanESM2) over 1900–2006 due to computational limitations. The inferred statistical downscaling is then applied to the other climate model output. This downscaling technique uses six major drainage basins of the Greenland ice sheet.
There are no direct observational time series of Greenland SMB in the twentieth century, but an observational estimate can be obtained using reanalysis data (instead of CMIP5 data) and the downscaling relation between temperature, precipitation, and SMB (Fettweis et al. 2013). Three reanalysis products are used to provide the observation-based estimate: ERA-20C (Poli et al. 2016), 20CRv2, and 20CRv2c (Compo et al. 2011).
The CMIP5-derived Greenland SMB contribution is at least an order of magnitude smaller than the contributions from glaciers and thermal expansion (Fig. 4). The ensemble mean contribution amounts to 0 ± 3 mm (1.65σ CMIP5 model spread) over the period 1901–20 to 1996–2015. For the period since 1990, both the CMIP5 results and the reanalyses show a sharp increase in the Greenland SMB contribution to GMSL, which is confirmed by observations and regional climate modeling (van Angelen et al. 2014). In the beginning of the twentieth century, the three reanalysis-based Greenland SMB estimates lead to a significantly larger contribution to GMSL than the CMIP5 models (Fig. 4), similar to the glacier contribution (section 3b). This is thought to be caused by an increase in air temperatures in and around Greenland around 1900–40, which led to increased melt in Greenland (Box et al. 2009; Bjørk et al. 2012; Fettweis et al. 2017) and surrounding glaciers (section 3b; Marzeion et al. 2012). This temperature increase is found in the reanalyses but not in the CMIP5 models. Similar to the approach for the glaciers, we propose a bias correction to account for this additional internal climate variability using the difference between the mean of the observational estimates and the CMIP5 ensemble mean (Fig. 4, dotted green line), which is 10 ± 4 mm from 1901–20 to 1996–2015 (1.65σ CMIP5 model spread) and set to zero from 1970 onward. Although this correction seems substantial, the corrected Greenland SMB contribution is a low-end estimate compared to independent observational evidence (Kjeldsen et al. 2015), suggesting a contribution of the whole of the Greenland ice sheet of 25 ± 9 mm (1σ) over the twentieth century. Note that the Kjeldsen estimate includes both SMB and ice dynamical changes but excludes glaciers peripheral to the Greenland ice sheet.
(ii) Antarctic surface mass balance
The contribution of Antarctic SMB changes to twentieth-century GMSL is estimated with two different methods, both using CMIP5 data. The first method uses the difference between precipitation (CMIP5 variable pr) and evaporation (CMIP5 variable evspsbl) on the Antarctic ice sheet [precipitation minus evaporation (PME)]. Although CMIP5 models usually have no correct representation of snowpack processes, the freshwater runoff on the Antarctic ice sheet is negligible for the twentieth century (Lenaerts et al. 2012), and thus PME is sufficient to represent the Antarctic SMB. The resulting SMB time series are scaled to fit the best estimate of the Antarctic SMB for the period 1985–2010 (merging historical and RCP8.5 experiments to go beyond 2005) from the regional climate model RACMO2.1 forced by ERA-Interim data (Lenaerts et al. 2012).
The second method to estimate Antarctic SMB changes uses CMIP5 surface temperature changes to estimate the change in precipitation over the Antarctic ice sheet (henceforth TAS), again with respect to the best estimate SMB from Lenaerts et al. (2012). As warmer air has a higher moisture holding capacity, the snowfall over Antarctica increases when temperatures increase, thus impacting the SMB of the ice sheet. The change in SMB is estimated at approximately 6% ± 0.7% K−1, which is consistent with observations (ice cores) and different global and regional climate models (Ligtenberg et al. 2013; Frieler et al. 2015).
The modeled SMB ensemble mean estimates from both methods show a small sea level fall over the twentieth century of −6 ± 3 mm for the TAS method and −6 ± 7 mm for the PME method (Fig. 5; 1.65σ model spread). The spread in the PME estimates is larger than in the TAS estimates (i.e., solid lines in Fig. 5 compared to dashed lines), probably because temperature is a more robust variable in the models than precipitation and evaporation. However, in the ensemble mean there is no significant difference for the 1900–2015 period. Both methods initially show only a minor Antarctic SMB contribution to GMSL (only 1 ± 3 mm up to 1970) but indicate an accelerated sea level fall from 1970 onward (−7 ± 5 mm for TAS and −8 ± 8 mm for PME between 1971 and 2015).
However, we note that although the relation between increasing temperature and increasing precipitation is robust, this small sea level fall is not yet significant when considering the large internal (snowfall) variability in Antarctic SMB change (on the order of 0.3–0.4 mm yr−1). A recent publication estimated that the signal is only projected to emerge from the noise in the first half of the twenty-first century (Previdi et al. 2016). This is in line with regional climate model results, which found no significant trend in Antarctic SMB for the period 1979–2010 (Lenaerts et al. 2012).
2) Ice sheet dynamics
Apart from changes in the surface mass balance, the Greenland and Antarctic ice sheets can also contribute to GMSL due to changes in dynamical discharge processes. Processes that could contribute to dynamical discharge include, for instance, the melt of floating ice shelves due to ocean warming or calving of marine-terminating glaciers (Church et al. 2013a, and references therein). However, these processes are not included in CMIP5 models and therefore have to be estimated in another way.
Here we use the time series from Shepherd et al. (2012), a community effort to present a complete ice sheet mass balance estimate that reconciled different methods and models. The total ice sheet mass balance estimate is available from 1993 to 2010 in four regions: West Antarctica, East Antarctica, the Antarctic Peninsula, and Greenland. The ice-dynamical estimate is computed by removing the SMB (based on RACMO2.1; Lenaerts et al. 2012; van Angelen et al. 2014) from the total ice sheet mass balance estimate. The ice dynamical time series are extended to 2015 using discharge estimates from Sutterley et al. (2014) and Enderlin et al. (2014), which assume that the West Antarctic discharge was slightly above the 2008–12 average and that the Greenland ice sheet discharge remained constant at the 2010 value. Discharge for East Antarctica and the Antarctic Peninsula is assumed to follow the 2001–10 average. The extrapolated values are in line with recent ice sheet mass balance studies (Velicogna et al. 2014; Harig and Simons 2015; McMillan et al. 2016).
The resulting total ice-dynamical contribution to GMSL over the period from 1993–95 to 2011–15 is 7 ± 4 mm from the Greenland ice sheet and 7 ± 8 mm from the Antarctic ice sheet (Fig. 6, Table 4). The 1σ uncertainties on the ice dynamics are derived by assuming they have the same uncertainty-to-signal ratios as the total mass balance change presented in Table 1 of Shepherd et al. (2012), which is 35% for Greenland and 75% for Antarctica.
Several studies have suggested that in the twentieth century there was a nonzero contribution to GMSL from long-term nonequilibrium changes on the ice sheets and/or in the deep ocean as a result of their relatively long response times to climate variations, which is not represented in the models (e.g., Gregory et al. 2013b; Church et al. 2011, 2013b). In Gregory et al. (2013a), this long-term change was suggested to range between 0 and 0.4 mm yr−1, the majority from Antarctica, which is supported by Masson-Delmotte et al. (2013), who find geological evidence for an ongoing long-term GMSL contribution of about 0.2 mm yr−1 in the last two millennia. This long-term change is not represented in the modeled SMB changes, which assumes a much faster response to climate changes. It is also important to keep in mind that in the climate models the dedrifting procedure removes long-term trends, which imposes an artificial equilibrium on the background trend. This will particularly affect the deep ocean, which takes much longer to respond to changes in forcing than the upper ocean. This is consistent with the thermal expansion estimates presented in Table 3, where the full-depth thermal expansion in the observations is underestimated by the modeled estimates, in contrast to the 0–700-m thermal expansion. We therefore propose to use a nonequilibrium constant to account for a combined ongoing ice sheet/deep ocean contribution in the twentieth century (Fig. 6). A value of 0.13 ± 0.35 mm yr−1 was derived in Slangen et al. (2016) using CMIP5 thermal expansion (CMIP5 variable zostoga) compared to observed thermal expansion (the same estimates as used in this study) for the period 1970–2005 (see their supplementary Fig. S2 for the derivation).
d. Terrestrial water storage change
Water is continuously being transferred between the ocean and land as part of the natural hydrological cycle and partly due to human intervention. While natural terrestrial water exchange, such as snowfall or surface water storage, is assumed to be in equilibrium on longer (from decadal to century) time scales (Church et al. 2013a; Wada et al. 2017), processes such as groundwater extraction and reservoir building can have a long-lasting effect on GMSL. People extract groundwater from deep reservoirs, to use for instance as drinking water, for industries, and to irrigate the land, bringing additional water into the active hydrological cycle and into the ocean. On the other hand, dams and reservoirs are built in rivers, leading to additional storage of water on land and a reduction of water flowing into the ocean.
However, these processes are not included in the climate models. Therefore, we use observation-based estimates of the twentieth-century contributions to GMSL from groundwater extraction and reservoir building. For water storage behind dams, we use the Chao et al. (2008) estimate (1900–2007), scaled by 0.85 of the nominal storage capacity as suggested in IPCC AR5 (Church et al. 2013b). We use two estimates for groundwater depletion. The first comes from Wada et al. (2012) for 1900–2005, corrected for recent findings in Wada et al. (2016), showing that only 80% of the depleted groundwater ends up in the ocean. This brings their estimate closer to earlier estimates that also reflected smaller contributions to GMSL from groundwater extraction (e.g., Konikow 2011). Wada et al. (2012, 2016) used a global hydrological model to compute groundwater depletion using country-specific groundwater abstraction data, but they did not model the reduction in base flow due to groundwater abstraction. The second estimate, from Döll et al. (2014), is available from 1902 to 2009 and presents a smaller contribution than that of Wada et al. (2012, 2016). They used a different hydrological model, which computes the change in groundwater storage from the difference between groundwater abstraction and recharge, including a decline in base flow as the storage decreases. Both groundwater depletion time series, as well as the reservoir storage time series, were extended to 2015 using the average rate in the last 5 years of available data.
The rates of change for reservoir storage and groundwater extraction are small before 1950 (Fig. 7). The reservoir storage then sharply increases, leading to sea level fall, and starts tapering off toward the end of the century because dam construction declined. Both groundwater estimates show a more gradual acceleration after 1950, with most of the increase toward the end of the period, contributing to sea level rise throughout the century. The reservoir storage (−25 mm from 1901–20 to 1996–2015) is larger than the groundwater contribution of 16 mm for the same period from Wada et al. (2012, 2016) and 11 mm for Döll et al. (2014). This study uses the average of the studies of Wada and Döll as the groundwater estimate, and the difference between the two groundwater studies (5 mm) is used as the 1σ uncertainty estimate for the landwater storage change component. Combined with the reservoir storage, this leads to an overall GMSL change of −11 ± 7 mm (1.65σ) from 1901–20 to 1996–2015.
4. Modeled sea level change compared to observations and reconstructions
a. The twentieth century to present day (1900–2015)
For the whole of the twentieth century, glacier mass loss and ocean thermal expansion are the largest contributors to GMSL (Fig. 8, Table 4). The proposed correction to the early-twentieth-century glacier contribution (section 3b) makes the glaciers the largest contributor in the beginning of the century, whereas from the middle of the century onward the thermal expansion accelerates faster than the glaciers. Both contributions show a change in pace around the 1990s (Fig. 8). From 1900 to 1990 the thermal expansion contributes approximately 30 mm in 90 years, which is the same amount as in the 25 years after 1990. This agrees with observational estimates of the increase in ocean heat content in the twentieth century (Gleckler et al. 2016). The change in pace is a result of the recovery after the Pinatubo eruption, as well as increasing greenhouse gas concentrations in tandem with decreasing anthropogenic aerosol concentrations (Gleckler et al. 2006; Slangen et al. 2014b). The glaciers contributed 50 mm (65 mm including the proposed corrections) in the first 90 years and about 25 mm in the last 25 years, doubling their rate of change. The landwater contribution causes a sea level fall before the 1950s due to the increasing reservoir storage, while closer to the present day the groundwater extraction increases, leading to a small and increasing positive contribution to GMSL from landwater storage changes. The ice sheet SMB contributions are relatively small compared to all the other contributions, apart from the proposed bias correction to the Greenland SMB component, which is of similar magnitude (but opposite sign) as the landwater storage change. Finally, the ice sheet dynamical component is initially dominated by the proposed nonequilibrium term, but by the end of the century the ice sheets start to show an increasing contribution from ice sheet dynamical discharge.
When all the contributions are combined, the models add up to a GMSL change of 92 ± 47 mm for the period from 1901–20 to 1996–2015 (Table 4, Fig. 9a). Compared to the average of the four reconstructed global mean time series for the overlapping period from 1901–20 to 1988–2007 (Table 5, Fig. 9a, the model simulations clearly underestimate the observed GMSL and explain only 50% ± 30% of the observed change (using ±1.65σ of the models to the mean of the observations). For the individual reconstructions, the explained percentages range from 43% (Jevrejeva et al. 2014) to 61% (Hay et al. 2015). These low percentages are mainly due to differences in the earlier part of the century, whereas the percentages significantly increase for later periods (see section 4b).
If the proposed corrections are added to the model simulations, the simulations explain a larger percentage of 75% ± 38% of the averaged reconstructed changes (Fig. 9b), with the lowest explained percentage being 65% (Jevrejeva et al. 2014) and the highest 92% (Hay et al. 2015). The increased percentages after adding the corrections indicate that the component biases (between models and observations of individual contributions) are consistent with the total sea level bias and that improvements are needed in the simulation of mass terms in the early twentieth century to be able to close the sea level budget. With the corrections, the twentieth-century simulated budget falls within 1.65σ uncertainty of the tide gauge observations (Table 5) and nearly closes the budget for the Hay et al. (2015) reconstruction.
Both the models and the observations show a more or less linear increase before 1950, followed by a short plateau, before starting to accelerate toward the present day (Fig. 9). However, the year-to-year variability in the observations is much larger than in the simulations, which show little variability apart from the response to volcanic eruptions. This may be explained by the fact that observations contain much more small-scale “noise.” The relatively large variability in the sparsely distributed observations tends to lead to an overestimate of the interannual variability in the global mean, as shown by comparisons between tide gauge reconstructions and satellite altimetry (Church and White 2011; Meyssignac et al. 2012). On the other hand, the variability in the climate model global mean may be underestimated because of processes that are either underestimated or absent from the climate models. Reasons for this could be the relatively coarse (ocean and atmospheric) grids of the models (both in latitude/longitude and in depth), which require parameterizations of subgrid-scale processes, such as eddies (e.g., Hallberg 2013) or turbulence (e.g., Melet et al. 2016).
b. The satellite era (1993–2015)
In the satellite era (1993–2015) the thermal expansion is the largest contribution, followed by the glaciers and the ice dynamic contribution (Fig. 10).This is different compared to the whole of the century (Fig. 8), where the ice sheet dynamics without the nonequilibrium correction were small: this contribution only starts to contribute from the 1990s onward. The landwater contribution in the satellite era is small, as the reservoir storage and the groundwater extraction largely cancel each other out. The Antarctic SMB contribution on its own would cause a sea level fall due to the increased snowfall on the ice sheet as a result of warmer air temperatures. Most contributions, apart from the Antarctic SMB, show a more or less linear increase for 1993–2015. As we saw in the previous section, both the thermal expansion and glacier contributions have increased significantly in these last 25 years compared to earlier in the twentieth century, doubling to tripling the contributions in recent times. The thermal expansion contribution shows a temporary response to the 1991 Pinatubo eruption for the two consecutive years but resumes a linear increase from 1993 onward.
Combined, the simulated GMSL change from 1993–97 to 2011–15 is 56 ± 18 mm (Table 6). The observations (tide gauges and satellites) and model estimates are much closer to each other over this more recent period than they were for the full twentieth century (cf. Figs. 11 and 9), agreeing well with a recent review study on the sea level budget during the altimetry era by Chambers et al. (2017). Also, the proposed corrections have a much smaller influence for this period than over the twentieth century, as the corrections due to early-twentieth-century warming around Greenland have no impact after 1970, and only the nonequilibrium contribution is still included. This leads to a modeled change of 58 ± 19 mm over the altimetry period (from 1993–97 to 2011–15), which is only 2 mm larger than the uncorrected model estimate. Since this is a relatively short period, and the averages are only taken over 5 years, one should keep in mind that internal climate variability may influence these values. The two satellite altimetry time series from ESA-CCI and CSIRO find nearly the same change from 1993–97 to 2011–15 (Table 6).
We also compare the models to the recently published satellite time series of Watson et al. (2015), who proposed to calibrate the satellite record with GPS-corrected tide gauges in an attempt to remove a possible bias drift in the first part of the altimeter record. This is in contrast to the treatment of the other two satellite time series above, which are completely independent from tide gauge estimates. The Watson et al. (2015) satellite time series has a smaller GMSL change than the other two estimates and less than the model ensemble means (Table 6) but is still within the model uncertainty range.
Over the satellite period, the model simulations explain 102% ± 33% of the satellite observations [Table 6; mean of the three satellite estimates, 95% for CSIRO, 95% for ESA-CCI, and 119% for Watson et al. (2015)]. When proposed corrections are included the percentages increase to 105% ± 35% [98% for CSIRO, 98% for ESA-CCI, and 123% for Watson et al. (2015)]. The percentages for both the uncorrected and the corrected simulations of the satellite period show a significant improvement compared to the whole of the twentieth century. Although there is a much smaller impact of the proposed corrections for this more recent period, they still improve the budget closure.
c. Increases in trends and changes in relative contributions to twentieth-century sea level rise
The yearly trends, averaged over 10 years (Fig. 12a), in simulated GMSL over the twentieth century are increasing from 0.7 mm yr−1 for 10 decades centered on 1904–13 to 3.2 mm yr−1 for 10 decades centered on 2001–10, and increase from 1.4 to 3.4 mm yr−1 for the corrected simulations. This is due to increasing trends in both the thermal expansion (from 0.2 to 1.4 mm yr−1) and the combined mass component of glaciers, ice sheets, and landwater (from 0.5 to 1.8 mm yr−1 for uncorrected and from 1.2 to 2.0 mm yr−1 corrected simulations). The increase in thermal expansion is consistent with Earth’s energy imbalance of the last decades, as a result of the long-term increase in radiative forcing associated with increased greenhouse gas emissions (Smith et al. 2015; von Schuckmann et al. 2016).
The 10-yr averaged trends in the observations also increase in the twentieth century, from 1.7 mm yr−1 (tide gauge trends 10-yr running means centered on 1904–13) to 2.9 mm yr−1 (satellite 10-yr running mean trends centered on 2001–10). The observations show 1 than the simulations (Fig. 12), but they fall within 1.65 standard deviations (1.65σ) of the model spread for most decades. The period of 1930–50 shows similar observed trends to the period 1990–2010, with the difference that in the more recent period the trends continue to increase rather than decrease as they did in 1950–60 (Fig. 12a). The difference is probably that in the beginning of the century the trends are driven by internal variability, whereas later in the century they are driven by external forcing (Slangen et al. 2016). However, part of the explanation could also be in the observations, which may have too much variability in the observations due to sparse spatial sampling (Church and White 2011; Meyssignac et al. 2012).
In the beginning of the twentieth century, the majority of the (simulated) GMSL comes from mass-driven contributions (i.e., glaciers, ice sheets, and landwater) rather than thermal expansion. This is enhanced when the mass bias corrections are included (Figs. 12b,c). These early-twentieth-century changes have been attributed to natural forcing and internal climate variability (e.g., the glaciers’ response to the end of the Little Ice Age and early-twentieth-century warming around Greenland) rather than anthropogenic forcing (Marzeion et al. 2014; Slangen et al. 2016). In contrast, for the late twentieth century (1991–2010), 69% ± 40% (±1.65σ) of the glacier contribution has been attributed to anthropogenic change (Marzeion et al. 2014), which clearly shows a shift from natural to anthropogenic forcing over the course of the twentieth century. Thermal expansion responds differently to external forcing: it starts to contribute positively to GMSL from 1910 onward and by 2015 is responsible for 46% of the cumulative simulated GMSL change since 1900 (34% of the corrected simulations). This contribution is mainly attributable to anthropogenic forcing (>80% since 1970; Marcos and Amores 2014; Slangen et al. 2014a). When the thermal expansion and mass contributions are combined, the majority of GMSL change since 1970 (69% ± 26%; ±1.65σ) is attributable to anthropogenic forcing (Slangen et al. 2016).
The mass term contributes positively to GMSL throughout the twentieth century (even though some individual mass contributions such as landwater storage can be negative contributors), but its relative contribution decreases as thermal expansion increases (Fig. 12c). Initially, the mass contribution almost entirely consists of the glacier contribution (and Greenland SMB when the bias corrections are included), but the relative contribution from glaciers decreases from 1993 onward. The relative ice dynamical contribution (included in the mass contribution) starts from 1993, accounts for 12% of the cumulative change by 2015 (18% when corrections are included), and is rapidly increasing over the period 1993–2015. Although both the mass terms and the thermal expansion are increasingly contributing to GMSL in an absolute sense (Figs. 12a,b), their relative contributions have stabilized from 2000 onward (Fig. 12c), indicating that the increase in the thermal expansion and the increase in the mass contributions are following the same pace for now. As the future changes in the ice sheets are still quite uncertain, in particularly the contribution of the Antarctic ice sheet (Church et al. 2013a; Slangen et al. 2014a; Kopp et al. 2014; Ritz et al. 2015; DeConto and Pollard 2016), the relative contributions may change if the mass contribution overtakes the thermal expansion contribution.
5. Discussion and conclusions
In this paper, we focused on evaluating twentieth-century observed GMSL change using model estimates. To get a consistent dataset, we have used the same 12 climate models’ data output to estimate the contributions of thermal expansion (assuming salinity change is negligible in the global mean), glaciers and ice sheet SMB, to GMSL. The simulated GMSL is completed by adding observation-based estimates of ice sheet dynamical processes and landwater storage changes. We use four sets of tide gauge reconstructions and three satellite altimetry time series to compare our model estimates to. With respect to previous comparisons (Church et al. 2013b; Slangen et al. 2016), we have extended the evaluation period up to 2015 using the RCP8.5 climate model scenario, and we use a consistent climate model ensemble, which allows us to evaluate individual climate model results.
Other differences with respect to IPCC AR5 and previous model–observation comparisons (Church et al. 2013a,b; Slangen et al. 2016) are in the individual contributions, where recent findings and model updates change some of the estimates, and different approaches were used to compute the ice sheet contributions. This leads to more consistent and robust estimates of the contributions to GMSL. The glacier contribution, for instance, contains an updated glacier inventory and improvements to the digital elevation model, which have caused a decrease in the CMIP5 estimated twentieth-century glacier contribution compared to the IPCC AR5 results (Marzeion et al. 2015).
The contributions of the Greenland and Antarctic ice sheets are a combination of increases in the runoff of surface meltwater, enhanced calving of solid ice, and melting of ice that is in contact with the ocean, with the former process dominating in Greenland and the latter two in Antarctica. Here, we use different methods to estimate these contributions, using CMIP5 data where possible to obtain a consistent dataset, an improvement with respect to previous work. At present, there are no models that we are aware of that simulate a robust and complete ice sheet in response to twentieth-century climate, but this is expected to change with the forthcoming CMIP6-endorsed Ice Sheet Model Intercomparison (ISMIP; Nowicki et al. 2016). As the future contribution of these ice sheets is expected to increase in a warming climate, such comprehensive simulations of the ice sheet contribution will be a significant contribution to the understanding of historical GMSL and evaluating our current understanding of GMSL and ice sheet models. Progress in these areas will lead to more robust projections of sea level change.
In addition, recent groundwater extraction estimates were used (e.g., Döll et al. 2014), which tend to be lower than the values included in IPCC AR5. This is in agreement with another recent publication, showing that only 80% of the extracted groundwater of Wada et al. (2012) ultimately reaches the ocean (Wada et al. 2016).
We furthermore suggest bias corrections for some of the GMSL contributions, based on differences between models and observations. For instance, observational evidence points to a much larger (SMB + ice dynamical) contribution from the Greenland ice sheet than previously thought (Kjeldsen et al. 2015). This is in line with the finding that the glacier contribution and the Greenland SMB contribution are much larger in the early twentieth century when they are computed using temperature reanalyses rather than CMIP5 model estimates (Slangen et al. 2016). In this study, we include a correction for the Greenland SMB contribution based on the differences between the CMIP5-driven results and the reanalysis driven results.
Following Slangen et al. (2016), we also explored the possibility that ice sheets and the deep ocean are not in equilibrium with twentieth-century climate by introducing a nonequilibrium constant, as their response time is likely to be on a century to millennia time scale. We use a constant of 0.13 ± 0.35 mm yr−1 as derived in Slangen et al. (2016).
Throughout this study, all the modeled and observational glacier estimates discussed have excluded the contribution from the Antarctic peripheral glaciers (PGs), following the approach taken in IPCC AR5, as this is an extremely uncertain contribution. Observations of these glaciers are sparse and the mechanics are poorly understood: it is unclear whether they should be included in the glacier estimate or the Antarctic SMB estimate. For a rough estimate of the PG contribution, Marzeion et al. (2012) suggested to scale the total area of the PGs to the total global glacier area (minus PGs) to derive a scaling factor that can be used to scale the global glacier contribution with. This scaling factor of 0.22 (132 867/ 593 925 m2) would result in an additional sea level contribution from the PGs of 12 mm from 1901–20 to 1996–2015, which would reduce the twentieth-century difference between observations and simulations. If the PG were assumed to follow the Antarctic SMB evolution the contribution would be much smaller, as the Antarctic contribution is small throughout the twentieth century. However, as the contribution is so poorly constrained, we have decided not to include this contribution in our simulated GMSL estimates.
When all the contributions are combined, there is still a large gap between the observations and the models, and only 50% ± 30% of the observations (mean of four tide gauge reconstructions) can be explained by the models for the period from 1901–20 to 1988–2007. The suggested bias corrections for Greenland SMB, glaciers, and deep ocean/ice sheet reduce the model–observation gap by construction, as they are based on model–observation differences, bringing the explained percentage to 75% ± 38% for the mean of the four reconstructions. Compared to the individual reconstructions, the bias-corrected simulations agree best with the Hay et al. (2015) reconstruction, explaining 92% of the observed change (Table 5).
For the more recent satellite period (from 1993 to 2015), the explained percentage (i.e., explained by the simulations) is 102% ± 33% (105% ± 35% when bias corrections are included), effectively closing the sea level budget for this period. In this later period, the uncertainties in the observations are smaller as the data resolution is higher, both spatially and temporally. The model uncertainties, mainly based on the spread between the different CMIP5 models (Table 4), are smaller in absolute sense for the later period, but the percentage remains the same, showing that the uncertainties in the models are relatively constant for the different periods.
The simulated GMSL time series shows increasing trends over the twentieth century due to increasing contributions from both thermal expansion and the mass components. Thermal expansion starts to contribute to GMSL from 1910 onward, and by 2015 accounts for 46% of the total simulated sea level rise. The mass contribution, which accounts for the remaining 54% in 2015, is dominated by the glacier contribution until the ice sheet dynamics start to play a role at the end of the twentieth century, accounting for 12% of total simulated GMSL in 2015.
In a warming climate, all of the contributions are expected to continue to increase in the coming century, with potentially the largest increase in the contributions from the ice sheets (Church et al. 2013a). Because of the delayed response in the deep ocean and ice sheets, sea level will continue to rise even if the emission of greenhouse gases were reduced today. It is therefore important to understand the response times and the magnitude of the different processes to better estimate future changes in sea level. This paper, presenting the GMSL budget, is a step forward in this process. The next step is to evaluate and understand regional patterns in sea level change as a result of these processes, which is done in the companion paper (M17).
This paper is an outcome of work with the ISSI international team on “Contemporary regional and global sea-level rise: assessment of satellite and in-situ observations and climate models.” We thank ISSI (www.issibern.ch) for facilitating the team meetings. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and leads development of software infrastructure in partnership with the Global Organization for Earth System Science Portal. We thank P. Leclercq (U. Oslo) for providing the updated glacier contribution based on glacier length records. A.S. received funding from a CSIRO Office of the Chief Executive postdoctoral fellowship and the NWO-Netherlands Polar Programme. M.P. and C.R. were supported by the Joint UK BEIS/Defra Met Office Hadley Centre Climate Programme (GA01101). In the early part of this study, J.C. was funded by CSIRO, the Australian Climate Change Science Program, and the University of New South Wales. S.L. was supported by an NWO ALW Veni Grant (865.15.023). K.R. was supported by the Austrian Science Fund (FP253620).