Joint inversion estimate of regional glacial isostatic adjustment in Antarctica considering a lateral varying Earth structure (ESA STSE Project REGINA)

A major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry, and to a lesser extent satellite altimetry, is the poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA). Although much progress has been made in consistently modelling the ice-sheet evolution throughout the last glacial cycle, as well as the induced bedrock deformation caused by these load changes, forward models of GIA remain ambiguous due to the lack of observational constraints on the ice sheet's past extent and thickness and mantle rheology beneath the continent. As an alternative to forward modelling GIA, we estimate GIA from multiple space-geodetic observations: GRACE, Envisat/ICESat and GPS. Making use of the different sensitivities of the respective satellite observations to current and past surface mass (ice mass) change and solid Earth processes, we estimate GIA based on viscoelastic response functions to disc load forcing. We calculate and distribute the viscoelastic response functions according to estimates of the variability of lithosphere thickness and mantle viscosity in Antarctica. We compare our GIA estimate with published GIA corrections and evaluate its impact in determining the ice mass balance in Antarctica from GRACE and satellite altimetry. Particular focus is applied to the Amundsen Sea Sector in West Antarctica, where uplift rates of several cm/yr have been measured by GPS. We show that most of this uplift is caused by the rapid viscoelastic response to recent ice-load changes, enabled by the presence of a low-viscosity upper mantle in West Antarctica. This paper presents the second and final contribution summarizing the work carried out within a European Space Agency funded study, REGINA, (www.regina-science.eu)

Altimetry, gravimetry, GPS and viscoelastic modeling data for the joint inversion for glacial isostatic adjustment in Antarctica (ESA STSE Project REGINA)

The poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA) is a major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry and to a lesser extent satellite altimetry. In the past decade, much progress has been made in consistently modeling ice sheet and solid Earth interactions; however, forward-modeling solutions of GIA in Antarctica remain uncertain due to the sparsity of constraints on the ice sheet evolution, as well as the Earth's rheological properties. An alternative approach towards estimating GIA is the joint inversion of multiple satellite data – namely, satellite gravimetry, satellite altimetry and GPS, which reflect, with different sensitivities, trends in recent glacial changes and GIA. Crucial to the success of this approach is the accuracy of the space-geodetic data sets. Here, we present reprocessed rates of surface-ice elevation change (Envisat/Ice, Cloud,and land Elevation Satellite, ICESat; 2003–2009), gravity field change (Gravity Recovery and Climate Experiment, GRACE; 2003–2009) and bedrock uplift (GPS; 1995–2013). The data analysis is complemented by the forward modeling of viscoelastic response functions to disc load forcing, allowing us to relate GIA-induced surface displacements with gravity changes for different rheological parameters of the solid Earth. The data and modeling results presented here are available in the PANGAEA database (https://doi.org/10.1594/PANGAEA.875745). The data sets are the input streams for the joint inversion estimate of present-day ice-mass change and GIA, focusing on Antarctica. However, the methods, code and data provided in this paper can be used to solve other problems, such as volume balances of the Antarctic ice sheet, or can be applied to other geographical regions in the case of the viscoelastic response functions. This paper presents the first of two contributions summarizing the work carried out within a European Space Agency funded study: Regional glacial isostatic adjustment and CryoSat elevation rate corrections in Antarctica (REGINA)

A comparison of GOCE and drifter-based estimates of the North Atlantic steady-state surface circulation

Over the last decade, due to the Gravity Recovery And Climate Experiment (GRACE) mission and, more recently, the Gravity and steady state Ocean Circulation Explorer (GOCE) mission, our ability to measure the ocean’s mean dynamic topography (MDT) from space has improved dramatically. Here we use GOCE to measure surface current speeds in the North Atlantic and compare our results with a range of independent estimates that use drifter data to improve small scales. We find that, with filtering, GOCE can recover 70% of the Gulf Steam strength relative to the best drifter-based estimates. In the subpolar gyre the boundary currents obtained from GOCE are close to the drifter-based estimates. Crucial to this result is careful filtering which is required to remove small-scale errors, or noise, in the computed surface. We show that our heuristic noise metric, used to determine the degree of filtering, compares well with the quadratic sum of mean sea surface and formal geoid errors obtained from the error variance–covariance matrix associated with the GOCE gravity model. At a resolution of 100 km the North Atlantic mean GOCE MDT error before filtering is 5 cm with almost all of this coming from the GOCE gravity model

High frequency oscillations in epileptic and non-epileptic human hippocampus during a cognitive task

Hippocampal high-frequency electrographic activity (HFOs) represents one of the major discoveries not only in epilepsy research but also in cognitive science over the past few decades. A fundamental challenge, however, has been the fact that physiological HFOs associated with normal brain function overlap in frequency with pathological HFOs. We investigated the impact of a cognitive task on HFOs with the aim of improving differentiation between epileptic and non-epileptic hippocampi in humans. Hippocampal activity was recorded with depth electrodes in 15 patients with focal epilepsy during a resting period and subsequently during a cognitive task. HFOs in ripple and fast ripple frequency ranges were evaluated in both conditions, and their rate, spectral entropy, relative amplitude and duration were compared in epileptic and non-epileptic hippocampi. The similarity of HFOs properties recorded at rest in epileptic and non-epileptic hippocampi suggests that they cannot be used alone to distinguish between hippocampi. However, both ripples and fast ripples were observed with higher rates, higher relative amplitudes and longer durations at rest as well as during a cognitive task in epileptic compared with non-epileptic hippocampi. Moreover, during a cognitive task, significant reductions of HFOs rates were found in epileptic hippocampi. These reductions were not observed in non-epileptic hippocampi. Our results indicate that although both hippocampi generate HFOs with similar features that probably reflect non-pathological phenomena, it is possible to differentiate between epileptic and non-epileptic hippocampi using a simple odd-ball task

Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. Here we show that the geodetic estimation problem has deep-seated connections to Slepian's spatiospectral localization problem on the sphere, which amounts to finding bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures as bias, variance and mean-square error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of this new method makes the use of spherical Slepian functions in earth and planetary geodesy practical.Comment: 14 figures, submitted to the Geophysical Journal Internationa

Constraints from orbital motions around the Earth of the environmental fifth-force hypothesis for the OPERA superluminal neutrino phenomenology

It has been recently suggested by Dvali and Vikman that the superluminal neutrino phenomenology of the OPERA experiment may be due to an environmental feature of the Earth, naturally yielding a long-range fifth force of gravitational origin whose coupling with the neutrino is set by the scale M_*, in units of reduced Planck mass. Its characteristic length lambda should not be smaller than one Earth's radius R_e, while its upper bound is expected to be slightly smaller than the Earth-Moon distance (60 R_e). We analytically work out some orbital effects of a Yukawa-type fifth force for a test particle moving in the modified field of a central body. Our results are quite general since they are not restricted to any particular size of lambda; moreover, they are valid for an arbitrary orbital configuration of the particle, i.e. for any value of its eccentricity $e$. We find that the dimensionless strength coupling parameter alpha is constrained to |alpha| <= 1 10^-10-4 10^-9 for 1 R_e <= lambda <= 10 R_e by the laser data of the Earth's artificial satellite LAGEOS II, corresponding to M_* >= 4 10^9 -1.6 10^10. The Moon perigee allows to obtain |alpha| <= 3 10^-11 for the Earth-Moon pair in the range 15 R_e <= lambda = 3 10^10 - 4.5 10^10. Our results are neither necessarily limited to the superluminal OPERA scenario nor to the Dvali-Vikman model, in which it is M_* = 10^-6 at lambda = 1 R_e, in contrast with our bounds: they generally extend to any theoretical scenario implying a fifth-force of Yukawa-type.Comment: LaTex2e, 18 pages, 4 figures, 1 table, 81 reference

Development of a Synthetic Earth Gravity Model by 3D mass optimisation based on forward modelling

Several previous Synthetic Earth Gravity Model (SEGM) simulations are based on existing information about the Earth’s internal mass distribution. However, currently available information is insufficient to model the Earth’s anomalous gravity field on a global scale. The low-frequency information is missing when modelling only topography, bathymetry and crust (including the Mohorovičić discontinuity), but the inclusion of information on the mantle and core does not seem to significantly improve this situation. This paper presents a method to determine a more realistic SEGM by considering simulated 3D mass distributions within the upper mantle as a proxy for all unmodelled masses within the Earth.The aim is to improve an initial SEGM based on forward gravity modelling of the topography, bathymetry and crust such that the missing low-frequency information is now included. The simulated 3D mass distribution has been derived through an interactive and iterative mass model optimisation algorithm, which minimises geoid height differences with respect to a degree-360 spherical harmonic expansion of the EGM2008 global external gravity field model. We present the developed optimisation algorithm by applying it to the development of a global SEGM that gives a reasonably close fit to EGM2008, and certainly closer than a SEGM based only on the topography, bathymetry and crust

Spectral estimation on a sphere in geophysics and cosmology

We address the problem of estimating the spherical-harmonic power spectrum of a statistically isotropic scalar signal from noise-contaminated data on a region of the unit sphere. Three different methods of spectral estimation are considered: (i) the spherical analogue of the one-dimensional (1-D) periodogram, (ii) the maximum likelihood method, and (iii) a spherical analogue of the 1-D multitaper method. The periodogram exhibits strong spectral leakage, especially for small regions of area $A\ll 4\pi$, and is generally unsuitable for spherical spectral analysis applications, just as it is in 1-D. The maximum likelihood method is particularly useful in the case of nearly-whole-sphere coverage, $A\approx 4\pi$, and has been widely used in cosmology to estimate the spectrum of the cosmic microwave background radiation from spacecraft observations. The spherical multitaper method affords easy control over the fundamental trade-off between spectral resolution and variance, and is easily implemented regardless of the region size, requiring neither non-linear iteration nor large-scale matrix inversion. As a result, the method is ideally suited for most applications in geophysics, geodesy or planetary science, where the objective is to obtain a spatially localized estimate of the spectrum of a signal from noisy data within a pre-selected and typically small region.Comment: Submitted to the Geophysical Journal Internationa

Evaluation of the third- and fourth-generation GOCE Earth gravity field models with Australian terrestrial gravity data in spherical harmonics

In March 2013 the fourth generation of ESA’s (European Space Agency) global gravity field models, DIR4 (Bruinsma et al, 2010b) and TIM4 (Pail et al, 2010), generated from the GOCE (Gravity field and steady-state Ocean Circulation Explorer) gravity observation satellite were released. We evaluate the models using an independent ground truth data set of gravity anomalies over Australia. Combined with GRACE (Gravity Recovery and Climate Experiment) satellite gravity, a new gravity model is obtained that is used to perform comparisons with GOCE models in spherical harmonics. Over Australia, the new gravity model proves to have significantly higher accuracy in the degrees below 120 as compared to EGM2008 and seems to be at least comparable to the accuracy of this model between degree 150 and degree 260. Comparisons in terms of residual quasi-geoid heights, gravity disturbances, and radial gravity gradients evaluated on the ellipsoid and at approximate GOCE mean satellite altitude (h=250 km) show both fourth generation models to improve significantly w.r.t. their predecessors.Relatively, we find a root-mean-square improvement of 39 % for the DIR4 and 23 % for TIM4 over the respective third release models at a spatial scale of 100 km (degree 200). In terms of absolute errors TIM4 is found to perform slightly better in the bands from degree 120 up to degree 160 and DIR4 is found to perform slightly better than TIM4 from degree 170 up to degree 250. Our analyses cannot confirm the DIR4 formal error of 1 cm geoid height (0.35 mGal in terms of gravity) at degree 200. The formal errors of TIM4, with 3.2 cm geoid height (0.9 mGal in terms of gravity) at degree 200, seem to be realistic. Due to combination with GRACE and SLR data, the DIR models, at satellite altitude, clearly show lower RMS values compared to TIM models in the long wavelength part of the spectrum (below degree and order 120). Our study shows different spectral sensitivity of different functionals at ground level and at GOCE satellite altitude and establishes the link among these findings and the Meissl scheme (Rummel and van Gelderen in Manuscripta Geodaetica 20:379–385, 1995)