Electromagnetic induction prediction of soil salinity and groundwater properties in a Tunisian Saharan oasis

Abstract

Electromagnetic induction measurements (EM) were taken in a saline gypsiferous soil of the Saharan-climate Fatnassa oasis (Tunisia) to predict the electrical conductivity of saturated soil extract (ECe) and shallow groundwater properties (depth, Dgw, and electrical conductivity, ECgw) using various models. The soil profile was sampled at 0.2 m depth intervals to 1.2 m for physical and chemical analysis. The best input to predict the log-transformed soil salinity (lnECe) in surface (0–0.2 m) soil was the EMh/EMv ratio. For the 0–0.6 m soil depth interval, the performance of multiple linear regression (MLR) models to predict lnECe was weaker using data collected over various seasons and years (R _a ² = 0.66 and MSE = 0.083 dS m^-1) as compared to those collected during the same period (R _a ² = 0.97, MSE = 0.007 dS m^-1). For similar seasonal conditions, for the Dgw–EMv relationship, R ² was 0.88 and the MSE was 0.02 m for Dgw prediction. For a validation subset, the R ² was 0.85 and the MSE was 0.03 m. Soil salinity was predicted more accurately when groundwater properties were used instead of soil moisture with EM variables as input in the MLR.

Editor D. Koutsoyiannis; Associate editor K. Heal

Citation Bouksila, F., Persson, M., Bahri, A., and Berndtsson, R., 2012. Electromagnetic induction predictions of soil salinity and groundwater properties in a Tunisian Saharan oasis. Hydrological Sciences Journal, 57 (7), 1473–1486.     

Comparison of an Electrochemical and Luminescence‐Based Oxygen Measuring System for Use in the Biodegradability Testing According to Closed Bottle Test (OECD 301D)

The effects of persistent organic pollutants on the aquatic environment still can hardly be assessed. Therefore, ready biodegradability is a desired key characteristic of chemicals. Accordingly, biodegradation testing is of utmost importance and is a vital part within the risk assessment for organic chemicals. In the closed bottle test (CBT; OECD 301D) the degree of biological degradation of a chemical is assessed by monitoring the oxygen consumption so far with a Clark electrode as standard method. However, this method is time‐consuming, laborious and led to frequent fluctuations of the test results, which limited their validity. Due to these shortcomings, an optical system that enables contact‐free measuring of the oxygen concentration in a closed bottle, based on dynamic luminescence quenching, was examined. The test results were compared to those obtained with the “classical” electrochemical technique (Clark electrode). At first, a basic standard for the handling of the new instrument in the laboratory was developed. Then possible influencing factors, reliability, comparability, and reproducibility of the test results were investigated by running the CBTs in parallel with the electrode and optode method. The findings proved the new optode method to be unambiguously superior to the electrode technique. The frequency of fluctuations of the test results and time and effort necessary for the test run could considerably be reduced. The degradation kinetics of the test substances could be followed easily in the CBT by increasing the measuring points without much additional effort.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

Geology of the epithermal Ag–Au Huevos Verdes vein system and San José district, Deseado massif, Patagonia, Argentina

The San José district is located in the northwest part of the Deseado massif and hosts a number of epithermal Ag–Au quartz veins of intermediate sulfidation style, including the Huevos Verdes vein system. Veins are hosted by andesitic rocks of the Bajo Pobre Formation and locally by rhyodacitic pyroclastic rocks of the Chon Aike Formation. New ⁴⁰Ar/³⁹Ar constraints on the age of host rocks and mineralization define Late Jurassic ages of 151.3 ± 0.7 Ma to 144.7 ± 0.1 Ma for volcanic rocks of the Bajo Pobre Formation and of 147.6 ± 1.1 Ma for the Chon Aike Formation. Illite ages of the Huevos Verdes vein system of 140.8 ± 0.2 and 140.5 ± 0.3 Ma are 4 m.y. younger than the volcanic host rock unit. These age dates are among the youngest reported for Jurassic volcanism in the Deseado massif and correlate well with the regional context of magmatic and hydrothermal activity. The Huevos Verdes vein system has a strike length of 2,000 m, with several ore shoots along strike. The vein consists of a pre-ore stage and three main ore stages. Early barren quartz and chalcedony are followed by a mottled quartz stage of coarse saccharoidal quartz with irregular streaks and discontinuous bands of sulfide-rich material. The banded quartz–sulfide stage consists of sulfide-rich bands alternating with bands of quartz and bands of chlorite ± illite. Late-stage sulfide-rich veinlets are associated with kaolinite gangue. Ore minerals are argentite and electrum, together with pyrite, sphalerite, galena, chalcopyrite, minor bornite, covellite, and ruby silver. Wall rock alteration is characterized by narrow (< 3 m) halos of illite and illite/smectite next to veins, grading outward into propylitic alteration. Gangue minerals are dominantly massive quartz intergrown with minor to accessory adularia. Epidote, illite, illite/smectite, and, preferentially at deeper levels, Fe-chlorite gangue indicate near-neutral pH hydrothermal fluids at temperatures of >220°C. Kaolinite occurring with the late sulfide-rich veinlet stage indicates pH < 4 and a temperature of <200°C. The Huevos Verdes system has an overall strike of 325°, dipping on average 65° NE. The orientations of individual ore shoots are controlled by vein strike and intersecting north-northwest-striking faults. We propose a structural model for the time of mineralization of the San José district, consisting of a conjugate shear pair of sinistral north-northwest- and dextral west-northwest-striking faults that correspond to R and R′ in the Riedel shear model and that are related to master faults (M) of north-northeast-strike. Veins of 315° strike can be interpreted as nearly pure extensional fractures (T). Variations in vein strike predict an induced sinistral shear component for strike directions of >315°, whereas strike directions of <315° are predicted with an induced dextral strike–slip movement. The components of the structural model appear to be present on a regional scale and are not restricted to the San José district.     

Automated mapping of rock slope geometry,kinematics and stability with RSS-GIS

A GIS-implemented, deterministic approach for the automated spatial evaluation of geometrical and kinematical properties of rock slope terrains is presented. Based on spatially distributed directional information on planar geological fabrics and DEM-derived topographic attribute data, the internal geometry of rock slopes can be characterized on a grid cell basis. For such computations, different approaches for the analysis and regionalization of available structural directional information applicable in specific tectonic settings are demonstrated and implemented in a GIS environment. Simple kinematical testing procedures based on feasibility criteria can be conducted on a pixel basis to determine which failure mechanisms are likely to occur at particular terrain locations. In combination with hydraulic and strength data on geological discontinuities, scenario-based rock slope stability evaluations can be performed. For conceptual investigations on rock slope failure processes, a GIS-based specification tool for a 2-D distinct element code (UDEC) was designed to operate with the GIS-encoded spatially distributed rock slope data. The concepts of the proposed methodology for rock slope hazard assessments are demonstrated at three different test sites in Germany.     

Phase relations of phlogopite with magnesite from 4 to 8 GPa

To evaluate the stability of phlogopite in the presence of carbonate in the Earth’s mantle, we conducted a series of experiments in the KMAS–H₂O–CO₂ system. A mixture consisting of synthetic phlogopite (phl) and natural magnesite (mag) was prepared (phl₉₀-mag₁₀; wt%) and run at pressures from 4 to 8 GPa at temperatures ranging from 1,150 to 1,550°C. We bracketed the solidus between 1,200 and 1,250°C at pressures of 4, 5 and 6 GPa and between 1,150 and 1,200°C at a pressure of 7 GPa. Below the solidus, phlogopite coexists with magnesite, pyrope and a fluid. At the solidus, magnesite is the first phase to react out, and enstatite and olivine appear. Phlogopite melts over a temperature range of ~150°C. The amount of garnet increases above solidus from ~10 to ~30 modal% to higher pressures and temperatures. A dramatic change in the composition of quench phlogopite is observed with increasing pressure from similar to primary phlogopite at 4 GPa to hypersilicic at pressures ≥5 GPa. Relative to CO₂-free systems, the solidus is lowered such, that, if carbonation reactions and phlogopite metasomatism take place above a subducting slab in a very hot (Cascadia-type) subduction environment, phlogopite will melt at a pressure of ~7.5 GPa. In a cold (40 mWm⁻²) subcontinental lithospheric mantle, phlogopite is stable to a depth of 200 km in the presence of carbonate and can coexist with a fluid that becomes Si-rich with increasing pressure. Ascending kimberlitic melts that are produced at greater depths could react with peridotite at the base of the subcontinental lithospheric mantle, crystallizing phlogopite and carbonate at a depth of 180–200 km.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

Simulation of DNAPL infiltration and spreading behaviour in the saturated zone at varying flow velocities and alternating subsurface geometries

The influence of varying groundwater flow velocities on DNAPL infiltration and spreading behaviour was investigated by multiphase modelling using TMVOC and PetraSim. The multiphase models were calibrated by results of previously conducted laboratory experiments for the complete spatio-temporal range of the experiments. The small scale 2D scenario modelling was applied to qualify and quantify changes in position, architecture, geometry and dissolution of a TCE body in a fully saturated homogeneous sandy medium. The applied flow velocities ranging from 0.05 up to 40.00 m/day exhibited that the DNAPL TCE is affected even at the lowest flow velocity in its position, its size and its architecture. Additionally, several impermeable lenses with simple geometry were assumed in the model, to investigate the influence of stratified subsoil. In the experimental set-ups, the DNAPL body reacts more sensitive to the applied groundwater flow velocities than to the geometrical set-up of the scenarios. A possible consequence can be the transportation and displacement of a DNAPL pool due to natural or anthropogenic induced high groundwater flow velocities, as by Pump and Treat facilities, complicating site investigation process and planning of remediation activities.     

Investigating sediments and rock structures beneath a river using underwater ERT

This article presents hydrogeophysical investigations performed in a well-developed, long-term hydrogeological gypsum karst research site where subsurface evaporite dissolution has led to the subsidence of a river dam and an adjacent highway; both constructed on gypsum-containing rock, southeast of Basel, Switzerland. An observation system was set up to improve the protection of surface and subsurface water resources during remedial construction measures of the highway and in order to understand the processes, as well as the temporal evolution, of rock water interaction (flow and dissolution). However, no detailed hydrogeological information beneath the river could be derived from the previous investigations. To supplement the basic knowledge on this area, underwater Electrical Resistivity Tomography (ERT) measurements were conducted in the river bed upstream of the dam. The ERT-data are interpreted together with drill-core information and a conceptual 3D-Model of the area behind the dam and beneath the river. Results help to delineate weathered zones, associated faults and the thickness of sediment deposits behind the dam, as well as to locate voids within the local karst system. The combination of the ERT and modeling allows the optimization of future site-specific remedial construction measures.     

Skyline-based registration of 3D laser scans

Acquisition and registration of terrestrial 3D laser scans is a fundamental task in mapping and modeling of cities in three dimensions. To automate this task marker-free registration methods are required. Based on the existence of skyline features, this paper proposes a novel method. The skyline features are extracted from panoramic 3D scans and encoded as strings enabling the use of string matching for merging the scans. Initial results of the proposed method in the old city center of Bremen are pre-sented.