Concepts and methods of quantitative assessment of risk to groundwater resources

This work presents a two-dimensional approach of risk assessment method based on the quantification of the probability of the occurrence of contaminant source terms, as well as the assessment of the resultant impacts. The risk is calculated using Monte Carlo simulation methods whereby synthetic contaminant source terms were generated to the same distribution as historically occurring pollution events or a priori potential probability distribution. The spatial and temporal distributions of the generated contaminant concentrations at pre-defined monitoring points within the aquifer were then simulated from repeated realisations using integrated mathematical models. The number of times when user defined ranges of concentration magnitudes were exceeded is quantified as risk. The utilities of the method were demonstrated using hypothetical scenarios, and the risk of pollution from a number of sources all occurring by chance together was evaluated. The results are presented in the form of charts and spatial maps. The generated risk maps show the risk of pollution at each observation borehole, as well as the trends within the study area. This capability to generate synthetic pollution events from numerous potential sources of pollution based on historical frequency of their occurrence proved to be a great asset to the method, and a large benefit over the contemporary methods.

Trichloroethylene fate and transport studies and biodegradation kinetics in the saturated zone

Permeable reactive barriers (PRB) are constructed from soil solid amendments to support the growth of bacteria that are capable of degrading organic contaminants. The objective of this study was to identify low-cost soil solid amendments that could retard the movement of trichloroethylene (TCE) while serving as long-lived carbon sources to foster its biodegradation in shallow groundwater through the use of a PRB. The natural amendments high in organic carbon content such as eucalyptus mulch, compost, wetland peat, organic humus were compared based on their geophysical characteristics, such as pHw, porosity and total organic carbon (TOC), and as well as TCE sorption potentials. The pHw values were within neutral range except for pine bark mulch and wetland peat. All other geophysical characteristics of the amendments showed suitability for use in a PRB. While the Freundlich model showed better fit for compost and pine bark mulch, the linear sorption model was adequate for eucalyptus mulch, wetland peat and Everglades muck within the concentration range studied (0.2-0.8 mg/L TCE). According to these results, two composts and eucalyptus mulch were selected for laboratory column experiments to evaluate their effectiveness at creating and maintaining conditions suitable for TCE anaerobic dechlorination. The columns were monitored for pH, ORP, TCE degradation, longevity of nutrients and soluble TOC to support TCE dechlorination. Native bacteria in the columns had the ability to convert TCE to DCEs; however, the inoculation with the TCE-degrading culture greatly increased the rate of biodegradation. This caused a significant increase in by-product concentration, mostly in the form of DCEs and VC followed by a slow degradation to ethylene. Of the tested amendments eucalyptus mulch was the most effective at supporting the TCE dechlorination. The experimental results of TCE sequential dechlorination took place in eucalyptus mulch and commercial compost from Savannah River Site columns were then simulated using the Hydrus-1D model. The simulations showed good fit with the experimental data. The results suggested that sorption and degradation were the dominant fate and transport mechanisms for TCE and DCEs in the column, supporting the use of these amendments in a permeable reactive barrier to remediate the TCE.

Experimental analysis of soil and dust stabilizers for controlling displacement of contaminated soils by wind tunnel experiments

During the remediation of burial grounds at the US Department of Energy's (DOE's) Hanford Site in Washington State, the dispersion of contaminated soil particles and dust is an issue that is faced by site workers on a daily basis. This contamination problem is even more of a concern when one takes into account the semi-arid characteristics of the region where the site is located. To mitigate this problem, workers at the site use a variety of engineered methods to minimize the dispersion of contaminated soil and dust (i.e. use of water and/or suppression agents that stabilizes the soil prior to soil excavation, segregation, and removal activities). A primary contributor to the dispersion of contaminated soil and dust is wind soil erosion. The erosion process occurs when the wind speed exceeds a certain threshold value which depends on a number of factors including wind force loading, particle size, surface soil moisture, and the geometry of the soil. Thus under these circumstances, the mobility of contaminated soil and generation and dispersion of particulate matter are significantly influenced by these parameters. This dependence of soil and dust movement on threshold shear velocity, fixative dilution and/or application rates, soil moisture content, and soil geometry were studied for Hanford's sandy soil through a series of wind tunnel experiments, laboratory experiments and theoretical analysis. In addition, the behavior of plutonium (Pu) powder contamination in the soil was studied by introducing a Pu simulant (cerium oxide). The results showed that soil dispersion and PM10 concentrations decreased with increasing soil moisture. Also, it was shown that the mobility of the soil was affected by increasing wind velocity. It was demonstrated that the use of fixative products greatly decreased the amount of soil and PM10 concentrations when exposed to varying wind conditions. In addition, it was shown that geometry of the soil sample affected the velocity profile and calculation of roughness surface coefficient when comparing round and flat soil samples. Finally, threshold shear velocities were calculated for soil with flat surface and their dependency on surface soil moisture was demonstrated. A theoretical framework was developed to explain these dependencies.

Investigating the origin of PM2.5 in Baltimore using highly time-resolved organic molecular markers measured during the Baltimore PM-supersite

The detailed organic composition of atmospheric fine particles with an aerodynamic diameter smaller than or equal to 2.5 micrometers (PM2.5) is an integral part of the knowledge needed in order to fully characterize its sources and transformation in the environment. For the study presented here, samples were collected at 3-hour intervals. This high time resolution allows gaining unique insights on the influence of short- and long-range transport phenomena, and dynamic atmospheric processes. A specially designed sequential sampler was deployed at the 2002-2003 Baltimore PM-Supersite to collect PM2.5 samples at a 3-hourly resolution for extended periods of consecutive days, during both summer and winter seasons. Established solvent-extraction and GC-MS techniques were used to extract and analyze the organic compounds in 119 samples from each season. Over 100 individual compounds were quantified in each sample. For primary organics, averaging the diurnal ambient concentrations over the sampled periods revealed ambient patterns that relate to diurnal emission patterns of major source classes. Several short-term releases of pollutants from local sources were detected, and local meteorological data was used to pinpoint possible source regions. Biogenic secondary organic compounds were detected as well, and possible mechanisms of formation were evaluated. The relationships between the observed continuous variations of the concentrations of selected organic markers and both the on-site meteorological measurements conducted parallel to the PM2.5 sampling, and the synoptic patterns of weather and wind conditions were also examined. Several one-to-two days episodes were identified from the sequential variation of the concentration observed for specific marker compounds and markers ratios. The influence of the meteorological events on the concentrations of the organic compounds during selected episodes was discussed. It was observed that during the summer, under conditions of pervasive influence of air masses originated from the west/northwest, some organic species displayed characteristics consistent with the measured PM2.5 being strongly influenced by the aged nature of these long-traveling background parcels. During the winter, intrusions from more regional air masses originating from the south and the southwest were more important.

An integrated flow and transport model to study the impact of mercury remediation strategies for East Fork Poplar Creek watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.