Search USGS⌕ Search

Geology topics

Claire R. Tiedeman

Publications and source records attributed to Claire R. Tiedeman.

At least 19 recordsLinked to original sources
USGS2022

Microbial community response to a bioaugmentation test to degrade trichloroethylene in a fractured rock aquifer, Trenton, N.J

Bioaugmentation is a promising strategy for enhancing trichloroethylene (TCE) degradation in fractured rock. However, slow or incomplete biodegradation can lead to stalling at degradation byproducts such as 1,2-dichloroethene ( cis -DCE) and vinyl chloride (VC). Over the course of 7 years, we examined the response of groundwater microbial populations in a bioaugmentation test where an emulsified vegetable oil solution (EOS ® ) and a dechlorinating consortium (KB-1 ® ), containing the established dechlorinator Dehalococcoides (DHC), were injected into a TCE-contaminated fractured rock aquifer. Indigenous microbial communities responded within 2 days to added substrate and outcompeted KB-1 ® , and over the years of monitoring, several other notable turnover events were observed. Concentrations of ethene, the end product in reductive dechlorination, had the strongest correlations ( P < .05) with members of Candidatus Colwellbacteria but their involvement in reductive dechlorination is unknown and warrants further investigation.DHC never exceeded 0.6% relative abundance of groundwater microbial communities, despite its previously presumed importance at the site. Increased concentrations of carbon dioxide, acetic acid, and methane were positively correlated with increasing ethene concentrations; however, concentrations of cis- DCE and VC remained high by the end of the monitoring period suggesting preferential enrichment of indigenous partial dechlorinators over bioaugmented complete dechlorinators. This study highlights the importance of characterizing in situ microbial populations to understand how they can potentially enhance or inhibit augmented TCE degradation.

New Jersey↗
USGS2022

Aqueous geochemistry of waters and hydrogeology of alluvial deposits, Pinnacles National Park, California

A cooperative study between the National Park Service (NPS) and the U.S. Geological Survey (USGS) characterized groundwater quality and hydrogeology in parts of Pinnacles National Park. The water-quality investigation assessed the geochemistry of springs, wells, surface water, and precipitation and analyzed geochemistry of rock formations that affect the water chemistry through water-rock interaction. The hydrogeology investigation used geophysical and groundwater level data to characterize groundwater-flow processes in the alluvial deposits of Bear Valley and the Chalone Creek watershed. Analysis of aqueous geochemical parameters in water samples from perennial springs, water-supply wells, and surface waters was conducted for samples collected after the dry season (autumnal) and after the wet season (vernal) to assess changes in geochemistry due to changes in groundwater levels or flow resulting from precipitation. The chemistry of bulk precipitation collected during the wet season was also analyzed. Bedrock samples were analyzed for geochemical parameters to help constrain groundwater sources, flow paths, and weathering. The geochemical investigations show a correspondence between the source rock and the spring-water chemistry that can be attributed to the mineralogy of the source rock. The narrow range of strontium isotopes in water samples, sourced in geochemically and mineralogically disparate rocks, indicates that the bedrock groundwater is relatively old and has reached quasi-steady state with respect to weathering of susceptible minerals. Groundwater-level monitoring indicated that the water table is shallow—from 0 to 10 meters (m) below land surface. In southern Bear Valley and in the Chalone Creek alluvium, water levels rose and declined by several meters over each annual cycle of this study. In northern Bear Valley, water levels rose modestly over two wet seasons but declined during a third wet season. In Bear Valley, groundwater/surface-water interaction occurs along the perennial reach of Sandy Creek. Groundwater discharges to the upstream part of the reach, becomes surface water and is partly consumed by evapotranspiration, and infiltrates farther downstream. In the Chalone Creek alluvium, runoff-generated surface-water flow in intermittent stream reaches is a major component of groundwater recharge. After the onset of significant streamflow, creek water rapidly recharges groundwater until water levels rise to nearly the creek level. Groundwater levels generally remain high throughout the wet season, then gradually decline after the creek becomes dry.

California↗
USGS2022

Contributing areas to domestic wells in dipping sedimentary rocks under extreme recharge events

We use particle tracking to determine contributing areas (CAs) to wells for transient flow models that simulate cyclic domestic pumping and extreme recharge events in a small synthetic watershed underlain by dipping sedimentary rocks. The CAs consist of strike-oriented bands at locations where the water table intersects high-hydraulic conductivity beds, and from which groundwater flows to the pumping well. Factors that affect the size and location of the CAs include topographic flow directions, rock dip direction, cross-bed fracture density, and position of the well relative to streams. For an effective fracture porosity (n e ) of 10 -4 , the fastest advective travel times from CAs to wells are only a few hours. These results indicate that wells in this type of geologic setting can be highly vulnerable to contaminants or pathogens flushed into the subsurface during extreme recharge events. Increasing n e to 10 -3 results in modestly smaller CAs and delayed well vulnerability due to slower travel times. CAs determined for steady-state models of the same setting, but with long-term average recharge and pumping rates, are smaller than CAs in the models with extreme recharge. Also, the earliest-arriving particles arrive at the wells later in the steady-state models than in the extreme-recharge models. The results highlight the importance of characterizing geologic structure, simulating plausible effective porosities, and simulating pumping and recharge transience when determining CAs in fractured rock aquifers to assess well vulnerability under extreme precipitation events.

Groundwater↗
USGS2021

Nuclear magnetic resonanance logs of fractured bedrock at the Hidden Lane Landfill site, Culpeper Basin, Virginia

In May 2018, the U.S. Geological Survey (USGS) in cooperation with the U.S. Environmental Protection Agency (EPA) collected borehole nuclear magnetic resonance (bNMR) logs in three boreholes completed in sandstone and siltstone of the Balls Bluff Member of the Bull Run Formation at a Superfund Site in Culpeper Basin, Virginia. The bNMR logs were used to aid in the evaluation of the aquifer by measuring the porosity, determining the mobile and immobile fractions of water, and estimating the hydraulic conductivity, to evaluate the potential storage and transport properties at the site. The bNMR method measures the transverse (T 2 ) decay in nuclear magnetism in response to radio-frequency pulses. The relaxation decay is related to water content and the size of the pores where the water resides. In addition, the relaxation decay parameters are used to estimate hydraulic conductivity. The results were compared to other borehole logs collected at the site and to regional groundwater investigations in similar rock formations.

Virginia↗
USGS2019

Hydraulic tomography: 3D hydraulic conductivity, fracture network, and connectivity in mudstone

We present the first demonstration of hydraulic tomography (HT) to estimate the three-dimensional (3D) hydraulic conductivity ( K ) distribution of a fractured aquifer at high-resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer-isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution ( K blocks <1 m 3 ), using no strict assumptions about K variability or fracture statistics. The resulting 3D K estimate ranges from approximately 0.1 (highest- K fractures) to approximately 10 −13 m/s (unfractured mudstone). Important estimated features include: (1) a highly fractured zone (HFZ) consisting of a sequence of high- K bedding-plane fractures; (2) a low- K zone that disrupts the HFZ; (3) several secondary fractures of limited extent; and (4) regions of very low- K rock matrix. The 3D K estimate explains complex drawdown behavior observed in the field. Drawdown tracing and particle tracking simulations reveal a 3D fracture network within the estimated K distribution, and connectivity routes through the network. Model fit is best in the shallower part of the wellfield, with high density of observations and tests. The capabilities of HT demonstrated for 3D fractured aquifer characterization at HRFS may support improved in situ remediation for contaminant source zones, and applications in mining, repository assessment, or geotechnical engineering.

New Jersey↗
USGS2019

The complex spatial distribution of trichloroethene and the probability of NAPL occurrence in the rock matrix of a mudstone aquifer

Methanol extractions for chloroethene analyses are conducted on rock samples from seven closely spaced coreholes in a mudstone aquifer that was subject to releases of the nonaqueous phase liquid (NAPL) form of trichloroethene (TCE) between the 1950's and 1990's. Although TCE concentration in the rock matrix over the length of coreholes is dictated by proximity to subhorizontal bedding planefractures, elevated TCE concentrations in the rock matrix are not continuous along the most permeable bedding plane fractures. A complex configuration of subvertical and subhorizontal fractures appears to be responsible for the TCE distribution from prior TCE releases at land surface. Phase partitioning calculations of TCE in the rock matrix show that most TCE is adsorbed to solid surfaces because of the large fraction of organic carbon ( f oc ) in the mudstone. Large TCE content in some cores indicate the likely presence of the NAPL form of TCE in the rock matrix. Using average values of porosity ( n ) and f oc in phase partitioning calculations identifies a number of locations of possible NAPL occurrence in the rock matrix. Samples of mudstone analyzed for n and f oc show variability in these properties over several orders of magnitude. Accounting for this variability in phase partitioning calculations identifies a probability of NAPL occurrence, P NAPL . The spatial variability of P NAPL along coreholes identifies a configuration that may be attributed to a TCE source zone that has evolved after emplacement due to NAPL dissolution, adsorption, and matrix diffusion.

New Jersey↗
USGS2017

Bioremediation in fractured rock: 1. Modeling to inform design, monitoring, and expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.

Groundwater↗
USGS2017

Bioremediation in fractured rock: 2. Mobilization of chloroethene compounds from the rock matrix

A mass balance is formulated to evaluate the mobilization of chlorinated ethene compounds (CE) from the rock matrix of a fractured mudstone aquifer under pre- and postbioremediation conditions. The analysis relies on a sparse number of monitoring locations and is constrained by a detailed description of the groundwater flow regime. Groundwater flow modeling developed under the site characterization identified groundwater fluxes to formulate the CE mass balance in the rock volume exposed to the injected remediation amendments. Differences in the CE fluxes into and out of the rock volume identify the total CE mobilized from diffusion, desorption, and nonaqueous phase liquid dissolution under pre- and postinjection conditions. The initial CE mass in the rock matrix prior to remediation is estimated using analyses of CE in rock core. The CE mass mobilized per year under preinjection conditions is small relative to the total CE mass in the rock, indicating that current pump-and-treat and natural attenuation conditions are likely to require hundreds of years to achieve groundwater concentrations that meet regulatory guidelines. The postinjection CE mobilization rate increased by approximately an order of magnitude over the 5 years of monitoring after the amendment injection. This rate is likely to decrease and additional remediation applications over several decades would still be needed to reduce CE mass in the rock matrix to levels where groundwater concentrations in fractures achieve regulatory standards.

New Jersey↗
USGS2015

Imaging pathways in fractured rock using three-dimensional electrical resistivity tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three-dimensional cross borehole electrical resistivity tomography (ERT) in a 9&thinsp;m (diameter)&thinsp;&times;&thinsp;15&thinsp;m well field to capture high-resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time-lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.

New Jersey↗
USGS2014

Integration of stable carbon isotope, microbial community, dissolved hydrogen gas, and 2 H H 2 O tracer data to assess bioaugmentation for chlorinated ethene degradation in fractured rocks

An in situ bioaugmentation (BA) experiment was conducted to understand processes controlling microbial dechlorination of trichloroethene (TCE) in groundwater at the Naval Air Warfare Center (NAWC), West Trenton, NJ. In the BA experiment, an electron donor (emulsified vegetable oil and sodium lactate) and a chloro-respiring microbial consortium were injected into a well in fractured mudstone of Triassic age. Water enriched in 2 H was also injected as a tracer of the BA solution, to monitor advective transport processes. The changes in concentration and the δ 13 C of TCE, cis-dichloroethene (cis-DCE), and vinyl chloride (VC); the δ 2 H of water; changes in the abundance of the microbial communities; and the concentration of dissolved H2 gas compared to pre- test conditions, provided multiple lines of evidence that enhanced biodegradation occurred in the injection well and in two downgradient wells. For those wells where the biodegradation was stimulated intensively, the sum of the molar chlorinated ethene (CE) concentrations in post-BA water was higher than that of the sum of the pre-BA background molar CE concentrations. The concentration ratios of TCE/(cis-DCE + VC) indicated that the increase in molar CE concentration may result from additional TCE mobilized from the rock matrix in response to the oil injection or due to desorption/diffusion. The stable carbon isotope mass-balance calculations show that the weighted average 13 C isotope of the CEs was enriched for around a year compared to the background value in a two year monitoring period, an effective indication that dechlorination of VC was occurring. Insights gained from this study can be applied to efforts to use BA in other fractured rock systems. The study demonstrates that a BA approach can substantially enhance in situ bioremediation not only in fractures connected to the injection well, but also in the rock matrix around the well due to processes such as diffusion and desorption. Because the effect of the BA was intensive only in wells where an amendment was distributed during injection, it is necessary to adequately distribute the amendments throughout the fractured rock to achieve substantial bioremediation. The slowdown in BA effect after a year is due to some extend to the decrease abundant of appropriate microbes, but more likely the decreased concentration of electron donor.

New Jersey↗
USGS2013

Effect of correlated observation error on parameters, predictions, and uncertainty

Correlations among observation errors are typically omitted when calculating observation weights for model calibration by inverse methods. We explore the effects of omitting these correlations on estimates of parameters, predictions, and uncertainties. First, we develop a new analytical expression for the difference in parameter variance estimated with and without error correlations for a simple one-parameter two-observation inverse model. Results indicate that omitting error correlations from both the weight matrix and the variance calculation can either increase or decrease the parameter variance, depending on the values of error correlation (ρ) and the ratio of dimensionless scaled sensitivities (r dss ). For small ρ, the difference in variance is always small, but for large ρ, the difference varies widely depending on the sign and magnitude of r dss . Next, we consider a groundwater reactive transport model of denitrification with four parameters and correlated geochemical observation errors that are computed by an error-propagation approach that is new for hydrogeologic studies. We compare parameter estimates, predictions, and uncertainties obtained with and without the error correlations. Omitting the correlations modestly to substantially changes parameter estimates, and causes both increases and decreases of parameter variances, consistent with the analytical expression. Differences in predictions for the models calibrated with and without error correlations can be greater than parameter differences when both are considered relative to their respective confidence intervals. These results indicate that including observation error correlations in weighting for nonlinear regression can have important effects on parameter estimates, predictions, and their respective uncertainties.

Water Resources Research↗
USGS2013

Knowledge, transparency, and refutability in groundwater models, an example from the Death Valley regional groundwater flow system

This work demonstrates how available knowledge can be used to build more transparent and refutable computer models of groundwater systems. The Death Valley regional groundwater flow system, which surrounds a proposed site for a high level nuclear waste repository of the United States of America, and the Nevada National Security Site (NNSS), where nuclear weapons were tested, is used to explore model adequacy, identify parameters important to (and informed by) observations, and identify existing old and potential new observations important to predictions. Model development is pursued using a set of fundamental questions addressed with carefully designed metrics. Critical methods include using a hydrogeologic model, managing model nonlinearity by designing models that are robust while maintaining realism, using error-based weighting to combine disparate types of data, and identifying important and unimportant parameters and observations and optimizing parameter values with computationally frugal schemes. The frugal schemes employed in this study require relatively few (10–1000 s), parallelizable model runs. This is beneficial because models able to approximate the complex site geology defensibly tend to have high computational cost. The issue of model defensibility is particularly important given the contentious political issues involved.

California, Nevada↗
USGS2012

Enhanced dichloroethene biodegradation in fractured rock under biostimulated and bioaugmented conditions

Significant microbial reductive dechlorination of [1,2 14 C] cis-dichloroethene (DCE) was observed in anoxic microcosms prepared with unamended, fractured rock aquifer materials, which were colonized in situ at multiple depths in two boreholes at the Naval Air Warfare Center (NAWC) in West Trenton, New Jersey. The lack of significant reductive dechlorination in corresponding water-only treatments indicated that chlororespiration activity in unamended, fractured rock treatments was primarily associated with colonized core material. In these unamended fractured rock microcosms, activity was highest in the shallow zones and generally decreased with increasing depth. Electron-donor amendment (biostimulation) enhanced chlororespiration in some but not all treatments. In contrast, combining electron-donor amendment with KB1 amendment (bioaugmentation) enhanced chlororespiration in all treatments and substantially reduced the variability in chlororespiration activity both within and between treatments. These results indicate (1) that a potential for chlororespiration-based bioremediation exists at NAWC Trenton but is limited under nonengineered conditions, (2) that the limitation on chlororespiration activity is not entirely due to electron-donor availability, and (3) that a bioaugmentation approach can substantially enhance in situ bioremediation if the requisite amendments can be adequately distributed throughout the fractured rock matrix.

New Jersey↗
USGS2012

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The east-central Florida transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration, runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into evapotranspiration, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average evapotranspiration (ET) over the 1995 to 2006 period was 34.5 inches per year, compared to the calculated average ET rate of 36.6 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.6 inches per year, compared with the calculated average of 3.2 inches per year from the model-independent waterbudget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 800 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida↗
USGS2012

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The East-Central Florida Transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration (ET), runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into ET, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average ET over the 1995 to 2006 period was 34.47 inches per year, compared to the calculated average ET rate of 36.39 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.58 inches per year, compared with the calculated average of 3.39 inches per year from the model-independent water-budget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 920 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida↗