Search USGS⌕ Search

SEARCH · Search USGS

Results for “Computers and Geology”

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 883 records · Page 49Linked to original sources

Statistical frequency analysis of flood records

The U.S. Geological Survey, like other Federal agencies, uses Hydrology Subcommittee Bulletin 17 for guidance in statistical frequency analysis of flood records. This paper describes the formal statistical and computational aspects of the Bulletin 17 methodology. The methodology includes provisions for dealing with high and low out-liers, historic peaks, and other anomalous flood data. If these options are inadequate, alternative procedures may be used if properly documented.

Conference Paper↗

Flood-inundation maps for the Green River in Colrain, Leyden, and Greenfield, Massachusetts, from U.S. Geological Survey streamgage 01170100 Green River near Colrain to the confluence with the Deerfield River

The U.S. Geological Survey developed flood elevations in cooperation with the Federal Emergency Management Agency for a 14.3-mile reach of the Green River in Colrain, Leyden, and Greenfield, Massachusetts, to assist landowners and emergency management workers to prepare for and recover from floods. The river reach extends from the U.S. Geological Survey Green River near Colrain, MA (01170100) streamgage downstream to the confluence with the Deerfield River. A series of seven digital flood inundation maps were developed for the upper 4.4 miles of the river reach downstream from the stream. Flood discharges corresponding to the 50-, 10-, 1-, and 0.2-percent annual exceedance probabilities were computed for the reach from updated flood-frequency analyses. These peak flows and the flood flows associated with the stages of 10.2, 12.4, and 14.4 feet (ft) at the Green River streamgage were routed through a one-dimensional step-backwater hydraulic model to obtain the corresponding peak water-surface elevations and to place the Tropical Storm Irene flood of August 28, 2011 (stage 13.97 ft), into historical context. The hydraulic model was calibrated by using the current (2015) stage-discharge relation at the U.S. Geological Survey Green River near Colrain, MA (01170100) streamgage and from documented high-water marks from the Tropical Storm Irene flood, which had a flow higher than a 0.2-percent annual exceedance probability flood discharge. The hydraulic model was used to compute water-surface profiles for flood stages referenced to the streamgage and ranging from the 50-percent annual exceedance probability (bankfull flow) at 7.6 ft (439.8 ft above the North American Vertical Datum of 1988 [NAVD 88]) to 14.4 ft (446.7 ft NAVD 88), which exceeds the maximum recorded water level of 13.97 ft (Tropical Storm Irene) at the streamgage. The mapped stages of 7.6 to 14.4 ft were selected to match the stages for bankfull; the 50-, 10-, 1-, and 0.2-percent annual exceedance probabilities; incremental stages of 10.2 and 12.4 ft; and the maximum stage of the stage-discharge rating curve. The simulated water-surface profiles were combined with a geographic information system digital elevation model derived from light detection and ranging (lidar) data having a 0.5-ft vertical accuracy to create a set of flood-inundation maps. The availability of the flood-inundation maps, combined with information regarding near real-time stage from U.S. Geological Survey Green River near Colrain, MA (01170100) streamgage, can provide emergency management personnel and residents with information that is critical for flood response activities, such as evacuations and road closures, and postflood recovery efforts. The flood-inundation maps are nonregulatory but provide Federal, State, and local agencies and the public with estimates of the potential extent of flooding during selected peak-flow events.

Massachusetts↗

Generalized altitude and configuration of the base of the High Plains regional aquifer, northwestern Oklahoma

During 1978, the U.S. Geological Survey began a 5-year study of the High Plains regional aquifer system to provide hydrologic information for evaluation of the effects of long-term development of the aquifer and to develop computer models for prediction of aquifer response to alternative changes in ground-water management (Weeks, 1978). This report is one of a series presenting hydrologic information of the High Plains aquifer in Oklahoma. This report consists of two maps showing the altitude and configuration of the eroded pre-Ogallala bedrock surface (sheets 1 and 2) that forms the base of the aquifer. Bedrock slopes generally from west to east and is composed of rocks of Permian, Triassic-Jurassic, and Cretaceous age. The subcrop extent of these rocks is shown by Morton (1973) for Oklahoma and by Weeks and Gutentag (1981) for the entire High Plains study area. Altitudes of the aquifer base were determined from drillers' logs provided by the Oklahoma Water Resources Board and from published and unpublished information in the files of the U.S. Geological Survey. Where data were sparse, the total well depth was used as the aquifer base under the assumption that wells generally are not drilled very deep into non-water yielding bedrock.

Open-File Report↗

Geologic interpretation of aeromagnetic map of the Seward and Blying Sound quadrangles, Alaska

An aeromagnetic survey of the Seward and northern part of the Blying Sound quadrangles was flown by Geometrics, Inc., during 1975-77 to provide magnetic data to aid in an appraisal of the mineral resources. Background details of the regional geology and major structures are described in a companion report by Tysdal and Case (1979), and interpretation of the gravity data are described in a report by Case and others (1979). Preliminary reports of a paleomagnetic investigation have been reported by Hillhouse and Gromme (1977), and a more complete analysis of the results is in preparation. The aeromagnetic survey was flown at a nominal height of 300 m above the surface and flight lines were spaced about 1.7 km apart. Total magnetic field was measured by a proton-precession magnetometer. The International Geomagnetic Reference Field (IGRF), updated to 1976, was removed, and the residual values were computer-contoured to produce the aeromagnetic map (sheet 1). Detailed aeromagnetic maps at scale 1:63,360 and a map at scale 1:250,000 have been released as Open File Reports 78-1080 through 1083 (U.S. Geological Survey, 1978). Magnetic anomalies in the Seward-Blying Sound region are somewhat less complex than in other parts of south-central Alaska, primarily because of the huge volume of relatively nonmagnetic flysch that constitutes the Valdez and Orca Groups. Several rather isolated groups of mafic-ultramafic rocks are significantly more magnetic than the flysch assemblages, and cause distinctive patterns of magnetic anomalies. A few of the granitic plutons are magnetic, but most appear to have little or no magnetic expression.

Alaska↗

A spatially referenced regression model (SPARROW) for suspended sediment in streams of the Conterminous U.S.

Suspended sediment has long been recognized as an important contaminant affecting water resources. Besides its direct role in determining water clarity, bridge scour and reservoir storage, sediment serves as a vehicle for the transport of many binding contaminants, including nutrients, trace metals, semi-volatile organic compounds, a nd numerous pesticides (U.S. Environmental Protection Agency, 2000a). Recent efforts to addr ess water-quality concerns through the Total Maximum Daily Load (TMDL) process have iden tified sediment as the single most prevalent cause of impairment in the Nation’s streams a nd rivers (U.S. Environmental Protection Agency, 2000b). Moreover, sediment has been identified as a medium for the tran sport and sequestration of organic carbon, playing a potentia lly important role in understa nding sources and sinks in the global carbon budget (Stallard, 1998). A comprehensive understanding of sediment fate a nd transport is considered essential to the design and implementation of effective plans for sediment management (Osterkamp and others, 1998, U.S. General Accounting Office, 1990). An exte nsive literature addr essing the problem of quantifying sediment transport has produced a nu mber of methods for estimating its flux (see Cohn, 1995, and Robertson and Roerish, 1999, for us eful surveys). The accuracy of these methods is compromised by uncertainty in the concentration measurements and by the highly episodic nature of sediment movement, particul arly when the methods are applied to smaller basins. However, for annual or decadal flux es timates, the methods are generally reliable if calibrated with extended periods of data (Robertson and Roerish, 1999). A substantial literature also supports the Universal Soil Loss Equation (U SLE) (Soil Conservation Service, 1983), an engineering method for estimating sheet and rill erosion, although the empirical credentials of the USLE have recently been questioned (Tri mble and Crosson, 2000). Conversely, relatively little direct evidence is available concerning the fate of sediment. The common practice of quantifying sediment fate with a sediment deliv ery ratio, estimated from a simple empirical relation with upstream basin area, does not artic ulate the relative importance of individual storage sites within a basin (Wolman, 1977). Rates of sediment deposition in reservoirs and flood plains can be determined from empirical measurement s , but only a limited number of sites have been monitored, and net rates of deposition or loss from other potential sinks and sources is largely unknown (Stallard, 1998). In particular, little is known about how much sediment loss from fields ultimately makes its way to stream channels, and how much sediment is subsequently stored in or lost from th e streambed (Meade and Parker, 1985, Trimble and Crosson, 2000). This paper reports on recent progress made to a ddress empirically the question of sediment fate and transport on a national scale. The model pres ented here is based on the SPAtially Referenced Regression On Watershed attr ibutes (SPARROW) methodology, fi rst used to estimate the distribution of nutrients in str eams and rivers of the United Stat es, and subsequently shown to describe land and stream processes affecting the delivery of nutrients (Smith and others, 1997, Alexander and others, 2000, Preston and Brakeb ill, 1999). The model makes use of numerous spatial datasets, available at the national level, to explain long-term sediment water-quality conditions in major streams and rivers throughou t the United States. Sediment sources are identified using sediment erosion rates from the National Resources I nventory (NRI) (Natural Resources Conservation Service, 2000) and apportioned over the landscape according to 30- meter resolution land-use information from th e National Land Cover Data set (NLCD) (U.S. Geological Survey, 2000a). More than 76,000 reservoirs from the National Inventory of Dams (NID) (U.S. Army Corps of Engin eers, 1996) are identified as pot ential sediment sinks. Other, non-anthropogenic sources and sinks are identified using soil in formation from the State Soil Survey Geographic (STATSGO) data base (Schwarz and Alexander, 1995) and spatial coverages representing surficial rock t ype and vegetative cover. The SPA RROW model empirically relates these diverse spatial datasets to estimates of long-term, mean annual sediment flux computed from concentration and flow measurements co llected over the period 1985 -95 from more than 400 monitoring stations maintained by the Na tional Stream Quality Accounting Network (Alexander and others, 1998), the National Wa ter Quality Assessment Program, and U.S. Geological Survey District offices (Turcios and Gray, in press). Th e calibrated model is used to estimate sediment flux for over 60,000 stream segments included in the River Reach File 1 (RF1) stream network (Alexander and others, 1999). SPARROW uses statis tical methods to calibrate a simple, structural model of riverine water quality, one that imposes mass ba lance in accounting for changes in contaminant flux. As applied here, the mass-balance approach facilitates the interpretation of model results in terms of physical processes affecting sediment transport, and makes possible the estimation of various rates of sediment generation and loss associated with stream channels and features of the landscape. The statistical approach provides a basi s for assessing the error of these inferred rates and of the error in extrapolated estimates of sediment flux made for streams in the RF1 network. An important implication of the holistic modeling approach adopted in this analysis is that estimates of sediment production and loss ar e based on, and therefore consistent with, measurements of in-stream flux. Other ancillary information, such as direct measurements of long-term sediment storage and release from rese rvoirs (Steffen, 1996), is incorporated into the analysis by specifying additional equations expl aining these ancillary variables. The imposition of cross-equation constraints affords this info rmation a statistically consistent weight in explaining in-stream sediment flux. Thus, the me thodology described here represents a general framework for synthesizing a wide spectrum of available information relevant to the understanding of sediment fate and transport.

Conterminous United States↗

Precise determination of δ 88 Sr in rocks, minerals, and waters by double-spike TIMS: A powerful tool in the study of chemical, geologic, hydrologic and biologic processes

We present strontium isotopic ( 88 Sr/ 86 Sr and 87 Sr/ 86 Sr) results obtained by 87 Sr– 84 Sr double spike thermal ionization mass-spectrometry (DS-TIMS) for several standards as well as natural water samples and mineral samples of abiogenic and biogenic origin. The detailed data reduction algorithm and a user-friendly Sr-specific stand-alone computer program used for the spike calibration and the data reduction are also presented. Accuracy and precision of our δ 88 Sr measurements, calculated as permil (‰) deviations from the NIST SRM-987 standard, were evaluated by analyzing the NASS-6 seawater standard, which yielded δ 88 Sr = 0.378 ± 0.009‰. The first DS-TIMS data for the NIST SRM-607 potassium feldspar standard and for several US Geological Survey carbonate, phosphate, and silicate standards (EN-1, MAPS-4, MAPS-5, G-3, BCR-2, and BHVO-2) are also reported. Data obtained during this work for Sr-bearing solids and natural waters show a range of δ 88 Sr values of about 2.4‰, the widest observed so far in terrestrial materials. This range is easily resolvable analytically because the demonstrated external error (±SD, standard deviation) for measured δ 88 Sr values is typically ≤0.02‰. It is shown that the “true” 87 Sr/ 86 Sr value obtained by the DS-TIMS or any other external normalization method combines radiogenic and mass-dependent mass-fractionation effects, which cannot be separated. Therefore, the “true” 87 Sr/ 86 Sr and the δ 87 Sr parameter derived from it are not useful isotope tracers. Data presented in this paper for a wide range of naturally occurring sample types demonstrate the potential of the δ 88 Sr isotope tracer in combination with the traditional radiogenic 87 Sr/ 86 Sr tracer for studying a variety of biological, hydrological, and geological processes.

Journal of Analytical Atomic Spectrometry↗

Estimation and comparison of 1-percent annual exceedance probability flood flows at Federal Emergency Management Agency flood insurance study flow locations across Pennsylvania

Flood-flow estimates were computed at over 5,000 Federal Emergency Management Agency (FEMA) flood insurance study (FIS) flow locations across Pennsylvania for the 1-percent annual exceedance probability flood event (1-percent AEP). Depending on a point of interest’s proximity to a streamgage, weighting techniques may be applied to obtain flood-flow estimates for ungaged flow locations using observed peak-flow data from a nearby streamgage. Following the U.S. Geological Survey’s (USGS) published guidance, stream segments were identified where the drainage-area ratio method could be leveraged. Using updated regional regression equations and recently published flood-flow estimates at USGS streamgage locations following USGS Bulletin 17C guidelines, weighted and transferred flood flows were computed, where appropriate. For locations not applicable for the drainage-area ratio method, regression equations were used to compute flood-flow estimates. These flood-flow estimates were then compared to FEMA FIS 1-percent AEP flood-flow estimates. Percentage-difference values were computed for 3,599 FIS flow locations determined to be suitable for analysis, finding that USGS-derived flood-flow estimates were consistently lower than FEMA FIS flood-flow estimates with a statewide median percentage difference of −10.1 percent. The dataset was normally distributed with a standard deviation of 45.7 percent. Allegheny County was found to have 74 FIS flow locations with percentage-difference values greater than or equal to 67 percent or less than or equal to −67 percent. The flood-flow region in which Allegheny County is contained, Region 2, had a median percentage-difference value of −39 percent. Although removed from the final analysis, flow locations with drainage-area values above the recommended threshold for regression-based estimation (about 1,000 square miles [mi 2 ]) were observed to have consistently higher percentage-difference values; a reminder of the limitations of use for regression-based flood-flow estimates. This report, the comparisons within, and a companion data release are intended to serve as tools to FEMA in assisting with the ongoing assessment of FIS flow locations across Pennsylvania.

Pennsylvania↗

National assessment of shoreline change—Summary statistics for updated vector shorelines and associated shoreline change data for the north coast of Alaska, U.S.-Canadian Border to Icy Cape

Long-term rates of shoreline change for the north coast of Alaska, from the U.S.-Canadian border to the Icy Cape region of northern Alaska, have been updated as part of the U.S. Geological Survey’s National Assessment of Shoreline Change Project. Short-term shoreline change rates are reported for the first time. Additional shoreline position data were used to compute rates where the previous rate-of-change assessment only included two shoreline positions at a given location. The calculation of uncertainty associated with the long-term average rates has also been updated to match refined methods used in other study regions of the National Assessment of Shoreline Change Project. The average rates of this report have a reduced amount of uncertainty compared to those presented in the first assessment for this region.

Alaska↗

A homogeneous stochastic model for earthquake occurrences

The objective of this study is to develop a probabilistic model for earthquake occurrences with temporal and spatial memory. Stochastic processes are used to characterize both the spatial and temporal dependencies of seismic occurrences along a fault. Currently, only homogeneous space and time transitions are considered. The resulting process however is evolutionary and depends on the specific sequence of events over a period of time. The model provides estimates on the cumulative activity of a fault over a future time period. In addition, probabilities of occurrences of individual events along a geologic fault at some specified future time are obtained when a descritized time scale is used. The information from these evaluations is particularly useful in engineering seismic hazard computations and for social and engineering seismic risk assessments.

Open-File Report↗

Compilation of ground water quality data in Pennsylvania

The U.S. Geological Survey's water quality file of 4,671 wells and springs in Pennsylvania provided ground-water-quality data for Pennsylvania. The data were assembled into computer-readable format and sorted into 15 major aquifer groups based on principal lithology, physiographic province, and age. Nineteen variables in each group were summarized by the Statistical Analysis Systems UNIVARIATE procedures to produce descriptive statistics including extreme values and quartiles. The bulk of the water-quality data are in the important aquifers in the unconsolidated Coastal Plain sediments, the Triassic sedimentary rocks, the igneous and metamorphic rocks, and the carbonate rocks. On the other hand, water-quality data for aquifer groups in the Appalachian Plateau and Valley and Ridge Province are sparse. Statewide, only six wells provide sufficient long-term data for trend analyses. Ground-water quality in Pennsylvania is highly diverse. High concentrations of dissolved solids, iron, manganese, sulfate, and nitrate are prevalent forms of natural and manmade contamination. The unconsolidated Coastal Plain aquifers have been most severely degraded. On the other hand, some of the best quality water is found in the quartzite, sandstone, and conglomerate rock units in the Cambrian and Precambrian rocks.

Open-File Report↗

Characterization of Suspended-Sediment Loading to and from John Redmond Reservoir, East-Central Kansas, 2007-2008

Storage capacity in John Redmond Reservoir is being lost to sedimentation more rapidly than in other federal impoundments in Kansas. The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, initiated a study to characterize suspended-sediment loading to and from John Redmond Reservoir from February 21, 2007, through February 21, 2008. Turbidity sensors were installed at two U.S. Geological Survey stream gages upstream (Neosho River near Americus and the Cottonwood River near Plymouth) and one stream gage downstream (Neosho River at Burlington) from the reservoir to compute continuous, real-time (15-minute) measurements of suspended-sediment concentration and loading. About 1,120,000 tons of suspended-sediment were transported to, and 100,700 tons were transported from John Redmond Reservoir during the study period. Dependent on the bulk density of sediment stored in the reservoir, 5.0 to 1.4 percent of the storage in the John Redmond conservation pool was lost during the study period, with an average deposition of 3.4 to 1.0 inches. Nearly all (98-99 percent) of the incoming sediment load was transported during 9 storms which occurred 25 to 27 percent of the time. The largest storm during the study period (peak-flow recurrence interval of about 4.6-4.9 years) transported about 37 percent of the sediment load to the reservoir. Suspended-sediment yield from the unregulated drainage area upstream from the Neosho River near Americus was 530 tons per square mile, compared to 400 tons per square mile upstream from the Cottonwood River near Plymouth. Comparison of historical (1964-78) to current (2007) sediment loading estimates indicate statistically insignificant (<90 percent confidence) differences at the Neosho River near Americus and the Cottonwood River near Plymouth, but a significant (>99 percent) decrease in sediment loading at the Neosho River at Burlington. Ninety-percent confidence intervals of streamflow-derived estimates of total sediment load were 7 to 21 times larger than turbidity-derived estimates. Results from this study can be used by natural resource managers to calibrate sediment models and estimate the ability of John Redmond Reservoir to support designated uses into the future.

Scientific Investigations Report↗

Estimating selected low-flow frequency statistics and harmonic-mean flows for ungaged, unregulated streams in Indiana

Information on low-flow characteristics of streams is essential for the management of water resources. This report provides equations for estimating the 1-, 7-, and 30-day mean low flows for a recurrence interval of 10 years and the harmonic-mean flow at ungaged, unregulated stream sites in Indiana. These equations were developed using the low-flow statistics and basin characteristics for 108 continuous-record streamgages in Indiana with at least 10 years of daily mean streamflow data through the 2011 climate year (April 1 through March 31). The equations were developed in cooperation with the Indiana Department of Environmental Management. Regression techniques were used to develop the equations for estimating low-flow frequency statistics and the harmonic-mean flows on the basis of drainage-basin characteristics. A geographic information system was used to measure basin characteristics for selected streamgages. A final set of 25 basin characteristics measured at all the streamgages were evaluated to choose the best predictors of the low-flow statistics. Logistic-regression equations applicable statewide are presented for estimating the probability that selected low-flow frequency statistics equal zero. These equations use the explanatory variables total drainage area, average transmissivity of the full thickness of the unconsolidated deposits within 1,000 feet of the stream network, and latitude of the basin outlet. The percentage of the streamgage low-flow statistics correctly classified as zero or nonzero using the logistic-regression equations ranged from 86.1 to 88.9 percent. Generalized-least-squares regression equations applicable statewide for estimating nonzero low-flow frequency statistics use total drainage area, the average hydraulic conductivity of the top 70 feet of unconsolidated deposits, the slope of the basin, and the index of permeability and thickness of the Quaternary surficial sediments as explanatory variables. The average standard error of prediction of these regression equations ranges from 55.7 to 61.5 percent. Regional weighted-least-squares regression equations were developed for estimating the harmonic-mean flows by dividing the State into three low-flow regions. The Northern region uses total drainage area and the average transmissivity of the entire thickness of unconsolidated deposits as explanatory variables. The Central region uses total drainage area, the average hydraulic conductivity of the entire thickness of unconsolidated deposits, and the index of permeability and thickness of the Quaternary surficial sediments. The Southern region uses total drainage area and the percent of the basin covered by forest. The average standard error of prediction for these equations ranges from 39.3 to 66.7 percent. The regional regression equations are applicable only to stream sites with low flows unaffected by regulation and to stream sites with drainage basin characteristic values within specified limits. Caution is advised when applying the equations for basins with characteristics near the applicable limits and for basins with karst drainage features and for urbanized basins. Extrapolations near and beyond the applicable basin characteristic limits will have unknown errors that may be large. Equations are presented for use in estimating the 90-percent prediction interval of the low-flow statistics estimated by use of the regression equations at a given stream site. The regression equations are to be incorporated into the U.S. Geological Survey StreamStats Web-based application for Indiana. StreamStats allows users to select a stream site on a map and automatically measure the needed basin characteristics and compute the estimated low-flow statistics and associated prediction intervals.

Indiana↗

Ground-water recharge in the arid and semiarid southwestern United States

Ground-water recharge in the arid and semiarid southwestern United States results from the complex interplay of climate, geology, and vegetation across widely ranging spatial and temporal scales. Present-day recharge tends to be narrowly focused in time and space. Widespread water-table declines accompanied agricultural development during the twentieth century, demonstrating that sustainable ground-water supplies are not guaranteed when part of the extracted resource represents paleorecharge. Climatic controls on ground-water recharge range from seasonal cycles of summer monsoonal and winter frontal storms to multimillennial cycles of glacial and interglacial periods. Precipitation patterns reflect global-scale interactions among the oceans, atmosphere, and continents. Large-scale climatic influences associated with El Niño and Pacific Decadal Oscillations strongly, but irregularly, control weather in the study area, so that year-to-year variations in precipitation and ground-water recharge are large and difficult to predict. Proxy data indicate geologically recent periods of naturally occurring multidecadal droughts unlike any in the modern instrumental record. Any anthropogenically induced climate change will likely reduce ground-water recharge through diminished snowpack at higher elevations. Future changes in El Niño and monsoonal patterns, both crucial to precipitation in the study area, are highly uncertain in current models. Current land-use modifications influence ground-water recharge through vegetation, irrigation, and impermeable area. High mountain ranges bounding the study area—the San Bernadino Mountains and Sierra Nevada to the west, and the Wasatch and southern Colorado Rocky Mountains to the east—provide external geologic controls on ground-water recharge. Internal geologic controls stem from tectonic processes that led to numerous, variably connected alluvial-filled basins, exposure of extensive Paleozoic aquifers in mountainous recharge areas, and distinct modes of recharge in the Colorado Plateau and Basin and Range subregions. The chapters in this professional paper present (first) an overview of climatic and hydrogeologic framework (chapter A), followed by a regional analysis of ground-water recharge across the entire study area (chapter B). These are followed by an overview of site-specific case studies representing different subareas of the geographically diverse arid and semiarid southwestern United States (chapter C); the case studies themselves follow in chapters D–K. The regional analysis includes detailed hydrologic modeling within the framework of a high-resolution geographic-information system (GIS). Results from the regional analysis are used to explore both the distribution of ground-water recharge for mean climatic conditions as well as the influence of two climatic patterns—the El Niño-Southern Oscillation and Pacific Decadal Oscillation—that impart a high degree of variability to the hydrologic cycle. Individual case studies employ a variety of geophysical and geochemical techniques to investigate recharge processes and relate the processes to local geologic and climatic conditions. All of the case studies made use of naturally occurring tracers to quantify recharge. Thermal and geophysical techniques that were developed in the course of the studies are presented in appendices. The quantification of ground-water recharge in arid settings is inherently difficult due to the generally low amount of recharge, its spatially and temporally spotty nature, and the absence of techniques for directly measuring fluxes entering the saturated zone from the unsaturated zone. Deep water tables in arid alluvial basins correspond to thick unsaturated zones that produce up to millennial time lags between changes in hydrologic conditions at the land surface and subsequent changes in recharge to underlying ground water. Recent advances in physical, chemical, isotopic, and modeling techniques have fostered new types of recharge assessments. Chemical and isotopic techniques include an increasing variety of environmental tracers that are useful and robust. Physically based techniques include the use of heat as a tracer and computationally intensive geophysical imaging tools for characterizing hydrologic conditions in the unsaturated zone. Modeling-based techniques include spatially distributed water-budget computations using high-resolution remotely sensed and ground-based geographic data. Application of these techniques to arid and semiarid settings in the southwestern United States reveals distinct patterns of recharge corresponding to geologic setting, climatic and vegetative history, and land use. Analysis of recharge patterns shows that large expanses of alluvial basin floors are drying out under current climatic conditions, with little to no recharge to underlying ground water. Ground-water recharge occurs mainly beneath upland catchments in which thin soils overlie permeable bedrock, ephemeral channels in which flow may average only several hours per year, and active agricultural areas. The chapters in this professional paper represent a coordinated attempt to develop a better understanding of one of the Nation's most critical yet difficult-to-quantify renewable resources.

Professional Paper↗

Flood-inundation maps for Indian Creek and Tomahawk Creek, Johnson County, Kansas, 2014

Digital flood-inundation maps for a 6.4-mile upper reach of Indian Creek from College Boulevard to the confluence with Tomahawk Creek, a 3.9-mile reach of Tomahawk Creek from 127th Street to the confluence with Indian Creek, and a 1.9-mile lower reach of Indian Creek from the confluence with Tomahawk Creek to just beyond the Kansas/Missouri border at State Line Road in Johnson County, Kansas, were created by the U.S. Geological Survey in cooperation with the city of Overland Park, Kansas. The flood-inundation maps, which can be accessed through the U.S. Geological Survey Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at the U.S. Geological Survey streamgages on Indian Creek at Overland Park, Kansas; Indian Creek at State Line Road, Leawood, Kansas; and Tomahawk Creek near Overland Park, Kansas. Near real time stages at these streamgages may be obtained on the Web from the U.S. Geological Survey National Water Information System at http://waterdata.usgs.gov/nwis or the National Weather Service Advanced Hydrologic Prediction Service at http://water.weather.gov/ahps/ , which also forecasts flood hydrographs at these sites. Flood profiles were computed for the stream reaches by means of a one-dimensional step-backwater model. The model was calibrated for each reach by using the most current stage-discharge relations at the streamgages. The hydraulic models were then used to determine 15 water-surface profiles for Indian Creek at Overland Park, Kansas; 17 water-surface profiles for Indian Creek at State Line Road, Leawood, Kansas; and 14 water-surface profiles for Tomahawk Creek near Overland Park, Kansas, for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from bankfull to the next interval above the 0.2-percent annual exceedance probability flood level (500-year recurrence interval). The simulated water-surface profiles were then combined in a geographic information system with a digital elevation model derived from light detection and ranging data (having a 0.429-foot vertical and 0.228-foot horizontal accuracy) to delineate the area flooded at each water level. The availability of these maps, along with Web information regarding current stage from the U.S. Geological Survey streamgages and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations, road closures, and postflood recovery efforts.

Kansas↗

Geotherm: the U.S. geological survey geothermal information system

GEOTHERM is a comprehensive system of public databases and software used to store, locate, and evaluate information on the geology, geochemistry, and hydrology of geothermal systems. Three main databases address the general characteristics of geothermal wells and fields, and the chemical properties of geothermal fluids; the last database is currently the most active. System tasks are divided into four areas: (1) data acquisition and entry, involving data entry via word processors and magnetic tape; (2) quality assurance, including the criteria and standards handbook and front-end data-screening programs; (3) operation, involving database backups and information extraction; and (4) user assistance, preparation of such items as application programs, and a quarterly newsletter. The principal task of GEOTHERM is to provide information and research support for the conduct of national geothermal-resource assessments. The principal users of GEOTHERM are those involved with the Geothermal Research Program of the U.S. Geological Survey. Information in the system is available to the public on request. ?? 1983.

Computers & Geosciences↗

Rapid estimation of minimum depth-to-bedrock from lidar leveraging deep-learning-derived surficial material maps

Previously glaciated landscapes often share similar surficial characteristics, including large areas of exposed bedrock, blankets of till deposits, and alluvium-floored valleys. These materials play significant roles in geologic and hydrologic resources, geohazards, and landscape evolution; however, the vast extents of many previously glaciated landscapes have rendered comprehensive, detailed field mapping difficult. While recent advances in remote sensing have facilitated mapping of surficial materials and landforms, manual map creation has remained a time-intensive task. The development of convolutional neural networks (CNNs) for image classification has provided a new opportunity for rapid characterization of digital elevation models, thus enabling efficient mapping of surficial materials and landforms. We have developed a methodology that leverages existing geologic maps and high-resolution (1–3 m) lidar data to train a U-Net CNN to classify alluvium and exposed bedrock in previously glaciated regions. Coupled with U.S. Geological Survey-developed geomorphometry tools capable of approximating stream incision depths, these classifications can be used to estimate the minimum thicknesses of stream-proximal hillslope sediments in areas where streams have undergone minimal incision into bedrock. We validate this approach in the context of the Neversink River watershed, a subbasin of the Delaware River Basin and significant water source for New York City. Evaluation of deep learning model performance demonstrates substantial agreement with manually drawn maps of alluvium and exposed bedrock. Validation of the minimum sediment thickness map using borehole data and passive seismic measurements shows the greatest performance for shallow materials and decreased performance in deep sediments, as well as in areas where bedrock exposures were too small to be resolved by lidar. To resolve these issues and create more accurate surficial maps, we are training new CNNs with additional geologic data and exploring advanced approaches for estimating depths of stream incision.

Delaware, New Jersey, New York, Pennsylvania↗

Application of the source-area concept of storm runoff to a small Arizona watershed

An attempt to demonstrate the source-area concept of storm runoff by analysis of the rainfall-runoff relation for the watershed of Queen Creek tributary in south-central Arizona was moderately successful. The demonstration was somewhat marred by the necessity to make several simplifying assumptions to eliminate some of the many basin variables of unknown magnitude. The percentage of watershed contributing storm runoff was related to total rainfall received that is, antecedent and storm rainfall on the assumption that the availability of rainfall excess from any subarea was dependent on the saturation, or near saturation, of a permeable upper layer of soil. An average unit hydrograph was used for all computed sizes of the runoff-contributing area.

Arizona↗

MODFLOW-2000, the U.S. Geological Survey modular ground-water model; user guide to the observation, sensitivity, and parameter-estimation processes and three post-processing programs

This report documents the Observation, Sensitivity, and Parameter-Estimation Processes of the ground-water modeling computer program MODFLOW-2000. The Observation Process generates model-calculated values for comparison with measured, or observed, quantities. A variety of statistics is calculated to quantify this comparison, including a weighted least-squares objective function. In addition, a number of files are produced that can be used to compare the values graphically. The Sensitivity Process calculates the sensitivity of hydraulic heads throughout the model with respect to specified parameters using the accurate sensitivity-equation method. These are called grid sensitivities. If the Observation Process is active, it uses the grid sensitivities to calculate sensitivities for the simulated values associated with the observations. These are called observation sensitivities. Observation sensitivities are used to calculate a number of statistics that can be used (1) to diagnose inadequate data, (2) to identify parameters that probably cannot be estimated by regression using the available observations, and (3) to evaluate the utility of proposed new data. The Parameter-Estimation Process uses a modified Gauss-Newton method to adjust values of user-selected input parameters in an iterative procedure to minimize the value of the weighted least-squares objective function. Statistics produced by the Parameter-Estimation Process can be used to evaluate estimated parameter values; statistics produced by the Observation Process and post-processing program RESAN-2000 can be used to evaluate how accurately the model represents the actual processes; statistics produced by post-processing program YCINT-2000 can be used to quantify the uncertainty of model simulated values. Parameters are defined in the Ground-Water Flow Process input files and can be used to calculate most model inputs, such as: for explicitly defined model layers, horizontal hydraulic conductivity, horizontal anisotropy, vertical hydraulic conductivity or vertical anisotropy, specific storage, and specific yield; and, for implicitly represented layers, vertical hydraulic conductivity. In addition, parameters can be defined to calculate the hydraulic conductance of the River, General-Head Boundary, and Drain Packages; areal recharge rates of the Recharge Package; maximum evapotranspiration of the Evapotranspiration Package; pumpage or the rate of flow at defined-flux boundaries of the Well Package; and the hydraulic head at constant-head boundaries. The spatial variation of model inputs produced using defined parameters is very flexible, including interpolated distributions that require the summation of contributions from different parameters. Observations can include measured hydraulic heads or temporal changes in hydraulic heads, measured gains and losses along head-dependent boundaries (such as streams), flows through constant-head boundaries, and advective transport through the system, which generally would be inferred from measured concentrations. MODFLOW-2000 is intended for use on any computer operating system. The program consists of algorithms programmed in Fortran 90, which efficiently performs numerical calculations and is fully compatible with the newer Fortran 95. The code is easily modified to be compatible with FORTRAN 77. Coordination for multiple processors is accommodated using Message Passing Interface (MPI) commands. The program is designed in a modular fashion that is intended to support inclusion of new capabilities.

Open-File Report↗