1 Introduction
Satellite-based estimates of the mass balance of the Antarctic ice sheet (AIS) allow one to quantify the response of the AIS to global warming. Projections show that solely the AIS may contribute between 4 and 34 cm to global mean sea level by 2100 . However, even recent mass loss estimates of the AIS vary over a wide range, e.g. from 94 to 202 Gt a over the time period 2010–2019 . This large spread reveals a lack of knowledge which propagates to projections.
The mass balance of an ice sheet – also referred to as the ice mass change (IMC) of an ice sheet – is the difference of the input mass flux, i.e. mainly the accumulation by precipitation, and the output flux, i.e. for the largest part ice discharge and meltwater runoff into the ocean. Commonly, three methods are applied to determine an ice sheet's mass balance using satellite data: (i) the gravimetric method deriving the mass balance from gravitational field changes measured by the Gravity Recovery and Climate Experiment (GRACE) and Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) satellite missions, (ii) the altimetric method deriving the mass balance from surface elevation changes measured by several radar and laser altimeter missions while assuming a volume-to-mass conversion, and (iii) the mass budget method deriving the mass balance by assessing the difference between the input and output mass fluxes. They are derived from regional climate modelling and from ice discharge estimates, which can be retrieved from remote-sensing satellite data and ice thickness data. All three methods have advantages but also limitations, which contribute to the large spread of estimates mentioned above, extensively documented elsewhere (e.g. ). To summarize, method (i) has the advantage of being directly sensitive towards mass changes, but there is the need to exclude all other sources of mass redistributions, which are superimposed onto the gravitational field changes. The present-day effect due to glacial isostatic adjustment (GIA) is the most relevant and most uncertain (e.g. ). The term present-day GIA effect refers to the presently observable effects resulting from the adjustment process to an isostatic state, which was induced by glacial mass changes in the past. This is to be distinguished from effects associated with the instantaneous elastic response to contemporaneous ice mass loading changes . GIA predictions for Antarctica differ by several tens of gigatonnes per year and disagree in their spatial patterns due to assumptions on rheology and ice loading history . Furthermore, IMC estimates derived from gravitational field changes usually only allow for a spatial resolution of a few hundred kilometres. Method (ii) has the advantage of capturing IMC with high spatial resolution (e.g. ), but the conversion from volume changes to mass changes is based on effective density hypotheses or needs to include auxiliary data, e.g. firn modelling results, where the uncertainties are largely unknown. Method (iii) has the advantage of aiming to resolve the full mass fluxes and not only the differential signal between input and output fluxes . However, as the mass balance amounts to only less than 10 % of the magnitude of the mass fluxes, even small errors have a strong impact on the result. The input flux and the output flux are each subject to an uncertainty that is on the order of magnitude of the AIS mass balance itself .
In addition to these mass balance estimation strategies, there are methods that combine data from satellite gravimetry and satellite altimetry to build on the advantages of both sensors. Some of the combination approaches aim to co-estimate the GIA effect rather than using GIA modelling results to account for the GIA signal . These regional approaches implement regional constraints to evaluate the global gravitational fields from satellite gravimetry in a regional domain. Except for , these approaches only allow a smoothed estimation of IMC with a spatial resolution comparable to GRACE and GRACE-FO-only estimates. In addition to these regional approaches, there are global inversion approaches that co-estimate the GIA signal by fitting prescribed GIA spatial patterns to data or by utilizing the signal covariance information. These approaches are limited in the sense that the estimated GIA depends on the applied modelling output. Hence, they fit presumably erroneous a priori information to the input data. The approaches according to and only allow for a basin-wise and a smoothed estimation of AIS IMC, respectively.
Here we extend the work of by applying a global inversion framework with a focus on Antarctica. We use data sets from satellite gravimetry, satellite altimetry, and regional climate and firn modelling. The latter are used to derive changes of the firn air content (FAC) based on the surface mass balance (SMB) from a regional climate model and firn thickness changes of a firn densification model (FDM). This approach aims to overcome some limitations of previous combination approaches that allow us to jointly determine IMC and GIA. First, it is a global framework, i.e. no regional constraints need to be implemented, and mass changes are parametrized across the globe. Second, the approach applies a GIA parametrization in Antarctica, utilizing local deglaciation impulse response patterns, which are globally consistent. In principle, this GIA parametrization allows us to spatially resolve GIA effects in Antarctica which have not been predicted by GIA forward modelling. Third, the IMC parametrization facilitates a spatial resolution of 50 km, which is useful to explore AIS IMC in more detail than results from previous combination studies allow. Fourth, the approach includes a parametrization for FAC changes, which allows us to circumvent the implementation of a firn density. Finally, the approach enables us to incorporate the error covariance information of all input data sets for rigorous accounting of input data quality limitations. The feasibility of the approach was demonstrated with simulation experiments .
We present and analyse results from applying this approach over the 10-year observation period from January 2011 to December 2020 (2011–2021), during which the following data sets are available at the same time: a satellite gravimetry data product from GRACE and GRACE-FO (ITSG-Grace2018; ), a satellite altimetry data product from CryoSat-2 , and changes of FAC derived from the RACMO2.3p2 SMB and IMAU-FDM v1.2A . We validate the results with independent global navigation satellite system (GNSS) data.
2 Material and methods
Our aim is to disentangle the IMC of the AIS, which is superimposed with other signals in global gravitational field changes observed by GRACE and GRACE-FO . From the perspective of Antarctica, gravitational field changes are caused by the global mass redistribution due to AIS IMC, GIA, and far-field effects from other mass redistributions in the Earth system. We refer the reader to recent review articles (e.g. ) for comprehensive background information on the physical processes related to mass changes of the AIS.
2.1 Global inversion framework
We apply an updated global inversion approach from , which builds upon the work from . The Supplement provides the updates we have made to the methodology described in in Sect. S1. In the following, we describe the applied inversion methodology. For a more extensive theoretical background and further details about the inversion setup, we refer the reader to and .
The global fingerprint inversion from enables one to partition observed sea level and to quantify the individual sea level budget components. For this purpose, globally consistent spatial patterns of the individual budget components are derived from a priori information. These spatial patterns serve as fingerprints in the inversion. Scaling factors for the individual fingerprints are then computed via a parameter estimation, utilizing observations from satellite altimetry over the ocean and from satellite gravimetry. The quality of the a priori information crucially affects the final result. found that the scaling factor of the Antarctic GIA fingerprint in particular was estimated too low, meaning that the GIA effect determined in Antarctica is likely unrealistic.
The inversion approach presented here is designed to resolve this aspect by co-estimating globally consistent AIS IMC and GIA from gravitational field changes. However, AIS IMC and GIA are superimposed in satellite gravimetry observations, i.e. a spatially resolved parametrization of these signals is strongly correlated, and a signal separation appears challenging. For this reason, we additionally introduce satellite observations from ice sheet altimetry over the AIS, which are sensitive to these signals as well. Furthermore, we make use of products from regional climate and firn modelling to account for ice sheet surface processes. Thus, this work is an advancement of the work from e.g. , , , and into a global framework. A globally consistent approach, however, requires a parametrization of far-field effects and cannot treat them by applying regional constraints . Relevant far-field effects, from the perspective of Antarctica, are the following global mass redistributions: (Northern Hemisphere) GIA , IMC of the Greenland ice sheet, glacier mass changes, and terrestrial hydrological mass changes.
We formulate the following observation equation including three observational groups, denoted by , and six parameter types, denoted by :
1
The design block matrices, denoted by , include the parametrization; i.e. they link the observational groups and parameter types indicated by subscripts and superscripts in Eq. ().
The three data sets, denoted by , are (1) spherical harmonic coefficients of surface density changes from gravimetry, ; (2) a grid with surface elevation changes of the AIS, ; and (3) a grid with changes of the FAC of the AIS, . The grid definition is chosen according to the grounded part of the AIS and peripheral glaciers from . The six parameter types are (1) GIA, ; (2) ice mass changes of the Antarctic ice sheet, ; (3) changes of the firn air content of the Antarctic ice sheet, ; (4) ice mass changes of the Greenland ice sheet, ; (5) ice mass changes of glaciers outside Antarctica and Greenland, ; and (6) mass changes of non-glaciated water on the continent (e.g. groundwater, surface water), here referred to as hydrological mass changes, .
Here, our intentional goal is the incorporation of error covariance information, which is a more rigorous approach to address the observational errors than minimizing error effects in the data sets by filtering. Furthermore, the large-scale fingerprints are not sensitive to small-scale errors, such as the typical GRACE and GRACE-FO stripe patterns. The error covariance matrix, , of the observations is 2
We introduce a GIA parametrization in and to account for GIA effects in Antarctica and for GIA effects outside Antarctica. In Antarctica, we implement globally consistent GIA fingerprints calculated in response to local deglaciation impulses. These impulse response patterns are generated using the GIA modelling software SELEN and enable us to capture GIA effects independently from GIA forward models with an effective spatial resolution of 450 km . This parametrization does not imply a particular spatial ice loading history that would predetermine the spatial occurrence of present-day GIA effects.
An exception to forward-model independent GIA parametrization is made on the Antarctic Peninsula. From our validation experiments, we found that we were not able to retrieve reasonable GIA results for the northern part of the Antarctic Peninsula (Graham Land). We attribute this mainly to the insufficient quality of surface elevation changes derived from radar altimetry here (e.g. ). In turn, the significant misfit between GRACE and GRACE-FO products and CryoSat-2 products is captured by an unphysical GIA signal. This is also the case for other inverse GIA estimates (e.g. ). To prevent an unphysical GIA, we decided not to co-estimate GIA in this particular region. We did not include local GIA patterns on the peninsula in our local GIA pattern parametrization. Instead, we approached the GIA effect here by a global GIA model result, which was then subtracted from the observations. We modelled the GIA effect with an ICE-6G ice history that solely exists in the Graham Land region by using SELEN . Figure S1 illustrates the modified GIA parametrization with the Antarctic Peninsula GIA pattern. Admittedly, this GIA pattern has strong limitations in representing the true GIA effect in this region. The upper-mantle viscosity is found to be low here . We therefore expect that GIA response timescales are similar to those in the Amundsen Sea Embayment region. This means that the applied pattern (Fig. S1) will only allow an incomplete representation of the actual GIA and will not resolve GIA effects induced by load changes over the last few centuries. Nevertheless, we argue that this methodological adjustment allows us to, at least, limit the bias to the entire Antarctic GIA estimate.
Outside Antarctica (“far-field”), we use four globally consistent GIA fingerprints generated with regionally tailored ice loading histories for Greenland, Laurentia, Fennoscandia, and other regions (e.g. Patagonia, Barents Sea, and Kara Sea) similarly to . Each fingerprint is modelled with SELEN using an ICE-6G ice history and VM5a rheology. Additionally, we include two GIA fingerprints to capture a potential residual rotational feedback signal given by GIA modelling limitations attributed to an erroneous lower-mantle viscosity .
The parametrization of Antarctic ice mass changes in the case of altimetry observations, , is realized by linking the IMC in a grid cell with its corresponding surface elevation change. This involves, on the one hand, the elevation change associated with a volume change of pure ice (assuming a density of 917 kg m) and, on the other hand, the elevation change caused by the elastic deformation due to the mass change. The latter is obtained while generating the IMC parametrization of the gravimetry observations. is created by assuming a point mass change in the centre of each grid cell. We represent this point mass change by a set of spherical harmonic coefficients . The sea level response of each (ice) point mass change is calculated by solving the sea level equation . The point mass change together with its sea level response is the globally consistent fingerprint of the ice mass change in each grid cell. The elastic deformation effect of this point mass change is co-calculated when solving the sea level equation and then implemented in .
Since the grid definitions of the observations and of the parameters are identical, and are identity matrices. The change in FAC is mapped one-to-one in the altimetry observations. Note that the surface elevation changes observed by satellite altimetry are parametrized in two parts: the first part is the elevation change associated with ice density and the related elastic deformation, and the second part is the change in FAC.
The gravitational field changes include, as mentioned above, far-field effects (far-field from the perspective of Antarctic mass changes). To take these into account, we introduce GIA parameters to account for GIA effects from outside Antarctica, as mentioned above. Furthermore, we introduce parameters for glacier mass changes, , and for continental hydrology mass changes, . The matrices , , and link observed surface density changes with IMC of the GIS, glaciers, and hydrology, respectively. Glacier mass changes and continental hydrology mass changes are parametrized with 68 and 60 globally consistent fingerprints, respectively (updated according to ). We extend the hydrology parametrization with a fingerprint to capture the residual hydrological mass change signal over the continents which is not resolved by the applied hydrology parametrization (see Sect. S1 in the Supplement). The parametrization of IMC in Greenland applies 16 fingerprints for the eight drainage basins of the Greenland ice sheet . For this purpose, each basin is divided into a part below and above 2000 m surface elevation, which leads to a total of 16 sub-basins. The mass change of each sub-basin is not assumed to be uniform. Instead, a more realistic mass change pattern within each sub-basin is chosen based on mean rates of surface elevation changes derived from CryoSat-2 satellite altimetry (updated according to ). The globally consistent mass change pattern (fingerprint) of each sub-basin is calculated by solving the sea level equation.
The parameters are estimated by generalized least-squares adjustment (e.g. ). We apply a variance component estimation to further optimize the estimates by a relative weighting of the uncertainty information of the three observational groups gravimetry, altimetry, and FAC changes (Sect. S2). Additionally, we implement a Tikhonov regularization of the Antarctic GIA parameters to prevent unphysical GIA results due to limitations of the error covariances of the input data sets (Sect. S3). We solve the following regularized normal equations: 3 where is the estimated parameters and is the corresponding covariance matrix. is the normal equation matrix, is the right-hand side, and is the regularization matrix . We determine the degree of regularization using the L-curve criterion . Sections S1–S3 of the Supplement include further details about the implemented parameter estimation strategy and closed-loop simulation results to justify the regularization. Furthermore, Sect. S3.2 in the Supplement provides information on how we choose the preferred inversion solution presented in Sect. .
2.2 Data setsThe observations, denoted by , are mean rates according to the time period from January 2011 until December 2020 (10 years). We use the gravitational field changes ITSG-Grace2018 , which are GRACE and GRACE-FO Level-2 products provided as monthly sets of spherical harmonic coefficients up to degree 96 without any GIA correction. We express the gravitational field changes as mass changes in a spherical layer, termed surface density changes . These Level-2 products have a low noise level compared to other products and at the same time almost completely retain the signal . The gravitational fields are complemented by degree-1 products derived according to . coefficients, and coefficients in the case of GRACE-FO and GRACE accelerometer failures, are replaced with satellite laser ranging products . We do not apply any filter to the gravitational fields. It should be noted that GRACE and GRACE-FO Level-3 products, e.g. mascon solutions, are not suitable for the investigation presented here due to the following reasons: (1) mascon solutions are already corrected for the GIA effect, i.e. this GIA correction would have to be back-processed; (2) the globally consistent parametrization cannot be applied to Level-3 data and would have to be completely re-developed; and (3) the rigorous propagation of covariance information would not be possible unless it were available along with the Level-3 products.
The surface elevation changes are derived from updated CryoSat-2 products according to . Finally, the FAC changes are derived from the RACMO2.3p2 SMB product and the IMAU-FDMv1.2A firn thickness change product . Surface elevation and FAC changes are resampled to a grid of 50 km 50 km by averaging over a grid cell.
Figure 1
Mean rates of data sets for the time period January 2011–December 2020. (a) GRACE and GRACE-FO-derived surface density rate using ITSG-Grace2018 monthly gravitational fields (no filter applied). (b) CryoSat-2-derived surface elevation rate updated according to . (c) Thickness change of firn air content (FAC) derived from RACMO2 SMB and IMAU-FDM .
[Figure omitted. See PDF]
The uncertainty characterization follows . is derived from normal equations provided along with the ITSG-Grace2018 gravitational field products. The degree-1 uncertainty is characterized analogously to based on a degree-1 ensemble over the period under investigation. In the case of and , we base the uncertainty characterization on ensembles of surface height rates and FAC rates. The ensemble of altimetry-derived surface elevation changes is based on surface elevation rates derived with altimetry processing techniques, commonly denoted as retrackers, “EWIDTH”, “ICE1”, “OCOG”, and “TFMRA”. We do not include retrackers in the ensemble that are known to lead to elevation changes with poor quality. The uncertainty characterization of FAC is based on mean rate differences of cumulated surface mass balance anomalies derived from the regional climate models RACMO2.3p2 and MARv3.11 . showed striking differences between products of these two climate models indicating systematic errors. To avoid artificial downweighting of the FAC information within the inversion by too-conservative uncertainties, we limited the difference between MARv3 and RACMO2. We specified that the magnitude of the mean rate difference between RACMO2 and MARv3 of 1 pixel should not be greater than the mean rate magnitude derived from RACMO2 for this pixel. Thus, we set a threshold of 100 % for the maximum deviation in the ensemble of mean rate differences.
2.3 Assessment methodsWe use three approaches to assess the quality and soundness of the results obtained with the described methodology: (1) we compare the resulting GIA estimates with independent data derived from GNSS observations, (2) we compare the Antarctic GIA result with results from an alternative inverse approach as well as from forward modelling, and (3) we perform sensitivity tests with alternative input data sets.
In Assessment (1), we use preliminary results from a consistent analysis of more than 270 GNSS stations distributed over the entirety of Antarctica. This analysis has been accomplished in the frame of the SCAR-endorsed Geodynamics In ANTarctica based on REprocessing GNSS dAta Initiative (GIANT-REGAIN; ). Figure S8 illustrates the locations of the GNSS sites. GIA-related bedrock motion rates from the GNSS data are compared with bedrock motion rates from the inversion results. The bedrock motion time series from the GNSS data are corrected for elastic deformation effects of the solid Earth due to the IMC during the observation period of each GNSS site. We determine the IMC using the surface elevation change time series derived from satellite altimetry and from the firn model IMAU-FDM . Based on these IMCs, we calculate the elastic-related (vertical) bedrock motion by using the Green's function approach in the spatial domain . We use load Love numbers derived from the preliminary reference Earth model (PREM; ). The elastic bedrock motion effect, which is part of the observed surface elevation changes , is accounted for by assuming it to be 1.5 % of the altimetry-derived surface elevation change . It has to be emphasized that the comparison of our GIA estimates with those inferred from GNSS serves only the purpose of an independent assessment and that this data set is not part of the inversion framework. We calculate the weighted root-mean-square difference (WRMSD) between the GNSS-derived rates and those from the inversion results (INV), as done by :
4 with the weight, , for each GNSS site, , as 5
indicates the uncertainty (standard deviation) of the estimated rate. We derive from (Eq. ). We obtain from the GNSS processing, and we additionally assume 10 % of the estimated elastic deformation as its uncertainty.
In Assessment (2), we compare the Antarctic GIA estimate to the modelling result from and to the regional inverse estimate from . modelled the present-day GIA effect based on assumptions on the ice loading history and the solid Earth's rheology. They applied a Bayesian inversion approach to find a best-fit GIA model with GNSS observations as well as relative sea level records. This framework includes GNSS observations also in Antarctica from . used a similar data-driven approach, as presented here, to isolate the GIA effect from satellite gravimetry and satellite altimetry. However, the approach is a regional approach, as it regionally constrains GIA by calibrating interim results over a low-precipitation zone. Furthermore, it differs from the result presented here, as it utilizes input data from GRACE and ICESat over the time period February 2003 to October 2009. incorporated GNSS observations within the estimation procedure to justify parametrization choices, thus making this approach not fully independent from GNSS observations.
For Assessment (3), we run sensitivity tests by using alternative input data sets for the inversion. We use surface elevation changes from as an alternative altimetry product. During the investigated time period, the input data to this alternative product are predominantly from CryoSat-2. The ICESat-2 observations included start only in October 2018. Envisat observations are not included for periods later than September 2010, that is, not included in the period considered here. Although dominated by the CryoSat-2 input, the alternative surface elevation change estimates by result from a different processing scheme and include data from one alternative altimetry mission to some extent. The uncertainty information provided along with the surface elevation time series from is less comprehensive than our error characterization for CryoSat-2 products. Since the data set is largely based on CryoSat-2 data during our investigation period, we consider it reasonable to assume the same uncertainty information as we use for the CryoSat-2-only data set (Sect. ). Moreover, we test the sensitivity of the results to an FAC variant by exchanging RACMO2.3p2 SMB with MARv3.11 SMB. Since our ensemble for error characterization was created from differences between MARv3 and RACMO2 SMB, we assume the same uncertainty information as we apply for FAC changes based on RACMO2.3p2 SMB. In the case of gravitational field products, we use solutions from the Centre for Space Research of the University of Texas at Austin, CSR RL06.1 products . For CSR RL06, the predecessor of CSR RL06.1, found the lowest noise level among the Science Data System (SDS) solutions and fair signal retainment. Nevertheless, the noise level of CSR RL06.1 is higher than in ITSG-Grace2018 gravitational fields (Fig. S9d). Full error covariance information, likewise with ITSG-Grace2018, is not provided as a standard product along with CSR RL06.1 products. Pragmatically, we assume ITSG-Grace2018 uncertainty information for CSR RL06.1 gravitational fields here by accepting that the sensitivity test is inconsistent to some extent. Figure S9 provides maps of the alternative input data sets and differences from the input data sets described in Sect. .
3 ResultsFigure illustrates the following results from the preferred inversion solution where we apply a regularization parameter of : the Antarctic GIA-related bedrock motion (Fig. a), the IMC of the AIS expressed in terms of surface density change (Fig. b), and the estimated thickness change of FAC (Fig. c). Figure d–f show the associated 2 uncertainties derived from (Eq. ). We intentionally show in Fig. a the present-day GIA effect in terms of surface density changes, rather than smoother geoid height changes as shown elsewhere , in order to demonstrate the limits of the spatial resolution of present-day GIA effects with the inverse approach applied here.
Figure 2
Maps of Antarctica illustrating the estimates of (a) vertical bedrock motion due to glacial isostatic adjustment (GIA), (b) surface density change due to ice mass change (IMC) of the Antarctic ice sheet (AIS), and (c) the thickness change of the firn air content (FAC) derived from the preferred inversion solution. Panels (d)–(f) show each of the uncertainties, respectively. Units indicated for (d)–(f) apply columnwise to (a)–(c).
[Figure omitted. See PDF]
The integrated values are 86 21, 144 27 Gt a, and 13 18 km a in the case of AIS GIA, AIS IMC, and AIS FAC, respectively (Table ). We quantify the apparent mass effect of GIA in order to demonstrate the effect that GIA-induced gravitational field changes would have on gravimetry-only ice mass balance estimates. This quantification depends on the adopted method to infer ice mass changes from gravitational field changes. Here we convert GIA-related gravitational field changes to equivalent surface mass density changes, as we do with gravimetry observations (Sect. ; ; ). Subsequently, we integrate over the ice sheet region extended by a 400 km buffer zone. Note that different methodologies of gravimetric IMC inferences would imply different ways of integrating the GIA-equivalent surface mass density change and that our adopted scheme is not used for our actual IMC estimates but exclusively for expressing our estimated GIA signals in terms of an integrated mass effect. Figure S7 illustrates the Antarctic integrals for all regularization parameters. The FAC change integrated over the AIS is 4 and 13 km a for the input data and the estimate of the preferred inversion solution, respectively. Thus, the input value is still within a uncertainty interval of the estimate ( 18 km a). In view of the uncertainty, it is not possible to conclude whether the mean FAC rate for the entire grounded AIS is positive or negative during the 10-year time interval.
Figure 3
Assessment (1): comparison of GIA-related bedrock motion at GNSS sites from the preferred inversion solution (a), from GNSS observations (b), and the difference between both (c).
[Figure omitted. See PDF]
In summary we find the following spatial features of IMC: prominent negative IMCs are evident in the Amundsen Sea region, in the Getz Ice Shelf region, and at the Totten and Denman glaciers (Wilkes Land). Positive IMCs were detected at Kamb Ice Stream, Ellsworth Land, Dronning Maud Land, Enderby Land, and to some extent Terre Adélie (cf. Fig. a for geographical names). Note that already the input data sets consistently reveal these spatial features for a large part (Fig. ), and the spatial pattern of the altimetry-derived mean rates (Fig. b) basically determines the spatial pattern of the determined AIS IMC (Fig. a). The spatial pattern of FAC change, which enters the inversion as an input data set, is for a large part identical to the pattern of FAC change that is estimated.
For Assessment (1) of these results (Sect. ), Fig. provides a spatial comparison of the estimated GIA bedrock motion of the preferred inversion solution and the GNSS-derived bedrock motion rates. The WRMSD (Eq. ) is 4.9 mm a. For Assessment (2), Fig. illustrates maps of the GIA estimate from the preferred inversion in comparison with the GIA results from and .
Figure 4
Assessment (2): the GIA-related surface density rate from the preferred inversion solution (a), from (b), and from (d). The deviation of the latter two from the preferred inversion solution is illustrated in (c) and (e), respectively. The integrated GIA mass effects are 86 Gt a (a), 146 Gt a (b), and 117 Gt a (d) using a 400 km offshore buffer zone . In (a) the following geographical names are labelled: Amundsen Sea region (ASR), Antarctic Peninsula (APIS), Dronning Maud Land (DML), Ellsworth Land (EL), Filchner–Ronne Ice Shelf (FRIS), Getz Ice Shelf (GIS), Graham Land (GL), Kamb Ice Stream (KIS), Ross Ice Shelf (RIS), Terre Adélie (TA), Transantarctic Mountains (TM), and Wilkes Land (WL).
[Figure omitted. See PDF]
In agreement between the GIA result of the preferred inversion solution and the alternative GIA inverse estimate from (Fig. b) and the optimized forward modelling result from (Fig. d) is the bedrock uplift in Ellsworth Land, and somewhat in the Ross Ice Shelf region, the Filchner–Ronne Ice Shelf region, and Wilkes Land (cf. Fig. a for geographical names). However, there are differences in the determined magnitudes (Fig. c and e). The comparison reveals some common features of the inverse estimates (Fig. b; ) that are not found in the forward modelling result (Fig. d; ). These are the bedrock uplift in the Amundsen Sea region, the bedrock uplift in the Transantarctic Mountains, and the bedrock subsidence in the Kamb Ice Stream area. However, between the two GIA inverse estimates the local magnitudes and the spatial assignments of features differ (Fig. c). The bedrock subsidence indicated by the preferred inversion solution in Dronning Maud Land and Terre Adélie is not part of the GIA results reported by and . identified GIA-induced bedrock subsidence at the Antarctic Peninsula and in the Getz Ice Shelf region that is not part of the GIA estimate of the preferred inversion solution and of the GIA forward modelling result .
Table 1Comparison of integrated GIA mass effect in Antarctica (Antarctic GIA), integrated ice mass change of the Antarctic ice sheet (AIS IMC), and firn air content volume change of the Antarctic ice sheet (AIS FAC) that result from the sensitivity tests using alternative input data sets. The “preferred solution” data set shows results from the preferred inversion solution based on ITSG-Grace2018 gravitational fields, CryoSat-2-derived surface elevation changes, SMB from RACMO2.3p2, and firn thickness changes from IMAU-FDM. The “alt_JPL” data set is based on surface elevation changes from . The “smb_MAR” data set uses FAC changes derived from MARv3.11 SMB . The “grav_CSR_ITSG_err” data set utilizes gravitational fields from CSR RL06.1 and includes the uncertainty information from ITSG-Grace2018. Values in parantheses show the deviation from the “preferred solution” result (first row). Indicated uncertainties are values derived from Eq. (). Figure provides maps of the sensitivity results.
Antarctic GIA | AIS IMC | AIS FAC | |
---|---|---|---|
in Gt a | in Gt a | in km a | |
preferred solution | 86 21 | 144 27 | 13 18 |
alt_JPL | 101 21 | 162 26 | 9 18 |
(15) | (19) | (4) | |
smb_MAR | 85 22 | 142 28 | 1 25 |
(1) | (1) | (14) | |
grav_CSR_ITSG_err | 61 32 | 113 45 | 9 25 |
(25) | (30) | (4) |
The WRMSD (Eq. ) between GNSS observations and the GIA estimate from the preferred inversion solution (Fig. ), from , and from is 5.3, 6.8, and 6.6 mm a, respectively. Note that we transferred the GIA-related surface density change from and to bedrock motion using a GIA density mask similar to that of . These three WRMSDs do not include any weight for GIA-related uncertainties but only include weight for GNSS-related uncertainties, because there is no consistent GIA uncertainty information for the three GIA estimates available. By not including the GIA uncertainty in the weights, we ensure that all three GIA models are treated equally for comparison with the GNSS. For this reason the WRMSD values differ from the WRMSD values given above and from the values illustrated in Fig. S5.
For Assessment (3), Fig. provide maps of the sensitivity test results in Antarctica. Differences between integrated results from the preferred inversion solution and integrated results from the sensitivity tests are smaller than the estimated uncertainties. Table summarizes the results of the sensitivity tests. Integrated values of the Antarctic GIA mass effect, the AIS IMC, and the change in FAC are compared. In each experiment, one of the three data sets has been substituted by an alternative data product.
Figure 5
Assessment (3): results from sensitivity tests using alternative input data sets (Fig. S9). The first row (a–c) shows the results from the preferred inversion result (Fig. 3). The second row (d–f) shows the differences between the inversion results and the first row where the surface elevation rate from (Fig. S9b) is used (test result minus preferred inversion solution). Similarly, the third row (g–i) and fourth row (j–l) show differences from the inversion results with MARv3.11 SMB data (; Fig. S9c) and CSR RL06.1 gravitational fields (; Fig. S9a), respectively.
[Figure omitted. See PDF]
4 Discussion4.1 Assessment and comparison
The AIS mass balance of 144 27 Gt a (2011–2021) from the preferred inversion solution is in the range from 94 to 202 Gt a (2010–2019) given in the Sixth Assessment Report of the IPCC . Our estimated AIS contribution to global mean sea level for 2011–2021 is 0.40 0.07 mm a. This is very close to the result from the most recent Ice sheet Mass Balance Inter-comparison Exercise (IMBIE) study at 0.40 0.09 mm a for the same decade. It is worth noting that the Antarctic result from is based on a total of 23 different ice mass balance estimates. The estimated integrated GIA effect of 86 21 Gt a is at the upper limit of integrated values presented by and but lower than the results by and . The stated uncertainties are also plausible. The differences of the sensitivity results from the reference result of the Antarctic-wide integrated values are always smaller than the uncertainties derived from the estimate (Table ). From this, we conclude that the accuracy of the presented integrated results is sound.
As described in Sect. , the comparison of the preferred solution with alternative GIA results reveals some similarities but also prominent differences. The GIA result of the preferred solution fits better to GNSS observations than the GIA solutions by and (Fig. ). However, GNSS-derived bedrock motion is only available for some parts of the Antarctic continent. The comparison is therefore always subject to the asymmetry given by the spatial coverage of GNSS data. In addition, the GIA result of the preferred solution significantly underestimates the bedrock uplift observed with GNSS in the Amundsen Sea region (Fig. c), which is presumably more realistically imaged in the result by .
Other data combination approaches, which aim to estimate present-day GIA effects, found significant GIA-induced bedrock subsidence in the Getz Ice Shelf region. Such subsidence is apparent in the GIA estimates according to , (somewhat offshore in Fig. b), and (clipped by the choice of the colour bar limits in Fig. 3a in , but visible in Fig. 2f in ). Before we implemented the IMC and FAC parametrization in the peripheral glacier regions, we also obtained this negative anomaly (Fig. 6.3 in ). This negative GIA anomaly vanishes by extending the IMC and FAC parametrization to include the peripheral glacier regions.
4.2 Methodological implications
In simulation experiments by , GIA could be spatially resolved without filtering or regularization. This requires profound knowledge about error covariances of the input data sets. We find that the used error covariances of the available input data sets (Fig. S12) are limited, and it is not useful to determine realistic Antarctic GIA effects (Fig. S6d and g) solely by relying on the error covariance information. In particular, Fig. S12e illustrates that the error characterization of the altimetry trends, based on a data processing ensemble, leads to the result that the errors are strongly correlated on a continental scale; i.e. they represent a bias. It is realistic that the error covariance information comprises biases, but whether we capture them sufficiently remains questionable. Note that it is not possible to account for locally or regionally limited errors by including these continental-scale error patterns in the parameter estimation.
However, with further simulations (Sect. S3.1) we demonstrated that a regularization can help to somewhat compensate this lack of knowledge and that it is possible to derive physically plausible results, especially in terms of integrated values (Fig. S3). Note that the found regularization optimum from the simulations cannot be applied to define the regularization parameter in the real-data case. This parameter is defined by the L-curve criterion discussed in the next paragraph. The simulation demonstrates that the noise level of the results is high, especially in the case of IMC (Fig. S4k). Applying the regularization to Antarctic GIA parameters and neglecting correlated errors of the input data sets are the main methodological limitations of the work presented here. In order to spatially resolve GIA independently from forward models, we need to accept for the moment an enhanced sensitivity towards input data errors.
For the results presented here, we can avoid classical filtering of the input data, such as Gaussian smoothing or decorrelation filtering (e.g. ) in the case of GRACE and GRACE-FO data, as it has been done in other data combination studies (e.g. ). It is not necessary to equalize the spatial resolution of the data sets prior to the joint inversion; i.e. we can avoid coarsening the spatial resolution of the input data sets. As discussed above, we implement a regularization of the Antarctic GIA parameters to prevent dominant spatial oscillations that otherwise appear in the GIA result and cannot be physically justified (Fig. S6). Like Gaussian filtering, regularization is methodologically less advanced than capturing errors by including the error covariance information. Nevertheless, we can justify the choice of a regularization parameter of (the amount of damping) with two arguments in the real-data case. First, it can be justified by the L-curve criterion (Fig. S5a). Second, the bend in the L-curve coincides with inversion solutions that show the smallest deviation from independent GNSS observations in terms of WRSMD (Eq. ). Despite the regularization, the GIA result of the preferred inversion solution shows spatial oscillations (Fig. a), which are anticorrelated to the IMC result to some degree (Figs. b, S14).
The implemented hydrological residual fingerprint (Fig. S2) allows us to capture possible far-field effects due to the limitations imposed by the imperfect hydrology parametrization. Evaluated over the Antarctic continent with a 400 km buffer zone, the integrated mass effect of the hydrological residual fingerprint amounts to 7.4 Gt a. Nevertheless, applying this fingerprint could be only an interim solution that we use here in an Antarctica-focused study. As soon as an improved globally consistent hydrology parametrization is available, this caveat can be remedied.
4.3 Interpretation
4.3.1 GIA estimate in East Antarctica
Limitations in spatially resolving GIA in Antarctica are indicated by the anticorrelation of some patterns of the IMC result and the GIA results (Figs. a and b, S13 and S14). For East Antarctica (with its rheology favouring GIA response times of millennia), we do not expect such anticorrelation for the actual signal of IMC and GIA, because such anticorrelation would require an associated correlation between patterns of deglaciation on millennial timescales and present-day IMC. Rather, the resolved GIA and IMC patterns in East Antarctica (Fig. a and b) indicate spatial error patterns propagated from the input data (Figs. a and b, S11e, ). The sensitivity experiments (Assessment 3) show that the GIA signal in Terre Adélie and Wilkes Land (Fig. ) depends on the choice of the altimetry product. In Wilkes Land it also depends significantly on the choice of the gravimetry product. This is obviously due to the non-consideration of correlated errors within the parameter estimation. Moreover, the simulation experiments (Sect. S3.1), where we test the regularization of Antarctic GIA, reveal that correlated altimetry errors are obviously reflected in the GIA and IMC result (Fig. S4j and k). This is also evident from the larger RMSE we find for AIS IMC than for the experiment where we have full knowledge of error covariance information. Nevertheless, the integral is very close to the simulated truth. From this we conclude that the preferred inversion solution still contains GIA and IMC patterns, which are artefacts due to data quality limitations rather than resolved physical GIA and IMC signals. This means one should be cautious when interpreting the short-scale spatial GIA patterns in East Antarctica by physical means. Nevertheless, predictions from GIA forward modelling disagree here, too, because there is a lack of knowledge in ice loading history and the rheological structure . Given this, it is currently challenging to ascertain how significant the identified spatial patterns of GIA-related bedrock motion in the interior of East Antarctica are in terms of physics. This implies that we cannot decide how useful our East Antarctic GIA estimate is as boundary information for testing glacial histories or rheological models. Nevertheless, the integrated values from the GIA estimate and thus the large-scale effects may hold some promise for this task. Measurements of the bedrock motion beneath the East Antarctic ice sheet would be helpful as independent information.
4.3.2 GIA estimate in West Antarctica
The spatial resolution capability of the chosen parametrization is, in a best case, based on reasonable physics. However, if the parametrization is at a finer resolution than the resolution capability of the data allows, this leads to overfitting in the inversion. demonstrated that the GIA parametrization applied here, which was chosen in agreement with the spatial resolving capability of solely GRACE and GRACE-FO data, is not able to resolve GIA effects associated with low upper-mantle viscosity and ice loading changes over the last few centuries. Such GIA effects, as postulated for the Antarctic Peninsula and the Amundsen Sea region , require a spatial resolution capability of 100 km (gravitational fields up to degree 200). The GNSS comparison with the GIA result of the preferred inversion solution illustrates the apparently limited GIA imaging capability within the Amundsen Sea region (Fig. ). Furthermore, the regularization dampens the GIA signal. In summary, with the inversion approach presented here, we are not able to fully spatially resolve GIA effects associated with low upper-mantle viscosity. What we present here is de facto a smoothed version of the true GIA signal. For comparative studies with forward modelling results that aim to represent the realistic rheological structure in West Antarctica, the comparison of smoothed results could at least help to constrain the parameter space. A high-resolution observation-based determination of GIA in this region remains a task for future work. This holds for the Antarctic Peninsula, too, where the GIA result presented here equals a classical GIA modelling result as described in Sect. .
4.3.3 IMC estimate
The spatial patterns of IMC are essentially determined by satellite altimetry (Sect. ), which enables the high spatial resolution based on the selected parametrization. It is noteworthy that this IMC estimate was determined consistently globally and reconciles GRACE and GRACE-FO and CryoSat-2 observations in a least-squares sense. Satellite gravimetry and altimetry are traditionally used separately to determine IMC and have differed significantly in IMBIE assessments . The result presented here is also in excellent agreement with the estimate from the statistical analysis of 23 different mass balances assessed in IMBIE. This lends confidence to our results and, hence, to the applied method. However, it is noteworthy that, as mentioned in Sect. , the spatial IMC features are partly anticorrelated with some of the found GIA features which we deem unphysical. This is also reflected in the results of the sensitivity experiments (Fig. ). This means that not all resolved spatial patterns can be interpreted as IMC. Based on the sensitivity experiments, we conclude that the stated uncertainty of 27 Gt a () is realistic. However, this uncertainty is still large considering that it amounts to almost 20 % of the magnitude.
4.3.4 FAC estimate in context of its uncertainty
We apply the FAC uncertainty information which assumes that differences between RACMO2.3p2 and MARv3.11 SMB products represent the true modelling error and can be used to characterize the SMB uncertainty. If we apply this empirical uncertainty information, this leads to unphysical GIA artefacts. In addition to ignoring correlations, we constrain the characterization of uncorrelated FAC errors. We presume that the empirical FAC uncertainty information is not fully sufficient to account for the true but unknown FAC error. MARv3.10 SMB especially shows a striking difference from the ensemble mean SMB in the (leeward of the) Transantarctic Mountain region (Fig. 6f; ), where we found unphysical GIA in preliminary results. As we use differences between RACMO2.3p2 and MARv3.11 SMB products to characterize the FAC uncertainty, this ends up in a large empirical uncertainty assumption here. In turn, the uncertainty assumption in this region allows unrealistic liberty within the inversion framework to explain the data. In fact, we presume that the spatial pattern of the differences (Fig. 6f; ) propagates to the GIA estimate presented here. For this reason, we constrain the mean rate ensemble from which we derive the FAC uncertainty as described in Sect. . Future studies may show the degree of improvement of FAC changes that can be achieved if a more sophisticated uncertainty characterization of FAC is available . Systematic SMB modelling errors, however, only explain part of the unphysical GIA effects of the preliminary results, as these also occur if the error covariance information of the other data sets is incorporated.
4.4 Outlook
If there is no improved error covariance information available, the spatial error patterns in the results could also be damped, for example, by applying Gaussian smoothing to the input data (instead of applying a regularization). This may lead to IMC results comparable to usual gravimetric mass balances but would, however, smooth the entire result. An alternative strategy may be to adjust the regularization depending on the region, e.g. by implementing a stronger regularization of the GIA effects within East Antarctica where they are presumably small . This needs to be justified based on additional information. Likewise a GIA parametrization, which is more oriented towards forward modelling results, may be used, but it has been the very intention here to make the estimation independent from possible GIA modelling errors.
The implementation of an extended IMC parametrization could further optimize the inversion result. A parametrization that allows for fine spatial resolution would only be desirable where it is justified by the signal-to-noise ratio (SNR) of the input data. The spatial resolution given by the parametrization may be fine where the SNR is large and coarse where the SNR is small. Further, fine spatial resolution is only needed where mass change processes occur on small spatial scales, i.e. in particular the ice sheet margin. The IMC estimate (Fig. b) indicates that the inversion is good at spatially separating large IMC amplitudes, e.g. at the ice sheet margin where ice dynamic flow changes govern IMC. This is probably not necessary in the East Antarctic interior, where small-scale IMC changes are less likely to be relevant. Such an adapted IMC parametrization could help to reduce the presumed error patterns in East Antarctica (Fig. S12e), which are currently erroneously assigned to a GIA effect in the inversion (Sect. ).
We expect a significant quality improvement of satellite-altimetry-derived surface elevation changes with new retracking methods in the case of radar altimetry and with the growing availability over time of laser altimetry products from the ICESat-2 mission. In terms of mean rates, the quality of the results in general will grow by investigating longer time periods. For instance, it is expected that the enhanced data products will facilitate the observation-based assessment of the GIA effect on the Antarctic Peninsula, which was not achieved in this study.
In a next step, the approach used here could be extended so that IMC changes can be resolved in time and not just as mean rates over defined time periods. According to the input data set availability, monthly IMC and FAC changes may be estimated. This requires characterizing uncertainties on the same temporal scales.
Moreover, the applied methodology is designed to serve as a complement to the global inversion of all sea level contributions and may allow us to resolve issues while co-estimating the GIA component.
5 Conclusions
We demonstrated the successful application of a joint global inversion approach with a focus on Antarctica. It combines the advantages of data sets derived from GRACE and GRACE-FO, CryoSat-2, and regional climate and firn modelling. We claim that the results presented have the following advantages over previous studies: the estimation procedure preserves global consistency in its representation of mass changes, because a global framework can avoid regional constraints as implemented in previous inverse GIA investigations. The inversion enables us to spatially resolve GIA effects in Antarctica largely detached from GIA forward modelling constraints. In addition, it enables us to determine high-resolution (50 km) IMC. FAC changes are implemented in the parameter estimation procedure instead of taking them into account deterministically only. Lastly, the estimation procedure uses a weighting based on realistic input data uncertainties and thus also allows sound uncertainty estimates of the results.
We estimate the following Antarctic-wide integrated values over the 10-year time interval from January 2011 until December 2020: a present-day GIA mass change effect of 86 21 Gt a, AIS IMC of 144 27 Gt a, and volume change of FAC of 13 18 km a. IMC and FAC integrals include peripheral glaciers. The GIA integral includes a 400 km offshore buffer zone. From the comparison with other published AIS IMC and GIA results, from the comparison with independent GNSS observations, and from sensitivity tests, we conclude that the presented results are sound and that the provided accuracy is reliable. We found better agreement of the GIA results from the preferred inversion solution with independent GNSS observations than GIA results from other inversion studies that even incorporate GNSS observations in their estimation procedure.
and this study demonstrated the relevance of having profound knowledge on error covariances of the input data sets available. So far, we were not able to completely eliminate spatial error patterns in the results which propagate from the input data. Moreover, we can attribute error sources only to some extent.
We see potential for improvement of the approach applied here by advancing the global hydrology parametrization, including improved error covariance information, and further developing the IMC and GIA parametrization.
Data availability
GRACE and GRACE-FO monthly gravitational fields can be obtained via 10.5880/ICGEM.2018.003 ().
CryoSat-2 data can be obtained from
RACMO2 SMB and IMAU-FDM are available on reasonable request (
MARv3.11 can be found at
The supplement related to this article is available online at:
Author contributions
Contributions are according to CRediT. Conceptualization: MOW and MH. Data curation: VH, EB, BU, MOW, and MS. Formal analysis: MOW. Funding acquisition: MH, JK, and MS. Investigation: MOW. Methodology: MOW and MH. Software: MOW. Supervision: MH. Validation: MOW. Visualization: MOW. Writing (original draft): MOW. Writing (review and editing): MH, EB, VH, BU, JK, MS, and MOW.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
We thank Yu Sun and the two anonymous referees for their very helpful feedback which improved the article. The work of Matthias O. Willen and Bernd Uebbing was funded by the grants HO 4232/4-2 and KU 1207/22-2, respectively, “Reconciling ocean mass change and GIA from satellite gravity and altimetry (OMCG)” (project no. 313917204) of the Deutsche Forschungsgemeinschaft (DFG) as part of the Special Priority Programme (SPP) 1889 “Regional Sea Level Change and Society” (SeaLevel). The work of Eric Buchta was funded by the grant SCHE 1426/26-1 and 2 (project number 404719077) of the Deutsche Forschungsgemeinschaft (DFG) as part of the SPP 1158 “Antarctic Research with Comparative Investigations in Arctic Ice Areas”. We kindly thank all colleagues and institutions who provided geodetic GNSS data in Antarctica to the SCAR-endorsed Geodynamics In ANTarctica based on REprocessing GNSS dAta Initiative (GIANT-REGAIN) led by Mirko Scheinert and Matt King (University of Tasmania, Hobart, Australia).
Financial support
This research has been supported by the Deutsche Forschungsgemeinschaft (project no. 313917204, grant nos. HO 4232/4-2 and KU 1207/22-2; project no. 404719077, grant no. SCHE 1426/26-1).This open access publication was financed by the Saxon State and University Library Dresden (SLUB).
Review statement
This paper was edited by Masashi Niwano and reviewed by Yu Sun and two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
A detailed understanding of how the Antarctic ice sheet (AIS) responds to a warming climate is needed because it will most likely increase the rate of global mean sea level rise. Time-variable satellite gravimetry, realized by the Gravity Recovery and Climate Experiment (GRACE) and Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) missions, is directly sensitive to AIS mass changes. However, gravimetric mass balances are subject to two major limitations. First, the usual correction of the glacial isostatic adjustment (GIA) effect by modelling results is a dominant source of uncertainty. Second, satellite gravimetry allows for a resolution of a few hundred kilometres only, which is insufficient to thoroughly explore causes of AIS imbalance. We have overcome both limitations by the first global inversion of data from GRACE and GRACE-FO, satellite altimetry (CryoSat-2), regional climate modelling (RACMO2), and firn densification modelling (IMAU-FDM). The inversion spatially resolves GIA in Antarctica independently from GIA modelling jointly with changes of ice mass and firn air content at 50 km resolution. We find an AIS mass balance of
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details




1 Institut für Planetare Geodäsie, Technische Universität Dresden, Dresden, Germany
2 Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany
3 Institute of Geodesy and Geoinformation, University of Bonn, Bonn, Germany