ISOEDIT : a program for interactive editing of seismic reflection velocity models computer file, Part A
Explore the source record for details and available documents.
SEARCH · Geology Reports
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Statistical operations termed model-adjustment procedures (MAP’s) can be used to incorporate local data into existing regression models to improve the prediction of urban-runoff quality. Each MAP is a form of regression analysis in which the local data base is used as a calibration data set. Regression coefficients are determined from the local data base, and the resulting ‘adjusted’ regression models can then be used to predict storm-runoff quality at unmonitored sites. The response variable in the regression analyses is the observed load or mean concentration of a constituent in storm runoff for a single storm. The set of explanatory variables used in the regression analyses is different for each MAP, but always includes the predicted value of load or mean concentration from a regional regression model. The four MAP’s examined in this study were: single-factor regression against the regional model prediction, P, (termed MAP-lF-P), regression against P,, (termed MAP-R-P), regression against P, and additional local variables (termed MAP-R-P+nV), and a weighted combination of P, and a local-regression prediction (termed MAP-W). The procedures were tested by means of split-sample analysis, using data from three cities included in the Nationwide Urban Runoff Program: Denver, Colorado; Bellevue, Washington; and Knoxville, Tennessee. The MAP that provided the greatest predictive accuracy for the verification data set differed among the three test data bases and among model types (MAP-W for Denver and Knoxville, MAP-lF-P and MAP-R-P for Bellevue load models, and MAP-R-P+nV for Bellevue concentration models) and, in many cases, was not clearly indicated by the values of standard error of estimate for the calibration data set. A scheme to guide MAP selection, based on exploratory data analysis of the calibration data set, is presented and tested. The MAP’s were tested for sensitivity to the size of a calibration data set. As expected, predictive accuracy of all MAP’s for the verification data set decreased as the calibration data-set size decreased, but predictive accuracy was not as sensitive for the MAP’s as it was for the local regression models.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
A daily water-budget model for estimating ground-water recharge, the Deep Percolation Model, was modularized for inclusion into the U.S. Geological Survey's Modular Modeling System. The model was modularized in order to facilitate estimation of ground-water recharge under a large range in climatic, landscape, and land-use and land-cover conditions. The model can be applied to areas as large as regions or as small as a field plot. An overview of the Modular Modeling System and the Deep Percolation Model is presented. Data requirements, parameters, and variables for the model are described. The modules that compose the Deep Percolation Model are documented.
The transient response of a hydrologic system can be of concern to water-resource managers, because it is often extreme relatively short-lived events, such as floods or droughts, that profoundly influence the management of the resource. The water available to a hydrologic system for stream flow and aquifer recharge is determined by the difference of precipitation and evapotranspiration (ET). As such, temporal variations in precipitation and ET determine the degree of influence each has on the transient response of the hydrologic system. Meteorological, ET, and hydrologic data collected from 1993 to 2003 and spanning 1- to 3 2/3 -year periods were used to develop a hydrologic model for each of five sites in central Florida. The sensitivities of simulated water levels and flows to simple approximations of ET were quantified and the adequacy of each ET approximation was assessed. ET was approximated by computing potential ET, using the Hargreaves and Priestley-Taylor equations, and applying vegetation coefficients to adjust the potential ET values to actual ET. The Hargreaves and Priestley-Taylor ET approximations were used in the calibrated hydrologic models while leaving all other model characteristics and parameter values unchanged. Two primary factors that influence how the temporal variability of ET affects hydrologic simulation in central Florida were identified: (1) stochastic character of precipitation and ET and (2) the ability of the local hydrologic system to attenuate variability in input stresses. Differences in the stochastic character of precipitation and ET, both the central location and spread of the data, result in substantial influence of precipitation on the quantity and timing of water available to the hydrologic system and a relatively small influence of ET. The temporal variability of ET was considerably less than that of precipitation at each site over a wide range of time scales (from daily to annual). However, when precipitation and ET are of similar magnitude, small errors in ET can produce relatively large errors in available water, and accurate estimates of actual ET are more important. Local hydrologic conditions can also be an important factor influencing the hydrologic response to ET variability. Various points along a flow path in a hydrologic system respond differently to temporal variations in ET. For example, soil moisture contents in the root zone are sensitive to daily variations in ET, whereas spring flow responds to only longer term variations in ET. Both the Hargreaves and Priestley-Taylor equations for potential ET, when applied with an annually invariant monthly vegetation coefficient derived from comparison of measured ET with computed potential ET values, can be used with a hydrologic model to produce reasonable predictions of water levels and flows. Baseline-adjusted modified coefficients of efficiency for simulated water levels ranged from 0.0, indicating that water levels were simulated equally as well with approximated ET as with actual ET values, to -0.6, indicating that water levels were simulated better with actual ET values. Simulations using the Hargreaves approximation consistently yielded larger absolute and relative errors than the Priestley-Taylor approximation. However, the differences between the Hargreaves and Priestley-Taylor simulations generally were much smaller than differences between these simulations and the simulations using actual ET. This suggests that the simpler Hargreaves equation may be an adequate substitute for the more complex Priestley-Taylor equation, depending on the level of accuracy required to satisfy the particular modeling objectives.
In 2014, the U.S. Geological Survey, in cooperation with the U.S. Department of Defense’s Strategic Environmental Research and Development Program, initiated a project to evaluate the potential impacts of projected climate-change on Department of Defense installations that rely on Guam’s water resources. A major task of that project was to develop a watershed model of southern Guam and a water-balance model for the Fena Valley Reservoir. The southern Guam watershed model provides a physically based tool to estimate surface-water availability in southern Guam. The U.S. Geological Survey’s Precipitation Runoff Modeling System, PRMS-IV, was used to construct the watershed model. The PRMS-IV code simulates different parts of the hydrologic cycle based on a set of user-defined modules. The southern Guam watershed model was constructed by updating a watershed model for the Fena Valley watersheds, and expanding the modeled area to include all of southern Guam. The Fena Valley watershed model was combined with a previously developed, but recently updated and recalibrated Fena Valley Reservoir water-balance model. Two important surface-water resources for the U.S. Navy and the citizens of Guam were modeled in this study; the extended model now includes the Ugum River watershed and improves upon the previous model of the Fena Valley watersheds. Surface water from the Ugum River watershed is diverted and treated for drinking water, and the Fena Valley watersheds feed the largest surface-water reservoir on Guam. The southern Guam watershed model performed “very good,” according to the criteria of Moriasi and others (2007), in the Ugum River watershed above Talofofo Falls with monthly Nash-Sutcliffe efficiency statistic values of 0.97 for the calibration period and 0.93 for the verification period (a value of 1.0 represents perfect model fit). In the Fena Valley watershed, monthly simulated streamflow volumes from the watershed model compared reasonably well with the measured values for the gaging stations on the Almagosa, Maulap, and Imong Rivers—tributaries to the Fena Valley Reservoir—with Nash-Sutcliffe efficiency values of 0.87 or higher. The southern Guam watershed model simulated the total volume of the critical dry season (January to May) streamflow for the entire simulation period within –0.54 percent at the Almagosa River, within 6.39 percent at the Maulap River, and within 6.06 percent at the Imong River. The recalibrated water-balance model of the Fena Valley Reservoir generally simulated monthly reservoir storage volume with reasonable accuracy. For the calibration and verification periods, errors in end-of-month reservoir-storage volume ranged from 6.04 percent (284.6 acre-feet or 92.7 million gallons) to –5.70 percent (–240.8 acre-feet or –78.5 million gallons). Monthly simulation bias ranged from –0.48 percent for the calibration period to 0.87 percent for the verification period; relative error ranged from –0.60 to 0.88 percent for the calibration and verification periods, respectively. The small bias indicated that the model did not consistently overestimate or underestimate reservoir storage volume. In the entirety of southern Guam, the watershed model has a “satisfactory” to “very good” rating when simulating monthly mean streamflow for all but one of the gaged watersheds during the verification period. The southern Guam watershed model uses a more sophisticated climate-distribution scheme than the older model to make use of the sparse climate data, as well as includes updated land-cover parameters and the capability to simulate closed depression areas. The new Fena Valley Reservoir water-balance model is useful as an updated tool to forecast short-term changes in the surface-water resources of Guam. Furthermore, the now spatially complete southern Guam watershed model can be used to evaluate changes in streamflow and recharge owing to climate or land-cover changes. These are substantial improvements to the previous models of the Fena Valley watershed and Reservoir. Datasets associated with this report are available as a U.S. Geological Survey data release (Rosa and Hay, 2017; DOI:10.5066/F7HH6HV4).
The San Antonio Creek Valley watershed (SACVW) is located in western Santa Barbara County, about 15 miles south of Santa Maria and 55 miles north of Santa Barbara, California. The SACVW is about 135 square miles and encompasses the San Antonio Creek Valley groundwater basin; the SACVW is separated from adjacent groundwater basins by the Casmalia and Solomon Hills to the north, and the Purisima Hills to the south. At the western, downstream part of the valley, uplifted, consolidated rocks cause groundwater to discharge at land surface at Barka Slough. Since the late 1800s, groundwater has been the primary source of water for agricultural, military, municipal, and domestic uses. Groundwater withdrawal by pumping exceeded the amount of water replenishing the aquifer system during water years 1948–2018, causing groundwater-level declines of more than 150 feet in parts of the valley and reducing base flow at Barka Slough. Reliance on groundwater for agricultural water use (primarily for the irrigation and frost protection of vineyards, and fruit and berry crops) continues to strain the sustainability of the groundwater system. Through a cooperative agreement, the Santa Barbara County Water Agency and Vandenberg Space Force Base invited the U.S. Geological Survey to address declines in groundwater levels, develop a better understanding of the hydrogeologic system, and provide tools to help evaluate and manage the effects of future development of the San Antonio Creek Valley groundwater basin within the encompassing San Antonio Creek Valley watershed (SACVW). The objectives of this study were to (1) refine the hydrogeologic framework of the San Antonio Creek Valley watershed, (2) quantify the hydrologic budget of the valley, and (3) develop hydrologic modeling tools to evaluate and aid in managing the groundwater resource. This report focuses on the first and second objectives to construct a hydrogeologic framework and characterize the historical and present-day hydrologic conditions of the SACVW during water years 1948–2018. As part of the second objective, work included quantifying the hydrologic budget and evaluating the hydrogeologic system using a combination of existing data and geologic and hydrologic data collected for this study. The groundwater-flow system in the SACVW consists of five hydrogeologic units. These separate water-bearing units were identified based on hydrogeologic properties, such as sediment grain size, vertical-head differences in multiple-depth, monitoring-well sites, long-term groundwater level responses to pumping and climate, and the chemical character of groundwater and groundwater age in the mostly semi-consolidated to unconsolidated basin-fill sediments. The hydrogeologic units that comprise the different aquifers vary in their lithologic composition. The upper and lower aquifers (upper Paso Robles Formation, and lower Paso Robles Formation and Careaga Sandstone, respectively) are relatively coarse grained and are comprised of sand, gravel, and clay; the middle confining unit (the middle Paso Robles Formation) is relatively fine grained and is comprised of primarily clay, silt, and sand. The Pezzoni-Casmalia and Los Alamos faults, which are inferred to transect the SACVW between the western and eastern areas of the valley floor, do not appear to substantially affect the groundwater system. Present-day recharge to the study area occurs primarily as infiltration from precipitation and streams in the upland areas of the Casmalia Hills and Solomon Hills, and along the main channel of San Antonio Creek. Reported estimates of annual natural recharge during water years 1948–2018 generally ranged from about 5,000 acre-feet to more than about 30,000 acre-feet. Stable and radioactive isotopes show that groundwater from the lower aquifer is old and probably was recharged as infiltration from precipitation and streams in the eastern upland areas of the Solomon Hills; however, the infiltration and recharge from these sources probably does not occur under present-day climatic conditions. Anthropogenic recharge, from sources such as return flow from agricultural irrigation, municipal water systems, and wastewater effluent, was estimated to range from about 600 acre-feet in 1948 to about 6,600 acre-feet in 2018. The average annual amount of groundwater removed from the SACVW by pumping during 1948–2018 was estimated to be about 17,200 acre-feet per year, increasing from about 3,000 acre-feet in 1948 to about 32,600 acre-feet in 2018. Estimates of annual pumpage generally exceeded estimates of annual recharge beginning in the mid-1970s and continuing through 2018. The predominant direction of groundwater flow under historical and present-day conditions was from the eastern uplands in the Solomon Hills to the west along San Antonio Creek to the discharge area in Barka Slough, and from the northern uplands in the Casmalia Hills south to San Antonio Creek. Pumpage since the early 1900s and the subsequent groundwater-level declines have substantially reduced the amount of natural groundwater discharge at Barka Slough. Estimates of base flow to San Antonio Creek at the western, downstream extent of the SACVW have varied over time in response to changes in groundwater pumpage and climate; however, there was an overall decline in base flow during water years 1956–2018, decreasing from an average of about 1,700 acre-feet per year during 1956–69, to about 300 acre-feet per year during 2016–18. The long-term extraction of groundwater correlates with a decrease in groundwater levels by more than about 150 feet since the early 1940s in the eastern part of the basin near Los Alamos, and as much as about 50 feet in the upland areas and in the western part of the basin. At Barka Slough, groundwater levels have declined below land surface in some places, altering native riparian vegetation in and around the slough. Surface-water quality in the SACVW varied depending on location and the time of year the samples were collected and on the amount of annual precipitation Most groundwater in the SACVW was calcium-bicarbonate-type water with total dissolved-solids concentrations of about 500–800 milligrams per liter generally representing water naturally recharged as infiltration from precipitation and streams. Total dissolved-solids concentrations in some wells ranged from 800 to 8,000 milligrams per liter, suggesting mixing of naturally recharged infiltrated water with water associated with oil-bearing geologic formations, agricultural products, or the evaporation of shallow groundwater. Concentrations of total dissolved solids and the chemical constituents chloride, nitrate plus nitrite (as nitrogen), calcium, and magnesium at selected wells generally increased during water years 1980–2018; increasing concentrations of these constituents may be associated with the expansion of agriculture in the watershed over time and the corresponding increase in the use of nitrates and calcium- and magnesium-based fertilizers and soil additives in modern agricultural practices. The predominant direction of groundwater flow during historical and present-day conditions was from the eastern uplands in the Solomon Hills to the west along San Antonio Creek toward Barka Slough, and from the western uplands in the Casmalia Hills south to San Antonio Creek. The age of groundwater in the SACVW was evaluated using radioactive isotopes, and the flow of groundwater within the SACVW was evaluated using radioactive and stable isotopes. Modern groundwater (recharged after 1952) was generally found adjacent to San Antonio Creek and its tributaries in wells with perforated depths that averaged about 270 feet below land surface. Pre-modern groundwater (recharged before 1952) was found in wells that had average perforation depths of about 540 ft below land surface. Pre-modern groundwater identified in wells in the eastern upland area is interpreted to have had long, slow travel times to the western part of the SACVW where it was eventually discharged as base flow at Barka Slough or extracted as groundwater pumpage.
This report describes modifications to a U.S. Geological Survey (USGS) three-dimensional solute-transport model (MODFLOW-GWT), which is incorporated into the USGS MODFLOW ground-water model as the Ground-Water Transport (GWT) Process. The modifications were made to create compatibility between the Multi-Node Well (MNW) Package for MODFLOW and the MODFLOW-GWT model. This compatibility improves the capability of MODFLOW-GWT to represent accurately solute transport in simulations that include multi-node wells because long-screen wells or long open boreholes that extend through multiple model layers can provide fast pathways for solutes to move from one location to another in a ground-water flow system. For nonpumping multi-node wells (used to simulate open boreholes or observation wells, for example), a simple routing and local mixing model was developed to calculate nodal concentrations within the borehole. A depth-averaged concentration is calculated for such nonpumping wells. For pumping multi-node wells (either withdrawal or injection) in which the flow between the well and the ground-water system is in the same direction at all nodes, the average concentration in the well is calculated as a flux-based mean assuming complete and instantaneous mixing in the wellbore of all inflows. For pumping multi-node wells (either withdrawal or injection) in which the flow between the well and the ground-water system is not unidirectional, the concentration distribution within the well is calculated using the same routing and local mixing model used for a nonpumping multi-node well, with the added assumption that the flux pumped in or out of the well is added or removed above the first well node.
Percolation of precipitation through unsaturated zones is important for recharge of ground water. Rain and snowmelt at land surface are partitioned into different pathways including runoff, infiltration, evapotranspiration, unsaturated-zone storage, and recharge. A new package for MODFLOW-2005 called the Unsaturated-Zone Flow (UZF1) Package was developed to simulate water flow and storage in the unsaturated zone and to partition flow into evapotranspiration and recharge. The package also accounts for land surface runoff to streams and lakes. A kinematic wave approximation to Richards? equation is solved by the method of characteristics to simulate vertical unsaturated flow. The approach assumes that unsaturated flow occurs in response to gravity potential gradients only and ignores negative potential gradients; the approach further assumes uniform hydraulic properties in the unsaturated zone for each vertical column of model cells. The Brooks-Corey function is used to define the relation between unsaturated hydraulic conductivity and water content. Variables used by the UZF1 Package include initial and saturated water contents, saturated vertical hydraulic conductivity, and an exponent in the Brooks-Corey function. Residual water content is calculated internally by the UZF1 Package on the basis of the difference between saturated water content and specific yield. The UZF1 Package is a substitution for the Recharge and Evapotranspiration Packages of MODFLOW-2005. The UZF1 Package differs from the Recharge Package in that an infiltration rate is applied at land surface instead of a specified recharge rate directly to ground water. The applied infiltration rate is further limited by the saturated vertical hydraulic conductivity. The UZF1 Package differs from the Evapotranspiration Package in that evapotranspiration losses are first removed from the unsaturated zone above the evapotranspiration extinction depth, and if the demand is not met, water can be removed directly from ground water whenever the depth to ground water is less than the extinction depth. The UZF1 Package also differs from the Evapotranspiration Package in that water is discharged directly to land surface whenever the altitude of the water table exceeds land surface. Water that is discharged to land surface, as well as applied infiltration in excess of the saturated vertical hydraulic conductivity, may be routed directly as inflow to specified streams or lakes if these packages are active; otherwise, this water is removed from the model. The UZF1 Package was tested against the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model for a vertical unsaturated flow problem that includes evapotranspiration losses. This report also includes an example in which MODFLOW-2005 with the UZF1 Package was used to simulate a realistic surface-water/ground-water flow problem that includes time and space variable infiltration, evapotranspiration, runoff, and ground-water discharge to land surface and to streams. Another simpler problem is presented so that the user may use the input files as templates for new problems and to verify proper code installation.
A database called the “Comprehensive Resource Database” (CRD) was prepared to support U.S. Geological Survey (USGS) assessments of technically recoverable hydrocarbons that might result from the injection of miscible or immiscible carbon dioxide (CO 2 ) for enhanced oil recovery (EOR). The CRD was designed by INTEK Inc., a consulting company under contract to the USGS. The CRD contains data on the location, key petrophysical properties, production, and well counts (number of wells) for the major oil and gas reservoirs in onshore areas and State waters of the conterminous United States and Alaska. The CRD includes proprietary data on petrophysical properties of fields and reservoirs from the “Significant Oil and Gas Fields of the United States Database,” prepared by Nehring Associates in 2012, and proprietary production and drilling data from the “Petroleum Information Data Model Relational U.S. Well Data,” prepared by IHS Inc. in 2012. This report describes the CRD and the computer algorithms used to (1) estimate missing reservoir property values in the Nehring Associates (2012) database, and to (2) generate values of additional properties used to characterize reservoirs suitable for miscible or immiscible CO 2 flooding for EOR. Because of the proprietary nature of the data and contractual obligations, the CRD and actual data from Nehring Associates (2012) and IHS Inc. (2012) cannot be presented in this report.
During 1997–98, the U.S. Geological Survey, in cooperation with the Houston-Galveston Area Council, collected stream-habitat and benthic macroinvertebrate data for 31 reaches on abovetidal streams in the Council service area near Houston, Texas. Stream-habitat, land-use and population, and benthic aquatic insect metrics were determined for the 31 reaches. Statistical analyses were used to determine the stream-habitat, land-use and population, and aquatic insect variables that are strongly intercorrelated and that explain the greatest amount of variation between the reaches. Comparison of stream-habitat and biological integrity scores computed for each of the 31 reaches indicated (1) reaches generally had larger stream-habitat integrity scores in drainage areas that were heavily forested and had fewer people per square mile, (2) larger biological integrity scores were significantly correlated with larger stream-habitat integrity scores, and (3) urban reaches generally had more simplified streamhabitat conditions and smaller biological integrity scores. Seven reaches in the study area were selected as reference reaches on the basis of high streamhabitat integrity and high biological integrity. The reference-reaches median biological integrity score was equaled or exceeded by three reaches (one on Spring Creek and two on Cypress Creek) that are on the State of Texas 303(d) list of threatened or impaired waters with respect to aquatic life. This indicates that direct measures of biological integrity could be used to supplement surrogatebased designations of biological integrity such as the State list. A statistically significant multipleregression model was developed that uses independent variables that can be obtained without fieldintensive studies to predict the biological integrity score for a reach. The deviation from the model’s predicted score with the score based on biological sampling can be used to interpret the degree of biological impairment in a reach. Data from reaches outside the group of reaches used in this study are needed to test the validity of the multipleregression model.
Tampa Bay depends on ground water for most of the water supply. Numerous wetlands and lakes in Pasco County have been impacted by the high demand for ground water. Central Pasco County, particularly the area within the Cypress Creek well field, has been greatly affected. Probable causes for the decline in surface-water levels are well-field pumpage and a decade-long drought. Efforts are underway to increase surface-water levels by developing alternative sources of water supply, thus reducing the quantity of well-field pumpage. Numerical ground-water flow simulations coupled with an optimization routine were used in a series of simulations to test the sensitivity of optimal pumpage to desired increases in surficial aquifer system heads in the Cypress Creek well field. The ground-water system was simulated using the central northern Tampa Bay ground-water flow model. Pumping solutions for 1987 equilibrium conditions and for a transient 6-month timeframe were determined for five test cases, each reflecting a range of desired target recovery heads at different head control sites in the surficial aquifer system. Results are presented in the form of curves relating average head recovery to total optimal pumpage. Pumping solutions are sensitive to the location of head control sites formulated in the optimization problem and as expected, total optimal pumpage decreased when desired target head increased. The distribution of optimal pumpage for individual production wells also was significantly affected by the location of head control sites. A pumping advantage was gained for test-case formulations where hydraulic heads were maximized in cells near the production wells, in cells within the steady-state pumping center cone of depression, and in cells within the area of the well field where confining-unit leakance is the highest. More water was pumped and the ratio of head recovery per unit decrease in optimal pumpage was more than double for test cases where hydraulic heads are maximized in cells located at or near the production wells. Additionally, the ratio of head recovery per unit decrease in pumpage was about three times more for the area where confining-unit leakance is the highest than for other leakance zone areas of the well field. For many head control sites, optimal heads corresponding to optimal pumpage deviated from the desired target recovery heads. Overall, pumping solutions were constrained by the limiting recovery values, initial head conditions, and by upper boundary conditions of the ground-water flow model.
The Prairie Island Indian Community in east-central Minnesota uses ground water from the Mt. Simon-Hinckley aquifer as its source of water supply. Tribal officials implemented a Source Water Protection Program to protect the quality of this water. Areas of contributing recharge were delineated for two community well fields. At well field A are two wells 325 m apart, and at well field B are two wells 25 m apart. A steady state single layer, two-dimensional ground-water flow model constructed with the computer program MODFLOW,combined with the particle-tracking computer program MODPATH, was used to track water particles (upgradient) from the two well fields. A withdrawal rate of 625 m 3 /d was simulated for each well field. The ground-water flow paths delineated areas of contributing recharge that are 0.38 and 0.65 km 2 based on 10- and 50-year travel times, respectively. The flow paths that define these areas extend for maximum distances of about 350 and 450 m, respectively, from the wells. At well field A the area of contributing recharge was delineated for each well as separate withdrawal points. At well field B the area of contributing recharge was delineated for the two wells as a single withdrawal point. Delineation of areas of contributing recharge to the well fields from land surface would require construction of a multi-layer ground-water flow model.