SEARCH · USGS Science

Results for “Quaternary”

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

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54Linked to original sources

USGS2015

California State Waters Map Series: Offshore of Refugio Beach, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Refugio Beach map lies within the western Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 0.5 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The coastal zone of this map area lies at the steep flank of the Santa Ynez Mountains. The crest of the range, which lies about 8 km from the shoreline (north of the map area), has a maximum elevation of about 780 m. This area is largely open space, partly used for livestock grazing, with no significant towns or population centers. Highway 101 crosses the map area, adjacent to and within a few hundred meters of the shoreline. The most significant developments are the recreational state beaches at El Capitan Beach and Refugio Beach. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Refugio Beach map area lies in the west-central part of the Santa Barbara littoral cell, which is characterized by west-to-east transport of sediment from Point Arguello on the northwest to Hueneme and Mugu Canyons on the southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central part of the littoral cell, including the map area, is mainly from relatively small coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Cañada del Capitan, Tajiquas Creek, Arroyo Hondo, Cañada del Molino, and several unnamed canyons and creeks. The Santa Ynez and Santa Maria Rivers, the mouths of which are 80 to 120 km northwest of the map area, are not considered significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 70 km to the southeast of the map area, are not sediment sources for the map area. The offshore part of the map area consists of relatively flat and shallow continental shelf, which dips gently seaward (about 0.8° to 1.0°) so that water depths at the shelf break, roughly coincident with the California’s State Waters limit, are about 80 to 100 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. In the map area, the shelf break is at depths of about 90 m and lies about 5.6 to 6.4 km offshore. Beyond the shelf break, the slope is steep (as much as about 7°) and unstable. Several submarine landslides, including the large (130 km 2 ) Goleta landslide complex, have been documented offshore of Goleta, a few kilometers east of the map area. This compound slump complex may have been initiated more than 200,000 year ago, but it also includes three recent failures that may have been generated 8,000 to 10,000 years ago. A local, 10-m-high tsunami may have been generated from these failure events. Small folds related to local faulting are superimposed on the homocline that makes up the south flank of the Santa Ynez Mountains. One of these superimposed anticlines hosts the Molino gas field, which was discovered in 1962 and subsequently developed through directional drilling from onshore wells. The Oligocene Sespe and Vaqueros Formations are the reservoirs in the Molino gas field, and the map area includes numerous seafloor hydrocarbon seeps. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft, unconsolidated sediment interspersed with isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Refugio Beach map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, primarily are composed of soft sediment interrupted by a few carbonate mounds. This homogeneous seafloor of sediment and low-relief bedrock provides promising habitat for groundfish, crabs, shrimp, and other marine benthic organisms.

California↗

USGS2017

Geologic map of the Providence Mountains in parts of the Fountain Peak and adjacent 7.5' quadrangles, San Bernardino County, California

Introduction The Providence Mountains are in the eastern Mojave Desert about 60 km southeast of Baker, San Bernardino County, California. This range, which is noted for its prominent cliffs of Paleozoic limestone, is part of a northeast-trending belt of mountainous terrain more than 100 km long that also includes the Granite Mountains, Mid Hills, and New York Mountains. Providence Mountains State Recreation Area encompasses part of the range, the remainder of which is within Mojave National Preserve, a large parcel of land administered by the National Park Service. Access to the Providence Mountains is by secondary roads leading south and north from Interstate Highways 15 and 40, respectively, which bound the main part of Mojave National Preserve. The geologic map presented here includes most of Providence Mountains State Recreation Area and land that surrounds it on the north, west, and south. This area covers most of the Fountain Peak 7.5′ quadrangle and small adjacent parts of the Hayden quadrangle to the north, the Columbia Mountain quadrangle to the northeast, and the Colton Well quadrangle to the east. The map area includes representative outcrops of most of the major geologic elements of the Providence Mountains, including gneissic Paleoproterozoic basement rocks, a thick overlying sequence of Neoproterozoic to Triassic sedimentary rocks, Jurassic rhyolite that intrudes and overlies the sedimentary rocks, Jurassic plutons and associated dikes, Miocene volcanic rocks, and a variety of Quaternary surficial deposits derived from local bedrock units. The purpose of the project was to map the area in detail, with primary emphasis on the pre-Quaternary units, to provide an improved stratigraphic, structural, and geochronologic framework for use in land management applications and scientific research.

California↗

USGS2019

Geologic map of the Poncha Pass area, Chaffee, Fremont, and Saguache Counties, Colorado

This report presents a 1:24,000-scale geologic map, cross sections, and descriptive and interpretative text for the Poncha Pass area in central Colorado. The map area is irregular in shape, covering all of one 7 ½' quadrangle (Poncha Pass) and parts of five others (Mount Ouray, Maysville, Salida West, Salida East, and Wellsville). The map boundaries were drawn to cover all of the “Poncha mountain block,” our designation for the approximately 15-kilometer-long northwestern end of the Sangre de Cristo Mountains. The map conveys the areal distribution of (1) Proterozoic basement rocks forming the core of the Poncha mountain block, (2) overlying Eocene and Oligocene volcanic rocks, (3) Miocene and younger basin-fill deposits, (4) Quaternary surficial glacial and alluvial deposits, and (5) faults and folds affecting all of the above units. The Poncha mountain block, which lies within the Rio Grande rift, is topographically and geologically distinctive. Generally, the Rio Grande rift is internally characterized by subsided structural basins or grabens and subdued, low-relief topography rather than elevated mountain blocks. The intrarift, topographically high Poncha mountain block spans the axial part of the rift and separates the low-lying basins of the west-tilted upper Arkansas River half graben and east-northeast-tilted San Luis half graben. These distinctive aspects of the Poncha mountain block were the primary motivations to conduct geologic mapping in the area. Important questions addressed by geologic mapping and related studies in the Poncha Pass area include (1) what were the structural controls and tectonic mechanism(s) that resulted in development of the Poncha mountain block in an intrarift environment; (2) did surface uplift of the Poncha block occur during rift development in the Neogene and Quaternary, and at what rate(s); (3) how was extensional strain accommodated and relayed across the Poncha block between the opposite-polarity rift basins and flanking mountain blocks; (4) is there a clear Laramide deformational signal in rocks of the map area; and (5) have earlier Laramide contractional structures, if they exist, influenced later rift-related extensional deformation through reactivation or strain localization. Prior to our mapping, the geology of much of the Poncha Pass area had only been mapped in reconnaissance fashion, reflecting the area’s poor bedrock exposures, poor access due to the rugged terrain, and geologic complexity. The map presented here provides new details of the geology of this difficult area and helps elucidate the development of the Poncha block and improves understanding of the geologic framework and geologic history of the area.

Colorado↗

USGS2026

Geologic map of the Emmons Lake volcanic center, Alaska

Introduction The Emmons Lake volcanic center is a spatially clustered group of stratovolcanoes and calderas in the southwestern part of the Alaska Peninsula, Alaska. The volcanic center is characterized by several ice- and snow-clad stratovolcanoes located within and along the margins of a nested-caldera complex that includes Emmons Lake. A shieldlike ancestral edifice (ancestral Mount Emmons) is truncated by the caldera complex and forms a broad volcanic platform around the center. The main stratovolcanoes of the Emmons Lake volcanic center are Pavlof Sister, Pavlof Volcano, Little Pavlof, Double Crater, Mount Hague, and Mount Emmons. Several small unnamed cinder cones and vents also are located within Emmons Lake volcanic center and on the east flank of Pavlof Volcano. Many of these cones and vents have been the source of the young lava flows that mantle the floor of the caldera. Pavlof Volcano, in the northeastern part of the Emmons Lake volcanic center, is one of the most historically (that is, the past about 300 years) active volcanoes in Alaska, and eruptions from Pavlof Volcano pose the greatest hazards to the region. Volcanic rocks of the Emmons Lake volcanic center overlie continental and marine sedimentary rocks of chiefly Late Jurassic to early Tertiary age. The oldest rocks in the area are those of the Naknek Formation, consisting of volcaniclastic sandstone, siltstone, and conglomerate of Late Jurassic age. The southern part of the area includes rocks of the Belkofski Formation, a thick sequence of volcaniclastic sandstone, siltstone, and conglomerate of middle Tertiary age. Lava flows, volcanic breccia, and fluvial volcaniclastic rocks of late Miocene age, which unconformably overlie the Belkofski Formation south of the Emmons Lake volcanic center, are primarily exposed on the islands just south of the Alaska Peninsula. The Emmons Lake volcanic center was affected multiple times by glaciation associated with the glacier expansion that characterized the Quaternary. Glaciation has played a key role in shaping the present-day landscape, and much of the eruptive history of the Emmons Lake volcanic center has involved interactions with glacier ice. Thus, a brief review of the Quaternary glacial history of the area is provided to establish the physical context for Emmons Lake volcanic center eruptive activity.

Alaska↗

USGS2004

Simulation of ground-water flow in the Cedar River alluvial aquifer flow system, Cedar Rapids, Iowa

The Cedar River alluvial aquifer is the primary source of municipal water in the Cedar Rapids, Iowa, area. Since 1992, the U.S. Geological Survey, in cooperation with the City of Cedar Rapids, has investigated the hydrogeology and water quality of the Cedar River alluvial aquifer. This report describes a detailed analysis of the ground-water flow system in the alluvial aquifer, particularly near well field areas. The ground-water flow system in the Cedar Rapids area consists of two main components, the unconsolidated Quaternary deposits and the underlying carbonate bedrock that has a variable fracture density. Quaternary deposits consist of eolian sand, loess, alluvium, and glacial till. Devonian and Silurian bedrock aquifers overlie the Maquoketa Shale (Formation) of Ordovician age, a regional confining unit. Ground-water and surface-water data were collected during the study to better define the hydrogeology of the Cedar River alluvial aquifer and Devonian and Silurian aquifers. Stream stage and discharge, ground-water levels, and estimates of aquifer hydraulic properties were used to develop a conceptual ground-water flow model and to construct and calibrate a model of the flow system. This model was used to quantify the movement of water between the various components of the aluvial aquifer flow system and provide an improved understanding of the hydrology of the alluvial aquifer. Ground-water flow was simulated for the Cedar River alluvial aquifer and the Devonian and Silurian aquifers using the three-dimensional finite-difference ground-water flow model MODFLOW. The model was discretized into 223 rows and 354 columns of cells. Areal cell sizes range from about 50 feet on a side near the Cedar River and the Cedar Rapids municipal wells to 1,500 feet on a side near the model boundaries and farthest away from the Cedar Rapids municipal well fields. The model is separated into five layers to account for the various hydrogeologic units in the model area. Model results indicate that the primary sources of inflow to the modeled area are infiltration from the Cedar River (53.0 percent) and regional flow in the glacial and bedrock materials (34.1 percent). The primary sources of outflow from the modeled area are discharge to the Cedar River (45.4 percent) and pumpage (44.8 percent). Current steady-state pumping rates have increased the flow of water from the Cedar River to the alluvial aquifer by 43.8 cubic feet per second. Steady-state and transient hypothetical pumpage scenarios were used to show the relation between changes in pumpage and changes in infiltration of water from the Cedar River. Results indicate that more than 99 percent of the water discharging from municipal wells infiltrates from the Cedar River, that the time required for induced river recharge to equilibrate with municipal pumpage may be 150 days or more, and that ground-water availability in the Cedar Rapids area will not be significantly affected by doubling current pumpage as long as there is sufficient flow in the Cedar River to provide recharge.

Iowa↗

USGS2004

Thrust-induced collapse of mountains-an example from the "Big Bend" region of the San Andreas Fault, western Transverse Ranges, California

Mount Pinos and Frazier Mountain are two prominent mountains just south of the San Andreas fault in the western Transverse Ranges of southern California, a region that has undergone rapid Quaternary contraction and uplift. Both mountains are underlain, at least in part, by thrusts that place granitic and gneissic rocks over sedimentary rocks as young as Pliocene. Broad profiles and nearly flat summits of each mountain have previously been interpreted as relicts of a raised erosion surface. However, several features bring this interpretation into question. First, lag or stream gravels do not mantle the summit surfaces. Second, extensive landslide deposits, mostly pre?Holocene and deeply incised, mantle the flanks of both mountains. Third, a pervasive fracture and crushed?rock network pervades the crystalline rocks underlying both mountains. The orientation of the fractures, prominent in roadcuts on Mount Pinos, is essentially random. 'Hill?and?saddle' morphology characterizes ridges radiating from the summits, especially on Mount Pinos; outcrops are sparse on the hills and are nonexistent in the saddles, suggesting fractures are concentrated in the saddles. Latest movement on the thrusts underlying the two mountain massifs is probably early Quaternary, during which the mountains were uplifted to considerably higher (although unknown) elevations than at present. A model proposes that during thrusting, ground accelerations in the hanging wall, particularly near thrust tips, were high enough to pervasively fracture the hanging?wall rocks, thereby weakening them and producing essentially an assemblage of loose blocks. Movement over flexures in the fault surface accentuated fracturing. The lowered shear stresses necessary for failure, coupled with deep dissection and ongoing seismic activity, reduced gravitational potential by spreading the mountain massifs, triggering flanking landslides and producing broad, flat?topped mountains. This study developed from mapping in the western Transverse Ranges as part of the U.S. Geological Survey's Southern California Areal Mapping Project (SCAMP).

California↗

USGS2006

Hydrogeology and simulation of ground-water flow in the Silurian-Devonian aquifer system, Johnson County, Iowa

Bedrock of Silurian and Devonian age (termed the “Silurian-Devonian aquifer system”) is the primary source of ground water for Johnson County in east-central Iowa. Population growth within municipal and suburban areas of the county has resulted in increased amounts of water withdrawn from this aquifer and water-level declines in some areas. A 3-year study of the hydrogeology of the Silurian-Devonian aquifer system in Johnson County was undertaken to provide a quantitative assessment of ground water resources and to construct a ground-water flow model that can be used by local governmental agencies as a management tool. Johnson County is underlain by unconsolidated deposits of Quaternary age and Paleozoic-age bedrock units. The bulk of the Quaternary deposits consists of weathered and unweathered glacial till; however, shallow alluvium and buried sand and gravel deposits also are present. Six bedrock hydrogeologic units are present in Johnson County (oldest to youngest): Maquoketa confining unit, Silurian aquifer, Wapsipinicon Group (aquifer and confining unit), Cedar Valley aquifer, Upper Devonian shale confining unit, and Cherokee confining unit. Although separate aquifers and confining units are described, the Silurian- and Devonian-age units are considered as a single aquifer system. The top of the Silurian-Devonian aquifer system is considered as the top of the Cedar Valley aquifer, where present, and the base of the aquifer system is considered as the top of the Maquoketa confining unit. The hydraulic properties of the rocks that comprise the Silurian-Devonian aquifer system are highly variable as a result of the variable composition of the rocks and the presence of solution features in some of the carbonate-rock units. For the combined Silurian-Devonian aquifer system, specific capacity averages 2.1 gallons per minute per foot of drawdown, transmissivity averages about 580 feet squared per day, and hydraulic conductivity averages 8.3 feet per day. Recharge to the Silurian-Devonian aquifer system in Johnson County is predominantly from infiltration of precipitation to the bedrock. Discharge from the aquifer is primarily to municipal, industrial, and private-development wells. Reliable measurements of the amount of recharge to or discharge from the ground-water system in Johnson County, however, are not available. Altitude of the 1996 potentiometric surface ranged from more than 750 feet above the North American Vertical Datum of 1988 (NAVD88) in northern Johnson County to less than 575 feet above NAVD88 in the central part of the county. A large cone of depression within the potentiometric surface is present in the central part of the county, between Coralville and Iowa City. A large limestone quarry is located near the center of this cone of depression. Ground water generally flows from the northern and western parts of Johnson County either toward the cone of depression in the center of the county or south out of the county. Ground water also flows toward the Cedar River in the northeastern part of the county. A ground-water divide in the northeastern part of the county roughly approximates the surface-water divide between the Iowa River and Cedar River drainages. A numerical ground-water-flow model of the Silurian-Devonian aquifer system in Johnson County was used to test concepts of ground-water flow, to assess the need for additional data, and to evaluate the potential effects of anticipated increased ground-water development and drought. The 1-layer model was calibrated to average 1996 ground-water conditions, which were assumed to approximate steady-state flow conditions. The model also was used to simulate steady-state conditions for 2004, steady-state conditions using anticipated pumping rates for 2025, and potential future drought conditions. The simulated potentiometric surface generally replicated the potentiometric surface for 1996 and 2004 conditions. The calculated root mean squared error values for the 1996 and 2004 simulations were 13.6 and 18.6 feet, respectively. The mean absolute differences between measured and simulated water levels for the 1996 and 2004 simulations were about 11 and 14 feet, respectively. Total model-calculated inflow to the ground-water system for the 1996 simulation was 19.6 million gallons per day (Mgal/d), and the largest model-calculated inflow component was areal recharge (15.1 Mgal/d). Total model-calculated outflow from the ground-water system was 19.7 Mgal/d, and the largest outflow component was discharge to wells (10.5 Mgal/d). Model-calculated water-budget components for the 2004 simulation were similar to the 1996 components. Potential future steady-state conditions were simulated using anticipated 2025 pumping rates. Pumpage both for existing wells and for assumed new wells, based on anticipated population growth in the northern part of the county and for the nearby municipalities, was included in the model. Simulated 2025 pumpage was about 1.5 Mgal/d greater than simulated 2004 pumpage. Simulated steady-state ground-water levels, using anticipated 2025 pumping rates, were lower than 2004 simulated levels throughout the county, and simulated water-level declines ranged from less than 1 foot near the county boundaries to about 11 feet. Potential future drought conditions were simulated by assuming that recharge to the Silurian-Devonian aquifer system is reduced by a factor of 0.75 and that water-supply pumpage is increased by a factor of 1.25 over the anticipated 2025 pumping rates. Overall, simulated water levels for future drought conditions were greater than 5 feet lower than simulated 2004 conditions and were a maximum of about 30 feet lower in the northeastern part of the county. The greatest limitation to the model is the lack of measured or estimated water-budget components for comparison to simulated water-budget components. Because the model is only calibrated to measured water levels, and not to water-budget components, the model results are nonunique. Other model limitations include the relatively coarse grid scale, lack of detailed information on pumpage from the quarry and from private developments and domestic wells, and the lack of separate water-level data for the Silurian- and Devonian-age rocks.

Iowa↗

USGS2006

Geohydrology and water chemistry of the Alexander Valley, Sonoma County, California

This study of the geohydrology and water chemistry of the Alexander Valley, California, was done to provide an improved scientific basis for addressing emerging water-management issues, including potential increases in water demand and changes in flows in the Russian River. The study tasks included (1) evaluation of existing geohydrological, geophysical, and geochemical data; (2) collection and analysis of new geohydrologic data, including subsurface lithologic data, ground-water levels, and streamflow records; and (3) collection and analysis of new water-chemistry data. The estimated total water use for the Alexander Valley for 1999 was approximately 15,800 acre-feet. About 13,500 acre-feet of this amount was for agricultural use, primarily vineyards, and about 2,300 acre-feet was for municipal/industrial use. Ground water is the main source of water supply for this area. The main sources of ground water in the Alexander Valley are the Quaternary alluvial deposits, the Glen Ellen Formation, and the Sonoma Volcanics. The alluvial units, where sufficiently thick and saturated, comprise the best aquifer in the study area. Average recharge to the Alexander Valley is estimated from a simple, basinwide water budget. On the basis of an estimated annual average of 298,000 acre-feet of precipitation, 160,000 acre-feet of runoff, and 113,000 to 133,000 acre-feet of evapotranspiration, about 5,000 to 25,000 acre-feet per year is available for ground-water recharge. Because this estimate is based on differences between large numbers, there is significant uncertainty in this recharge estimate. Long-term changes in ground-water levels are evident in parts of the study area, but because of the sparse network and lack of data on well construction and lithology, it is uncertain if any significant changes have occurred in the northern part of the study area since 1980. In the southern half of the study area, ground-water levels generally were lower at the end of the 2002 irrigation season than at the end of the 1980 season, which suggests that a greater amount of ground water is being pumped in the southern half of the study area in recent years compared with that pumped in the early 1980s. Water-chemistry data for samples collected from 11 wells during 2002-04 indicate that water quality in the study area generally is acceptable for potable use. Two wells, however, each contained one constituent (241 ?g/L of manganese and 1,350 ?g/L of boron) in excess of the recommended standards for drinking water (50 ?g/L and 1,000 ?g/L, respectively). The chemical composition of water from most of the wells sampled for major ions plot as a mixed cation-bicarbonate, magnesium-bicarbonate, or calcium-bicarbonate type water. The ionic composition of the historical and recent samples from wells in the Alexander Valley is similar to that of the historical surface-water samples collected from the Russian River near Healdsburg. This suggests a similar source of water, particularly for wells that are less than 200 feet total depth and perforated in Quaternary alluvial deposits. Water from deeper, non-alluvial wells may contain slightly higher concentrations of sodium as a result of cation exchange. Water samples collected from several wells over an approximately 30-year period suggest a progressive change in water chemistry over time. Samples from the southern part of the valley show a trend towards higher ionic concentrations and increasing concentrations of particular constituents such as sulfate.

California↗

USGS2006

Geology, Water, and Wind in the Lower Helmand Basin, Southern Afghanistan

This report presents an overview of the geology, hydrology, and climate of the lower Helmand Basin, a large, closed, arid basin in southern Afghanistan. The basin is drained by the Helmand River, the only perennial desert stream between the Indus and Tigris-Euphrates Rivers. The Helmand River is the lifeblood of southern Afghanistan and has supported desert civilizations in the Sistan depression for over 6,000 years. The Helmand Basin is a structurally closed basin that began to form during the middle Tertiary as a consequence of the collision of several Gondwanaland fragments. Aeromagnetic studies indicate the basin is 3-5 kilometers deep over basement rocks. Continued subsidence along basin-bounding faults in Iran and Pakistan throughout the Neogene has formed the Sistan depression in the southwest corner of the basin. Lacustrine, eolian, and fluvial deposits are commonly exposed in the basin and were intruded by latest Miocene-middle Quaternary volcanoes, which indicates that depositional environments in the lower Helmand Basin have not substantially changed for nearly 10 million years. Lakes expanded in the Sistan depression during the Quaternary; however, the size and extent of these pluvial lakes are unknown. Climate conditions in the lower Helmand Basin likely mirrored climate changes in the Rajasthan Desert to the east and in Middle Eastern deserts to the west: greater aridity during global episodes of colder temperatures and increased available moisture during episodes of warmer temperatures. Eolian processes are unusually dominant in shaping the landscape in the basin. A strong wind blows for 120 days each summer, scouring dry lakebeds and creating dune fields from annual flood deposits. Nearly one-third of the basin is mantled with active or stabilized dunes. Blowing winds combined with summer temperatures over 50? Celsius and voluminous insect populations hatched from the deltaic wetlands create an environment referred to as the 'most odious place on earth' by 19th century visitors. During dry years, large plumes of dust originating from Sistan are recorded by weather satellites. The Helmand River drains about 40 percent of Afghanistan and receives most of its moisture from melting snow and spring storms. Similar to many desert streams, the Helmand and its main tributary, the Arghandab River, are characterized by large fluctuations in monthly and annual discharges. Water from the Helmand accumulates in several hamuns (shallow lakes) in the Sistan depression. The wetlands surrounding these hamuns are the largest in western Asia and are directly affected by droughts and floods on the Helmand. Average annual discharge on the Helmand is about 6.12 million megaliters (million cubic meters), and the annual discharge varies by a factor of five. In 2005, the region was just beginning to recover from the longest drought (1998-2005) of record back to 1830. Annual peak discharges range from less than 80 cubic meters per second in 1971 to nearly 19,000 cubic meters per second in 1885. Large floods fill each hamun to overflowing to create one large lake that overflows into the normally dry Gaud-i Zirreh basin. The interaction of flooding, active subsidence, and wind erosion causes frequent channel changes on the Helmand delta. A major development effort on the Helmand River was initiated after World War II with substantial aid from the United States. Two dams and several major canals were completed in the 1950s; however, poor drainage conditions on the newly prepared agricultural fields caused extensive waterlogging and salinization. New drains were installed and improved agricultural methods were implemented in the 1970s, and some lands became more productive. Since 1980, Afghanistan has endured almost constant war and civil and political strife. In 2005, the country was on a path to rebuild much of its technical infrastructure. Revitalization of agricultural lands in the lower Helmand Basin and improved managem

Scientific Investigations Report↗

USGS2008

Substrate Geochemistry and Soil Development in Boreal Forest and Tundra Ecosystems in the Yukon-Tanana Upland and Seward Peninsula, Alaska

We report on soil development as a function of bedrock type and the presence of loess in two high latitude ecosystems (boreal forest and tundra) and from two regions in Alaska?the Yukon-Tanana Upland (YTU, east-central Alaska) and the Seward Peninsula (SP, far-west coastal Alaska). This approach to the study of 'cold soils' is fundamental to the quantification of regional geochemical landscape patterns. Of the five state factors in this study, bedrock and biota (ecosystem; vegetation zone) vary whereas climate (within each area) and topography are controlled. The influence of time is assumed to be controlled, as these soils are thousands of years old (late Quaternary to Holocene). The primary minerals in soils from YTU, developed over loess and crystalline bedrock (metamorphic and intrusive), are quartz, plagioclase, and 2:1 clays; whereas in the SP, where loess and metasedimentary bedrock (schist and quartzite) predominate, they are quartz and muscovite. The A horizon of both regions is rich in peat. Examination of the ratio of mobile (K2O, CaO, and Fe2O3) to immobile (TiO2) major oxides, within each region, shows that very little difference exists in the chemical weathering of soils developed between the two ecosystems examined. Differences were observed between tundra soils developed in the two regions. These differences are most probably due to the dissimilarity in the geochemical importance of both loess and bedrock. A minimal loss of cadmium with soil depth is seen for soils developed over YTU crystalline bedrock in the boreal forest environments. This trend is related to the mobility of cadmium in these soils as well as to its biogenic cycling. Major differences were observed in the proportion of cadmium and zinc among the A, B, and C horizon material sequestered in various soil fractions as measured by sequential soil extractions. These trends followed such variables as the decrease with depth in organic matter, the change in clay minerals, and the change in the proportion of oxides/hydroxides. An analysis of the bulk soil mineralogy and the relation between CaO and MgO and Al2O3 and Fe2O3 indicates that the silty textured soils of the YTU are predominantly eolian (that is, of late Tertiary or Quaternary age) but not broad-regional in origin. Their composition instead is probably the result of locally derived dusts as well as input from long-term, in-place bedrock weathering.

Scientific Investigations Report↗

USGS2008

Statistical Stationarity of Sediment Interbed Thicknesses in a Basalt Aquifer, Idaho National Laboratory, Eastern Snake River Plain, Idaho

The statistical stationarity of distributions of sedimentary interbed thicknesses within the southwestern part of the Idaho National Laboratory (INL) was evaluated within the stratigraphic framework of Quaternary sediments and basalts at the INL site, eastern Snake River Plain, Idaho. The thicknesses of 122 sedimentary interbeds observed in 11 coreholes were documented from lithologic logs and independently inferred from natural-gamma logs. Lithologic information was grouped into composite time-stratigraphic units based on correlations with existing composite-unit stratigraphy near these holes. The assignment of lithologic units to an existing chronostratigraphy on the basis of nearby composite stratigraphic units may introduce error where correlations with nearby holes are ambiguous or the distance between holes is great, but we consider this the best technique for grouping stratigraphic information in this geologic environment at this time. Nonparametric tests of similarity were used to evaluate temporal and spatial stationarity in the distributions of sediment thickness. The following statistical tests were applied to the data: (1) the Kolmogorov-Smirnov (K-S) two-sample test to compare distribution shape, (2) the Mann-Whitney (M-W) test for similarity of two medians, (3) the Kruskal-Wallis (K-W) test for similarity of multiple medians, and (4) Levene's (L) test for the similarity of two variances. Results of these analyses corroborate previous work that concluded the thickness distributions of Quaternary sedimentary interbeds are locally stationary in space and time. The data set used in this study was relatively small, so the results presented should be considered preliminary, pending incorporation of data from more coreholes. Statistical tests also demonstrated that natural-gamma logs consistently fail to detect interbeds less than about 2-3 ft thick, although these interbeds are observable in lithologic logs. This should be taken into consideration when modeling aquifer lithology or hydraulic properties based on lithology.

Scientific Investigations Report↗

USGS2009

Hydrogeologic Framework of the Yakima River Basin Aquifer System, Washington

The Yakima River basin aquifer system underlies about 6,200 square miles in south-central Washington. The aquifer system consists of basin-fill deposits occurring in six structural-sedimentary basins, the Columbia River Basalt Group (CRBG), and generally older bedrock. The basin-fill deposits were divided into 19 hydrogeologic units, the CRBG was divided into three units separated by two interbed units, and the bedrock was divided into four units (the Paleozoic, the Mesozoic, the Tertiary, and the Quaternary bedrock units). The thickness of the basin-fill units and the depth to the top of each unit and interbed of the CRBG were mapped. Only the surficial extent of the bedrock units was mapped due to insufficient data. Average mapped thickness of the different units ranged from 10 to 600 feet. Lateral hydraulic conductivity (Kh) of the units varies widely indicating the heterogeneity of the aquifer system. Average or effective Kh values of the water-producing zones of the basin-fill units are on the order of 1 to 800 ft/d and are about 1 to 10 ft/d for the CRBG units as a whole. Effective or average Kh values for the different rock types of the Paleozoic, Mesozoic, and Tertiary units appear to be about 0.0001 to 3 ft/d. The more permeable Quaternary bedrock unit may have Kh values that range from 1 to 7,000 ft/d. Vertical hydraulic conductivity (Kv) of the units is largely unknown. Kv values have been estimated to range from about 0.009 to 2 ft/d for the basin-fill units and Kv values for the clay-to-shale parts of the units may be as small as 10-10 to 10-7 ft/d. Reported Kv values for the CRBG units ranged from 4x10-7 to 4 ft/d. Variations in the concentrations of geochemical solutes and the concentrations and ratios of the isotopes of hydrogen, oxygen, and carbon in groundwater provided information on the hydrogeologic framework and groundwater movement. Stable isotope ratios of water (deuterium and oxygen-18) indicated dispersed sources of groundwater recharge to the CRBG and basin-fill units and that the source of surface and groundwater is derived from atmospheric precipitation. The concentrations of dissolved methane were larger than could be attributable to atmospheric sources in more than 80 percent of wells with measured methane concentrations. The concentrations of the stable isotope of carbon-13 of methane were indicative of a thermogenic source of methane. Most of the occurrences of methane were at locations several miles distant from mapped structural fault features, suggesting the upward vertical movement of thermogenic methane from the underlying bedrock may be more widespread than previously assumed or there may be a more general occurrence of unmapped (buried) fault structures. Carbon and tritium isotope data and the concentrations of dissolved constituents indicate a complex groundwater flow system with multiple contributing zones to groundwater wells and relative groundwater residence time on the order of a few tens to many thousands of years. Potential mean annual recharge for water years 1950-2003 was estimated to be about 15.6 in. or 7,149 ft3/s (5.2 million acre-ft) and includes affects of human activities such as irrigation of croplands. If there had been no human activities (predevelopment conditions) during that time period, estimated recharge would have been about 11.9 in. or 5,450 ft3/s (3.9 million acre-ft). Estimated mean annual recharge ranges from virtually zero in the dry parts of the lower basin to more than 100 in. in the humid uplands, where annual precipitation is more than 120 in. Groundwater in the different hydrogeologic units occurs under perched, unconfined, semiconfined, and confined conditions. Groundwater moves from topographic highs in the uplands to topographic low areas along the streams. The flow system in the basin-fill units is compartmentalized due to topography and geologic structure. The flow system also is compartmentalized for the CRBG units but not to as large

Scientific Investigations Report↗

USGS2010

Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model area in Michigan and Indiana. Overall water use for the selected areas in Wisconsin, Michigan, Indiana, and Illinois was less during early time intervals than during more recent intervals, with large increases beginning around the 1960s. Total estimated groundwater withdrawals for model input range from 18.01 million gallons per day (Mgal/d) for interval 1 (1864-1900) to 1,280.25 Mgal/d for interval 12 (2001-5). Withdrawals for the public-supply category make up the majority of the withdrawals in each of the four states. In Wisconsin and Michigan, the second largest withdrawals are for the irrigation category; in Indiana and Illinois, industrial withdrawals account for the second largest withdrawal amounts. The smallest withdrawals are for miscellaneous uses in Wisconsin and irrigation uses in Indiana and Illinois. Estimated groundwater withdrawals in the Southern Lower Peninsula of Michigan, Northeastern Illinois, and the farfield model area are generally larger than in the other model subregions. Withdrawals in Michigan and Indiana are predominantly from the Quaternary aquifer system, whereas withdrawals in Illinois are predominantly from the Cambrian-Ordovician aquifer systems. Withdrawals in Wisconsin are about equal from the Quaternary and Cambrian-Ordovician aquifer systems. Estimated groundwater withdrawals in Michigan and Indiana are predominantly from the unconfined unconsolidated aquifer type. Withdrawals in Illinois are largely from the deep confined bedrock aquifer type, although they decreased considerably in more recent time intervals. Wisconsin withdrawals are about equal from unconfined unconsolidated and deep confined bedrock aquifer types. Groundwater-withdrawal estimates in Wisconsin were compiled for the 47 easternmost counties within the boundary of the Lake Michigan Basin model, of which 32 counties, though not entirely contained, are at least partly within the Lake Michigan Basin. Overall, 6,457 withdrawal locations were estima

Scientific Investigations Report↗

USGS2010

Characterization of geologic deposits in the vicinity of US Ecology, Amargosa Basin, southern Nevada

Multiple approaches have been applied to better understand the characteristics of geologic units exposed at the surface and buried at depth in the vicinity of US Ecology (USE), a low-level commercial waste site in the northern Amargosa Desert, Nevada. Techniques include surficial geologic mapping and interpretation of the subsurface using borehole data. Dated deposits at depth were used to estimate rates of sediment accumulation. The subsurface lithologies have been modeled in three dimensions. Lithologic cross sections have been created from the three-dimensional model and have been compared to resistivity data at the same location. Where deposits appear offset, a fault was suspected. Global Positioning System elevation transects were measured and trenches were excavated to locate a strand of the Carrara Fault. The presence of the fault helps to better understand the shape of the potentiometric surface. These data will be used to better understand the hydrologic parameters controlling the containment of the waste at US Ecology. Quaternary geologic units exposed at the surface, in the vicinity of US Ecology, are derived from the alluvium shed off the adjacent range front and the Amargosa River. These deposits vary from modern to early Pleistocene in age. At depth, heterogeneous sands and gravel occur. Observed in deep trenches and boreholes, the subsurface deposits are characterized as fining-upward sequence of sediment from 5- to 8-meters thick. No volcanic units or fine-grained playa deposits were described in the boreholes to a depth of 200 meters. Based on Infrared Stimulated Luminescence dated core samples, short-term rates of sediment accumulation (<70,000 years) are an average of 2.7 millimeters per year, however, long-term rates (<3,900,000 years) are orders of magnitude less. Resistivity data, when compared to lithologic cross sections, generally are consistent with lithology grain size and probable soil carbonate accumulations. Surface resistivity displays a fining-upward sequence of sediments at the surface with a soil carbonate imprint. Finally, trenching north of US Ecology successfully exposed offset Quaternary deposits on a splay of the Carrara Fault. Holocene deposits do not appear to be faulted, however, a fault zone does intersect middle and late Pleistocene aged units.

Nevada↗

USGS2012

Geologic framework, regional aquifer properties (1940s-2009), and spring, creek, and seep properties (2009-10) of the upper San Mateo Creek Basin near Mount Taylor, New Mexico

The U.S. Geological Survey, in cooperation with the U.S. Forest Service, examined the geologic framework, regional aquifer properties, and spring, creek, and seep properties of the upper San Mateo Creek Basin near Mount Taylor, which contains areas proposed for exploratory drilling and possible uranium mining on U.S. Forest Service land. The geologic structure of the region was formed from uplift of the Zuni Mountains during the Laramide Orogeny and the Neogene volcanism associated with the Mount Taylor Volcanic Field. Within this structural context, numerous aquifers are present in various Paleozoic and Mesozoic sedimentary formations and the Quaternary alluvium. The distribution of the aquifers is spatially variable because of the dip of the formations and erosion that produced the current landscape configuration where older formations have been exhumed closer to the Zuni Mountains. Many of the alluvial deposits and formations that contain groundwater likely are hydraulically connected because of the solid-matrix properties, such as substantive porosity, but shale layers such as those found in the Mancos Formation and Chinle Group likely restrict vertical flow. Existing water-level data indicate topologically downgradient flow in the Quaternary alluvium and indiscernible general flow patterns in the lower aquifers. According to previously published material and the geologic structure of the aquifers, the flow direction in the lower aquifers likely is in the opposite direction compared to the alluvium aquifer. Groundwater within the Chinle Group is known to be confined, which may allow upward migration of water into the Morrison Formation; however, confining layers within the Chinle Group likely retard upward leakage. Groundwater was sodium-bicarbonate/sulfate dominant or mixed cation-mixed anion with some calcium/bicarbonate water in the study area. The presence of the reduction/oxidation-sensitive elements iron and manganese in groundwater indicates reducing conditions at some time or in some location(s) in most aquifers. Frequent detections of zinc in the alluvium aquifer may represent anthropogenic influences such as mining. Along the mesas in the upper San Mateo Creek Basin, springs that form various creeks, including El Rito and San Mateo Creeks, discharge from the basalt-cap layer and the upper Cretaceous sedimentary layers. Streamflow in El Rito and San Mateo Creeks flows down steep gradients near the mesas sustained by groundwater discharges, and this streamflow transitions to shallow groundwater contained within the valley alluvium through infiltration where the subsequent groundwater is restricted from downward migration by the shaly Menefee Formation. This shallow groundwater reemerges at seeps where the land surface has been eroded below the groundwater level. Spring- and creek-water samples contained small amounts of dissolved solutes, and seep water contained substantially larger amounts of dissolved solutes. The pH of water within the creeks was neutral to alkaline, and all locations exhibited well-oxygenated conditions, although typically at substantially less than saturated levels. Changes in the stable-isotope ratios of water between spring and summer samples indicate differences in source-water inputs that likely pertain to seasonal recharge sources. Results of the water-isotope analysis and geochemical modeling indicate little evaporation and chemical weathering at the spring and creek sites but stronger evaporation and chemical weathering by the time the water reaches the seep locations in the center of the upper San Mateo Creek Basin.

New Mexico↗

USGS2013

Groundwater and surface-water interaction within the upper Smith River Watershed, Montana 2006-2010

The 125-mile long Smith River, a tributary of the Missouri River, is highly valued as an agricultural resource and for its many recreational uses. During a drought starting in about 1999, streamflow was insufficient to meet all of the irrigation demands, much less maintain streamflow needed for boating and viable fish habitat. In 2006, the U.S. Geological Survey, in cooperation with the Meagher County Conservation District, initiated a multi-year hydrologic investigation of the Smith River watershed. This investigation was designed to increase understanding of the water resources of the upper Smith River watershed and develop a detailed description of groundwater and surface-water interactions. A combination of methods, including miscellaneous and continuous groundwater-level, stream-stage, water-temperature, and streamflow monitoring was used to assess the hydrologic system and the spatial and temporal variability of groundwater and surface-water interactions. Collectively, data are in agreement and show: (1) the hydraulic connectedness of groundwater and surface water, (2) the presence of both losing and gaining stream reaches, (3) dynamic changes in direction and magnitude of water flow between the stream and groundwater with time, (4) the effects of local flood irrigation on groundwater levels and gradients in the watershed, and (5) evidence and timing of irrigation return flows to area streams. Groundwater flow within the alluvium and older (Tertiary) basin-fill sediments generally followed land-surface topography from the uplands to the axis of alluvial valleys of the Smith River and its tributaries. Groundwater levels were typically highest in the monitoring wells located within and adjacent to streams in late spring or early summer, likely affected by recharge from snowmelt and local precipitation, leakage from losing streams and canals, and recharge from local flood irrigation. The effects of flood irrigation resulted in increased hydraulic gradients (increased groundwater levels relative to stream stage) or even reversed gradient direction at several monitoring sites coincident with the onset of nearby flood irrigation. Groundwater-level declines in mid-summer were due to groundwater withdrawals and reduced recharge from decreased precipitation, increased evapotranspiration, and reduced leakage in some area streams during periods of low flow. Groundwater levels typically rebounded in late summer, a result of decreased evapotranspiration, decreased groundwater use for irrigation, increased flow in losing streams, and the onset of late-season flood irrigation at some sites. The effect of groundwater and surface-water interactions is most apparent along the North and South Forks of the Smith River where the magnitude of streamflow losses and gains can be greater than the magnitude of flow within the stream. Net gains consistently occurred over the lower 15 miles of the South Fork Smith River. A monitoring site near the mouth of the South Fork Smith River gained (flow from the groundwater to the stream) during all seasons, with head gradients towards the stream. Two upstream sites on the South Fork Smith River exhibited variable conditions that ranged from gaining during the spring, losing (flowing from the stream to the groundwater) during most of the summer as groundwater levels declined, and then approached or returned to gaining conditions in late summer. Parts of the South Fork Smith River became dry during periods of losing conditions, thus classifying this tributary as intermittent. The North Fork Smith River is highly managed at times through reservoir releases. The North Fork Smith River was perennial throughout the study period although irrigation diversions removed a large percentage of streamflow at times and losing conditions persisted along a lower reach. The lowermost reach of the North Fork Smith River near its mouth transitioned from a losing reach to a gaining reach throughout the study period. Groundwater and surface-water interactions occur downstream from the confluence of the North and South Fork Smith Rivers, but are less discernible compared to the overall magnitude of the main-stem streamflow. The Smith River was perennial throughout the study. Monitoring sites along the Smith River generally displayed small head gradients between the stream and the groundwater, while one site consistently showed strongly gaining conditions. Synoptic streamflow measurements during periods of limited irrigation diversion in 2007 and 2008 consistently showed gains over the upper 41.4 river miles of the main stem Smith River where net gains ranged from 13.0 to 28.9 cubic feet per second. Continuous streamflow data indicated net groundwater discharge and small-scale tributary inflow contributions of around 25 cubic feet per second along the upper 10-mile reach of the Smith River for most of the 2010 record. A period of intense irrigation withdrawal during the last two weeks in May was followed by a period (early June 2010 to mid-July 2010) with the largest net increase (an average of 71.1 cubic feet per second) in streamflow along this reach of the Smith River. This observation is likely due to increased groundwater discharge to the Smith River resulting from irrigation return flow. By late July, the apparent effects of return flows receded, and the net increase in streamflow returned to about 25 cubic feet per second. Two-dimensional heat and solute transport VS2DH models representing selected stream cross sections were used to constrain the hydraulic properties of the Quaternary alluvium and estimate temporal water-flux values through model boundaries. Hydraulic conductivity of the Quaternary alluvium of the modeled sections ranged from 3x10-6 to 4x10-5 feet per second. The models showed reasonable approximations of the streambed and shallow aquifer environment, and the dynamic changes in water flux between the stream and the groundwater through different model boundaries.

Montana↗

USGS2014

Aquifers of Arkansas: protection, management, and hydrologic and geochemical characteristics of groundwater resources in Arkansas

Sixteen aquifers in Arkansas that currently serve or have served as sources of water supply are described with respect to existing groundwater protection and management programs, geology, hydrologic characteristics, water use, water levels, deductive analysis, projections of hydrologic conditions, and water quality. State and Federal protection and management programs are described according to regulatory oversight, management strategies, and ambient groundwater-monitoring programs that currently (2013) are in place for assessing and protecting groundwater resources throughout the State. Physical attributes, groundwater geochemistry, and groundwater quality are described for each of the 16 aquifers of the State. Information in regard to the hydrology and geochemistry of each of the aquifers is summarized from about 550 historical and recent publications. Additionally, more than 8,000 sites with groundwater-quality data were obtained from the U.S. Geological Survey National Water Information System and the Arkansas Department of Environmental Quality databases and entered into a spatial database to investigate distribution and trends in chemical constituents for each of the aquifers. The 16 aquifers of the State were divided into two major physiographic regions of the State: the Coastal Plain Province (referred to as Coastal Plain) of eastern and southern Arkansas, which includes 11 of the 16 aquifers, and the Interior Highlands Division (referred to as Interior Highlands) of western Arkansas, which includes the remaining 5 aquifers. The 11 aquifers in the Coastal Plain consist of various geologic units that are Cenozoic in age and consist primarily of Cretaceous, Tertiary, and Quaternary sands, gravels, silts, and clays. Groundwater in the Coastal Plain represents one of the most valuable natural resources in the State, driving the economic engines of agriculture, while also supplying abundant water for commercial, industrial, and public-supply use. In terms of age from youngest to oldest, the aquifers of the Coastal Plain include Quaternary alluvial aquifers, including the Mississippi River Valley alluvial aquifer (the most important aquifer in Arkansas in terms of volume of use and economic benefits), the Jackson Group (a regional confining unit that served for decades as an important source of domestic supply), and the Cockfield, Sparta, Cane River, Carrizo, Wilcox, Nacatoch, Ozan, Tokio, and Trinity aquifers. The Mississippi River Valley alluvial aquifer accounts for approximately 94 percent of all groundwater used in the State, and the aquifer is used primarily for irrigation purposes. The Sparta aquifer is the second most important aquifer in terms of use, and the aquifer was used in the past dominantly as a source of public and industrial supply, although increasing irrigation use is occurring because of critically declining water levels in the Mississippi River Valley alluvial aquifer. Other aquifers of the Coastal Plain generally are used as important local sources of domestic, industrial, and public supply, in addition to other minor uses. Water quality generally is good for all aquifers of the Coastal Plain, except for elevated iron concentrations and localized areas of high salinity. The high salinity results from intrusion from underlying formations, evapotranspiration processes in areas of low recharge, and inadequate flushing in downgradient areas of residual salinity from deposition in marine environments. Trends in the spatial distribution of individual chemical constituents are related to position along the flow path for most aquifers of the Coastal Plain. These trends include elevated iron and nitrate concentrations with lower pH values and dissolved solids in groundwater from the outcrop areas, transitioning to lower iron and nitrate (related to changes in redox) and higher pH and dissolved solids (dominantly from the dissolution of carbonate minerals) in groundwater downgradient from outcrop areas. Groundwater generally trended from a calcium- to a sodium-bicarbonate water type with increasing cation exchange along the flow path. The Interior Highlands of western Arkansas has less reported groundwater use than other areas of the State, reflecting a combination of factors. These factors include prevalent and increasing use of surface water, less intensive agricultural uses, lower population and industry densities, lesser potential yield of the resource, and lack of detailed reporting. The overall low yields of aquifers of the Interior Highlands result in domestic supply as the dominant use, with minor industrial, public, and commercial-supply use. Where greater volumes are required for growth of population and industry, surface water is the greatest supplier of water needs in the Interior Highlands. The various aquifers of the Interior Highlands generally occur in shallow, fractured, well-indurated, structurally modified bedrock of this mountainous region of the State, as compared to the relatively flat-lying, unconsolidated sediments of the Coastal Plain. In terms of age from youngest to oldest, the aquifers of the Interior Highlands include: the Arkansas River Valley alluvial aquifer, the Ouachita Mountains aquifer, the Western Interior Plains confining system, the Springfield Plateau aquifer, and the Ozark aquifer. Spatial trends in groundwater geochemistry in the Interior Highlands differ greatly from trends noted for aquifers of the Coastal Plain. In the Coastal Plain, the prevalence of long regional flow paths results in regionally predictable and mappable geochemical changes along the flow paths. In the Interior Highlands, short, topographically controlled flow paths (from hilltops to valleys) within small watersheds represent the predominant groundwater-flow system. As such, dense data coverage from numerous wells would be required to effectively characterize these groundwater basins and define small-scale geochemical changes along any given flow path for aquifers of the Interior Highlands. Changes in geochemistry generally were related to rock type and residence time along individual flow paths. Dominant changes in geochemistry for the Ouachita Mountains aquifer and the Western Interior Plains confining system are attributed to rock/water interaction and changes in redox zonation along the flow path. In these areas, groundwater evolves along flow paths from a calcium- to a sodium-bicarbonate water type with increasing reducing conditions resulting in denitrification, elevated iron and manganese concentrations, and production of methane in the more geochemically evolved and strongest reducing conditions. In the Ozark and Springfield Plateau aquifers, rapid influx of surface-derived contaminants, especially nitrogen, coupled with few to no attenuation processes was attributed to the karst landscape developed on Mississippian- and Ordovician-age carbonate rocks of the Ozark Plateaus. Increasing nitrate concentrations are related to increasing agricultural land use, and areas of mature karst development result in higher nitrate concentrations than areas with less karst features.

Arkasas↗

USGS2016

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗