Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

A systematic algorithm for obtaining SODAR wind profiles

Summary Current methods of obtaining wind profiles from SODARs generally employ ad hoc strategies for Doppler spectrum peak selection, profile smoothing, and consistency checks. This paper addresses the problem of systematically obtaining vertical profiles of wind vectors from a SODAR. A new approach is to regard every spectrum frequency as a possible estimator of the radial wind component, but to weight these Doppler estimates according to the spectral power at each frequency. This allows the peak detection, smoothing, and assimilation of a priori information to be combined into a constrained linear inversion methodology. The inversion process also allows an estimation of signal-to-noise ratio, degrees of freedom, and information content versus height. Examples are given of profiles obtained using this new inverse method. The relationship to the physical limitations of the SODAR instrument and the atmosphere is discussed.     

Petrogenesis of orthopyroxene- and amphibole-bearing andesites, Mustique, Grenadine Islands, Lesser Antilles Arc: Isotope, trace element and physical constraints

Abstract On the island of Mustique, fresh and propylitized olivine–plagioclase–clinopyroxene basalt, plagioclase–clinopyroxene–orthopyroxene and plagioclase–clinopyroxene–amphibole andesite lavas and minor intrusions are interbedded with Oligocene pyroclastic and epiclastic rocks. Chemical data show that two isotopically identical, but chemically different, suites of lava are present: (i) the OPXS (⁸⁷Sr/⁸⁶Sr 0.70403–0.70454; ¹⁴³Nd/¹⁴⁴Nd 0.512952–0.512986; δ¹⁸O_cpx 5.49 and 5.61), comprising basalts and orthopyroxene‐bearing andesites; and (ii) the AMPHS (⁸⁷Sr/⁸⁶Sr 0.70401–0.70457; ¹⁴³Nd/¹⁴⁴Nd 0.512981–0.513037; δ¹⁸O_cpx 5.54), made up of basalts and amphibole‐bearing andesites. The OPXS has higher contents of TiO₂, P₂O₅, light rare earth elements, Sm, Pb, Th, U, Zr, Y and Nb, and higher La/Yb ratios than the AMPHS. The isotopic data suggest that both suites formed from melts derived from the same subduction‐modified depleted mantle source as the volcanic rocks of nearby St Vincent and Bequia, and the northern islands of the Lesser Antilles Arc. The immobile trace element contents, and La/Yb ratios, of the OPXS are indicative of ～10% partial melting of the source, whereas those of the AMPHS are indicative of ～25% partial melting. The within‐suite chemical variation of the OPXS is consistent with ～45% fractional crystallization of its intratelluric mineral assemblages, and that of the AMPHS is consistent with the removal of ～65% of its intratelluric assemblages. Experimental evidence suggests that both suites of basalt crystallized at pressures <8 kbar from melts containing 1–2 wt% water. After extensive fractional crystallization, the andesites crystallized at pressures between approximately 5 and 2 kbar. The OPXS magmas appear to have lost more of their water content than the AMPHS magmas. Thus, the OPXS andesites formed from melts with an estimated water content of 2–3 wt%, whereas the AMPHS andesites formed from melts containing at least 4.5 wt% water.     

Influence of water level fluctuation on the spatial dynamic of the water chemistry at Lake Fertõ/Neusiedler See

Lake Fertõ/Neusiedler See is a steppe lake situated on the Hungarian-Austrian border with a surface area of 309 km². It is an extremely shallow lake with regulated outflow. 54% of the whole lake and 85% of the Hungarian part (75 km²) is covered by reed. There are numerous reedless areas (inner ponds) of variable size within the reed belt, which is enmeshed with canals connecting the inner ponds with the open water areas. Fundamental differences were recorded in the water chemistry of the three water type areas, open water, canals and inner ponds. The degree of these differences is basically determined by the distance from the open water areas and the water level of the actual years. An extremely high, an average and an extremely low water level year was examined, in order to present the spatial and temporal differences in the water chemistry of the investigated water bodies. Differences due to the water level fluctuation between years in the three different water types are discussed on the basis of electrical conductivity, cation, and anion concentration changes of the water. Our results support the concept of high instability of Lake Fertõ and describe habitats with remarkable temporal and spatial variability. The high influence of water level on habitat conditions at different sites of the lake recommend a carefull water level regulation strategy.     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.     

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.     

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).