USGS Science⌕ Search

SEARCH · USGS Science

Results for “Facies”

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,189 records · Page 66Linked to original sources
USGS2022

Bedrock geologic map of the Crown Point quadrangle, Essex County, New York, and Addison County, Vermont

The bedrock geology of the 7.5-minute Crown Point quadrangle consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of granitic orthogneiss at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the host rocks. Subsequent D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into previously deformed lithologies. The LMG crystallized at approximately 1.15 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation during D3 produced F3 folds exhibited in regional domes and basins, such as the Keeney Mountain synform, local reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. The Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. Cambrian-Ordovician stratigraphy records an approximately 1-km-thick section and a transition from synrift clastics to passive margin peritidal carbonate buildups to gradually deeper water subtidal to shelf carbonates during foreland basin development associated with the Taconic orogeny. The Paleozoic rocks are weakly folded and block faulted. Large areas of the Champlain Valley are covered by undifferentiated glacial deposits, some of which contain mapped landslides. The map also shows waste rock piles and tailings from historical mining operations and large areas of artificial fill. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondacks, provide a context for historical iron mines in the eastern Adirondacks, and update the stratigraphy of the Champlain Valley in New York and Vermont. This Scientific Investigations Map of the Crown Point 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, three cross sections, and a simplified surficial geologic map that includes lidar percent slope. The explanatory pamphlet includes a discussion of the geology. The bedrock geologic map on the map sheet is multi-layered and has been designed to enable the user to turn off the surficial map layer to view the concealed bedrock map units.

New York, Vermont↗
USGS2026

Bedrock geologic map of the Eagle Lake quadrangle, Essex County, New York

The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.

New York↗
USGS2008

The Yampa Bed — A regionally extensive tonstein in the Williams Fork Formation, northwestern Piceance Creek and southern Sand Wash Basins, Colorado

A regionally persistent and distinctive unit of Upper Cretaceous age is here formally named the Yampa Bed of the Williams Fork Formation for exposures in the Yampa, Danforth Hills, and Grand Hogback coal fields, Moffat and Routt Counties, northwest Colorado; the name is derived from the Yampa River valley. The type section was measured in the NE¼ SW¼ sec. 6, T. 5 N., R. 91 W., about 8 miles south of Craig, Colo., where the bed is 38 inches thick and lies within the C-D coal bed in the lower part of the Williams Fork Formation, about 165 feet above the Trout Creek Sandstone Member of the Iles Formation. The Yampa Bed is dated at 72.2 ± .1 mega-annum using the K-Ar method. Regionally, the Yampa Bed is a 0.5- to 5-ft-thick, regionally persistent tonstein that can be readily identified in several different lithofacies in the lower part of the Williams Fork Formation. The unit is useful as a regional datum in the correlation of facies within the Williams Fork, and it is easily recognized on geophysical logs by its low resistivity response. Evidence suggests that it is a diagenetically altered airfall ash.

Colorado↗
USGS2008

Geologic model for oil and gas assessment of the Kemik-Thomson Play, central North Slope, Alaska

A geologic model was developed to assess undiscovered oil and gas resources in the Kemik-Thomson Play of the Central North Slope, Alaska. In this model, regional erosion during the Early Cretaceous produced an incised valley system on the flanks and crest of the Mikkelsen High and formed the Lower Cretaceous unconformity. Locally derived, coarse-grained siliciclastic and carbonate detritus from eroded Franklinian-age basement rocks, Carboniferous Kekiktuk Conglomerate (of the Endicott Group), Lisburne Group, and Permian-Triassic Sadlerochit Group may have accumulated in the incised valleys during lowstand and transgression, forming potential reservoirs in the Lower Cretaceous Kemik Sandstone and Thomson sandstone (informal term). Continued transgression resulted in the deposition of the mudstones of the over-lying Cretaceous pebble shale unit and Hue Shale, which form top seals to the potential reservoirs. Petroleum from thermally mature facies of the Triassic Shublik Formation, Jurassic Kingak Shale, Hue Shale (and pebble shale unit), and the Cretaceous-Tertiary Canning Formation might have charged Thomson and Kemik sandstone reservoirs in this play during the Tertiary. The success of this play depends largely upon the presence of reservoir-quality units in the Kemik Sandstone and Thomson sandstone.

Alaska↗
USGS2017

Mineral-deposit model for lithium-cesium-tantalum pegmatites

Lithium-cesium-tantalum (LCT) pegmatites comprise a compositionally defined subset of granitic pegmatites. The major minerals are quartz, potassium feldspar, albite, and muscovite; typical accessory minerals include biotite, garnet, tourmaline, and apatite. The principal lithium ore minerals are spodumene, petalite, and lepidolite; cesium mostly comes from pollucite; and tantalum mostly comes from columbite-tantalite. Tin ore as cassiterite and beryllium ore as beryl also occur in LCT pegmatites, as do a number of gemstones and high-value museum specimens of rare minerals. Individual crystals in LCT pegmatites can be enormous: the largest spodumene was 14 meters long, the largest beryl was 18 meters long, and the largest potassium feldspar was 49 meters long. Lithium-cesium-tantalum pegmatites account for about one-fourth of the world’s lithium production, most of the tantalum production, and all of the cesium production. Giant deposits include Tanco in Canada, Greenbushes in Australia, and Bikita in Zimbabwe. The largest lithium pegmatite in the United States, at King’s Mountain, North Carolina, is no longer being mined although large reserves of lithium remain. Depending on size and attitude of the pegmatite, a variety of mining techniques are used, including artisanal surface mining, open-pit surface mining, small underground workings, and large underground operations using room-and-pillar design. In favorable circumstances, what would otherwise be gangue minerals (quartz, potassium feldspar, albite, and muscovite) can be mined along with lithium and (or) tantalum as coproducts. Most LCT pegmatites are hosted in metamorphosed supracrustal rocks in the upper greenschist to lower amphibolite facies. Lithium-cesium-tantalum pegmatite intrusions generally are emplaced late during orogeny, with emplacement being controlled by pre-existing structures. Typically, they crop out near evolved, peraluminous granites and leucogranites from which they are inferred to be derived by fractional crystallization. In cases where a parental granite pluton is not exposed, one is inferred to lie at depth. Lithium-cesium-tantalum LCT pegmatite melts are enriched in fluxing components including H2O, F, P, and B, which depress the solidus temperature, lower the density, and increase rates of ionic diffusion. This, in turn, enables pegmatites to form thin dikes and massive crystals despite having a felsic composition and temperatures that are significantly lower than ordinary granitic melts. Lithium-cesium-tantalum pegmatites crystallized at remarkably low temperatures (about 350–550 °C) in a remarkably short time (days to years). Lithium-cesium-tantalum pegmatites form in orogenic hinterlands as products of plate convergence. Most formed during collisional orogeny (for example, Kings Mountain district, North Carolina). Specific causes of LCT pegmatite-related magmatism could include: ordinary arc processes; over thickening of continental crust during collision or subduction; slab breakoff during or after collision; slab delamination before, during, or after collision; and late collisional extensional collapse and consequent decompression melting. Lithium-cesium-tantalum pegmatite deposits are present in all continents including Antarctica and in rocks spanning 3 billion years of Earth history. The global age distribution of LCT pegmatites is similar to those of common pegmatites, orogenic granites, and detrital zircons. Peak times of LCT pegmatite genesis at about 2640, 1800, 960, 485, and 310 Ma (million years before present) correspond to times of collisional orogeny and supercontinent assembly. Between these pulses were long intervals when few or no LCT pegmatites formed. These minima overlap with supercontinent tenures at ca. 2450–2225, 1625–1000, 875–725, and 250–200 Ma. Exploration and assessment for LCT pegmatites are guided by a number of observations. In frontier areas where exploration has been minimal at best, the key first-order criteria are an orogenic hinterland setting, appropriate regional metamorphic grades, and the presence of evolved granites and common granitic pegmatites. New LCT pegmatites are most likely to be found near known deposits. Pegmatites tend to show a regional mineralogical and geochemical zoning pattern with respect to the inferred parental granite, with the greatest enrichment in the more distal pegmatites. Mineral-chemical trends in common pegmatites that can point toward an evolved LCT pegmatite include: increasing rubidium in potassium feldspar, increasing lithium in white mica, increasing manganese in garnet, and increasing tantalum and manganese in columbite-tantalite. Most LCT pegmatite bodies show a distinctive internal zonation featuring four zones: border, wall, intermediate (where lithium, cesium, and tantalum are generally concentrated), and core. This zonation is expressed both in cross section and map view; thus, what may appear to be a common pegmatite may instead be the edge of a mineralized body. Neither lithium-cesium-tantalum pegmatites nor their parental granites are likely to cause serious environmental concerns. Soils and country rock surrounding a LCT pegmatite, as well as waste from mining operations, may be enriched in characteristic elements relative to global average soil and bedrock values. These elements may include lithium, cesium, tantalum, beryllium, boron, fluorine, phosphorus, manganese, gallium, rubidium, niobium, tin, and hafnium. Among this suite of elements, however, the only ones that might present a concern for environmental health are beryllium and fluorine, which are included in the U.S. Environmental Protection Agency drinking-water regulations with maximum contaminant levels of 4 micrograms per liter and 4 milligrams per liter, respectively.

Scientific Investigations Report↗
USGS2011

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↗
USGS2014

Geochronology and correlation of Tertiary volcanic and intrusive rocks in part of the southern Toquima Range, Nye County, Nevada

Extensive volcanic and intrusive igneous activity, partly localized along regional structural zones, characterized the southern Toquima Range, Nevada, in the late Eocene, Oligocene, and Miocene. The general chronology of igneous activity has been defined previously. This major episode of Tertiary magmatism began with emplacement of a variety of intrusive rocks, followed by formation of nine major calderas and associated with voluminous extrusive and additional intrusive activity. Emplacement of volcanic eruptive and collapse megabreccias accompanied formation of some calderas. Penecontemporaneous volcanism in central Nevada resulted in deposition of distally derived outflow facies ash-flow tuff units that are interleaved in the Toquima Range with proximally derived ash-flow tuffs. Eruption of the Northumberland Tuff in the north part of the southern Toquima Range and collapse of the Northumberland caldera occurred about 32.3 million years ago. The poorly defined Corcoran Canyon caldera farther to the southeast formed following eruption of the tuff of Corcoran Canyon about 27.2 million years ago. The Big Ten Peak caldera in the south part of the southern Toquima Range Tertiary volcanic complex formed about 27 million years ago during eruption of the tuff of Big Ten Peak and associated air-fall tuffs. The inferred Ryecroft Canyon caldera formed in the south end of the Monitor Valley adjacent to the southern Toquima Range and just north of the Big Ten Peak caldera in response to eruption of the tuff of Ryecroft Canyon about 27 million years ago, and the Moores Creek caldera just south of the Northumberland caldera developed at about the same time. Eruption of the tuff of Mount Jefferson about 26.8 million years ago was accompanied by collapse of the Mount Jefferson caldera in the central part of the southern Toquima Range. An inferred caldera, mostly buried beneath alluvium of Big Smoky Valley southwest of the Mount Jefferson caldera, formed about 26.5 million years ago with eruption of the tuff of Round Mountain. The Manhattan caldera south of the Mount Jefferson caldera and northwest of the Big Ten Peak caldera formed in association with eruption of a series of tuffs, principally the Round Rock Formation, mostly ash-flow tuff, about 24.4 million years ago. Extensive 40 Ar/ 39 Ar dating of about 60 samples that represent many of the Tertiary extrusive and intrusive rocks in the southern Toquima Range provides precise ages that refine the chronology of previously dated units. New geochronologic data indicate that the petrogenetically related Corcoran Canyon, Ryecroft Canyon, and Mount Jefferson calderas formed during a period of about 560,000 years. Electron microprobe analyses of phenocrysts from 20 samples of six dated units underscore inferred petrogenetic relations among some of these units. In particular, compositions of augite, hornblende, and biotite in tuffs erupted from the Corcoran Canyon, Ryecroft Canyon, and Mount Jefferson calderas are similar, which suggests that magmas represented by these tuffs have similar petrogenetic histories. The unique occurrence of hypersthene in Isom-type tuff confirms its derivation from a source beyond the southern Toquima Range.

Nevada↗
USGS2014

Quality of groundwater in the Denver Basin aquifer system, Colorado, 2003-5

Groundwater resources from alluvial and bedrock aquifers of the Denver Basin are critical for municipal, domestic, and agricultural uses in Colorado along the eastern front of the Rocky Mountains. Rapid and widespread urban development, primarily along the western boundary of the Denver Basin, has approximately doubled the population since about 1970, and much of the population depends on groundwater for water supply. As part of the National Water-Quality Assessment Program, the U.S. Geological Survey conducted groundwater-quality studies during 2003–5 in the Denver Basin aquifer system to characterize water quality of shallow groundwater at the water table and of the bedrock aquifers, which are important drinking-water resources. For the Denver Basin, water-quality constituents of concern for human health or because they might otherwise limit use of water include total dissolved solids, fluoride, sulfate, nitrate, iron, manganese, selenium, radon, uranium, arsenic, pesticides, and volatile organic compounds. For the water-table studies, two monitoring-well networks were installed and sampled beneath agricultural (31 wells) and urban (29 wells) land uses at or just below the water table in either alluvial material or near-surface bedrock. For the bedrock-aquifer studies, domestic- and municipal-supply wells completed in the bedrock aquifers were sampled. The bedrock aquifers, stratigraphically from youngest (shallowest) to oldest (deepest), are the Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers. The extensive dataset collected from wells completed in the bedrock aquifers (79 samples) provides the opportunity to evaluate factors and processes affecting water quality and to establish a baseline that can be used to characterize future changes in groundwater quality. Groundwater samples were analyzed for inorganic, organic, isotopic, and age-dating constituents and tracers. This report discusses spatial and statistical distributions of chemical constituents and evaluates natural and human-related processes that affect water quality. Findings are synthesized to assess the vulnerability of the Denver Basin aquifer system to groundwater contamination. The chemistry of groundwater samples collected from the water-table wells was generally different from that of samples collected from the bedrock-aquifer wells. Samples from the water-table wells tended to have higher concentrations of total dissolved solids and most major ions. Concentrations of several constituents with potential human-health concerns, including nitrate, selenium, uranium, and arsenic, decreased with depth and were highest in samples from the water-table wells. Exceedances of drinking-water standards and water-quality benchmarks were more frequently associated with shallow groundwater samples; concentrations of total dissolved solids and sulfate exceeded water-quality benchmarks for about half or more of samples from the water-table wells. The sediments and rocks of the Denver Basin are natural sources of the trace elements selenium, uranium, and arsenic, which affect their concentrations in groundwater. Detections of organic contaminants, which are typically indicative of human sources of contamination to groundwater, were more frequent in samples from the water-table wells. Pesticide compounds and volatile organic compounds were detected in 33 and 62 percent, respectively, of water-table well samples. Detected organic contaminant concentrations were much less than the associated drinking-water standards. Samples collected from the bedrock aquifers had lower concentrations of total dissolved solids than did samples collected from the water-table wells, although within the bedrock-aquifer samples, concentrations increased from the Dawson to Denver to Arapahoe to Laramie-Fox Hills aquifers. Concentrations of total dissolved solids and many constituents varied spatially and with depth in the bedrock aquifers, likely as a result of ion-exchange and oxidation-reduction reactions, which are important processes affecting water quality. Major-ion chemistry generally evolved from a calcium-bicarbonate to calcium-sulfate composition, with some sodium-bicarbonate and sodium-sulfate facies in the deeper bedrock aquifers, likely resulting from longer residence times and more extensive water-rock interaction. Oxidation-reduction conditions generally evolved from oxic at the water table to anoxic with increasing depth in the bedrock aquifers. Most samples from the bedrock aquifers were anoxic. Exceedances of drinking-water standards and water-quality benchmarks for the bedrock aquifers occurred in 1 percent or less of samples for nitrate, selenium, or arsenic; there were no exceedances for uranium. Exceedances for total dissolved solids, sulfate, manganese, and iron were generally between about 10 and 20 percent for the bedrock-aquifer samples. Radon concentrations, which were only measured in samples collected from two of the bedrock aquifers, exceeded the lower proposed drinking-water standard for more than 90 percent of samples but exceeded the higher alternative standard for less than 5 percent of samples. Pesticide compounds and volatile organic compounds were detected in 3 and 22 percent, respectively, of bedrock-aquifer samples, all at concentrations that were that were much less than drinking-water standards. Water-quality data were synthesized to evaluate factors that affect spatial and depth variability in water quality and to assess aquifer vulnerability to contaminants from geologic materials and those of human origin. The quality of shallow groundwater in the alluvial aquifer and shallow bedrock aquifer system has been adversely affected by development of agricultural and urban areas. Land use has altered the pattern and composition of recharge. Increased recharge from irrigation water has mobilized dissolved constituents and increased concentrations in the shallow groundwater. Concentrations of most constituents associated with poor or degraded water quality in shallow groundwater decreased with depth; many of these constituents are not geochemically conservative and are affected by geochemical reactions such as oxidation-reduction reactions. Groundwater age tracers provide additional insight into aquifer vulnerability and help determine if young groundwater of potentially poor quality has migrated to deeper parts of the bedrock aquifers used for drinking-water supply. Age-tracer results were used to group samples into categories of young, mixed, and old groundwater. Groundwater ages transitioned from mostly young in the water-table wells to mostly mixed in the shallowest bedrock aquifer, the Dawson aquifer, to mostly old in the deeper bedrock aquifers. Although the bedrock aquifers are mostly old groundwater of good water quality, several lines of evidence indicate that young, contaminant-bearing recharge has reached shallow to moderate depths in some areas of the bedrock aquifers. The Dawson aquifer is the most vulnerable of the bedrock aquifers to contamination, but results indicate that the older (deeper) bedrock aquifers are also vulnerable to groundwater contamination and that mixing with young recharge has occurred in some areas. Heavy pumping has caused water-level declines in the bedrock aquifers in some parts of the Denver Basin, which has the potential to enhance the transport of contaminants from overlying units. Results of this study are consistent with the existing conceptual understanding of aquifer processes and groundwater issues in the Denver Basin and add new insight into the vulnerability of the bedrock aquifers to groundwater contamination.

Colorado↗
USGS2017

Three-dimensional hydrogeologic framework model of the Rio Grande transboundary region of New Mexico and Texas, USA, and northern Chihuahua, Mexico

As part of a U.S. Geological Survey study in cooperation with the Bureau of Reclamation, a digital three-dimensional hydrogeologic framework model was constructed for the Rio Grande transboundary region of New Mexico and Texas, USA, and northern Chihuahua, Mexico. This model was constructed to define the aquifer system geometry and subsurface lithologic characteristics and distribution for use in a regional numerical hydrologic model. The model includes five hydrostratigraphic units: river channel alluvium, three informal subdivisions of Santa Fe Group basin fill, and an undivided pre-Santa Fe Group bedrock unit. Model input data were compiled from published cross sections, well data, structure contour maps, selected geophysical data, and contiguous compilations of surficial geology and structural features in the study area. These data were used to construct faulted surfaces that represent the upper and lower subsurface hydrostratigraphic unit boundaries. The digital three-dimensional hydrogeologic framework model is constructed through combining faults, the elevation of the tops of each hydrostratigraphic unit, and boundary lines depicting the subsurface extent of each hydrostratigraphic unit. The framework also compiles a digital representation of the distribution of sedimentary facies within each hydrostratigraphic unit. The digital three-dimensional hydrogeologic model reproduces with reasonable accuracy the previously published subsurface hydrogeologic conceptualization of the aquifer system and represents the large-scale geometry of the subsurface aquifers. The model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

Chihuahua, New Mexico, Texas↗
USGS2018

Postglacial eruptive history and geochemistry of Semisopochnoi volcano, western Aleutian Islands, Alaska

Semisopochnoi Island, located in the Rat Islands group of the western Aleutian Islands and Aleutian volcanic arc, is a roughly circular island composed of scattered volcanic vents, the prominent caldera of Semisopochnoi volcano, and older, ancestral volcanic rocks. The oldest rocks on the island are gently radially dipping lavas that are the remnants of a shield volcano and of Ragged Top, which is an eroded stratocone southeast of the current caldera. None of these oldest rocks have been dated, but they all are likely Pleistocene in age. Anvil Peak, to the caldera’s north, has the morphology of a young stratocone and is latest Pleistocene to early Holocene in age. The oldest recognized Holocene deposits are those of the caldera-forming eruption, which produced the 7- by 6-km caldera in the center of the island, left nonwelded ignimbrite in valleys below the edifice, and left welded ignimbrite high on its flanks. The caldera-forming eruption produced rocks showing a range of intermediate whole-rock compositions throughout the eruption sequence, although a majority of clasts analyzed form a fairly tight cluster on SiO 2 -variation diagrams at 62.9 to 63.4 weight percent SiO 2 . This clustering of compositions at about 63 weight percent SiO 2 includes black, dense, obsidian-like clasts, as well as tan, variably oxidized, highly inflated pumice clasts. The best estimate for the timing of the eruption is from a soil dated at 6,920±60 14 C years before present underlying a thin facies of the ignimbrite deposit on the island’s north coast. Shortly after the caldera-forming eruption, two scoria cones on the northwest flank of the volcano outside the caldera, Ringworm crater and Threequarter Cone, simultaneously erupted small volumes of andesite. The oldest intracaldera lavas, on the floor of the caldera, are andesitic to dacitic, but are mostly covered by younger lavas and tephras. These intracaldera lavas include the basaltic andesites of small Windy cone, as well as the more voluminous basaltic andesites of three-peaked Mount Cerberus, which takes up most of the west half of the caldera and has erupted lavas that flowed to the sea on the southwestern coast of the island. Apparently active at the same time as Mount Cerberus, extracaldera Sugarloaf Peak at the southern point of the island has exclusively erupted basalts. Its young satellite peak, Sugarloaf Head, has erupted morphologically young lavas and cinder cones and may be the source of the last historical eruption in 1987. Several tephra sections on the east half of the island record as many as 50 tephras, mostly from Mount Cerberus, Sugarloaf Peak, and Sugarloaf Head, over the past several thousand years. Eruptive products of Semisopochnoi Island show an overall compositional range of basalt to dacite, though basaltic andesite and andesite constitute the largest proportions of rock types. They are tholeiitic, low to medium K, and have geochemical characteristics typical of magmatic arcs. The earliest Pleistocene lavas are mostly basalts that show the greatest geochemical diversity, as illustrated by, for example, La N /Yb N ratios of 1.9 to 3.5, suggesting fluctuations in the magma source region over the hundreds of thousands of years recorded by these older lavas. The Holocene rocks, in contrast, follow arrays in compositional space that suggest crystallization differentiation from discrete, subtly different batches of magma under varying pressure and temperature conditions. Increasingly negative Eu anomalies and an only modestly increasing alumina saturation index value with differentiation suggest that plagioclase and mafic silicates (amphibole and pyroxene) were involved to varying degrees in fractional crystallization to produce Semisopochnoi’s magmatic diversity. The crystal-poor, andesitic magmas that erupted during caldera formation likely separated from a plagioclase-, amphibole-, and clinopyroxene-dominated crystal residue in the upper crust at less than 900 °C, possibly following a period of decreased magmatic flux. During the Holocene, basaltic Sugarloaf Peak appears to bypass any upper crustal magmatic storage region and erupt crystal-rich basalts. Recent seismic swarms and long-lived warm springs attest to ongoing magmatic activity. The Holocene eruptive record at Semisopochnoi volcano is one of diverse eruptive styles as well as frequent eruptions from multiple vents located within and outside the caldera. The number and diversity of postcaldera vents means that the sites of future eruptions cannot be predicted with certainty. Future eruptions of ash similar in magnitude to the VEI 3 or less eruptions recorded in the documented tephra deposits would pose a hazard to aircraft in the region.

Alaska↗
USGS2022

Revisions to the Virginia Coastal Plain hydrogeologic framework southwest of the James River

New drilling information reveals that altitudes of some hydrogeologic units of the Virginia Coastal Plain aquifer system differ by as much as 50 feet (ft) from those previously known, namely the Aquia and Potomac aquifers, the Potomac confining zone, and the Nanjemoy-Marlboro and Saint Marys confining units. In addition, the lateral margins of some hydrogeologic units are located as much as several miles from previously estimated locations. The largest revisions to unit margins were for the Aquia aquifer and the Nanjemoy-Marlboro and Saint Marys confining units. Interpretation of new geophysical logs, sediment core, and cuttings as well as revised interpretations to existing data indicate channels and embayments are also preserved on eroded top surfaces of the shallowest hydrogeologic units including the Yorktown confining zone, Yorktown-Eastover aquifer, Saint Marys confining unit, Potomac confining zone, and Potomac aquifer. Enhanced details on the configuration of part of the aquifer system southwest of the James River are provided by sediment cores and cuttings as well as geophysical logs from 36 recently drilled boreholes. These, along with reinterpretation of data from 93 preexisting boreholes, form the basis for revised top-surface altitudes and margins of hydrogeologic units beneath parts of Prince George, Surry, Sussex, Isle of Wight, and Southampton Counties and the cities of Franklin and Suffolk. Groundwater withdrawals in the Virginia Coastal Plain cause widespread water-level declines, create the potential for saltwater intrusion, and contribute to regionwide land subsidence. A description of the aquifer system, termed a hydrogeologic framework, was developed by the U.S. Geological Survey in 2006 and provides information needed to base withdrawal-permitting decisions by the Virginia Department of Environmental Quality. This revision of part of the hydrogeologic framework southwest of the James River is based on interpretations of both new and previously analyzed borehole data. The revision is strictly confined to the study area extent and hydrogeologic units not found within the study area were not revised and are not discussed in this report. The newly determined hydrogeologic-unit altitudes and margins have implications for groundwater-withdrawal permitting. New interpretations have found that the Yorktown Eastover aquifer is absent in the southwestern part of the City of Suffolk, owing to what is most likely an isolated area of sediment-texture facies change. Most notably, the top-surface altitudes of the Aquia and Potomac aquifers have been lowered by as much as 50 ft from previous interpretations. This means that wells previously believed to be screened in the top of the Potomac aquifer could, based on these new interpretations, be screened in the bottom of the Aquia aquifer. These changes to aquifers in which wells are screened means that there is potentially more room in the groundwater withdrawal permitting for the Potomac aquifer, the largest and most productive aquifer in Virginia, and overpumping occurring in the Aquia aquifer.

Virginia↗
USGS1956

Uranium-bearing coal in the central part of the Great Divide basin, Sweetwater County, Wyoming

Field work leading to this report was done by the U.S. Geological Survey for the Division of Raw Materials of the U.S. Atomic Energy Commission. Nearly 24 townships were mapped in the central part of the Great Divide Basin, Sweetwater County, Wyoming. Fourteen of these townships contain outcrops of uranium-bearing coal. Thirty coal beds were mapped, but only seven of them have uranium-bearing coal reserves as defined in this report. Coal beds 2.5 or more feet thick are considered in calculating coal reserves, and of these, only beds containing 0.003 or more percent uranium are considered in calculating reserves of uranium in coal. Reserves of uranium in coal ash include those beds 2.5 or more feet thick that contain 0.015 or more percent uranium in coal ash. Measured and indicated coal reserves total about 700,000,000 short tons which contain about 2,600 short tons of uranium in the coal, or about 2,400 short tons of uranium in the coal ash. Strippable reserves, defined as reserves in beds beneath 60 or less feet of overburden, are about 250,000,000 short tons of coal containing about 1,100 short tons of uranium in coal, or about 600 tons of uranium in coal ash. The thickest coal beds underlie a relatively narrow belt that trends northwest and coincides approximately with the axis of the Red Desert syncline. The coal beds contain the most uranium on the east flank of the syncline near the southwesternmost edge of the Battle Spring formation (new). This formation is of early and middle Eocene age and consists predominantly of very coarse-grained arkosic sandstone which is highly permeable. It intertongues southwestward with the Tess permeable Green River and Wasatch formations. The Green River formation consists from youngest to oldest of the Morrow Creek and Laney shale members and the Tipton and Luman (new) tongues. The Wasatch formation interfingers with the Green River formation and consists from youngest to oldest of the Cathedral Bluffs, Niland, and Red Desert tongues. The latter two are here recognized for the first time and contain all the coal beds in the Wasatch formation. The Morrow Creek member of the Green River formation is of middle Eocene age. The Laney member and Cathedral Bluffs tongue of the Green River and Wasatch formation, respectively intertongue and are in part equivalent, and are of early or middle Eocene age. The other units are of early Eocene age. A broad gentle arch which is either the eastern extension of the Wamsutter arch, or a separate arch en echelon to it, separates the Washakie Basin the southeastern part of the map area from the Red Desert syncline in the northeastern part. The Red Desert syncline plunges gently northwest into the nearly circular structural Niland basin. The south flank of the Wamsutter(?) arch dips southeastward at an average rate of about 230 feet per mile. The northern flank of the arch dips northeastward at an average rate of about 140 feet per mile. North of the Niland basin, the structure is more complex and is dominated by a northwestwardtrending graben which parallels and includes the Cyclone Rim in the northern part of die map area. Schroeckingerite deposits in the northeastern part of the map area lie within this graben in a sequence of arkosic sandstone and clay shale stratigraphically equivalent to the lower part of the Cathedral Bluffs tongue of the Wasatch formation and to the uppermost part of the Tipton tongue of the Green River formation. Weakly radioactive tuffaceous sandstone beds of the Browns Park formation that probably once blanketed the entire Great Divide Basin are preserved a short distance north and northwest of the schroeckingerite deposits. The geologic settings of the schroeckingerite and uranium-bearing coal are similar and their source of uranium is probably the same. The uranium probably was leached from tuffaceous beds in the Browns Park formation and carried to its present site of deposition by groundwater whose circulation was guided by structure and changes in facies.

Wyoming↗
USGS1955

Radioactivity and uranium content of the Sharon Springs member of the Pierre shale and associated rocks in western Kansas and eastern Colorado

As a part of the Geological Survey's program of investigating uranium-bearing carbonaceous rocks on behalf of the Division of Raw Materials of the U.S. Atomic Energy Commission, a reconnaissance of the Sharon Springs member of the Pierre shale in western Kansas and eastern Colorado was conducted during 1954. The Sharon Springs member of the Pierre shale and its lateral equivalents ranges from 155 to about 500 feet in thickness and generally contains about 0.001 percent uranium, but some beds contain larger amounts. A 6-foot thick shale bed in Cheyenne County, Colo., contains about 0.006 percent uranium, a 4 1/2-foot thick sequence of beds in Crowley County, Colo., is estimated to contain between 0.004 and 0.005 percent uranium, and a 3 1/2-foot thick sequence of beds in Kiowa County, Colo., contains about 0.004 percent uranium. At several outcrop localities, sequences of beds as much as 9 1/2 feet thick contain about 0.003 percent uranium. Data from wells indicate that the 4 1/2-foot thick sequence of beds in Crowley County, Colo., may have a lateral extent of at least 5 1/2 miles. A gamma-ray log of a well in Yuma County, Colo., indicates the presence of a sequence of beds 66 feet thick which contains 0.005 to 0.010 percent equivalent uranium. No definite pattern of areal distribution of radioactivity and uranium content in the Sharon Springs is indicated by available data. Lateral variation in uranium content of individual beds was not noted in outcrops, which seldom extend more than 150 feet, but subsurface data from gamma-ray logs of wells indicate that both the maximum radioactivity and the thickness of radioactive beds are variable within distances of a few miles. Vertical variation in radioactivity and uranium content of the more radioactive beds is usually abrupt, but in the rocks as a whole the range of uranium content is so small that large variations in content are absent. In most of the gamma-ray logs examined there is only part of the sequence of rocks comprising the Pierre shale and Niobrara formation that exhibits radioactivity in excess of the average radioactivity of the two formations. Comparison of features of gamma-ray logs of wells in north-eastern Colorado suggests that the most radioactive part referred to above is a laterally correlatable sequence of beds. The stratigraphic position of the radioactive unit relative to the Pierre shale-Niobrara formation contact in oil industry scout reports, as identified from electric logs and wells, is variable within short distances. This may indicate that some of the Pierre-Niobrara contacts picked from electric logs may not correspond to the boundary that would be selected by examination of the rocks themselves, or may indicate that there is a facies relationship between teh lowermost part of the Pierre shale and the uppermost part of the Niobrara formation.

Colorado;Kansas↗
USGS1957

The Gas Hills uranium district and some probable controls for ore deposition

Uranium deposits occur in the upper coarse-grained facies of the Wind River formation of Eocene age in the Gas Hills district of the southern part of the Wind River Basin. Some of the principal deposits lie below the water table in the unoxidized zone and consist of uraninite and coffinite occurring as interstitial fillings in irregular blanket-like bodies. In the near-surface deposits that lie above the water table, the common yellow uranium minerals consist of uranium phosphates, silicates, and hydrous oxides. The black unoxidized uraninite -coffinite ores show enrichment of molybdenum, arsenic, and selenium when compared to the barren sandstone. Probable geologic controls for ore deposits include: 1) permeable sediments that allowed passage of ore-bearing solutions; 2) numerous faults that acted as impermeable barriers impounding the ore -bearing solutions; 3) locally abundant pyrite, carbonaceous material, and natuial gas containing hydrogen sulfide that might provide a favorable environment for precipitation of uranium. Field and laboratory evidence indicate that the uranium deposits in the Gas Hills district are very young and related to the post-Miocene to Pleistocene regional tilting to the south associated with the collapse of the Granite Mountains fault block. This may have stopped or reversed ground water movement from a northward (basinward) direction and alkaline ground water rich in carbonate could have carried the uranium into the favorable environment that induced precipitation.

Trace Elements Memorandum↗
USGS2000

Geologic framework and hydrogeologic characteristics of the Edwards Aquifer outcrop, Medina County, Texas

The hydrogeologic subdivisions of the Edwards aquifer outcrop in Medina County generally are porous and permeable. The most porous and permeable appear to be hydrogeologic subdivision VI, the Kirschberg evaporite member of the Kainer Formation; and hydrogeologic subdivision III, the leached and collapsed members, undivided, of the Person Formation. The most porous and permeable rocks of the Devils River Formation in Medina County appear to be in the top layer. The upper member of the Glen Rose Limestone, the lower confining unit, has much less porosity and permeability than that observed in the Edwards aquifer. The Edwards aquifer has relatively large porosity and permeability resulting, in part, from the development or redistribution of secondary porosity. Lithology, stratigraphy, diagenesis, and karstification account for the effective porosity and permeability in the Edwards aquifer outcrop. Karst features that can greatly enhance effective porosity and permeability in the Edwards aquifer outcrop include sinkholes, dolines, and caves. The Edwards aquifer rocks in Medina County change from the eight-member Edwards Group to the essentially indivisible Devils River Formation. The facies change occurs along a line extending northwestward from just south of Medina Lake.

Texas↗
USGS2001

Stratigraphy and vertical hydraulic conductivity of the St. Francois confining unit in townships 25-27 N. and ranges 01-02 W., southeastern Missouri

The St. Francois confining unit (DerbyDoerun Dolomite and Davis Formation) lies beneath the Ozark aquifer (Jefferson City Dolomite to the Potosi Dolomite) and impedes the circulation of water between the overlying Ozark aquifer and the underlying St. Francois aquifer (Bonneterre Formation and Lamotte Sandstone). The Bonneterre Formation is the potential host formation for lead-zinc deposits in the area. There is concern that mine dewatering in the Bonneterre Formation could lower water levels in the Ozark aquifer. To address this concern, the vertical hydraulic conductivity of the St. Francois confining unit in six townships (T. 25-27 N. and R. 01- 02 W.) of Oregon, Carter, and Ripley Counties of southeastern Missouri was evaluated by describing the stratigraphy and measuring the vertical hydraulic conductivity of core samples. The Davis Formation is an intrashelf basin facies consisting of a series of shales interbedded with shaley limestones, shale-free limestones, and local dolostones, and ranges from 24 to 320ft (feet) thick, but typically the thickness is 100 to 200 ft. Shale-dominant sequences can be tens of feet thick, and contain as much as 90 percent shale. Carbonate-dominant zones may be 70 ft thick or greater. The top of the Davis Formation (based on 56 data points) ranges from 620 to 2,022 ft deep and ranges in altitude from 40ft below sea level in the northern part of the study area to 1, 182 ft below sea level in the southern part of the study area. The Derby-Doerun Dolomite represents a pair of superimposed carbonate ramp cycles. Where present, the basal shaley sequence represents a transition with the Davis Formation. The formation (based on 50 data points) ranges from 50 to 386ft thick, but typically is 120 to 180ft thick in the study area. The top of the DerbyDoerun Dolomite ranges from 495 to 2,020 ft deep (based on 53 data points), and ranges in altitude from 85 ft above sea level to 94 7 ft below sea level. The St. Francois confining unit is thickest in the central and southern parts of the study area. The thickness, as determined by 51 core logs that completely penetrate the unit, ranges from less than 200ft in the northwestern and east-central parts of the study area to 411 ft in the central part, but typically ranges from 270 to 340 ft. The net shale thickness of the confining unit (based on 29 data points) ranges from 1. 7 ft in the east -central part of the study area to 89 ft in the southwest part. These net shale thickness values include the cumulative shale thickness of rock from the top of the Derby-Doerun Dolomite to the base of the False Davis. Vertical hydraulic conductivities of 35 rock core samples from the St. Francois confining unit in the study area range from 7.6 x 10 -15 to 2.1 x 10 -10 ft/s (foot per second). The logarithmic transformed vertical hydraulic conductivities of the Derby-Doerun Dolomite and Davis Formation are similar (p-value = 0.073) using the statistical twosample t-test; however, this p-value approaches the level of significance value of 0.05. The vertical hydraulic conductivity of the Derby-Doerun Dolomite is larger and less variable than the Davis Formation. When grouped by rock type, the vertical hydraulic conductivity of samples that contain carbonate, shale, or both carbonate and shale, are similar. A comparison on the ranked data using the Mann-Whitney test shows the confining unit in the study area is statistically different (p-value = 0.020) from the confining unit in the prospecting area (west and adjacent to the study area). The median value of the vertical hydraulic conductivity data from the study area (6.7 x 10 -13 ft/s) is three times larger than the median vertical hydraulic conductivity value for the prospecting area (2.2 x 10 -13 ft/s ). The interquartile range shows that the variability of the study area data spans one order of magnitude (2.0 x 10 -13 to 2.2 x 10 -12 ft/s) and that the corresponding data from the prospecting area spans nearly two orders of magnitude (3.2 x 10 -14 to 1.1 x 10 -12 ft/s). The ranked vertical hydraulic conductivities of the Derby-Doerun Dolomite in the two areas are statistically similar (p-value = 0.514). The median vertical hydraulic conductivity of the study area data ( 1.2 X 10 -12 ft/s) is about three times greater than the median value of the prospecting area data (4.4 x 10 -13 ft/s). The variability of the data, as shown by the interquartile range, is less in the study area (5.5 x 10 -13 to 2.2 x 10 -12 ft/s; spanning less than one order of magnitude) as compared to the prospecting area (3.2 x 10 -14 to 6.3 x 10 -10 ft/s; spanning over four orders of magnitude). The ranked vertical hydraulic conductivities of the Davis Formation in the two areas show these data sets are statistically similar (p-value = 0.076). The median vertical hydraulic conductivity value of study area samples ( 4.5 x 10 -13 ft/s) is three times greater than the median value of the prospecting area data (1.6 x 10 -13 ft/s). The interquartile range of the study area data spans one order of magnitude (1.2 x 10 -13 to 1.4 x 10 -12 ft/s) and the corresponding data from the prospecting area spans nearly 1.5 orders of magnitude (3.2 x 10 -14 to 7.4 x 10 -13 ft/s). The Mann-Whitney test shows the ranked vertical hydraulic conductivities of each rock type from the study area are statistically similar to the same rock type in the prospecting area [carbonates (p-value = 0.225), shales (p-value = 0.668), and carbonates and shales (p-value = 0.227)]. However, in each of the three cases the study area samples have larger median values and less variability than the prospecting area samples. Because the vertical hydraulic conductivity of the various rock types of the confining unit in the study area are statistically similar, the entire carbonate-shale thickness is the primary factor determining the effectiveness of the confining unit. The range of effective vertical hydraulic conductivity of the St. Francois confining unit in the study area using appropriate minimum and maximum thickness, net shale thickness, and vertical hydraulic conductivities is 3 X 10 -13 to 2 X 10 -12 ft/s. The vertical hydraulic conductivity of the confining unit is small, and the confining unit effectively impedes the ground-water flow between the Ozark aquifer and the St. Francois aquifer, unless preferred- path secondary permeability has developed along faults and fractures that extend through the confining unit.

Missouri↗
USGS2002

Geologic framework of the regional ground-water flow system in the Upper Deschutes Basin, Oregon

Ground water is increasingly relied upon to satisfy the needs of a growing population in the upper Deschutes Basin, Oregon. Hydrogeologic studies are being undertaken to aid in management of the ground-water resource. An understanding of the geologic factors influencing ground-water flow is basic to those investigations. The geology of the area has a direct effect on the occurrence and movement of ground water. The permeability and storage properties of rock material are influenced by the proportion, size, and degree of interconnection of open spaces the rocks contain. These properties are the result of primary geologic processes such as volcanism and sedimentation, as well as subsequent processes such as faulting, weathering, or hydrothermal alteration. The geologic landscape in the study area evolved during about 30 million years of volcanic activity related to a north-south trending volcanic arc, the current manifestation of which are today’s Cascade Range volcanoes. The oldest rock unit in the upper Deschutes Basin study area, the John Day Formation, is a sequence of upper Eocene to lower Miocene volcanic and sedimentary rocks. Weathering and alteration of the rocks has resulted in very low permeability; consequently, the unit forms the hydrologic basement for the regional ground- water flow system throughout much of the area. The Deschutes Formation and age-equivalent deposits that overlie the John Day Formation, in contrast, are highly permeable and are the most widely used ground-water-bearing units in the study area. The Deschutes Formation consists of a variety of volcanic and sedimentary deposits ranging in age from late Miocene to Pliocene (approximately 7.5 to 4.0 million years). Three distinct depositional environments previously described for the formation provide useful hydrogeologic subdivisions. The ancestral Deschutes River deposits and some units within the arc- adjacent alluvial-plain region are among the highest yielding units within the Deschutes Formation, with some wells producing up to a few thousand gallons per minute. Opal Springs basalt, Pelton basalt, and the rhyodacite dome complex near Steelhead Falls are particularly productive subunits within the Deschutes Formation and provide tens to hundreds of cubic feet per second of ground-water discharge to the Deschutes and Crooked Rivers, upstream of Round Butte Dam. Most ground-water recharge in the upper Deschutes Basin occurs in Quaternary deposits of the Cascade Range and Newberry Volcano. These deposits are highly permeable, and the fractured character of the lava flows facilitates rapid infiltration of precipitation and snowmelt, as well as movement of ground water to lower elevations. Additional recharge from canal leakage occurs along sections of unlined canals near Bend, constructed on lava flows from Newberry Volcano. Hydrothermal alteration and secondary mineralization at depth beneath the Cascade Range and Newberry Volcano has drastically reduced the permeability of the material in those regions, effectively restricting most ground water to the strata above the altered rocks. The top of the hydrothermally altered region is considered the base of the regional ground-water system beneath the Cascade Range and Newberry Volcano. Structural features influence ground-water flow within the upper Deschutes Basin mainly by juxtaposing materials with contrasting permeability. This juxtaposition can be caused by fault movement or by the influence of a fault on subsequent deposition. Several depositional centers have formed along the base of fault-line scarps or in grabens within the study area, and the infilling sedimentary deposits have permeability that differs from the surrounding rocks. The effects of faults on ground-water flow may be masked in some areas. For example, the water-table gradient changes slope in the vicinity of the Sisters fault zone, but the slope change also corresponds with a major precipitation gradient change; therefore, any influence of the fault zone is unclear. Geologic units in the Deschutes Basin were divided into several distinct hydrogeologic units. In some instances the units correspond to existing stratigraphic divisions. In other instances, hydrogeologic units correspond to different facies within a single stratigraphic unit or formation. The hydrogeologic units include Quaternary sediment, deposits of the Cascade Range and Newberry Volcano, four zones within the Deschutes Formation and age-equivalent rocks that roughly correspond with depositional environments, and pre-Deschutes-age strata.

Water-Resources Investigations Report↗
USGS2003

Development and calibration of a ground-water flow model for the Sparta Aquifer of southeastern Arkansas and north-central Louisiana and simulated response to withdrawals, 1998-2027

The Sparta aquifer, which consists of the Sparta Sand, in southeastern Arkansas and north-central Louisiana is a major water resource and provides water for municipal, industrial, and agricultural uses. In recent years, the demand in some areas has resulted in withdrawals from the Sparta aquifer that substantially exceed replenishment of the aquifer. Considerable drawdown has occurred in the potentiometric surface forming regional cones of depression as water is removed from storage by withdrawals. These cones of depression are centered beneath the Grand Prairie area and the cities of Pine Bluff and El Dorado in Arkansas, and Monroe in Louisiana. The rate of decline for hydraulic heads in the aquifer has been greater than 1 foot per year for more than a decade in much of southern Arkansas and northern Louisiana where hydraulic heads are now below the top of the Sparta Sand. Continued hydraulic-head declines have caused water users and managers alike to question the ability of the aquifer to supply water for the long term. Concern over protecting the Sparta aquifer as a sustainable resource has resulted in a continued, cooperative effort by the Arkansas Soil and Water Conservation Commission, U.S. Army Corps of Engineers, and the U.S. Geological Survey to develop, maintain, and utilize numerical ground-water flow models to manage and further analyze the ground-water system. The work presented in this report describes the development and calibration of a ground-water flow model representing the Sparta aquifer to simulate observed hydraulic heads, documents major differences in the current Sparta model compared to the previous Sparta model calibrated in the mid-1980's, and presents the results of three hypothetical future withdrawal scenarios. The current Sparta model-a regional scale, three-dimensional numerical ground-water flow model-was constructed and calibrated using available hydrogeologic, hydraulic, and water-use data from 1898 to 1997. Significant changes from the previous model include grid rediscretization of the aquifer, extension of the active model area northward beyond the Cane River Formation facies change, and representation of model boundaries. The current model was calibrated with the aid of parameter estimation, a nonlinear regression technique, combined with trial and error parameter adjustment using a total of 795 observations from 316 wells over 4 different years-1970, 1985, 1990, and 1997. The calibration data set provides broad spatial and temporal coverage of aquifer conditions. Analysis of the residual statistics, spatial distribution of residuals, simulated compared to observed hydrographs, and simulated compared to observed potentiometric surfaces were used to analyze the ability of the calibrated model to simulate aquifer conditions within acceptable error. The calibrated model has a root mean square error of 18 feet for all observations, an improvement of more than 12 feet from the previous model. The current Sparta model was used to predict the effects of three hypothetical withdrawal scenarios on hydraulic heads over the period 1998-2027 with one of those extended indefinitely until equilibrium conditions were attained, or steady state. In scenario 1a, withdrawals representing the time period from 1990 to 1997 was held constant for 30 years from 1998 to 2027. Hydraulic heads in the middle of the cone of depression centered on El Dorado decreased by 10 feet from the 1997 simulation to 222 feet below NGVD of 1929 in 2027. Hydraulic heads in the Pine Bluff cone of depression showed a greater decline from 61 feet below NGVD of 1929 to 78 feet below NGVD of 1929 in the center of the cone. With these same withdrawals extended to steady state (scenario 1b), hydraulic heads in the Pine Bluff cone of depression center declined an 2 Development and Calibration of a Ground-Water Flow Model for the Sparta Aquifer of Southeastern Arkansas and North-Central Louisiana and Simulated Response to Withdrawa

Water-Resources Investigations Report↗