The 1988–2003 Greenland ice sheet melt extent using passive microwave satellite data and a regional climate model

Measurements from ETH-Camp and JAR1 AWS (West Greenland) as well as coupled atmosphere-snow regional climate simulations have highlighted flaws in the cross-polarized gradient ratio (XPGR) technique used to identify melt from passive microwave satellite data. It was found that dense clouds (causing notably rainfall) on the ice sheet severely perturb the XPGR melt signal. Therefore, the original XPGR melt detection algorithm has been adapted to better incorporate atmospheric variability over the ice sheet and an updated melt trend for the 1988–2003 period has been calculated. Compared to the original algorithm, the melt zone area increase is eight times higher (from 0.2 to 1.7% year⁻¹). The increase is higher with the improved XPGR technique because rainfall also increased during this period. It is correlated to higher atmospheric temperatures. Finally, the model shows that the total ice sheet runoff is directly proportional to the melt extent surface detected by satellites. These results are important for the understanding of the effect of Greenland melting on the stability of the thermohaline circulation.     

Interaction between an eddy and a zonal jet : Part II. Two-and-a-half-layer model

In a two-and-a-half-layer quasi-geostrophic model, a process study is conducted on the interaction between a vortex and a zonal jet, both with constant potential vorticity. The vortex is a stable anticyclone, initially located north of the eastward jet. The potential vorticity of the jet is allowed to have various vertical structures, while the vortex is concentrated in only one layer. The flow parameters are set to values characteristic of the Azores region.First, the jet is stable. Weak vortices steadily drift north of the jet without crossing it while strong vortices can cross the jet and tear off a cyclone with which they pair as a heton (baroclinic dipole). This heton often breaks later in the shear exerted by the jet; the two vortices finally drift apart. When crossed by deep anticyclones, the jet develops meanders with 375 km wavelength. These results exhibit a noticeable similarity with the one-and-a-half-layer case studied in Part I.Secondly, the jet is allowed to be linearly unstable. In the absence of the vortex, it develops meanders with 175 km wavelength and 25-day e-folding time on the β-plane. For various vertical structures of the jet, baroclinic instability is shown to barely affect jet–vortex interaction if the linear growth rate of unstable waves is smaller than 1/(14 days). Further simulations with a linearly unstable, nonlinearly equilibrated jet evidence its strong temporal variability when crossed by a deep vortex on the β-plane. In particular, long waves can dominate the spectrum for a few months after jet crossing by the vortex. Again in this process, the deep vortex couples with a surface cyclone and both drift southwestward.     

On matrix diffusion: formulations, solution methods and qualitative effects

 Matrix diffusion has become widely recognized as an important transport mechanism. Unfortunately, accounting for matrix diffusion complicates solute-transport simulations. This problem has led to simplified formulations, partly motivated by the solution method. As a result, some confusion has been generated about how to properly pose the problem. One of the objectives of this work is to find some unity among existing formulations and solution methods. In doing so, some asymptotic properties of matrix diffusion are derived. Specifically, early-time behavior (short tests) depends only on φ _m ² R _m D _m  / L _m ², whereas late-time behavior (long tracer tests) depends only on φ _m R _m , and not on matrix diffusion coefficient or block size and shape. The latter is always true for mean arrival time. These properties help in: (a) analyzing the qualitative behavior of matrix diffusion; (b) explaining one paradox of solute transport through fractured rocks (the apparent dependence of porosity on travel time); (c) discriminating between matrix diffusion and other problems (such as kinetic sorption or heterogeneity); and (d) describing identifiability problems and ways to overcome them. Received, October 1997 · Revised, November 1997 · Accepted, December 1997     

Source spectral phases of surface waves recovered from observed spectral amplitudes by using Hilbert transforms

Spectral phases of surface waves radiated at the seismic source are recovered by means of the Hilbert transform of spectral amplitudes. Through numerical simulations as well as the study of an earthquake with well known focal mechanism, it is shown that the recovered phases satisfy a criterion of coherency, so that they can be used in an inversion process for the retrieval of the seismic moment tensor.     

Generation of earth stratified models compatible with both ellipticity and phase velocity observations of Rayleigh waves

In the eigenvalue eigenvector decomposition used in inversion techniques, we can sometimes obtain eigenvalues close to zero. Related to these eigenvalues close to zero, a set of stratified earth models is generated. These are compatible with the theoretical observations according to a given model. The method is applied to ellipticities and phase velocities of Rayleigh waves, and supply proof of the advantage of working with both sets of data together.     

Thrusting of young lithosphere in subduction zones with special reference to structures in ophiolitic peridotites

Two contrasted types of structures have been recognized in peridotites from ophiolites and from the oceanic environment. The first one, typical of high-temperature/moderate-stress conditions, is observed in the upper part of ophiolitic peridotites and has been ascribed to plastic flow in an oceanic ridge environment. The second one, typical of moderate-temperature/high-stress conditions, is more specially dealt with here. It is printed in the peridotites above the basal metamorphic aureole found in many ophiolites. The strain increases downward over 1–2 km to produce peridotite mylonites at the contact with the metamorphic aureole. Similarities with rocks from trench and island arc environments suggest ascribing this deformation in ophiolites to a trench environment. We propose that shear fracturing in a young oceanic lithosphere is initiated by the compressive elastic stress in its lower part which is produced by bending of the subducted plate. An externally applied compressive stress is responsible for subsequent overthrusting of the fractured lithosphere. This interpretation is in good agreement with the available geophysical data on young subducted plates and with the physical data on ophiolitic peridotites.     

A new tool for in situ monitoring of Fe-mobilization in soils

The aim of this study was to design and test a new tool for (i) the quantitative in situ monitoring of Fe(III) reduction in soils and (ii) the tracking of the potential mineralogical changes of Fe-oxides. The tool consists of small (2 × 2 × 0.2 cm) striated polymer plates coated with synthetic pure ferrihydrite or As-doped ferrihydrite (Fh–As). These slides were then inserted within two different horizons (organo-mineral and albic) located in a wetland soil with alternating redox conditions. Dissolution was quantified by X-ray fluorescence (XRF) analyses of total metal contents before and after insertion into the soil. The crystallographic evolution of Fe-oxides was characterized by scanning electron microscope equipped with an energy-dispersive spectrometer (SEM–EDS). Over the months, the ferrihydrite progressively disappeared, at rates comparable to those previously measured in laboratory studies, i.e. in the 1–10 × 10⁻¹² mol Fe m⁻² s⁻¹ range. SEM observations indicate that the supports were highly colonized by bacteria and biofilms in the organo-mineral horizon, suggesting a biological-mediated process, while the albic horizon appeared to be characterized by a mostly chemical-mediated process. In the albic horizon, Fe-sulphide and other micro-precipitates were formed after 7 months of incubation in balance with a quasi dissolution of initial Fe-oxides.     

A spectral approach to simulating intrinsic random fields with power and spline generalized covariances

This article presents a variant of the spectral turning bands method that allows fast and accurate simulation of intrinsic random fields with power, spline, or logarithmic generalized covariances. The method is applicable in any workspace dimension and is not restricted in the number and configuration of the locations where the random field is simulated; in particular, it does not require these locations to be regularly spaced. On the basis of the central limit and Berry–Esséen theorems, an upper bound is derived for the Kolmogorov distance between the distributions of generalized increments of the simulated random fields and the normal distribution.     

Regionalized Classification of Geochemical Data with Filtering of Measurement Noises for Predictive Lithological Mapping

Natural Resources Research - A method for predictive lithological mapping is proposed, which combines geostatistical simulation of geochemical concentrations with coregionalization analysis and...