Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration

Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use.     

From subduction to transform motion: a seabeam survey of the Hellenic trench system

Preliminary results of a multi-narrow beam survey of the Hellenic trench system, in the Eastern Mediterranean, are presented. The southwestern Ionian branch is divided in small basins, partly filled with Pleistocene sediments. The morphology suggests that the basins are deformed by a compressional stress acting roughly perpendicularly to the trench along N50°E. This direction is the direction of the regional slip vector of the shallow thrust-type earthquakes. The structure of the southeastern Pliny-Strabo branch is quite different. Narrow en-e´chelon slots, oriented N40°E, have been mapped within the main troughs oriented N60°E. The regional earthquake slip vector is also oriented along N40°E. We conclude that the Hellenic trench system is an active subduction system, dominated by thrust along the Ionian branch and by transform motion along the Pliny-Strabo branch.     

Formation of Zn–Ca phyllomanganate nanoparticles in grass roots

It is now well established that a number of terrestrial and aquatic microorganisms have the capacity to oxidize and precipitate Mn as phyllomanganate. However, this biomineralization has never been shown to occur in plant tissues, nor has the structure of a natural Mn(IV) biooxide been characterized in detail. We show that the graminaceous plant Festuca rubra (red fescue) produces a Zn-rich phyllomanganate with constant Zn:Mn and Ca:Mn atomic ratios (0.46 and 0.38, respectively) when grown on a contaminated sediment. This new phase is so far the Zn-richest manganate known to form in nature (chalcophanite has a Zn:Mn ratio of 0.33) and has no synthetic equivalent. Visual examination of root fragments under a microscope shows black precipitates about ten to several tens of microns in size, and their imaging with backscattered and secondary electrons demonstrates that they are located in the root epidermis. In situ measurements by Mn and Zn K-edge extended X-ray absorption fine structure (EXAFS) spectroscopy and X-ray diffraction (XRD) with a micro-focused beam can be quantitatively described by a single-phase model consisting of Mn(IV) octahedral layers with 22% vacant sites capped with tetrahedral and octahedral Zn in proportions of 3:1. The layer charge deficit is also partly balanced by interlayer Mn and Ca. Diffracting crystallites have a domain radius of 33 Å in the ab plane and contain only 1.2 layers (8.6 Å) on average. Since this biogenic Mn oxide consists mostly of isolated layers, basal 00l reflections are essentially absent despite its lamellar structure. Individual Mn layers are probably held together in the Mn–Zn precipitates by stabilizing organic molecules. Zinc biomineralization by plants likely is a defense mechanism against toxicity induced by excess concentrations of this metal in the rhizosphere.     

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

A Top-Cut Model for Deposits with Heavy-Tailed Grade Distribution

In some ore deposits, the grade distribution is heavy-tailed and high values make the inference of first- and second-order statistics nonrobust. In the case of gold data, high values are usually cut and the block estimation is performed using truncated grades. With this method, a part of the deposit is omitted, resulting in a potential bias on resources figures. Ad-hoc corrections can be applied on the final figures to take into account the tail of the grade distribution, but no theoretical guidelines are available. A geostatistical model has been developed to handle high values based on the assumption that for high grade zones, the only tangible information is the geometry. The grade variable z can be split into three parts: the truncated grade ( $\operatorname{Min} (z, z_{\mathrm{e}})$ where z _e is the top-cut grade), a weighted indicator above top-cut grade (1{z≥z _e}), and a residual. Within this framework, the residual is poorly structured, and in most cases is a pure nugget-effect. Moreover, it has no spatial correlation with the truncated grade and the indicator above top-cut grade. This decomposition makes the variographic study more robust because variables (indicator and truncated grade) do not keep high grade values. The underlying hypotheses can be tested on experimental indicator variograms and the top-cut grade can be justified. Finally, the estimation is based on a truncated grade and indicator cokriging. The model is applied to blast holes from a gold deposit and on a simulated example. Both cases illustrate the benefits of keeping the high values in the estimation process via an appropriate mathematical model.     

U-series and oxygen isotope chronology of the mid-Pleistocene Lake Amora (Dead Sea basin)

This study establishes for the first time the chronology and limnological history of Lake Amora (Dead Sea basin, Israel), whose deposits (the Amora Formation) comprise one of the longest exposed lacustrine records of the Pleistocene time. The Amora Formation consists of sequences of laminated primary aragonite and silty-detritus, Ca-sulfate minerals, halite and clastic units. This sedimentary sequence was uplifted and tilted by the rising Sedom salt diapir, exposing ∼320 m of sediments on the eastern flanks of Mt. Sedom (the Arubotaim Cave (AC) section).The chronology of the AC section is based on U-disequilibrium dating (²³⁰Th-²³⁴U and ²³⁴U-²³⁸U ages) combined with floating δ¹⁸O stratigraphy and paleomagnetic constraints. The determination of the ²³⁰Th-²³⁴U ages required significant corrections to account for detrital Th and U. These corrections were performed on individual samples and on suites of samples from several stratigraphic horizons. The most reliable corrected ages were used to construct an age-elevation model that was further tuned to the oxygen isotope record of east Mediterranean foraminifers (based on the long-term similarity between the sea and lake oxygen isotope archives).The combined U-series-δ¹⁸O age-elevation model indicates that the (exposed) Amora sequence was deposited between ∼740 and 70 ka, covering seven glacial-interglacial cycles (Marine Isotope Stages (MIS) 18 to 5).Taking the last glacial Lake Lisan and the Holocene Dead Sea lacustrine systems as analogs of the depositional-limnological environment of Lake Amora, the latter oscillated between wet (glacial) and more arid (interglacial) conditions, represented by sequences of primary evaporites (aragonite and gypsum that require enhanced supply of freshwater to the lakes) and clastic sediments, respectively. The lake evolved from a stage of rapid shifts between high and low-stand conditions during ∼740 to 550 ka to a sabkha-like environment that existed (at the AC site) between 550 and 420 ka. This stage was terminated by a dry spell represented by massive halite deposition at 420 ka (MIS12-11). During MIS10-6 the lake fluctuated between lower and higher stands reaching its highest stand conditions at the late glacial MIS6, after which a significant lake level decline corresponds to the transition to the last interglacial (MIS5) low-stand lake, represented by the uppermost part of the Formation.δ¹⁸O values in the primary aragonite range between 6.0 and −1.3‰, shifting cyclically between glacial and interglacial intervals. The lowest δ¹⁸O values are observed during interglacial stages and may reflect short and intense humid episodes that intermittently interrupted the overall arid conditions. These humid episodes, expressed also by enhanced deposition of travertines and speleothems, seem to characterize the Negev Desert, and in contrast to the overall dominance of the Atlantic-Mediterranean system of rain patterns in the Dead Sea basin, some humid episodes during interglacials may be traced to southern sources.     

Middle Carboniferous crustal melting in the Variscan Belt: New insights from U–Th–Pbtot. monazite and U–Pb zircon ages of the Montagne Noire Axial Zone (southern French Massif Central)

In France, the Devonian–Carboniferous Variscan orogeny developed at the expense of continental crust belonging to the northern margin of Gondwana. A Visean–Serpukhovian crustal melting has been recently documented in several massifs. However, in the Montagne Noire of the Variscan French Massif Central, which is the largest area involved in this partial melting episode, the age of migmatization was not clearly settled. Eleven U–Th–Pb_tot. ages on monazite and three U–Pb ages on associated zircon are reported from migmatites (La Salvetat, Ourtigas), anatectic granitoids (Laouzas, Montalet) and post-migmatitic granites (Anglès, Vialais, Soulié) from the Montagne Noire Axial Zone are presented here for the first time. Migmatization and emplacement of anatectic granitoids took place around 333–326 Ma (Visean) and late granitoids emplaced around 325–318 Ma (Serpukhovian). Inherited zircons and monazite date the orthogneiss source rock of the Late Visean melts between 560 Ma and 480 Ma. In migmatites and anatectic granites, inherited crystals dominate the zircon populations. The migmatitization is the middle crust expression of a pervasive Visean crustal melting event also represented by the “Tufs anthracifères” volcanism in the northern Massif Central. This crustal melting is widespread in the French Variscan belt, though it is restricted to the upper plate of the collision belt. A mantle input appears as a likely mechanism to release the heat necessary to trigger the melting of the Variscan middle crust at a continental scale.