Self-potentials in partially saturated media: the importance of explicit modeling of electrode effects

Self-potential (SP) data are of interest to vadose zone hydrology because of their direct sensitivity to water flow and ionic transport. There is unfortunately little consensus in the literature about how to best model SP data under partially saturated conditions, and different approaches (often supported by one laboratory data set alone) have been proposed. We argue that this lack of agreement can largely be traced to electrode effects that have not been properly taken into account. A series of drainage and imbibition experiments were considered in which we found that previously proposed approaches to remove electrode effects were unlikely to provide adequate corrections. Instead, we explicitly modeled the electrode effects together with classical SP contributions using a flow and transport model. The simulated data agreed overall with the observed SP signals and allowed decomposing the different signal contributions to analyze them separately. After reviewing other published experimental data, we suggest that most of them include electrode effects that have not been properly taken into account. Our results suggest that previously presented SP theory works well when considering the modeling uncertainties presently associated with electrode effects. Additional work is warranted to not only develop suitable electrodes for laboratory experiments but also to assure that associated electrode effects that appear inevitable in longer term experiments are predictable, so that they can be incorporated into the modeling framework.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Long-term and long-range migration of radioactive fallout in a Karst system.

Mountainous areas are often covered by little evolved soils from which deposited radionuclides can potentially leak into the vadose zone. In the Swiss Jura mountains, we observed unusual isotopic ratios of nuclear weapon test (NWTs) fallout with an apparent loss of NWTs plutonium relative to &supl;³⁷Cs of Chernobyl origin in thinner soils. Here, we studied the karstic watershed of a vauclusian spring to determine the residence times of plutonium, ²⁴&supl;Am, and ⁹⁰Sr deposited by global fallout and their respective mobility in carbonaceous soils. The results show that ⁹⁰Sr is washed most efficiently from the watershed with a residence time of several hundred years. The estimated plutonium residence time is more than 10 times higher (in the range of 5000-10,000 years), and the ²⁴&supl;Am residence time is double that of plutonium. The spring water ²⁴&supl;Am/²³⁹+²⁴⁰Pu isotopic ratio is lower (0.12 - 0.28) than found in watershed soils (0.382 ± 0.077). Similar differences are found in aquatic mosses (²⁴&supl;Am/²³⁹+²⁴⁰Pu isotopic ratio 0.05-0.12), which are permanently submerged in spring waters. In contrast to plutonium, ⁹⁰Sr is leached from these mosses with 0.5M HCl, demonstrating that strontium is probably associated with calcium carbonate precipitations on the mosses. The higher plutonium to americium isotopic ratio found in the samples of spring water and mosses at the outlet of the karst shows that plutonium mobility is enhanced.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Groundwater recharge estimate at Alto Rio Grande - MG watershed

Springs are outcrops of aquifers surface, and the water cycle in this environment pass through the recharge, generally defined as the amount of water added to the aquifer, which may occur locally from rainwater infiltration. This study uses the Water Table Fluctuation (WTF) method to estimate the direct recharge and a groundwater balance to estimate the deep recharge on unconfined aquifers. The WTF method employs data of the aquifer water levels and its specific yield to estimate the direct recharge. The groundwater balance considers the direct recharge estimated by the WTF method, as the water input in the system and outputs as the base flow and deep recharge. The recharge was estimated at four areas at the watershed of Alto Rio Grande city, Minas Gerais (MG) state, in Brazil. The direct recharge estimate was 121.11; 64.62; 83.99; 152.46 (mm/year) for the L1, L2, M1 and M2 areas. The effect of the presence of forest in the recharge area can prevail over slope of relief, allowing more direct recharge, even in sources with steeper relief. The runoff from the springs in the study period exceeded the direct recharge, indicating a situation in which the saturated zone feeds the vadose zone. The annual flow was above the direct recharge pointing a situation of over exploitation of the aquifer, a non sustainable situation. The specific yield of the aquifers could also have been underestimated.

Petrography, Geochemistry and Diagenesis of Coral Deposits of Kavaratti and Minicoy Islands, Lakshadweep, India

In the present thesis the petrographic, geochemical and digenetic variability of the sediments in the islands of Kavaratti and Minicoy has been investigated .The beach profile studies show that in the Kavaratti lagoon beach, the slope is steeper in the southern and south central part than in the northern end’s marginal deposition is taking place in the northern end of the Kavaratti island, whereas a marginal erosion is observed at the southern end. In Minicoy the slope of the lagoon beach is gentle in the south and is slightly steeper at the northern part of the beach. The southern and northern beach sections show a marginal deposition. Based on the mineralogical and geochemical studies it is concluded that the main digenetic changes observed is the transformation of aragonite to LMC.The transformation takes place mainly in the vadose zone and is caused by the abundance of fresh water infiltration.

PALEOCLIMATIC CONTROLS ON STABLE OXYGEN AND CARBON ISOTOPES IN CALICHE OF THE ABO FORMATION (PERMIAN), SOUTH-CENTRAL NEW-MEXICO, USA

The stable oxygen and carbon isotopic composition of caliche in fluvial and supratidal rocks of the Abo Formation (Permian), south-central New Mexico, is controlled by palecoclimate and depositional environment. Fluvial caliche consists of low-Mg calcite nodules and vertically oriented tubules that display stage II texture. Micrite matrix support, brecciation, ooids/pisoliths, aveolar-septal texture, and peloids are common in the fluvial caliche and, along with red color and slickensides in the host shale, indicate pedogenesis in a well-oxidized vadose zone. In contrast, periodic waterlogging of the supratidal paleosols, probably due to high water table, is indicated by drab colors, carbonaceous flecks, horizontal rhizoliths, and the paucity of vadose textures in the stage II caliche nodules.Stable oxygen isotopes are similar in the fluvial and supratidal caliches and range from 21.6 to 30.5 parts per thousand (SMOW). The data exhibit a crude bimodality and delta-O-18 enrichment with a decrease in age (higher in the section). Consideration of these data in the context of delta-temperature relations suggests that 1) surface waters responsible for caliche formation increased in delta-O-18 (from roughly -8 to + 1 parts per thousand) over the 18 m.y. time interval that separated the lowest stratigraphic nodule horizon from the highest, 2) the increasing delta-O-18 values also reflect a warming trend (approximately 15-degrees to nearly 30-degrees-C) in the mean monthly temperature over this same time period, with perhaps an associated increase in Permian ocean temperatures, and 3) the significant variation in delta-O-18 from oldest to youngest caliche was probably enhanced by the amount effect, such that as the temperature increased, the amount of precipitation decreased, resulting in high delta-O-18 values.Caliches in the Abo are enriched in heavy carbon (-7.2 to -1.5 part per thousand PDB) compared to that of soil carbonate derived exclusively from C3 plants (-12 part per thousand PDB), and the supratidal caliches contain somewhat heavier carbon compared to the fluvial caliche. The delta-C-13 values for both environments increase with a decrease in caliche age. These results indicate that as the temperature increased and rainfall decreased with time, the level of C3 plant productivity apparently declined, allowing a greater influx of atmospheric CO2 into the soil. This can only occur when soil respiration rates are quite low or at very shallow depths (less than 10 cm), or both. Atmospheric CO2 seems to have invaded the supratidal soils to a somewhat greater extent than the fluvial soils.

Management Impacts on Soil Organic Matter of Tropical Soils

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A mathematical view of water table fluctuations in a shallow aquifer in Brazil

Detailed monitoring of the groundwater table can provide important data about both short- and long-term aquifer processes, including information useful for estimating recharge and facilitating groundwater modeling and remediation efforts. In this paper, we presents results of 4 years (2002 to 2005) of monitoring groundwater water levels in the Rio Claro Aquifer using observation wells drilled at the Rio Claro campus of São Paulo State University in Brazil. The data were used to follow natural periodic fluctuations in the water table, specifically those resulting from earth tides and seasonal recharge cycles. Statistical analyses included methods of time-series analysis using Fourier analysis, cross-correlation, and R/S analysis. Relationships could be established between rainfall and well recovery, as well as the persistence and degree of autocorrelation of the water table variations. We further used numerical solutions of the Richards equation to obtain estimates of the recharge rate and seasonable groundwater fluctuations. Seasonable soil moisture transit times through the vadose zone obtained with the numerical solution were very close to those obtained with the cross-correlation analysis. We also employed a little-used deep drainage boundary condition to obtain estimates of seasonable water table fluctuations, which were found to be consistent with observed transient groundwater levels during the period of study.

Ausbreitung von Nitroaromaten über das Sickerwasser der ehemaligen Espagit AG in Hallschlag

Die in den Jahren 1915-1923 erfolgte Produktion von Sprengstoffen und Vernichtung von Munition am Rüstungsaltstandort Hallschlag (Eifel) führte zu einer Nitroaromatenbelastung des Bodens, die bis heute nachgewiesen werden kann. Im Rahmen vorliegender Studie wurde das Sickerwasser außerhalb des Schadensherdes über einen Zeitraum von rund zwei Jahren mit Saugkerzensystemen beprobt und auf Nitroaromaten untersucht. Weiterhin wurde die klimatische Wasserbilanz für den Standort berechnet und Pegeldaten einer im Rahmen der Standortsicherung errichteten Sickerwassersammelleitung ausgewertet. Mit Hilfe dieses Datenmaterials konnte die mittlere Nitroaromatenkonzentration im Sickerwasser innerhalb und außerhalb des Schadensherdes ermittelt und die im hydrologischen Jahr 2004 ausgetragene Nitroaromatenmenge abgeschätzt werden.

The Influence of Fractures on Vulnerability to Pollution of a Carbonate Aquifer: Monte Conero (Italy)

The carbonate outcrops of the anticline of Monte Conero (Italy) were studied in order to characterize the geometry of the fractures and to establish their influence on the petrophysical properties (hydraulic conductivity) and on the vulnerability to pollution. The outcrops form an analog for a fractured aquifer and belong to the Maiolica Fm. and the Scaglia Rossa Fm. The geometrical properties of fractures such as orientation, length, spacing and aperture were collected and statistically analyzed. Five types of mechanical fractures were observed: veins, joints, stylolites, breccias and faults. The types of fractures are arranged in different sets and geometric assemblages which form fracture networks. In addition, the fractures were analyzed at the microscale using thin sections. The fracture age-relationships resulted similar to those observed at the outcrop scale, indicating that at least three geological episodes have occurred in Monte Conero. A conceptual model for fault development was based on the observations of veins and stylolites. The fracture sets were modelled by the code FracSim3D to generate fracture network models. The permeability of a breccia zone was estimated at microscale by and point counting and binary image methods, whereas at the outcrop scale with Oda’s method. Microstructure analysis revealed that only faults and breccias are potential pathways for fluid flow since all veins observed are filled with calcite. According this, three scenarios were designed to asses the vulnerability to pollution of the analogue aquifer: the first scenario considers the Monte Conero without fractures, second scenario with all observed systematic fractures and the third scenario with open veins, joints and faults/breccias. The fractures influence the carbonate aquifer by increasing its porosity and hydraulic conductivity. The vulnerability to pollution depends also on the presence of karst zones, detric zones and the material of the vadose zone.