Combined MODFLOW‐FRACTRAN Application to Assess Chlorinated Solvent Transport and Remediation in Fractured Sedimentary Rock

Detailed field investigations and numerical modeling were conducted to evaluate transport and fate of chlorinated solvent contamination in a fractured sedimentary bedrock aquifer (sandstone/siltstone/mudstone) at a Superfund site in central New Jersey. Field investigations provided information on the fractured rock system hydrogeology, including hydraulic gradients, bulk hydraulic conductivity, fracture network, and rock matrix, and on depth discrete contaminant distribution in fractures (via groundwater sampling) and matrix (via detailed subsampling of continuous cores). The numerical modeling endeavor involved application of both an equivalent porous media (EPM) model for flow and a discrete fracture network (DFN) model for transport. This combination of complementary models, informed by appropriate field data, allowed a quantitative representation of the conceptual site model (CSM) to assess relative importance of various processes, and to examine efficacy of remedial alternatives. Modeling progressed in two stages: first a large‐scale (20 km x 25 km domain) 3‐D EPM flow model (MODFLOW) was used to evaluate the bulk groundwater flow system and contaminant transport pathways under historic and current aquifer stress conditions and current stresses. Then, results of the flow model informed a 2‐D DFN transport model (FRACTRAN) to evaluate transport along a 1,000‐m flowpath from the source represented as a 2‐D vertical cross‐section. The combined model results were used to interpret and estimate the current and potential future extent of rock matrix and aqueous‐phase contaminant conditions and evaluate remedial strategies. Results of this study show strong effects of matrix diffusion and other processes on attenuating the plume such that future impacts on downgradient well fields under the hydraulic stresses modeled should be negligible. Results also showed futility of source remediation efforts in the fractured rock, and supported a technical impracticability (TI) waiver for the site. © 2013 Wiley Periodicals, Inc.     

Effects of hydrodynamic chromatography on colloid-facilitated migration of radionuclides in the fractured rock

The performance assessment of high level radioactive waste disposal has emphasized the role of colloids in the migration of radionuclides in the geosphere. Previous literature [Nagasaki S, Tanaka S, Suzuki A. Fast transport of colloidal particles through quartz-packed columns. J. Nucl. Sci. Technol. 1975;30(11):1136] indicates that owing to hydrodynamic chromatography the colloid velocity may not be equal to that of groundwater. Using hydrodynamic chromatography, this work investigates the effects of the size of colloidal particles on the radionuclide migration facilitated by colloids in a single fractured porous rock. Also, a methodology is proposed to develop a predictive model to assess transport within the fracture rock as well as various other phenomenological coefficients, particularly the size of colloidal particles. In addition, a fully developed concentration profile for non-reactive colloids in the fracture is developed to elucidate hydrodynamic chromatography of colloids in geological media. The external forces acting on colloidal particles hypothesized in the model proposed herein include inertial force, van der Waals attractive force, double layer force as well as gravitational force. The dispersion coefficient of colloids and the distribution coefficient for radionuclides with colloids are also considered as they pertain to the size of the colloid. In addition, the size distributions of colloids are utilized to investigate the effects of polydispersed colloids.     

Textural characterization of media composed of compacted pieces of cardboard and polyethylene using a gas tracer method

The aim of this work is the experimental determination of effective transport properties of porous media consisting of compacted pieces of cardboard and polyethylene (PE). The proposed method itself is more general and can be applied to many different materials and contexts. Three major transport properties were determined: porosity, tortuosity factor and permeability. Three parameters characterizing the media were varied over a wide range: the bulk density, the size of the elements entering the mix, and the proportion of cardboard and PE in the mix. The properties were measured by means of a specially designed experimental device based on miscible gas tracing. The porosity and tortuosity factor were simultaneously determined by parametric identification, based on the experimental sample output response to an inlet gas concentration step change compared to the results of a direct numerical model. Permeability was calculated in the standard way from the measurement of the pressure drop across the sample. The reproducibility of the measurements was very good. It was found that changing the material density of the medium significantly affects all three structural properties. When the bulk density is varied between 300 and 900 kg m(-3), the tortuosity factor varies in a range as large as 18-8 and the permeability decreases by a ratio of 2-3. The tortuosity factor shows unusual variation, characterized by a decrease when density is increased above 500 kg m(-3). The size of the elements does not significantly affect the structural properties of the medium in the range of parameters studied.     

Fracture Characterization Using Borehole Radar: Numerical Modeling

Forward modeling of borehole radar data for a series of synthetic discrete-fracture network (DFN) models provides a conceptual framework for interpreting experimental field data at fractured rock sites. A finite-difference time-domain (FDTD) radar wave propagation model was developed for this purpose. Synthetic examples demonstrate the utility of single-hole reflection-mode and cross-hole transmission-mode borehole radar for (1) identification of fracture location and orientation, and (2) identification of fracture pore-fluid properties, which might change as a result of tracer tests or flooding exercises in support of resource development or site remediation. A two-dimensional, synthetic DFN was generated statistically based on hypothetical distributions of fracture length, orientation, aperture, permeability, and inter-connectivity. The DFN includes a zone of permeable fractures embedded within a network of lower-permeability fractures and a low-permeability rock matrix. We modeled the unconnected and non-permeable fractures as being filled with freshwater. To simulate tracer experiments, contaminant releases, or engineered-remediation processes, we considered alternately the inter-connected, permeable fractures to be filled with freshwater air, or saline water (tracer). Synthetic radar data sets for both single-hole reflection and cross-hole transmission modes were generated. The features in synthetic radargrams were then examined and compared to the DFN model to evaluate the likelihood of identifying fracture location, orientation, and pore fluid in field situations. This comparison demonstrates that (1) the replacement of freshwater with saline water in permeable fractures generally increases the amplitude of reflections from permeable, connected fractures; and (2) in general, radar reflection-mode data contains more information about fracture properties than transmission-mode data.     

Contribution of 3-D time-lapse ERT to the study of leachate recirculation in a landfill

Leachate recirculation is a key process in the operation of municipal waste landfills as bioreactors. It aims at increasing the moisture content to optimise the biodegradation. Because waste is a very heterogeneous and anisotropic porous media, the geometry of the leachate plume recirculation is difficult to delineate from the surface at the scale of the bioreactor site. In this study, 3-D time-lapse electrical resistivity tomography (ERT) was used to obtain useful information for understanding leachate recirculation hydrodynamics. The ERT inversion methodology and the electrode arrays were optimised using numerical modelling simulating a 3-D leachate injection scenario. Time-lapse ERT was subsequently applied at the field scale during an experimental injection. We compared ERT images with injected volumes to evaluate the sensitivity of time-lapse ERT to delineate the plume migration. The results show that time-lapse ERT can accomplish the following: (i) accurately locate the injection plume, delineating its depth and lateral extension; (ii) be used to estimate some hydraulic properties of waste.     

Comparison of Microscopic and Macroscopic Modeling Approaches for Subsurface Contaminant Transport

The practice of contaminant transport and remediation has shown significant progress in recent years. However, despite the significant progress made, remediation efforts are often delayed by extremely long breakthrough curve tails that render efforts to bring the level of contaminants below the regulatory standards inefficient. One hypothesis is that these long tails are due to the reservoir-like slow diffusive processes in soil micropore zones. This study compares the effects of micropores at macroscopic and microscopic levels and establishes a link between these approaches for validation and calibration purposes. The link between macroscopic and microscopic levels is established through comparisons and testing of the two models while incorporating appropriate scale and boundary effects. Despite the differences in conceptual approaches and simulation time, the two approaches rendered meaningful results. The link helps forecast the effects of micropore zone transport processes in the subsurface efficiently and thus allows development of numerical tools that could contribute towards more efficient remediation design.     

Numerical modelling of multiphase flow and transport processes in landfills.

Waste material in municipal landfills can be described as heterogeneous porous media, where flow and transport processes of gases and liquids are combined with local material degradation. This paper deals with the basic formulation of a multiphase flow and transport model applicable to the numerical analysis of coupled transport and reaction processes inside landfills. The transport model treats landfills within the framework of continuum mechanics, where flow and transport processes are described on a macroscopic level. The composition of organic and inorganic matter in the solid phase and its degradation are modelled on a microscopic scale. The degradation model captures the different reaction schemes of various microbial activities. Subsequently, transport and reaction processes have to be coupled, since emissions at the surface and from the drainage layer depend on the flow of leachate and gas, the transport of various substances and heat, and the biodegradation of organic matter. The theoretical considerations presented here are fundamental to the development of numerical models for the simulation of multiphase flow and transport processes inside landfills coupled with biochemical reactions and heat generation. The implicit modelling of leachate and gas flows including growth and decay of micro-organisms are innovative contributions to landfill modelling     

Two‐dimensional transport of lactate‐modified nanoscale iron particles in porous media

This study investigates the two‐dimensional transport of nanoscale iron particles (NIP) and lactate‐modified NIP (LMNIP) in homogeneous and heterogeneous porous media under typical pressurized groundwater flow conditions. A two‐dimensional bench‐scale test setup was developed and a series of experiments was conducted simulating homogeneous sand profile and two‐layer profile with two different sands. NIP and LMNIP at a concentration of 4 g/L were prepared in electrolyte simulating groundwater conditions and were injected at the inlet of the test setup under different pressure gradients (0.5. 0.8, 1, and 2 pounds per square inch). During the testing, effluent was collected and its volume and nanoiron concentrations were measured. At the end of the testing, soil cores were obtained at different distances from the inlet and were used to measure nanoiron concentrations and magnetic susceptibility values. Results showed that the transport of NIP and LMNIP was enhanced by increased pressure gradient. LMNIP transport occurred more uniformly as compared to bare NIP. The iron concentrations decreased with distance from the inlet to the outlet and increased from the top to the bottom of the test cell. The data indicate that, as the particles were transported, they underwent aggregation and sedimentation, which resulted in the observed non‐uniform spatial distribution of iron. The NIP and LMNIP transported through the high‐porosity and high‐permeability soil layer in the heterogeneous soil profile, implying that the transport occurred predominantly along the path of least resistance for the flow. Magnetic susceptibility values are found to have good correlation with the iron content in the soil and are helpful to characterize the transport of NIP and LMNIP. Overall, this study shows that the non‐uniform distribution of NIP and LMNIP occurs under two‐dimensional transport conditions and the soil heterogeneities can significantly impact the transport of NIP and LMNIP. The design of field delivery systems should consider such conditions and optimize the pressurized injection systems. © 2011 Wiley Periodicals, Inc.     

Application of Biobarriers for Groundwater Containment at Fractured Bedrock Sites

Biological barriers are a beneficial application of biofilms that aim at reducing the hydraulic conductivity (K) in geological formations. Several studies have shown the potential benefits of creating such barriers either by stimulating the indigenous microbial community (biostimulation) or injecting bacteria (bioaugmentation). For example, laboratory experiments show that groundwater microorganisms attached to a ceramic surface and generated a biofilm as thick as 1,100 μm. In a limestone fracture, this bacterial community clogged a single fracture up to 99.2 percent within 22 days. At the field scale, applications in porous aquifers led to a five‐fold decrease in K after 2.5 days of biostimulation, and a bioaugmentation with a starved, adapted bacterial culture decreased K by 99.4 percent. One promising development of the biobarrier concept is a field application at a fractured bedrock site. Using a multidisciplinary approach and focusing on a well‐characterized fracture system, a field trial was undertaken in Southern Ontario to measure the extent of bioclogging and the stability over time. This article focuses on the literature pertinent to the preparation of this field trial and presents the innovative approach selected to monitor the bioclogging in such a challenging environment.     

Stochastic modelling of landfill processes incorporating waste heterogeneity and data uncertainty

A landfill is a very complex heterogeneous environment and as such it presents many modelling challenges. Attempts to develop models that reproduce these complexities generally involve the use of large numbers of spatially dependent parameters that cannot be properly characterised in the face of data uncertainty. An alternative method is presented, which couples a simplified microbial degradation model with a stochastic hydrological and contaminant transport model. This provides a framework for incorporating the complex effects of spatial heterogeneity within the landfill in a simplified manner, along with other key variables. A methodology for handling data uncertainty is also integrated into the model structure. Illustrative examples of the model's output are presented to demonstrate effects of data uncertainty on leachate composition and gas volume prediction.     

Performing contaminant mass balances for remedy assessments

Contaminant mass‐balance assessments are useful tools to help quantify various mass transport and removal mechanisms that may be active in a remedial system setting. This article presents the basics of performing a mass balance and illustrates the utility of using the information derived to support project management decisions. It is important to understand the partitioning of contaminant mass into various environmental media and physical forms, as well as the relationships among the partitions. Contaminant partitioning tends toward an equilibrium state, so natural or engineered mass transfer into or out of one partition will affect the others. Mass balances are exercises that quantify, to the extent possible, the contaminant mass in the various environmental partitions and the transfer and transformation processes that affect contaminant distribution. Understanding mass partitioning and transfer mechanisms helps remediation practitioners to engineer and optimize those mechanisms that contribute to risk reduction at a contaminated site. Such knowledge can inform risk managers when natural mechanisms may dominate engineered approaches and help identify uncertainties in contaminant fate and transport. © 2009 Wiley Periodicals, Inc.     

Visualization of the Malleability of the Rubber Core of Rubber Particles from Parthenium argentatum Gray and Other Rubber-Producing Species Under Extremely Cold Temperatures

Rubber particles from Parthenium argentatum Gray (guayule) were frozen in liquid nitrogen (–196°C), fractured, and visualized using cryo-scanning electron microscopy. We observed that the rubber polymer core of the rubber particles was still malleable at this extremely cold temperature, and the core stretched substantially during separation of the fracture planes. This malleability was observed in situ in tissue sections, as well as in purified rubber particles, and was found to be independent of purification procedure, guayule line, tissue age, or season. The malleability or stretching phenomenon suggests that P. argentatum rubber has some unique properties because rubber particles from Hevea brasiliensis Müll. Arg. and Ficus elastica Roxb. were brittle at this temperature, fractured cleanly, or showed only tiny threads of material pulling out of the core.     

Waste ecocompatibility in storage and reuse scenarios: global methodology and detailed presentation of the impact study on the recipient environments

In 1995, the ADEME launched a research program called "Waste Ecocompatibility" in order to define a reliable methodology for measuring the impact of waste in storage or reuse scenarios. The French concept of "Ecocompatibility" is defined as the situation where the pollutant flux from waste disposed of or used in specified conditions is compatible with the environmental acceptance of the receiving environments. The chief feature of this definition is to integrate the evaluation of the three following terms: pollutants emission from the waste, transport of the pollutants from the waste to the receptor cells and the environmental acceptance of the receiving environments. The "Waste Ecocompatibility" program consisted of a literature survey and an experimental part. The literature study aimed to determine factors and waste characteristics to be considered for a reliable ecocompatility assessment, to provide an overview of the available tools for measuring those factors and characteristics and to propose a first approach of the methodology. In the framework of the experimental program, this approach was then applied to three theoretical scenarios to validate the laboratory tools (comparative study of laboratory and field results) and to calibrate the global methodology. This paper deals with the results of the experimental program concerning the impact study on receiving environments: impact on plants and microorganisms living in soil, impacts on soil fauna and aquatic fauna. In other papers we intend to present the operational methodology for the assessment of waste ecocompatibility. It includes bio-assays at laboratory scale (microcosms), pilot scale (mesocosms) and in situ experiments (experimental prairie). To limit the use of in situ experiments other research works are necessary to validate bio-assays at laboratory or pilot scale.     

Nanotechnology for site remediation

As a remediation tool, nanotechnology holds promise for cleaning up hazardous waste sites cost‐effectively and addressing challenging site conditions, such as the presence of dense nonaqueous phase liquids (DNAPLs). Some nanoparticles, such as nanoscale zero‐valent iron (nZVI) are already in use in full‐scale projects with encouraging success. Ongoing research at the bench and pilot scale is investigating particles such as self‐assembled monolayers on mesoporous supports (SAMMS™), dendrimers, carbon nanotubes, and metalloporphyrinogens to determine how to apply their unique chemical and physical properties for full‐scale remediation. There are many unanswered questions regarding nanotechnology. Further research is needed to understand the fate and transport of free nanoparticles in the environment, whether they are persistent, and whether they have toxicological effects on biological systems. In October 2008, the U.S. Environmental Protection Agency's Office of Superfund Remediation and Technology Innovation (OSRTI) prepared a fact sheet entitled “Nanotechnology for Site Remediation,” and an accompanying list of contaminated sites where nanotechnology has been tested. The fact sheet contains information that may assist site project managers in understanding the potential applications of this group of technologies. This article provides a synopsis of the US EPA fact sheet, available at http://clu‐in.org/542F08009 , and includes background information on nanotechnology; its use in site remediation; issues related to fate, transport, and toxicity; and a discussion of performance and cost data for field tests. The site list is available at http://clu‐in.org/products/nanozvi . © 2008 Wiley Periodicals, Inc.     

Migration of radionuclides in porous rock in the presence of colloids: effects of kinetic interactions

This work investigates the colloid-facilitated migration of radionuclides with radioactive decay in porous media. The sorption processes for radionuclides with both the solid matrix and colloids are treated as equilibrium or nonequilibrium. An analytical solution is obtained from a simplified linear equilibrium interaction mechanism. In addition, the adsorption processes for radionuclides with colloids and porous rock can be assumed as nonequilibrium and modeled by the linear kinetic adsorption. The numerical method is employed to solve the coupled colloid and radionuclide transport equations under nonequilibrium sorption assumption. Moreover, the reaction rates of the adsorption processes for radionuclides with the solid matrix and colloids affect the transport characteristics of radionuclides. The fast reaction rate of radionuclides with colloids causes a higher concentration of radionuclides adsorbed on colloids in a dispersed phase and enlarges acceleration caused by colloids. However, the fast reaction rate for radionuclides with solid matrix increases the retardation effect caused by the solid matrix. This work developed a predictive model for the transport of colloid-facilitated radionuclides in porous medium and to assess the importance of various phenomenological coefficients, particularly parameters for the adsorption interactions.     

Borehole‐scale testing of matrix diffusion for contaminated‐rock aquifers

A new method was developed to assess the effect of matrix diffusion on contaminant transport and remediation of groundwater in fractured rock. This method utilizes monitoring wells constructed of open boreholes in the fractured rock to conduct backward diffusion experiments on chlorinated volatile organic compounds (CVOCs) in groundwater. The experiments are performed on relatively unfractured zones (called test zones) of the open boreholes over short intervals (approximately 1 meter) by physical isolation using straddle packers. The test zones were identified with a combination of borehole geophysical logging and chemical profiling of CVOCs with passive samplers in the open boreholes. To confirm the test zones are within inactive flow zones, they are subjected to a series of hydraulic tests. Afterward, the test zones are air sparged with argon to volatilize the CVOCs from aqueous to air phase. Backward diffusion is then measured by periodic passive‐sampling of water in the test zone to identify rebound. The passive (nonhydraulically stressed) sampling negates the need to extract water and potentially dewater the test zone. The authors also monitor active flowing zones of the borehole to assess trends in concentrations in other parts of the fractured rock by purge and passive sampling methods. The testing was performed at the former Pease Air Force Base (PAFB) in Portsmouth, New Hampshire. Bedrock at the former PAFB consists of fractured metasedimentary rocks where the authors investigated back diffusion of cis‐1,2‐dichloroethylene (cis‐1,2‐DCE), a CVOC. Postsparging concentrations of cis‐1,2‐DCE showed initial rebounding followed by declines, excluding an episodic spike in concentrations from a groundwater recharge event. The authors theorize that there are three processes that controlled concentration responses in the test zones postsparging. First, the limited back diffusion of CVOCs from a halo or thin zone of rock around the borehole contributes to the initial rebounding. Second, aerobic degradation of cis‐1,2‐DCE occurred causing declines in concentrations in the test zone. Third, microflow from microfractures contributed to the episodic spike in concentrations following the groundwater recharge event. In active flow zones, the latter two processes are not measurable due to equilibration from groundwater transport between the borehole and active flowing fractures.     

Technical and environmental long-term properties of industrial residues--summary of field and laboratory investigations

In Sweden, use of industrial residues is still hindered by concern for their long-term properties. A three-year research project was therefore initiated aiming to (1) identify the crucial processes of ageing related to the usefulness of residues in roads; (2) investigate the consequences of these processes for technical and environmental properties of the residues, and (3) propose a method for accelerated ageing to predict the long-term properties. This paper gives an overview of the project methodology, a summary of the test results and references to papers where further details are given.The project, running through 2006-2008, compared naturally aged samples of two residues used as sub-bases in existing asphalt paved roads with samples of fresh residues from producers’ piles. Steel slag of electric arc furnace (EAF) type and municipal solid waste incinerator (MSWI) bottom ash were chosen. The samples were thoroughly characterised in order to identify which ageing processes had been crucial.The results showed that:

-

Bottom ash from the pavement edge was more aged than bottom ash from the road centre. However, no difference in pH was found, instead the differences were caused by differences in water exposure.

-

Steel slag from the pavement edge showed traces of carbonation and leaching processes, whereas slag from the road centre was identical to fresh slag.

-

Water exposure to the subbase materials after ten years in an asphalt paved road was calculated to less than 0.1-0.5 litres per kg.

-

Ageing reactions in steel slag and MSWI bottom ash, ready for use, were too small to be verified by laboratory measurement of deformation properties under loaded conditions.

An accelerated ageing test for steel slag was set up to achieve the carbonation (decrease in pH) and leaching that was observed in the pavement edge material.An accelerated ageing test for bottom ash was set up to achieve the pozzolan reactions that were observed in SEM analyses of in situ specimens.It is recommended to use uncrushed particles when properties of aged material are studied, in order to preserve the original particle surfaces.     

The Development of a GIS-based Inventory of Standing Waters in Great Britain together with a Risk-based Prioritisation Protocol

An inventory of standing waters (freshwater lakes and lochs) wasderived from Ordnance Survey digital map data at a scale of 1:50 000 and represents the most comprehensive survey of its kind for Great Britain. The inventory includes 43 738 water bodies in England, Scotland, Wales and the Isle of Man and contains basic physical data such as location, surface area, perimeter and altitude. Catchment areas were computed for water bodies with a surface area larger than 1 ha from a digital terrain model (DTM) using customised routines in a geographical information system (GIS). The resulting polygons were then used to derive catchment-related information from a variety of national datasets including population density, livestock density, land cover, solid and drift geology, meteorological data, freshwater sensitivity status, acid deposition and conservation status. Using data derived from the inventory a risk-based prioritisation protocol was developed to identify standing waters at risk of harm from acidification and eutrophication. This information is required by the Environment Agency, Scottish Environmental Protection Agency and the U.K. statutory conservation bodies to co-ordinate actions and monitor change under international, European and national legislation.