Search USGSSearch

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 73 records · Page 4Linked to original sources

Documentation of computer programs to compute and display pathlines using results from the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model

A particle tracking post-processing package was developed to compute three-dimensional path lines based on output from steady-state simulations obtained with the U.S. Geological Survey modular 3-dimensional finite difference groundwater flow model. The package consists of two FORTRAN 77 computer programs: (1) MODPATH, which calculates pathlines, and (2) MODPATH-PLOT, which presents results graphically. MODPATH uses a semi-analytical particle tracking scheme. The method is based on the assumption that each directional velocity component varies linearly within a grid cell in its own coordinate direction. This assumption allows an analytical expression to be obtained describing the flow path within a grid cell. Given the initial position of a particle anywhere in a cell, the coordinates of any other point along its path line within the cell, and the time of travel between them, can be computed directly. Data is input to MODPATH and MODPATH-PLOT through a combination of files and interactive dialogue. Examples of how to use MODPATH and MODPATH-PLOT are provided for a sample problem. Listings of the computer codes and detailed descriptions of input data format and program options are also presented. (Author 's abstract)

Open-File Report

Evaluation of factors that influence estimated zones of transport for six municipal wells in Clark County, Washington

A ground-water flow model was used in conjunction with particle tracking to estimate zones of transport for six municipal well sites in Clark County, Washington. A zone of transport for a well is a three-dimensional volume within a ground-water system that contains all of the ground water that will discharge from that well within a specified time period. All of the zones of transport for a well compose the zone of contribution for the well. Zones of transport and contribution are important considerations in the delineation of wellhead-protection areas. Hydrogeologic factors, such as hydraulic conductivity and porosity, influence the shape and size of the zones of transport, and, therefore, uncertainty in these and other factors can lead to uncertainty in the delineation of the zones of transport. The sensitivity of the zones of transport to uncertainty in selected hydrogeologic factors was evaluated for the six wells. Estimates of the zones of transport were delineated by the U.S. Geological Survey program MODTOOLS from three-dimensional pathlines computed by the U.S. Geological Survey program MODPATH. Input to MODPATH came from steady-state simulations calculated by the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model, MODFLOW. Three-dimensional modeling is the best method for delineating zones of transport within stratigraphically complex, heterogeneous, anisotropic aquifers that have complex boundary conditions such as streams and multiple, simultaneously discharging wells. In this study, zones of transport were delineated by using simulated particle locations computed from the results of a three-dimensional steady-state regional model for 0-0.5, 0.5-1, 1-5, 5-10, 10-20, and 20-50 year travel times to the selected wells. Zones of transport for a well were delineated by tracking particles along pathlines in the reverse direction of ground-water flow. Sensitivity of the zones of transport to change in the discharge rate of the selected well, porosity, and hydraulic conductivity, as well as to the presence or absence of interfering wells, was evaluated at six well sites to evaluate the effect of uncertainties in these factors on the size and shape of zones of transport. Uncertainty in porosity contributed the most to the uncertainty in delineating the zones of transport. Uncertainty in other factors, such as well discharge rate and horizontal hydraulic conductivity, had measurable effects on the zones of transport, but errors introduced through these factors were less significant. Insight into the causes of the changes in the size and shape of the zones of transport to varying conditions was gained by evaluating the simulated water budget and ground-water levels in the vicinity of the well. Changes in the simulated water budget and ground-water levels provided information to better understand the effects of uncertainties in the data on simulation results.The results of this study suggest that ground-water velocity is the underlying control on the size of the zones of transport. The regional hydraulic gradient is the most significant factor controlling the shape and orientation of the zones of transport. Spatial variation in recharge, discharge, and hydraulic properties can also affect the shape of the zones of transport, however. Underestimation of porosity or overestimation of horizontal hydraulic conductivity leads to overestimation of ground-water velocity and overestimation of the size of zones of transport. Overestimation of porosity or underestimation of horizontal hydraulic conductivity leads to underestimation of ground-water velocity and underestimation of the size of zones of transport. Well discharge rate affects ground-water velocities near the well. Underestimation of discharge (and therefore velocities) will result in underestimation of the size of the zones of transport. The sensitivity of estimated zones of transport to uncertainty in parameters such as porosity and horizontal hydraulic conductivity is a function of the well discharge rate and the proximity of the well to boundaries, such as streams and rivers.

Washington

Computer programs for obtaining and analyzing daily mean streamflow data from the U.S. Geological Survey National Water Information System web site

Five computer programs were developed for obtaining and analyzing streamflow from the National Water Information System (NWISWeb). The programs were developed as part of a study by the U.S. Geological Survey, in cooperation with the Federal Highway Administration, to develop a stochastic empirical loading and dilution model. The programs were developed because reliable, efficient, and repeatable methods are needed to access and process streamflow information and data. The first program is designed to facilitate the downloading and reformatting of NWISWeb streamflow data. The second program is designed to facilitate graphical analysis of streamflow data. The third program is designed to facilitate streamflow-record extension and augmentation to help develop long-term statistical estimates for sites with limited data. The fourth program is designed to facilitate statistical analysis of streamflow data. The fifth program is a preprocessor to create batch input files for the U.S. Environmental Protection Agency DFLOW3 program for calculating low-flow statistics. These computer programs were developed to facilitate the analysis of daily mean streamflow data for planning-level water-quality analyses but also are useful for many other applications pertaining to streamflow data and statistics. These programs and the associated documentation are included on the CD-ROM accompanying this report. This report and the appendixes on the CD-ROM describe the implementation and use of the programs and the interpretation of results from the programs. The body of this report provides an overview of the five programs included on this CD-ROM. The appendixes are the software manuals for each program. These manuals describe statistical and numerical methods used to implement each program, input-file formats, output-file formats, installation of the programs, and use of the programs. Each appendix is written as a self-contained manual because each program may have many uses alone or in tandem with other programs on the CD-ROM. Each of these programs uses graphical user interface that follows standard Microsoft Windows interface conventions.

Open-File Report

Nonstationary flood frequency analysis using regression in the north-central United States

Traditional flood frequency methods assume that the statistical properties of peak streamflow do not change with time and may not be appropriate for many areas in the north-central United States. This study examines a nonstationary flood frequency analysis method that uses ordinary least squares linear regression to estimate flood magnitudes at U.S. Geological Survey streamgages that exhibit trends and change points in a nine-State region including Montana, North Dakota, South Dakota, Minnesota, Illinois, Iowa, Wisconsin, Missouri, and Michigan. Additionally, an extension of this method is introduced, which enables nonstationary flood frequency based on a statistical relation with a stochastic climate predictor. Estimates of the 1-percent annual exceedance probability flood using regression equations to adjust for conditions in 2020 were computed at U.S. Geological Survey streamgages across the study area. Regression equations used either a time index or a climate variable as the explanatory variable for changes in peak streamflow. Of 153 candidate streamgages, the assumptions of time-adjusted analyses were met at 137 streamgages. Climate-adjusted flood frequency analyses were applicable at 98 streamgages based on annual precipitation, annual temperature, or annual snowfall. Time- and climate-adjusted methods produced similar estimates of the 1-percent annual exceedance probability flood magnitude at streamgages where both methods were applicable. Nonstationary estimates of the 1-percent annual exceedance probability flood were primarily greater than stationary estimates in eastern North and South Dakota, Minnesota, Iowa, Illinois, and parts of Missouri and less than stationary estimates in Montana, western North and South Dakota, and Wisconsin. The largest differences between stationary and nonstationary flood estimates were in North and South Dakota and Minnesota.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon

Ground-water quality atlas of Wisconsin

This report summarizes data on ground-water quality stored in the U.S. Geological Survey's computer system (WATSTORE). The summary includes water quality data for 2,443 single-aquifer wells, which tap one of the State's three major aquifers (sand and gravel, Silurian dolomite, and sandstone). Data for dissolved solids, hardness, alkalinity, calcium, magnesium, sodium, potassium, iron, manganese, sulfate, chloride, fluoride, and nitrate are summarized by aquifer and by county, and locations of wells for which data are available 1 are shown for each aquifer. Calcium, magnesium, and bicarbonate (the principal component of alkalinity) are the major dissolved constituents in Wisconsin's ground water. High iron concentrations and hardness cause ground-water quality problems in much of the State. Statewide ,summaries of trace constituent (selected trace metals; arsenic, boron, and organic carbon) concentrations show that these constituents impair water quality in only a few isolated wells.

Wisconsin

Index of surface-water records, part 3, Ohio River basin, to September 30, 1948

This report summarizes ground-water-quality data, for Wisconsin, stored in the U.S. Geological Survey 's computer system (WATSTORE). The summary includes water-quality data for 2,443 wells which tap one of the State 's three major aquifers (sand and gravel, Silurian dolomite, and sandstone). Data for dissolved solids, hardness, alkalinity, calcium, magnesium, sodium, potassium, iron, manganese, sulfate, chloride, fluoride, and nitrate are summarized by aquifer and by county, and locations of wells where data are available are shown for each aquifer. Calcium, magnestium, and bicarbonate (the principal component of alkalinity) are the major dissolved constituents present in Wisconsin 's ground water. High iron concentrations and hardness cause ground-water-quality problems in much of the State. Statewide summaries of trace constituent (selected trace metals, arsenic, boron, and organic carbon) concentrations show that they impair water quality in only a few isolated wells. A bibliography of published reports that describe or summarize various aspects of ground-water quality in Wisconsin is included. (USGS)

Wisconsin

Graphic and algebraic solutions of the discordant lead-uranium age problem

Uranium-bearing minerals that give lead-uranium and lead—lead ages that are essentially in agreement, i.e. concordant, generally are considered to have had a relatively simple geologic history and to have been unaltered since their deposition. The concordant ages obtained on such materials are, therefore, assumed to approach closely the actual age of the minerals. Many uranium-bearing samples, particularly uranium ores, give the following discordant age sequences; Pb 206 U 238 &lt; Pb 207 U 235 &#x2AA1; Pb 207 Pb 206 "> Pb 206 U 238 <Pb 207 U 235 ⪡Pb 207 Pb 206 or, less frequently, Pb 207 Pb 206 &#x2AA1; Pb 207 U 235 &lt; Pb 206 U 238 "> Pb 207 Pb 206 ⪡Pb 207 U 235 <Pb 206 U 238 . These discordant age sequences have been attributed most often to uncertainties in the common lead correction, selective loss of radio-active daughter products, loss or gain of lead or uranium, or contamination by an older generation of radiogenic lead. The evaluation of discordant lead isotope age data may be separated into two operations. The first operation, with which this paper is concerned, is mechanical in nature and involves the calculation of the different possible concordant ages corresponding to the various processes assumed to have produced the discordant ages. The second operation is more difficult to define and requires, in part, some personal judgement. It includes a synthesis of the possible concordant age solutions with other independent geologic and isotopic evidence. The concordant age ultimately chosen as most acceptable should be consistent not only with the known events in the geologic history of the area, the age relations of the enclosing rocks, and the mineralogic and paragenetic evidence, but also with other independent age measurements and the isotopic data obtained on the lead in related or associated non-radioactive minerals. The calculation of the possible concordant ages from discordant age data has been greatly simplified by Wetherill's graphical method of plotting the mole ratios of radiogenic Pb 206 U 238 "> Pb 206 U 238 ( N 206 N 238 "> N 206 N 238 ) vs. radiogenic Pb 207 U 235 "> Pb 207 U 235 ( N 207 N 235 "> N 207 N 235 ) after correcting for the contaminating common Pb 206 and Pb 207 . The linear relationships noted in this graphical procedure have been extended to plots of the mole ratios of total Pb 206 U 238 "> Pb 206 U 238 ( t N 206 N 238 "> t N 206 N 238 ) vs. total Pb 207 U 235 "> Pb 207 U 235 ( t N 207 N 235 "> t N 207 N 235 ). This modification permits the calculation of concordant ages for unaltered samples using only the Pb 207 Pb 206 "> Pb 207 Pb 206 ratio of the contaminating common lead. If isotopic data are available for two samples of the same age, x and y , from the same or related deposits or outcrops, graphs of the normalized difference ratios [ ( N 206 N 204 )x &#x2212; ( N 206 N 204 )y ( N 238 N 204 )x &#x2212;( N 238 N 204 )y ] vs. [ ( N 207 N 204 )x &#x2212; ( N 207 N 204 )y ( N 235 N 204 )x &#x2212;( N 235 N 204 )y ] "> [(N 206 N 204 )x − (N 206 N 204 )y(N 238 N 204 )x −(N 238 N 204 )y] vs. [(N 207 N 204 )x − (N 207 N 204 )y(N 235 N 204 )x −(N 235 N 204 )y] can give concordant ages corrected for unknown amounts of a common lead with an unknown Pb 207 / Pb 206 ratio. (If thorium is absent the difference ratios may be normalized with the more abundant index isotope, Pb 208 .) Similar plots of tho normalized, difference ratios for three genetically related samples ( x − y ) and( x − z ), will give concordant ages corrected, in addition, for either one unknown period of past alteration or initial contamination by an older generation of radiogenic lead of unknown Pb 207 /Pb 206 ratio. Practical numerical solutions for many of tho concordant age calculations are not currently available. However, the algebraic equivalents of these new graphical methods give equations which may be programmed for computing machines. For geologically probable parameters the equations of higher order have two positive real roots that rapidly converge on the exact concordant ages corrected for original radiogenic lead and for loss or gain of lead or uranium. Modifications of these general age equations expanded only to the second degree have been derived for use with desk calculators. These graphical and algebraic methods clearly suggest both the type and minimum number of samples necessary for adequate mathematical analysis of discordant lead isotope age data. This mathematical treatment also makes it clear that discordant lead isotope data alone cannot provide the basis for the choice of one of the possible concordant age solutions. The new equations, in particular, provide an incentive to improve our physical constants, analytical techniques and sampling methods in order that we may derive all of the useful geologic information that is available in a comprehensive lead isotope age study.

Geochimica et Cosmochimica Acta

Catalog of US GeoData

The development of geographic information systems (GIS) is a rapidly growing industry that supports natural resources, studies, land management, environmental analysis, and urban and transporation planning. The increasing use of computers for storing and analyzing earth science information has greatly expanded the demand for digital cartographic and geographic data. Digital cartography involves the collection, storage, processing, analysis, and display of map data with the aid of computers. The U.S. Geological Survey (USGS), the Nation's largest earth science research agency, through its National Mapping Program, has expanded digital cartography operations to include the collection of elevation, planimetric, land use and land cover, and geographic names information in digital form. This digital information is available on 9-track magnetic tapes and, in the case of 1:2,000,000-scale planimetric digital line graph data, in Compact Disc Read Only Memory (CD-ROM) format. Digital information can be used with all types of geographic and land information systems.

General Information Product

US GeoData: Digital cartographic and geographic data

The increasing use of computers for storing and analyzing earth science information has sparked a growth in the demand for various types of cartographic data in digital form. The production of map data in computerized form is called digital cartography, and it involves the collection, storage, processing, analysis, and display of map data with the aid of computers. The U.S. Geological Survey, the Nation's largest earth science research agency, has expanded its national mapping program to incorporate operations associated with digital cartography, including the collection of planimetric, elevation, and geographic names information in digital form. This digital information is available for use in meeting the multipurpose needs and applications of the map user community.

General Information Product

Sources of water to the Rio Grande upstream from San Marcial, New Mexico

The Rio Grande watershed is a complex hydrologic system that includes numerous tributaries, inflow from transmountain diversions, irrigation diversions, agricultural return lows, reservoirs, and ground-water inflows and outflows. Many people depend on and are affected by the Rio Grande, which is the largest river of the surface-water system draining the Rio Grande watershed. To provide information to further the understanding of the Rio Grande surface-water system, the U.S. Geological Survey, in cooperation with the City of Albuquerque, evaluated streamflow data computed by the U.S. Geological Survey at 60 streamflowgaging stations. This fact sheet examines the sources of water to the Rio Grande and the "water balance" of the Rio Grande surface-water system upstream from San Marcial, New Mexico, by comparing the mean annual discharge (streamflow) at different stations.

Fact Sheet

Flood Magnitude and Frequency of the Delaware River in New Jersey, New York, and Pennsylvania

From September 2004 to June 2006, the Delaware River in New Jersey, New York, and Pennsylvania experienced three major floods that caused extensive damage. The Federal Emergency Management Agency (FEMA) needed updated information on the flood magnitude and frequency for the eight active streamflow-gaging stations along the main stem Delaware River in New Jersey, New York, and Pennsylvania that included the three recent floods in order to update its flood insurance studies. Therefore, the U.S. Geological Survey (USGS) computed updated flood magnitude and frequency values following the guidelines published by the Interagency Advisory Committee on Water Data in its Bulletin 17B. The updated flood-frequency values indicate that the recurrence interval of the September 2004 flood ranged from 20 to 35 years, the recurrence interval of the April 2005 flood ranged from 40 to 70 years, and the recurrence interval of the June 2006 flood ranged from 70 to greater than 100 years. Examination of trends in flood discharges indicate no statistically significant trends in peak flows during the period of record for any of the eight streamflow-gaging stations.

Open-File Report

Floods of August 21-24, 2007, in Northwestern and North-Central Ohio

Heavy rains in northwestern and north-central Ohio on August 19-22, 2007, caused severe flooding and widespread damages to residential, public, and commercial structures in the communities of Bluffton, Bucyrus, Carey, Columbus Grove, Crestline, Findlay, Mansfield, Ottawa, and Shelby. On August 27, 2007, the Federal Emergency Management Agency (FEMA) issued a notice of a Presidential declaration of a major disaster affecting Allen, Crawford, Hancock, Hardin, Putnam, Richland, Seneca, and Wyandot Counties as a result of the severe flooding. Rainfall totals for most of the flooded area were 3 to 5 in., with some locations reporting as much as 8 to 10 in. Three National Weather Service (NWS) gages in the area indicated a rainfall recurrence interval of greater than 1,000 years, and two indicated a recurrence interval between 500 and 1,000 years. Total damages are estimated at approximately $290 million, with 8,205 residences registering for financial assistance. The U.S. Geological Survey (USGS) computed flood recurrence intervals for peak streamflows at 22 streamgages and 8 ungaged sites in and around the area of major flooding. The peak streamflows at Sandusky River near Bucyrus streamgage and at seven of the eight ungaged sites had estimated recurrence intervals of greater than 500 years. The USGS located and surveyed 421 high-water marks and plotted high-water profiles for approximately 44.5 miles of streams throughout the nine communities.

Open-File Report

Simulation of groundwater flow in a volatile organic compound-contaminated area near Bethpage, Nassau County, New York: A discussion of modeling considerations

The 2010 Bethpage groundwater-flow model (ARCADIS, 2010) was based on a steady state assumption. Although it is widely acknowledged that significant water-level changes have occurred in the past, the reviewed model does not represent changing water levels. The steady state approach limits the effectiveness of the following: 1. identification of sources of contamination, 2. analysis of model accuracy, 3. model calibration, and 4. simulations of future scenarios. Future plume movement was simulated in an incomplete manner through an unchanging groundwater-flow field. Available time-series information on temporal variation of factors affecting groundwater-flow dynamics includes: 1. public-supply pumping, 2. groundwater discharges from systems remediating volatile organic compound (VOC) plumes, 3. recharge and precipitation rates, and 4. water levels and streamflows. Transient phenomena that might be useful in future hypothetical simulations include pumping variations, redirection of containment-system waters for industrial use, and climate-change scenarios. Public-domain computer programs, U.S. Geological Survey guidance reports on transient-state calibration and uncertainty methods (Doherty and Hunt, 2010), and additional local and regional datasets are available to provide additional confidence in model evaluations and allow better evaluation of their limitations.

New York

Construction data and retrieval procedures for selected wells drilled from 1985 through 1987 at Oak Ridge National Laboratory, Tennessee

Twenty-eight wells were constructed by the U. S. Geological Survey for use in describing the groundwater flow system in Melton Valley, at the Oak Ridge National Laboratory in eastern Tennessee. The wells were installed at 18 locations in Melton Valley and along the Clinch River during the period 1985 through 1987. During the same period, 19 wells were constructed by Oak Ridge National Laboratory at 7 locations in or near radioactive-waste burial grounds in Melton Valley. Construction data for all 47 wells are in the U.S. Geological Survey Groundwater Site Inventory data system, where information is also stored for 450 wells that were completed at the laboratory in earlier years. The data can be electronically retrieved by personnel who have access to the U.S. Geological Survey Prime computer located in Nashville, Tennessee, and retrieval procedures are given in the report. (USGS)

Tennessee

Hydrogeology, distribution, and volume of saline groundwater in the southern midcontinent and adjacent areas of the United States

The hydrogeology, distribution, and volume of saline water in 22 aquifers in the southern midcontinent of the United States were evaluated to provide information about saline groundwater resources that may be used to reduce dependency on freshwater resources. Those aquifers underlie six States in the southern midcontinent—Arkansas, Kansas, Louisiana, Missouri, Oklahoma, and Texas—and adjacent areas including all or parts of Alabama, Colorado, Florida, Illinois, Kentucky, Mississippi, Nebraska, New Mexico, South Dakota, Tennessee, and Wyoming and some offshore areas of the Gulf of Mexico. Saline waters of the aquifers were evaluated by defining salinity zones; digitizing data, primarily from the Regional Aquifer-System Analysis Program of the U.S. Geological Survey; and computing the volume of saline water in storage. The distribution of saline groundwater in the southern midcontinent is substantially affected by the hydrogeology and groundwater-flow systems of the aquifers. Many of the aquifers in the southern midcontinent are underlain by one or more aquifers, resulting in vertically stacked aquifers containing groundwaters of varying salinity. Saline groundwater is affected by past and present hydrogeologic conditions. Spatial variation of groundwater salinity in the southern midcontinent is controlled primarily by locations of recharge and discharge areas, groundwater-flow paths and residence time, mixing of freshwater and saline water, and interactions with aquifer rocks and sediments. The volume calculations made for the evaluated aquifers in the southern midcontinent indicate that about 39,900 million acre-feet (acre-ft) of saline water is in storage. About 21,600 million acre-ft of the water in storage is slightly to moderately saline (1,000–10,000 milligrams per liter [mg/L] dissolved solids), and about 18,300 million acre-ft is very saline (10,000–35,000 mg/L dissolved solids). The largest volumes of saline water are in the coastal lowlands (about 16,300 million acre-ft), Mississippi embayment and Texas coastal uplands (about 12,000 million acre-ft), and Great Plains (about 8,170 million acre-ft) aquifer systems. Of the 22 aquifers evaluated in this report, the Maha aquifer in the Great Plains aquifer system contains both the largest total volume of saline water (about 6,280 million acre-ft) and the largest volume of slightly to moderately saline water (about 5,150 million acre-ft).

Alabama;Arkansas;Colorado;Florida;Illinois;Kansas;

Traveltime, longitudinal-dispersion, reaeration, and basin characteristics of the White River, Colorado and Utah

Traveltime and longitudinal-dispersion characteristics were measured for the White River using dye tracers. Discharges ranged from 281 to 1,840 cubic feet per second and velocities ranged from 1.26 to 3.17 miles per hour. Traveltimes were determined for discharges other than measured discharges by a graphical method and a linear-regression method. Longitudinal-dispersion coefficients ranged from 284 square feet per second at a discharge of 539 cubic feet per second to 5,430 square feet per second at a discharge of 1,580 cubic feet per second. Reaeration was measured in four reaches of the White River during a medium-flow period in August 1982. Reaeration coefficients at 20 degrees Celsius ranged from 5.3 to 25.3 per day. The results of a comparison with measured reaeration coefficients and reaeration coefficients predicted using empirical equations showed that the most accurate equations were by Bennett and Rathbun (1972) and Isaacs and Gaudy (1968). Basin characteristics were computed using U.S. Geological Survey topographic maps, precipitation data from the National Weather Service, and aerial photographs taken on September 11, 1981. (USGS)

Water-Resources Investigations Report

Quantity and quality of urban runoff from the Chester Creek basin Anchorage, Alaska

Urbanization has affected both the flow characteristics and water quality of streams in the Chester Creek basin, of Anchorage, Alaska. Peak flows are higher in the urban rather than rural parts of the basin, and the percent of effective impervious area has a significant effect on storm runoff volumes and peaks. Water quality in the Chester Creek basin varies according to season and flow conditions. During low or base-flow conditions, concentrations of most water quality constituents measured are within State of Alaska drinking water standards, except for fecal coliform bacteria. During periods of high flow due to snowmelt or rainfall, concentrations of trace metal lead usually exceed recommended maximum levels. The primary sources of trace metal lead and suspended sediments are commercial areas , while the primary source of nutrients and fecal coliform bacteria is residential areas. Streamflow and water quality data collected at five sites representing different land-use categories were used to calibrate and verify three U.S. Geological Survey computer-based models: the Distributed Routing Rainfall-Runoff Model-Version II (DR3M-II), the Multi-Event Urban Runoff Quality Model (DR3M-QUAL), and the Precipitation Runoff Modeling Systems (PRMS). The PRMS can be used to simulate the effects of increased urbanization on daily flows. The DR3M-II can be used to simulate storm effects on small basins of < 40 acres. The DR3M-QUAL can be used to estimate seasonal loads of suspended sediment from basins of < 40 acres. (Lantz-PTT)

Alaska

Hydrogeology and ground-water flow in the carbonate rocks of the Little Lehigh Creek basin, Lehigh County, Pennsylvania

The Little Lehigh Creek basin is underlain mainly by a complex assemblage of highly-deformed Cambrian and Ordovician carbonate rocks. The Leithsville Formation, Allentown Dolomite, Beekmantown Group, and Jacksonburg Limestone act as a single hydrologic unit. Ground water moves through fractures and other secondary openings and generally is under water-table conditions. Median annual ground-water discharge (base flow) to Little Lehigh Creek near Allentown (station 01451500) during 1946-86 was 12.97 inches or 82 percent of streamflow. Average annual recharge for 1975-83 was 21.75 inches. Groundwater and surface-water divides do not coincide in the basin. Ground-water underflow from the Little Lehigh Creek basin to the Cedar Creek basin in 1987 was 4 inches per year. A double-mass curve analysis of the relation of cumulative precipitation at Allentown to the flow of Schantz Spring for 1956-84 showed that cessation of quarry pumping and development of ground water for public supply in the Schantz Spring basin did not affect the flow of Schantz Spring. Ground-water flow in the Little Lehigh Creek basin was simulated using a finite-difference, two-dimensional computer model. The geologic units in the modeled area were simulated as a single water-table aquifer. The 134-squaremile area of carbonate rocks between the Lehigh River and Sacony Creek was modeled to include the natural hydrologic boundaries of the ground-water-flow system. The ground-water-flow model was calibrated under steady-state conditions using 1975-83 average recharge, evapotranspiration, and pumping rates. Each geologic unit was assigned a different hydraulic conductivity. Initial aquifer hydraulic conductivity was estimated from specific-capacity data. The average (1975-83) water budget for the Little Lehigh Creek basin was simulated. The simulated base flow from the carbonate rocks of the Little Lehigh Creek basin above gaging station 01451500 is 11.85 inches per year. The simulated ground-water underflow from the Little Lehigh Creek basin to the Cedar Creek basin is 4.04 inches per year. For steady-state calibration, the root-mean-squared difference between observed and simulated heads was 21.19 feet. The effects of increased ground-water development on base flow and underflow out of the Little Lehigh Creek basin for average and drought conditions were simulated by locating a hypothetical well field in different parts of the basin. Steady-state simulations were used to represent equilibrium conditions, which would be the maximum expected long-term effect. Increased ground-water development was simulated as hypothetical well fields pumping at the rate of 15, 25, and 45 million gallons per day in addition to existing ground-water withdrawals. Four hypothetical well fields were located near and away from Little Lehigh Creek in upstream and downstream areas. The effects of pumping a well field in different parts of the Little Lehigh Creek basin were compared. Pumping a well field located near the headwaters of Little Lehigh Creek and away from the stream would have greatest effect on inducing underflow from the Sacony Greek basin and the least effect on reducing base flow and underflow to the Ceda^r Creek basin. Pumping a well field located near the headwaters of Little Leh|igh Creek near the stream would have less impact on inducing underflow from|the Sacony Creek basin and a greater impact on reducing the base flow of Little Lehigh Creek because more of the pumpage would come from diverted base flow. Pumping a well field located in the downstream area of the Little Lehigh Creek basin away from the stream would have the greatest effect on the underflow to the Cedar Creek basin. Pumping a well field located in the downstream area of the Little Lehigh Creek basin near the stream would have the greatest effect on reducing the base flow of Little Lehigh Cteek. Model simulations show that groundwater withdrawals do not cause a proportional reduction in base flow. Under average conditions, ground-water withdrawals are equal to 48 to 70 percent of simulated base-flow reductions; under drought conditions, ground-water withdrawals are equal to 35 to 73 percent of simulated base-flow reductions. The hydraulic effects of pumping largely depend on well location. In the Little Lehigh basin, surface-water and ground-water divides do not coincide, and ground-water development, especially near surface-water divides, can cause ground-water divides to shift and induce ground-water underflow from adjacent basins. Large-scale ground-water pumping in a basin may not produce expected reductions of base flow in that basin because of shifts in the ground-water divide; however, such shifts can reduce base flow in adjacent surface-water basins.

Water-Resources Investigations Report