Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interpretation”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70Linked to original sources
USGS2018

Hydrogeologic controls and geochemical indicators of groundwater movement in the Niles Cone and southern East Bay Plain groundwater subbasins, Alameda County, California

Beginning in the 1970s, Alameda County Water District began infiltrating imported water through ponds in repurposed gravel quarries at the Quarry Lakes Regional Park, in the Niles Cone groundwater subbasin, to recharge groundwater and to minimize intrusion of saline, San Francisco Bay water into freshwater aquifers. Hydraulic connection between distinct aquifers underlying Quarry Lakes allows water to recharge the upper aquifer system to depths of 400 feet below land surface, and the Deep aquifer to depths of more than 650 feet. Previous studies of the Niles Cone and southern East Bay Plain groundwater subbasins suggested that these two subbasins may be hydraulically connected. Characterization of storage capacities and hydraulic properties of the complex aquifers and the structural and stratigraphic controls on groundwater movement aids in optimal storage and recovery of recharged water and provides information on the ability of aquifers shared by different water management agencies to fulfill competing storage and extraction demands. The movement of recharge water through the Niles Cone groundwater subbasin from Quarry Lakes and the possible hydraulic connection between the Niles Cone and the southern East Bay Plain groundwater subbasins were investigated using interferometric synthetic aperture radar (InSAR), water-chemistry, and isotopic data, including tritium/helium-3, helium-4, and carbon-14 age-dating techniques. InSAR data collected during refilling of the Quarry Lakes recharge ponds show corresponding ground-surface displacement. Maximum uplift was about 0.8 inches, reasonable for elastic expansion of sedimentary materials experiencing an increase in hydraulic head that resulted from pond refilling. Sodium concentrations increase while calcium and magnesium concentrations in groundwater decrease along groundwater flowpaths from the Niles Cone groundwater subbasin through the Deep aquifer to the northwest toward the southern East Bay Plain groundwater subbasin. Residual effects of pre-1970s intrusion of saline water from San Francisco Bay, including high chloride concentrations in groundwater, are evident in parts of the Niles Cone subbasin. Noble gas recharge temperatures indicate two primary recharge sources (Quarry Lakes and Alameda Creek) in the Niles Cone groundwater subbasin. Although recharge at Quarry Lakes affects hydraulic heads as far as the transition zone between the Niles Cone and East Bay Plain groundwater subbasins (about 5 miles), the effect of recharged water on water quality is only apparent in wells near (less than 2 miles) recharge sources. Groundwater chemistry from upper aquifer system wells near Quarry Lakes showed an evaporated signal (less negative oxygen and hydrogen isotopic values) relative to surrounding groundwater and a tritium concentration (2 tritium units) consistent with recently recharged water from a surface-water impoundment. Uncorrected carbon-14 activities measured in water sampled from wells in the Niles Cone groundwater subbasin range from 16 to 100 percent modern carbon (pmC). The geochemical reaction modeling software NETPATH was used to interpret carbon-14 ages along a flowpath from Quarry Lakes toward the East Bay Plain groundwater subbasin. Model results indicate that changes in groundwater chemistry are controlled by cation exchange on clay minerals and weathering of primary silicate minerals. Old groundwater (lower carbon-14 activities) is characterized by high dissolved silica and pH. Interpreted carbon-14 ages ranged from 830 to more than 7,000 years before present and are less than helium-4 ages that range from 2,000 to greater than 11,000 years before present. The average horizontal groundwater velocity along the studied flowpath, as calculated using interpreted carbon-14 ages, through the Deep aquifer of the Niles Cone groundwater subbasin is between 3 and 12 feet per year. The groundwater velocity decreases near the boundary of the transition zone to the southern East Bay Plain groundwater subbasin to about 0.5 feet per year. These changes may result from water recharged from different sources converging in flowpaths north of the transition zone, or a boundary to flow between the Niles Cone and southern East Bay Plain groundwater subbasins, likely owing to changes in lithology caused by depositional patterns.

California↗
USGS2018

Data analysis considerations for pesticides determined by National Water Quality Laboratory schedule 2437

In 2013, the U.S. Geological Survey National Water Quality Laboratory (NWQL) made a new method available for the analysis of pesticides in filtered water samples: laboratory schedule 2437. Schedule 2437 is an improvement on previous analytical methods because it determines the concentrations of 225 fungicides, herbicides, insecticides, and associated degradates in one method at similar or lower concentrations than previously available methods. Additionally, the pesticides included in schedule 2437 were strategically identified in a prioritization analysis that assessed likelihood of occurrence, prevalence of use, and potential toxicity. When the NWQL reports pesticide concentrations for analytes in schedule 2437, the laboratory also provides supplemental information useful to data users for assessing method performance and understanding data quality. That supplemental information is discussed in this report, along with an initial analysis of analytical recovery of pesticides in water-quality samples analyzed by schedule 2437 during 2013–2015. A total of 523 field matrix spike samples and their paired environmental samples and 277 laboratory reagent spike samples were analyzed for this report (1,323 samples total). These samples were collected in the field as part of the U.S. Geological Survey National Water-Quality Assessment groundwater and surface-water studies and as part of the NWQL quality-control program. This report reviews how pesticide samples are processed by the NWQL, addresses how to obtain all the data necessary to interpret pesticide concentrations, explains the circumstances that result in a reporting level change or the occurrence of a raised reporting level, and describes the calculation and assessment of recovery. This report also discusses reasons why a data user might choose to exclude data in an interpretive analysis and outlines the approach used to identify the potential for decreased data quality in the assessment of method recovery. The information provided in this report is essential to understanding pesticide data determined by schedule 2437 and should be reviewed before interpretation of these data.

Scientific Investigations Report↗
USGS2018

Documentation of single-well aquifer tests and integrated borehole analyses, Pahute Mesa and Vicinity, Nevada

Single-well aquifer testing has been carried out at Pahute Mesa in southern Nevada since 1962. These tests include single-well pumping and slug tests to estimate geologic formation hydraulic properties. Initially, aquifer tests focused on identifying low-permeability rocks suitable for testing large-yield nuclear devices, whereas later hydrologic investigations focused on potential subsurface transport of radionuclides away from Pahute Mesa. The number of single-well aquifer tests analyzed for Pahute Mesa and vicinity was much greater than the number of actual tests because each response to a water-level displacement was interpreted multiple times using different analytical models by different investigators, which were reported as distinct, individual aquifer tests. Some aquifer tests also were reinterpreted multiple times because the wells had not been defined uniquely. Previous databases and reports identified wells by different borehole names, reported various depths to top and bottom of open intervals, and used inconsistent wetted-aquifer thicknesses, which were coupled with inappropriate analytical methods to interpret aquifer-test data. This work standardized the reporting of single-well aquifer-test results at Pahute Mesa and vicinity. A hydraulic-parameter database was created that uniquely identifies tested wells, reports all replicated aquifer tests in each well, and identifies the best transmissivity estimate for each well. Datasets, interpreted models, and results from 1,459 analyses are presented and evaluated for 360 unique wells. Integrated borehole analyses were used to estimate hydraulic conductivity by depth in Pahute Mesa because vertical variations in hydraulic conductivity greatly affect groundwater velocities, which directly affect subsurface transport rates and directions. The integrated analyses reconciled differences between transmissivity estimated from pumping tests and the summed transmissivity from straddle-packer slug tests done at different depth intervals in a borehole. An integration borehole database was created that presents integrated aquifer-test results from 17 boreholes.

Nevada↗
USGS2018

Chemical and isotopic characteristics of methane in groundwater of Ohio, 2016

In 2016, the U.S. Geological Survey, in cooperation with the Ohio Water Development Authority, investigated the hydrogeologic setting, chemical and isotopic characteristics, and origin of methane in groundwater of Ohio. Understanding the occurrence and distribution of methane in groundwater is important in terms of public safety because methane in water wells can pose a risk of explosion. In addition, documenting the chemical and isotopic characteristics of methane in groundwater can make an important contribution to future stray gas investigations. Water samples were collected from 15 domestic water wells known to produce methane, which were in 12 counties in diverse parts of Ohio. The wells were 75–345 feet deep and tapped a range of aquifer types, including glacial deposits and bedrock of Upper Ordovician, Upper Devonian, Lower Mississippian, and Pennsylvanian ages. Although the hydrogeologic settings were varied, there was a broad similarity among the well sites in that the bedrock was predominantly shale and the glacial deposits were predominantly clay. The wells were sampled for dissolved inorganic constituents; dissolved organic carbon; methane and other dissolved gases; stable isotopes (carbon, hydrogen, and oxygen) of methane, water, and dissolved inorganic carbon; and carbon-14 of methane. Gas composition and stable isotopes of methane were used to differentiate thermogenic and microbial methane. The degree of fractionation of hydrogen and carbon isotopes was used to evaluate the pathway of microbial methanogenesis (carbon dioxide [CO 2 ] reduction or acetate fermentation) and the effects of secondary processes such as oxidation, mixing, and migration. The concentration of carbon-14 of methane was used to evaluate the relative age of the carbon source. The quality of water from the 15 wells differed greatly; water types ranged from CaMgHCO 3 to NaCl, and total dissolved solids concentrations ranged from 318 to 2,940 milligrams per liter (mg/L). Methane concentrations ranged from 1.2 to 120 mg/L. Of the 15 samples, 12 had methane concentrations greater than 28 mg/L, the level that can pose a risk of explosion. Of the 15 samples, 12 had chemical and isotopic characteristics or "signatures" consistent with microbial methane formed by CO 2 reduction. CO 2 reduction is commonly associated with microbial degradation of organic matter in anaerobic aquifers and with the formation of microbial shale gas and coalbed methane along margins of sedimentary basins. Two of 15 samples were interpreted as having a component of thermogenic methane based on the δ 13 C of methane (−50.96 and −47.74 parts per thousand [per mil]) and gas dryness (28 and 5). One of 15 samples (from the shallowest well) had chemical and isotopic characteristics consistent with methane oxidation by sulfate reduction based on light δ 13 C of dissolved inorganic carbon (−31.6 per mil) and evidence of sulfate reduction in terms of the odor and appearance of the water. For the 12 samples interpreted as microbial methane formed by CO 2 reduction, the δ 13 C of methane varied from −75 to −56 per mil. Multiple samples from the same aquifer demonstrated a general trend of increasing δ 13 C of methane with depth. Samples with lighter δ 13 C of methane (−75 to −62 per mil) were from shallower wells (or wells with shallow open intervals), and the isotopic signature of the water was consistent with modern or postglacial groundwater recharge. Three samples with heavier δ 13 C of methane (−61 to −56 per mil) were from deeper wells or more confined aquifers where the isotopic signature of water was consistent with older (glacial) recharge. In addition, δ 13 C of dissolved inorganic carbon was enriched (+12 to +18.9 per mil), and carbon-14 of methane was consistent with carbon associated with Paleozoic bedrock or older glacial deposits. These observations are generally consistent with increased Rayleigh-type fractionation at greater depths; however, other interpretations are possible. Isotopic signatures can be ambiguous, especially in areas with complex geologic histories that include multiple episodes of migration, mixing, and (or) oxidation. Many of the wells were in proximity to multiple potential natural and anthropogenic pathways of methane migration; however, it is not possible to determine if the methane in any of the wells is related to human activities based on the chemical and isotopic data collected for this study.

Ohio↗
USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗
USGS2020

Quality of data from the U.S. Geological Survey National Water Quality Network for water years 2013–17

Water samples from 122 sites in the U.S. Geological Survey National Water Quality Network were collected in 2013–17 to document ambient water-quality conditions in surface water of the United States and to determine status and trends of loads and concentrations for nutrients, contaminants, and sediment to estuaries and streams. Quality-control (QC) samples collected in the field with environmental samples were combined with QC samples from laboratory processing to provide information and documentation about the quality of the environmental data. Quality assurance for inorganic and organic compounds assessed in the National Water Quality Network includes collection of field blanks to determine contamination bias and field replicates to determine variability bias. No contamination bias was found for 6 of the 13 nutrient compounds analyzed, and some potential contamination bias for some years was found for the other 7 nutrient compounds. Contamination bias was not found for carbon compounds or ultraviolet-absorbance measurements and was not assessed for sediment. All major ions and trace elements except potassium and lithium showed moderate contamination bias for at least 1 water year; generally, this bias was not at environmentally relevant concentrations. All compounds in the nutrient, carbon, and sediment group and in the major ions and trace elements group had low variability both in detection frequency and in concentration. Exceptions to this low variability were total particulate inorganic carbon and sediment for 2015, both of which are particulate substances with intrinsically high sampling variability. The risk of contamination bias for pesticides in National Water Quality Network samples was low, as indicated by very few detections in field blanks. Sixteen pesticide compounds showed potential contamination bias based on unexpected detections in third-party blind spikes (false-positive results for compounds that are not included in the spike mixture of a sample, where the identity as a QC sample is unknown to the analyst), and 47 different compounds (out of 225 pesticide compounds) showed potential contamination bias from laboratory blanks. However, when timing and relative magnitudes of detections in blank samples, environmental samples, and benchmark concentrations are considered, most of this potential contamination is not relevant to interpretation of published pesticide results. Overall variability in detection frequency for pesticides from field replicates was low or moderate. Also based on field replicates, 55 pesticides had overall high variability in concentrations for at least 1 water year, although these assessments likely overestimate high variability. At least 1 QC issue was found for 87 pesticides; however, most of the QC issues had no or little effect on the interpretation of environmental results because the U.S. Geological Survey National Water Quality Laboratory addressed the QC issue before publishing the environmental results, environmental results were almost entirely nondetections, concentrations of environmental results were higher than potential contamination bias, or benchmark concentrations were orders of magnitude higher than all environmental results. Eight compounds affected by two QC issues had a benchmark less than 100 nanograms per liter and warranted careful consideration of timing and magnitude of QC results in relation to surface-water results before interpretive use.

Scientific Investigations Report↗
USGS2023

Compressional-wave seismic velocity, bulk density, and their empirical relations for geophysical modeling of the Midcontinent Rift System in the Lake Superior region

Compressional-wave seismic velocity (velocity) and bulk density (density) data were compiled from published sources for rock suites and earth materials that are significant for geophysical modeling of the Mesoproterozoic Midcontinent Rift System in the Lake Superior region. The data include laboratory measurements of outcrop and drill core samples, seismic refraction studies, and a sonic log from a 1.5-kilometer-deep exploration well. Rock suites of the Midcontinent Rift System include basalts of the Mesoproterozoic Keweenawan Supergroup, Oronto Group sedimentary rocks (divided into arenaceous versus argillaceous units), and several sedimentary formations overlying the Oronto Group that have been correlated across the area. Intrusive units include diabase, gabbro, and felsic igneous rocks. Other geologic units important for geophysical modeling in the Lake Superior region include Archean crystalline crust, Paleoproterozoic metasedimentary and crystalline rocks, lower Mesoproterozoic sedimentary rocks, and Holocene to Pleistocene surficial deposits. Empirical velocity-density relations for each rock suite were determined by comparing the compiled data to published relations, such as the Nafe-Drake curve, Gardner’s relation, and best-fit equations developed for different rock types from laboratory studies. Graphical representations of these velocity-density relations provide a way to easily understand how velocity and density differ between tectonic settings and by rock type. Overlaps in velocity and density ranges for different geologic units are significant and have especially important implications for geologic interpretation of seismic data. Important examples include similar velocities but differing densities for argillaceous Oronto Group versus units overlying the Oronto Group and arenaceous Oronto Group versus basalt of the Keweenawan Supergroup. Similar densities but differing velocities were found for diabase versus gabbro. In addition, expected velocity ranges by rock type show that high-velocity intervals (6.9–7.1 kilometers per second) interpreted as basalt in previous seismic-reflection studies more likely indicate diabase or gabbro, suggesting that these interpretations may warrant additional consideration.

Lake Superior↗
USGS1959

Comparison of thickness, grade, and depth of radioactive layers as determined by gamma-ray logging and by core sampling

Thickness, grade, and depth data were obtained by analyzing gamma-ray logs and core samples from 56 diamond drill holes penetrating uranium deposits in the Colorado Plateau. The data from the two methods were compared to determine variations found in gamma-ray log interpretation and chemical and radiometric analyses of the drill core. Correlations within each parameter varied among the drilling areas analyzed. Gamma-ray interpretations of grade compared to chemical analyses were within the range of -10 to +25 percent. Most depth measurements determined by gamma-ray log interpretation compared to drill core measurement were within 0.5 percent. Results of the study indicate a need for better thickness definition in both gamma-ray logging and core scanning equipment.

Trace Elements Investigations↗
USGS2023

White-Nose Syndrome Diagnostic Laboratory Network handbook

When responding to a wildlife disease outbreak, managers depend on consistent and clear data to make decisions. However, diagnostic methods for detecting pathogens of wildlife often lack the level of procedural and interpretational standardization that occurs in the investigation of human and domestic animal diseases. This lack of standardization can hamper diagnostic reliability in two ways. First is the inappropriate application of tests to new species or in situations that are outside of the original (in other words, validated) purpose. Second is the use of laboratory-specific modifications or analytical parameters without thorough investigation of how those changes affect result comparisons across institutions or the ability to make broader conclusions about pathogen or disease. White-nose syndrome (WNS) is a disease caused by the fungal pathogen Pseudogymnoascus destructans ( Pd ), which has spread rapidly and is causing population-level declines in some species of North American bats. During the last decade, quantitative polymerase chain reaction (qPCR) has become the most common method of testing for Pd because of qPCR’s speed, accuracy, and simplicity across a wide range of invasive and noninvasive sample types. Its widespread use by many State, Federal, Provincial, and academic institutions has inevitably led to variations in methodology and interpretation among laboratories. The progressive geographic spread of fungus and disease has also led to sampling contexts and strategies that differ from those for which the qPCR assay was originally developed and validated. These factors have resulted in inconsistencies among results tested in different laboratories and, subsequently, confusion for managers and decision makers. To address these challenges, the WNS National Response Team Diagnostic Working Group launched a project congruent with increased calls for the harmonization of wildlife disease diagnostic results, and reporting standards across disparate methodologies and laboratories. Beginning in 2019, interlaboratory testing was done to better understand how variations to Pd qPCR methodology affect diagnostic consistency and to reassess the assay’s fit for purpose in new testing contexts. This information led to expanded conversations within the Diagnostic Working Group related to best practices in Pd qPCR diagnostic testing, the development of common interpretation language for classifying test results, and the incorporation of that language into an updated WNS case definition. This handbook is the resulting product and is intended to help further harmonize Pd qPCR diagnostic testing by establishing recommendations related to voluntary participation in a WNS Diagnostic Laboratory Network, documenting the currently (2022) practiced Pd qPCR methodologies, discussing general best practices for molecular diagnostics and laboratory networks, and elaborating on the epidemiologic and diagnostic basis of the agreed-upon classification language for Pd qPCR results. Through this voluntary, consensus-based approach to diagnostic harmonization, this work aims to improve the confidence of management agencies in reported Pd qPCR results and can serve as an example of national diagnostic coordination for other unregulated wildlife diseases.

Techniques and Methods↗
USGS2002

An integrated surface-geophysical investigation of the University of Connecticut landfill, Storrs, Connecticut: 2000

A surface-geophysical investigation to characterize the hydrogeology and contaminant distribution of the former landfill area at the University of Connecticut in Storrs, Connecticut, was conducted in 2000 to supplement the preliminary hydrogeologic assessment of the contamination of soil, surface water, and ground water at the site. A geophysical-toolbox approach was used to characterize the hydrogeology and contaminant distribution of the former landfill. Two-dimensional direct-current resistivity, inductive terrain-conductivity, and seismic-refraction surface-geophysical data were collected and interpreted in an iterative manner with exploratory drilling, borehole geophysics, and hydraulic testing. In this investigation, a geophysical-toolbox approach was used to 1) further define previously identified conductive anomalies and leachate plumes; 2) identify additional leachate plumes, possible fracture zones, and (or) conductive lithologic layers in the bedrock; and 3) delineate bedrock-surface topography in the drainage valleys north and south of the landfill. Resistivity and terrain-conductivity surveys were used to further delineate previously identified geophysical anomalies to the north and southwest of the landfill. A conductive anomaly identified in the terrain-conductivity survey to the north of the landfill in 2000 had a similar location and magnitude as an anomaly identified in terrain-conductivity surveys conducted in 1998 and 1999. Collectively, these surveys indicated that the magnitude of the conductive anomaly decreased with depth and with distance from the landfill. These anomalies indicated landfill leachate in the overburden and shallow bedrock. Results of previous surface-geophysical investigations southwest of the landfill indicated a shallow conductive anomaly in the overburden that extended into the fractured-bedrock aquifer. This conductive anomaly had a sheet-like geometry that had a north-south strike, dipped to the west, and terminated abruptly about 450 feet southwest of the landfill. The sheet-like conductive anomaly was interpreted as a fractured, conductive lithologic feature filled with conductive fluids. To further delineate this anomaly, two two-dimensional resistivity profiles were collected west of the sheet-like conductive anomaly to assess the possibility that the sheet-like conductive anomaly continued to the west in its down-dip direction. Each of the north-south oriented resistivity profiles showed bullet-shaped rather than linear-shaped anomalies, with a relatively smaller magnitude of conductivity than the sheet-like conductive anomaly to the east. If these bullet-like features are spatially connected, they may represent a linear, or pipe-like, conductive anomaly in the bedrock with a trend of N290°E and a plunge of 12°. Additional surveys were conducted to assess the apparent southern termination of the sheet-like conductive feature. Terrain-conductivity surveys indicated the sheet-like feature was not continuous to the south. A two-dimensional resistivity line and a coincident terrain-conductivity profile indicated the presence of a steep, eastward dipping, low magnitude, electrically conductive anomaly on the eastern end of the profile. Although the sheet-like conductive anomaly apparently did not continue to the south, the survey conducted in 2000 identified an isolated, weak conductive anomaly south of the previously identified anomaly. Inductive terrain-conductivity surveys performed north of the sheet-like conductive anomaly and west of the landfill indicated the anomaly did not extend to the north into the area of the former chemical-waste disposal pits. No conductive plumes or conductive features were observed in the subsurface bedrock west of the landfill. A conductive anomaly was identified in the southern section of the new terrain-conductivity grid. The magnitude and distribution of the apparent conductivity of this anomaly was identified as a nearly vertical sheet-like conductive feature. The anomaly extended north-south for 150 feet, and the depth and (or) magnitude of the conductive anomaly decreased towards the south. Based on the location, orientation, and dip, this feature was interpreted as a separate anomaly rather than as a continuation of the sheet-like, conductive bedrock feature that was previously identified southwest of the landfill. No other new conductive anomalies were identified south or west of the landfill. Seismic-refraction surveys were used to delineate the depth to the water table and the depth and topography of the bedrock surface. The seismic-refraction surveys at the northern and southern ends of the landfill confirm the presence of shallow bedrock at 25 feet below land surface. Seismic-refraction surveys conducted southwest of the landfill in a minor topographic valley indicate the bedrock is about 10 feet below land surface.

Connecticut↗
USGS2003

Quality-control results for ground-water and surface-water data, Sacramento River Basin, California, National Water-Quality Assessment, 1996-1998

Evaluating the extent that bias and variability affect the interpretation of ground- and surface-water data is necessary to meet the objectives of the National Water-Quality Assessment (NAWQA) Program. Quality-control samples used to evaluate the bias and variability include annual equipment blanks, field blanks, field matrix spikes, surrogates, and replicates. This report contains quality-control results for the constituents critical to the ground- and surface-water components of the Sacramento River Basin study unit of the NAWQA Program. A critical constituent is one that was detected frequently (more than 50 percent of the time in blank samples), was detected at amounts exceeding water-quality standards or goals, or was important for the interpretation of water-quality data. Quality-control samples were collected along with ground- and surface-water samples during the high intensity phase (cycle 1) of the Sacramento River Basin NAWQA beginning early in 1996 and ending in 1998. Ground-water field blanks indicated contamination of varying levels of significance when compared with concentrations detected in environmental ground-water samples for ammonia, dissolved organic carbon, aluminum, and copper. Concentrations of aluminum in surface-water field blanks were significant when compared with environmental samples. Field blank samples collected for pesticide and volatile organic compound analyses revealed no contamination in either ground- or surface-water samples that would effect the interpretation of environmental data, with the possible exception of the volatile organic compound trichloromethane (chloroform) in ground water. Replicate samples for ground water and surface water indicate that variability resulting from sample collection, processing, and analysis was generally low. Some of the larger maximum relative percentage differences calculated for replicate samples occurred between samples having lowest absolute concentration differences and(or) values near the reporting limit. Surrogate recoveries for pesticides analyzed by gas chromatography/mass spectrometry (GC/MS), pesticides analyzed by high performance liquid chromatography (HPLC), and volatile organic compounds in ground- and surface-water samples were within the acceptable limits of 70 to 130 percent and median recovery values between 82 and 113 percent. The recovery percentages for surrogate compounds analyzed by HPLC had the highest standard deviation, 20 percent for ground-water samples and 16 percent for surface-water samples, and the lowest median values, 82 percent for ground-water samples and 91 percent for surface-water samples. Results were consistent with the recovery results described for the analytical methods. Field matrix spike recoveries for pesticide compounds analyzed using GC/MS in ground- and surface-water samples were comparable with published recovery data. Recoveries of carbofuran, a critical constituent in ground- and surface-water studies, and desethyl atrazine, a critical constituent in the ground-water study, could not be calculated because of problems with the analytical method. Recoveries of pesticides analyzed using HPLC in ground- and surface-water samples were generally low and comparable with published recovery data. Other methodological problems for HPLC analytes included nondetection of the spike compounds and estimated values of spike concentrations. Recovery of field matrix spikes for volatile organic compounds generally were within the acceptable range, 70 and 130 percent for both ground- and surface-water samples, and median recoveries from 62 to 127 percent. High or low recoveries could be related to errors in the field, such as double spiking or using spike solution past its expiration date, rather than problems during analysis. The methodological changes in the field spike protocol during the course of the Sacramento River Basin study, which included decreasing the amount of spike solu

Water-Resources Investigations Report↗
USGS2003

Water resources of Monroe County, New York, water years 1997-99, with emphasis on water quality in the Irondequoit Creek basin—Atmospheric deposition, ground water, streamflow, trends in water quality, and chemical loads to Irondequoit Bay

Irondequoit Creek drains 169 square miles in the eastern part of Monroe County. Over time, nutrients transported by Irondequoit Creek to Irondequoit Bay on Lake Ontario have contributed to the eutrophication of the bay. Sewage-treatment-plant effluent, a major source of nutrients to the creek and its tributaries, was eliminated from the basin in 1979 by diversion to a regional wastewater-treatment facility, but sediment and contaminants from nonpoint sources continue to enter the creek and Irondequoit Bay. This report, the fourth in a series of reports that present interpretive analyses of the hydrologic data collected in Monroe County since 1984, interprets data from four surface-water monitoring sites in the Irondequoit Creek basin—Irondequoit Creek at Railroad Mills, East Branch Allen Creek at Pittsford, Allen Creek near Rochester, and Irondequoit Creek at Blossom Road. It also interprets data from three sites in the the Genesee River basin—Oatka Creek at Garbutt, Honeoye Creek at Honeoye Falls, and Black Creek at Churchville—as well as the Genesee River at Charlotte Pump Station, and also from a site on Northrup Creek at North Greece. The Northrup Creek site drains a 23.5-square-mile basin in western Monroe County, and provides information on surface-water quality in streams west of the Genesee River and on loads of nutrients delivered to Long Pond, a small eutrophic embayment of Lake Ontario. The report also includes water-level and water-quality data from nine observation wells in Ellison Park, and atmospheric-deposition data from a collection site at Mendon Ponds County Park. Average annual loads of some chemical constituents in atmospheric deposition for 1997–99 differed considerably from those for the long-term period 1984–96. Ammonia and potassium loads for 1997-99 were 144 and 118 percent greater, respectively, than for the previous period. Sodium and ammonia + organic nitrogen loads were 87 and 60 percent greater, respectively. Average annual loads of sulfate and orthophosphate for 1997-99 were 36 and 30 percent lower, respectively, than for the previous period. Loads of all nutrients deposited on the Irondequoit basin from atmospheric sources during 1997–99 greatly exceeded those transported by Irondequoit Creek. The ammonia load deposited on the basin was 139 times the load transported at Blossom Road (the most downstream site); the ammonia + organic nitrogen load was 6.3 times greater, orthophosphate 7.5 times greater, total phosphorus 1.3 times greater and nitrite + nitrate 1.5 times greater. Average yields of dissolved chloride and dissolved sulfate from atmospheric sources were much smaller than those transported by streamflow at Blossom Road.chloride was about 2 percent and sulfate about 8 percent of the amount transported. Trends in concentration of chemical constituents in surface water generally can be attributed to changes in land use, annual and seasonal variations in streamflow, and annual variations in the application of road salt to county highways and roads. Concentrations of several constituents in streams of the Irondequoit Creek basin showed statistically significant (α=0.05) trends from the beginning of their period of record through 1999. The constituent with the greatest number of significant trends was ammonia + organic nitrogen, with downward trends ranging from 4.1 to 5.6 percent per year at Allen Creek, Irondequoit Creek at Blossom Road, and East Branch Allen Creek. Orthophosphate showed an upward trend of 4.1 percent per year at Irondequoit Creek at Railroad Mills (the most upstream site). Dissolved chloride showed upward trends at Railroad Mills, Allen Creek, and Blossom Road. No trends in volatile suspended solids were noted at any of the four Irondequoit basin sites. Northrup Creek showed significant downward trends in concentrations of ammonia + organic nitrogen (3.3 percent per year), total phosphorus (3.4 percent per year), and orthophosphate (5.5 percent per year), and an upward trend for dissolved sulfate (1.8 percent per year). The Genesee River at Charlotte Pump Station showed downward trends of 6.1 percent per year for ammonia + organic nitrogen and 0.1 percent per year for chloride, and upward trends of 1.7 percent per year for total phosphorus and 6.6 percent per year for orthophosphate. Mean annual yields (mass per unit area) of most constituents at the Irondequoit Creek basin sites were similar to those noted for the previous report period (1994–96). East Branch Allen Creek showed lower yields of all constituents during 1997–99 than previously, even though runoff during 1997–99 was greater. These lower yields are attributed to the construction of an upstream detention basin on East Branch Allen Creek in 1995. Statistical analysis of long-term (greater than 12 years) streamflow records for unregulated streams in Monroe County indicated that annual mean flows for water years 1997–99 were in the normal range (75th to 25th percentile), although Allen Creek continues to show a significant downward trend in mean monthly streamflow during the 1984–99 water years.

New York↗
USGS2004

Characterization of aquifer heterogeneity using cyclostratigraphy and geophysical methods in the upper part of the Karstic Biscayne Aquifer, Southeastern Florida

This report identifies and characterizes candidate ground-water flow zones in the upper part of the shallow, eogenetic karst limestone of the Biscayne aquifer in the Lake Belt area of north-central Miami-Dade County using cyclostratigraphy, ground-penetrating radar (GPR), borehole geophysical logs, and continuously drilled cores. About 60 miles of GPR profiles were used to calculate depths to shallow geologic contacts and hydrogeologic units, image karst features, and produce qualitative views of the porosity distribution. Descriptions of the lithology, rock fabrics, and cyclostratigraphy, and interpretation of depositional environments of 50 test coreholes were linked to the geophysical interpretations to provide an accurate hydrogeologic framework. Molluscan and benthic foraminiferal paleontologic constraints guided interpretation of depositional environments represented by rockfabric facies. Digital borehole images were used to characterize and quantify large-scale vuggy porosity. Preliminary heat-pulse flowmeter data were coupled with the digital borehole image data to identify candidate ground-water flow zones. Combined results show that the porosity and permeability of the karst limestone of the Biscayne aquifer have a highly heterogeneous and anisotropic distribution that is mostly related to secondary porosity overprinting vertical stacking of rock-fabric facies within high-frequency cycles (HFCs). This distribution of porosity produces a dual-porosity system consisting of diffuse-carbonate and conduit flow zones. The nonuniform ground-water flow in the upper part of the Biscayne aquifer is mostly localized through secondary permeability, the result of solution-enlarged carbonate grains, depositional textures, bedding planes, cracks, root molds, and paleokarst surfaces. Many of the resulting pore types are classified as touching vugs. GPR, borehole geophysical logs, and whole-core analyses show that there is an empirical relation between formation porosity, permeability, formation electrical conductivity, and GPR reflection amplitudes? as porosity and permeability increase, formation electrical conductivity increases and reflection amplitude decreases. This relation was observed throughout the entire vertical and lateral section of the upper part of the Biscayne aquifer in the study area. Further, upward-shallowing brackish- or freshwatercapped cycles of the upper part of the Fort Thompson Formation show low-amplitude reflections near their base that correspond to relatively higher porosity and permeability. This distribution is related to a systematic vertical stacking of rock-fabric facies within the cycle. Inferred flow characteristics of the porosity distribution within the upper part of the Biscayne aquifer were used to identify four ground-water flow classes, with each characterized by a discrete pore system that affects vertical and horizontal groundwater flow: (1) a low-permeability peat, muck, and marl ground-water flow class; (2) a horizontal conduit ground-water flow class; (3) a leaky, low-permeability ground-water flow class; and (4) a diffuse-carbonate ground-water flow class. At the top of the Biscayne aquifer, peat, muck, and marl can combine to form a relatively low-permeability layer of Holocene sediment that water moves through slowly. Most horizontal conduit flow is inferred to occur along touching vugs in portions of the following rock-fabric facies: (1) touchingvug pelecypod floatstone and rudstone, (2) sandy touching-vug pelecypod floatstone and rudstone, (3) vuggy wackestone and packstone, (4) laminated peloid grainstone and packstone, (5) peloid grainstone and packstone, and (6) peloid wackestone and packstone. Gastropod floatstone and rudstone, mudstone and wackestone, and pedogenic limestone rock-fabric facies are the main hosts for leaky, low-permeability units.

Florida↗
USGS2003

Possible Extent and Depth of Salt Contamination in Ground Water Using Geophysical Techniques, Red River Aluminum Site, Stamps, Arkansas, April 2003

A surface-geophysical investigation of the Red River Aluminum site at Stamps, Arkansas, was conducted in cooperation with the Arkansas Department of Environmental Quality to determine the possible extent and depth of saltwater contamination. Water-level measurements indicate the distance to water level below land surface ranges from about 1.2 to 3.9 feet (0.37 to 1.19 meters) in shallow monitor wells and about 10.5 to 17.1 feet (3.20 to 5.21 meters) in deeper monitoring wells. The two-dimensional, direct-current resistivity method identified resistivities less than 5 ohm-meters which indicated possible areas of salt contamination occurring in near-surface or deep subsurface ground water along four resistivity lines within the site. One line located east of the site yielded data that demonstrated no effect of salt contamination. Sections from two of the five data sets were modeled. The input model grids were created on the basis of the known geology and the results and interpretations of borehole geophysical data. The clay-rich Cook Mountain Formation is modeled as 25 ohm-meters and extends from 21 meters (68.9 feet) below land surface to the bottom of the model (about 52 meters (170.6 feet)). The models were used to refine interpretation of the resistivity data and to determine extent of saltwater contamination and depth to the Cook Mountain Formation. Data from the resistivity lines indicate both near-surface and subsurface saltwater contamination. The near-surface contamination appears as low resistivity (less than 5 ohm-meters) on four of the five resistivity lines, extending up to 775 meters (2,542.8 feet) horizontally in a line that traverses the entire site south to north. Model resistivity data indicate that the total depth of saltwater contamination is about 18 meters (59 feet) below land surface. Data from four resistivity lines identified areas containing low resistivity anomalies interpreted as possible salt contamination. A fifth line located just east of the site showed no saltwater contamination.

Water-Resources Investigations Report↗
USGS1985

Management of ground water and evolving hydrogeologic studies in New Jersey : a heavily urbanized and industrialized state in the northeastern United States

New Jersey is the most densely populated and one of the most industrialized states in the United States. An abundance of freshwater and proximity to major northeastern metropolitan centers has facilitated this development. Pumpage of freshwater from all aquifers in the State in 1980 was 730 million gallons per day (2.76 million cubic meters per day). Management and efficient development of the ground-water resources of the State are the responsibility of the New Jersey Department of Environmental Protection. Laws have been enacted and updated by the State legislature to manage water allocation and to control the disposal of hazardous wastes. Present resource management is guided by the New Jersey Water-Supply Master Plan of 1981. Funding for management activities is partially derived from the sale of state-approved bonds. Effective planning and regional management require accurate and up-to-date hydrologic information and analyses. The U.S. Geological Survey, in cooperation with the New Jersey Geological Survey, is conducting three intensive ground-water studies involving the collection and interpretation of hydrologic data to meet the urgent water-management needs of New Jersey. These studies are part of a long-term cooperative program and are funded through the Water-Supply Bond Act of 1981. They began in 1983 and are scheduled to be completed in 1988. The project areas are situated in the New Jersey part of the Atlantic Coastal Plain in and near Atlantic City, Camden, and South River. They range in size from 400 to 1,200 mil (1,040 to 3,120 km2). The studies are designed to define the geology, hydrology, and geochemistry of the local ground-water systems. The results of these studies will enable the State to address more effectively major problems in these areas such as declining water levels, overpumping, saltwater intrusion, and ground-water contamination resulting from the improper disposal of hazardous wastes. Specific objectives of these studies by the U.S. Geological Survey are to (1) develop an accurate and up-to-date hydrogeologic data base, (2) design and implement a data-collection program and establish a computerized information management system, (3) refine the conceptualization of the ground-water flow system, and (4) define the geochemistry of the aquifer system by conducting a water-quality appraisal. The objectives are accomplished by standard hydrogeologic methods. Information concerning hydrogeologic framework, ground-water levels, water use, hydraulic characteristics, and water quality in the study areas is compiled from all available sources. Additional data needed are collected through well inventories, surface geophysical surveys, water-quality samplings, water-level measurements, and a well-drilling program. Interpretation of the flow system is based on the use of standard analytical techniques and digital flow modeling. Calibrated flow models will provide ground-water managers with a mechanism to develop and test regional water-supply strategies. Definition of the geochemistry of the aquifer system is accomplished through a variety of methods which depend on the problems and available data in the particular study area. The approach includes statistical analysis of water-quality data, reaction-path modeling, and determination of the movement of chemical constituents using analytical and numerical modeling techniques. A combined staff of 25 to 30 professionals and technicians from the New Jersey District office of the U.S. Geological Survey is committed to the three studies. The staff has specialists in geohydrology, numerical modeling, geochemistry, geophysics, and computer science. The findings of these studies will be published in data reports, interpretive reports, instructional manuals and journal articles.

New Jersey↗
USGS1993

Application of ground-penetrating radar methods in determining hydrogeologic conditions in a karst area, west-central Florida

Ground-penetrating radar (GPR) is useful as a surface geophysical method for exploring geology and subsurface features in karst settings. Interpretation of GPR data was used to infer lithology and hydrogeologic conditions in west-central Florida. This study demonstrates how GPR methods can be used to investigate the hydrogeology of an area. GPR transmits radio- frequency electromagnetic waves into the ground and receives reflected energy waves from subsurface interfaces. Subsurface profiles showing sediment thickness, depth to water table and clay beds, karst development, buried objects, and lake-bottom structure were produced from GPR traverses obtained during December 1987 and March 1990 in Pinellas, Hillsborough, and Hardee Counties in west-central Florida. Performance of the GPR method is site specific, and data collected are principally affected by the sediment and pore fluids, conductances and dielectric constants. Effective exploration depths of the GPR surveys through predominately unsaturated and saturated sand and clay sediments at five study sites ranged from a few feet to greater than 50 feet below land surface. Exploration depths were limited when high conductivity clay was encountered, whereas greater exploration depths were possible in material composed of sand. Application of GPR is useful in profiling subsurface conditions, but proper interpretation depends upon the user's knowledge of the equipment and the local hydrogeological setting, as well as the ability to interpret the graphic profile.

Water-Resources Investigations Report↗
USGS1995

Analysis of steady-state flow and advective transport in the eastern Snake River Plain aquifer system, Idaho

Quantitative estimates of ground-water flow directions and traveltimes for advective flow were developed for the regional aquifer system of the eastern Snake River Plain, Idaho. The work included: (1) descriptions of compartments in the aquifer that function as intermediate and regional flow systems, (2) descriptions of pathlines for flow originating at or near the water table, and (3) quantitative estimates of traveltimes for advective transport originating at or near the water table. A particle-tracking postprocessing program was used to compute pathlines on the basis of output from an existing three-dimensional steady-state flow model. The flow model uses 1980 conditions to approximate average annual conditions for 1950-80. The advective transport model required additional information about the nature of flow across model boundaries, aquifer thickness, and porosity. Porosity of two types of basalt strata has been reported for more than 1,500 individual cores from test holes, wells, and outcrops near the south side of the Idaho National Engineering Laboratory. The central 80 percent of samples had porosities of 0.08 to 0.25, the central 50 percent of samples, O. 11 to 0.21. Calibration of the model involved choosing a value for porosity that yielded the best solution. Two radiologic contaminants, iodine-129 and tritium, both introduced to the flow system about 40 years ago, are relatively conservative tracers. Iodine- 129 was considered to be more useful because of a lower analytical detection limit, longer half-life, and longer flow path. The calibration value for porosity was 0.21. Most flow in the aquifer is contained within a regional-scale compartment and follows paths that discharge to the Snake River downstream from Milner Dam. Two intermediate-scale compartments exist along the southeast side of the aquifer and near Mud Lake.One intermediate-scale compartment along the southeast side of the aquifer discharges to the Snake River near American Fails Reservoir and covers an area of nearly 1,000 square miles. This compartment, which receives recharge from an area of intensive surface-water irrigation, is apparently fairly stable. The other intermediate-scale compartment near Mud Lake covers an area of 300 square miles. The stability and size of this compartment are uncertain, but are assumed to be in a state of change. Traveltimes for advective flow from the water table to discharge points in the regional compartment ranged from 12 to 350 years for 80 percent of the particles; in the intermediate-scale flow compartment near American Falls Reservoir, from 7 to 60 years for 80 percent of the particles; and in the intermediate-scale compartment near Mud Lake, from 25 to 100 years for 80 percent of the particles. Traveltimes are sensitive to porosity and assumptions regarding the importance of the strength of internal sinks, which represent ground-water pumpage. A decrease in porosity results in shorter traveltimes but not a uniform decrease in traveltime, because the porosity and thickness is different in each model layer. Most flow was horizontal and occurred in the top 500 feet of the aquifer. An important limitation of the model is the assumption of steady-state flow. The most recent trend in the flow system has been a decrease in recharge since 1987 because of an extended drought and changes in land use. A decrease in flow through the system will result in longer traveltimes than those predicted for a greater flow. Because the interpretation of the model was limited to flow on a larger scale, and did not consider individual wells or well fields, the interpretations were not seriously limited by the discretization of well discharge. The interpretations made from this model also were limited by the discretization of the major discharge areas. Near discharge areas, pathlines might not be representative at the resolution of the grid. Most improvement in the estimates of ground-waterflow directions and travelt

Water-Resources Investigations Report↗
USGS1997

Review of selected features of the natural system model, and suggestions for applications in south Florida

A study was conducted to review selected features of the Natural System Model, version 4.3. The Natural System Model is a regional-scale model that uses recent climatic data and estimates of historic vegetation and topography to simulate pre-canal-drainage hydrologic response in south Florida. Equations used to represent the hydrologic system and the numerical solution of these equations in the model were documented and reviewed. Convergence testing was performed using 1965 input data, and selected other aspects of the model were evaluated. Some conclusions from the evaluation of the Natural System Model include the following observations. Simulations were generally insensitive to the temporal resolution used in the model. However, reduction of the computational cell size from 2-mile by 2-mile to 2/3-mile by 2/3-mile resulted in a decrease in spatial mean ponding depths for October of 0.35 foot for a 3-hour time step. Review of the computer code indicated that there is no limit on the amount of water that can be transferred from the river system to the overland flow system, on the amount of seepage from the river to the ground-water system, on evaporation from the river system, or on evapotranspiration from the overland-flow system. Oscillations of 0.2 foot or less in simulated river stage were identified and attributed to a volume limiting function which is applied in solution of the overland-flow equations. The computation of the resistance coefficient is not consistent with the computation of overland-flow velocity. Ground-water boundary conditions do not always ensure a no-flow condition at the boundary. These inconsistencies had varying degrees of effects on model simulations, and it is likely that simulations longer than 1 year are needed to fully identify effects. However, inconsistencies in model formulations should not be ignored, even if the effects of such errors on model results appear to be small or have not been clearly defined. The Natural System Model can be a very useful tool for estimating pre-drainage hydrologic response in south Florida. The model includes all of the important physical processes needed to simulate a water balance. With a few exceptions, these hydrologic processes are represented in a reasonable manner using empirical, semiempirical, and mechanistic relations. The data sets that have been assembled to represent physical features, and hydrologic and meteorological conditions are quite extensive in their scope. Some suggestions for model application were made. Simulation results from the Natural System Model need to be interpreted on a regional basis, rather than cell by cell. The available evidence suggests that simulated water levels should be interpreted with about a plus or minus 1 foot uncertainty. It is probably not appropriate to use the Natural System Model to estimate pre-drainage discharges (as opposed to hydroperiods and water levels) at a particular location or across a set of adjacent computational cells. All simulated results for computational cells within about 10 miles of the model boundaries have a higher degree of uncertainty than results for the interior of the model domain. It is most appropriate to interpret the Natural System Model simulation results in connection with other available information. Stronger linkages between hydrologic inputs to the Everglades and the ecological response of the system would enhance restoration efforts.

Florida↗