USGS ScienceSearch

SEARCH · USGS Science

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 253 records · Page 14Linked to original sources

Computer-model analysis of ground-water flow and simulated effects of contaminant remediation at Naval Weapons Industrial Reserve Plant, Dallas, Texas

In June 1993, the Department of the Navy, Southern Division Naval Facilities Engineering Command (SOUTHDIV), began a Resource Conservation and Recovery Act (RCRA) Facility Investigation (RFI) of the Naval Weapons Industrial Reserve Plant (NWIRP) in north-central Texas. The RFI has found trichloroethene, dichloroethene, vinyl chloride, as well as chromium, lead, and other metallic residuum in the shallow alluvial aquifer underlying NWIRP. These findings and the possibility of on-site or off-site migration of contaminants prompted the need for a ground-water-flow model of the NWIRP area. The resulting U.S. Geological Survey (USGS) model: (1) defines aquifer properties, (2) computes water budgets, (3) delineates major flowpaths, and (4) simulates hydrologic effects of remediation activity. In addition to assisting with particle-tracking analyses, the calibrated model could support solute-transport modeling as well as help evaluate the effects of potential corrective action. The USGS model simulates steadystate and transient conditions of ground-water flow within a single model layer. The alluvial aquifer is within fluvial terrace deposits of Pleistocene age, which unconformably overlie the relatively impermeable Eagle Ford Shale of Late Cretaceous age. Over small distances and short periods, finer grained parts of the aquifer are separated hydraulically; however, most of the aquifer is connected circuitously through randomly distributed coarser grained sediments. The top of the underlying Eagle Ford Shale, a regional confining unit, is assumed to be the effective lower limit of ground-water circulation and chemical contamination. The calibrated steady-state model reproduces long-term average water levels within +5.1 or –3.5 feet of those observed; the standard error of the estimate is 1.07 feet with a mean residual of 0.02 foot. Hydraulic conductivity values range from 0.75 to 7.5 feet per day, and average about 4 feet per day. Specific yield values range from 0.005 to 0.15 and average about 0.08. Simulated infiltration rates range from 0 to 2.5 inches per year, depending mostly on local patterns of ground cover. Computer simulation indicates that, as of December 31, 1998, remediation systems at NWIRP were removing 7,375 cubic feet of water per day from the alluvial aquifer, with 3,050 cubic feet per day coming from aquifer storage. The resulting drawdown prevented 1,800 cubic feet per day of ground water from discharging into Cottonwood Bay, as well as inducing another 1,325 cubic feet per day into the aquifer from the bay. An additional 1,200 cubic feet of water per day (compared to pre-remediation conditions) was prevented from discharging into the west lagoon, east lagoon, Mountain Creek Lake, and Mountain Creek swale. Particle-tracking simulations, assuming an aquifer porosity of 0.15, were made to delineate flowpath patterns, or contaminant “capture zones,” resulting from 2.5- and 5-year periods of remediation activity at NWIRP. The resulting flowlines indicate three such zones, or areas from which ground water is simulated to have been removed during July 1996–December 1998, as well as extended areas from which ground water would be removed during the next 2.5 years (January 1999– June 2001). Simulation indicates that, as of December 31, 1998, the recovery trench was intercepting about 827 cubic feet per day of ground water that—without the trench—would have discharged into Cottonwood Bay. During this time, the trench is simulated to have removed about 3,221 cubic feet per day of water from the aquifer, with about 934 cubic feet per day (29 percent) coming from the south (Cottonwood Bay) side of the trench.

Texas

Characterization and occurrence of confined and unconfined aquifers in Quaternary sediments in the glaciated conterminous United States

The glacial aquifer system, which is a collection of aquifers within Quaternary sediments in the glaciated conterminous United States, is a principal aquifer that supplies groundwater that serves about 42 million people and accounts for about 5 percent of the Nation’s drinking water. This aquifer system (the area of maximum glacial advance) underlies parts of 25 States and covers 1.87×10 6 square kilometers. A hydrogeologic framework is presented that divides the glaciated United States into 17 distinct hydrogeologic terranes using a geologic approach based on previous mapping. Each hydrogeologic terrane contains Quaternary sediment that is derived from a common depositional history and can be characterized by similar texture and thickness. Characteristics of Quaternary sediments are described using attributes computed from a lithologic database of well logs compiled from 24 States (excluding Kentucky). The hydrogeologic framework presents a nationwide picture of the glacial aquifer system and provides generalizations concerning the nature of aquifers within it (for example, whether the aquifers are shallow or deep, and unconfined or confined). In this way insights can be gained from understanding the similarities and differences in distinct parts of the glacial aquifer system and how they relate to water use and quality and to aquifer vulnerability. Delineation of hydrogeologic terranes was based on an interpretation of existing geologic mapping of Quaternary sediments and the thickness of unconsolidated material. Overall thickness of Quaternary sediment was used to qualitatively rank the generalized complexity of the hydrogeologic framework in each terrane: “lower” complexity (assigned a terrane code 1), “moderate” complexity (terrane code 2), and “higher” complexity (terrane code 3). Letter designations appended to the terrane codes (for example, 1A, 1B, or 1C) differentiate terranes of similar complexity. Two unique areas, where thick, stratified, coarse-grained sediment dominates, were assigned terrane code 4. Elements of this hydrogeologic framework include a glacial environments and surficial sediments geodatabase, which includes lithologic, geomorphic, and stratigraphic characterization of Quaternary sediments based on previous mapping; a gridded database of sediment and aquifer characteristics computed from lithologic logs obtained from water-well driller records; a water-use database with information on public-water supply systems and sources of groundwater; and estimated recharge computed from a geologically based soil-water balance model. A generalized map of the bedrock geology based on previous State-level mapping is included as well. Quaternary sediment in the glaciated United States includes glacial, postglacial (Holocene) and nonglacial sediments. At land surface, 60 percent of the glacial sediment is till. Large areas of outwash and ice contact sediments are extensive throughout the Midwest but generally are confined to valleys in the Northeast and the Northwest. Lacustrine sediments were deposited in proglacial lakes adjacent to the present Great Lakes and in glacial Lake Agassiz in the eastern Dakotas and northwestern Minnesota. The median thickness of Quaternary sediment ranges from 6 to 45 meters across the 17 hydrogeologic terranes, but the maximum thickness is more than 500 meters in some areas. Quaternary sediments generally contain less than 10 percent coarse material; the median range is near zero percent under till to about 50 percent under ice contact and outwash sediments. About 80 percent of the coarse material lies within 25 to 40 meters of land surface. In most of the glaciated United States, there is a small likelihood of penetrating an aquifer-material interval containing coarse material at least 3 meters thick. A single aquifer-material interval was recorded in about 44 percent of lithologic logs, whereas about 11 percent of the logs penetrated multiple intervals. About 44 percent of water wells in the lithologic database are completed in Quaternary sediment, and many of these Quaternary water wells (42 percent) are confined by at least 7.5 meters of fine materials. About 33 percent of these Quaternary water wells are unconfined—the remainder are where only thin layers (less than 3 meters) of coarse material are present. The median depths of Quaternary water wells range from 13 to 40 meters among the 17 hydrogeologic terranes. Recharge ranges from more than 400 millimeters per year in the Northeast to 11 millimeters or less per year in the Dakotas and Montana (median value of 136 millimeters per year). Annual groundwater withdrawals compiled by county range on an areal basis from less than 1 to 370 millimeters per year, and the mean is 7.4 millimeters per year. About 36 percent of the withdrawals are for public-water supply, of which 70 percent are derived from Quaternary sediments. Groundwater withdrawals are less than 10 percent of recharge throughout most of the glaciated conterminous United States but are a larger proportion of recharge near urban areas in the Northeast and the Midwest, and in counties throughout drier parts of the Midwest. The salient characteristics of the 17 hydrogeologic terranes are presented through maps and a set of descriptive plots to facilitate visual comparisons between selected sediment and aquifer characteristics. The thickness of Quaternary sediment generally increases from the lower complexity terranes through the higher complexity terranes, consistent with their delineation. Median proportions of coarse material in Quaternary sediment and depths to aquifer-material intervals are highly variable (less than 10 to 50 percent, and 0 to 30 meters, respectively). Median thicknesses of aquifer-material intervals generally fall within a narrow range (10 to 20 meters), except in two terranes that contain thick coarse-grained sediment (30 to 35 meters). The source of water in wells varies from mostly bedrock wells in the lower complexity terranes to mostly Quaternary wells in the higher complexity terranes where the sediment is thickest. A tree diagram compiled from a hierarchical cluster analysis of a matrix composed of metrics based on sediment and aquifer characteristics, and the distribution of water wells in each terrane, indicates some groups of terranes that can be treated as comparable when analyzing groundwater flow and quality. Aquifer-material intervals indicated on maps prepared from the lithologic logs, including unconfined and confined conditions, correlate well with aquifer systems delineated on state maps for Illinois, Indiana, and North Dakota. The large scale of the study limits the resolution at which the maps can be interpreted, however, and alluvial units are not mapped correctly for some valleys in the Northeast and the Northwest. Lithologic logs used in the study are biased toward shallow depths because not all logs penetrate the entire thickness of Quaternary sediment, but this bias should not limit the utility of the sediment and aquifer descriptions because shallow depths are commonly exploited for water supply. The hydrogeologic framework will support ongoing studies of groundwater flow and quality in the U.S. Geological Survey National Water Quality Assessment program for the glaciated United States.

Scientific Investigations Report

Selected low-flow frequency statistics for continuous-record streamgage locations in Maryland, 2010

According to a 2008 report by the Governor's Advisory Committee on the Management and Protection of the State's Water Resources, Maryland's population grew by 35 percent between 1970 and 2000, and is expected to increase by an additional 27 percent between 2000 and 2030. Because domestic water demand generally increases in proportion to population growth, Maryland will be facing increased pressure on water resources over the next 20 years. Water-resources decisions should be based on sound, comprehensive, long-term data and low-flow frequency statistics from all available streamgage locations with unregulated streamflow and adequate record lengths. To provide the Maryland Department of the Environment with tools for making future water-resources decisions, the U.S. Geological Survey initiated a study in October 2009 to compute low-flow frequency statistics for selected streamgage locations in Maryland with 10 or more years of continuous streamflow records. This report presents low-flow frequency statistics for 114 continuous-record streamgage locations in Maryland. The computed statistics presented for each streamgage location include the mean 7-, 14-, and 30-consecutive day minimum daily low-flow dischages for recurrence intervals of 2, 10, and 20 years, and are based on approved streamflow records that include a minimum of 10 complete climatic years of record as of June 2010. Descriptive information for each of these streamgage locations, including the station number, station name, latitude, longitude, county, physiographic province, and drainage area, also is presented. The statistics are planned for incorporation into StreamStats, which is a U.S. Geological Survey Web application for obtaining stream information, and is being used by water-resource managers and decision makers in Maryland to address water-supply planning and management, water-use appropriation and permitting, wastewater and industrial discharge permitting, and setting minimum required streamflows to protect freshwater biota and ecosystems.

Open-File Report

Simulation of the interaction of karstic lakes Magnolia and Brooklyn with the upper Floridan Aquifer, southwestern Clay County, Florida

The stage of Lake Brooklyn, in southwestern Clay County, Florida, has varied over a range of 27 feet since measurements by the U.S. Geological Survey began in July 1957. The large stage changes have been attributed to the relation between highly transient surface-water inflow to the lake and subsurface conduits of karstic origin that permit a high rate of leakage from the lake to the Upper Floridan aquifer. After the most recent and severe stage decline (1990-1994), the U.S. Geological Survey began a study that entailed the use of numerical ground-water flow models to simulate the interaction of the lake with the Upper Floridan aquifer and the large fluctuations of stage that were a part of that process. A package (set of computer programs) designed to represent lake/aquifer interaction in the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model (MODFLOW-96) and the Three-Dimensional Method-of-Characteristics Solute-Transport Model (MOC3D) simulators was prepared as part of this study, and a demonstration of its capability was a primary objective of the study. (Although the official names are Brooklyn Lake and Magnolia Lake (Florida Geographic Names), in this report the local names, Lake Brooklyn and Lake Magnolia, are used.) In the simulator of lake/aquifer interaction used in this investigation, the stage of each lake in a simulation is updated in successive time steps by a budget process that takes into account ground-water seepage, precipitation upon and evaporation from the lake surface, stream inflows and outflows, overland runoff inflows, and augmentation or depletion by artificial means. The simulator was given the capability to simulate both the division of a lake into separate pools as lake stage falls and the coalescence of several pools into a single lake as the stage rises. This representational capability was required to simulate Lake Brooklyn, which can divide into as many as 10 separate pools at sufficiently low stage. In the first of two calibrated models, recharge to the water table, specified as a monthly rate, was set equal to 40 percent of the monthly rainfall rate. The specified rate of inflow to the uppermost stream segment was set equal to outflows from Lake Lowry estimated from lake stage and the 1994-97 rating table. Leakage to the intermediate and Upper Floridan aquifers was assumed to occur from the surficial aquifer system through the confining layers directly beneath deeper parts of the lake bottom. A leakance coefficient value of 0.001 feet per day per foot of thickness was used beneath Lake Magnolia, and a value of 0.005 feet per day per foot of thickness was used beneath most of Lake Brooklyn. With these values, the conductance through the confining layers beneath Lake Brooklyn was about 19 times that beneath Lake Magnolia. The simulated stages of Lake Brooklyn matched the measured stages reasonably well in the early (1957-72) and later (1990-98) parts of the simulation time period, but the match was unsatisfactory in an intermediate time period (1973-89). To resolve this discrepancy, the hypothesis was proposed that undocumented losses of water from Alligator Creek upstream from Lake Brooklyn or from the lake itself occurred between 1973 and 1989 when there was sufficient streamflow. The resulting simulation of lake stages matched the measured lake stages accurately during the entire simulation time period. The model was then revised to incorporate the assumption that only 20 percent of precipitation recharged the water table (the second calibrated model). Recalibration of the model required that leakance values for the confining units under deeper parts of the lakes also be reduced by nearly 50 percent. The stages simulated with the new parameter assumptions, but retaining the assumption of surface-water losses, were an excellent match of the measured values. The stage of Lake Magnolia was also simulated accurately. The results of sensitivity analyses show that simulated streamflow between Lakes Magnolia and Brooklyn tends to be water-budget controlled, and is not appreciably affected by the specified outflow altitude or channel characteristics of the receiving stream. To match heads measured in observation wells of the surficial aquifer network, the assigned hydraulic conductivity values were zoned, and ranged from a minimum of 4 feet per day to a maximum of 400 feet per day in the first calibrated model. These values were reduced by about 50 percent in the second calibrated model. Differences between observation wells were noted in the abruptness of changes of measured head values, and in the relation of the timing of peak measured heads and simulated peak heads. These differences seemed to be correlated with the depth of the water table below land surface. Spatially uniform values of transmissivity were specified for the intermediate (10,000 feet squared per day) and Upper Floridan (100,000 feet squared per day) aquifers. Simulated heads in the Upper Floridan aquifer layer follow the trend of the heads measured in a long-term observation well with data beginning in 1960. This result suggests that the observed head decline could be explained entirely in terms of the stage decline in Lake Brooklyn and may not indicate a regional trend.

Florida

User guide for the drawdown-limited, multi-node well (MNW) package for the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model, versions MODFLOW-96 and MODFLOW-2000

A computer program called the drawdown-limited, Multi-Node Well (MNW) Package was developed for the U.S. Geological Survey three-dimensional finite-difference modular ground-water flow model, commonly referred to as MODFLOW. The MNW Package allows MODFLOW users to simulate wells that extend beyond a single model node. Multi-node wells can simulate wells that are completed in multiple aquifers or in a single heterogeneous aquifer, partially penetrating wells, and horizontal wells. Multi-aquifer wells dynamically distribute flow between nodes under pumping, recharging, or unpumped conditions. Variations in intraborehole flow can be simulated with the MNW Package, which is limited by how finely an aquifer system has been discretized vertically. Simulated discharge from single-node and multi-node wells also can be drawdown limited, which is user specified for pumping or recharging conditions. The MNW Package also has the ability to track potential mixes of a water-quality attribute. Simulated wellbore flow can be compared with measured wellbore flow, which provides another constraint for model calibration.

Open-File Report

Liquefaction hazard maps for three earthquake scenarios for the communities of San Jose, Campbell, Cupertino, Los Altos, Los Gatos, Milpitas, Mountain View, Palo Alto, Santa Clara, Saratoga, and Sunnyvale, Northern Santa Clara County, California

Maps showing the probability of surface manifestations of liquefaction in the northern Santa Clara Valley were prepared with liquefaction probability curves. The area includes the communities of San Jose, Campbell, Cupertino, Los Altos, Los Gatos Milpitas, Mountain View, Palo Alto, Santa Clara, Saratoga, and Sunnyvale. The probability curves were based on complementary cumulative frequency distributions of the liquefaction potential index (LPI) for surficial geologic units in the study area. LPI values were computed with extensive cone penetration test soundings. Maps were developed for three earthquake scenarios, an M7.8 on the San Andreas Fault comparable to the 1906 event, an M6.7 on the Hayward Fault comparable to the 1868 event, and an M6.9 on the Calaveras Fault. Ground motions were estimated with the Boore and Atkinson (2008) attenuation relation. Liquefaction is predicted for all three events in young Holocene levee deposits along the major creeks. Liquefaction probabilities are highest for the M7.8 earthquake, ranging from 0.33 to 0.37 if a 1.5-m deep water table is assumed, and 0.10 to 0.14 if a 5-m deep water table is assumed. Liquefaction probabilities of the other surficial geologic units are less than 0.05. Probabilities for the scenario earthquakes are generally consistent with observations during historical earthquakes.

Open-File Report

The digital geologic map of Colorado in ARC/INFO format: Part A. Documentation

This geologic map was prepared as a part of a study of digital methods and techniques as applied to complex geologic maps. The geologic map was digitized from the original scribe sheets used to prepare the published Geologic Map of Colorado (Tweto 1979). Consequently the digital version is at 1:500,000 scale using the Lambert Conformal Conic map projection parameters of the state base map. Stable base contact prints of the scribe sheets were scanned on a Tektronix 4991 digital scanner. The scanner automatically converts the scanned image to an ASCII vector format. These vectors were transferred to a VAX minicomputer, where they were then loaded into ARC/INFO. Each vector and polygon was given attributes derived from the original 1979 geologic map. This database was developed on a MicroVAX computer system using VAX V 5.4 nd ARC/INFO 5.0 software. UPDATE: April 1995, The update was done solely for the purpose of adding the abilitly to plot to an HP650c plotter. Two new ARC/INFO plot AMLs along with a lineset and shadeset for the HP650C design jet printer have been included. These new files are COLORADO.650, INDEX.650, TWETOLIN.E00 and TWETOSHD.E00. These files were created on a UNIX platform with ARC/INFO 6.1.2. Updated versions of INDEX.E00, CONTACT.E00, LINE.E00, DECO.E00 and BORDER.E00 files that included the newly defined HP650c items are also included. * Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Descriptors: The Digital Geologic Map of Colorado in ARC/INFO Format Open-File Report 92-050

Colorado

The digital geologic map of Colorado in ARC/INFO format

This geologic map was prepared as a part of a study of digital methods and techniques as applied to complex geologic maps. The geologic map was digitized from the original scribe sheets used to prepare the published Geologic Map of Colorado (Tweto 1979). Consequently the digital version is at 1:500,000 scale using the Lambert Conformal Conic map projection parameters of the state base map. Stable base contact prints of the scribe sheets were scanned on a Tektronix 4991 digital scanner. The scanner automatically converts the scanned image to an ASCII vector format. These vectors were transferred to a VAX minicomputer, where they were then loaded into ARC/INFO. Each vector and polygon was given attributes derived from the original 1979 geologic map. This database was developed on a MicroVAX computer system using VAX V 5.4 nd ARC/INFO 5.0 software. UPDATE: April 1995, The update was done solely for the purpose of adding the abilitly to plot to an HP650c plotter. Two new ARC/INFO plot AMLs along with a lineset and shadeset for the HP650C design jet printer have been included. These new files are COLORADO.650, INDEX.650, TWETOLIN.E00 and TWETOSHD.E00. These files were created on a UNIX platform with ARC/INFO 6.1.2. Updated versions of INDEX.E00, CONTACT.E00, LINE.E00, DECO.E00 and BORDER.E00 files that included the newly defined HP650c items are also included. * Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Descriptors: The Digital Geologic Map of Colorado in ARC/INFO Format Open-File Report 92-050

Colorado

Map showing characteristics of lode gold in the Medford 1 degree by 2 degrees Quadrangle, Oregon-California

About 500 lode gold occurrences, prospects, and mines have been reported from the Medford 1 o x 2 o quadrangle, Oregon-California. Summary descriptions of individual lode gold deposits and districts which include those occurrences in the Medford 1 o x 2 o quadrangle have been published by Brooks and Ramp (1968) and Hotz (1971) and by others cited in the reference lists in these two publications. As part of the Medford CUSMAP project between 1974 and 1980, some of these deposits were examined in the field during geologic mapping and geochemical sampling and most of the available information on individual deposits was compiled into a CRIB computer file (Computerized Resource Information Bank). This report is a summary of the more pertinent geologic and economic characteristics of the lode gold deposits for use in the resource evaluation part of the Medford CUSMAP project. The significant features described include the distribution of deposits and their relation to gross geologic and lithologic features, production information, mineralogy of the ore and gangue, gold-silver ratios, and speculations on possible genetic models for some types of deposits.

California, Oregon

Evaluation of the Acoustic Doppler Velocity Meter for Computation of Discharge Records at Three Sites in Colorado, 2004-2005

The U.S. Geological Survey (USGS), in cooperation with the Colorado Water Conservation Board, conducted a study in 2004-2005 at three sites in Colorado: Bear Creek at Morrison, Clear Creek near Empire, and Redlands Canal near Grand Junction. The study was done to evaluate acoustic Doppler velocity meter (ADVM) technology in different hydrologic settings that are characteristic of many Colorado streamflow-gaging sites. ADVMs have been tested and used extensively in many parts of the United States by USGS but not in Colorado where relatively small, shallow, clear, coarse-bed streams that ice up in the winter may affect the ADVM suitability. In this study, ADVM instrumentation was successfully used and discharge computations compared favorably, generally within 5 to 10 percent, with conventional USGS stage/discharge methods at the three Colorado sites. However, two factors, encountered in this study, may adversely affect the use of ADVM technology in Colorado. First, for some streams, the depth required (about 1.5 feet for a side-looking instrument) cannot be met during low-flow periods of the year. Second, cold temperatures and freezing-thawing cycles can produce ice effects that could prevent collection of usable ADVM (and stage) data.

Scientific Investigations Report

Assessment of selected inorganic constituents in streams in the Central Arizona Basins study area, Arizona and northern Mexico, through 1998

Stream properties and water-chemistry constituent concentrations from data collected by the National Water-Quality Assessment and other U.S. Geological Survey water-quality programs were analyzed to (1) assess water quality, (2) determine natural and human factors affecting water quality, and (3) compute stream loads for the surface-water resources in the Central Arizona Basins study area. Stream temperature, pH, dissolved-oxygen concentration and percent saturation, and dissolved-solids, suspended-sediment, and nutrient concentration data collected at 41 stream-water quality monitoring stations through water year 1998 were used in this assessment. Water-quality standards applicable to the stream properties and water-chemistry constituent concentration data for the stations investigated in this study generally were met, although there were some exceedences. In a few samples from the White River, the Black River, and the Salt River below Stewart Mountain Dam, the pH in reaches designated as a domestic drinking water source was higher than the State of Arizona standard. More than half of the samples from the Salt River below Stewart Mountain Dam and almost all of the samples from the stations on the Central Arizona Project Canal—two of the three most important surface-water sources used for drinking water in the Central Arizona Basins study area—exceeded the U.S. Environmental Protection Agency drinking water Secondary Maximum Contaminant Level for dissolved solids. Two reach-specific standards for nutrients established by the State of Arizona were exceeded many times: (1) the annual mean concentration of total phosphorus was exceeded during several years at stations on the main stems of the Salt and Verde Rivers, and (2) the annual mean concentration of total nitrogen was exceeded during several years at the Salt River near Roosevelt and at the Salt River below Stewart Mountain Dam. Stream properties and water-chemistry constituent concentrations were related to streamflow, season, water management, stream permanence, and land and water use. Dissolved-oxygen percent saturation, pH, and nutrient concentrations were dependent on stream regulation, stream permanence, and upstream disposal of wastewater. Seasonality and correlation with streamflow were dependant on stream regulation, stream permanence, and upstream disposal of wastewater. Temporal trends in streamflow, stream properties, and water-chemistry constituent concentrations were common in streams in the Central Arizona Basins study area. Temporal trends in the streamflow of unregulated perennial reaches in the Central Highlands tended to be higher from 1900 through the 1930s, lower from the 1940s through the 1970s, and high again after the 1970s. This is similar to the pattern observed for the mean annual precipitation for the Southwestern United States and indicates long-term trends in flow of streams draining the Central Highlands were driven by long-term trends in climate. Streamflow increased over the period of record at stations on effluent-dependent reaches as a result of the increase in the urban population and associated wastewater returns to the Salt and Gila Rivers in the Phoenix metropolitan area and the Santa Cruz River in the Tucson metropolitan area. Concentrations of dissolved solids decreased in the Salt River below Stewart Mountain Dam and in the Verde River below Bartlett Dam. This decrease represents an improvement in the water quality and resulted from a concurrent increase in the amount of runoff entering the reservoirs. Stream loads of water-chemistry constituents were compared at different locations along the streams with one another, and stream loads were compared to upstream inputs of the constituent from natural and anthropogenic sources to determine the relative importance of different sources and to determine the fate of the water-chemistry constituent. Of the dissolved solids transported into the Basin and Range Lowlands each year from the Central Arizona Project Canal and from streams draining the Central Highlands, about 1.2 billion kilograms accumulated in the soil, unsaturated zone, and aquifers in agricultural and urban areas as a result of irrigating crops and urban vegetation. Stream loads of phosphorus decreased from the 91st Avenue Wastewater-Treatment Plant downstream to the Gila River at Gillespie Dam, probably as a result of adsorption of phosphorus to the streambed sediments. In this same reach, stream loads of nitrogen increased, most likely because of inputs from fertilizers. The annual mass of nitrogen and phosphorus input to developed basins from quantifiable sources was much larger than the mass input to basins that had little or no municipal or agricul-tural development. These computed inputs exclude the mass of nitrogen and phosphorus from sources such as geologic formations and soils that could not be quantified. The quantifiable annual inputs of nitrogen and phosphorus for the upper Salt River Basin and the upper Verde River Basin were similar to those for the West Clear Creek Basin. This similarity suggests that the small amount of municipal and agricultural development in the upper Salt River and the upper Verde River Basins did not greatly change the basin input flux. For basins with minimal urban and agricultural development, the largest quantifiable source of nitrogen was precipitation, and the largest source of phosphorus was human bodily waste treated by sewer and septic systems. This was in contrast to developed basins, for which fertilizer was the largest quantifiable source of both nutrients. For most basins examined, quantifiable inputs of nitrogen and phosphorus from nonpoint sources were greater than inputs from point sources. This relation emphasizes the importance of land- and water-management policies that protect surface-water resources from nonpoint sources of nutrients as well as from point sources. The amount of nitrogen and phosphorus transpor-ted out of basins was a small fraction of the total for the quantifiable inputs. This result indicated that most of the nutrients input to basins were not transported out of the basins in surface water, but rather were transported to the subsurface (the soil, unsaturated zone, or aquifer), released to the atmosphere (such as volatilized ammonia), or incorporated into the biomass.

Central Arizona Basins study area

Overview of selected surrogate technologies for high-temporal resolution suspended-sediment monitoring

Traditional methods for characterizing selected properties of suspended sediments in rivers are being augmented and in some cases replaced by cost-effective surrogate instruments and methods that produce a temporally dense time series of quantifiably accurate data for use primarily in sediment-flux computations. Turbidity is the most common such surrogate technology, and the first to be sanctioned by the U.S. Geological Survey for use in producing data used in concert with water-discharge data to compute sediment concentrations and fluxes for storage in the National Water Information System. Other technologies, including laser-diffraction, digital photo-optic, acoustic-attenuation and backscatter, and pressure-difference techniques are being evaluated for producing reliable sediment concentration and, in some cases, particle-size distribution data. Each technology addresses a niche for sediment monitoring. Their performances range from compelling to disappointing. Some of these technologies have the potential to revolutionize fluvial-sediment data collection, analysis, and availability.

Conference Paper

Coal-seismic, desktop computer programs in BASIC; Part 5, Perform X-square T-square analysis and plot normal moveout lines on seismogram overlay

Processing of data taken with the U.S. Geological Survey's coal-seismic system is done with a desktop, stand-alone computer. Programs for this computer are written in the extended BASIC language used by the Tektronix 4051 Graphic System. This report presents computer programs to perform X-square/T-square analyses and to plot normal moveout lines on a seismogram overlay.

Open-File Report

Resident research associateships, postdoctoral research awards 1989: opportunities for research at the U.S. Geological Survey, U.S. Department of the Interior

The scientists of the U.S. Geological Survey are engaged in a wide range of geologic, geophysical, geochemical, hydrologic, and cartographic programs, including the application of computer science to them. These programs offer exciting possibilities for scientific achievement and professional growth to young scientists through participation as Research Associates.

Monograph

Sources of groundwater and characteristics of surface-water recharge at Bell, White, and Suwannee Springs, Florida, 2012–13

Discharge from springs in Florida is sourced from aquifers, such as the Upper Floridan aquifer, which is overlain by an upper confining unit that locally can have properties of an aquifer. Water levels in aquifers are affected by several factors, such as precipitation, recharge, and groundwater withdrawals, which in turn can affect discharge from springs. Therefore, identifying groundwater sources and recharge characteristics can be important in assessing how these factors might affect flows and water levels in springs and can be informative in broader applications such as groundwater modeling. Recharge characteristics include the residence time of water at the surface, apparent age of recharge, and recharge water temperature. The groundwater sources and recharge characteristics of three springs that discharge from the banks of the Suwannee River in northern Florida were assessed for this study: Bell Springs, White Springs, and Suwannee Springs. Sources of groundwater were also assessed for a 150-foot-deep well finished within the Upper Floridan aquifer, hereafter referred to as the UFA well. Water samples were collected for geochemical analyses in November 2012 and October 2013 from the three springs and the UFA well. Samples were analyzed for a suite of major ions, dissolved gases, and isotopes of sulfur, strontium, oxygen, and hydrogen. Daily means of water level and specific conductance at White Springs were continuously recorded from October 2012 through December 2013 by the Suwannee River Water Management District. Suwannee River stage at White Springs was computed on the basis of stage at a U.S. Geological Survey streamgage about 2.4 miles upstream. Water levels in two wells, located about 2.5 miles northwest and 13 miles southeast of White Springs, were also used in the analyses. Major ion concentrations were used to differentiate water from the springs and Upper Floridan aquifer into three groups: Bell Springs, UFA well, and White and Suwannee Springs. When considered together, evidence from water-level, specific conductance, major-ion concentration, and isotope data indicated that groundwater at Bell Springs and the UFA well was a mixture of surface water and groundwater from the upper confining unit, and that groundwater at White and Suwannee Springs was a mixture of surface water, groundwater from the upper confining unit, and groundwater from the Upper Floridan aquifer. Higher concentrations of magnesium in groundwater samples at the UFA well than in samples at Bell Springs might indicate less mixing with surface water at the UFA well than at Bell Springs. Characteristics of surface-water recharge, such as residence time at the surface, apparent age, and recharge water temperature, were estimated on the basis of isotopic ratios, and dissolved concentrations of gases such as argon, tritium, and sulfur hexafluoride. Oxygen and deuterium isotopic ratios were consistent with rapid recharge by rainwater for samples collected in 2012, and longer residence time at the surface (ponding) for samples collected in 2013. Apparent ages of groundwater samples, computed on the basis of tritium activity and sulfur hexafluoride concentration, indicated groundwater recharge occurred after the late 1980s; however, the estimated apparent ages likely represent the average of ages of multiple sources. Recharge since the 1980s is consistent with groundwater from shallow sources, such as the upper confining unit and Upper Floridan aquifer. Recharge water temperature computed for the three springs and UFA well averaged 20.1 degrees Celsius, which is similar to the mean annual air temperature of 20.6 degrees Celsius at a nearby weather station for 1960–2014.

Florida

Documentation of a computer program to simulate stream-aquifer relations using a modular, finite-difference, ground-water flow model

Computer models are widely used to simulate groundwater flow for evaluating and managing the groundwater resource of many aquifers, but few are designed to also account for surface flow in streams. A computer program was written for use in the US Geological Survey modular finite difference groundwater flow model to account for the amount of flow in streams and to simulate the interaction between surface streams and groundwater. The new program is called the Streamflow-Routing Package. The Streamflow-Routing Package is not a true surface water flow model, but rather is an accounting program that tracks the flow in one or more streams which interact with groundwater. The program limits the amount of groundwater recharge to the available streamflow. It permits two or more streams to merge into one with flow in the merged stream equal to the sum of the tributary flows. The program also permits diversions from streams. The groundwater flow model with the Streamflow-Routing Package has an advantage over the analytical solution in simulating the interaction between aquifer and stream because it can be used to simulate complex systems that cannot be readily solved analytically. The Streamflow-Routing Package does not include a time function for streamflow but rather streamflow entering the modeled area is assumed to be instantly available to downstream reaches during each time period. This assumption is generally reasonable because of the relatively slow rate of groundwater flow. Another assumption is that leakage between streams and aquifers is instantaneous. This assumption may not be reasonable if the streams and aquifers are separated by a thick unsaturated zone. Documentation of the Streamflow-Routing Package includes data input instructions; flow charts, narratives, and listings of the computer program for each of four modules; and input data sets and printed results for two test problems, and one example problem. (Lantz-PTT)

Open-File Report

Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model

Removal of ground water by pumping from aquifers may result in compaction of compressible fine-grained beds that are within or adjacent to the aquifers. Compaction of the sediments and resulting land subsidence may be permanent if the head declines result in vertical stresses beyond the previous maximum stress. The process of permanent compaction is not routinely included in simulations of ground-water flow. To simulate storage changes from both elastic and inelastic compaction, a computer program was written for use with the U.S. Geological Survey modular finite-difference ground- water flow model. The new program, the Interbed-Storage Package, is designed to be incorporated into this model. In the Interbed-Storage Package, elastic compaction or expansion is assumed to be proportional to change in head. The constant of proportionality is the product of the skeletal component of elastic specific storage and the thickness of the sediments. Similarly, inelastic compaction is assumed to be proportional to decline in head. The constant of proportionality is the product of the skeletal component of inelastic specific storage and the thickness of the sediments. Storage changes are incorporated into the ground-water flow model by adding an additional term to the right-hand side of the flow equation. Within a model time step, the package appropriately apportions storage changes between elastic and inelastic components on the basis of the relation of simulated head to the previous minimum (preconsolidation) head. Two tests were performed to verify that the package works correctly. The first test compared model-calculated storage and compaction changes to hand-calculated values for a three-dimensional simulation. Model and hand-calculated values were essentially equal. The second test was performed to compare the results of the Interbed-Storage Package with results of the one-dimensional Helm compaction model. This test problem simulated compaction in doubly draining confining beds stressed by head changes in adjacent aquifers. The Interbed-Storage Package and the Helm model computed essentially equal values of compaction. Documentation of the Interbed-Storage Package includes data input instructions, flow charts, narratives, and listings for each of the five modules included in the package. The documentation also includes an appendix describing input instructions and a listing of a computer program for time-variant specified-head boundaries. That package was developed to reduce the amount of data input and output associated with one of the Interbed-Storage Package test problems.

Techniques of Water-Resources Investigations

A computer program for estimating instream travel times and concentrations of a potential contaminant in the Yellowstone River, Montana

The Yellowstone River is very important in a variety of ways to the residents of southeastern Montana; however, it is especially vulnerable to spilled contaminants. In 2004, the U.S. Geological Survey, in cooperation with Montana Department of Environmental Quality, initiated a study to develop a computer program to rapidly estimate instream travel times and concentrations of a potential contaminant in the Yellowstone River using regression equations developed in 1999 by the U.S. Geological Survey. The purpose of this report is to describe these equations and their limitations, describe the development of a computer program to apply the equations to the Yellowstone River, and provide detailed instructions on how to use the program. This program is available online at [http://pubs.water.usgs.gov/sir2006-5057/includes/ytot.xls]. The regression equations provide estimates of instream travel times and concentrations in rivers where little or no contaminant-transport data are available. Equations were developed and presented for the most probable flow velocity and the maximum probable flow velocity. These velocity estimates can then be used to calculate instream travel times and concentrations of a potential contaminant. The computer program was developed so estimation equations for instream travel times and concentrations can be solved quickly for sites along the Yellowstone River between Corwin Springs and Sidney, Montana. The basic types of data needed to run the program are spill data, streamflow data, and data for locations of interest along the Yellowstone River. Data output from the program includes spill location, river mileage at specified locations, instantaneous discharge, mean-annual discharge, drainage area, and channel slope. Travel times and concentrations are provided for estimates of the most probable velocity of the peak concentration and the maximum probable velocity of the peak concentration. Verification of estimates of instream travel times and concentrations for the Yellowstone River requires information about the flow velocity throughout the 520 mi of river in the study area. Dye-tracer studies would provide the best data about flow velocities and would provide the best verification of instream travel times and concentrations estimated from this computer program; however, data from such studies does not currently (2006) exist and new studies would be expensive and time-consuming. An alternative approach used in this study for verification of instream travel times is based on the use of flood-wave velocities determined from recorded streamflow hydrographs at selected mainstem streamflow-gaging stations along the Yellowstone River. The ratios of flood-wave velocity to the most probable velocity for the base flow estimated from the computer program are within the accepted range of 2.5 to 4.0 and indicate that flow velocities estimated from the computer program are reasonable for the Yellowstone River. The ratios of flood-wave velocity to the maximum probable velocity are within a range of 1.9 to 2.8 and indicate that the maximum probable flow velocities estimated from the computer program, which corresponds to the shortest travel times and maximum probable concentrations, are conservative and reasonable for the Yellowstone River.

Scientific Investigations Report