Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Assessing Initial Conditions for Chloride Transport Across Low- permeability Argillaceous Rocks, Wellenberg, Switzerland

Information about fluid evolution and solute transport in a low-permeability metamorphic rock sequence has been obtained by comparing chloride concentrations and chlorine isotope ratios of pore water, groundwater, and fluid inclusions. The similarity of d37Cl values in fluid inclusions and groundwater suggests a closed-system evolution during the metamorphic overprint, and signatures established at this time appear to form the initial conditions for chloride transport after exhumation of the rock sequence.

Isolation and Characterization of a Microorganism from Groundwater that Reduces Arsenate

This is an investigation into the microbially mediated processes involved in the transformation of arsenic. With the recent change in the Federal Maximum Contaminant Level for arsenic in drinking water, an increasing amount of resources are being devoted to understanding the mechanisms involved in the movement of arsenic. Arsenic in drinking water typically comes from natural sources, but the triggers that result in increased release of arsenic from parent material are poorly understood. Knowledge of these processes is necessary in order to make sound engineering decisions regarding drinking water management practices. Recent years have brought forth the idea that bacteria play a significant role in arsenic cycling. Groundwater is a major source of potable water in this and many other countries. To date, no reports have been made indicating the presence and activity of arsenate reducing bacteria in groundwater settings, which may increase dissolved arsenic concentrations. This research was designed to address this question and has shown that these bacteria are present in Maine groundwater. Two Maine wells were sampled in order to culture resident bacteria that are capable of dissimilatory arsenate reduction. Samples were collected using anaerobic techniques fiom wells in Northport and Green Lake. These samples were amended with specific compounds to enrich the resident population of arsenate utilizing bacteria. These cultures were monitored over time to establish rates of arsenate reduction. Cultures fiom both sites exhibited arsenate reduction in initial enrichment cultures. Isolates obtained fiom the Green Lake enrichments, however, did not reduce arsenate. This indicates either that a symbiotic relationship was required for the observed arsenate reduction or that fast-growing fermentative organisms that could survive in high arsenate media were picked in the isolation procedure. The Northport cultures exhibited continued arsenate reduction after isolation and successive transfers into fiesh media. The cultured bacteria reduced the majority of 1 a arsenate solutions in less than one week, accompanied by a corresponding oxidation of lactate. The 16s rRNA fiom the isolate was arnplifled and sequenced. The results of the DNA sequence analysis indicate that the rRNA sequence of the bacteria isolated at the Northport site is unique. This means that this strain of bacteria has not been reported before. It is in the same taxonomic subgroup as two previously described arsenate respirers. The implications of this study are significant. The fact that resident bacteria are capable of reducing arsenate has implications for water management practices. Reduction of arsenate to arsenite increases the mobility of the compound, as well as the toxicity. An understanding of the activity of these types of organisms is necessary in order to understand the contribution they are making to arsenic concentrations in drinking water. The next step in this work would be to quantitj the actual loading of dissolved arsenic present in aquifers because of these organisms.

Fate and Transport of Naproxen in a Sandy Aquifer Material: Saturated Column Studies and Model Evaluation

Naproxen-C14H14O3 is a nonsteroidal anti-inflammatory drug which has been found at detectable concentrations in wastewater, surface water, and groundwater. Naproxen is relatively hydrophilic and is in anionic form at pH between 6 and 8. In this study, column experiments were performed using an unconsolidated aquifer material from an area near Barcelona (Spain) to assess transport and reaction mechanisms of Naproxen in the aquifer matrix under different pore water fluxes. Results were evaluated using HYDRUS-1D, which was used to estimate transport parameters. Batch sorption isotherms for Naproxen conformed with the linear model with a sorption coefficient of 0.42 (cm3 g−1), suggesting a low sorption affinity. Naproxen breakthrough curves (BTCs) measured in soil columns under steady-state, saturated water flow conditions displayed similar behavior, with no apparent hysteresis in sorption or dependence of retardation (R, 3.85-4.24) on pore water velocities. Soil sorption did not show any significant decrease for increasing flow rates, as observed from Naproxen recovery in the effluent. Sorption parameters estimated by the model suggest that Naproxen has a low sorption affinity to aquifer matrix. Most sorption of Naproxen occurred on the instantaneous sorption sites, with the kinetic sorption sites representing only about 10 to 40% of total sorption.

Vinasse application to sugar cane fields. Effect on the unsaturated zone and groundwater at Valle del Cauca (Colombia)

Extensive application of vinasse, a subproduct from sugar cane plantations for bioethanol production, is currently taking place as a source of nutrients that forms part of agricultural management in different agroclimatic regions. Liquid vinasse composition is characterised by high variability of organic compounds and major ions, acid pH (4.7), high TDS concentration (117,416–599,400 mg L− 1) and elevated EC (14,350–64,099 μS cm− 1). A large-scale sugar cane field application is taking place in Valle del Cauca (Colombia), where monitoring of soil, unsaturated zone and the aquifer underneath has been made since 2006 to evaluate possible impacts on three experimental plots. For this assessment, monitoring wells and piezometers were installed to determine groundwater flow and water samples were collected for chemical analysis. In the unsaturated zone, tensiometers were installed at different depths to determine flow patterns, while suction lysimeters were used for water sample chemical determinations. The findings show that in the sandy loam plot (Hacienda Real), the unsaturated zone is characterised by low water retention, showing a high transport capacity, while the other two plots of silty composition presented temporal saturation due to La Niña event (2010–2011). The strong La Niña effect on aquifer recharge which would dilute the infiltrated water during the monitoring period and, on the other hand dissolution of possible precipitated salts bringing them back into solution may occur. A slight increase in the concentration of major ions was observed in groundwater (~ 5% of TDS), which can be attributed to a combination of factors: vinasse dilution produced by water input and hydrochemical processes along with nutrient removal produced by sugar cane uptake. This fact may make the aquifer vulnerable to contamination.

Transport and fate of contaminants in the subsurface.

Mode of access: Internet.

Processes affecting subsurface transport of leaking underground tank fluids /

"CR 810052."

A plan for the comprehensive evaluation of the occurrence, transport, and fate of ground-water contaminants in the Lake Calumet area of southeast Chicago /

"HWRIC TN88-010."

Low-flow transport models for conservative and sorbed solutes, Uvas Creek, near Morgan Hill, California /

Bibliography: p. 73-74.

A Study of Shallow Groundwater Migration through Anthropogenic Fill along a Seawall

This report provides the findings and opinions of a historical document review, hydraulic balance calculation, and proposed additional study for a property that was historically used as a bulk petroleum storage and distribution facility. The property lies along the base, west, of a heavily vegetated bluff with a tidally influenced body of water west-adjacent to the property. The western portion of the property is bounded by a seawall spanning approximately 3,200 linear feet trending north-south. The seawall’s construction details are not known, save for a 225-foot section of driven sheet pile wall located within the northern portion of the property’s seawall. Due to the presence of petroleum hydrocarbons in soil and groundwater at the property, a cleanup action for the property will likely be overseen by the state regulatory agency. The property is currently undergoing remedial investigation in an effort to identify the lateral and vertical extent in which contaminants at the property have come to be located, also known as the “site” as defined by the Model Toxics Control Act (MTCA). The majority of the property bounded within the seawall area has been characterized; however, the shoreline sediments located immediately west-adjacent of the seawall have not been properly delineated. Identifying the bounds of the site to the west within sediment is pivotal for the purposes of the remedial investigation. Since the west-adjacent shoreline is so extensive, conducting a complete sediment sampling event along the entire shoreline would be cost-prohibitive due to analytical costs and logistical issues at the property. Because of the extensive nature of the shoreline, it would greatly benefit the client and project to focus sampling efforts at areas of greater risk for contaminants along the shoreline by identifying potential preferential pathways for contaminants to migrate off of the property and into adjacent shoreline sediments. The review of historical studies of the property yielded some useful information; however much of the findings included within the historical studies were lacking original raw data, therefore limiting the information obtained. The calculated hydraulic balance for the property yielded a relatively large surplus of recharge to groundwater after precipitation events, reinforcing the concept that contaminant have potentially historically, and currently, been migrating into the adjacent shoreline through preferential pathways along the seawall. Due to the limitation within the historical studies for the property as well as the groundwater recharge identified in the hydraulic balance, an additional study was proposed in an effort to provide additional aquifer characteristics along the seawall, and the ability to observe flow propagation at and proximate to the seawall in two-dimensions through time without the need to piece separate studies together. This proposed study includes a single simultaneous tidal study which comprises select monitoring points along the seawall. This report has identified the need for additional data that can be collected through available avenues for the property based upon the client’s desires and project needs. Ultimately, the proposed additional study is suggested based upon its relatively low capital investment, and ability meet the requirements relevant to the specific project needs and scope. Assuming preferential pathways are identified through the additional study proposed within this report, a representative and cost-effective sediment sampling plan can then be put in place in an effort to define the site.

Ultrastructure, aggregation-state, and crystal growth of biogenic nanocrystalline sphalerite and wurtzite

In this study, we investigated the size, submicrometer-scale structure, and aggregation state of ZnS formed by sulfate-reducing bacteria (SRB) in a SRB-dominated biofilm growing on degraded wood in cold (Tsimilar to8degreesC), circumneutral-pH (7.2-8.5) waters draining from an abandoned, carbonate-hosted Pb-Zn mine. High-resolution transmission electron microscope (HRTEM) data reveal that the earliest biologically induced precipitates are crystalline ZnS nanoparticles 1-5 nm in diameter. Although most nanocrystals have the sphalerite structure, nanocrystals of wurtzite are also present, consistent with a predicted size dependence for ZnS phase stability. Nearly all the nanocrystals are concentrated into 1-5 mum diameter spheroidal aggregates that display concentric banding patterns indicative of episodic precipitation and flocculation. Abundant disordered stacking sequences and faceted, porous crystal-aggregate morphologies are consistent with aggregation-driven growth of ZnS nanocrystals prior to and/or during spheroid formation. Spheroids are typically coated by organic polymers or associated with microbial cellular surfaces, and are concentrated roughly into layers within the biofilm. Size, shape, structure, degree of crystallinity, and polymer associations will all impact ZnS solubility, aggregation and coarsening behavior, transport in groundwater, and potential for deposition by sedimentation. Results presented here reveal nanometer- to micrometer-scale attributes of biologically induced ZnS formation likely to be relevant to sequestration via bacterial sulfate reduction (BSR) of other potential contaminant metal(loid)s, such as Pb2+, Cd2+, As3+ and Hg2+, into metal sulfides. The results highlight the importance of basic mineralogical information for accurate prediction and monitoring of long-term contaminant metal mobility and bioavailability in natural and constructed bioremediation systems. Our observations also provoke interesting questions regarding the role of size-dependent phase stability in biomineralization and provide new insights into the origin of submicrometer- to millimeter-scale petrographic features observed in low-temperature sedimentary sulfide ore deposits.