Space-time mapping of soil salinity using probabilistic bayesian maximum entropy

The mapping of saline soils is the first task before any reclamation effort. Reclamation is based on the knowledge of soil salinity in space and how it evolves with time. Soil salinity is traditionally determined by soil sampling and laboratory analysis. Recently, it became possible to complement these hard data with soft secondary data made available using field sensors like electrode probes. In this study, we had two data sets. The first includes measurements of field salinity (EC_a) at 413 locations and 19 time instants. The second, which is a subset of the first (13 to 20 locations), contains, in addition to EC_a, salinity determined in the laboratory (EC_2.5). Based on a procedure of cross-validation, we compared the prediction performance in the space-time domain of 3 methods: kriging using either only hard data (HK) or hard and mid interval soft data (HMIK), and Bayesian maximum entropy (BME) using probabilistic soft data. We found that BME was less biased, more accurate and giving estimates, which were better correlated with the observed values than the two kriging techniques. In addition, BME allowed one to delineate with better detail saline from non-saline areas.     

Managed water levels and the expansion of emergent vegetation along a lakeshore

Water-level fluctuations may be used to promote the expansion of emergent vegetation along lakeshores. We present the case of the lake Volkerak-Zoommeer in the Netherlands, a fresh-water lake created in 1987 after the enclosure of an estuary. Using an experimental area in which the water level could be manipulated, it was shown that partial summer drawdown of the shoreline created suitable conditions for germination and growth of tall emergent species (in particular Phragmites australis). Plant survival and growth depended on subsequent water-level fluctuations and grazing by waterbirds. Based on the experiment and empirical data, a model was developed to predict the effects of the water-level regime on potential reed bed development. The model was applied for four hydrological scenarios that have been considered for the water-level management of the lake.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Decoupling of As and Fe release to Bangladesh groundwater under reducing conditions. Part I: Evidence from sediment profiles

This study reexamines the notion that extensive As mobilization in anoxic groundwater of Bangladesh is intimately linked to the dissolution of Fe oxyhydroxides on the basis of analyses performed on a suite of freshly collected samples of aquifer material. Detailed sediment profiles extending to 40 to 70 m depth below the surface were obtained at six sites where local groundwater As concentrations were known to span a wide range. The sediment properties that were measured include (1) the proportion of Fe(II) in the Fe fraction leached in hot 1.2 N HCl, (2) diffuse spectral reflectance, and (3) magnetic susceptibility.In parallel with local concentrations of dissolved As ranging from <5 to 600 μg/L, Fe(II)/Fe ratios in shallow (gray) Holocene sands tended to gradually increase with depth from values of 0.3 to 0.5 to up to 0.9. In deeper (orange) aquifers of presumed Pleistocene age that were separated from shallow sands by a clay layer and contained <5 μg/L dissolved As, leachable Fe(II)/Fe ratios averaged ∼0.2. There was no consistent relation between sediment Fe(II)/Fe and dissolved Fe concentrations in groundwater in nearby wells. The reflectance measurements indicate a systematic linear relation (R² of 0.66; n = 151) between the first derivative transform of the reflectance at 520 nm and Fe(II)/Fe. The magnetic susceptibility of the shallow aquifer sands ranged from 200 to 3600 (x 10⁻⁹ m³/kg SI) and was linearly related (R² of 0.75; n = 29) to the concentrations of minerals that could be magnetically separated (0.03 to 0.79% dry weight). No systematic depth trends in magnetic susceptibility were observed within the shallow sands, although the susceptibility of deeper low-As aquifers was low (up to ∼200 × 10⁻⁹ m³/kg SI).This set of observations, complemented by incubation results described in a companion paper by van Geen et al. (this volume), suggests that the release of As is linked to the transformation of predominantly Fe (III) oxyhydroxide coatings on sand particles to Fe(II) or mixed Fe(II/III) solid phases with a flatter reflectance spectrum such as siderite, vivianite, or magnetite, without necessarily resulting in the release of Fe to groundwater. The very low As/Fe ratio of magnetically separated minerals compared to the As/Fe of bulk acid leachate (2 vs. 40 10⁻⁶, respectively) suggests that such a transformation could be accompanied by a significant redistribution of As to a mobilizable phase on the surface of aquifer particles.     

Particularités de la contamination crustale des phonolites : exemple du Velay oriental (Massif central)

Sr and Nd isotopic compositions of one trachyte, eight phonolites and five basalts have been measured. The isotopic characteristics of the trachyte can be explained by a combined assimilation–fractional crystallization process within an upper crustal magmatic chamber. Some phonolites display isotopic signatures identical to basalts, suggesting that they have been protected against any crustal assimilation during their formation. Some others have low Sr contents, whereas they are enriched in radiogenic Sr (0.70451<⁸⁷Sr/⁸⁶Sr_i<0.71192), and display basaltic ¹⁴³Nd/¹⁴⁴Nd ratios. Both observations could be explained by very strong alkali feldspar fractionation and by subsequent very low assimilation of surrounding rocks (between 0.3 and 4%) during intrusion. To cite this article: J.-M. Dautria et al., C. R. Geoscience 336 (2004).     

Crustal thickening in an extensional regime: application to the mid-Norwegian Vøring margin

The structure of the mid-Norwegian volcanic Vøring margin at the onset of the Maastrichtian–Paleocene extension phase reflects the cumulative effect of earlier consecutive rifting events. Lateral structural differences present on the margin at that time are a consequence of migration of the location of maximum extension in time between Norway and Greenland. The most important imprints (Moho depth, thermal structure) of these events on the lithosphere are incorporated in a numerical simulation of the final extension phase. We focus on a possible mechanism of formation of the Vøring Marginal High and address the relationship between spatial and temporal evolution of crustal thinning and thickening, uplift of the surface and strength of the lithosphere.It is found that the Vøring Basin formed the strongest part of the margin which explains why the Maastrichtian–Paleocene rift axis was not located here but instead jumped westward with respect to the earlier rift axes locations. The modeling study predicts that local crustal thickening during extension can be expected when large lateral thermal variations are present in the lithosphere at the onset of extension. Negative buoyancy induced by lateral temperature differences increases downwelling adjacent to the rifting zone; convergence of material at the particular part of the margin is mainly taken up by the lower crust. The model shows that during the final phase of extension, the crust in the Vøring Marginal High area was thickened and the surface uplifted. It is likely that this dynamic process and the effects of magmatic intrusions both acted in concert to form the Marginal High.     

The influence of particle composition on thorium scavenging in the NE Atlantic ocean (POMME experiment)

²³⁰Th, ²³²Th and ²³⁴Th were analyzed in sinking particles collected by moored and drifting sediment traps in the NE Atlantic Ocean (POMME experiment) in order to constrain the phase(s) carrying Th isotopes in the water column. It reveals a contrasted behaviour between ²³⁴Th and ²³⁰Th. ²³⁴Th is correlated to the particulate organic carbon suggesting that it is primarily scavenged by organic compounds in the surface waters. ²³⁰Th_xs is correlated with Mn, Ba and the lithogenic fraction that are enriched in small suspended particles and incorporated in the sinking particulate flux throughout the water column. The lack of correlation between ²³⁰Th_xs and CaCO₃ or biogenic silica (bSi) indicates that CaCO₃ and bSi are not responsible for ²³⁰Th scavenging in the deep waters of this oceanic region. ²³⁰Th is generally correlated with the lithogenic content of the trapped material but this correlation disappears in winter during strong atmospheric dust inputs suggesting that lithogenic matter is not directly responsible for ²³⁰Th scavenging in the deep waters or that sufficient time is required to achieve particle–solution equilibration. MnO₂ could be the prevalent ²³⁰Th_xs-bearing phase. The narrow range of K_{d_MnO2}^Th obtained for very contrasted oceanic environments supports a global control of ²³⁰Th_xs scavenging by MnO₂ and raises the possibility that the ²³⁰Th–²³¹Pa fractionation is controlled by the amount of colloidal MnO₂ in seawater.