Decarbonation reaction of magnesite in subducting slabs at the lower mantle

High-pressure and temperature experiments (28–62 GPa, and 1,490–2,000 K, corresponding to approximately 770–1,500 km depth in the mantle) have been conducted on a MgCO₃ + SiO₂ mixture using a laser-heated diamond anvil cell combined with analytical transmission electron microscope observation of the product phases to constrain the fate of carbonates carried on the subducting basalt into the lower mantle. At these conditions, the decarbonation reaction MgCO₃ (magnesite) + SiO₂ (stishovite) → MgSiO₃ (perovskite) + CO₂ (solid) has been recognized. This indicates that above reaction takes place as a candidate for decarbonation of the carbonated subducting mid ocean ridge basalts in the Earth’s lower mantle.     

Synchrotron IR study of hydrous ringwoodite (γ-Mg2SiO4) up to 30 GPa

High-pressure synchrotron infrared (IR) absorption spectra were collected between 650 and 4,000 cm⁻¹ at ambient temperature for hydrous Mg-ringwoodite (γ-Mg₂SiO₄) up to 30 GPa. The main feature in the OH⁻ stretching region is an extremely broad band centred at 3,150 cm⁻¹. The hydrogen bond is strong for most protons and the most probable site for protonation is the tetrahedral edge. With increasing pressure, this band shifts downward while decreasing its integrated intensity until disappearance at a pressure of 25 GPa. Only one band at 2,450 cm⁻¹ and an absorption plateau persist with a maximum wavenumber of 3,800 cm⁻¹. This behaviour is reversible upon pressure release. We interpret this as a second-order phase transition occurring in hydrated Mg-ringwoodite at high pressure (beyond ∼ 25 GPa). This result is compatible with the observation by Kleppe et al. (Phys Chem Miner 29:473–476, 2002a) who suggested the presence of Si–O–Si linkages and/or partial increase in the coordination of Si. Beyond the phase transition, the protons are delocalized and their environment on the ringwoodite structure is probably quite different from that at low pressure. Data obtained in situ at high pressures and temperatures are needed to better understand the effect of protonation on the structure and to better constrain this phase transition.     

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997–2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the “truth” signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.     

Spatial process and data models: Toward integration of agent-based models and GIS

The use of object-orientation for both spatial data and spatial process models facilitates their integration, which can allow exploration and explanation of spatial-temporal phenomena. In order to better understand how tight coupling might proceed and to evaluate the possible functional and efficiency gains from such a tight coupling, we identify four key relationships affecting how geographic data (fields and objects) and agent-based process models can interact: identity, causal, temporal and topological. We discuss approaches to implementing tight integration, focusing on a middleware approach that links existing GIS and ABM development platforms, and illustrate the need and approaches with example agent-based models.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

GGM02 – An improved Earth gravity field model from GRACE

A new generation of Earth gravity field models called GGM02 are derived using approximately 14 months of data spanning from April 2002 to December 2003 from the Gravity Recovery And Climate Experiment (GRACE). Relative to the preceding generation, GGM01, there have been improvements to the data products, the gravity estimation methods and the background models. Based on the calibrated covariances, GGM02 (both the GRACE-only model GGM02S and the combination model GGM02C) represents an improvement greater than a factor of two over the previous GGM01 models. Error estimates indicate a cumulative error less than 1 cm geoid height to spherical harmonic degree 70, which can be said to have met the GRACE minimum mission goals. Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Investigating the propagation mechanism of unmodelled systematic errors on coordinate time series estimated using least squares

The propagation of unmodelled systematic errors into coordinate time series computed using least squares is investigated, to improve the understanding of unexplained signals and apparent noise in geodetic (especially GPS) coordinate time series. Such coordinate time series are invariably based on a functional model linearised using only zero and first-order terms of a (Taylor) series expansion about the approximate coordinates of the unknown point. The effect of such truncation errors is investigated through the derivation of a generalised systematic error model for the simple case of range observations from a single known reference point to a point which is assumed to be at rest by the least squares model but is in fact in motion. The systematic error function for a one pseudo-satellite two-dimensional case, designed to be as simple but as analogous to GPS positioning as possible, is quantified. It is shown that the combination of a moving reference point and unmodelled periodic displacement at the unknown point of interest, due to ocean tide loading, for example, results in an output coordinate time series containing many periodic terms when only zero and first-order expansion terms are used in the linearisation of the functional model. The amplitude, phase and period of these terms is dependent on the input amplitude, the locations of the unknown point and reference point, and the period of the reference point's motion. The dominant output signals that arise due to truncation errors match those found in coordinate time series obtained from both simulated data and real three-dimensional GPS data.