Mapping material distribution in a heterogeneous sand tank by image analysis

Water flow and solute transport through soils are strongly influenced by the spatial arrangement of soil materials with different hydraulic and chemical properties. Knowing the specific or statistical arrangement of these materials is considered as a key toward improved predictions of solute transport. Our aim was to obtain two-dimensional material maps from photographs of exposed profiles. We developed a segmentation and classification procedure and applied it to the images of a very heterogeneous sand tank, which was used for a series of flow and transport experiments. The segmentation was based on thresholds of soil color, estimated from local median gray values, and of soil texture, estimated from local coefficients of variation of gray values. Important steps were the correction of inhomogeneous illumination and reflection, and the incorporation of prior knowledge in filters used to extract the image features and to smooth the results morphologically. We could check and confirm the success of our mapping by comparing the estimated with the designed sand distribution in the tank. The resulting material map was used later as input to model flow and transport through the sand tank. Similar segmentation procedures may be applied to any high-density raster data, including photographs or spectral scans of field profiles.

Matrix Porewater In Crystalline Rocks: Extraction and Analysis

The chemical and isotopic characterization of porewater residing in the inter- and intragranular pore space of the low-permeability rock matrix is an important component with respect to the site characterization and safety assessment of potential host rocks for a radioactive waste disposal. The chemical and isotopic composition of porewater in such low permeability rocks has to be derived by indirect extraction techniques applied to naturally saturated rock material. In most of such indirect extraction techniques – especially in case of rocks of a porosity below about 2 vol.% – the original porewater concentrations are diluted and need to be back-calculated to in-situ concentrations. This requires a well-defined value for the connected porosity – accessible to different solutes under in-situ conditions. The derivation of such porosity values, as well as solute concentrations, is subject to various perturbations during drilling, core sampling, storage and experiments in the laboratory. The present study aims to demonstrate the feasibility of a variety of these techniques to charac-terize porewater and solute transport in crystalline rocks. The methods, which have been de-veloped during multiple porewater studies in crystalline environments, were applied on four core samples from the deep borehole DH-GAP04, drilled in the Kangerlussuaq area, Southwest Greenland, as part of the joint NWMO–Posiva–SKB Greenland Analogue Project (GAP). Potential artefacts that can influence the estimation of in situ porewater chemistry and isotopes, as well as their controls, are described in detail in this report, using specific examples from borehole DH-GAP04

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Pore water chemistry of sediment core GeoB15101-7 and water chemistry of CTD casts GeoB15101-1 and GeoB15101-6

Submarine brine lakes feature sharp and persistent concentration gradients between seawater and brine, though these should be smoothed out by free diffusion in open ocean settings. The anoxic Urania basin of the Eastern Mediterranean contains an ultra sulfidic, hypersaline brine of Messinian origin above a thick layer of suspended sediments. With a dual modeling approach we reconstruct its contemporary stratification by geochemical solute transport fundamentals, and show that thermal convection is required to maintain mixing in the brine and mud layer. The origin of the Urania basin stratification was dated to 1650 years before present, which may be linked to a major earthquake in the region. The persistence of the chemoclines may be key to the development of diverse and specialized microbial communities. Ongoing thermal convection in the fluid mud layer may have important, yet unresolved consequences for sedimentological and geochemical processes, also in similar environments.

Fractal Behaviour and Irregularity of Wetting Fronts in Heterogeneous Porous Media

The wetting front is the zone where water invades and advances into an initially dry porous material and it plays a crucial role in solute transport through the unsaturated zone. Water is an essential part of the physiological process of all plants. Through water, necessary minerals are moved from the roots to the parts of the plants that require them. Water moves chemicals from one part of the plant to another. It is also required for photosynthesis, for metabolism and for transpiration. The leaching of chemicals by wetting fronts is influenced by two major factors, namely: the irregularity of the fronts and heterogeneity in the distribution of chemicals, both of which have been described by using fractal techniques. Soil structure can significantly modify infiltration rates and flow pathways in soils. Relations between features of soil structure and features of infiltration could be elucidated from the velocities and the structure of wetting fronts. When rainwater falls onto soil, it doesn?t just pool on surfaces. Water ?or another fluid- acts differently on porous surfaces. If the surface is permeable (porous) it seeps down through layers of soil, filling that layer to capacity. Once that layer is filled, it moves down into the next layer. In sandy soil, water moves quickly, while it moves much slower through clay soil. The movement of water through soil layers is called the the wetting front. Our research concerns the motion of a liquid into an initially dry porous medium. Our work presents a theoretical framework for studying the physical interplay between a stationary wetting front of fractal dimension D with different porous materials. The aim was to model the mass geometry interplay by using the fractal dimension D of a stationary wetting front. The plane corresponding to the image is divided in several squares (the minimum correspond to the pixel size) of size length ". We acknowledge the help of Prof. M. García Velarde and the facilities offered by the Pluri-Disciplinary Institute of the Complutense University of Madrid. We also acknowledge the help of European Community under project Multi-scale complex fluid flows and interfacial phenomena (PITN-GA-2008-214919). Thanks are also due to ERCOFTAC (PELNoT, SIG 14)

Calibration of WAVE in irrigated maize: fallow vs. cover crops.

Nitrate leaching decreases crop available N and increases water contamination. Replacing fallow by cover crops (CC) is an alternative to reduce nitrate contamination, because it reduces overall drainage and soil mineral N accumulation. A study of the soil N and nitrate leaching was conducted during 5 years in a semi-arid irrigated agricultural area of Central Spain. Three treatments were studied during the intercropping period of maize (Zea mays L.): barley (Hordeum vulgare L.), vetch (Vicia villosa L.), and fallow. Cover crops, sown in October, were killed by glyphosate application in March, allowing direct seeding of maize in April. All treatments were irrigated and fertilised following the same procedure. Soil water content was measured using capacity probes. Soil Nmin accumulation was determined along the soil profile before sowing and after harvesting maize. Soil analysis was conducted at six depths every 0.20m in each plot in samples from 0 to 1.2-m depth. The mechanistic water balance model WAVE was applied in order to calculate drainage and plant growth of the different treatments, and apply them to the N balance. We evaluated the water balance of this model using the daily soil water content measurements of this field trial. A new Matlab version of the model was evaluated as well. In this new version improvements were made in the solute transport module and crop module. In addition, this new version is more compatible with external modules for data processing, inverse calibration and uncertainty analysis than the previous Fortran version. The model showed that drainage during the irrigated period was minimized in all treatments, because irrigation water was adjusted to crop needs, leading to nitrate accumulation on the upper layers after maize harvest. Then, during the intercrop period, most of the nitrate leaching occurred. Cover crops usually led to a shorter drainage period, lower drainage water amount and lower nitrate leaching than the treatment with fallow. These effects resulted in larger nitrate accumulation in the upper layers of the soil after CC treatments.

Cotransport of water by the Na+/glucose cotransporter

Water is transported across epithelial membranes in the absence of any hydrostatic or osmotic gradients. A prime example is the small intestine, where 10 liters of water are absorbed each day. Although water absorption is secondary to active solute transport, the coupling mechanism between solute and water flow is not understood. We have tested the hypothesis that water transport is directly linked to solute transport by cotransport proteins such as the brush border Na+/glucose cotransporter. The Na+/glucose cotransporter was expressed in Xenopus oocytes, and the changes in cell volume were measured under sugar-transporting and nontransporting conditions. We demonstrate that 260 water molecules are directly coupled to each sugar molecule transported and estimate that in the human intestine this accounts for 5 liters of water absorption per day. Other animal and plant cotransporters such as the Na+/Cl−/γ-aminobutyric acid, Na+/iodide and H+/amino acid transporters are also able to transport water and this suggests that cotransporters play an important role in water homeostasis.

High Expression of the Tonoplast Aquaporin ZmTIP1 in Epidermal and Conducting Tissues of Maize1

Aquaporins are integral membrane proteins of the tonoplast and the plasma membrane that facilitate the passage of water through these membranes. Because of their potentially important role in regulating water flow in plants, studies documenting aquaporin gene expression in specialized tissues involved in water and solute transport are important. We used in situ hybridization to examine the expression pattern of the tonoplast aquaporin ZmTIP1 in different organs of maize (Zea mays L.). This tonoplast water channel is highly expressed in the root epidermis, the root endodermis, the small parenchyma cells surrounding mature xylem vessels in the root and the stem, phloem companion cells and a ring of cells around the phloem strand in the stem and the leaf sheath, and the basal endosperm transfer cells in developing kernels. We postulate that the high level of expression of ZmTIP1 in these tissues facilitates rapid flow of water through the tonoplast to permit osmotic equilibration between the cytosol and the vacuolar content, and to permit rapid transcellular water flow through living cells when required.

Convergent and parallel activation of low-conductance potassium channels by calcium and cAMP-dependent protein kinase.

K+ channels, which have been linked to regulation of electrogenic solute transport as well as Ca2+ influx, represent a locus in hepatocytes for the concerted actions of hormones that employ Ca2+ and cAMP as intracellular messengers. Despite considerable study, the single-channel basis for synergistic effects of Ca2+ and cAMP on hepatocellular K+ conductance is not well understood. To address this question, patch-clamp recording techniques were applied to a model liver cell line, HTC hepatoma cells. Increasing the cytosolic Ca2+ concentration ([Ca2+]i) in HTC cells, either by activation of purinergic receptors with ATP or by inhibition of intracellular Ca2+ sequestration with thapsigargin, activated low-conductance (9-pS) K+ channels. Studies with excised membrane patches suggested that these channels were directly activated by Ca2+. Exposure of HTC cells to a permeant cAMP analog, 8-(4-chlorophenylthio)-cAMP, also activated 9-pS K+ channels but did not change [Ca2+]i. In excised membrane patches, cAMP-dependent protein kinase (the downstream effector of cAMP) activated K+ channels with conductance and selectivity identical to those of channels activated by Ca2+. In addition, cAMP-dependent protein kinase activated a distinct K+ channel type (5 pS). These data represent the differential regulation of low-conductance K+ channels by signaling pathways mediated by Ca2+ and cAMP. Moreover, since low-conductance Ca(2+)-activated K+ channels have been identified in a variety of cell types, these findings suggest that differential regulation of K+ channels by hormones with distinct signaling pathways may provide a mechanism for hormonal control of solute transport and Ca(2+)-dependent cellular functions in the liver as well as other nonexcitable tissues.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Migração de solutos em basalto fraturado: quantificação experimental em laboratório e validação matemática

A avaliação do risco a contaminação e a escolha de técnicas de remediação de poluentes em aquíferos fraturados depende da quantificação dos fenômenos envolvidos no transporte de solutos. A geometria da fratura, usualmente caracterizada pela abertura, é o principal parâmetro que indiretamente controla o transporte nos aquíferos fraturados. A simplificação mais comum desse problema é assumir que as fraturas são um par de placas planas e paralelas, isto é, com uma abertura constante. No entanto, por causa do limitado número de trabalhos experimentais, não está esclarecida a adequabilidade do uso de uma abertura constante para simular o transporte conservativo em fraturas do Aquífero Serra Geral (ASG), Brasil. O objetivo deste trabalho é avaliar a influência da abertura de uma fratura natural do Aquífero Serra Geral sob o transporte conservativo de solutos. Uma amostra natural de basalto fraturado foi usada em um experimento hidráulico e de transporte de um traçador conservativo (escala de laboratório). O campo de abertura foi medido usando a técnica avançada, de alta resolução e tridimensional, chamada microtomografia computadorizada de raios-X. A concentração de traçador medida foi utilizada para validar uma solução analítica unidimensional da Equação de Advecção-dispersão (ADE). O desemprenho do ajuste da ADE às curvas de passagem experimentais foi avaliado para quatro diferentes tipos de aberturas constantes. Os resultados mostraram que o escoamento de água e o transporte de contaminantes pode ocorrer através de fraturas micrométricas, ocasionando, eventualmente, a contaminação do ASG. A abertura de balanço de massa é a única que pode ser chamada propriamente de \"abertura equivalente\". O uso de aberturas constantes na ADE não permitiu representar completamente o formato das curvas de passagem porque o campo de velocidade não é uniforme e intrinsicamente bidimensional. Portanto, na simulação do transporte deve-se incorporar a heterogeneidade da abertura da fratura.

Dispersivity Determination Through a Modeling Approach From a Tracer Test Based on Total Br Concentration in Soil Samples

Determination of reliable solute transport parameters is an essential aspect for the characterization of the mechanisms and processes involved in solute transport (e.g., pesticides, fertilizers, contaminants) through the unsaturated zone. A rapid inexpensive method to estimate the dispersivity parameter at the field scale is presented herein. It is based on the quantification by the X-ray fluorescence solid-state technique of total bromine in soil, along with an inverse numerical modeling approach. The results show that this methodology is a good alternative to the classic Br− determination in soil water by ion chromatography. A good agreement between the observed and simulated total soil Br is reported. The results highlight the potential applicability of both combined techniques to infer readily solute transport parameters under field conditions.

Numerical investigation of seawater intrusion at Gooburrum, Bundaberg, Queensland, Australia

Seawater intrusion in coastal agricultural areas due to groundwater abstraction is a major environmental problem along the northeastern coast of Australia. Management options are being explored using numerical modelling, however, questions remain concerning the appropriate level of sophistication in models, choice of seaward boundary conditions, and how to accommodate heterogeneity and data uncertainty. The choice of seaward boundary condition is important since it affects the amount of salt transported into the aquifers and forms the focus of the present study. The impact of this boundary condition is illustrated for the seawater-intrusion problem in the Gooburrum aquifers, which occur within Tertiary sedimentary strata. A two-dimensional variable-density groundwater and solute-transport model was constructed using the computer code 2DFEMFAT (Cheng et al. 1998). The code was tested against an experiment for a steady-state freshwater-saltwater interface and against the Elder (Elder 1967) free-convection problem. Numerical simulations show that the imposition of the commonly-used equivalent hydrostatic freshwater heads, combined with a constant salt concentration at the seaward boundary, results in overestimated seawater intrusion in the lower Gooburrum aquifer. Since the imposition of this boundary condition allows water flow across the boundary, which subsequently takes salt into the aquifer, a careful check is essential to estimate whether too much mass of salt is introduced.

Modelling the fate of chromated copper arsenate in a sandy soil

A pulse of chromated copper arsenate (CCA, a timber preservative) was applied in irrigation water to an undisturbed field soil in a laboratory column. Concentrations of various elements in the leachate from the column were measured during the experiment. Also, the remnants within the soil were measured at the end of the experiment. The geochemical modelling package, PHREEQC-2, was used to simulate the experimental data. Processes included in the CCA transport modelling were advection, dispersion, non-specific adsorption (cation exchange) and specific adsorption by clay minerals and organic matter, as well as other possible chemical reactions such as precipitation/dissolution. The modelling effort highlighted the possible complexities in CCA transport and reaction experiments. For example, the uneven dosing of CCA as well as incomplete knowledge of the soil properties resulted in simulations that gave only partial, although reasonable, agreement with the experimental data. Both the experimental data and simulations show that As and Cu are strongly adsorbed and therefore, will mostly remain at the top of the soil profile, with a small proportion appearing in leachate. On the other hand, Cr is more mobile and thus it is present in the soil column leachate. Further simulations show that both the quantity of CCA added to the soil and the pH of the irrigation water will influence CCA transport. Simulations suggest that application of larger doses of CCA to the soil will result in higher leachate concentrations, especially for Cu and As. Irrigation water with a lower pH will dramatically increase leaching of Cu. These results indicate that acidic rainfall or significant accidental spillage of CCA will increase the risk of groundwater pollution.

Quantifying effects of soil heterogeneity on groundwater pollution in four sites in the USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction lysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater. The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.