ALBERT

Description: The geocentre motion is the motion of the centre of mass of the entire Earth, considered an isolated system, in a terrestrial system of reference. We first derive a formula relating the harmonic degree-1 Lagrangian variation of the gravity at a station to both the harmonic degree-1 vertical displacement of the station and the displacement of the whole Earth's centre of mass. The relationship is independent of the nature of the Earth deformation and is valid for any source of deformation. We impose no constraint on the system of reference, except that its origin must initially coincide with the centre of mass of the spherically symmetric Earth model. Next, we consider the geocentre motion caused by surface loading. In a system of reference whose origin is the centre of mass of the solid Earth, we obtain a specific relationship between the gravity variation at the surface, the geocentre displacement and the load Love number $h^{\prime }_1$ , which demands the Earth's structure and rheological behaviour be known. For various networks of real or fictitious stations, we invert synthetic signals of surface gravity variations caused by atmospheric loading to retrieve the degree-1 variation of gravity. We then select six well-distributed stations of the Global Geodynamics Project, which is a world network of superconducting gravimeters, to invert actual gravity data for the degree-1 variations and determine the geocentre displacement between the end of 2004 and the beginning of 2012, assuming it to be due to surface loading. We find annual and semi-annual displacements with amplitude 0.5–2.3 mm.

Description: We review the theory of the Earth's elastic and gravitational response to a surface disk load. The solutions for displacement of the surface and the geoid are developed using expansions of Legendre polynomials, their derivatives and the load Love numbers. We provide a matlab  function called diskload that computes the solutions for both uncompensated and compensated disk loads. In order to numerically implement the Legendre expansions, it is necessary to choose a harmonic degree, n max , at which to truncate the series used to construct the solutions. We present a rule of thumb (ROT) for choosing an appropriate value of n max , describe the consequences of truncating the expansions prematurely and provide a means to judiciously violate the ROT when that becomes a practical necessity.

Description: Two types of signals are clearly visible in continuous GPS (cGPS) time-series in Iceland, in particular in the vertical component. The first one is a yearly seasonal cycle, usually sinusoid-like with a minimum in the spring and a maximum in the fall. The second one is a trend of uplift, with higher values the closer the cGPS stations are to the centre of Iceland and ice caps. Here, we study the seasonal cycle signal by deriving its average at 71 GPS sites in Iceland. We estimate the annual and semi-annual components of the cycle in their horizontal and vertical components using a least-squares adjustment. The peak-to-peak amplitude of the cycle of the vertical component at the studied sites ranges from 4 mm near the coastline up to 27 mm at the centre of the Vatnajökull, the largest ice cap in Iceland. The minimum of the seasonal cycle occurs earlier in low lying areas than in the central part of Iceland, consistent with snow load having a large influence on seasonal deformation. Modelling shows that the seasonal cycle is well explained by accounting for elastically induced surface displacements due to snow, atmosphere, reservoir lake and ocean variations. Model displacement fields are derived considering surface loads on a multilayered isotropic spherical Earth. Through forward and inverse modelling, we were able to reproduce a priori information on the average seasonal cycle of known loads (atmosphere, snow in non-glaciated areas and lake reservoir) and get an estimation of other loads (glacier mass balance and ocean). The seasonal glacier mass balance cycle in glaciated areas and snow load in non-glaciated areas are the main contributions to the seasonal deformation. For these loads, induced seasonal vertical displacements range from a few millimetres far from the loads in Iceland, to more than 20 mm at their centres. Lake reservoir load also has to be taken into account on local scale as it can generate up to 20 mm of vertical deformation. Atmosphere load and ocean load are observable and generate vertical displacements in the order of a few millimetres. Inversion results also shows that the Iceland crust is less rigid than the world average. Interannual deviation from the GPS seasonal cycle can occur and are caused by unusual weather conditions over extended period of time.

Description: Seismic waves produced by fault ruptures give rise to gravity perturbations. So far, these perturbations have either been modelled as permanent coseismic gravity change in a half-space or spherical Earth model, or as full time-domain model in infinite space. In this paper, we present the explicit solution of gravity perturbations in time domain produced by a double-couple buried in a homogeneous half-space. This result is especially suited to study gravity perturbations up to a few hundreds of kilometres from the epicentre. It facilitates detailed parametric studies of gravity perturbations from fault rupture, and predicts gravity perturbations of real earthquakes with greatly improved accuracy. The results may serve to develop first designs of gravity-assisted earthquake early-warning systems, made possible by a new generation of ultrasensitive gravity gradiometers, which is currently under development.

Description: Continuous gravimetric observations have been made with three successive generations of superconducting gravimeter over 20 yr at Syowa Station ( $39.6\deg$ E, $69.0\deg$ S), East Antarctica. The third-generation instrument, OSG#058, was installed in January 2010 and was calibrated by an absolute gravimeter during January and February, 2010. The estimated scale factor was –73.823 ± 0.053 μGal V –1 (1 μGal = 10 –8 m s –2 ). The first 5 yr of OSG#058 data from 2010 January 7 to 2015 January 10 were decomposed into tidal waves (M3 to Ssa) and other non-tidal components by applying the Bayesian tidal analysis program BAYTAP. Long-term non-tidal gravity residuals, which were obtained by subtracting annual and 18.6 year tidal waves and the predicted gravity response to the Earth's variable rotation, showed significant correlation with the accumulated snow depth measured at Syowa Station. The greatest correlation occurred when the gravity variations lagged the accumulated snow depth by 21 d. To estimate the gravitational effect of the accumulated snow mass, we inferred a conversion factor of 3.13 ± 0.08 μGal m –1 from this relation. The accumulated snow depth at Syowa Station was found to represent an extensive terrestrial water storage (the snow accumulation) around Syowa Station, which was estimated from the Gravity Recovery and Climate Experiment satellite gravity data. The snow accumulation around Syowa Station was detectable by the superconducting gravimeter.

Description: The new release AIUB-RL02 of monthly gravity models from GRACE GPS and K-Band range-rate data is based on reprocessed satellite orbits referring to the reference frame IGb08. The release is consistent with the IERS2010 conventions. Improvements with respect to its predecessor AIUB-RL01 include the use of reprocessed (RL02) GRACE observations, new atmosphere and ocean dealiasing products (RL05), an upgraded ocean tide model (EOT11A), and the interpolation of shallow ocean tides (admittances). The stochastic parametrization of AIUB-RL02 was adapted to include daily accelerometer scale factors, which drastically reduces spurious signal at the 161 d period in C 20 and at other low degree and order gravity field coefficients. Moreover, the correlation between the noise in the monthly gravity models and solar activity is considerably reduced in the new release. The signal and the noise content of the new AIUB-RL02 monthly gravity fields are studied and calibrated errors are derived from their non-secular and non-seasonal variability. The short-period time-variable signal over the oceans, mostly representing noise, is reduced by 50 per cent with respect to AIUB-RL01. Compared to the official GFZ-RL05a and CSR-RL05 monthly models, the AIUB-RL02 stands out by its low noise at high degrees, a fact emerging from the estimation of seasonal variations for selected river basins and of mass trends in polar regions. Two versions of the monthly AIUB-RL02 gravity models, with spherical harmonics resolution of degree and order 60 and 90, respectively, are available for the time period from March 2003 to March 2014 at the International Center for Global Earth Models or from ftp://ftp.unibe.ch/aiub/GRAVITY/GRACE (last accessed 22 March 2016).

Description: We document two kinds of traveling ionospheric disturbances, namely, CTIDs (Co-tsunami-Traveling-Ionospheric-disturbances) and ATIDs (Ahead-of-Tsunami-Traveling-Ionospheric-disturbances) related to the Tohoku-Oki tsunami of 2011 March 11. They are referred to the disturbances that remain behind and ahead of the principal tsunami wave front, respectively. We first note their presence in a numerical experiment performed using a simulation code coupling the tsunami, atmosphere and ionosphere. This code uses the tsunami wavefield as an input and simulates acoustic-gravity waves (AGWs) in the atmosphere and TIDs, in the form of total electron content (TEC) disturbance, in the ionosphere. The simulated TEC reveals the excitation of CTIDs (at about 2 TECU) and ATIDs (at about 1 TECU), representing up to 5 per cent disturbance over the ambient electron density, and they arise from the dissipation of AGWs in the thermosphere. A novel outcome is that during the tsunami passage between ~6° and 12° of epicentral distance, strong ATIDs arrive ~20–60 min ahead of the tsunami wave front covering ~3°–10° of distance from the tsunami location. Simulation results are compared with the far-field observations using GNSS satellites and confirm that ATIDs are the first detected TEC maximum, occurring 20–60 min ahead of the tsunami arrival. Our simulation also confirms the presence of largest TEC maximum representing CTIDs, 10–20 min after the first tsunami wave. ATIDs reported in this study have characteristics that can be potentially used for the early warning of the tsunami.

Description: Apparent acceleration in Gravity Recovery and Climate Experiment (GRACE) Antarctic ice mass time-series may reflect both ice discharge and surface mass balance contributions. However, a recent study suggests there is also contamination from errors in atmospheric pressure de-aliasing fields [European Center for Medium-Range Weather Forecast (ECMWF) operational products] used during GRACE data processing. To further examine this question, we compare GRACE atmospheric pressure de-aliasing (GAA) fields with in situ surface pressure data from coastal and inland stations. Differences between the two are likely due to GAA errors, and provide a measure of error in GRACE solutions. Time-series of differences at individual weather stations are fit to four presumed error components: annual sinusoids, a linear trend, an acceleration term and jumps at times of known ECMWF model changes. Using data from inland stations, we estimate that atmospheric pressure error causes an acceleration error of about +7.0 Gt yr –2 , which is large relative to prior GRACE estimates of Antarctic ice mass acceleration in the range of –12 to –14 Gt yr –2 . We also estimate apparent acceleration rates from other barometric pressure (reanalysis) fields, including ERA-Interim, MERRA and NCEP/DOE. When integrated over East Antarctica, the four mass acceleration estimates (from GAA and the three reanalysis fields) vary considerably (by ~2–16 Gt yr –2 ). This shows the need for further effort to improve atmospheric mass estimates in this region of sparse in situ observations, in order to use GRACE observations to measure ice mass acceleration and related sea level change.

Description: The long-wavelength gravity field contains information about processes in the sublithospheric mantle. As satellite-derived gravity models now provide the long to medium-wavelength gravity field at unprecedented accuracy, techniques used to process gravity data need to be updated. We show that when determining these long-wavelengths, the treatment of topographic-isostatic effect (TIE) and isostatic effects (IE) is a likely source of error. We constructed a global isostatic model and calculated global TIE and IE. These calculations were done for ground stations as well as stations at satellite height. We considered both gravity and gravity gradients. Using these results, we determined how much of the gravity signal comes from distant sources. We find that a significant long-wavelength bias is introduced if far-field effects on the topographic effect are neglected. However, due to isostatic compensation far-field effects of the topographic effect are to a large degree compensated by the far-field IE. This means that far-field effects can be reduced effectively by always considering topographic masses together with their compensating isostatic masses. We show that to correctly represent the ultra-long wavelengths, a global background model should be used. This is demonstrated both globally and for a continental-scale case area in North America. In the case of regional modelling, where the ultra-long wavelengths are not of prime importance, gravity gradients can be used to help minimize correction errors caused by far-field effects.

Description: Constraining laterally varying structures in planetary interiors is important for understanding both the composition and the internal dynamics of a planet. Recognizing that seismic imaging technique is currently only viable for studying the Earth's interior structures, methods that can be supported by advanced space geodetic techniques may become alternatives to ‘image’ the interiors of other planets. The method of tidal tomography is one possibility, and it relies on high precision measurement of the response of a planet to its body tide. However, it is essential to develop an efficient analytical tool that computes the dependence of tidal response to 3-D interior structures. In this paper, we present a complete formulation of such an analytical tool, which calculates to high accuracy the tidal response of a terrestrial planet with lateral heterogeneities in its elastic and density structures. We treat the lateral heterogeneities as small perturbations and derive the governing equations based on the perturbation theory. In a spherical harmonic representation, equations at each order of perturbation are reduced into multiple matrix equations at harmonics that are allowed by mode couplings, and the total response equals the sum of all those single-harmonic responses, which can be solved semi-analytically. We test our perturbation method by applying it to the Moon with a harmonic degree-1 mantle structure for which the perturbation solutions of the tidal response are compared with those from a fully numerical method. The remarkable agreement between results from these two methods validates the perturbation method. As an example, we then use the perturbation method to evaluate the impact of lunar crustal thickness variations on tidal response of the Moon. We find that lunar crust produces much smaller degree-3 tidal responses than a relatively weak degree-1 structure in the deep lunar mantle. Our calculations show that degree-3 tidal response measurements may hold key constraints on possible degree-1 mantle structure of the Moon, as suggested from previous modelling results.