A mountain-scale thermal-hydrologic model for simulating fluid flow and heat transfer in unsaturated fractured rock

A multidimensional, mountain-scale, thermal-hydrologic (TH) numerical model is presented for investigating unsaturated flow behavior in response to decay heat from the proposed radioactive waste repository in the Yucca Mountain unsaturated zone (UZ), The model, consisting of both two-dimensional (2-D) and three-dimensional (3-D) representations of the UZ repository system, is based on the current repository design, drift layout, thermal loading scenario, and estimated current and future climate conditions. This mountain-scale TH model evaluates the coupled TH processes related to mountain-scale UZ flow. It also simulates the impact of radioactive waste heat release on the natural hydrogeological system, including heat-driven processes occurring near and far away from the emplacement tunnels or drifts. The model simulates predict thermally perturbed liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature elevations, as well as the changes in water flux driven by evaporation/condensation processes and drainage between drifts. These simulations provide insights into mountain-scale thermally perturbed flow fields under thermal loading conditions.     

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21 m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3 mx4 m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten alpha parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns.     

A stochastic multi-channel model for solute transport--analysis of tracer tests in fractured rock

Some of the basic assumptions of the advection-dispersion model (AD-model) are revisited. This model assumes a continuous mixing along the flowpath similar to Fickian diffusion. This implies that there is a constant dispersion length irrespective of observation distance. This is contrary to most field observations. The properties of an alternative model based on the assumption that individual water packages can retain their identity over long distances are investigated. The latter model is called the multi-channel model (MCh-model). Inherent in the latter model is that if the waters in the different pathways are collected and mixed, the "dispersion length" is proportional to distance. The conditions for when non-mixing between adjacent streams can be assumed are explored. The MCh- and AD-models are found to have very similar residence time distributions (RTD) for Peclet numbers larger than 3. A generalized relation between flowrate and residence time is developed, including the so-called cubic law and constant aperture assumptions. The two models extrapolate very differently when there is strong matrix interaction. The AD-model could severely underestimate the effluent concentration of a tracer pulse and overestimate the residence time. The conditions are explored for when in-filling particles in the fracture will not be equilibrated but will act as if there was seemingly a much larger flow wetted surface (FWS). It is found that for strongly sorbing tracers, relatively small particles can act in this way for systems and conditions that are typical of many tracer tests. The assumption that the tracer residence time found by cautiously injecting a small stream of traced water represents the residence time in the whole fracture is explored. It is found that the traced stream can potentially sample a much larger fraction of the fracture than the ratio between the traced flowrate and the total pumped flowrate. The MCh-model was used to simulate some recent tracer tests in what is assumed to be a single fracture at the Asp? Hard rock laboratory in Sweden. Non-sorbing tracers, HTO and Uranin were used to determine the mean residence time and its variance. Laboratory data on diffusion and sorption properties were used to "predict" the RTD of the sorbing tracers. At least 30 times larger FWS or 1000 times larger diffusion or sorption coefficients would be needed to explain the observed BTCs. Some possible reasons for such behavior are also explored.     

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.     

Analytical model of solute transport by unsteady unsaturated gravitational infiltration

Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils.     

Rapid transport from the surface to wells in fractured rock: a unique infiltration tracer experiment

A unique infiltration tracer experiment was performed whereby a fluorescent dye was applied to the land surface in an agricultural field, near Perth, Ontario, Canada, to simulate the transport of solutes to two pumped monitoring wells drilled into the granitic gneiss aquifer. This experiment, interpreted using the discrete-fracture capability of the numerical model HydroGeoSphere, showed that solute transport from the surface through thin soil (less than 2m) to wells in fractured bedrock can be extremely rapid (on the order of hours). Also, it was demonstrated that maximum concentrations of contaminants originating from the ground surface will not necessarily be the highest in the shallow aquifer horizon. These are important considerations for both private and government-owned drinking water systems that draw water from shallow fractured bedrock aquifers. This research illustrates the extreme importance of protecting drinking water at the source.     

Modeling tracer diffusion in fractured and unfractured, unsaturated, porous media

A proposed tracer diffusion test for the Exploratory Shaft Facility at Yucca Mountain, NV, is modeled. For the proposed test, a solution containing conservative tracers will be introduced into a borehole in the geologic medium of interest. The tracers will diffuse and advect from the saturated source region into the unsaturated matrix in the surrounding tuff. After some time, the borehole is to be overcored, and tracer concentrations in the fluid will be measured in the core as a function of distance from emplacement. The data will be used to evaluate diffusive behavior and to derive effective diffusion coefficients for the tracers in the specific tuff. Numerical simulations are used to study the effects of effective diffusion coefficient, porosity, saturation, and fracturing on tracer transport. Results are reported for numerical simulations of tests in the Topopah Spring Member and the Tuff of Calico Hills, which have significantly different porosities and saturations. The simulations make the following predictions: The spread of tracer during the test will be sensitive to the effective diffusion coefficient of the tracer. Tracer will diffuse farther in the Topopah Spring Member than in the Tuff of Calico Hills because of the former's lower porosity and saturation. Tracer transport by advection into the Topopah Spring Member will be greater than that into the Tuff of Calico Hills because of capillary effects. While advection will be a significant mechanism for tracer penetration into the Topopah Spring tuff, it will be less significant for tracer penetration into the Calico Hills tuff. The proximity of a single vertical fracture to the source region determines its effects on tracer transport, especially if the fracture diverts fluid flowing from the source region into the matrix.     

A numerical model for transient-hysteretic flow and solute transport in unsaturated porous media

A two-dimensional flow and transport model was developed for simulating transient water flow and nonreactive solute transport in heterogeneous, unsaturated porous media containing air and water. The model is composed of a unique combination of robust and accurate numerical algorithms for solving the Richards', Darcy flux, and advection-dispersion equations. The mixed form of Richards' equation is solved using a finite-element formulation and a modified Picard iteration scheme. Mass lumping is employed to improve solution convergence and stability behavior. The flow algorithm accounts for hysteresis in the pressure head-water content relationship. Darcy fluxes are approximated with a Galerkin and Petrov-Galerkin finite-element method developed for random heterogeneous porous media. The transport equation is solved using an Eulerian-Lagrangian method. A multi-step, fourth-order Runge-Kutta, reverse particle tracking technique and a quadratic-linear interpolation scheme are shown to be superior for determining the advective concentration. A Galerkin finite-element method is used for approximating the dispersive flux. The unsaturated flow and transport model was applied to a variety of rigorous problems and was found to produce accurate, mass-conserving solutions when compared to analytical solutions and published numerical results.     

Pilot-scale feasibility study for the stabilization of coal tailings via microbially induced calcite precipitation

Environmental Science and Pollution Research - Sustainable long-term solutions to managing tailings storage facilities (TSFs) are integral for mines to operate in a safe and environmentally...     

A remediation performance model for enhanced metabolic reductive dechlorination of chloroethenes in fractured clay till

A numerical model of metabolic reductive dechlorination is used to describe the performance of enhanced bioremediation in fractured clay till. The model is developed to simulate field observations of a full scale bioremediation scheme in a fractured clay till and thereby to assess remediation efficiency and timeframe. A relatively simple approach is used to link the fermentation of the electron donor soybean oil to the sequential dechlorination of trichloroethene (TCE) while considering redox conditions and the heterogeneous clay till system (clay till matrix, fractures and sand stringers). The model is tested on lab batch experiments and applied to describe sediment core samples from a TCE-contaminated site. Model simulations compare favorably to field observations and demonstrate that dechlorination may be limited to narrow bioactive zones in the clay matrix around fractures and sand stringers. Field scale simulations show that the injected donor is expected to be depleted after 5 years, and that without donor re-injection contaminant rebound will occur in the high permeability zones and the mass removal will stall at 18%. Long remediation timeframes, if dechlorination is limited to narrow bioactive zones, and the need for additional donor injections to maintain dechlorination activity may limit the efficiency of ERD in low-permeability media. Future work should address the dynamics of the bioactive zones, which is essential to understand for predictions of long term mass removal.     

A comparison of the low frequency electrical signatures of iron oxide versus calcite precipitation in granular zero valent iron columns

Geophysical methods have been proposed as technologies for non-invasively monitoring geochemical alteration in permeable reactive barriers (PRBs). We conducted column experiments to investigate the effect of mineralogy on the electrical signatures resulting from iron corrosion and mineral precipitation in Fe0 columns using (a) Na2SO4, and (b) NaHCO3 plus CaCl2 mixture, solutions. At the influent interface where the reactions were most severe, a contrasting time-lapse electrical response was observed between the two columns. Solid phase analysis confirmed the formation of corrosion halos and increased mineralogical complexity in the corroded sections of the columns compared to the minimal/non-corroded sections. We attribute the contrasting time-lapse signatures to the differences in the electrical properties of the mineral phases formed within the two columns. While newly precipitated/transformed polarizable and semi-conductive iron oxides (mostly magnetite and green rust) increase the polarization and conductivity of the sulfate column, the decrease of both parameters in the bicarbonate column is attributed to the precipitation of non-polarizable and non-conductive calcite. Our results show that precipitate mineralogy is an important factor influencing the electrical properties of the corroded iron cores and must be considered if electrical geophysical methods are to be developed to monitor PRB barrier corrosion processes in situ.     

Numerical simulations of pyrite oxidation and acid mine drainage in unsaturated waste rock piles

Numerical simulations of layered, sulphide-bearing unsaturated waste rock piles are presented to illustrate the effect of coupled processes on the generation of acid mine drainage (AMD). The conceptual 2D systems were simulated using the HYDRUS model for flow and the POLYMIN model for reactive transport. The simulations generated low-pH AMD which was buffered by sequential mineral dissolution and precipitation. Sulphide oxidation rates throughout the pile varied by about two orders of magnitude (0.004-0.4 kg m-3 year-1) due to small changes in moisture content and grain size. In the fine-grained layers, the high reactive surface area induced high oxidation rates, even though capillary forces kept the local moisture content relatively high. In waste rock piles with horizontal layers, most of the acidity discharged through vertical preferential flow channels while with inclined fine grained layers, capillary diversion channeled the AMD to the outer slope boundary, keeping the pile interior relatively dry. The simulation approach will be useful for helping evaluate design strategies for controlling AMD from waste rock.     

Effects of solution conditions on the precipitation of phosphate for recovery. A thermodynamic evaluation

To understand the effects of solution conditions on the precipitation of calcium phosphates from wastewater for recovery, a computer programme PHREEQC was employed to calculate the speciation and saturation-index (SI) with respect to hydroxyapatite of a chemically defined precipitation system, which contains phosphate of 1–200 mg P/l, with Ca/P molar ratios of one to 10 times of the stoichiometric calcium to phosphorus molar ratio of hydroxyapatite, at a pH range of 7.0–11.0. The results show that the SI is respectively the logarithmic function of the phosphate concentration and the calcium concentration, increasing with the increase of either of them; the SI is a polynomial function of the solution pH value and increases with its increase, and the effect of solution pH value is due to its influence on base uptake of the precipitation reaction and the speciation of phosphate and calcium ions; the SI is also a logarithmic function of the solution ionic strength but decreases with its increase; at the temperature range of 5–30 °C the SI increases linearly with solution temperature and the effect of temperature is also due to its influence on the speciation of phosphate and calcium ions.     

An active region model for capturing fractal flow patterns in unsaturated soils: model development

Preferential flow commonly observed in unsaturated soils allows rapid movement of solute from the soil surface or vadose zone to the groundwater, bypassing a significant volume of unsaturated soil and increasing the risk of groundwater contamination. A variety of evidence indicates that complex preferential patterns observed from fields are fractals. In this study, we developed a relatively simple active region model to incorporate the fractal flow pattern into the continuum approach. In the model, the flow domain is divided into active and inactive regions. Flow occurs preferentially in the active region (characterized by fractals), and inactive region is simply bypassed. A new constitutive relationship (the portion of the active region as a function of saturation) was derived. The validity of the proposed model is demonstrated by the consistency between field observations and the new constitutive relationship.     

Non-linear regression model for spatial variation in precipitation chemistry for South India

Chemical composition of rainwater changes from sea to inland under the influence of several major factors – topographic location of area, its distance from sea, annual rainfall. A model is developed here to quantify the variation in precipitation chemistry under the influence of inland distance and rainfall amount. Various sites in India categorized as ‘urban’, ‘suburban’ and ‘rural’ have been considered for model development. pH, HCO₃, NO₃ and Mg do not change much from coast to inland while, SO₄ and Ca change is subjected to local emissions. Cl and Na originate solely from sea salinity and are the chemistry parameters in the model.Non-linear multiple regressions performed for the various categories revealed that both rainfall amount and precipitation chemistry obeyed a power law reduction with distance from sea. Cl and Na decrease rapidly for the first 100 km distance from sea, then decrease marginally for the next 100 km, and later stabilize. Regression parameters estimated for different cases were found to be consistent (R² ～ 0.8). Variation in one of the parameters accounted for urbanization. Model was validated using data points from the southern peninsular region of the country. Estimates are found to be within 99.9% confidence interval.Finally, this relationship between the three parameters – rainfall amount, coastline distance, and concentration (in terms of Cl and Na) was validated with experiments conducted in a small experimental watershed in the south-west India. Chemistry estimated using the model was in good correlation with observed values with a relative error of ～5%. Monthly variation in the chemistry is predicted from a downscaling model and then compared with the observed data. Hence, the model developed for rain chemistry is useful in estimating the concentrations at different spatio-temporal scales and is especially applicable for south-west region of India.     

Responses of grassland ecosystem carbon fluxes to precipitation and their environmental factors in the Badain Jaran Desert

Environmental Science and Pollution Research - Studying the effects of precipitation on carbon exchange in grassland ecosystems is critical for revealing the mechanisms of the carbon cycle. In this...     

Geochemical evaluation of different groundwater-host rock systems for radioactive waste disposal

The geochemical suitability of a deep bedrock repository for radioactive waste disposal is determined by the composition of geomatrix and groundwater. Both influence radionuclide solubility, chemical buffer capacity and radionuclide retention. They also determine the chemical compatibility of waste forms, containers and backfill materials. Evaluation of different groundwater-host rock systems is performed by modeling the geochemical environments and the resulting radionuclide concentrations. In order to demonstrate the evaluation method, model calculations are applied to data sets available for various geological formations such as granite, clay and rocksalt. The saturation state of the groundwater-geomatrix system is found to be fundamental for the evaluation process. Hence, calculations are performed to determine if groundwater is in equilibrium with mineral phases of the geological formation. In addition, corrosion of waste forms in different groundwater is examined by means of reaction path modeling. The corrosion reactions change the solution compositions and pH, resulting in significant changes of radionuclide solubilities. The results demonstrate that geochemical modeling of saturation state and compatibility of the host formation environment with the radioactive waste proves to be a feasible tool for evaluation of various sites considered as deep underground repositories.     

A risk assessment tool for contaminated sites in low-permeability fractured media

A risk assessment tool for contaminated sites in low-permeability fractured media is developed, based on simple transient and steady-state analytical solutions. The discrete fracture (DF) tool, which explicitly accounts for the transport along fractures, covers different source geometries and history (including secondary sources) and can be applied to a wide range of compounds. The tool successfully simulates published data from short duration column and field experiments. The use for risk assessment is illustrated by three typical risk assessment case studies, involving pesticides, chlorinated solvents, benzene and MTBE. The model is compared with field data and with results from a simpler approach based on an Equivalent Porous Media (EPM). Risk assessment conclusions of the DF and EPM approaches are very different due to the early breakthrough, long term tailing, and lower attenuation due to degradation associated with fractured media. While the DF tool simulates the field data, it is difficult to conclude that the DF model is superior to an EPM model because of a lack of long term monitoring data. However, better agreement with existing field data by the DF model using observed physical fracture parameters favors the use of this model over the EPM model for risk assessments.     

A composite medium approximation for unsaturated flow in layered sediments

Saturated-unsaturated flow in strictly layered sediments proceeds via conductors in parallel in the direction parallel to bedding, and via resistors in series in the direction perpendicular to bedding. On sufficiently small scales of space and time, flow in such media will be subject to approximate capillary equilibrium locally, which provides a basis for approximating the effective hydraulic conductivity of a composite multi-layer medium in terms of the conductivities of the individual layers. Equations for the hydraulic conductivity tensor in "composite medium approximation" (COMA) are given in a coordinate system aligned with bedding. Hydraulic conductivity parallel to bedding is generally larger than in the perpendicular direction. The anisotropy depends on the spread of the conductivity distribution, and tends to increase for dryer conditions. The COMA model was implemented in a multi-phase flow simulator and tested by comparison with high-resolution simulations in which all layering heterogeneity is resolved explicitly. Under favorable conditions, COMA is found to accurately represent sub-grid scale flow and transport processes, providing a practical method for simulating field-scale flow and transport in layered media. The approximation improves when layers are thinner, and when flow rates are smaller.     

A model study of smoke-haze influence on clouds and warm precipitation formation in Indonesia 1997/1998

In the last few decades, fire and smoke-haze occurrence increased in Indonesia by intentionally set land clearing fires and higher fire susceptibility of disturbed forests. Particularly, during El Niño years with prolonged droughts in Indonesia, land clearing fires become uncontrolled wildfires and produce large amounts of gaseous and particulate emissions. This paper investigates the influence of smoke-haze aerosols from such fires on clouds and precipitation over Indonesia during the El Niño event 1997/1998 by numerical modelling. Warm precipitation formation in both layered and convective clouds is calculated dependent on the atmospheric aerosol concentration. In the smoke-haze affected regions of Indonesia, aerosol–cloud interactions induce events with both precipitation suppression and increase compared to a reference simulation without aerosol–cloud interactions. The effect of precipitation suppression is found to dominate with about 2/3 of all precipitation modification events pointing to a prolongation of smoke-haze episodes. The corresponding convective cloud top height of shallow clouds is increased whereas distinct lower deep convective cloud top heights are found. The remaining about 1/3 events are characterised by increased precipitation and cloud liquid water content, accompanied by lower convective cloud top heights of shallow clouds and higher deep convective clouds.