Present-day secular variations in the low-degree harmonics of the geopotential: Sensitivity analysis on spherically symmetric Earth models

Size: px
Start display at page:

Download "Present-day secular variations in the low-degree harmonics of the geopotential: Sensitivity analysis on spherically symmetric Earth models"

Transcription

1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B12, 2378, doi: /2001jb000696, 2002 Present-day secular variations in the low-degree harmonics of the geopotential: Sensitivity analysis on spherically symmetric Earth models M. E. Tamisiea and J. X. Mitrovica Department of Physics, University of Toronto, Toronto, Ontario, Canada J. Tromp Seismological Laboratory, California Institute of Technology, Pasadena, California, USA G. A. Milne Department of Geological Sciences, University of Durham, Durham, UK Received 24 May 2001; revised 30 May 2002; accepted 18 July 2002; published 31 December [1] We predict present-day rates of change of the long-wavelength components of the Earth s geopotential due to the late Pleistocene glacial cycles (hereinafter referred to as glacial isostatic adjustment, or GIA). These predictions are generated using spherically symmetric, self-gravitating, (Maxwell) viscoelastic Earth models. Previous studies have generally considered only zonal (i.e., azimuthally independent) harmonics, or the so-called _J coefficients. We extend these efforts to focus on the nonzonal harmonics and explore the sensitivity of our predictions to changes in a variety of model inputs. As an example,we examine the influence on the GIA predictions of ignoring rotational perturbations and assuming an elastically incompressible model and find that these assumptions can have a significant impact on the degree-2, order-1 coefficients, while the other harmonics are generally affected by <10%. Though the GIA predictions are generally insensitive to variations in lithospheric thickness and upper mantle viscosity, reasonable variations in either lower mantle viscosity or the late Pleistocene ice history will have a more significant effect. Mindful of upcoming dedicated gravity missions, we show that predictions of the GIA-induced changes of the geoid can vary by at least 0.5 mm yr 1 over a plausible range of Earth and ice models. INDEX TERMS: 1213 Geodesy and Gravity: Earth s interior dynamics (8115, 8120); 1236 Geodesy and Gravity: Rheology of the lithosphere and mantle (8160); 8162 Tectonophysics: Evolution of the Earth: Rheology mantle; KEYWORDS: secular, geoid, GIA, Stokes coefficients Citation: Tamisiea, M. E., J. X. Mitrovica, J. Tromp, and G. A. Milne, Present-day secular variations in the low-degree harmonics of the geopotential: Sensitivity analysis on spherically symmetric Earth models, J. Geophys. Res., 107(B12), 2378, doi: / 2001JB000696, Introduction [2] Most previous studies of the present-day secular change of the Earth s geopotential (or geoid height) have focused on the low-degree zonal coefficients of its spherical harmonic expansion, with particular emphasis on the degree two case (i.e., _J 2 ). This focus reflects both the limited nature of estimates available from previous gravity missions [e.g., Yoder et al., 1983; Eanes and Bettadpur, 1996; Cazenave et al., 1996; Nerem and Klosko, 1996; Cheng et al., 1997; Klosko and Chao, 1998], and the mathematical connection between load-induced _J 2 and length-of-day variations [e.g., Wu and Peltier, 1984; Mitrovica and Milne, 1998]. Copyright 2002 by the American Geophysical Union /02/2001JB000696$09.00 [3] Glacial isostatic adjustment (GIA) refers to the response of the Earth system to loading associated with the late Pleistocene glacial cycles. The contribution of GIA to secular variations in the zonal harmonics of the Earth s geopotential has been the focus of nearly two decades of study [e.g., Peltier, 1983; Wu and Peltier, 1984; Yuen and Sabadini, 1985; Ivins et al., 1993; Mitrovica and Peltier, 1993; Vermeersen et al., 1998; Devoti et al., 2001]. Several analyses have revealed a detailed depth-dependent sensitivity of numerical predictions to variations in the radial profile of mantle viscosity [e.g., Mitrovica and Peltier, 1993; Devoti et al., 2001]; in the case of the lowest-degree harmonics this sensitivity extends into the deepest regions of the mantle. The sensitivity of the predictions to variations in the ice load has also been explored. Wu and Peltier [1984], for example, showed that the GIA-induced _J 2 signal was primarily sensitive to net melting of late Pleistocene ice sheets (or the ETG 18-1

2 ETG 18-2 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS equivalent eustatic sea level change associated with the deglaciation). Mitrovica and Peltier [1993] and Ivins et al. [1993] extended this sensitivity analysis to higher degrees and noted a strong sensitivity of predicted odd-degree zonal harmonics (particularly at degree 3) to variations in the ice model. In contrast to degree two, the sign of the Legendre polynomial basis function at degree three is different at the two poles of the Earth. As a consequence, the _J 3 prediction is sensitive to small changes in the balance between mass loss from northern and southern ice complexes. [4] The impact of present-day glacial forcings on the zonal harmonics has also been an area of active interest [e.g., Yoder and Ivins, 1985; Sabadini et al., 1988; Mitrovica and Peltier, 1993; James and Ivins, 1997]. In this case, the Earth s response is primarily elastic and hence uncertainties in the radial viscosity profile are largely irrelevant. At degree two the predictions are sensitive to the net ongoing mass flux from the global ice reservoirs, and thus an accurate prediction of this signal can be obtained without knowledge of the individual contributions to the net flux [Mitrovica and Peltier, 1993]. For the reason discussed above, this simplified requirement does not hold for odd degree zonal harmonics [Mitrovica and Peltier, 1993; James and Ivins, 1997]. [5] Nonzonal secular variations in the geopotential, induced by either GIA or ongoing ice mass changes, have received less attention. Yuen et al. [1987] predicted zonal and nonzonal harmonics for a range of different lower mantle viscosities and both Maxwell and Burger s body rheologies. Their GIA calculations did not vary the upper mantle viscosity or lithospheric thickness, and they adopted simplified disk load models of the Laurentian and Fennoscandian ice complexes. Ivins et al. [1993] extended this work by considering a set of global ice models, as well as regional subsets of each, and a suite of both layered and continuously (radially) varying viscosity profiles. Both studies assumed an incompressible rheology and uniform meltwater redistribution. [6] With the initiation of a new set of gravity missions, namely, CHAMP and GRACE, the accuracy of constraints on the spherical harmonic (zonal and nonzonal Stokes) coefficients of the secular change in the Earth s geopotential may be significantly improved [e.g., Dickey et al., 1997]. These impending constraints motivated the present study, in which we focus on a suite of sensitivity analyses related to the GIA signal. Recent GIA analyses have extended traditional numerical formalisms to include effects such as rotation-induced feedback [e.g., Milne and Mitrovica, 1996, 1998; Milne et al., 1999; Peltier, 1998, 1999]. We first explore the sensitivity of our predictions of the Stokes coefficients to several of these extensions. We next perform a more traditional sensitivity analysis in which we vary both the radial viscosity profile of the Earth model and the adopted late Pleistocene ice sheet history. This traditional analysis will involve both the individual Stokes coefficients and predictions of global maps of the secular geoid change. A companion paper extends the calculations performed here to consider Earth models that include lateral variations in mantle viscosity. [7] In practice, there are two possible applications of the predictions presented here. First, if all other signals associated with secular variations in the geopotential were known, then the predictions might be used to constrain rheology or ice history. Alternatively, one might use the predictions to correct observations of the secular changes in the geopotential for GIA and use the residual signal in some other geophysical analysis (e.g., constraining mass variations associated with ongoing hydrological, oceanographic and cryospheric processes). Our bias is toward the latter philosophy as long as the sensitivity analyses establish a plausible range for the correction, and thus a plausible uncertainty in the residual signal. 2. Secular Geoid Changes Due to GIA [8] In this paper we focus on perturbations to the geopotential and geoid (or sea surface) height. These perturbations are related via Brun s formula through a simple multiplicative factor equal to the surface gravitational acceleration [Lambeck, 1988]. Predictions of the geoid height anomaly due to GIA involve two distinct numerical steps. First, for a given viscoelastic Earth model and late Pleistocene ice history, the sea level equation is solved to determine a gravitationally self-consistent ocean load. In this study we adopt elements of the recent sea level theory outlined by Milne et al. [1999] for this purpose. In the second step, we convolve the total (ice plus ocean) surface load with the Green s function for the geoid height anomaly derived by Mitrovica and Peltier [1989]. This step also incorporates the influence of rotation perturbations on geoid height [Mitrovica et al., 2001]. All calculations are performed using spectral approaches with a truncation at spherical harmonic degree and order 64. [9] Let us denote the geoid height anomaly computed using the above procedure as G. In order to calculate the present-day (t = t p ) secular change of the geoid, _G, we use a finite difference method 1 _G q; f; t p ¼ G q; f; t p þ t G q; f; tp t ; ð1þ 2t as in the work by Mitrovica and Peltier [1993], where t was chosen to be 100 years and q and f are the colatitude and east longitude, respectively. _G can be expanded in terms of normalized Legendre functions ~P m [e.g., Chao and Gross, 1987] X _G 1 X q; f; t p ¼ a ¼0 m¼0 ~P m ðcos qþ : _C m t p cos mf þ _S m t p sin mf ; ð2þ where a is the radius of the Earth. The constants _C m and _S m are the Stokes coefficients, where is the degree and m is the order of the Legendre function. Past studies have concentrated on the values of the zonal harmonics, _J, which can be related to the values of _C 0 by p _J t p ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffi 2 þ 1 _C 0 t p : Also, _S 0 0. [10] Numerical predictions for the Stokes coefficients require the input of both an ice and Earth model. For the Earth model, we use elastic parameters given by PREM ð3þ

3 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS ETG 18-3 [Dziewonski and Anderson, 1981] and divide the mantle into two isoviscous layers, separated at the 670 km seismic discontinuity. The upper and lower mantle viscosities are denoted by n UM and n LM, respectively. In addition, the model includes an elastic lithosphere of thickness, h LT.In general, we will adopt a glacial cycle model based on the ICE-3G deglaciation history [Tushingham and Peltier, 1991] but modified to include a glaciation phase of duration 90 kyr. However, we also present predictions using a similarly modified version of the ICE1 deglaciation history [Peltier and Andrews, 1976]. Unless otherwise noted, two glacial cycles are used in all calculations. [11] In the sensitivity analyses described below it is impossible to present the complete set of Stokes coefficients for any reasonably broad degree range. (There are, for example, 49 such coefficients for spherical harmonic degrees up to 6.) As a result, in any specific figure we choose a subset (which will vary from figure to figure) that illustrates the main point of the calculations Effects of Common Simplifications [12] In Figure 1 we show predictions of various Stokes coefficients as a function of the lower mantle viscosity. The prediction of _C 20 has a form which is characteristic of numerous past predictions of the zonal harmonics when the upper mantle is fixed to Pa s (see also the _C 31 case in Figure 1). At low viscosities the predicted signal is small because these models will have relaxed significantly, nearly reaching equilibrium, in the time (5 kyr) since the end of the model deglaciation. Models with high viscosity have progressively more sluggish relaxation throughout the postglaciation phase, and hence these models also yield progressively lower present-day signals. The predicted signals are thus highest at some intermediate value of the lower mantle viscosity, and the result is the profile shown in the top frame of Figure 1. For _C 20, the peak predicted amplitude is obtained for n LM Pa s. [13] The results for the remaining, nonzonal harmonics in Figure 1 exhibit a variety of forms. The rather complicated sensitivity of the _C 21 and _S 21 predictions to variations in n LM is due to rotational effects. Changes in the rotation vector perturb the centrifugal potential, and this perturbation (known as the rotational potential) acts to load the planet. The rotational potential can have noticable effects on global geoid maps [Peltier, 1999; Mitrovica et al., 2001]. Because the rotational potential is dominated by true polar wander (TPW), it has a geometry governed by a spherical harmonic with = 2 and m = 1 [e.g., Han and Wahr, 1989]. Thus, one would expect the Stokes coefficients _C 21 and _S 21 to be perturbed the greatest due to this effect. This expectation is confirmed in Figure 1. Predictions of _C 21 and _S 21 without the rotational signal included (dashed lines) show the usual form associated with GIA results for nonrotating Earth models. The signal due to rotation is large enough to change the sign of the prediction relative to the nonrotating case. Furthermore, the sensitivity of the predictions to variation in n LM is similar to GIA-induced TPW speed predictions as a function of the same suite of Earth models [e.g., Yuen et al., 1986]. [14] Rotational effects are generally negligible for Stokes coefficients other than those at degree two and order one. Figure 1. Effect of rotational perturbations on predictions of various Stokes coefficients. Predictions are shown for a suite of lower mantle viscosities (n LM, on the abscissa scale), and in all cases the lithospheric thickness, h LT,is71 km and upper mantle viscosity, n UM,is Pa s. The elastic parameters are taken from PREM and the ice model includes two cycles of a modified version of ICE-3G (see text). The solid lines are predictions based on the theory for a rotating planet, while the dashed lines assume a nonrotating Earth (and thus ignore effects associate with perturbations in Earth rotation). The dotted line on the bottom frame is analogous to the solid with the exception that ocean loading effects are not included.

4 ETG 18-4 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS scale). Ocean loading effects contribute significantly in this case due to the geometry of the ocean basin. As a consequence, rotational perturbations to the ocean load also yield a nonnegligible contribution to the predicted signal. [15] A second simplification commonly employed in GIA studies is to assume that the Earth is elastically incompressible [e.g., Spada et al., 1992; Vermeersen and Sabadini, 1996]. This assumption has a particularly large effect on the _C 21 and _S 21 coefficients (Figure 2) and yields predictions Figure 2. The influence of assuming an elastically incompressible model on predictions of various Stokes coefficients. Predictions are shown for a suite of lower mantle viscosities (n LM, on the abscissa scale), and in all cases the lithospheric thickness, h LT, is 120 km and upper mantle viscosity, n UM,is10 21 Pa s. The elastic parameters are taken from PREM and the ice model includes one cycle of a modified version of ICE-3G (see text). The solid and dashed lines are predictions based on elastically compressible and incompressible Earth models, respectively. However, there are exceptions, and an example is given in the bottom frame of Figure 1. The dotted line on this frame shows a prediction for a rotating Earth in which all of the effects associated with the ocean load are removed. For _C 33 the total signal is very small (notice the ordinate Figure 3. Variation of selected Stokes coefficients with lithospheric thickness, h LT. Results are generated using the Earth models and ice history described in Figure 1 (n UM = Pa s), with the exception that we include the three cases: h LT = 71, 96 and 120 km (as labeled).

5 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS ETG 18-5 Figure 4. Sensitivity of _C m and _S m predictions up to degree and order 5 (as labeled) to variations in n UM and n LM. In all cases the lithospheric thickness, h LT, is 71 km. The elastic parameters are taken from PREM and the ice model includes two cycles of a modified version of ICE-3G. Each label at the top of the frame is proceeded by a minus sign if the prediction is negative and is followed by the multiplicative factor for the contour levels in parenthesis. The units are yr 1. The shaded overlay on each frame reflects the range of values for n UM and n LM inferred in recent analyses of GIA data (see text). that are over a factor of two too large at low values of n LM. Once again, this effect originates from the signal due to the rotational perturbation. Mitrovica and Milne [1998] have shown that TPW speeds for models with relatively weak lower mantle viscosity are significantly overestimated when an assumption of elastic incompressibility is introduced. This assumption affects the prediction of many of the other Stokes coefficients, though generally less than 10%. Indeed, the errors observed in the bottom three frames of Figure 2 are representative of the error introduced in many of the coefficients Earth and Ice Model Variations [16] Next, we examine the sensitivity of the predicted Stokes coefficients to variations in the Earth and ice model adopted in the GIA calculation. Previous studies have found that low degree zonal harmonics are generally insensitive to variations in lithospheric thickness [e.g., Mitrovica and

6 ETG 18-6 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS Figure 4. (continued) Peltier, 1993]. Figure 3 shows the prediction of selected Stokes coefficients as a function of n LM for three different values of h LT. Once again, _C 21 and _S 21 present a special case where the sensitivity evident on the figure originates from rotational effects. Predictions of GIA-induced TPW are sensitive to variations in the lithospheric thickness [e.g., Yuen et al., 1983; Peltier and Wu, 1983; Wu and Peltier, 1984; Mitrovica and Milne, 1998], particularly as the lower mantle viscosity of the Earth model is reduced [Mitrovica and Milne, 1998]. Specifically, thicker values of h LT increase the magnitude of the predicted TPW speed. For the other coefficients, the sensitivity to variations in h LT increases with an increase in lower mantle viscosity. As the degree of the Stokes coefficient increases, this trend becomes evident at lower values of n LM. These trends are straightforward to explain. As the degree of the coefficient or the lower mantle viscosity is increased, GIA-induced deformations become confined to progressively shallower depths. [17] In Figure 4 we provide contour plots showing predictions of Stokes coefficients for degrees two to five for a suite of Earth models in which both n UM and n LM are varied by nearly 2 orders of magnitude. For an upper mantle viscosity of Pa s and for harmonics other than _C 21 and _S 21, the predicted signal peaks for an intermediate value of n LM and decreases relatively smoothly for smaller and greater values of the lower mantle viscosity. This simply confirms a general form evident in many prior studies. As the degree increases, the predicted peak amplitude is obtained for progressively lower values of n LM and the sensitivity to

7 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS ETG 18-7 n UM tends to increase. The sensitivity to upper mantle viscosity also increases as n LM is increased and n UM is weakened (and thus as flow is more confined to the shallow mantle). As an example, when n UM < Pa s, the sensitivity of many coefficients (e.g., _C 51 ) to variations in lower mantle viscosity disappears for n LM >10 22 Pa s, and in these cases the double root structure of the n LM sensitivity evident for the n UM =10 21 Pa s case also disappears. [18] Some consensus appears to be emerging in recent inferences of the average upper and lower mantle viscosity based upon GIA data [e.g., Mitrovica, 1996; Lambeck et al., 1998]. These studies suggest a value of n UM somewhat less than Pa s, but greater than Pa s, and n LM between Pa s and Pa s. We overlay, on Figure 4, a shading that reflects this range of current uncertainty. The relatively strong sensitivity to variations in n LM below Pa s, evident for nearly all of the Stokes coefficients, is clearly avoided by limiting the viscosity model space in this manner. For most of the lowdegree harmonics this sensitivity is due to the sharp reduction in the present-day disequilibrium predicted for models in which n LM is reduced toward Pa s. However, in the case of _C 21 and _S 21, the predictions increase as n LM is reduced below Pa s as a consequence of a large increase in the predicted TPW signal. The highest amplitudes in Figure 4 are achieved for the zonal Stokes coefficients. However, we note that a set of the nonzonal harmonics reach relatively large values (e.g., _S 31, _S 41, etc.) and these may be strong candidates for detailed consideration using satellitederived constraints. [19] Uncertainties in the adopted late Pleistocene ice model also have potentially large effects on the predictions (for a detailed discussion, see Ivins et al. [1993]). In Figure 5 we show a set of predictions generated using both the ICE- 3G model and an ice model (ICE1B) modified from ICE1 [Peltier and Andrews, 1976] by the addition of an Antarctica loading cycle. The latter contribution was originally derived from Hughes et al. [1981] by Wu and Peltier [1983]. As we discussed in section 1, the large variation in the solution for _C 30 (and other odd zonal harmonics) is due to the sensitivity of this coefficient to the relative balance between the ice loads in the Northern Hemisphere and the Southern Hemisphere [Mitrovica and Peltier, 1993; Ivins et al., 1993]. Figure 5 indicates that predictions of the nonzonal harmonics are also sensitive to plausible variations in the ice model. In this regard, the differences apparent in Figure 5 are characteristic of variations we predicted for a much broader suite of Stokes coefficients. [20] In addition to providing a sensitivity analysis, Figures 3 5 also provide bounds on the GIA-induced contribution to low degree Stokes coefficients governing the secular variation of the geoid. These bounds will be useful for rigorous analysis of constraints emerging from the new set of satellite missions. These missions will constrain not only the Stokes coefficients but also regional variations of the geoid change, and so we turn next to predictions of the latter. [21] We computed global maps of the GIA-induced geoid change for a large suite of Earth models in which all combinations of the following values were adopted: h LT = 71, 96, or 120 km; n UM =1,2,3,5,8,or Pa s; and n LM = 2, 3, 5, 8, 10, 20, or Pa s. In addition, Figure 5. Differences in various Stokes coefficients predicted using two cycles of ICE-3G and ICE1B (see text). Results are generated using the Earth models described in Figure 1 (n UM = Pa s, h LT = 71 km). for each of these Earth models, we computed results for both the ICE-3G and ICE1B ice histories. To proceed, we then culled this set of 252 predictions to include only those which satisfy a recent constraint on the present-day uplift rate in Hudson Bay (10.7 ± 1.4 mm yr 1 [Larson and van Dam, 2000]). This procedure yielded a set of 66 predictions, the simple mean of which is shown in Figure 6a. The geographic variation along 90 W longitude for the individual predictions is given in Figure 6b. The latter profile passes through the center of the Laurentide deglaciation,

8 ETG 18-8 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS Figure 6. Regional variations of secular geoid rates computed from a large suite of Earth and ice models. (a) Average secular geoid change, through degree and order 64, of predictions from 66 combinations of h LT, n UM, n LM, and ice model (see text). (b) Profiles of the predictions for each of the 66 models along 90 W longitude, as indicated by the arrows in Figure 6a. Predictions generated using ICE- 3G are shown in blue, while those using ICE1B are shown in green. The red line is the average value of the profiles, as shown in Figure 6a. and we adopt different colors to distinguish predictions based on the two ice histories. [22] The global map is characterized by regions of upwarping over the main centers of late Pleistocene deglaciation, including Laurentia, Antarctica and Fennoscandia. In the first two of these regions the average predicted signal exceeds 1 mm yr 1. In addition, there is a distinct Y 21 pattern associated with rotational (TPW) effects that is evident in the far-field geoid signal. In the absence of the TPW signal there would be a broad subsidence of the geoid in the far-field. The rotation signal adds to this subsidence in the vicinity of Central America and southern Australia and the net effect is an average peak subsidence that exceeds 0.5 mm yr 1. Clearly, the GIA-induced geoid rate

9 TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS ETG 18-9 is nonnegligible even in the far field of the late Pleistocene ice cover [Mitrovica and Peltier, 1991; Peltier, 1999; Mitrovica et al., 2001]. [23] Predictions of the geoid rate over Laurentia and Antarctica are clearly sensitive to variations in the adopted ice and Earth models (Figure 6b). However, this sensitivity remains significant in the far field, where predictions have a range of 0.5 mm yr 1 (or 0.25 mm yr 1 around the average). When regional maps of geoid rates become available, efforts will be made to correct these maps using some standard model for GIA. The results in Figure 6b indicate that this correction may be subject to a relatively large uncertainty, even within the far-field. 3. Summary [24] Previous studies have shown that the predicted zonal harmonics of the secular geoid change due to GIA are sensitive to variations in the adopted Earth and ice model. We have extended this conclusion to nonzonal harmonics and have provided a detailed sensitivity analysis involving variations in both the main input parameters (mantle viscosity, lithospheric thickness, ice history) and elements of the forward theory (Earth rotation, elastic compressibility). Some Stokes coefficients show sensitivities akin to those identified in earlier zonal analyses, while others show distinct sensitivities originating, for example, from rotational effects (e.g., _C 21 and _S 21 ) and ocean loading (e.g., _C 33 ). Regardless of the source of the sensitivities, our analysis provides bounds on the variation of these coefficients which will be important as new values are derived from upcoming satellite missions. [25] The same reasoning motivated our detailed analysis of regional predictions of the secular geoid change. We find that these predictions are also highly sensitive to variations in the Earth and ice model. Indeed, even in the far field of the late Pleistocene ice sheets this variation is 0.5 mm yr 1 or greater. Because long time series will be required to constrain GIA from (for example) GRACE data, it will be difficult to robustly infer a GIA Earth/ice model combination from the data before the end of the mission. Thus it is more likely, at least in the short term, that the analysis of this data will proceed by correcting the data using some independently derived standard model. In this case, the ranges we have established may serve as a measure of the uncertainty in this correction. [26] There are at least two areas of this research that warrant further study. The first of these concerns predictions that involve a nonnegligible contribution from TPW effects (in particular, the coefficients _C 21 and _S 21, and the Y 21 imprinting evident on maps such as Figure 6a). Recent GIA analyses indicate that predictions of polar wander driven by surface mass loading are sensitive to the treatment of the density discontinuity at 670 km depth (nonadiabatic, as in this study, versus adiabatic (e.g., Mitrovica and Milne [1998] and, especially, Johnston and Lambeck [1999, Figure B1]) and the adoption of realistic viscoelastic stratification within the lithosphere [Nakada, 2000]. In regard to the latter, the validity (and accuracy) of our assumption of a global elastic lithosphere for TPW predictions on GIA timescales is uncertain. [27] The second concern involves the adoption of spherically symmetric Earth models. Upper mantle (including lithospheric) asymmetry between ocean and continent environments, as well as thermochemical variations driven by mantle convection, imply significant lateral variations in mantle viscosity. In a future paper we will examine the sensitivity of the predictions to these lateral variations. The calculations will be based on both the forward theory developed by Tromp and Mitrovica [1999] and finite element numerical simulations. Furthermore, the geometry of the variations will be prescribed using recent models from seismic tomography. [28] Acknowledgments. We thank Erik Ivins and Tom James for their thorough reviews of the original manuscript. We acknowledge the support of the Natural Sciences and Engineering Research Council, the Government of Ontario (via a PREA Award) and the Canadian Institute for Advanced Research. References Cazenave, A. A., P. Gegout, G. Ferhat, and R. Biancale, Temporal variations of the gravity field from LAGEOS 1 and LAGEOS 2 observations, in Global Gravity Field and Its Temporal Variations, edited by R. H. Rapp, A. A. Cazenave, and R. S. Nerem, pp , Springer-Verlag, New York, Chao, B. F., and R. S. Gross, Changes in the Earth s rotation and lowdegree gravitational field induced by earthquakes, Geophys. J. R. Astron. Soc., 91, , Cheng, M. K., C. K. Shum, and B. D. Tapley, Determination of long-term changes in the Earth s gravity field from satellite laser ranging observations, J. Geophys. Res., 102, 22,377 22,390, Devoti, R., V. Luceri, C. Sciarretta, G. Bianco, G. Di Donato, L. L. A. Vermeersen, and R. Sabadini, The SLR secular gravity variations and their impact on the inference of mantle rheology and lithospheric thickness, Geophys. Res. Lett., 28, , Dickey, J. O., et al., Satellite Gravity and the Geosphere: Contributions to the Study of the Solid Earth and Its Fluid Envelope, 112 pp., Natl. Acad. Press, Washington, D. C., Dziewonski, A. M., and D. L. Anderson, Preliminary reference Earth model (PREM), Phys. Earth Planet. Inter., 25, , Eanes, R. J., and S. V. Bettadpur, Temporal variability of Earth s gravitational field from satellite laser ranging, in Global Gravity Field and Its Temporal Variations, edited by R. H. Rapp, A. A. Cazenave, and R. S. Nerem, pp , Springer-Verlag, New York, Han, D., and J. Wahr, Post-glacial rebound analysis for a rotating Earth, in Slow Deformations and Transmission of Stress in the Earth, Geophys. Monogr. Ser., vol. 49, edited by S. Cohen and P. Vanicek, pp. 1 6, AGU, Washington, D. C., Hughes, T. J., G. H. Denton, B. G. Andersen, D. H. Schilling, J. L. Fastook, and C. S. Lingle, The last great ice sheets: a global view, in The Last Great Ice Sheets, edited by G. H. Denton and T. J. Hughes, pp , John Wiley, New York, Ivins, E. R., C. G. Sammis, and C. F. Yoder, Deep mantle viscous structure with prior estimate and satellite constraint, J. Geophys. Res., 98, , James, T. S., and E. R. Ivins, Global geodetic signatures of the Antarctic ice sheet, J. Geophys. Res., 102, , Johnston, P., and K. Lambeck, Postglacial rebound and sea level contributions to changes in the geoid and the Earth s rotation axis, Geophys. J. Int., 136, , Klosko, S., and B. Chao, Secular variations of the zonal gravity field, global sea level, and polar motion as geophysical constraints, Phys. Chem. Earth, 23, , Lambeck, K., Geophysical Geodesy, Oxford Univ. Press, New York, Lambeck, K., C. Smither, and P. Johnston, Sea-level change, glacial rebound and mantle viscosity for northern Europe, Geophys. J. Int., 134, , Larson, K. M., and T. van Dam, Measuring postglacial rebound with GPS and absolute gravity, Geophys. Res. Lett., 27, , Milne, G. A., and J. X. Mitrovica, Postglacial sea-level change on a rotating Earth: First results from a gravitationally self-consistent sea-level equation, Geophys. J. Int., 126, F13 F20, Milne, G. A., and J. X. Mitrovica, Postglacial sea-level change on a rotating Earth, Geophys. J. Int., 133, 1 10, 1998.

10 ETG TAMISIEA ET AL.: SENSITIVITY OF SECULAR GEOID PREDICTIONS Milne, G. A., J. X. Mitrovica, and J. L. Davis, Near-field hydro-isostasy: The implementation of a revised sea-level equation, Geophys. J. Int., 139, , Mitrovica, J. X., Haskell [1935] Revisited, J. Geophys. Res., 101, , Mitrovica, J. X., and G. A. Milne, Glaciation-induced perturbations in the Earth s rotation: A new appraisal, J. Geophys. Res., 103, , Mitrovica, J. X., and W. R. Peltier, Pleistocene deglaciation and the global gravity field, J. Geophys. Res., 94, 13,651 13,671, Mitrovica, J. X., and W. R. Peltier, On postglacial geoid subsidence over the equatorial oceans, J. Geophys. Res., 96, 20,053 20,071, Mitrovica, J. X., and W. R. Peltier, Present-day secular variations in the zonal harmonics of the Earth s geopotential, J. Geophys. Res., 98, , Mitrovica, J. X., G. A. Milne, and J. L. Davis, Glacial isostatic adjustment of a rotating Earth, Geophys. J. Int., 147, , Nakada, M., Effects of the viscoelastic lithosphere on polar wander speed caused by the late Pleistocene glacial cycles, Geophys. J. Int., 143, , Nerem, R. S., and S. M. Klosko, Secular variations of the zonal harmonics and polar motion as geophysical constraints, in Global Gravity Field and Its Temporal Variations, edited by R. H. Rapp, A. A. Cazenave, and R. S. Nerem, pp , Springer-Verlag, New York, Peltier, W. R., Constraint on deep mantle viscosity from LAGEOS acceleration data, Nature, 304, , Peltier, W. R., Postglacial variations in the level of the sea: Implications for climate dynamics and solid-earth geophysics, Rev. Geophys., 36, , Peltier, W. R., Global sea level rise and glacial isostatic adjustment, Global Planet. Change, 20, , Peltier, W. R., and J. T. Andrews, Glacial-isostatic adjustment, I, The forward problem, Geophys. J. R. Astron. Soc., 46, , Peltier, W. R., and P. Wu, Continental lithospheric thickness and deglaciation induced true polar wander, Geophys. Res. Lett., 10, , Sabadini, R., D. A. Yuen, and P. Gasperini, Mantle rheology and satellite signatures from present-day glacial forcings, J. Geophys. Res., 93, , Spada, G., R. Sabadini, D. A. Yuen, and Y. Ricard, Effects on post-glacial rebound from the hard rheology in the transition zone, Geophys. J. Int., 109, , Tromp, J., and J. X. Mitrovica, Surface loading of a viscoelastic Earth, I, General theory, Geophys. J. Int., 137, , Tushingham, A. M., and W. R. Peltier, ICE-3G: A new global model of late Pleistocene deglaciation based upon geophysical predictions of postglacial relative sea level change, J. Geophys. Res., 96, , Vermeersen, L. L. A., and R. Sabadini, Significance of the fundamental mantle relaxation mode in polar wander simulations, Geophys. J. Int., 127, F5 F9, Vermeersen, L. L. A., R. Sabadini, R. Devoti, V. Luceri, P. Rutigliano, C. Sciarretta, and G. Bianco, Mantle viscosity from joint inversions of Pleistocene deglaciation-induced changes in geopotential with a new SLR analysis and polar wander, Geophys. Res. Lett., 25, , Wu, P., and W. R. Peltier, Glacial isostatic adjustment and the free air gravity anomaly as a constraint on deep mantle viscosity, Geophys. J. R. Astron. Soc., 74, , Wu, P., and W. R. Peltier, Pleistocene deglaciation and the Earth s rotation: A new analysis, Geophys. J. R. Astron. Soc., 76, , Yoder, C. F., and E. R. Ivins, Changes in the Earth s gravity field from Pleistocene deglaciation and present-day glacial melting, Eos Trans. AGU, 66(18), 245, Yoder, C., J. Williams, J. Dickey, B. Schutz, R. Eanes, and B. Tapley, Secular variations of the earth s gravitational harmonic J 2 coefficient from LAGEOS and nontidal acceleration of earth rotation, Nature, 303, , Yuen, D. A., and R. Sabadini, Viscosity stratification of the lower mantle as inferred from the _J 2 observation, Ann. Geophys., 3, , Yuen, D. A., R. Sabadini, and E. Boschi, The dynamical equations of polar wander and the global characteristics of the lithosphere as extracted from rotational data, Phys. Earth Planet. Inter., 33, , Yuen, D. A., R. Sabadini, P. Gasperini, and E. Boschi, On transient rheology and glacial isostasy, J. Geophys. Res., 91, 11,420 11,438, Yuen, D. A., P. Gasperini, R. Sabadini, and E. Boschi, Azimuthal dependence in the gravity field induced by recent and past cryospheric forcings, Geophys. Res. Lett., 8, , G. A. Milne, Department of Geological Sciences, University of Durham, Science Laboratories, South Road, Durham, DH1 3LE, UK. (g.a.milne@ durham.ac.uk) J. X. Mitrovica and M. E. Tamisiea, Department of Physics, University of Toronto, 60 St. George St., Toronto, Ontario, Canada M5S 1A7. ( jxm@physics.utoronto.ca; tamisiea@physics.utoronto.ca) J. Tromp, California Institute of Technology, Seismological Laboratory, Mail Code , 1200 E. California Blvd., Pasadena, CA 91125, USA. ( jtromp@gps.caltech.edu)

Glacial isostatic adjustment on a rotating earth

Glacial isostatic adjustment on a rotating earth Geophys. J. Int. (2001) 147, 562 578 Glacial isostatic adjustment on a rotating earth Jerry X. Mitrovica, 1 Glenn A. Milne, 2 and James L. Davis 3 1 Department of Physics, University of Toronto, 60 St.

More information

High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation

High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation L.L.A. Vermeersen (1), H. Schotman (1), M.-W. Jansen (1), R. Riva (1) and R. Sabadini (2) (1) DEOS,

More information

The rotational stability of a triaxial ice age Earth

The rotational stability of a triaxial ice age Earth Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006564, 2010 The rotational stability of a triaxial ice age Earth I. Matsuyama, 1 J. X. Mitrovica, 2 A. Daradich,

More information

is the coefficient of degree 2, order 0 of the non-dimensional spherical harmonic

is the coefficient of degree 2, order 0 of the non-dimensional spherical harmonic Materials and Methods J is the coefficient of degree, order 0 of the non-dimensional spherical harmonic representation of the mass distribution of the Earth system. It is directly related to the diagonal

More information

The rotational stability of an ice-age earth

The rotational stability of an ice-age earth Geophys. J. Int. (005) 161, 491 506 doi: 10.1111/j.1365-46X5.0609.x he rotational stability of an ice-age earth Jerry X. Mitrovica, 1 John Wahr, Isamu Matsuyama 3 and Archie Paulson 1 Department of Physics,

More information

Significance of the fundamental mantle rotational relaxation mode in polar wander simulations

Significance of the fundamental mantle rotational relaxation mode in polar wander simulations Geophys. J. Int. (1996) 127, F5-F9 FAST-TRACK PAPER Significance of the fundamental mantle rotational relaxation mode in polar wander simulations L. L. A. Vermeersen* and R. Sabadini Dipartimento di Scienze

More information

Relationship between glacial isostatic adjustment and gravity perturbations observed by GRACE

Relationship between glacial isostatic adjustment and gravity perturbations observed by GRACE GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl048624, 2011 Relationship between glacial isostatic adjustment and gravity perturbations observed by GRACE A. Purcell, 1 A. Dehecq, 1,2 P. Tregoning,

More information

ISSN Scientific Technical Report STR 07/05. DOI: /GFZ.b GeoForschungsZentrum Potsdam

ISSN Scientific Technical Report STR 07/05. DOI: /GFZ.b GeoForschungsZentrum Potsdam ISSN 1610-0956 Ingo Sasgen, Robert Mulvaney, Volker Klemann and Detlef Wolf Glacial-isostatic adjustment and sea-level change near Berkner Island, Antarctica Ingo Sasgen Dep. 1: Geodesy and Remote Sensing

More information

Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth?

Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth? Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth? BARRETO L., FORTIN M.-A., IREDALE A. In this simple analysis, we compare the historical record

More information

Constraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data

Constraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data Constraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data Hugo Schotman 1,2, Bert Vermeersen 2, Pieter Visser 2 1 2 3 rd International GOCE User Workshop, ESA Esrin, 7 November 2006 glacial-isostatic

More information

Toward global inverse solutions for current and past ice mass variations: Contribution of secular satellite gravity and topography change measurements

Toward global inverse solutions for current and past ice mass variations: Contribution of secular satellite gravity and topography change measurements JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B11, 2291, doi:10.1029/2001jb000543, 2002 Toward global inverse solutions for current and past ice mass variations: Contribution of secular satellite gravity

More information

Space-Geodetic Constraints on Glacial Isostatic Adjustment in Fennoscandia

Space-Geodetic Constraints on Glacial Isostatic Adjustment in Fennoscandia Space-Geodetic Constraints on Glacial Isostatic Adjustment in Fennoscandia G. A. Milne, 1 * J. L. Davis, 2 Jerry X. Mitrovica, 3 H.-G. Scherneck, 4 J. M. Johansson, 4 M. Vermeer, 5 H. Koivula 5 Analysis

More information

Isostasy. Introduction

Isostasy. Introduction SEA LEVEL STUDIES/Isostasy 3043 Fleming, K., Johnston, P., Zwartz, D., Yokoyama, Y., Lambeck, K., and Chappell, J. (1998). Refining the eustatic sea-level curve since the Last Glacial Maximum using far-

More information

Inference of mantle viscosity from GRACE and relative sea level data

Inference of mantle viscosity from GRACE and relative sea level data Geophys. J. Int. (7) 7, 97 58 doi:./j.365-6x.7.3556.x FAST TRACK PAPER Inference of mantle viscosity from GRACE and relative sea level data Archie Paulson, Shijie Zhong and John Wahr Department of Physics,

More information

Vertical Motion from Satellite Altimetry and Tide gauges

Vertical Motion from Satellite Altimetry and Tide gauges Vertical Motion from Satellite Altimetry and Tide gauges Alexander Braun, C.K. Shum and C.-Y. Kuo Byrd Polar Research Center and Laboratory for Space Geodesy and Remote Sensing, The Ohio State University,

More information

Anomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment

Anomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment ANNALS OF GEOPHYSICS, 57, 4, 2014, S0432; doi:10.4401/ag-6548 Anomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment Giorgio Spada 1,*, Marco Olivieri 2, Gaia Galassi

More information

Isolating the PGR signal in the GRACE data: impact on mass balance estimates in Antarctica and Greenland

Isolating the PGR signal in the GRACE data: impact on mass balance estimates in Antarctica and Greenland Geophys. J. Int. (28) 172, 18 3 doi: 1.1111/j.1365-246X.27.363.x GJI Geodesy, potential field and applied geophysics Isolating the PGR signal in the GRACE data: impact on mass balance estimates in Antarctica

More information

Constraints on lithosphere and mantle rheology from in-situ observations

Constraints on lithosphere and mantle rheology from in-situ observations Constraints on lithosphere and mantle rheology from in-situ observations Shijie Zhong Department of Physics University of Colorado at Boulder Collaborators: Archie Paulson and John Wahr on post-glacial

More information

Constraints on Mantle Structure from Surface Observables

Constraints on Mantle Structure from Surface Observables MYRES I: Heat, Helium & Whole Mantle Convection Constraints on Mantle Structure from Surface Observables Magali Billen University of California, Davis Department of Geology The Goal Use observations of

More information

Kollo, Karin; Spada, Giorgio; Vermeer, Martin Studying earth rheology using GNSS permanent stations and GIA modelling tools

Kollo, Karin; Spada, Giorgio; Vermeer, Martin Studying earth rheology using GNSS permanent stations and GIA modelling tools Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Kollo, Karin; Spada, Giorgio; Vermeer,

More information

Limitations on the inversion for mantle viscosity from postglacial rebound

Limitations on the inversion for mantle viscosity from postglacial rebound Geophys. J. Int. (27) 168, 119 129 doi: 1.1111/j.136-246X.26.3222.x Limitations on the inversion for mantle viscosity from postglacial rebound Archie Paulson, Shijie Zhong and John Wahr Department of Physics,

More information

Postseismic signature of the 2004 Sumatra earthquake on low-degree gravity harmonics

Postseismic signature of the 2004 Sumatra earthquake on low-degree gravity harmonics Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005296, 2008 Postseismic signature of the 2004 Sumatra earthquake on low-degree gravity harmonics V. Cannelli,

More information

How important are elastic deflections in the Fennoscandian postglacial uplift?

How important are elastic deflections in the Fennoscandian postglacial uplift? How important are elastic deflections in the Fennoscandian postglacial uplift? WILLY FJELDSKAAR Fjeldskaar, W. How important are elastic deftections in the Fennoscandian postglacial uplift? Norsk Geologisk

More information

Rheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals

Rheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals (part 1): Deformation mechanisms and flow rules of mantle minerals What is rheology? Rheology is the physical property that characterizes deformation behavior of a material (solid, fluid, etc) solid mechanics

More information

Three-dimensional finite-element modelling of Earth s viscoelastic deformation: effects of lateral variations in lithospheric thickness

Three-dimensional finite-element modelling of Earth s viscoelastic deformation: effects of lateral variations in lithospheric thickness Geophys. J. Int. (23) 1, 679 695 Three-dimensional finite-element modelling of Earth s viscoelastic deformation: effects of lateral variations in lithospheric thickness Shijie Zhong, Archie Paulson and

More information

Course Business. Today: isostasy and Earth rheology, paper discussion

Course Business. Today: isostasy and Earth rheology, paper discussion Course Business Today: isostasy and Earth rheology, paper discussion Next week: sea level and glacial isostatic adjustment Email did you get my email today? Class notes, website Your presentations: November

More information

GEOLOGY Vol. II - Mantle Dynamics and Plate Kinematics - Carlo Doglioni, Roberto Sabadini

GEOLOGY Vol. II - Mantle Dynamics and Plate Kinematics - Carlo Doglioni, Roberto Sabadini MANTLE DYNAMICS AND PLATE KINEMATICS Carlo Doglioni La Sapienza University, Rome, Italy Roberto Sabadini University of Milan, Italy Keywords: lithosphere, asthenosphere, mantle, core, transition zone,

More information

Low degree gravity changes from GRACE, Earth rotation, geophysical models, and satellite laser ranging

Low degree gravity changes from GRACE, Earth rotation, geophysical models, and satellite laser ranging JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005397, 2008 Low degree gravity changes from GRACE, Earth rotation, geophysical models, and satellite laser ranging J. L. Chen 1 and C. R.

More information

Estimating geocenter variations from a combination of GRACE and ocean model output

Estimating geocenter variations from a combination of GRACE and ocean model output JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 3,, doi:10.1029/2007jb005338, 2008 Estimating geocenter variations from a combination of GRACE and ocean model output Sean Swenson, 1 Don Chambers, 2 and John Wahr

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION Departures from eustasy in Pliocene sea-level records Supplementary Information Maureen E. Raymo 1 * #, Jerry X. Mitrovica 2#, Michael J. O'Leary 3, Robert M. DeConto 4 and Paul

More information

The Rotational and Gravitational Signature of Recent Great Earthquakes

The Rotational and Gravitational Signature of Recent Great Earthquakes The Rotational and Gravitational Signature of Recent Great Earthquakes Richard S. Gross Jet Propulsion Laboratory California Institute of Technology Pasadena, CA 91109 8099, USA 7th IVS General Meeting

More information

Geodetic Observing Systems: tools in observing the Glacial Isostatic Adjustment. Markku Poutanen Finnish Geodetic Institute

Geodetic Observing Systems: tools in observing the Glacial Isostatic Adjustment. Markku Poutanen Finnish Geodetic Institute Geodetic Observing Systems: tools in observing the Glacial Isostatic Adjustment Markku Poutanen Finnish Geodetic Institute Glacial Isostatic Adjustment Land uplift is just one consequence of the physical

More information

Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data

Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data Chapter 2 Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data W.Sun,G.Fu,andSh.Okubo Abstract Dislocation theories were developed conventionally for a deformed earth

More information

Seismotectonics of intraplate oceanic regions. Thermal model Strength envelopes Plate forces Seismicity distributions

Seismotectonics of intraplate oceanic regions. Thermal model Strength envelopes Plate forces Seismicity distributions Seismotectonics of intraplate oceanic regions Thermal model Strength envelopes Plate forces Seismicity distributions Cooling of oceanic lithosphere also increases rock strength and seismic velocity. Thus

More information

Edge Driven Convection and Iceland

Edge Driven Convection and Iceland Edge Driven Convection and Iceland Scott D. King Department of Earth and Atmospheric Sciences Purdue University, West Lafayette, Indiana One of the alternative hypotheses for hotspot volcanism is Edge-Driven

More information

Earth gravity field recovery using GPS, GLONASS, and SLR satellites

Earth gravity field recovery using GPS, GLONASS, and SLR satellites 13-01-08 Earth gravity field recovery using GPS, GLONASS, and SLR satellites Krzysztof Sośnica (1), Adrian Jäggi (1), Daniela Thaller (2), Ulrich Meyer (1), Christian Baumann (1), Rolf Dach (1), Gerhard

More information

Recent tectonic plate decelerations driven by mantle convection

Recent tectonic plate decelerations driven by mantle convection Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L23301, doi:10.1029/2009gl040224, 2009 Recent tectonic plate decelerations driven by mantle convection A. M. Forte, 1 R. Moucha, 1 D.

More information

GRACE impact in geodesy and geophysics. R. Biancale (GRGS-CNES Toulouse), M. Diament (IPG Paris)

GRACE impact in geodesy and geophysics. R. Biancale (GRGS-CNES Toulouse), M. Diament (IPG Paris) GRACE impact in geodesy and geophysics R. Biancale (GRGS-CNES Toulouse), M. Diament (IPG Paris) Improvement of gravity models Since 2002 the GRACE mission has changed some goals in geodesy. It has become

More information

Sea level change recent past, present, future

Sea level change recent past, present, future Sea level change recent past, present, future Anny Cazenave LEGOS-CNES Toulouse, France CCI_Colocation meeting, ESA/ESRIN Global mean sea level rise during the 20 th century (red : tide gauges; blue :

More information

Journal of Geodynamics

Journal of Geodynamics Journal of Geodynamics 72 (2013) 59 66 Contents lists available at ScienceDirect Journal of Geodynamics journal homepage: http://www.elsevier.com/locate/jog Estimate of glacial isostatic adjustment uplift

More information

True polar wander associated with continental drift on a hypothetical Earth

True polar wander associated with continental drift on a hypothetical Earth Earth Planets Space, 59, 513 522, 27 True polar wander associated with continental drift on a hypothetical Earth Masao Nakada Department of Earth and Planetary Sciences, Faculty of Science, Kyushu University,

More information

topography half is commonly represented in maps, and how it interacts with

topography half is commonly represented in maps, and how it interacts with Topography T. Perron 12.001 We ll spend a large fraction of the second half of the course discussing Earth s surface. Today we ll do two things: First, we ll discuss the ways topography is commonly represented

More information

Shape of thermal plumes in a compressible mantle with depth-dependent viscosity

Shape of thermal plumes in a compressible mantle with depth-dependent viscosity GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl050959, 2012 Shape of thermal plumes in a compressible mantle with depth-dependent viscosity Wei Leng 1 and Michael Gurnis 1 Received 18 January

More information

Spatial variability of sea level rise due to water impoundment behind dams

Spatial variability of sea level rise due to water impoundment behind dams Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl043462, 2010 Spatial variability of sea level rise due to water impoundment behind dams Julia W. Fiedler 1 and Clinton

More information

B6 Isostacy. B6.1 Airy and Pratt hypotheses. Geophysics 210 September 2008

B6 Isostacy. B6.1 Airy and Pratt hypotheses. Geophysics 210 September 2008 B6 Isostacy B6.1 Airy and Pratt hypotheses Himalayan peaks on the Tibet-Bhutan border In the 19 th century surveyors used plumblines and theodolites to map India. A plumb line was used when measuring the

More information

7. Post Glacial Rebound. Ge 163 4/16/14-

7. Post Glacial Rebound. Ge 163 4/16/14- 7. Post Glacial Rebound Ge 163 4/16/14- Outline Overview Order of magnitude estimate of mantle viscosity Essentials of fluid mechanics Viscosity Stokes Flow Biharmonic equation Half-space model Channel

More information

RELATIVE SEA LEVEL VARIATIONS CAUSED BY SUBDUCTION

RELATIVE SEA LEVEL VARIATIONS CAUSED BY SUBDUCTION RELATIVE SEA LEVEL VARIATIONS CAUSED BY SUBDUCTION CLAUDIA PIROMALLO 1, GIORGIO SPADA 2, ROBERTO SABADINI 3 and YANICK RICARD 4 1 Istituto Nazionale di Geofisica, 00143 Roma, Italy; 2 Dipartimento di Fisica

More information

Long-term polar motion excited by ocean thermal expansion

Long-term polar motion excited by ocean thermal expansion Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L17603, doi:10.1029/2009gl039692, 2009 Long-term polar motion excited by ocean thermal expansion Felix W. Landerer, 1,2 Johann H. Jungclaus,

More information

Postglacial Rebound Modeling with Power-Law Rheology

Postglacial Rebound Modeling with Power-Law Rheology Postglacial Rebound Modeling with Power-Law Rheology Patrick Wu Dept. of Geology & Geophysics, University of Calgary, Calgary, Alberta T2N-1N4, Canada, ppwu@acs.ucalgary.ca Keywords: rheology, power law

More information

Physics Graduate Theses & Dissertations

Physics Graduate Theses & Dissertations University of Colorado, Boulder CU Scholar Physics Graduate Theses & Dissertations Physics Spring 1-1-2013 Viscoelastic Response to Surface and Tidal Loading Applications to Glacial Isostatic Adjustment

More information

Gravitational constraints

Gravitational constraints Gravitational constraints Reading: Fowler p172 187 Gravity anomalies Free-air anomaly: g F = g g( λ ) + δg obs F Corrected for expected variations due to the spheroid elevation above the spheroid Bouguer

More information

Accuracy of the International Terrestrial Reference Frame origin and Earth expansion

Accuracy of the International Terrestrial Reference Frame origin and Earth expansion GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl047450, 2011 Accuracy of the International Terrestrial Reference Frame origin and Earth expansion X. Wu, 1 X. Collilieux, 2 Z. Altamimi, 2 B. L.

More information

Satellite Gravimetry and its Application to Glaciology by Anthony Arendt for the UAF Summer School in Glaciology, June post-glacial rebound

Satellite Gravimetry and its Application to Glaciology by Anthony Arendt for the UAF Summer School in Glaciology, June post-glacial rebound Satellite Gravimetry and its Application to Glaciology by Anthony Arendt for the UAF Summer School in Glaciology, June 2010 1 Overview The Earth is a dynamic system in which components of the core, surface

More information

Global Inverse for Surface Mass Variations, Geocenter Motion, and Earth Rheology

Global Inverse for Surface Mass Variations, Geocenter Motion, and Earth Rheology Global Inverse for Surface Mass Variations, Geocenter Motion, and Earth Rheology Xiaoping Frank Wu Objectives Seasonal to Interannual Mass Variation: Elastic Earth - Multi-satellite data combination to

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO1521 Possible links between long-term geomagnetic variations and whole-mantle convection processes Biggin, A.J., Steinberger, B., Aubert, J., Suttie, N., Holme, R., Torsvik, T.H., van der

More information

New satellite mission for improving the Terrestrial Reference Frame: means and impacts

New satellite mission for improving the Terrestrial Reference Frame: means and impacts Fourth Swarm science meeting and geodetic missions workshop ESA, 20-24 March 2017, Banff, Alberta, Canada New satellite mission for improving the Terrestrial Reference Frame: means and impacts Richard

More information

Geophysics Departmental Exam: 2004 Part 1

Geophysics Departmental Exam: 2004 Part 1 2004 Geophysics Departmental Exam: 2004 Part 1 This section is 90 minutes, closed book, and consists of questions designed to test your knowledge of facts and figures in the geosciences. The focus will

More information

Contributions of Geodesy to Oceanography

Contributions of Geodesy to Oceanography Contributions of Geodesy to Oceanography B. Tapley and J. Ries Center for Space Research, The University of Texas at Austin Dynamic Planet 2005 Cairns, Australia August 22-26, 2005 August 22-26, 2005 Dynam

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

Glacial isostatic adjustment and nonstationary signals observed by GRACE

Glacial isostatic adjustment and nonstationary signals observed by GRACE Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb006161, 2009 Glacial isostatic adjustment and nonstationary signals observed by GRACE Paul Tregoning, 1 Guillaume

More information

Dynamic Subsidence and Uplift of the Colorado Plateau. Supplementary Material

Dynamic Subsidence and Uplift of the Colorado Plateau. Supplementary Material GSA DATA REPOSITORY 2010177 Liu and Gurnis Dynamic Subsidence and Uplift of the Colorado Plateau Supplementary Material Lijun Liu and Michael Gurnis Seismological Laboratory California Institute of Technology

More information

Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE)

Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE) WATER RESOURCES RESEARCH, VOL. 39, NO. 8, 1223, doi:10.1029/2002wr001808, 2003 Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE)

More information

Hydrology signal from GRACE gravity data in the Nelson River basin, Canada: a comparison of two approaches

Hydrology signal from GRACE gravity data in the Nelson River basin, Canada: a comparison of two approaches Li et al. Earth, Planets and Space (2) 7: https://doi.org/./s2---x FULL PAPER Hydrology signal from GRACE gravity data in the Nelson River basin, Canada: a comparison of two approaches Tanghua Li, Patrick

More information

Lower GRACE estimates of Antarctic sea-level contribution

Lower GRACE estimates of Antarctic sea-level contribution Lower GRACE estimates of Antarctic sea-level contribution Matt King 1,2, Rory Bingham 1, Philip Moore 1, Pippa Whitehouse 3, Mike Bentley 3, Glenn Milne 4 1 Newcastle University, 2 University of Tasmania,

More information

Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula

Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula Geophys. J. Int. (2008) 174, 672 695 doi:.1111/j.1365-246x.2008.03847.x Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula D. Melini, 1 V.

More information

Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula

Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula Geophys. J. Int. (2008) 174, 672 695 doi:.1111/j.1365-246x.2008.03847.x Post-seismic rebound of a spherical Earth: new insights from the application of the Post Widder inversion formula D. Melini, 1 V.

More information

ARTICLE IN PRESS. Quaternary Science Reviews

ARTICLE IN PRESS. Quaternary Science Reviews Quaternary Science Reviews xxx (2009) 1 17 Contents lists available at ScienceDirect Quaternary Science Reviews journal homepage: www.elsevier.com/locate/quascirev Closure of the budget of global sea level

More information

Recent postglacial rebound, gravity change and mantle flow in Fennoscandia

Recent postglacial rebound, gravity change and mantle flow in Fennoscandia Geophys. J. int. (1996) 126,229-234 Recent postglacial rebound, gravity change and mantle flow in Fennoscandia Martin Ekmanl and Jaakko Makinen2 'National Land Survey, Division of Geodetic Research, S-801

More information

GEOID UNDULATION DIFFERENCES BETWEEN GEOPOTENTIAL. RICHARD H. RAPP and YAN MING WANG

GEOID UNDULATION DIFFERENCES BETWEEN GEOPOTENTIAL. RICHARD H. RAPP and YAN MING WANG GEOID UNDULATION DIFFERENCES BETWEEN GEOPOTENTIAL MODELS RICHARD H. RAPP and YAN MING WANG Department of Geodetic Science and Surveying, The Ohio State University, Columbus, Ohio, U.S.A. (Received 15 September,

More information

Observation of glacial isostatic adjustment in stable North America with GPS

Observation of glacial isostatic adjustment in stable North America with GPS Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L02306, doi:10.1029/2006gl027081, 2007 Observation of glacial isostatic adjustment in stable North America with GPS Giovanni F. Sella,

More information

GLOBAL GLACIAL ISOSTASY AND THE SURFACE

GLOBAL GLACIAL ISOSTASY AND THE SURFACE Annu. Rev. Earth Planet. Sci. 2004. 32:111 49 doi: 10.1146/annurev.earth.32.082503.144359 Copyright c 2004 by Annual Reviews. All rights reserved First published online as a Review in Advance on February

More information

Continental Hydrology, Rapid Climate Change, and the Intensity of the Atlantic MOC: Insights from Paleoclimatology

Continental Hydrology, Rapid Climate Change, and the Intensity of the Atlantic MOC: Insights from Paleoclimatology Continental Hydrology, Rapid Climate Change, and the Intensity of the Atlantic MOC: Insights from Paleoclimatology W.R. Peltier Department of Physics University of Toronto WOCE derived N-S salinity section

More information

Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates

Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L06619, doi:10.1029/2003gl019285, 2004 Impact of short period, non-tidal, temporal mass variability on GRACE gravity estimates P. F. Thompson, S. V. Bettadpur, and

More information

Wiener optimal combination and evaluation of the Gravity Recovery and Climate Experiment (GRACE) gravity fields over Antarctica

Wiener optimal combination and evaluation of the Gravity Recovery and Climate Experiment (GRACE) gravity fields over Antarctica JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004605, 2007 Wiener optimal combination and evaluation of the Gravity Recovery and Climate Experiment (GRACE) gravity fields over Antarctica

More information

Physics of the Earth

Physics of the Earth Physics of the Earth Fourth edition Frank D Stacey CSIRO Exploration and Mining, Brisbane, Australia Paul M Davis Department of Earth and Space Sciences, University of California, Los Angeles, USA CAMBRIDGE

More information

Originally published as:

Originally published as: Originally published as: Lorenzo Martín, F., Wang, R., Roth, F. (2002): The effect of input parameters on visco-elastic models of crustal deformation. - Física de la Tierra, 14, 33-54 The effect of input

More information

Constraints on upper mantle viscosity from the flow induced pressure gradient across the Australian continental keel

Constraints on upper mantle viscosity from the flow induced pressure gradient across the Australian continental keel Article Volume 11, Number 6 4 June 2010 Q06004, doi:10.1029/2010gc003038 ISSN: 1525 2027 Click Here for Full Article Constraints on upper mantle viscosity from the flow induced pressure gradient across

More information

Love Numbers and Gravimetric Factor for Diurnal Tides. Piravonu M. Mathews

Love Numbers and Gravimetric Factor for Diurnal Tides. Piravonu M. Mathews Journal of the Geodetic Society of Japan Vol. 46, No. 4, (2001), pp. 231-236 Love Numbers and Gravimetric Factor for Diurnal Tides Piravonu M. Mathews Department of Theoretical Physics, University of Madras

More information

Deformation cycles of great subduction earthquakes in a viscoelastic Earth

Deformation cycles of great subduction earthquakes in a viscoelastic Earth Deformation cycles of great subduction earthquakes in a viscoelastic Earth Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada School of Earth and Ocean Science, University of Victoria????

More information

Determination of long-term changes in the Earth's gravity field from satellite laser ranging observations

Determination of long-term changes in the Earth's gravity field from satellite laser ranging observations JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 102, NO. B10, PAGES 22,377-22,390, OCTOBER 10, 1997 Determination of long-term changes in the Earth's gravity field from satellite laser ranging observations M. K.

More information

The geoid constraint in global geodynamics: viscosity structure, mantle heterogeneity models and boundary conditions

The geoid constraint in global geodynamics: viscosity structure, mantle heterogeneity models and boundary conditions Geophys. J. Znt. (1997) 131, 1-8 The geoid constraint in global geodynamics: viscosity structure, mantle heterogeneity models and boundary conditions Catherine Thoraval and Mark A. Richards Department

More information

Effects of ice melting on GRACE observations of ocean mass trends

Effects of ice melting on GRACE observations of ocean mass trends Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L05610, doi:10.1029/2006gl029171, 2007 Effects of ice melting on GRACE observations of ocean mass trends Don P. Chambers, 1 Mark E. Tamisiea,

More information

Modification of the lithospheric stress field by lateral variations in plate-mantle coupling

Modification of the lithospheric stress field by lateral variations in plate-mantle coupling GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L22307, doi:10.1029/2009gl040484, 2009 Modification of the lithospheric stress field by lateral variations in plate-mantle coupling J. B. Naliboff, 1 C. P. Conrad,

More information

The influence of short wavelength variations in viscosity on subduction dynamics

The influence of short wavelength variations in viscosity on subduction dynamics 1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.

More information

Neogene Uplift of The Barents Sea

Neogene Uplift of The Barents Sea Neogene Uplift of The Barents Sea W. Fjeldskaar A. Amantov Tectonor/UiS, Stavanger, Norway FORCE seminar April 4, 2013 The project (2010-2012) Funding companies Flat Objective The objective of the work

More information

The relation between gravity rate of change and vertical displacement in previously glaciated areas.

The relation between gravity rate of change and vertical displacement in previously glaciated areas. The relation between gravity rate of change and vertical displacement in previously glaciated areas. Per-Anders Olsson 1,2, Glenn Milne 3, Hans-Georg Scherneck 2, Jonas Ågren 1 Submitted to Journal of

More information

Time-variable gravity from SLR and DORIS tracking

Time-variable gravity from SLR and DORIS tracking Time-variable gravity from SLR and DORIS tracking Frank G. Lemoine 1, Steven M. Klosko 2, Christopher M. Cox 3, Thomas J. Johnson 4 1. Planetary Geodynamics Laboratory, NASA Goddard Space Flight Center,

More information

Application of Satellite Laser Ranging for Long- Wavelength Gravity Field Determination

Application of Satellite Laser Ranging for Long- Wavelength Gravity Field Determination Application of Satellite Laser Ranging for Long- Wavelength Gravity Field Determination J. C. Ries Center for Space Research The University of Texas at Austin Low Degree Gravity Variations from SLR GRACE

More information

C3.4.1 Vertical (radial) variations in mantle structure

C3.4.1 Vertical (radial) variations in mantle structure C3.4 Mantle structure Mantle behaves as a solid on short time scales (seismic waves travel through it and this requires elastic behaviour). Over geological time scales the mantle behaves as a very viscous

More information

Modelling of the GIA-induced surface gravity change over Fennoscandia

Modelling of the GIA-induced surface gravity change over Fennoscandia Modelling of the GIA-induced surface gravity change over Fennoscandia Per-Anders Olsson 1, Jonas Ågren 1, Hans-Georg Scherneck 2 June 15, 212 1 Geodetic Research Division, Lantmäteriet (Swedish mapping,

More information

12 Gravity and Topography

12 Gravity and Topography 12. Gravity and Topography 126 12 Gravity and Topography In the limit of perfect hydrostatic equilibrium, there is an exact relationship between gravity and topography... and no new information emerges:

More information

Physical Dynamics of the Coastal Zone in the Mediterranean on Annual to Decadal Scales

Physical Dynamics of the Coastal Zone in the Mediterranean on Annual to Decadal Scales Physical Dynamics of the Coastal Zone in the Mediterranean on Annual to Decadal Scales Hans-Peter Plag Nevada Bureau of Mines and Geology and Seismological Laboratory University of Nevada, Reno, Nevada,

More information

Radial profile of mantle viscosity: Results from the joint inversion of convection and postglacial

Radial profile of mantle viscosity: Results from the joint inversion of convection and postglacial JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 102, NO. B2, PAGES 2751-2769, FEBRUARY 10, 1997 Radial profile of mantle viscosity: Results from the joint inversion of convection and postglacial rebound observables

More information

Mars without the equilibrium rotational figure, Tharsis, and the remnant rotational figure

Mars without the equilibrium rotational figure, Tharsis, and the remnant rotational figure JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.109/010je003686, 010 Mars without the equilibrium rotational figure, Tharsis, and the remnant rotational figure I. Matsuyama 1 and M. Manga 1 Received

More information

Copyright 2004 American Geophysical Union. Further reproduction or electronic distribution is not permitted.

Copyright 2004 American Geophysical Union. Further reproduction or electronic distribution is not permitted. Copyright 2004 American Geophysical Union. Further reproduction or electronic distribution is not permitted. Citation: Thompson, P. F., S. V. Bettadpur, and B. D. Tapley (2004), Impact of short period,

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION Online Supplementary Information SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1829 (Chen et al., Contribution of polar and mountain glacier melt to recent sea level rise) 1. Steric Sea Level Changes From

More information

Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle

Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2733 Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle Edward J. Garnero 1, Allen K. McNamara 1, and Sang-Heon D. Shim 1

More information

A sea-level test for inertial interchange true polar wander events

A sea-level test for inertial interchange true polar wander events Geophys. J. Int. (1999) 136, F5-F10 FAST-TRACK PAPER A sea-level test for inertial interchange true polar wander events Jon E. Mound, 1 Jerry X. Mitrovica, 1 David A. D. Evans 2,* and Joseph L. Kirschvink

More information

The Earth s time-variable gravity field observed by GOCE

The Earth s time-variable gravity field observed by GOCE The Earth s time-variable gravity field observed by GOCE GOCE+ Time-Variations, part of STSE (Theme 4, Innovative Feasibility Studies) J. Bouman, M. Fuchs, C. Haberkorn, V. Lieb, M. Schmidt T. Broerse,

More information

Special Bureau for Loading: Potential Role Versus Reality

Special Bureau for Loading: Potential Role Versus Reality Special Bureau for Loading: Potential Role Versus Reality Hans-Peter Plag Nevada Bureau of Mines and Geology and Seismological Laboratory University of Nevada, Reno Special Bureau for Loading: Potential

More information