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.

Integrated assessment of water conservation practices for sustainable management strategies

Miami-Dade County implemented a series of water conservation programs, which included rebate/exchange incentives to encourage the use of high efficiency aerators (AR), showerheads (SH), toilets (HET) and clothes washers (HEW), to respond to the environmental sustainability issue in urban areas. This study first used panel data analysis of water consumption to evaluate the performance and actual water savings of individual programs. Integrated water demand model has also been developed for incorporating property’s physical characteristics into the water consumption profiles. Life cycle assessment (with emphasis on end-use stage in water system) of water intense appliances was conducted to determine the environmental impacts brought by each practice. Approximately 6 to 10 % of water has been saved in the first and second year of implementation of high efficiency appliances, and with continuing savings in the third and fourth years. Water savings (gallons per household per day) for water efficiency appliances were observed at 28 (11.1%) for SH, 34.7 (13.3%) for HET, and 39.7 (14.5%) for HEW. Furthermore, the estimated contributions of high efficiency appliances for reducing water demand in the integrated water demand model were between 5 and 19% (highest in the AR program). Results indicated that adoption of more than one type of water efficiency appliance could significantly reduce residential water demand. For the sustainable water management strategies, the appropriate water conservation rate was projected to be 1 to 2 million gallons per day (MGD) through 2030. With 2 MGD of water savings, the estimated per capita water use (GPCD) could be reduced from approximately 140 to 122 GPCD. Additional efforts are needed to reduce the water demand to US EPA’s “Water Sense” conservation levels of 70 GPCD by 2030. Life cycle assessment results showed that environmental impacts (water and energy demands and greenhouse gas emissions) from end-use and demand phases are most significant within the water system, particularly due to water heating (73% for clothes washer and 93% for showerhead). Estimations of optimal lifespan for appliances (8 to 21 years) implied that earlier replacement with efficiency models is encouraged in order to minimize the environmental impacts brought by current practice.

Evaluation of wind-induced internal pressure in low-rise buildings: A multi scale experimental and numerical approach

Hurricane is one of the most destructive and costly natural hazard to the built environment and its impact on low-rise buildings, particularity, is beyond acceptable. The major objective of this research was to perform a parametric evaluation of internal pressure (IP) for wind-resistant design of low-rise buildings and wind-driven natural ventilation applications. For this purpose, a multi-scale experimental, i.e. full-scale at Wall of Wind (WoW) and small-scale at Boundary Layer Wind Tunnel (BLWT), and a Computational Fluid Dynamics (CFD) approach was adopted. This provided new capability to assess wind pressures realistically on internal volumes ranging from small spaces formed between roof tiles and its deck to attic to room partitions. Effects of sudden breaching, existing dominant openings on building envelopes as well as compartmentalization of building interior on the IP were systematically investigated. Results of this research indicated: (i) for sudden breaching of dominant openings, the transient overshooting response was lower than the subsequent steady state peak IP and internal volume correction for low-wind-speed testing facilities was necessary. For example a building without volume correction experienced a response four times faster and exhibited 30–40% lower mean and peak IP; (ii) for existing openings, vent openings uniformly distributed along the roof alleviated, whereas one sided openings aggravated the IP; (iii) larger dominant openings exhibited a higher IP on the building envelope, and an off-center opening on the wall exhibited (30–40%) higher IP than center located openings; (iv) compartmentalization amplified the intensity of IP and; (v) significant underneath pressure was measured for field tiles, warranting its consideration during net pressure evaluations. The study aimed at wind driven natural ventilation indicated: (i) the IP due to cross ventilation was 1.5 to 2.5 times higher for Ainlet/Aoutlet>1 compared to cases where Ainlet/Aoutlet<1, this in effect reduced the mixing of air inside the building and hence the ventilation effectiveness; (ii) the presence of multi-room partitioning increased the pressure differential and consequently the air exchange rate. Overall good agreement was found between the observed large-scale, small-scale and CFD based IP responses. Comparisons with ASCE 7-10 consistently demonstrated that the code underestimated peak positive and suction IP.

Development of an integrated surface and subsurface model of Everglades National Park

An integrated surface-subsurface hydrological model of Everglades National Park (ENP) was developed using MIKE SHE and MIKE 11 modeling software. The model has a resolution of 400 meters, covers approximately 1050 square miles of ENP, includes 110 miles of drainage canals with a variety of hydraulic structures, and processes hydrological information, such as evapotranspiration, precipitation, groundwater levels, canal discharges and levels, and operational schedules. Calibration was based on time series and probability of exceedance for water levels and discharges in the years 1987 through 1997. Model verification was then completed for the period of 1998 through 2005. Parameter sensitivity in uncertainty analysis showed that the model was most sensitive to the hydraulic conductivity of the regional Surficial Aquifer System, the Manning's roughness coefficient, and the leakage coefficient, which defines the canal-subsurface interaction. The model offers an enhanced predictive capability, compared to other models currently available, to simulate the flow regime in ENP and to forecast the impact of topography, water flows, and modifying operation schedules.

A methodology to estimate time varying user responses to travel time and travel time reliability in a road pricing environment

Road pricing has emerged as an effective means of managing road traffic demand while simultaneously raising additional revenues to transportation agencies. Research on the factors that govern travel decisions has shown that user preferences may be a function of the demographic characteristics of the individuals and the perceived trip attributes. However, it is not clear what are the actual trip attributes considered in the travel decision- making process, how these attributes are perceived by travelers, and how the set of trip attributes change as a function of the time of the day or from day to day. In this study, operational Intelligent Transportation Systems (ITS) archives are mined and the aggregated preferences for a priced system are extracted at a fine time aggregation level for an extended number of days. The resulting information is related to corresponding time-varying trip attributes such as travel time, travel time reliability, charged toll, and other parameters. The time-varying user preferences and trip attributes are linked together by means of a binary choice model (Logit) with a linear utility function on trip attributes. The trip attributes weights in the utility function are then dynamically estimated for each time of day by means of an adaptive, limited-memory discrete Kalman filter (ALMF). The relationship between traveler choices and travel time is assessed using different rules to capture the logic that best represents the traveler perception and the effect of the real-time information on the observed preferences. The impact of travel time reliability on traveler choices is investigated considering its multiple definitions. It can be concluded based on the results that using the ALMF algorithm allows a robust estimation of time-varying weights in the utility function at fine time aggregation levels. The high correlations among the trip attributes severely constrain the simultaneous estimation of their weights in the utility function. Despite the data limitations, it is found that, the ALMF algorithm can provide stable estimates of the choice parameters for some periods of the day. Finally, it is found that the daily variation of the user sensitivities for different periods of the day resembles a well-defined normal distribution.

Development of an enhanced hydro-geochemical model to address mercury-speciation fate and transport in aquatic environments

An awareness of mercury (Hg) contamination in various aquatic environments around the world has increased over the past decade, mostly due to its ability to concentrate in the biota. Because the presence and distribution of Hg in aquatic systems depend on many factors (e.g., pe, pH, salinity, temperature, organic and inorganic ligands, sorbents, etc.), it is crucial to understand its fate and transport in the presence of complexing constituents and natural sorbents, under those different factors. An improved understanding of the subject will support the selection of monitoring, remediation, and restoration technologies. The coupling of equilibrium chemical reactions with transport processes in the model PHREEQC offers an advantage in simulating and predicting the fate and transport of aqueous chemical species of interest. Thus, a great variety of reactive transport problems could be addressed in aquatic systems with boundary conditions of specific interest. Nevertheless, PHREEQC lacks a comprehensive thermodynamic database for Hg. Therefore, in order to use PHREEQC to address the fate and transport of Hg in aquatic environments, it is necessary to expand its thermodynamic database, confirm it and then evaluate it in applications where potential exists for its calibration and continued validation. The objectives of this study were twofold: 1) to develop, expand, and confirm the Hg database of the hydrogeochemical PHREEQC to enhance its capability to simulate the fate of Hg species in the presence of complexing constituents and natural sorbents under different conditions of pH, redox, salinity and temperature; and 2) to apply and evaluate the new database in flow and transport scenarios, at two field test beds: Oak Ridge Reservation, Oak Ridge, TN and Everglades National Park, FL, where Hg is present and is of much concern. Overall, this research enhanced the capability of the PHREEQC model to simulate the coupling of the Hg reactions in transport conditions. It also demonstrated its usefulness when applied to field situations.

Developing brownfields ranking models using decision analytic methods

The redevelopment of Brownfields has taken off in the 1990s, supported by federal and state incentives, and largely accomplished by local initiatives. Brownfields redevelopment has several associated benefits. These include the revitalization of inner-city neighborhoods, creation of jobs, stimulation of tax revenues, greater protection of public health and natural resources, the renewal and reuse existing civil infrastructure and Greenfields protection. While these benefits are numerous, the obstacles to Brownfields redevelopment are also very much alive. Redevelopment issues typically embrace a host of financial and legal liability concerns, technical and economic constraints, competing objectives, and uncertainties arising from inadequate site information. Because the resources for Brownfields redevelopment are usually limited, local programs will require creativity in addressing these existing obstacles in a manner that extends their limited resources for returning Brownfields to productive uses. Such programs may benefit from a structured and defensible decision framework to prioritize sites for redevelopment: one that incorporates the desired objectives, corresponding variables and uncertainties associated with Brownfields redevelopment. This thesis demonstrates the use of a decision analytic tool, Bayesian Influence Diagrams, and related decision analytic tools in developing quantitative decision models to evaluate and rank Brownfields sites on the basis of their redevelopment potential.

Flood risk analysis of the White Oak Creek watershed in Oak Ridge, Tennessee using MIKE SHE and MIKE 11 software

The purpose of this study was to determine the flooding potential of contaminated areas within the White Oak Creek watershed in the Oak Ridge Reservation in Tennessee. The watershed was analyzed with an integrated surface and subsurface numerical model based on MIKE SHE/MIKE 11 software. The model was calibrated and validated using five decades of historical data. A series of simulations were conducted to determine the watershed response to 25 year, 100 year and 500 year precipitation forecasts; flooding maps were generated for those events. Predicted flood events were compared to Log Pearson III flood flow frequency values for validation. This investigation also provides an improved understanding of the water fluxes between the surface and subsurface subdomains as they affect flood frequencies. In sum, this study presents crucial information to further assess the environmental risks of potential mobilization of contaminants of concern during extreme precipitation events.