Time Variations In Gravity And Inferences On The Earth's Structure And Dynamics

This paper reviews how the study of the surface gravity changes is able to provide useful information on the Earth's structure and global dynamics. The spectral range which is observable with superconducting gravimeters is broad and goes from the seismic frequency band to periods longer than one year. We first investigate the seismic and sub-seismic bands with a special attention paid to the gravity detection of core modes in the liquid core and to the Slichter mode of translation of the solid inner core. In the tidal bands, we show how accurate measurements allow us to infer constraints on various phenomena such as mantle (an-)elasticity, as well as ocean and atmospheric loading. The observation of the Free Core Nutation resonance in the diurnal frequency band is reviewed and indirectly suggests an increase in the ellipticity of the core-mantle boundary with respect to its hydrostatic value. A similar resonance is also theoretically predicted in the diurnal band for the rotation of the solid inner core (Free Inner Core Nutation) but we show that its detection is much more difficult because of the small amplitude and lack of a nearby tidal frequency. Oceanic and atmospheric loading mechanisms induce gravity changes over a wide spectral range and we present some recent progress in this field. Finally, because superconducting gravimeters have high calibration stability and small long-term instrumental drift, they can easily measure longperiod gravity variations due to polar motion and hydrogeology.     

A modification to the Soil Conservation Service curve number method for steep slopes in the Loess Plateau of China

The Soil Conservation Service curve number (CN) method is widely used for predicting direct runoff from rainfall. However, despite the extent of cultivation on hillslope areas, very few attempts have been made to incorporate a slope factor into the CN method. The objectives of this study were (1) to evaluate existing approaches integrating slope in the CN method, and (2) to develop an equation incorporating a slope factor into the CN method for application in the steep slope areas of the Loess Plateau of China. The dataset consisted of 11 years of rainfall and runoff measurements from two experimental sites with slopes ranging from 14 to 140%. The results indicated that the standard CN method underestimated large runoff events and overestimated small events. For our experimental conditions, the optimized and non‐optimized forms of the slope‐modified CN method of the Erosion Productivity Impact Calculator model improved runoff prediction for steep slopes, but large runoff events were still underestimated and small ones overpredicted. Based on relationships between slope and the observed and theoretical CN values, an equation was developed that better predicted runoff depths with an R² of 0·822 and a linear regression slope of 0·807. This slope‐adjusted CN equation appears to be the most appropriate for runoff prediction in the steep areas of the Loess Plateau of China. Copyright © 2005 John Wiley & Sons, Ltd.     

Comparison of reduced-order, sequential and variational data assimilation methods in the tropical Pacific Ocean

This paper presents a comparison of two reduced-order, sequential, and variational data assimilation methods: the singular evolutive extended Kalman filter (SEEK) and the reduced 4D-Var (R-4D-Var). A hybridization of the two, combining the variational framework and the sequential evolution of covariance matrices, is also preliminarily investigated and assessed in the same experimental conditions. The comparison is performed using the twin-experiment approach on a model of the tropical Pacific domain. The assimilated data are simulated temperature profiles at the locations of the TAO/TRITON array moorings. It is shown that, in a quasilinear regime, both methods produce similarly good results. However, the hybrid approach provides slightly better results and thus appears as potentially fruitful. In a more nonlinear regime, when tropical instability waves develop, the global nature of the variational approach helps control model dynamics better than the sequential approach of the SEEK filter. This aspect is probably enhanced by the context of the experiments in that there is a limited amount of assimilated data and no model error.     

General database for ground water site information

In most cases, analysis and modeling of flow and transport dynamics in ground water systems require long-term, high-quality, and multisource data sets. This paper discusses the structure of a multisite database (the H+ database) developed within the scope of the ERO program (French Environmental Research Observatory, http://www.ore.fr). The database provides an interface between field experimentalists and modelers, which can be used on a daily basis. The database structure enables the storage of a large number of data and data types collected from a given site or multiple-site network. The database is well suited to the integration, backup, and retrieval of data for flow and transport modeling in heterogeneous aquifers. It relies on the definition of standards and uses a templated structure, such that any type of geolocalized data obtained from wells, hydrological stations, and meteorological stations can be handled. New types of platforms other than wells, hydrological stations, and meteorological stations, and new types of experiments and/or parameters could easily be added without modifying the database structure. Thus, we propose that the database structure could be used as a template for designing databases for complex sites. An example application is the H+ database, which gathers data collected from a network of hydrogeological sites associated with the French Environmental Research Observatory.     

Single-station estimates of the seismic moment of the 1960 Chilean and 1964 Alaskan earthquakes,using the mantle magnitudeM m

Measurements are taken of the mantle magnitudeM _m , developed and introduced in previous papers, in the case of the 1960 Chilean and 1964 Alaskan earthquakes, by far the largest events ever recorded instrumentally. We show that theM _m algorithm recovers the seismic moment of these gigantic earthquakes with an accuracy (typically 0.2 to 0.3 units of magnitude, or a factor of 1.5 to 2 on the seismic moment) comparable to that achieved on modern, digital, datasets. In particular, this study proves that the mantle magnitudeM _m does not saturate for large events, as do standard magnitude scales, but rather keeps growing with seismic moment, even for the very largest earthquakes. We further prove that the algorithm can be applied in unfavorable experimental conditions, such as instruments with poor response at mantle periods, seismograms clipped due to limited recording dynamics, or even on microbarograph records of air coupled Rayleigh waves.In addition, we show that it is feasible to use acoustic-gravity air waves generated by those very largest earthquakes, to obtain an estimate of the seismic moment of the event along the general philosophy of the magnitude concept: a single-station measurement ignoring the details of the earthquake's focal mechanism and exact depth.     

Volcanological and petrological evolution of Nyiragongo volcano,Virunga volcanic field,Zaire

Three major phases are distinguished during the growth of Nyiragongo, an active volcano at the western limit of the Virunga Range, Zaire. Lavas erupted during phase 1 are strongly undersaturated melilitites characterized by the presence of kalsilite phenocrysts, perovskite, and the abundance of calcite in the matrix. Such lavas crop out mainly on the inner crater wall and progressively evolve toward more aphyric melilite nephelinites well represented on the flanks of the volcano. Adventive vents lying at the base of the cone developed along radial fracture systems and erupted olivine and/or clinopyroxene – rich melilitites or nephelinites. Stage 2 lavas are melilite-free nephelinites. Clinopyroxene is the main phenocryst and feldspathoids are abundant in the lavas exposed on the crater wall. These flows result from periodic overflowing of a magma column from an open crater. Extensive fissure flows which erupted from the base of the cone at the end of this stage are related to widespread draining out of magma which in turn induces the formation of the summit pit crater. Magmas erupted during stage 3 are relatively aphyric melilite nephelinites and the main volcanological characteristic is the permanent lava lake observed into the pit crater until the 1977 eruption. Fluctuations of the level of the lava lake was responsible for the development of the inner terraces. Periodic overflowing of the lava lake from the central pit formed the nepheline aggregate lava flows. Petrography and major element geochemistry allow the determination of the principal petrogenetic processes. Melilitites and nephelinites erupted from the summit crater are lavas derived, via clinopyroxene fractionation, from a more primitive melt. The abundance of feldspathoids in these lavas is in keeping with nepheline flotation. Aphyric melilite nephelinites covering the flanks and the extensive fissure flows have a homogeneous chemical composition; rocks from the historical lava lake are slightly more evolved. All these lavas differentiated in a shallow reservoir. Lavas erupted from the parasitic vents are mainly olivine and/or clinopyroxene-phyric rocks. Rushayite and picrites from Muja cone are peculiar high-magnesium lavas resulting from the addition of olivine xenocrysts to melilitic or nephelinitic melts. Fluid and melt inclusions in olivine and clinopyroxene phenocrysts indicate a crystallization depth of 10–14 km. A model involving two reservoirs located at different depths and periodically connected is proposed to explain the petrography of the lavas; this hypothesis is in accordance with geophysical data. Received: July 8, 1993/Accepted: September 10, 1993     

Achievements of the Earth orientation parameters prediction comparison campaign

Precise transformations between the international celestial and terrestrial reference frames are needed for many advanced geodetic and astronomical tasks including positioning and navigation on Earth and in space. To perform this transformation at the time of observation, that is for real-time applications, accurate predictions of the Earth orientation parameters (EOP) are needed. The Earth orientation parameters prediction comparison campaign (EOP PCC) that started in October 2005 was organized for the purpose of assessing the accuracy of EOP predictions. This paper summarizes the results of the EOP PCC after nearly two and a half years of operational activity. The ultra short-term (predictions to 10 days into the future), short-term (30 days), and medium-term (500 days) EOP predictions submitted by the participants were evaluated by the same statistical technique based on the mean absolute prediction error using the IERS EOP 05 C04 series as a reference. A combined series of EOP predictions computed as a weighted mean of all submissions available at a given prediction epoch was also evaluated. The combined series is shown to perform very well, as do some of the individual series, especially those using atmospheric angular momentum forecasts. A main conclusion of the EOP PCC is that no single prediction technique performs the best for all EOP components and all prediction intervals.     

Combination of temporal gravity variations resulting from superconducting gravimeter (SG) recordings, GRACE satellite observations and global hydrology models

Gravity recovery and climate experiment (GRACE)-derived temporal gravity variations can be resolved within the μgal (10^?8 m/s ²) range, if we restrict the spatial resolution to a half-wavelength of about 1,500 km and the temporal resolution to 1 month. For independent validations, a comparison with ground gravity measurements is of fundamental interest. For this purpose, data from selected superconducting gravimeter (SG) stations forming the Global Geodynamics Project (GGP) network are used. For comparison, GRACE and SG data sets are reduced for the same known gravity effects due to Earth and ocean tides, pole tide and atmosphere. In contrast to GRACE, the SG also measures gravity changes due to load-induced height variations, whereas the satellite-derived models do not contain this effect. For a solid spherical harmonic decomposition of the gravity field, this load effect can be modelled using degree-dependent load Love numbers, and this effect is added to the satellite-derived models. After reduction of the known gravity effects from both data sets, the remaining part can mainly be assumed to represent mass changes in terrestrial water storage. Therefore, gravity variations derived from global hydrological models are applied to verify the SG and GRACE results. Conversely, the hydrology models can be checked by gravity variations determined from GRACE and SG observations. Such a comparison shows quite a good agreement between gravity variation derived from SG, GRACE and hydrology models, which lie within their estimated error limits for most of the studied SG locations. It is shown that the SG gravity variations (point measurements) are representative for a large area within the accuracy, if local gravity effects are removed. The individual discrepancies between SG, GRACE and hydrology models may give hints for further investigations of each data series.     

Different mechanisms of the ionospheric influence on GPS occultation signals

A local mechanism for strong ionospheric effects on radio occultation (RO) global positioning satellite system (GPS) signals is described. Peculiar zones centered at the critical points (the tangent points) in the ionosphere, where the gradient of the electron density is perpendicular to the RO ray trajectory, strongly influence the amplitude and phase of RO signals. It follows from the analytical model of local ionospheric effects that the positions of the critical points depend on the RO geometry and the structure of the ionospheric disturbances. Centers of strong ionospheric influence on RO signals can exist, for example, in the sporadic E-layers, which are inclined by 3–6° relative to the local horizontal direction. Also, intense F2 layer irregularities can contribute to the RO signal variations. A classification of the ionospheric influence on the GPS RO signals is introduced using the amplitude data, which indicates different mechanisms (local, diffraction, etc.) for radio waves propagation. The existence of regular mechanisms (e.g., local mechanism) indicates a potential for separating the regular and random parts in the ionospheric influence on the RO signals.     

Focal mechanism solutions and nature of plate movements in Pakistan

From the seismic point of view, the territory of Pakistan which lies between latitude 23°–37° N and longitude 61°–75° E is one of the most active zones in the world. The importance of this area lies in terms of movements of the Indian plate with respect to Eurasia on the west. Seismicity, as well as focal mechanism- solutions, throws a considerable light on the nature of forces acting in the area. All the available solutions, along with 12 new ones, have been considered for the present study. Their relationship to major faults in the area is discussed. The majority of the solutions in the central and northern parts show strike-slip faulting with a left-lateral sense of motion, followed by thrust faulting; few show normal faulting. This suggests that the Indian plate is moving with respect to the Eurasian plate along the Chaman fault, Quetta transverse zone, Sulaiman Ranges and the Hazara thrusts region joining the Hazara/Kashmir syntaxis. The orientations of P and T axes have been studied. It is seen that in a large number of cases compressive stress is acting nearly in NNW-SSE to N-S directions. The Hazara thrust region appears to be the most complex. Here, the influence of the Himalayan thrust front is evident to a large extent.The nature of faulting along the Chaman fault and Quetta transverse zone is to some extent similar to that of the San Andreas fault system of California. So far as the energy release is concerned, the maximum energy is being released in the form of strike-slip movements close to the Chaman fault and Quetta transverse ranges.