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

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

Response of the salt-freshwater interface in a coastal aquifer to a wave-induced groundwater pulse: field observations and modelling

Field observations on an unconfined coastal aquifer showed that a groundwater pulse, generated by it moderate (significant wave height, H-sig similar to 4.5 m) wave/storm event, induced significant oscillations in the salt-freshwater interface of the order of several metres in the horizontal direction. A dynamic sharp-interface model is developed to quantify the mechanism of these interface oscillations. The model uses the 50% seawater salinity contour as the location of the equivalent sharp-interface. The model was calibrated against the observed groundwater table fluctuations. It predicted reasonably well the interface oscillations with a slight over-prediction of the oscillation magnitude and a steepening of the interface. The neglect of mixing in the salt-freshwater mixing zone by the sharp-interface model is suggested as a possible contributor to the discrepancies between the model predictions and observations. In contrast with the significant wave effects, there was no observable response of the interface to diurnal or semidiurnal tides. (C) 2004 Elsevier Ltd. All rights reserved.

Response of coastal groundwater table to offshore storms

Large groundwater table fluctuations were observed in a coastal aquifer during an offshore storm. The storm induced significant changes of the mean shoreline elevation, characterized by a pulse-like oscillation. This pulse propagated in the aquifer, resulting in the water table fluctuations. A general analytical solution is derived to quantify this new mechanism of water table fluctuation. The solution is applied to field observations and is found to be able to predict reasonably well the observed storm-induced water table fluctuations. Based on the analytical solution, the damping characteristics and phase shift of the oscillation as it propagates inland are examined.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Quantitative risk assessment of groundwater quality utilizing GIS technology and coupled groundwater models

The thesis presents a two-dimensional Risk Assessment Method (RAM) where the assessment of risk to the groundwater resources incorporates both the quantification of the probability of the occurrence of contaminant source terms, as well as the assessment of the resultant impacts. The approach emphasizes the need for a greater dependency on the potential pollution sources, rather than the traditional approach where assessment is based mainly on the intrinsic geo-hydrologic parameters. The risk is calculated using Monte Carlo simulation methods whereby random pollution events were generated to the same distribution as historically occurring events or a priori potential probability distribution. Integrated mathematical models then simulate contaminant concentrations at the predefined monitoring points within the aquifer. The spatial and temporal distributions of the concentrations were calculated from repeated realisations, and the number of times when a user defined concentration magnitude was exceeded is quantified as a risk. The method was setup by integrating MODFLOW-2000, MT3DMS and a FORTRAN coded risk model, and automated, using a DOS batch processing file. GIS software was employed in producing the input files and for the presentation of the results. The functionalities of the method, as well as its sensitivities to the model grid sizes, contaminant loading rates, length of stress periods, and the historical frequencies of occurrence of pollution events were evaluated using hypothetical scenarios and a case study. Chloride-related pollution sources were compiled and used as indicative potential contaminant sources for the case study. At any active model cell, if a random generated number is less than the probability of pollution occurrence, then the risk model will generate synthetic contaminant source term as an input into the transport model. The results of the applications of the method are presented in the form of tables, graphs and spatial maps. Varying the model grid sizes indicates no significant effects on the simulated groundwater head. The simulated frequency of daily occurrence of pollution incidents is also independent of the model dimensions. However, the simulated total contaminant mass generated within the aquifer, and the associated volumetric numerical error appear to increase with the increasing grid sizes. Also, the migration of contaminant plume advances faster with the coarse grid sizes as compared to the finer grid sizes. The number of daily contaminant source terms generated and consequently the total mass of contaminant within the aquifer increases in a non linear proportion to the increasing frequency of occurrence of pollution events. The risk of pollution from a number of sources all occurring by chance together was evaluated, and quantitatively presented as risk maps. This capability to combine the risk to a groundwater feature from numerous potential sources of pollution proved to be a great asset to the method, and a large benefit over the contemporary risk and vulnerability methods.

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.

Quantifying the effects of temperature and concentration on variable-density flow in numerical modeling of groundwater systems: Implications for predictive uncertainty and data collection

Groundwater systems of different densities are often mathematically modeled to understand and predict environmental behavior such as seawater intrusion or submarine groundwater discharge. Additional data collection may be justified if it will cost-effectively aid in reducing the uncertainty of a model's prediction. The collection of salinity, as well as, temperature data could aid in reducing predictive uncertainty in a variable-density model. However, before numerical models can be created, rigorous testing of the modeling code needs to be completed. This research documents the benchmark testing of a new modeling code, SEAWAT Version 4. The benchmark problems include various combinations of density-dependent flow resulting from variations in concentration and temperature. The verified code, SEAWAT, was then applied to two different hydrological analyses to explore the capacity of a variable-density model to guide data collection. ^ The first analysis tested a linear method to guide data collection by quantifying the contribution of different data types and locations toward reducing predictive uncertainty in a nonlinear variable-density flow and transport model. The relative contributions of temperature and concentration measurements, at different locations within a simulated carbonate platform, for predicting movement of the saltwater interface were assessed. Results from the method showed that concentration data had greater worth than temperature data in reducing predictive uncertainty in this case. Results also indicated that a linear method could be used to quantify data worth in a nonlinear model. ^ The second hydrological analysis utilized a model to identify the transient response of the salinity, temperature, age, and amount of submarine groundwater discharge to changes in tidal ocean stage, seasonal temperature variations, and different types of geology. The model was compared to multiple kinds of data to (1) calibrate and verify the model, and (2) explore the potential for the model to be used to guide the collection of data using techniques such as electromagnetic resistivity, thermal imagery, and seepage meters. Results indicated that the model can be used to give insight to submarine groundwater discharge and be used to guide data collection. ^

Comparative study of dissolved organic matter from groundwater and surface water in the Florida coastal Everglades using multi-dimensional spectrofluorometry combined with multivariate statistics

Dissolved organic matter (DOM) in groundwater and surface water samples from the Florida coastal Everglades were studied using excitation–emission matrix fluorescence modeled through parallel factor analysis (EEM-PARAFAC). DOM in both surface and groundwater from the eastern Everglades S332 basin reflected a terrestrial-derived fingerprint through dominantly higher abundances of humic-like PARAFAC components. In contrast, surface water DOM from northeastern Florida Bay featured a microbial-derived DOM signature based on the higher abundance of microbial humic-like and protein-like components consistent with its marine source. Surprisingly, groundwater DOM from northeastern Florida Bay reflected a terrestrial-derived source except for samples from central Florida Bay well, which mirrored a combination of terrestrial and marine end-member origin. Furthermore, surface water and groundwater displayed effects of different degradation pathways such as photodegradation and biodegradation as exemplified by two PARAFAC components seemingly indicative of such degradation processes. Finally, Principal Component Analysis of the EEM-PARAFAC data was able to distinguish and classify most of the samples according to DOM origins and degradation processes experienced, except for a small overlap of S332 surface water and groundwater, implying rather active surface-to-ground water interaction in some sites particularly during the rainy season. This study highlights that EEM-PARAFAC could be used successfully to trace and differentiate DOM from diverse sources across both horizontal and vertical flow profiles, and as such could be a convenient and useful tool for the better understanding of hydrological interactions and carbon biogeochemical cycling.

Hydrochemical and microbiological properties of groundwater from pristine aquifer sampled at the Lower Rhine Embayment

The hydrochemistry and the microbial diversity of a pristine aquifer system near Garzweiler, Germany next to the open-pit lignite mine Garzweiler 1, were characterized. Hydrogeochemical and isotopic data indicate a recent activity of sulfate-reducing bacteria in the Tertiary marine sands. The community structure in the aquifer was studied by fluorescence in situ hybridization (FISH). Up to 7.3 x 10**5 cells/ml were detected by DAPIstaining. Bacteria (identified by the probe EUB338) were dominant, representing 51.9% of the total cell number (DAPI). Another 25.7% of total cell were affiliated with the domain Archaea as identified by the probe ARCH915. Within the domain Bacteria, the beta-Proteobacteria were most abundant (21.0% of total cell counts). Using genusspecific probes for sulfate-reducing bacteria (SRB), 2.5% of the total cells were identified as members of the genus Desulfotomaculum. This reflects the predominant role these microorganisms have been found to play in sulfatereducing zones of aquifers at other sites. Previously, all SRB cultured from this site were from the spore-forming genera Desulfotomaculum and Desulfosporosinus. Samples were taken after pumping for >= 40 min and after parameters such as temperature, pH, redox potential, oxygen and conductivity of the groundwater had remained stable for >= 15 min due to recharge of aquifer water. Hybridization and microscopy counts of hybridized and 4',6'-diamidino-2-phenylindole (DAPI)- stained cells were performed as described in Snaidr et al., (1997, http://aem.asm.org/content/63/7/2884.full.pdf). Means were calculated from 10 to 20 randomly chosen fields on each filter section, corresponding to 800-1000 DAPI stained cells. Counting results were always corrected by subtracting signals observed with the probe NON338. Formamide concentrations and oligonucleotide probes used please see further details.

Hydrochemical properties of groundwater samples in Semarang area, Java Island, Indonesia (1992, 1993, 2003, 2006, and 2007)

In Semarang City, groundwater has been exploited as a natural resource since 1841. The groundwater exploited in deep wells is concentrated in confined aquifers. The previous hydrogeological model was developed in one unit of aquifer and refined then by using several hydrostratigraphical units following a regional hydrogeological map without any further analysis. At present, there is a lack of precise hydrogeological model which integrates geological and hydrogeological data, in particular for multiple aquifers in Semarang. Thus, the aim of this paper is to develop a hydrogeological model for the multiple aquifers in Semarang using an integrated data approach. Groundwater samples in the confined aquifers have been analyzed to define the water type and its lateral distribution. Two hydrogeological cross sections were then created based on several borelog data to define a hydrostratigraphical unit (HSU). The HSU result indicates the hydrogeological model of Semarang consists of two aquifers, three aquitards, and one aquiclude. Aquifer 1 is unconfined, while Aquifer 2 is confined. Aquifer 2 is classified into three groups (2a, 2b, and 2c) based on analyses of major ion content and hydrostratigraphical cross sections.

(Table 1) Water chemistry of groundwater and snow samples from three sites in Ontario, Canada

Snow samples collected from hand-dug pits at two sites in Simcoe County, Ontario, Canada were analysed for major and trace elements using the clean lab methods established for polar ice. Potentially toxic, chalcophile elements are highly enriched in snow, relative to their natural abundance in crustal rocks, with enrichment factor (EF) values (calculated using Sc) in the range 107 to 1081 for Ag, As, Bi, Cd, Cu, Mo, Pb, Sb, Te, and Zn. Relative to M/Sc ratios in snow, water samples collected at two artesian flows in this area are significantly depleted in Ag, Al, Be, Bi, Cd, Cr, Cu, Ni, Pb, Sb, Tl, V, and Zn at both sites, and in Co, Th and Tl at one of the sites. The removal from the waters of these elements is presumably due to such processes as physical retention (filtration) of metal-bearing atmospheric aerosols by organic and mineral soil components as well as adsorption and surface complexation of ionic species onto organic, metal oxyhydroxide and clay mineral surfaces. In the case of Pb, the removal processes are so effective that apparently ''natural'' ratios of Pb to Sc are found in the groundwaters. Tritium measurements show that the groundwater at one of the sites is modern (ie not more than 30 years old) meaning that the inputs of Pb and other trace elements to the groundwaters may originally have been much higher than they are today; the M/Sc ratios measured in the groundwaters today, therefore, represent a conservative estimate of the extent of metal removal along the flow path. Lithogenic elements significantly enriched in the groundwaters at both sites include Ba, Ca, Li, Mg, Mn, Na, Rb, S, Si, Sr, and Ti. The abundance of these elements can largely be explained in terms of weathering of the dominant silicate (plagioclase, potassium feldspar, amphibole and biotite) and carbonate minerals (calcite, dolomite and ankerite) in the soils and sediments of the watershed. Arsenic, Mo, Te, and especially U are also highly enriched in the groundwaters, due to chemical weathering: these could easily be explained if there are small amounts of sulfides (As, Mo, Te) and apatite (U) in the soils of the source area. Elements neither significantly enriched nor depleted at both sites include Fe, Ga, Ge, and P.

Modeling long-term variability and change of soil moisture and groundwater level - from catchment to global scale

The water stored in and flowing through the subsurface is fundamental for sustaining human activities and needs, feeding water and its constituents to surface water bodies and supporting the functioning of their ecosystems. Quantifying the changes that affect the subsurface water is crucial for our understanding of its dynamics and changes driven by climate change and other changes in the landscape, such as in land-use and water-use. It is inherently difficult to directly measure soil moisture and groundwater levels over large spatial scales and long times. Models are therefore needed to capture the soil moisture and groundwater level dynamics over such large spatiotemporal scales. This thesis develops a modeling framework that allows for long-term catchment-scale screening of soil moisture and groundwater level changes. The novelty in this development resides in an explicit link drawn between catchment-scale hydroclimatic and soil hydraulics conditions, using observed runoff data as an approximation of soil water flux and accounting for the effects of snow storage-melting dynamics on that flux. Both past and future relative changes can be assessed by use of this modeling framework, with future change projections based on common climate model outputs. By direct model-observation comparison, the thesis shows that the developed modeling framework can reproduce the temporal variability of large-scale changes in soil water storage, as obtained from the GRACE satellite product, for most of 25 large study catchments around the world. Also compared with locally measured soil water content and groundwater level in 10 U.S. catchments, the modeling approach can reasonably well reproduce relative seasonal fluctuations around long-term average values. The developed modeling framework is further used to project soil moisture changes due to expected future climate change for 81 catchments around the world. The future soil moisture changes depend on the considered radiative forcing scenario (RCP) but are overall large for the occurrence frequency of dry and wet events and the inter-annual variability of seasonal soil moisture. These changes tend to be higher for the dry events and the dry season, respectively, than for the corresponding wet quantities, indicating increased drought risk for some parts of the world.

Development of an Observation Well Network in the Mahomet Aquifer of East-Central Illinois

Illinois State Water Survey

Vegetables production in NW Portugal, agriculture practices and groundwater quality

In NW of Portugal, the region between Esposende and Vila do Conde was established into a unique system called masseiras; masseiras are small fields that were created at the beginning of XX century by removing sand from dunes. Since the 80’s sand between closed masseiras has been removed giving way for larger fields, the “enlarged masseira with greenhouses”. This region is part of vulnerable zone nº1 (VZ1). Nowadays, the masseiras have almost disappeared; they had been replaced by “enlarged masseiras with plastic greenhouses” exploited by a new generation of growers. As water quality is an indicator of agriculture practices, 8 wells were selected and water samples were collected on four dates from September 2009 to May 2010. Laboratorial analysis was performed on nitrates and 28 pesticides. The results confirm that water in VZ1 is contaminated with nitrates (as it was expected) and also with pesticides. The water was contaminated in both dug wells and drilled wells all year round. Masseiras’ wells have more nitrates than greenhouses’ wells, probably due to poor agricultural practices used by growers who still employ high quantities of fertilizers and water. Pesticides detected in water samples include some that have been sprayed in the field and others that have not. Some forbidden substances were found. Our results confirm that, in VZ1, there is a diffusion contamination, which origin is unknown. In this communication, we put in picture the history of masseiras fields: how masseiras fields have been created and how they have disappeared; what trends exist in pest management strategies and the use of fertilizers, pesticides and water resources. These results reinforce the need to change the agricultural system in this region of traditional vegetable production.