Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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

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

At least 1,567 records · Page 87Linked to original sources

Land Disturbance Associated with Oil and Gas Development and Effects of Development-Related Land Disturbance on Dissolved-Solids Loads in Streams in the Upper Colorado River Basin, 1991, 2007, and 2025

Oil and gas resource development in the Upper Colorado River Basin (UCRB) has increased substantially since the year 2000. The UCRB encompasses several significant oil and gas producing areas that have the potential for continued oil and gas resource development. Land disturbance associated with oil and gas resource development is caused by activities related to constructing drill pads to contain drilling and well maintenance equipment and roads to access the drill pad. Land disturbed by oil and gas development has the potential to cause increased erosion, stream degradation, habitat fragmentation and alteration, and increase public use of areas that may be environmentally sensitive. Land disturbance resulting from oil and gas resource development has not been monitored and mapped on a regional scale in the UCRB. However, information on the location and age of oil and gas wells in the UCRB is available. These data combined with geographic data analysis and modeling techniques were used to estimate the total area of disturbed land associated with oil and gas resource development in 1991 and in 2007 in the UCRB. Additional information about anticipated oil and gas development in the UCRB was used to project land disturbance to the year 2025. Results of the analysis indicate that approximately 117,500 acres (183 mi2) of total land disturbance was associated with drill pads and related roads in the UCRB in 1991. The estimated area of disturbed land associated with oil and gas development increased 53 percent to 179,400 acres (280 mi2) in 2007. Projecting oil and gas development through 2025 results in a potential near doubling of the land surface disturbance to approximately 319,300 acres (500 mi2). Estimated land disturbance for 1991 and 2007 were input to a contaminant transport model developed for the UCRB to assess the statistical significance of energy-related land disturbance to contributing dissolved solids to basin streams. The statistical assessment was an observational study based on an existing model and available water-quality monitoring data for the basin. No new data were collected for the analysis. The source coefficient calibrated for the disturbed lands associated with oil and gas development in 2007 was zero, which indicated that estimated land disturbance from oil and gas development is not statistically significant in explaining dissolved solids in UCRB streams. The lack of significance in the contaminant transport modeling framework may be due to the amount of available monitoring data, the spatial distribution of monitoring sites with respect to land disturbance, or the overall quantity of land disturbance associated with oil and gas development basin wide. Finally, dissolved-solids loads derived from natural landscapes may be similar to loads derived from lands disturbed by oil and gas resource development. The model recalibration done for this study confirms calibration results from Kenney and others (2009): the most significant contributor to dissolved solids in the UCRB is irrigated agricultural land, which covers an area substantially larger than the estimated area disturbed by oil and gas development and is subjected to artificially applied water.

Scientific Investigations Report↗

Ground-Water System in the Chimacum Creek Basin and Surface Water/Ground Water Interaction in Chimacum and Tarboo Creeks and the Big and Little Quilcene Rivers, Eastern Jefferson County, Washington

A detailed study of the ground-water system in the unconsolidated glacial deposits in the Chimacum Creek Basin and the interactions between surface water and ground water in four main drainage basins was conducted in eastern Jefferson County, Washington. The study will assist local watershed planners in assessing the status of the water resources and the potential effects of ground-water development on surface-water systems. A new surficial geologic map of the Chimacum Creek Basin and a series of hydrogeologic sections were developed by incorporating LIDAR imagery, existing map sources, and drillers' logs from 110 inventoried wells. The hydrogeologic framework outlined in the study will help characterize the occurrence of ground water in the unconsolidated glacial deposits and how it interacts with the surface-water system. Water levels measured throughout the study show that the altitude of the water table parallels the surface topography and ranges from 0 to 400 feet above the North American Vertical Datum of 1988 across the basin, and seasonal variations in precipitation due to natural cycles generally are on the order of 2 to 3 feet. Synoptic stream-discharge measurements and instream mini-piezometers and piezometers with nested temperature sensors provided additional data to refine the positions of gaining and losing reaches and delineate seasonal variations. Chimacum Creek generally gains water from the shallow ground-water system, except near the community of Chimacum where localized losses occur. In the lower portions of Chimacum Creek, gaining conditions dominate in the summer when creek stages are low and ground-water levels are high, and losing conditions dominate in the winter when creek stages are high relative to ground-water levels. In the Quilcene Bay area, three drainage basins were studied specifically to assess surface water/ground water interactions. The upper reaches of Tarboo Creek generally gain water from the shallow ground-water system throughout most of the year and the lower reaches have little or no gains. The Big Quilcene River generally gains water from the shallow ground-water system after it emerges from a bedrock canyon and loses water from the town of Quilcene to the mouth of the river in Quilcene Bay. The Little Quilcene River generally loses water to the shallow ground-water system, although two localized areas were found to have gaining conditions. The Big Quilcene and Little Quilcene Rivers incur significant losses on the alluvial plain at the head of Quilcene Bay. Each of the creeks examined had a unique pattern of gaining and losing reaches, owing to the hydraulic conductivity of the streambed material and the relative altitude of the surrounding water table. Although the magnitudes of gains and losses varied seasonally, the spatial distribution did not vary greatly, suggesting that patterns of gains and losses in surface-water systems depend greatly on the geology underlying the streambed.

Scientific Investigations Report↗

Flood recovery maps for the White River in Bethel, Stockbridge, and Rochester, Vermont, and the Tweed River in Stockbridge and Pittsfield, Vermont, 2014

From August 28 to 29, 2011, Tropical Storm Irene delivered rainfall ranging from about 4 inches to more than 7 inches in the White River Basin. The rainfall resulted in severe flooding throughout the basin and significant damage along the White River and Tweed River. In response to the flooding, the U.S. Geological Survey, in cooperation with the Federal Emergency Management Agency, conducted a new flood study to aid in the flood recovery and restoration. This flood study includes a 20.7-mile reach of the White River from the downstream end at about 2,000 feet downstream from the State Route 107 bridge in the Village of Bethel, Vermont, to the upstream end at about 1,000 feet upstream from the River Brook Drive bridge in the Village of Rochester, Vt., and a 7.9-mile reach of the Tweed River from its mouth in Stockbridge, Vt., to the confluence of the West and South Branches of the Tweed River and continuing upstream on the South Branch Tweed River to the Pittsfield, Vt., town line. This report presents water-surface elevations determined for the study reaches using the U.S. Army Corps of Engineers one-dimensional step-backwater Hydrologic Engineering Center River Analysis System model, also known as HEC– RAS. The water-surface elevations were determined for floods having a 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probability (AEP) and for the floodway. Eighteen high-water marks from Tropical Storm Irene were available along the studied reaches. The discharges in the Tropical Storm Irene HEC–RAS model were adjusted so that the resulting water-surface elevations matched the high-water mark elevations along the study reaches. This allowed for an estimation of the water-surface profile throughout the study area resulting from Tropical Storm Irene. From a comparison of the estimated water-surface profile of Tropical Storm Irene to the water-surface profiles of the 1- and 0.2-percent AEP floods, it was determined that the high-water elevations resulting from Tropical Storm Irene exceeded the estimated 1-percent AEP flood throughout the White River and Tweed River study reaches and exceeded the estimated 0.2-percent AEP flood in 16.7 of the 28.6 study reach miles. The simulated water-surface profiles were then combined with a geographic information system digital elevation model derived from light detection and ranging (lidar) data having a 18.2-centimeter vertical accuracy at the 95-percent confidence level and 1-meter horizontal resolution to delineate the area flooded for each water-surface profile.

Vermont↗

Processes of Terrace Formation on the Piedmont of the Santa Cruz River Valley During Quaternary Time, Green Valley-Tubac Area, Southeastern Arizona

In this report we describe a series of stepped Quaternary terraces on some piedmont tributaries of the Santa Cruz River valley in southeastern Arizona. These terraces began to form in early Pleistocene time, after major basin-and-range faulting ceased, with lateral planation of basin fill and deposition of thin fans of alluvium. At the end of this cycle of erosion and deposition, tributaries of the Santa Cruz River began the process of dissection and terrace formation that continues to the present. Vertical cutting alternated with periods of equilibrium, during which streams cut laterally and left thin deposits of channel fill. The distribution of terraces was mapped and compiled with adjacent mapping to produce a regional picture of piedmont stream history in the middle part of the Santa Cruz River valley. For selected tributaries, the thickness of terrace fill was measured, particle size and lithology of gravel were determined, and sedimentary features were photographed and described. Mapping of terrace stratigraphy revealed that on two tributaries, Madera Canyon Wash and Montosa Canyon Wash, stream piracy has played an important role in piedmont landscape development. On two other tributaries, Cottonwood Canyon Wash and Josephine Canyon Wash, rapid downcutting preempted piracy. Two types of terraces are recognized: erosional and depositional. Gravel in thin erosional terraces has Trask sorting coefficients and sedimentary structures typical of streamflood deposits, replete with bar-and-swale surface topography on young terraces. Erosional-terrace fill represents the channel fill of the stream that cuts the terrace; the thickness of the fill indicates the depth of channel scour. In contrast to erosional terraces, depositional terraces show evidence of repeated deposition and net aggradation, as indicated by their thickness (as much as 20+ m) and weakly bedded structure. Depositional terraces are common below mountain-front canyon mouths where streams drop their load in response to abrupt flattening of gradients and expansion of channel banks, and they extend down the piedmont along Josephine Canyon Wash. Gravel in depositional terraces also has sorting coefficients typical of streamflood deposits. Sedimentary features in both types of terraces are consistent with deposition by flash floods in ephemeral streams, suggesting the climate was arid. Bedding and clast armor are weakly developed, clast clusters and imbrication are common, and crossbedding is generally absent. Debris-flow deposits, even near the mountain front, are surprisingly rare. On the tectonically stable piedmont of southeastern Arizona, stream piracy and climate change are the most likely agents of terrace formation. Both piracy and climate change can cause rapid changes in discharge and sediment supply, which initiate cycles of incision, lateral cutting, and aggradation. Increased stream discharge initiates downcutting, but increased sediment supply interrupts downcutting and causes streams to cut laterally and aggrade. At times, on Madera Canyon Wash and Montosa Canyon Wash, stream piracy affected stream discharge and sediment supply, but on Cottonwood Canyon Wash and Josephine Canyon Wash, only climate change could have initiated terrace cutting. Terraces probably formed during extended arid intervals when sparse vegetation and flashy stream discharge combined to increase sediment supply. In most cases, sediment supply was sufficient to promote lateral cutting but not long-term aggradation. Thus, most streams formed erosional terraces. The middle Pleistocene Josephine Canyon Wash formed a depositional terrace because it had a source of abundant unconsolidated sediment.

Scientific Investigations Report↗

Guide to Surficial Geology and River-Bluff Exposures, Noatak National Preserve, Northwestern Alaska

From its origin in rugged granitic highlands of the central Brooks Range, the Noatak River flows westward between the De Long Mountains and the Baird Mountains before turning south to enter Kotzebue Sound. Glaciers of middle and late Pleistocene age entered the Noatak River valley from the east, north, and south. Glaciers flowed down the upper Noatak River valley from the rugged peaks at its head, merging with tributary glaciers that issued from cirque-headed valleys along its south flank. Farther downvalley, small glaciers flowed northward from the Baird Mountains and much larger glaciers issued from the De Long Mountains. The De Long Mountains glaciers expanded southward to cover parts of the Noatak valley floor; they dammed the Noatak River during successive advances, creating a series of glacial lakes. The more extensive glacial advances dammed huge lakes that filled the Aniuk Lowland to overflowing. At various times, overflow waters spilled northward through Howard Pass, southward via Hunt River into the Kobuk River system, and westward down a series of channelways that skirted south of the glacier margins. Prominent bluffs along the Noatak River and its principal tributaries reveal glacial, glaciolacustrine, fluvial, and eolian sediments. More than 120 measured bluff exposures are described and illustrated in this report. These are dated by 92 radiocarbon age determinations and by the presence of the old Crow tephra, which was deposited about 130,000-140,000 years ago. Six geologic base maps, which cover sections of the Noatak River valley from east to west, show the locations of the river bluffs in relation to the glacial, glaciolacustrine, and fluvial deposits that cover the valley floor. The upper Noatak River valley is dominated by a bulky end moraine near Douglas Creek that was deposited during the last glacial maximum about 25,000-15,000 14C yr BP (termed the Itkillik II phase in the central Brooks Range glacial succession). Bluffs along this section of the Noatak River reveal thick till that underlies the moraine and interfingers downvalley with outwash and upvalley with moraine-dammed lake deposits. Remnants of older river gravels that underlie the set of glacial deposits contain wood fragments that are dated at about 35,000-30,000 14C yr BP. The Aniuk River area, which includes much of the eastern Aniuk Lowland, contains older moraines derived from headwaters of the Noatak valley that lie downvalley from the Douglas Creek moraine. These older moraines are assigned to the Itkillik IA and IB advances of the central Brooks Range glacial succession. Their deposits are seldom visible in river bluffs, but associated outwash and glaciolacustrine sediments are commonly exposed. More ancient end moraines farther downvalley are buried beneath lake deposits of the Aniuk Lowland, but are traceable as subdued arcuate drainage divides and as boulder concentrations in river bluffs or along their bases. The Cutler River area was occupied by glacial lakes assignable to three separate glacial phases. The oldest of these was probably dammed by the Cutler moraine, which crosses the Noatak valley floor near the mouth of Cutler River. The younger two are correlated with Itkillik-age deposits in the Aniuk Lowland. Glaciers of Itkillik age also flowed northwestward down the Cutler and Imelyak valley systems from cirques along the Noatak-Kobuk divide, but they did not reach the Noatak valley floor. The western Aniuk Lowland, which extends westward from the Cutler River mouth to the lower course of Nimiuktuk River, is dominated by a series of large end moraines deposited by glaciers from the De Long Mountains and that flowed southeastward down the Nimiuktuk valley system and then up the Noatak River valley. The Cutler moraine is the most extensive of these deposits. Following the Cutler glaciation, less extensive glacial advances built end moraines near the present-day mouths of Makpik Creek and Anisak River during in

Scientific Investigations Report↗

Potential effects of energy development on environmental resources of the Williston Basin in Montana, North Dakota, and South Dakota—Water resources

The Williston Basin has been a leading oil and gas producing area for more than 50 years. While oil production initially peaked within the Williston Basin in the mid-1980s, production rapidly increased in the mid-2000s, largely because of improved horizontal (directional) drilling and hydraulic fracturing methods. In 2012, energy development associated with the Bakken Formation was identified as a priority requiring collaboration toward improved timeliness of issuing permits for new wells combined with reasonable measures to maintain environmental quality. Shortly thereafter, the Bakken Federal Executive Group was created to address common challenges associated with energy development. The Bakken Federal Executive Group partner agencies identified a gap in current understanding of the cumulative environmental challenges attributed to energy development throughout the area, resulting in an effort to aggregate scientific data and identify additional research and information needs related to natural resources within areas of energy development in the Williston Basin. As part of this effort, water resources in the area (including groundwater; streams and rivers; and lakes, reservoirs, and wetlands) were characterized and described in terms of physical occurrence, flow characteristics, recharge, water quality, and water use. Similarly, waters produced during energy-development activities also were characterized even though these waters are not considered usable resources within the area. Groundwater resources were characterized by the major hydrogeologic units, or aquifers, identifying the units that supply most groundwater used for domestic, stock, agricultural, and industrial purposes. The groundwater characterization included other deeper hydrogeologic units in the Williston Basin that may be a useable source of water with treatment, have utility as a reservoir for reinjection of produced waters, or be a source of minerals and energy resources. A generalized groundwater budget and flow system identifying the sources of recharge (stream infiltration, precipitation, and movement [leakage] from other aquifers) and the general groundwater flow direction is included for each of the major hydrogeologic units. Rivers and streams within the Williston Basin with 10 or more years of continuous streamflow data were identified. For a subset of these sites, streamflow characteristics, including the monthly and annual mean flow, were generated to identify seasonal and interannual changes in streamflow and thus provide information on the drivers and reliability of streamflow at the seasonal or multiyear scale. Daily streamflow and annual extreme flows (peak and low flow) also were estimated for the subset of sites. The daily streamflow and annual extreme flow values provide information on short-term or extreme events that are relevant to infrastructure design and evaluating spills, leaks, or accidental discharges of water or petroleum products. Surface-water features (lakes, ponds, and wetlands) were classified using the Cowardin system and identified on the National Wetlands Inventory maps generated by the U.S. Fish and Wildlife Service. The spatial distribution of the surface-water features was analyzed by State, county, and specifically in comparison to the Prairie Pothole Region. The proximity of the surface-water features to energy development infrastructure (specifically oil or gas well pads) was evaluated. It was determined that, although oil or gas wells are often near a surface-water feature, most surface-water features do not have wells nearby, with the exception of wells in the Prairie Pothole Region. Water-quality data were aggregated from two data sources: (1) the Water-Quality Portal, sponsored by the U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (EPA), and National Water Quality Monitoring Council; and (2) a data compilation completed as part of the USGS National Water-Quality Assessment project. The Water-Quality Portal integrates publicly available water-quality data from databases maintained by the USGS, EPA, and U.S. Department of Agriculture, including water-quality data from Tribal, State, and local databases. Water-quality data for 15 commonly measured water-quality constituents were aggregated for groundwater, rivers and streams, and lakes and reservoirs. For each aggregated dataset (groundwater, rivers and streams, and lakes and reservoirs), analyses of the water-quality data included summary statistics, maps of spatial distribution of constituent values, boxplots of constituent values by timeframe or hydrogeologic unit, spatial comparisons of site locations and constituent values to petroleum well density, and comparisons of the constituent values measured to EPA drinking-water standards/guidelines. Produced water includes all fluids brought to the surface along with the targeted hydrocarbons as part of the oil and gas exploration and extraction processes. These fluids may include formation water (waters that co-exist with rock/oil/gas), hydraulic fracturing fluids, and other combinations of water and chemicals used during oil and gas well drilling, development, treatments, recompletions, and workovers. Produced water datasets were aggregated from two sources: the USGS National Produced Waters Geochemical database (ver. 2.1) and a series of projects focused specifically on sampling produced water in the Williston Basin from 2010 to 2014. The National Produced Waters Geochemical database was useful for a general understanding of produced-water chemistry. Produced waters are characterized by extreme salinity and contain elevated concentrations of other constituents (including arsenic, barium, cadmium, lead, zinc, radium-226/radium-228, and ammonium) that could negatively affect water and aquatic resources if released. Produced waters also have a generally unique chemical (isotopic) signature that may be useful in tracking water from different geologic units; for example, the oxygen/deuterium and strontium ratio values measured in brine waters from the Bakken Formation are distinct from brines collected from other geologic units in the Williston Basin. Water-use information related to energy production in the area also was aggregated and summarized. The summary of water use is not limited to oil and gas production but includes water used to produce all types of energy resources in the Williston Basin, including coal/lignite, thermoelectric power, oil and gas, hydropower, biomass and biofuels, wind, geothermal, and solar. Each State has its own methods for regulating and reporting water usage within its jurisdiction. These methods can introduce problems when examining water use from sources, such as the Missouri River or Fox Hills aquifer, that are shared across political boundaries. Without the one-to-one match for usage types and amounts used from a water source, it is difficult to develop a comprehensive water budget for the water source being evaluated. A large amount of freshwater is required to prepare a well for oil and gas well production; in some cases, 3 to 7 million gallons of water are needed per well. The EPA estimates that hydraulic fracturing in the Williston Basin uses between 70 to 140 billion gallons per year. Water also is used for myriad other purposes related to ancillary oil and gas extraction. In addition to water used for immediate energy development, the expanded human workforce migrating into the area and other support staff who have moved into the area during the development also use water. Research and information needs were identified that could be relevant in the evaluation of the effects of energy development on water resources. Information needs related to the evaluation of groundwater resources include the following: improved potentiometric-surface maps for glacial units; availability of a uniform stream network digital geographic coverage that spans the international boundary with Canada; enhanced surface-water use information with regards to the gain and loss of streamflow to shallow groundwater, which would increase understanding groundwater and surface-water interactions; and expanded geophysical assessments. Gaps in the availability of streamflow data include the lack of information on ice-jam flooding despite potential for effects to infrastructure (pipelines, roads, and facilities) and an understanding of the cumulative effects of largely undocumented stock and diversion dams. Although this study resulted in the aggregation of a large quantity of water-quality data, the availability of consistently collected, systematically processed and reported data over large parts of the Williston Basin is sparse. Few samples have been analyzed for constituents that may indicate the effect of energy development on water resources. Constituents that could be considered include boron, chloride, bromide, iodine, fluoride, manganese, lithium, radium, strontium isotopes, volatile organic compounds, and isotopes of inorganic ions (such as hydrogen and carbon). Collaboration between Tribal, Federal, State, and local entities to identify a common study design, common monitoring constituents, and consistent sampling locations would generate datasets with broad utility and would likely result in overall cost savings for monitoring over time. Similarly, there is a need for standardized sample collection, processing, laboratory analytical methods, and the collection of ancillary data for produced waters sampling. Additional characterization of the range of chemical, microbial, and isotopic compositions and quantities of “end-member” produced waters, and the collection of time-series datasets to document the changes in produced waters during and after well development also were needs identified during this study. Water-use estimates would be improved through the implementation of comprehensive studies of water use from groundwater and surface-water sources using consistent methodologies across the Williston Basin. The submission of chemical and water data related to hydraulic fracturing collected by the oil and gas industry would add to the quantity of available data. Consistent implementation of regulations and monitoring controls across political boundaries (State, county, and international) would further improve the consistency of data available for the estimates of water use.

Montana, North Dakota, South Dakota↗

Areal distribution and concentration of contaminants of concern in surficial streambed and lakebed sediments, Lake St. Clair and tributaries, Michigan, 1990-2003

As part of the Lake St. Clair Regional Monitoring Project, the U.S. Geological Survey evaluated data collected from surficial streambed and lakebed sediments in the Lake Erie-Lake St. Clair drainages. This study incorporates data collected from 1990 through 2003 and focuses primarily on the U.S. part of the Lake St. Clair Basin, including Lake St. Clair, the St. Clair River, and tributaries to Lake St. Clair. Comparable data from the Canadian part of the study area are included where available. The data are compiled into 4 chemical classes and consist of 21 compounds. The data are compared to effects-based sediment-quality guidelines, where the Threshold Effect Level and Lowest Effect Level represent concentrations below which adverse effects on biota are not expected and the Probable Effect Level and Severe Effect Level represent concentrations above which adverse effects on biota are expected to be frequent. Maps in the report show the spatial distribution of the sampling locations and illustrate the concentrations relative to the selected sediment-quality guidelines. These maps indicate that sediment samples from certain areas routinely had contaminant concentrations greater than the Threshold Effect Concentration or Lowest Effect Level. These locations are the upper reach of the St. Clair River, the main stem and mouth of the Clinton River, Big Beaver Creek, Red Run, and Paint Creek. Maps also indicated areas that routinely contained sediment contaminant concentrations that were greater than the Probable Effect Concentration or Severe Effect Level. These locations include the upper reach of the St. Clair River, the main stem and mouth of the Clinton River, Red Run, within direct tributaries along Lake St. Clair and in marinas within the lake, and within the Clinton River headwaters in Oakland County. Although most samples collected within Lake St. Clair were from sites adjacent to the mouths of its tributaries, samples analyzed for trace-element concentrations were collected throughout the lake. The distribution of trace-element concentrations corresponded well with the results of a two-dimensional hydrodynamic model of flow patterns from the Clinton River into Lake St. Clair. The model was developed independent from the bed sediment analysis described in this report; yet it showed a zone of deposition for outflow from the Clinton River into Lake St. Clair that corresponded well with the spatial distribution of trace-element concentrations. This zone runs along the western shoreline of Lake St. Clair from L'Anse Creuse Bay to St. Clair Shores, Michigan and is reflected in the samples analyzed for mercury and cadmium. Statistical summaries of the concentration data are presented for most contaminants, and selected statistics are compared to effects-based sediment-quality guidelines. Summaries were not computed for dieldrin, chlordane, hexachlorocyclohexane, lindane, and mirex because insufficient data are available for these contaminants. A statistical comparison showed that the median concentration for hexachlorobenzene, anthracene, benz[ a ]anthracene, chrysene, and pyrene are greater than the Threshold Effect Concentration or Lowest Effect Level. Probable Effect Concentration Quotients provide a mechanism for comparing the concentrations of contaminant mixtures against effects-based biota data. Probable Effect Concentration Quotients were calculated for individual samples and compared to effects-based toxicity ranges. The toxicity-range categories used in this study were nontoxic (quotients < 0.5) and toxic (quotients > 0.5). Of the 546 individual samples for which Probable Effect Concentration Quotients were calculated, 469 (86 percent) were categorized as being nontoxic and 77 (14 percent) were categorized as being toxic. Bed-sediment samples with toxic Probable Effect Concentration Quotients were collected from Paint Creek, Galloway Creek, the main stem of the Clinton River, Big Beaver Creek, Red Run, Clinton River towards the mouth, Lake St. Clair along the western shore, and the St. Clair River near Sarnia.

Michigan↗

Distribution and sources of polychlorinated biphenyls in Woods Inlet, Lake Worth, Fort Worth, Texas, 2003

Woods Inlet is a flooded stream channel on the southern shore of Lake Worth along the western boundary of Air Force Plant 4 in Fort Worth, Texas, where elevated polychlorinated biphenyl (PCB) concentrations in sediment were detected in a previous study. In response, the U.S. Geological Survey, in cooperation with the U.S. Air Force, conducted a study in 2003 to map the extent of elevated PCB concentrations in Woods Inlet and to identify possible sources (or more specifically, source areas) of PCBs in the watershed of Woods Inlet. Three gravity cores (penetration to pre-reservoir sediment at three sites) and 17 box cores (surficial bottom sediment samples) were collected in Woods Inlet. Suspended sediment in stormwater runoff and streambed sediment were sampled in tributaries to Woods Inlet following storms. Assemblages of PCB congeners in surficial inlet sediments and suspended and streambed sediments were analyzed to indicate sources of PCBs in the inlet sediments on the basis of chemical signatures of PCBs. Woods Inlet receives runoff primarily from three tributaries: (1) Gruggs Park Creek, (2) the small unnamed creek that drains a Texas National Guard maintenance facility, called TNG Creek for this report, and (3) Meandering Road Creek. Twenty-seven of 209 possible PCB congeners were analyzed. The sum of the congeners was used as a measure of total PCB. The spatial distribution of total PCB concentrations in the inlet indicates that most PCBs are originating in the Meandering Road Creek watershed. Peak total PCB concentrations in the three gravity cores occurred at depths corresponding to sediment deposition dates of about 1960 for two of the cores and about 1980 for the third core. The magnitudes of peak total PCB concentrations in the gravity cores followed a spatial distribution generally similar to that of surficial bottom sediment concentrations. Total PCB concentrations in suspended and streambed sediment varied greatly between sites and indicated a likely source of PCBs associated with a sampling site that receives runoff from Air Force Plant 4. Three approaches to the analyses of congener assemblages indicate that PCBs in surficial bottom sediment of Woods Inlet primarily enter Lake Worth from Meandering Road Creek and that runoff from Air Force Plant 4 is a source of the PCBs in Meandering Road Creek. Although current (2003) transport of PCBs from Air Force Plant 4 to the creek is occurring, large decreases in PCB concentrations with decreasing age in two cores indicate that PCB loading to the inlet has decreased greatly since the 1960s. Because runoff entering Meandering Road Creek from some parts of Air Force Plant 4 was not measured or sampled in this study, it cannot be said with certainty that the Air Force Plant 4 site sampled is the only source of PCBs to Meandering Road Creek.

Texas↗

Geology, ground-water hydrology, geochemistry, and ground-water simulation of the Beaumont and Banning Storage Units, San Gorgonio Pass area, Riverside County, California

Ground water has been the only source of potable water supply for residential, industrial, and agricultural users in the Beaumont and Banning storage units of the San Gorgonio Pass area, Riverside County, California. Ground-water levels in the Beaumont area have declined as much as 100 feet between the early 1920s and early 2000s, and numerous natural springs have stopped flowing. In 1961, the San Gorgonio Pass Water Agency (SGPWA) entered into a contract with the California State Department of Water Resources to receive 17,300 acre-feet per year of water to be delivered by the California State Water Project (SWP) to supplement natural recharge. Currently (2005), a pipeline is delivering SWP water into the area, and the SGPWA is artificially recharging the ground-water system using recharge ponds located along Little San Gorgonio Creek in Cherry Valley with the SWP water. In addition to artificial recharge, SGPWA is considering the direct delivery of SWP water for the irrigation of local golf courses and for agricultural supply in lieu of ground-water pumpage. To better understand the potential hydrologic effects of different water-management alternatives on ground-water levels and movement in the Beaumont and Banning storage units, existing geohydrologic and geochemical data were compiled, new data from a basin-wide ground-water level and water-quality monitoring network were collected, monitoring wells were installed near the Little San Gorgonio Creek recharge ponds, geohydrologic and geochemical analyses were completed, and a ground-water flow simulation model was developed. The San Gorgonio Pass area was divided into several storage units on the basis of mapped or inferred faults. This study addresses primarily the Beaumont and Banning storage units. The geologic units in the study area were generalized into crystalline basement rocks and sedimentary deposits. The younger sedimentary deposits and the surficial deposits are the main water-bearing deposits in the San Gorgonio Pass area. The water-bearing deposits were divided into three aquifers: (1) the perched aquifer, (2) the upper aquifer, and (3) the lower aquifer based on lithologic and downhole geophysical logs. Natural recharge in the San Gorgonio Pass area was estimated using INFILv3, a deterministic distributed- parameter precipitation-runoff model. The INFILv3 model simulated that the potential recharge of precipitation and runoff in the Beaumont and Banning storage units was about 3,710 acre-feet per year and that the potential recharge in 28 sub-drainage basins upstream of the storage units was about 6,180 acre-feet per year. The water supply for the Beaumont and Banning storage units is supplied by pumping ground water from wells in the Canyon (Edgar and Banning Canyons), Banning Bench, Beaumont, and Banning storage units. Total annual pumpage from the Beaumont and Banning storage units ranged from about 1,630 acre-feet in 1936 to about 20,000 acre-feet in 2003. Ground-water levels declined by as much as 100 feet in the Beaumont storage unit from 1926-2003 in response to ground-water pumping of about 450,160 acre-feet during this period. Since ground-water development began in the San Gorgonio Pass area, there have been several sources of artificial recharge to the basin including return flow from applied water on crops, golf courses, and landscape; septic-tank seepage; and infiltration of storm runoff diversions and imported water into recharge ponds. Return flow from applied water and septic-tank seepage was estimated to reach a maximum of about 8,100 acre-feet per year in 2003. Owing to the great depth of water in much of study area (in excess of 150 feet), the return flow and septic-tank seepage takes years to decades to reach the water table. Stable-isotope data indicate that the source of ground-water recharge was precipitation from storms passing through the San Gorgonio Pass as opposed to runoff from the higher altitudes of the San Bernar

Scientific Investigations Report↗

Death Valley regional ground-water flow system, Nevada and California -- hydrogeologic framework and transient ground-water flow model

A numerical three-dimensional (3D) transient ground-water flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the ground-water flow system and previous less extensive ground-water flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect ground-water flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley Regional Ground-Water Flow System (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the ground-water flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural ground-water discharge occurring through evapotranspiration and spring flow; the history of ground-water pumping from 1913 through 1998; ground-water recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional ground-water flow system. Ground-water flow in the Death Valley region is composed of several interconnected, complex ground-water flow systems. Ground-water flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional ground-water flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and evapotranspiration (ET) are the dominant natural ground-water discharge processes. Ground water also is withdrawn for agricultural, commercial, and domestic uses. Ground-water flow in the DVRFS was simulated using MODFLOW-2000, a 3D finite-difference modular ground-water flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 m on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins, also simulated by drains. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model gener

Scientific Investigations Report↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River near Kansas City, Missouri, August 8–9, 2023

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 8 bridge crossings of the Missouri River near Kansas City, Missouri, on August 8–9, 2023. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches that extended about 1,550 to 1,640 feet longitudinally and generally extended laterally across the active channel from bank to bank during low floodflow to nonflood conditions. These surveys provided the channel geometry and hydraulic conditions of the river at the time of the surveys and provided characteristics of scour holes, which may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. The data collected from the surveys may also be useful to the Missouri Department of Transportation as a record of low floodflow conditions in regards to the stability and integrity of the bridges with respect to bridge scour. Bathymetric data were collected around every in-channel pier. Scour holes were at most piers where bathymetry could be obtained, except for those piers on banks or surrounded by riprap. All the bridge sites in this study were surveyed and documented in previous studies. The average difference between the bathymetric surfaces ranged from 0.07 to 4.16 feet higher in 2023 than 2019, which indicates overall deposition between the survey dates, as might be expected based purely on streamflow at the time of the survey. However, the average difference between the bathymetric surfaces ranged from 1.44 feet higher to 1.88 feet lower in 2023 than 2015, which indicates a dynamic equilibrium of scour and deposition overall between those surveys, despite the lower flow conditions in 2023. Similarly, the average difference between the bathymetric surfaces ranged from 3.18 feet higher to 5.19 feet lower in 2023 than 2011, which indicates a relative equilibrium between scour and deposition overall, albeit the trend was toward scour as might be expected because of the substantial flood event in 2011. Riprap blankets and alignment to flow had a substantial effect on the size of the scour hole for a given pier. Piers that were partially or fully surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers with no rock or riprap and effectively mitigated the scour holes historically observed at these piers. Several of the structures had piers that were skewed to primary approach flow. At most of the structures, the scour hole was deeper and longer on the side of the pier with impinging flow than the leeward side, with some amount of deposition on the leeward side, as typically observed at piers skewed to approach flow.

Kansas, Missouri↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, April-May, 2013

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, in the vicinity of 10 bridges at 9 highway crossings of the Missouri River between Lexington and Washington, Missouri, from April 22 through May 2, 2013. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches ranging from 1,640 to 1,840 feet longitudinally and extending laterally across the active channel between banks and spur dikes in the Missouri River during low- to moderate-flow conditions. These bathymetric surveys indicate the channel conditions at the time of the surveys and provide characteristics of scour holes that may be useful in the development of predictive guidelines or equations for scour holes. These data also may be useful to the Missouri Department of Transportation to assess the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every pier that was in water, except those at the edge of water or in very shallow water (less than about 6 feet). Scour holes were present at most piers for which bathymetry could be obtained, except at piers on channel banks, near or embedded in lateral or longitudinal spur dikes, and on exposed bedrock outcrops. Scour holes observed at the surveyed bridges were examined with respect to depth and shape. Although exposure of parts of foundational support elements was observed at several piers, at most sites the exposure likely can be considered minimal compared to the overall substructure that remains buried in channel-bed material; however, there were several notable exceptions where the bed material thickness between the bottom of the scour hole and bedrock was less than 6 feet. Such substantial exposure of usually buried substructural elements may warrant special observation in future flood events. Previous bathymetric surveys had been done at all of the sites in this study during the flood of 2011. Comparisons between bathymetric surfaces from the previous surveys and those of this study generally indicate a consistent increase in the elevation of the bed and decrease in the size of scour holes at these sites, both likely caused by a substantial decrease in discharge and water-surface elevation compared to the 2011 surveys at most sites. However, multiple surveys at one of the sites indicate that the flow condition is not the sole variable in the determination of the size of scour holes at sites with a dual bridge configuration. Furthermore, another site had a smaller and shallower scour hole even though the discharge in this study was slightly greater than in 2011. Pier size, nose shape, and alignment to flow also had a substantial effect on the size of the scour hole observed.

Missouri↗

Potentiometric Surfaces and Water-Level Trends in the Cockfield and Wilcox Aquifers of Southern and Northeastern Arkansas, 2006

The Cockfield Formation of Claiborne Group and the Wilcox Group contain aquifers that provide sources of ground water in southern and northeastern Arkansas. In 2000, about 9.9 million gallons per day was withdrawn from the Cockfield Formation of Claiborne Group and about 22.2 million gallons per day was withdrawn from the Wilcox Group. Major withdrawals from the aquifers were for industrial and public water supplies. A study was conducted by the U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission and the Arkansas Geological Survey to determine the water level associated with the aquifers in the Cockfield Formation of Claiborne Group and the Wilcox Group in southern and northeastern Arkansas. During February and March 2006, 56 water-level measurements were made in wells completed in the Cockfield aquifer and 59 water-level measurements were made in wells completed in the Wilcox aquifer, 16 in southwestern and 43 in northeastern Arkansas. This report presents the results as potentiometric-surface maps and as long-term water-level hydrographs. The regional direction of ground-water flow in the Cockfield Formation of Claiborne Group generally is towards the east and southeast, away from the outcrop, except in areas of intense ground-water withdrawals, such as western Drew County, southeastern Lincoln County, southwestern Calhoun County, and near Crossett in Ashley County. There are three cones of depression indicated by relatively low water-level altitudes in southeastern Lincoln County, southwestern Calhoun County, and near Crossett in Ashley County. The lowest water-level altitude measured was 44 feet above the National Geodetic Vertical Datum of 1929 in Lincoln County; the highest water-level altitude measured was 346 feet above the National Geodetic Vertical Datum of 1929 in Columbia County at the outcrop area. Hydrographs from 40 wells with historical water levels from 1986 to 2006 were evaluated using linear regression to calculate the annual rise or decline. Calhoun and Cleveland Counties have mean annual rises from 0.01 to 0.07 feet per year. Arkansas, Ashley, Bradley, Chicot, Columbia, Drew, Lincoln, and Union Counties have mean annual declines from 0.4 to 0.55 feet per year. Desha County has a mean annual decline of about 1.35 feet per year. The direction of ground-water flow in the southwestern study area of the Wilcox Group generally is south and east. The lowest water-level altitude measured in southwestern Arkansas was 147 feet above the National Geodetic Vertical Datum of 1929 near the Ouachita River in Clark County; the highest water-level altitude measured was 397 feet above the National Geodetic Vertical Datum of 1929 in the outcrop area of Hempstead County. The direction of ground-water flow in the northeastern study area of the Wilcox Group generally is south and east. The lowest water-level altitude measured in northeastern Arkansas was 120 feet above the National Geodetic Vertical Datum of 1929 near West Memphis in Crittenden County; the highest water-level altitude measured was 368 feet above the National Geodetic Vertical Datum of 1929 on Crowleys Ridge in Clay County. Hydrographs from 28 wells with historical water levels from 1986 to 2006 were evaluated using linear regression to calculate the annual rise or decline. All 28 wells showed an annual decline from 1986 to 2006. Craighead, Greene, Mississippi, and Poinsett Counties have mean annual declines from 0.27 to 1.00 feet per year. Crittenden, Lee, and St. Francis Counties have mean annual declines from 1.39 to 1.64 feet per year.

Scientific Investigations Report↗

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is &ldquo;subhorizontal&rdquo; in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗

Preliminary geohydrologic assessment of Buenos Aires National Wildlife Refuge, Altar Valley, southeastern Arizona

The Buenos Aires National Wildlife Refuge is located in the southern part of Altar Valley, southwest of Tucson in southeastern Arizona. The primary water-supply well at the Buenos Aires National Wildlife Refuge has experienced a two-decade decrease in groundwater levels in the well, as have other wells in the southern part of Altar Valley. In part to understand this trend, a study was undertaken by the U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, to summarize what is known about the geohydrologic system on the refuge and analyze groundwater-level trends and precipitation-groundwater correlations. In addition, available data were compiled where possible on the climate, land cover, soils, geology, and hydrology to provide a foundation for future modeling of the system. Altar Valley is a sedimentary basin bounded by a mixture of Paleozoic to Tertiary sedimentary, volcanic, granitic, and metamorphic rocks. The valley fill is undifferentiated Tertiary to Quaternary sediments underlain by middle Miocene to Pliocene rocks that consist of moderately to strongly consolidated conglomerate and sandstone. Surface water, when present in the predominantly ephemeral streams of the valley, flows from south to north. Arivaca Creek has a cienega (or wetland) where groundwater surfaces before it flows as a short perennial reach out of Arivaca Basin. Groundwater maps compiled between 1934 and 2016 showed groundwater flowing from south to north. Before the 1980s, temporal patterns of groundwater levels in wells in Altar Valley varied substantially from one well to another. In the mid-1980s, comparatively high levels of precipitation occurred: the 1980s median value was 15.3 inches, whereas the median for the period of record was 13.2 inches. In addition, apparently corresponding groundwater level increases were seen in nearly all wells studied. After this initial increase, two different groundwater-level trends began to be observed in two spatially distinct sets of wells: in the northern part, groundwater levels were relatively steady, whereas in the southern part, groundwater levels declined from 10 to 20 feet between 1990 and 2019. Annual groundwater pumpage declined substantially in the northern part of the valley beginning in the early 1980s, but it began to increase again in the 1990s. Pumpage in the southern part has remained low and relatively steady compared to the northern part. Although the precise reasons for the declining groundwater levels in the southern part remain unclear, groundwater levels may be affected by factors such as climate cycles, long-term drought, and temperature-induced declines in recharge, resulting in increased evapotranspiration. Preliminary analyses of two wells, one selected from each part of the valley, using linear regression and lag correlation to investigate correlation between annual precipitation and groundwater levels, showed a maximum correlation at a lag of about 17 years in the southern part of the valley and about 25 years in the northern part, indicating that, although variable sources and traveltimes of recharged water may be needed to propagate to each location, the strongest correlation at each well is with precipitation that was recharged 17 and 25 years prior to the groundwater response in that well. Assuming a constant flow of groundwater from the southern to the northern part of the valley, a decrease in recharge is expected to lead to a decrease in aquifer storage. As to the comparatively stable groundwater levels in the northern part, pumpage is still only about one-half what it was in the early 1980s, even though pumpage has increased there since the 1990s. Water levels in most wells in the northern part were drawn down prior to the decrease in pumping in the early 1980s, possibly owing to a combination of pumping and the nearly 20-year midcentury drought that occurred between 1940 and 1960. Water levels were in the process of recovering when the increase in pumping occurred in the 1990s. Because the water levels were recovering (increasing) instead of remaining static, the increased pumping may have only limited the recovery rather than causing a decrease in water levels, as a new quasi-equilibrium state may have been reached. Additional possible causes for the stable groundwater levels include (1) upgradient aquifer transmissivity that was high enough to offset pumping, (2) a low-permeability barrier, such as bedrock or clay, at the north end of the valley that caused groundwater pooling, (3) higher lateral inflow of groundwater in the northern part of the valley, (4) a delay in the effect of storage declines propagating from the south, or (5) some combination thereof.

Arizona↗

Porphyry copper assessment of East and Southeast Asia: Philippines, Taiwan (Republic of China), Republic of Korea (South Korea), and Japan

The U.S. Geological Survey collaborated with member countries of the Coordinating Committee for Geoscience Programmes in East and Southeast Asia (CCOP) on an assessment of the porphyry copper resources of East and Southeast Asia as part of a global mineral resource assessment. The assessment covers the Philippines in Southeast Asia, and the Republic of Korea (South Korea), Taiwan (Province of China), and Japan in East Asia. The Philippines host world class porphyry copper deposits, such as the Tampakan and Atlas deposits. No porphyry copper deposits have been discovered in the Republic of Korea (South Korea), Taiwan (Province of China), or Japan. Thirteen geographic areas were delineated as tracts that are permissive for porphyry copper deposits in the assessed area. Individual tracts range from about 3,000 to 100,000 square kilometers in area. Permissive tracts are delineated on the basis of mapped distributions of igneous rocks of specific age ranges that define subduction-related magmatic arcs or magmatic belts that might contain porphyry copper deposits. Most of these magmatic arcs are subduction related, although some porphyry deposits and prospects are present in back-arc or poorly understood tectonic settings. Maps at various scales were used in the compilation; however, the final tract boundaries are intended for use at a scale of 1:1,000,000. Numbers of undiscovered deposits were estimated at different levels of confidence for 10 permissive tracts in the Philippines including one area that extends to eastern Taiwan (Republic of China); permissive tracts in South Korea and Japan are discussed qualitatively. Estimates of numbers of undiscovered deposits were combined with grade and tonnage models using Monte Carlo simulation to estimate amounts of undiscovered resources. Grades and tonnages of known porphyry copper deposits in the study area were compared with global grade and tonnage models to determine the appropriate model for simulation of undiscovered resources. Most of the known deposits are best described as copper-gold subtypes of porphyry copper deposits. For some permissive tracts, a general porphyry copper-gold-molybdenum model was used. Thirty-eight porphyry copper deposits are known in the Philippines; the mean number of undiscovered deposits was estimated to be 28. Mean (arithmetic) resources that could be associated with the undiscovered deposits are 90 million metric tons of copper and 5,800 metric tons of gold, as well as byproduct molybdenum and silver. Additional resources that could be discovered in extensions to known deposits were not evaluated. Assessment results, presented in tables and graphs, indicate expected amounts of total contained metal and mineralized rock in undiscovered deposits at different quantile levels, as well as the arithmetic mean for each tract. The Philippines have a long history of porphyry exploration cycles and mine development, interrupted at times by political and social unrest, environmental concerns, and natural disasters. Changes in mining laws within the region and the recent high price of gold on the world market have prompted renewed interest in porphyry copper deposits in the region. South Korea and Japan have been thoroughly explored for many types of mineral deposits. Available data suggest that the permissive rocks in South Korea typically are too deeply eroded to preserve porphyry copper deposits. Porphyry copper systems may be present in Japan, but are likely to lie at depths greater than the 1 kilometer from the surface protocol adopted for this study. Descriptions of the geologic basis for delineating each tract, the data used, the geologic criteria and rationale for the assessment, and results of the assessment are included in appendixes along with the description of a geographic information system (GIS) that includes tract boundaries, known porphyry copper deposits and significant prospects, and assessment results.

Scientific Investigations Report↗

Hydrogeology and simulation of groundwater flow in fractured rock in the Newark basin, Rockland County, New York

Groundwater in the Newark basin aquifer flows primarily through discrete water-bearing zones parallel to the strike and dip of bedding, whereas flow perpendicular to the strike is restricted, thereby imparting anisotropy to the groundwater flow field. The finite-element model SUTRA was used to represent bedrock structure in the aquifer by spatially varying the orientation of the hydraulic conductivity tensor to reflect variations in the strike and dip of the bedding. Directions of maximum and medium hydraulic conductivity were oriented parallel to the bedding, and the direction of minimum hydraulic conductivity was oriented perpendicular to the bedding. Groundwater flow models were prepared to simulate local flow in the vicinity of the Spring Valley well field and regional flow through the Newark basin aquifer. The Newark basin contains sedimentary rocks deposited as alluvium during the Late Triassic and is one of a series of basins that developed when Mesozoic rifting of the super continent Pangea created the Atlantic Ocean. The westward-dipping basin is filled with interbedded facies of coarse-grained to fine-grained rocks that were intruded by diabase associated with Jurassic volcanism. The Newark basin aquifer is bounded to the north and east by the Palisades sill and to the west by the Ramapo Fault. Although the general dip of bedding is toward the fault, mapping of conglomerate beds indicates the rocks are folded into broad anticlines and synclines. An alternative, more uniform pattern of regional structure, based on interpolated strike and dip measurements from a number of sources, has also been proposed. Two groundwater flow models (A for the former type of bedrock structure and B for the latter type) were developed to represent these alternative depictions of bedrock structure. Transient simulations were calibrated to reproduce measured water-level recoveries in a 9.3 mi&sup2 area surrounding the Spring Valley well field during a 5-day aquifer test in 1992. The models represented a 330-ft thick rock mass divided vertically into 10 equally spaced layers and were calibrated through nonlinear regression. Results of model B best matched the observed water-level recoveries with an estimated hydraulic conductivity of 9.5 ft/day, specific storage of 7.6 x 10 -6 ft -1, and K max : K min anisotropy ratio (hydraulic conductivity parallel to bedding: perpendicular to bedding) of 72:1. Model error was 50 percent greater in model A because the assumed structure did not match the actual strike of bedding in this area. Steady-state simulations of regional flow through the 85.4-mi2 modeled extent of the Newark basin aquifer represented both the alluvial aquifer beneath the Mawah River and the fractured bedrock. The rock mass was divided into two aquifer units: an upper 500-ft thick unit divided into 10 equally spaced layers through which most groundwater is assumed to flow and a lower unit divided into 7 layers with increasing thickness. Models were calibrated through nonlinear regression to average water levels measured in 140 wells from August 2005 through April 2007. Water levels simulated using the two models were similar and generally matched those observed, and the average recharge rate estimated using both models was 19 inches/year for the simulated period. Estimated transmissivity parallel to the strike of bedding (1,100 ft&sup2/d) was uniform in two transmissivity (T) zones in model A, but in model B the transmissivity of a high T zone (1,600 ft&sup2/d), delineated on the basis of aquifer test data, was slightly greater than in a low T zone (1,300 ft&sup2/d). The K max : K min anisotropy was estimated to be 58:1 in model A and 410:1 in model B, so the proportion of flow perpendicular to bedding is less in model B than in model A. Distributions of groundwater age simulated with models A and B are similar and indicate that most shallow ground-water (225 ft below the bedrock surface) is 5 t

Scientific Investigations Report↗

Significance of the precambrian basement and late Cretaceous thrust nappes on the location of tertiary ore deposits in the Oquirrh Mountains, Utah

The Oquirrh Mountains are located in north central Utah, in the easternmost part of the Basin and Range physiographic province, immediately south of the Great Salt Lake. The range consists of a northerly trending alignment of peaks 56 km long. Tooele and Rush Valleys flank the Oquirrh Mountains on the western side and Salt Lake and Cedar Valleys lie on the eastern side. The world class Bingham mine in the central part of the range hosts disseminated copper-bearing porphyry, skarn, base-and precious-metal vein and replacement ore deposits. The district includes the outlying Barneys Canyon disseminated-gold deposits. Disseminated gold in the Mercur mining district in the southern part of the range has become exhausted. The Ophir and Stockton base- and precious-metal mining districts in the range north of Mercur also are inactive. A geologic map of the range (Tooker and Roberts, 1998), available at a scale of 1:50,000, is a summation of U.S. Geological Survey (USGS) studies. Information about the range and its mining areas is scattered. This report summarizes map locations, new stratigraphic and structural data, and reexamined data from an extensive published record. Unresolved controversial geological interpretations are considered, and, for the first time, the complete geological evidence provides a consistent regional basis for the location of the ore deposits in the range. The geological setting and the siting of mineral deposits in the Oquirrh Mountains began with the formation of a Precambrian craton. Exposures of folded Proterozoic basement rocks of the craton, in the Wasatch Mountains east of Salt Lake City, were accreted and folded onto an Archean crystalline rock terrane. The accretion suture lies along the north flank of the Uinta Mountains. The western part of the accreted block was offset to northern Utah along a north-trending fault lying approximately along the Wasatch Front (Nelson and others, 2002), thereby creating a prominant basement barrier or buttress east of the Salt Lake area. The accretion suture along the north flank of the Uinta Anticline overlaps an earlier Precambrian east-west mobile zone, the Uinta trend (Erickson, 1976, Bryant and Nichols, 1988 and John, 1989), which extends westward across western Utah and into Nevada. A trace of the trend underlies the middle part of the Oquirrh Mountains. Its structure is recognized by disrupted Paleozoic stratigraphic units and fold and fault evidence of thrust faulting, intermittent local uplift and erosion, the alignment of Tertiary intrusives and associated ore deposits. Geologic readjustments along the trend continued intermittently through the Paleozoic, Cenozoic, Tertiary, and the development of clastic deposits along the shores of Pleistocene Lake Bonneville. Paleozoic sedimentary rocks were deposited on the craton platform shelf in westernmost Utah and eastern Nevada as the shelf subsided gradually and differentially. Debris was shed into two basins separated by the uplifted Uinta trend, the Oquirrh Basin on the south and Sublette Basin on the north. Sediments were derived from the craton to the east, the Antler orogenic zone on the west (Roberts, 1964), and locally from uplifted parts of the trend itself. Thick accumulations of clastic calcareous quartzite, shale, limestone, and dolomite of Lower and Upper Paleozoic ages are now exposed in the Oquirrh Mountains, the result of thrust faults. Evidence of decollement thrust faults in in the Wasatch Mountains during the Late Cretaceous Sevier orogeny, recognized by Baker and others (1949) and Crittenden (1961, is also recognized in the Oquirrh Mountains by Roberts and others (1965). During the late Cretaceous Sevier Orogeny, nappes were thrust sequentially along different paths from their western hinterland to the foreland. Five distinct nappes converged over the Uinta trend onto an uplifted west-plunging basement buttress east of the Oquirrh Mountains area: the Pass Canyon, Bingham,

Scientific Investigations Report↗