Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Effet du vieillissement des céramiques poreuses sur leur capacité à évaluer la concentration de pesticide en solution

The soil solution sampling by ceramic cups allows pesticide transfer monitoring in the soil during long times. The ageing of material involves a bias in the sampling results. In laboratory, the comparison of two types of ceramic suction cups, new and installed in situ during four years, shows a modification of the hydrodynamic properties and a possible evolution of the adsorption capacity of the matrix. The passage rate, as well qualitative as quantitative, is better for the old material. Recommendations about site management are finally exposed. To cite this article: N. Domange et al., C. R. Geoscience 336 (2004).     

Nouvelle interprétation structurale de la « faille Nord- Pyrénéenne » en vallée d'Aspe (Pyrénées-Atlantiques). Remise en question d'un plutonisme ophitique danien dans le secteur de Bedous

In the Aspe Valley (western Pyrenees), the Europe/Iberia boundary corresponds to a complex fracturing zone, called the ‘Bielle–Accous Wrench-Faulting Corridor’, which represents the classical ‘North-Pyrenean Fault’. Located between the High Primary Range and the North-Pyrenean Zone, the BAWC shows different south-verging sheets mainly composed of Triassic materials. The Bedous ophite, associated with Muschelkalk and Keuper sediments, is also Triassic in age and involved in the same Pyrenean thrusting structures. So, contrary to a recent interpretation, this magmatic rock cannot be related to a supposed Danian plutonism inducing metamorphic processes in the surrounding Mesozoic formations. To cite this article: J. Canérot et al., C. R. Geoscience 336 (2004).     

Nickel speciation in Sebertia acuminata,a plant growing on a lateritic soil of New Caledonia

Nickel speciation in a nickel hyperaccumulating plant (Sebertia acuminata) and its associated soil of southern New Caledonia was studied using various analytical methods. The soil is formed of iron oxides (goethite, hematite), which contain almost all the nickel. The available nickel is probably linked to the organic matter in the litter. Sebertia acuminata, acts as a nickel pump, and concentrates the metal in its leaves. It partitions nickel and silica; nickel is concentrated in the cells (probably in the vacuoles) as organometallic complexes, whereas silica forms the framework of the cells, and the phytolithes. A thorough study of these plants seems essential in order to define the soil–plant relations, and to propose appropriate ways for ecological restoration. To cite this article: N. Perrier et al., C. R. Geoscience 336 (2004).     

Dérive à la surface de l'océan sous l'effet des vagues

We model the drift velocity near the ocean surface separating the motion induced by the local current, itself influenced by winds and waves, and the motion induced by the waves, which are generated by local and remote winds. Application to the drift of ‘tar balls’, following the sinking of the oil tanker Prestige-Nassau in November 2002, shows that waves contribute at least one third of the drift for pollutants floating 1 m below the surface, with a mean direction about 30° to the right of the wind-sea direction. Although not new, this result was previously obtained with specific models, whereas the formalism used here combines classical wave and circulation forecasting models. To cite this article: F. Ardhuin et al., C. R. Geoscience 336 (2004).     

Effet de la teneur en eau d'un sol sur la réduction bactérienne d'oxydes de fer

Iron-reducing activity of autochthonous bacteria from two temporary hydromorphic soils is evaluated by the study of iron reductive dissolution, as a function of water content. The release of ferrous iron in solution is coupled to the mineralization of soil organic carbon. Water soil saturation is not necessary for iron reductive dissolution, since the highest dissolution is obtained for a wet, but not water-saturated soil (100% of water holding capacity WHC), and dissolution is also very high in a soil at 75% WHC. To cite this article: S.J. Stemmler et al., C. R. Geoscience 336 (2004).     

Mise en évidence d'un nouveau front de chevauchement dans l'Atlas tunisien oriental de Tunisie par sismique réflexion. Contexte structural régional et rôle du Trias salifère

Structural interpretations of newly acquired seismic lines in northeastern Tunisia allow us to highlight a new thrust front for the Atlasic range of Tunisia, in contrast to the previously Zaghouan fault thrust Dorsale zone. This new thrust front takes place on weakness tectonic zones, materialized by inherited faults anchored on the pre-Triassic basement. This front seems to be a paleogeographic trend controlling structural style and basin fill with a synsedimentary activity. The front is expressed by reverse faults, thrust faults, back thrusting, and decollement structures. To cite this article: S. Khomsi et al., C. R. Geoscience 336 (2004).     

Seismic signature of a phreatic explosion: hydrofracturing damage at Karthala volcano, Grande Comore Island, Indian Ocean

Karthala volcano is a basaltic shield volcano with an active hydrothermal system that forms the southern two-thirds of the Grande Comore Island, off the east coat of Africa, northwest of Madagascar. Since the start of volcano monitoring by the local volcano observatory in 1988, the July 11th, 1991 phreatic eruption was the first volcanic event seismically recorded on this volcano, and a rare example of a monitored basaltic shield. From 1991 to 1995 the VT locations, 0.5<M_l<4.3, show a crack shaped pattern (3 km long, 1 km wide) within the summit caldera extending at depth from –2 km to +2 km relative to sea level. This N-S elongated pattern coincides with the direction of the regional maximum horizontal stress as deduced from regional focal mechanism solutions. This brittle signature of the damage associated with the 1991 phreatic eruption is a typical pattern of the seismicity induced by controlled fluid injections such as those applied at geothermal fields, in oil and gas recovery, or for stress measurements. It suggests the 1991 phreatic eruption was driven by hydraulic fracturing induced by forced fluid flow. We propose that the extremely high LP and VT seismicity rates, relative to other effusive volcanoes, during the climax of the 1991 phreatic explosion, are due to the activation of the whole hydrothermal system, as roughly sized by the distribution of VT hypocenters. The seismicity rate in 1995 was still higher than the pre-eruption seismicity rate, and disagrees with the time pattern of thermo-elastic stress readjustment induced by single magma intrusions at basaltic volcanoes. We propose that it corresponds to the still ongoing relaxation of pressure heterogeneity within the hydrothermal system as suggested by the few LP events that still occurred in 1995.Editorial responsibility: H Shinohara