ALBERT

Description: Sensitive instruments like strainmeters and tiltmeters are necessary for measuring slowly varying low amplitude Earth deformations. Nonetheless, laser and fibre interferometers are particularly suitable for interrogating such instruments due to their extreme precision and accuracy. In this paper, a practical design of a simple pendulum borehole tiltmeter based on laser fibre interferometric displacement sensors is presented. A prototype instrument has been constructed using welded borosilicate with a pendulum length of 0.85 m resulting in a main resonance frequency of 0.6 Hz. By implementing three coplanar extrinsic fibre Fabry-Perot interferometric probes and appropriate signal filtering, our instrument provides tilt measurements that are insensitive to parasitic deformations caused by temperature and pressure variations. This prototype has been installed in an underground facility (Rustrel, France) where results show accurate measurements of Earth strains derived from Earth and ocean tides, local hydrologic effects, as well as local and remote earthquakes. The large dynamic range and the high sensitivity of this tiltmeter render it an invaluable tool for numerous geophysical applications such as transient fault motion, volcanic strain and reservoir monitoring.

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: In the present work we illustrate a new local inversion algorithm to retrieve the Moho depth from GOCE (Gravity field and steady-state Ocean Circulation Explorer) gravity field. In details the proposed procedure can be divided into two main steps: the first one consists in recognizing and isolating the different geological provinces in the study area by exploiting information coming from the GOCE global gravity field model. Once the main geological provinces are defined, a function relating the crust density of each province with depth is built and used to reduce the data. The gravitational effects of sediments, topography, bathymetry and upper mantle are also removed. In the second step the residual gravitational field is inverted to retrieve the Moho depth and some information on the crustal density. In particular, the clustering of geological province is performed by means of an automatic Bayesian classification algorithm while the inversion of GOCE residual field is performed by adapting the global algorithm developed in the framework of the GEMMA project to the local scale. The procedure, based on an iterative Wiener filter, allows to compute the Moho depth considering lateral as well as radial variations of crustal density. The algorithm has been applied to the fifth release of GOCE time-wise global gravity field model to infer information on the crustal structure in the Western Balkan area, that is, the region laying between Bulgaria and the Adriatic Sea. This region is one of the most complex and active, from the tectonic point of view, in the whole Europe and it is characterized by the presence of the Alpine-Himalayan orogenic belt, formed by the collision between the African and Eurasian plates, and by the opening of the Pannonian Basin. Results show a good agreement between the obtained geological provinces with the actual knowledge on the region. The resulting Moho depth ranges between about 20 km beneath the Adriatic Sea and 45 km in the Dinarides. Comparisons with available seismic data show differences smaller than 1 km (standard deviation).

Description: Normal mode treatments of the Earth's body tide response were developed in the 1980s to account for the effects of Earth rotation, ellipticity, anelasticity and resonant excitation within the diurnal band. Recent space-geodetic measurements of the Earth's crustal displacement in response to luni-solar tidal forcings have revealed geographical variations that are indicative of aspherical deep mantle structure, thus providing a novel data set for constraining deep mantle elastic and density structure. In light of this, we make use of advances in seismic free oscillation literature to develop a new, generalized normal mode theory for the tidal response within the semi-diurnal and long-period tidal band. Our theory involves a perturbation method that permits an efficient calculation of the impact of aspherical structure on the tidal response. In addition, we introduce a normal mode treatment of anelasticity that is distinct from both earlier work in body tides and the approach adopted in free oscillation seismology. We present several simple numerical applications of the new theory. First, we compute the tidal response of a spherically symmetric, non-rotating, elastic and isotropic Earth model and demonstrate that our predictions match those based on standard Love number theory. Second, we compute perturbations to this response associated with mantle anelasticity and demonstrate that the usual set of seismic modes adopted for this purpose must be augmented by a family of relaxation modes to accurately capture the full effect of anelasticity on the body tide response. Finally, we explore aspherical effects including rotation and we benchmark results from several illustrative case studies of aspherical Earth structure against independent finite-volume numerical calculations of the semi-diurnal body tide response. These tests confirm the accuracy of the normal mode methodology to at least the level of numerical error in the finite-volume predictions. They also demonstrate that full coupling of normal modes, rather than group coupling, is necessary for accurate predictions of the body tide response.

Description: We have extended backwards from 2001 to 1979 the current release 05 (RL05) of the Gravity Recovery and Climate Experiment (GRACE) Atmospheric and Oceanic De-aliasing Level-1B (AOD1B) product and studied the impact of this and a previous release 04 (RL04) of the AOD1B product on precise orbits of five altimetry satellites (ERS-1, ERS-2, TOPEX/Poseidon, Envisat and Jason-1) for the time span 1991–2012, as compared to the case when no AOD1B product is used. We have found that using AOD1B RL05 product reduces root mean square (RMS) fits of satellite laser ranging (SLR) observations by about 1.0–6.4 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by about 1.3–12.0, 0.3–10.0 and 2.0–10.0 per cent, respectively, for various satellites tested, as compared to the case without AOD1B product. Using AOD1B RL05 product instead of RL04 one reduces SLR RMS fits by 0.1–0.7 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by 0.1–0.6, 0.1–1.3 and 0.2–1.2 per cent, respectively, for the satellite orbits tested. The multi-mission crossover analysis shows that the application of an AOD1B product reduces the scatter of radial errors by 0.4–2.8 per cent for the satellite missions studied. At the regions with the most pronounced changes the use of the AOD1B products improves the consistency between the sea level as measured by the TOPEX and ERS-2 missions and by the Jason-1 and Envisat missions by 5 to 10 per cent (globally by about 2 per cent). The results of our study show that extended AOD1B RL05 product performs better than the AOD1B RL04 and improves orbits of altimetry satellites and consistency of sea level products.

Description: Extracting geophysical signals from Global Positioning System (GPS) coordinate time-series is a well-established practice that has led to great insights into how the Earth deforms. Often small discontinuities are found in such time-series and are traceable to either broad-scale deformation (i.e. earthquakes) or discontinuities due to equipment changes and/or failures. Estimating these offsets accurately enables the identification of coseismic deformation estimates in the former case, and the removal of unwanted signals in the latter case which then allows tectonic rates to be estimated more accurately. We develop a method to estimate accurately discontinuities in time series of GPS positions at specified epochs, based on a so-called ‘offset series’. The offset series are obtained by varying the amount of GPS data before and after an event while estimating the offset. Two methods, a mean and a weighted mean method, are then investigated to produce the estimated discontinuity from the offset series. The mean method estimates coseismic offsets without making assumptions about geophysical processes that may be present in the data (i.e. tectonic rate, seasonal variations), whereas the weighted mean method includes estimating coseismic offsets with a model of these processes. We investigate which approach is the most appropriate given certain lengths of available data and noise within the time-series themselves. For the Sumatra–Andaman event, with 4.5 yr of pre-event data, we show that between 2 and 3 yr of post-event data are required to produce accurate offset estimates with the weighted mean method. With less data, the mean method should be used, but the uncertainties of the estimated discontinuity are larger.

Description: Traditional processing of Global Navigation Satellite System (GNSS) data using dedicated scientific software has provided the highest levels of positional accuracy, and has been used extensively in geophysical deformation studies. To achieve these accuracies a significant level of understanding and training is required, limiting their availability to the general scientific community. Various online GNSS processing services, now freely available, address some of these difficulties and allow users to easily process their own GNSS data and potentially obtain high quality results. Previous research into these services has focused on Continually Operating Reference Station (CORS) GNSS data. Less research exists on the results achievable with these services using large campaign GNSS data sets, which are inherently noisier than CORS data. Even less research exists on the quality of velocity fields derived from campaign GNSS data processed through online precise point positioning services. Particularly, whether they are suitable for geodynamic and deformation studies where precise and reliable velocities are needed. In this research, we process a very large campaign GPS data set (spanning 10 yr) with the online Jet Propulsion Laboratory Automated Precise Positioning Service. This data set is taken from a GNSS network specifically designed and surveyed to measure deformation through the central North Island of New Zealand. This includes regional CORS stations. We then use these coordinates to derive a horizontal and vertical velocity field. This is the first time that a large campaign GPS data set has been processed solely using an online service and the solutions used to determine a horizontal and vertical velocity field. We compared this velocity field to that of another well utilized GNSS scientific software package. The results show a good agreement between the CORS positions and campaign station velocities obtained from the two approaches. We discuss the implications of these results for how future GNSS campaign field surveys might be conducted and how their data might be processed.

Description: In a pioneering study, Wahr & Bergen developed the widely adopted, pseudo-normal mode framework for predicting the impact of anelastic effects on the Earth's body tides. Lau et al. have recently derived an extended normal mode treatment of the problem (as well as a minor variant of the theory known as the direct solution method) that makes full use of theoretical developments in free oscillation seismology spanning the last quarter century and that avoids a series of assumptions and approximations adopted in the traditional theory for predicting anelastic effects. There are two noteworthy differences between these two theories: (1) the traditional theory only considers perturbations to the eigenmodes of an elastic Earth, whereas the new theory augments this set of modes to include the relaxation modes that arise in anelastic behaviour; and (2) the traditional theory approximates the complex perturbation to the tidal Love number as a scaled version of the complex perturbation to the elastic moduli, whereas the new theory computes the full complex perturbation to each eigenmode. In this study, we highlight the above differences using a series of synthetic calculations, and demonstrate that the traditional theory can introduce significant error in predictions of the complex perturbation to the Love numbers due to anelasticity and the related predictions of tidal lag angles. For the simplified Earth models we adopt, the computed lag angles differ by ~20 per cent. The assumptions in the traditional theory have important implications for previous studies that use model predictions to correct observables for body tide signals or that analyse observations of body tide deformation to infer mantle anelastic structure. Finally, we also highlight the fundamental difference between apparent attenuation (i.e. attenuation inferred from observations or predicted using the above theories) and intrinsic attenuation (i.e. the material property investigated through experiments), where both are often expressed in terms of lag angles or Q –1 . In particular, we demonstrate the potentially significant (factor of two or more) bias introduced in estimates of Q –1 and its frequency dependence in studies that have treated Q –1 determined from tidal phase lags or measured experimentally as being equal. The observed or theoretically predicted lag angle (or apparent Q –1 ) differs from the intrinsic, material property due to inertia, self-gravity and effects associated with the energy budget. By accounting for these differences we derive, for a special case, an expression that accurately maps apparent attenuation predicted using the extended normal mode formalism of Lau et al. into intrinsic attenuation. The theory allows for more generalized mappings which may be used to robustly connect observations and predictions of tidal lag angles to results from laboratory experiments of mantle materials.

Description: We consider a new approach to both the forward and inverse problems in post-seismic deformation. We present a method for forward modelling post-seismic deformation in a self-gravitating, heterogeneous and compressible earth with a variety of linear and nonlinear rheologies. We further demonstrate how the adjoint method can be applied to the inverse problem both to invert for rheological structure and to calculate the sensitivity of a given surface measurement to changes in rheology or time-dependence of the source. Both the forward and inverse aspects are illustrated with several numerical examples implemented in a spherically symmetric earth model.

Description: Moment accumulation rate (also referred to as moment deficit rate) is a fundamental quantity for evaluating seismic hazard. The conventional approach for evaluating moment accumulation rate of creeping faults is to invert for the slip distribution from geodetic measurements, although even with perfect data these slip-rate inversions are non-unique. In this study, we show that the slip-rate versus depth inversion is not needed because moment accumulation rate can be estimated directly from surface geodetic data. We propose an integral approach that uses dense geodetic observations from Interferometric Synthetic Aperture Radar (InSAR) and the Global Positioning System (GPS) to constrain the moment accumulation rate. The moment accumulation rate is related to the integral of the product of the along-strike velocity and the distance from the fault. We demonstrate our methods by studying the Creeping Section of the San Andreas fault observed by GPS and radar interferometry onboard the ERS and ALOS satellites. Along-strike variation of the moment accumulation rate is derived in order to investigate the degree of partial locking of the Creeping Section. The central Creeping Segment has a moment accumulation rate of 0.25–3.1  x  10 15 Nm yr –1 km –1 . The upper and lower bounds of the moment accumulation rates are derived based on the statistics of the noise. Our best-fitting model indicates that the central portion of the Creeping Section is accumulating seismic moment at rates that are about 5 per cent to 23 per cent of the fully locked Carrizo segment that will eventually be released seismically. A cumulative moment budget calculation with the historical earthquake catalogue ( M  〉 5.5) since 1857 shows that the net moment deficit at present is equivalent to a M w 6.3 earthquake.