Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Clear Waters”

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,045 records · Page 58Linked to original sources

Downstream effects of dams on channel geometry and bottomland vegetation: Regional patterns in the Great Plains

The response of rivers and riparian forests to upstream dams shows a regional pattern related to physiographic and climatic factors that influence channel geometry. We carried out a spatial analysis of the response of channel geometry to 35 dams in the Great Plains and Central Lowlands, USA. The principal response of a braided channel to an upstream dam is channel-narrowing, and the principal response of a meandering channel is a reduction in channel migration rate. Prior to water management, braided channels were most common in the southwestern Plains where sand is abundant, whereas meandering channels were most common in the northern and eastern Plains. The dominant response to upstream dams has been channel-narrowing in the southwestern Plains (e.g., six of nine cases in the High Plains) and reduction in migration rate in the north and east (e.g., all of twelve cases in the Missouri Plateau and Western Lake Regions). Channel-narrowing is associated with a burst of establishment of native and exotic woody riparian pioneer species on the former channel bed. In contrast, reduction in channel migration rate is associated with a decrease in reproduction of woody riparian pioneers. Thus, riparian pioneer forests along large rivers in the southwestern Plains have temporarily increased following dam construction while such forests in the north and east have decreased. These patterns explain apparent contradictions in conclusions of studies that focused on single rivers or small regions and provide a framework for predicting effects of dams on large rivers in the Great Plains and elsewhere. These conclusions are valid only for large rivers. A spatial analysis of channel width along 286 streams ranging in mean annual discharge from 0.004 to 1370 cubic meters per second did not produce the same clear regional pattern, in part because the channel geometries of small and large streams are affected differently by a sandy watershed.

Great Plains↗

Optimizing selection of training and auxiliary data for operational land cover classification for the LCMAP initiative

The U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) initiative is a new end-to-end capability to continuously track and characterize changes in land cover, use, and condition to better support research and applications relevant to resource management and environmental change. Among the LCMAP product suite are annual land cover maps that will be available to the public. This paper describes an approach to optimize the selection of training and auxiliary data for deriving the thematic land cover maps based on all available clear observations from Landsats 4–8. Training data were selected from map products of the U.S. Geological Survey’s Land Cover Trends project. The Random Forest classifier was applied for different classification scenarios based on the Continuous Change Detection and Classification (CCDC) algorithm. We found that extracting training data proportionally to the occurrence of land cover classes was superior to an equal distribution of training data per class, and suggest using a total of 20,000 training pixels to classify an area about the size of a Landsat scene. The problem of unbalanced training data was alleviated by extracting a minimum of 600 training pixels and a maximum of 8000 training pixels per class. We additionally explored removing outliers contained within the training data based on their spectral and spatial criteria, but observed no significant improvement in classification results. We also tested the importance of different types of auxiliary data that were available for the conterminous United States, including: (a) five variables used by the National Land Cover Database, (b) three variables from the cloud screening ‘‘Function of mask” (Fmask) statistics, and (c) two variables from the change detection results of CCDC. We found that auxiliary variables such as a Digital Elevation Model and its derivatives (aspect, position index, and slope), potential wetland index, water probability, snow probability, and cloud probability improved the accuracy of land cover classification. Compared to the original strategy of the CCDC algorithm (500 pixels per class), the use of the optimal strategy improved the classification accuracies substantially (15-percentage point increase in overall accuracy and 4-percentage point increase in minimum accuracy).

ISPRS Journal of Photogrammetry and Remote Sensing↗

Geologic Map of the Saint Helens Quadrangle, Columbia County, Oregon, and Clark and Cowlitz Counties, Washington

The Saint Helens 7.5' quadrangle is situated in the Puget-Willamette Lowland approximately 35 km north Portland, Oregon. The lowland, which extends from Puget Sound into west-central Oregon, is a complex structural and topographic trough that lies between the Coast Range and the Cascade Range. Since late Eocene time, Cascade Range has been the locus of a discontinuously active volcanic arc associated with underthrusting of oceanic lithosphere beneath the North American continent along the Cascadia Subduction Zone. The Coast Range occupies the forearc position within the Cascadia arc-trench system and consists of a complex assemblage of Eocene to Miocene volcanic and marine sedimentary rocks. The Saint Helens quadrangle lies in the northern part of the Portland Basin, a roughly 2000-km2 topographic and structural depression. It is the northernmost of several sediment-filled structural basins that collectively constitute the Willamette Valley segment of the Puget-Willamette Lowland (Beeson and others, 1989; Swanson and others, 1993; Yeats and others, 1996). The rhomboidal basin is approximately 70 km long and 30 km wide, with its long dimension oriented northwest. The Columbia River flows west and north through the Portland Basin at an elevation near sea level and exits through a confined bedrock valley less than 2.5 km wide about 16 km north of Saint Helens. The flanks of the basin consist of Eocene through Miocene volcanic and sedimentary rocks that rise to elevations exceeding 2000 ft (610 m). Seismic-reflection profiles (L.M. Liberty, written commun., 2003) and lithologic logs of water wells (Swanson and others, 1993; Mabey and Madin, 1995) indicate that as much as 550 m of late Miocene and younger sediments have accumulated in the deepest part of the basin near Vancouver. Most of this basin-fill material was carried in from the east by the Columbia River but contributions from streams draining the adjacent highlands are locally important. The Portland Basin has been interpreted as a pull-apart basin located in the releasing stepover between two echelon, northwest-striking, right-lateral fault zones (Beeson and others, 1985, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). These fault zones are thought to reflect regional transpression and dextral shear within the forearc in response to oblique subduction along the Cascadia Subduction Zone Pezzopane and Weldon, 1993; Wells and others, 1998). The southwestern margin of the Portland Basin is a well-defined topographic break along the base of the Tualatin Mountains, an asymmetric anticlinal ridge that is bounded its northeast flank by the Portland Hills Fault Zone (Balsillie and Benson, 1971; Beeson and others, 1989; Blakely and others, 1995), which is probably an active structure (Wong and others, 2001; Liberty and others, 2003). The nature of the corresponding northeastern margin of the basin is less clear, but a poorly defined and partially buried dextral extensional fault zone has been hypothesized from topography, microseismicity, potential fieldanomalies, and reconnaissance geologic mapping (Beeson and others, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). Another dextral structure, the Kalama Structural Zone of Evarts (2002), may underlie the north-northwest-trending reach of the Columbia River north of Woodland (Blakely and others, 1995). This map is a contribution to a U.S. Geological Survey (USGS) program designed to improve the geologic database for the Portland Basin region of the Pacific Northwest urban corridor, the populated forearc region of western Washington and Oregon. Better and more detailed information on the bedrock and surficial geology of the basin and its surrounding area is needed to refine assessments of seismic risk (Yelin and Patton, 1991; Bott and Wong, 1993), ground-failure hazards (Madin and Wang, 1999; Wegmann and Walsh, 2001) and resource availability in this rapid

Scientific Investigations Map↗

Reinterpreting the importance of oxygen-based biodegradation in chloroethene-contaminated groundwater

Chlororespiration is common in shallow aquifer systems under conditions nominally identified as anoxic. Consequently, chlororespiration is a key component of remediation at many chloroethene-contaminated sites. In some instances, limited accumulation of reductive dechlorination daughter products is interpreted as evidence that natural attenuation is not adequate for site remediation. This conclusion is justified when evidence for parent compound (tetrachloroethene, PCE, or trichloroethene, TCE) degradation is lacking. For many chloroethene-contaminated shallow aquifer systems, however, nonconservative losses of the parent compounds are clear but the mass balance between parent compound attenuation and accumulation of reductive dechlorination daughter products is incomplete. Incomplete mass balance indicates a failure to account for important contaminant attenuation mechanisms and is consistent with contaminant degradation to nondiagnostic mineralization products like CO 2 . While anoxic mineralization of chloroethene compounds has been proposed previously, recent results suggest that oxygen-based mineralization of chloroethenes also can be significant at dissolved oxygen concentrations below the currently accepted field standard for nominally anoxic conditions. Thus, reassessment of the role and potential importance of low concentrations of oxygen in chloroethene biodegradation are needed, because mischaracterization of operant biodegradation processes can lead to expensive and ineffective remedial actions. A modified interpretive framework is provided for assessing the potential for chloroethene biodegradation under different redox conditions and the probable role of oxygen in chloroethene biodegradation.

Ground Water Monitoring and Remediation↗

Geology and ground-water resources of the island of Kauai, Hawaii

Kauai is one of the oldest, and is structurally the most complicated, of the Hawaiian Islands. Like the others, it consists principally of a huge shield volcano, built up from the sea floor by many thousands of thin flows of basaltic lava. The volume of the Kauai shield was on the order of 1,000 cubic miles. Through much of its growth it must have resembled rather closely the presently active shield volcano Mauna Loa, on the island of Hawaii. When the Kauai volcano started its growth is not known with certainty, but it is believed that activity started late in the Tertiary period, possibly in the early or middle part of the Pliocene epoch. Growth of the shield was rapid and probably was completed before the end of the Pliocene. Toward the end of the growth of the shield, its summit collapsed to form a broad caldera, the largest that has been found in the Hawaiian Islands. Like the calderas of Kilauea and Mauna Loa, that of Kauai volcano had boundaries that were, in part, rather indefinite. The principal depression was bordered by less depressed fault blocks, some of which merged imperceptibly with the outer slopes of the volcano. Elsewhere the caldera rim was low, and flows spilled over it onto the outer slopes. The well-defined central depression of the Kauai caldera was approximately 10 to 12 miles across. At about the same time as the formation of the major caldera, another, smaller caldera was formed by collapse around a minor eruptive center on the southeastern side of the Kauai shield. Lavas accumulated in the calderas, gradually filling them and burying banks of talus that formed along the foot of the boundary cliffs. The caldera-filling lavas differed from those that built the major portion of the shield in being much thicker and more massive as a result of ponding in the depressions. The petrographic types for the most part are the same throughout. Both the flank flows that built most of the shield and the flows that filled the calderas are predominantly olivine basalt. Picrite-basalt (oceanite), containing very abundant large phenocrysts of olivine, and basalt containing little or no olivine are present but together comprise less than 10 percent of the whole. Late in the period of filling of the major caldera a small amount of basaltic andesine andesite was extruded. Near the end of the period of filling of the major caldera further collapse occurred, forming a large graben on the southwestern side of the shield. Lava flows erupting within the caldera poured southwestward over the cliff bounding the graben and spread over the gently sloping graben floor. Near the present Waimea Canyon their advance was obstructed by the fault scarp at the west edge of the graben. The cliff along the northeast edge of the graben eventually was buried by lava flows from within the caldera, but that along the west edge continued to stand above the level of the flows in the graben. The flows that accumulated in the graben are of the same types as those that filled the caldera, and like them are mostly thick and massive because of ponding by the graben walls and of the gentle slopes of the graben floor over which they spread. The rocks of the major Kauai shield volcano are known as the Waimea Canyon volcanic series. The thin flows that accumulated on the flanks of the shield, which compose the major portion of the volcanic edifice, are named the Napali formation of the Waimea Canyon volcanic series. The rocks that accumulated in the big summit caldera are named the Olokele formation, and those that filled the small caldera on the southeast flank of the shield are named the Haupu formation. The volcanic rocks accumulated in the graben on the southwestern side of the shield are named the Makaweli formation of the Waimea Canyon volcanic series, and sedimentary rocks interbedded with them are known as the Mokuone member of the Makaweli formation. Few vents of the Waimea Canyon volcanic series have been recognized, probably because most of them have been destroyed by erosion or are buried by later lavas. Large numbers of dikes cut the lavas of the Napali formation along Waimea Canyon and the Napali Coast and along the east edge of the Waialeale massif. Fewer dikes are found in the other members of the series. Some tendency toward radial arrangement of the dikes is present, but the dominant trend all over the island is east-northeastward. Another great collapse took place on the eastern flank of the volcano at about the time the major shield became extinct, or shortly afterward. A subcircular graben 6 or 7 miles across sank several thousand feet, forming a broad depression between the Waialeale massif on the west and Kalepa and Nonou ridges on the east. This collapsed structure cannot be as clearly demonstrated as the Makaweli graben on the southwest side of the shield, because its walls have been greatly eroded and its floor is deeply buried by lavas of the later Koloa volcanic series. It appears, however, to be the only reasonable explanation of the physiography of the eastern side of the island. After the completion of the great Kauai shield came a long period of erosion during which no volcanic activity occurred. Waves cut high sea cliffs around the island, and streams cut canyons as much as 3,000 feet deep. Thick soil formed over much of the mountain. Then volcanism was renewed. Eruption occurred from a series of minor vents arranged in nearly north-south and northeast-southwest lines across the eastern two-thirds of the island. The lavas, cinder cones, and ash beds of this period of volcanism are known as the Koloa volcanic series. Lavas of the Koloa volcanic series include olivine basalt, picrite-basalt (mimosite) with few phenocrysts of olivine, basanite, nepheline basalt, melilite-nepheline basalt, and ankaratrite (nepheline basalt very rich in pyroxene and olivine). Inclusions of dunite, composed almost entirely of olivine, are common in flows of the Koloa. Just before and during the eruption of the Koloa volcanic series, voluminous landslides and mudflows brought down a large amount of rock debris and soil from the steep slopes of the mountainous central upland and deposited it as breccias at the foot of the steep slopes in valley heads and along the border of the marginal lowland. Streams distributed part of the material across the lowland. The breccias and conglomerates thus formed, and later buried by lavas of the Koloa volcanic series, are named the Palikea formation of the Koloa volcanic series. The structures formed at Koloa vents include cinder cones, one tuff cone, and lava cones. The latter are miniature shields resembling the major shield volcano, formed by repeated outpourings of fluid lava. The tuff cone, at the west side of Kilauea Bay, was formed by phreatomagmatic explosions caused by rising magma coming in contact with water-saturated rocks. Volcanism during Koloa time continued for a long period but was not continuous over the entire area. Locally, long periods of quiet occurred, allowing streams to re-excavate some of the canyons filled by earlier flows of the Koloa volcanic series, and weathering to form soils later buried by new flows. Some of the canyons thus formed during the time when the Koloa was being deposited were several hundred feet deep. Volcanism probably continued throughout most of the Pleistocene epoch. The latest flow of the Koloa volcanic series appears very recent, and rests on lithified calcareous dunes formed during one of the Pleistocene low stands of the sea. During the Pleistocene epoch stream valleys and sea cliffs were eroded to base levels governed by one or more stands of the sea more than 100 feet below present sea level. Beaches of calcareous sand were formed, and the sand blown inland to form calcareous dunes, now lithified. A test boring near Moloaa penetrated calcareous sand 160 feet below sea level, at the foot of a high sea cliff. Coral reef also was built around part or all of the island, and in part buried by lavas of the Koloa volcanic series. The explosions that built the tuff cone at Kilauea Bay threw up fragments of limestone from a buried reef. Much of the apron of lavas of the Kalna series around the northeastern side of the island probably rests on a platform formed below present sea level by wave erosion and the growth of coral reef. As the sea rose around the island, the valley mouths were alluviated. Several levels of the sea higher than the present one probably are represented. Some stream terraces may be graded to a stand of the sea as high as 260 feet above present sea level, but no positive evidence for stands higher than 25 feet have been found. Well-preserved shorelines are recognized approximately 25 and 5 feet above sea level. Much of the present coral reef appears to have been formed when the sea stood about 5 feet higher than now, and reduced to its present level by solutional weathering and wave erosion. The lavas of the Napali formation of the Waimea Canyon volcanic series are highly permeable. They carry basal water over much of the island, and yield it freely to wells. This water is fresh everywhere except very close to the coast on the leeward side of the island. In some areas they may contain water confined at high levels between dikes. The lavas of the Olokele and Haupu formations are moderately to poorly permeable. They probably contain fresh water at sea level, but would not yield it readily to wells. Locally, ash beds perch small bodies of fresh water at high levels in the lavas of the Olokele formation, but these are of no economic importance. The lavas of the Makaweli formation also arc moderately to poorly permeable. They carry fresh or brackish water at sea level. In general, they yield water to wells less readily than the lavas of the Napali formation, but more readily than the lavas of the Olokele. The conglomerates and breccias of the Mokuone member are poorly permeable, but are not known to perch more than a slight amount of water in the overlying lavas, The lava flows of the Koloa volcanic series are poorly to moderately permeable. They carry fresh or brackish water at sea level, but generally yield it slowly to wells. Locally, small bodies of fresh water are perched at high levels in the lavas of the Koloa by beds of ash and soil and by breccia and conglomerate of the Palikea formation. Both the older and the younger alluvium generally are poorly permeable, but contain small amounts of fresh or brackish water. The lithified calcareous dunes are permeable, but they appear to contain only brackish water. Lagoon deposits on the Mana plain are poorly to moderately permeable and yield brackish water to wells.

Hawaii↗

Landscape level effects of invasive plants and animals on water infiltration through Hawaiian tropical forests

Watershed degradation due to invasion threatens downstream water flows and associated ecosystem services. While this topic has been studied across landscapes that have undergone invasive-driven state changes (e.g., native forest to invaded grassland), it is less well understood in ecosystems experiencing within-system invasion (e.g. native forest to invaded forest). To address this subject, we conducted an integrated ecological and ecohydrological study in tropical forests impacted by invasive plants and animals. We measured soil infiltration capacity in multiple fenced (i.e., ungulate-free)/unfenced and native/invaded forest site pairs along moisture and substrate age gradients across Hawaii to explore the effects of invasion on hydrological processes within tropical forests. We also characterized forest composition, structure and soil characteristics at these sites to assess the direct and vegetation-mediated impacts of invasive species on infiltration capacity. Our models show that invasive ungulates negatively affect soil infiltration capacity consistently across the wide moisture and substrate age gradients considered. Additionally, several soil characteristics known to be affected by invasive ungulates were associated with local infiltration rates, indicating that the long-term secondary effects of high ungulate densities in tropical forests may be stronger than effects observed in this study. The effect of invasive plants on infiltration was complex and likely to depend on their physiognomy within existing forest community structure. These results provide clear evidence for managers that invasive ungulate control efforts can improve ecohydrological function of mesic and wet forest systems critical to protecting downstream and nearshore resources and maintaining groundwater recharge.

Hawaii↗

DRIFteRS: A dataset of drift invertebrate densities in streams and rivers across western North America, 1997–2024

Prey availability is among the most influential and highly variable determinants of fish growth and freshwater habitat carrying capacity, yet it remains understudied compared to physical habitat variables (Ouellet et al., 2025; Rosenfeld et al., 2014; Weber et al., 2017). We often lack a clear understanding of how much food is available to fishes, how it varies spatially and temporally, and how it influences responses to restoration (Ouellet et al., 2025; Rossi et al., 2024; Wipfli et al., 2010). Drift invertebrates—the primary food source for juvenile salmonids and other drift-foraging fishes—play a pivotal role in these dynamics. To better understand the spatiotemporal variability of drift invertebrate abundance and biomass across the freshwater range of drift-feeding salmonids in western North America, we compiled the DRift Invertebrates For salmonids in River Systems (DRIFteRS) dataset. The dataset encompasses 6125 samples of drift invertebrates, and, for a subset of drift samples, associated benthic invertebrate density data, collected from 1360 reaches on 459 unique rivers and streams spanning 55 river basins considered hydrologically independent (i.e., not nested within the same larger watershed) across British Columbia, Canada, and the U.S. states of Alaska, Arizona, California, Colorado, Idaho, Nevada, New Mexico, Oregon, Utah, Washington, and Wyoming. Sample sites represent a diverse array of river and stream habitats (e.g., headwater, mainstem, side channel), in watersheds with diverse land uses (e.g., urban, wilderness, agricultural), and disturbance histories (e.g., fire, restoration). Collected between 1997 and 2024, the data span the full calendar year and capture daily and seasonal patterns in drift abundance and biomass densities. When paired with water quality and quantity data as well as remotely sensed environmental landscape data, such as land use/land cover, climate, and disturbance history, channel morphology, and riparian vegetation composition, the DRIFteRS dataset can aid in identifying key drivers of drift invertebrate densities and mean body size and support predictive modeling in unsampled locations and times. The dataset may also be used to analyze aquatic-terrestrial resource flows, derive prey-encounter rates and profitability (mean prey size), and inform broader investigations of sit-and-wait foraging ecology, especially when paired with data on drift-foraging predators. For salmonid-focused applications, the dataset can be integrated into habitat evaluation models, including bioenergetic (e.g., Naman et al., 2019) and life cycle models (e.g., Beechie et al., 2023), to improve estimates of habitat capacity and population dynamics for river- and stream-rearing salmonids. Understanding prey availability dynamics is increasingly important, because rising water temperatures increase salmonid metabolic demands (Crozier et al., 2010). Flow regime transitions (i.e., snow or glacier dominated to rain dominated; Beechie et al., 2013), wildfire frequency and intensity (Hessburg et al., 2021), as well as plant community and phenology shifts (Cleland et al., 2007; Franklin et al., 2016) are all predicted to change with rising temperatures and are potential drivers of terrestrial and aquatic invertebrate prey quality and availability in lotic systems. These insights can ultimately inform restoration prioritization and design, helping managers consider food resource implications when evaluating restoration priorities and match habitat improvement to food supply. The data are released under the Creative Commons Attribution 4.0 International license and may be reused with attribution.

Alaska, Arizona, British Columbia, California, Col↗

Changes in streamflow and the flux of nutrients in the Mississippi-Atchafalaya River Basin, USA, 1980-2007

Nutrients and freshwater delivered by the Mississippi and Atchafalaya Rivers drive algal production in the northern Gulf of Mexico, which eventually results in the widespread occurrence of hypoxic bottom waters along the Louisiana and Texas coast. Researchers have demonstrated a relation between the extent of the hypoxic zone and the magnitude of streamflow, nutrient fluxes, and nutrient concentrations in the Mississippi River, with springtime streamflows and fluxes being the most predictive. In 1999 the U.S. Geological Survey (USGS) estimated the flux of nitrogen, phosphorus, and silica at selected sites in the Mississippi Basin and to the Gulf of Mexico for 1980-1996. These flux estimates provided the baseline information used by the Mississippi River/Gulf of Mexico Watershed Nutrient Task Force to develop an Action Plan for reducing hypoxia in the northern Gulf of Mexico. The primary goal of the Action Plan was to achieve a reduction in the size (areal extent) of the hypoxic zone from an average of approximately 14,000 square kilometers in 1996-2000 to a 5-year moving average of less than 5,000 square kilometers by 2015. Improved statistical models and adjusted maximum likelihood estimation using USGS Load Estimator (LOADEST) software were used to estimate annual and seasonal nutrient fluxes for 1980-2007 at selected sites on the Mississippi River and its tributaries. These data provide a means to evaluate the influence of natural and anthropogenic effects on delivery of water and nutrients to the Gulf of Mexico; to define subbasins that are the most important contributors of nutrients to the gulf; and to investigate the relations among streamflow, nutrient fluxes, and the size and duration of the Gulf of Mexico hypoxic zone. A comparative analysis between the baseline period of 1980-1996 and 5-year moving averages thereafter indicate that the average annual streamflow and fluxes of total nitrogen, nitrate, orthophosphate, and silica to the Gulf of Mexico have decreased. However, the flux of total phosphorus between the baseline period and subsequent 5-year periods has increased. The average spring (April, May, and June) streamflow and fluxes of silica, total nitrogen, nitrate, and orthophosphate to the Gulf of Mexico also decreased, whereas the spring flux of total phosphorus has increased. Similar changes in streamflow and nutrient flux were observed at many sites Buxtonwithin the basin. The inputs of water, total nitrogen, and total phosphorus from the major subbasins of the Mississippi-Atchafalaya River Basin as a percentage of the to-the-gulf totals have increased from the Ohio River Basin, decreased from the Missouri River Basin, and remained relatively unchanged from the Upper Mississippi, Red, and Arkansas River Basins. Changes in streamflow and nutrient fluxes are related, but short-term variations in sources of streamflow and nutrients complicate the interpretation of factors that affect nutrient delivery to the Gulf of Mexico. Parametric time-series models are used to try and separate natural variability in nutrient flux from changes due to other causes. Results indicate that the decrease in annual nutrient fluxes that has occurred between the 1980-1996 baseline period and more recent years can be largely attributed to natural causes (climate and streamflow) and not management actions or other human controlled activities in the Mississippi-Atchafalaya River Basin. The downward trends in total nitrogen, nitrate, ammonium, and orthophosphate that were detected at either the Mississippi River near St. Francisville, La., or the Atchafalaya River at Melville, La., occurred prior to 1995. In spite of the general decrease in nutrient flux, the average size of the Gulf of Mexico hypoxic zone has increased between 1997 and 2007. The reasons for this are not clear but could be due to the type or nature of nutrient delivery. Whereas the annual flux of total nitrogen to the Gulf of Mexico has decreased, the proporti

Scientific Investigations Report↗

Geologic Map of the Woodland Quadrangle, Clark and Cowlitz Counties, Washington

The Woodland 7.5' quadrangle is situated in the Puget-Willamette Lowland approximately 50 km north of Portland, Oregon (fig. 1). The lowland, which extends from Puget Sound into west-central Oregon, is a complex structural and topographic trough that lies between the Coast Range and the Cascade Range. Since late Eocene time, the Cascade Range has been the locus of an active volcanic arc associated with underthrusting of oceanic lithosphere beneath the North American continent along the Cascadia Subduction Zone. The Coast Range occupies the forearc position within the Cascadia arc-trench system and consists of a complex assemblage of Eocene to Miocene volcanic and marine sedimentary rocks. The Woodland quadrangle lies at the northern edge of the Portland Basin, a roughly 2000-km2 topographic and structural depression that is the northernmost of several sediment-filled structural basins, which collectively constitute the Willamette Valley segment of the Puget-Willamette Lowland (Beeson and others, 1989; Swanson and others, 1993; Yeats and others, 1996). The Portland Basin is approximately 70 km long and 30 km wide; its long dimension is oriented northwest. Its northern boundary coincides, in part, with the lower Lewis River, which flows westward through the center of the quadrangle. The Lewis drains a large area in the southern Washington Cascade Range, including the southern flank of Mount St. Helens approximately 25 km upstream from the quadrangle, and joins the Columbia River about 6 km south of Woodland (fig. 1). Northwest of Woodland, the Columbia River exits the broad floodplain of the Portland Basin and flows northward through a relatively narrow bedrock valley at an elevation near sea level. The flanks of the Portland Basin consist of Eocene through Miocene volcanic and sedimentary rocks that rise to elevations exceeding 2000 ft (610 m). Seismic-reflection profiles (L.M. Liberty, written commun., 2003) and lithologic logs of water wells (Swanson and others, 1993; Mabey and Madin, 1995) indicate that as much as 550 m of late Miocene and younger sediments have accumulated in the deepest part of the basin near Vancouver. Most of this basin-fill material was carried in from the east by the Columbia River but sediment deposited by streams draining the adjacent highlands are locally important. The Portland Basin has been interpreted as a pull-apart basin located in the releasing stepover between two en echelon, northwest-striking, right-lateral fault zones (Beeson and others, 1985, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). These fault zones are thought to reflect regional transpression and dextral shear within the forearc in response to oblique subduction of the Pacific Plate along the Cascadia Subduction Zone (Pezzopane and Weldon, 1993; Wells and others, 1998). The southwestern margin of the Portland Basin is a well-defined topographic break along the base of the Tualatin Mountains, an asymmetric anticlinal ridge that is bounded on its northeast flank by the Portland Hills Fault Zone (Balsillie and Benson, 1971; Beeson and others, 1989; Blakely and others, 1995), which is probably an active structure (Wong and others, 2001; Liberty and others, 2003). The nature of the corresponding northeastern margin of the basin is less clear, but a poorly defined and partially buried dextral extensional fault zone has been hypothesized from topography, microseismicity, potential field-anomalies, and reconnaissance geologic mapping (Beeson and others, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). Another dextral structure may control the north-northwest-trending reach of the Columbia River between Portland and Longview (Blakely and others, 1995; Evarts, 2002; Evarts and others, 2002). This map is a contribution to a U.S. Geological Survey program designed to improve the geologic database for the Portland Basin part of the Pacific Northwest urban corridor,

Scientific Investigations Map↗

Water-quality and geochemical variability in the Little Arkansas River and Equus aquifer, south-central Kansas, 2001–16

The city of Wichita’s water supply currently (2019) comes from two primary sources: Cheney Reservoir and the Equus Beds aquifer. The Equus Beds aquifer storage and recovery project was developed to help the city of Wichita meet increasing future water demands. Source water for artificial recharge comes from the Little Arkansas River during above-base-flow conditions, is treated using National Primary Drinking Water Regulations as a guideline, and is injected into the Equus Beds aquifer through recharge wells or surface spreading basins for later use. The Equus Beds aquifer storage and recovery project currently (2019) consists of two coexisting phases. Phase I began in 2007 and captures Little Arkansas River water and indirect streambank diversion well water for aquifer recharge using 4 wells and 2 recharge basins. Phase II began in 2013 and currently (2019) includes a surface-water treatment facility, a river intake facility, eight recharge injection wells, and a third recharge basin. The U.S. Geological Survey, in cooperation with the City of Wichita, completed this study to summarize water-quality and geochemical variability of the Equus Beds aquifer. Data in this report can be used to establish baseline conditions before implementing artificial aquifer recharge further, document groundwater quality, evaluate changing conditions, identify environmental factors affecting groundwater, provide science-based information for decision making, and help meet regulatory monitoring requirements. Physicochemical properties were measured and water-quality data were collected from 2 Little Arkansas River surface-water sites and 63 Equus Beds aquifer groundwater sites, including 38 areal assessment index wells (IWs) during 2001 through 2016. Data collection included discrete samples and additional continuous measurements at selected sites. Discretely collected samples were analyzed for physicochemical properties, dissolved solids, primary ions, nutrients (nitrogen and phosphorus species), organic carbon, indicator bacteria, trace elements, arsenic species, organic compounds, and radioactivity. This report focuses discussion on aquifer water quality. Federal drinking-water criteria were used to evaluate aquifer water quality. Primary drinking-water criteria are those that are enforceable for public drinking water. Secondary criteria are those that can cause aesthetics or tastes that are unpleasant. Continuously collected data at a subset of sites included streamflow, groundwater levels, water temperature, specific conductance, pH, oxidation-reduction potential (ORP), dissolved oxygen, turbidity, nitrate plus nitrite, and fluorescent dissolved organic matter. Continuous measurement of physicochemical properties in near-real time allowed characterization of Little Arkansas River surface water and Equus Beds aquifer groundwater during conditions and time scales that would not have been possible otherwise and served as a complement to discrete water-quality sampling. During 2001 through 2016, less than 1 percent of chloride and nitrate plus nitrite, 7 percent of dissolved iron, 48 percent of dissolved manganese, 12 percent of dissolved arsenic, and 39 percent of atrazine detections in surface-water samples exceeded their respective Federal primary or secondary drinking-water criteria. None of the surface-water samples collected exceeded the Federal sulfate criterion, and every sample had detections of total coliform bacteria during the study. Constituents of concern in the Equus Beds aquifer exceeded their respective Federal criteria throughout the study period and included chloride, sulfate, nitrate plus nitrite, Escherichia coli ( E. coli ), total coliforms, and dissolved iron and arsenic species. About 5 percent of shallow (less than 80 feet) and 7 percent of deep (greater than 80 feet) IW chloride sample concentrations exceeded the secondary Federal criterion of 250 milligrams per liter (mg/L). Chloride tended to exceed its criterion in shallow and deep wells along the Arkansas River and near Burrton, Kansas, an area with past oil and gas activities. Chloride concentrations near Burrton were larger in the deep parts of the aquifer. About 18 percent of shallow and 13 percent of deep IW sulfate sample concentrations exceeded the secondary Federal criterion of 250 mg/L. Mean sulfate concentrations tended to exceed the criterion in the central part of the study area. Shallow IW mean nitrate plus nitrite (hereafter referred to as “nitrate”) was substantially larger than mean deep IW nitrate. Geochemical conditions in the deeper aquifer reduced forms of nitrogen to species such as ammonia. About 15 percent of shallow and less than 1 percent of deep IW nitrate sample concentrations exceeded the Federal criterion of 10 mg/L. Mean shallow IW nitrate concentrations exceeded the criterion in the northeastern and southeastern parts of the study area; on average, deep IW nitrate concentrations did not exceed the criterion. E. coli and fecal coliform bacteria detections were usually at or near the detection limit. E. coli was detected in 3 percent of shallow and deep IWs, and fecal coliform bacteria were detected in 8 percent of shallow and 6 percent of deep IWs. Total coliforms were detected in 24 percent of shallow and 12 percent of deep IWs. E. coli coliphage was detected in two shallow IW samples (1 percent of samples) at the detection limit and was not detected in deep IW samples. Dissolved iron was detected in 51 percent of shallow and 62 percent of deep IW samples. Dissolved iron concentrations exceeded the secondary Federal criterion of 0.3 mg/L in 38 percent of shallow and 46 percent of deep IW samples. Mean dissolved iron concentrations were largest mostly in the central and northwest part of the study area corresponding to an area of the aquifer where aquifer material is more clay-rich. The distribution of large dissolved iron concentrations was similar to that of large sulfate concentrations. About 55 percent of shallow and 92 percent of deep IW dissolved manganese samples exceeded the secondary Federal criterion of 0.05 mg/L. Almost all samples from the central and northern parts of the study area had mean dissolved manganese concentrations that exceeded the Federal criterion in the shallow part of the aquifer. Mean dissolved manganese concentrations in the shallow part of the aquifer were substantially large (greater than 1,000 micrograms per liter [μg/L]) in wells near the Little Arkansas River and in the central part of the study area because of chemically reducing conditions in the aquifer that likely related to larger percentages of clay in the aquifer material. Concentrations of dissolved arsenic species generally were larger in the deep parts of the aquifer. Arsenite was the dominant form of arsenic on average in shallow (52 percent) and deep (55 percent) IWs. About 12 percent of shallow and 34 percent of deep IW dissolved arsenic sample concentrations exceeded the Federal primary drinking criterion of 10 μg/L. Shallow IW dissolved arsenic concentrations were larger near the Little Arkansas River and the center of the study area; large shallow IW dissolved arsenic concentrations (10–50 μg/L) in the center of the study area correspond to areas that have had the most water-level recovery since the historical low in 1993. Mean ORP in shallow IWs generally decreased with increasing water-level depths and were inversely related to mean dissolved arsenic concentrations because of more reducing conditions (smaller ORP) at larger depths below the land surface. Larger dissolved arsenic concentrations in the shallow parts of the aquifer were associated with decreases in water levels and a subsequent decrease in ORP and thus more reducing conditions. Atrazine was detected in about 58 percent of shallow and 28 percent of deep IWs and did not exceed the primary Federal criterion of 3 μg/L in any groundwater samples. Atrazine concentrations in shallow IWs generally were largest in the northwest part of the study area near the North Branch Kisiwa Creek, and atrazine concentrations in deep IWs generally were largest most often in the southern part of the study area. Gross α radioactivity concentrations exceeded the primary Federal criterion of 15 picocuries per liter in 4 percent of shallow IW samples. Gross α and gross β radioactivity concentrations generally were larger in the southern third of the aquifer. Most groundwater-sample-simulated minerals saturation indices (SIs) were consistently negative (undersaturated). Minerals that had SI values that were consistently or typically positive (oversaturated) included iron oxide, hydroxide, and quartz-group minerals. Several SI values for arsenic- and manganese-bearing minerals were consistently negative. Some manganese-bearing mineral SI values ranged from undersaturated to oversaturated in shallow and deep IWs during the study. Several carbonate minerals in shallow and deep IWs varied across their equilibrium state. Calcite SI values were larger more often in the deep parts of the aquifer and did not show a clear distributional pattern. Mean and median calcite SI values for shallow and deep IWs were negative (undersaturated) indicating the potential for calcite dissolution if calcite is present for a substantial part of the study period. However, some individual calcite SI values in this study indicated saturation and subsequent calcite precipitation may occur in the study area, potentially resulting in formation of calcite mineral deposits that may reduce efficiency of injection wells. SI values with respect to iron hydroxide varied across their equilibrium states. Mean and median SI values with respect to iron hydroxide were undersaturated in shallow and deep IWs; however, some samples had positive SI values indicating there is potential for iron hydroxide precipitation, possibly caused by leaching and oxidation of iron-containing minerals, like pyrite, in the aquifer material.

Kansas↗

Weathering, landscape equilibrium, and carbon in four watersheds in eastern Puerto Rico: Chapter H in Water quality and landscape processes of four watersheds in eastern Puerto Rico

The U.S. Geological Survey's Water, Energy, and Biogeochemical Budgets (WEBB) program research in eastern Puerto Rico involves a double pair-wise comparison of four montane river basins, two on granitic bedrock and two on fine-grained volcaniclastic bedrock; for each rock type, one is forested and the other is developed. A confounding factor in this comparison is that the developed watersheds are substantially drier than the forested (runoff of 900–1,600 millimeters per year compared with 2,800–3,700 millimeters per year). To reduce the effects of contrasting runoff, the relation between annual runoff and annual constituent yield were used to estimate mean-annual yields at a common, intermediate mean-annual runoff of 1,860 millimeters per year. Upon projection to this intermediate runoff, the ranges of mean-annual yields among all watersheds became more compact or did not substantially change for dissolved bedrock, sodium, silica, chloride, dissolved organic carbon, and calcium. These constituents are the primary indicators of chemical weathering, biological activity on the landscape, and atmospheric inputs; the narrow ranges indicate little preferential influence by either geology or land cover. The projected yields of biologically active constituents (potassium, nitrate, ammonium ion, phosphate), and particulate constituents (suspended bedrock and particulate organic carbon) were considerably greater for developed landscapes compared with forested watersheds, consistent with the known effects of land clearing and human waste inputs. Equilibrium rates of combined chemical and physical weathering were estimated by using a method based on concentrations of silicon and sodium in bedrock, river-borne solids, and river-borne solutes. The observed rates of landscape denudation greatly exceed rates expected for a dynamic equilibrium, except possibly for the forested watershed on volcaniclastic rock. Deforestation and agriculture can explain the accelerated physical erosion in the two developed watersheds. Because there has been no appreciable deforestation, something else, possibly climate or forest-quality change, must explain the accelerated erosion in the forested watersheds on granitic rocks. Particulate organic carbon yields are closely linked to sediment yields. This relation implies that much of the particulate organic carbon transport in the four rivers is being caused by this enhanced erosion aided by landslides and fast carbon recovery. The increase in particulate organic carbon yields over equilibrium is estimated to range from 300 kilomoles per square kilometer per year (6 metric tons carbon per square kilometer per year) to 1,700 kilomoles per square kilometer per year (22 metric tons carbon per square kilometer per year) and is consistent with human-accelerated particulate-organic-carbon erosion and burial observed globally. There is no strong evidence of human perturbation of silicate weathering in the four study watersheds, and differences in dissolved inorganic carbon are consistent with watershed geology. Although dissolved organic carbon is slightly elevated in the developed watersheds, that elevation is not enough to unambiguously demonstrate human causes; more work is needed. Accordingly, the dissolved organic carbon and dissolved inorganic carbon yields of tropical rivers, although large, are of secondary importance in the study of the anthropgenically perturbed carbon cycle.

Professional Paper↗

Floods of Ohio and Mississippi Rivers, January-February 1937, with a section on the Flood deposits of the Ohio River, January-February 1937

In January and February 1937 the Ohio and mid-Mississippi Rivers experienced floods which, over reaches many hundreds of miles in length, exceeded all previously recorded stages. When measured by the loss of life and property, extent of damage, and general disruption of human activities, these floods constituted a major catastrophe. The floods were caused by a succession of heavy rainstorms that began late in December 1936 and continued nearly to the end of the following January. Although the storms covered a considerable part of the lower Mississippi River Basin and almost the entire Ohio River Basin, the center of heaviest precipitation was in the middle and lower portions of the Ohio River Valley. The total storm period can be subdivided into several individual storms, which were more or less clearly demarked by short intervening periods of little or no precipitation. Although the individual storm periods were the same or nearly the same over wide areas, their subdivisions were somewhat different in the most widely separated parts of the affected areas, with intermediate gradations in the intervening areas. The heaviest rainfall--that of January 20 to 25--was centered in the lower Ohio Valley, and, falling as it did upon a region with soil saturated and waterways already running full, it had the effect of producing extreme floods. The small quantity of snow on the ground over the higher eastern parts of the area at the beginning of the storm period disappeared in a short time. Some of the precipitation occurred in the form of snow, but this snow and the associated cold weather were much less significant in their influence on the floods than in the misery and discomfort they caused to ill-sheltered flood refugees and flood-bound people. Sequence and time of the storms were such that in the upper and smaller tributary basins the associated flood rises tended to clear to a notable degree before the-next flood rises came; hence many of these tributaries were at no time in extreme flood. In the lower reaches of the largest tributaries, and especially on the middle and lower reaches of the Ohio River, there were extreme and almost continuously increasing accumulations of run-off, which culminated in the region of Louisville, Ky., in stages 10 or 11 feet higher than any previously known. The precipitation was heaviest in the Ohio River Basin, and the flood in the Mississippi River, like other notable floods of the past, was caused largely by the extraordinary contributions from the Ohio River. The river stages exceeded those previously recorded for the lower 700 miles on the Ohio River and for 250 miles .on the Mississippi River below the Ohio. At Cairo, Ill., at the mouth Of the Ohio River, the river stage was higher for a period of 19 days, from January 24 to February 11, than at any previous time on record. The height above previous flood stages diminished materially as the flood progressed down the Mississippi. The mean precipitation ever the Ohio River Basin during the storm period was. 12.85 inches. The snow on the ground at the beginning of the period is estimated to have been equivalent to a mean depth of 0.10 inch of water over the basin. Out of the total precipitation 8.9 inches appeared as flood flow. On January 26 the computed volume of water in the stream channels of the Ohio River Basin was 56,000,000 acre-feet, equivalent to a depth of 5.1 inches over the drainage basin. The maximum discharge of the Ohio River at its mouth was 1,880,000 second-feet on February 1. On February 2, the day of the crest stage at the mouth of the Ohio, the computed volume of water on the surface channel system was equivalent to a depth of 3.7 inches over the drainage basin, of which 2.4 inches was in the 337-mile reach of the Ohio River between Louisville, Ky., and the mouth. This water-supply paper presents records of stage and discharge for the period including the floods at about 250 measurement stations, records of stage and discharge for the period including the floods at about 250 measurement stations, records of storage in many reservoirs, a summary of peak discharges with comparative data for other floods at about 470 measurement points, and tables showing crest stages along an aggregate length of stream channel for 5,000 miles. The report also includes basic information in regard to the weather associated with the floods, results of detailed studies of the rainfall and run-off, analyses of the volume of flood waters'in the surface channel systems during the progress of the floods, and many other kinds of flood information. Following the main flood report is a brief report entitled "Flood" deposits of the Ohio River, January-February 1937, a study of sedimentation." An abstract of that report is presented on page 693.

Ohio River, Mississippi River↗

Channel morphology of Cottonwood Creek near Cottonwood, California, from 1940 to 1985

Proposed construction of two dams on Cottonwood Creek California , has caused concern that resulting streamflow modification may alter downstream channel morphology. Baseline information on Cottonwood Creek channel characteristics from 1982-83 field surveys and 1940-84 aerial photographs indicates an alluvial channel that consists of a braided inner main channel within a broader flood channel, with no clear topographic break between the main and flood channels. The braided channel is subject to large and rapid position shifts in the flood channel, which meanders within the valley fill; however, the position of the flood channel has remained relatively stable since 1940. Mean slope of Cottonwood Creek is 0.0017 and that of South Fork Cottonwood Creek is 0.0020. Fluctuations of mean bed elevation appear to be random with no apparent long-term trend of aggradation or degradation. Mean bed-material size ranged from 3 to 82 millimeters. Low-flow channel sinuosity ranged from 1.04 to 1.47 from 1940 to 1984. Cumulative lateral migration for Cottonwood Creek decreased upstream, while for South Fork, it remained relatively constant. Net lateral migration was toward the right bank on Cottonwood Creek whereas no trend in net lateral migration is apparent for the South Fork. (USGS)

Water-Resources Investigations Report↗

Tectonostratigraphic terranes and their Paleozoic boundaries in the central and southern Appalachians

Parts of the central and southern Appalachian orogen appear to have evolved away from Proterozoic North America (Laurentia) and to have been accreted to it during the Paleozoic orogenies that collectively formed the orogen. Identifying each tectonostratigraphic terrane is a necessary step in understanding the evolution of the orogen. The terranes in the central and southern Appalachians are delineated, interpreted, and classified with varying degrees of confidence as: (1) Laurentian native terranes , (2) internal continental terranes of the Appalachian orogen, (3) disrupted terranes , (4) possible oceanic crustal remnants, (5) volcanic-arc terranes , (6) a continental terrane of Gondwanaland affinity, and (7) metamorphic complexes of undetermined affinity. The Laurentian native terranes consist of external massifs of Laurentian basement (Grenvillian and older), their rift- and shelf-facies cover rocks, and slope-rise prism deposits. External massifs are present in the Blue Ridge tectonic province, Reading Prong, and Honey Brook Upland. Rocks of the Talladega block are stratigraphically tied to Laurentia and, with the possible exception of the Hillabee greenstone, are also considered native. Offshore, deep-water, post-rift deposits of the Hamburg and Westminster terranes have no direct stratigraphic ties to Laurentia and are considered discrete native (not suspect) terranes . The internal continental terranes of the Appalachian orogen are isolated massifs of Middle Proterozoic (Grenvillian) continental basement and their cover sequences that occur within the metamorphic core of the orogen. These terranes , the Baltimore, Sauratown, and Pine Mountain terranes , could be either structurally isolated outliers of Laurentia or microcontinental fragments of Laurentian crust displaced by rifting or transcurrent faulting and later reassembled. Disrupted terranes in the central and southern Appalachians contain mélange complexes as well as more coherent terrane fragments (volcanic, ophiolitic, or continental) intermingled with the mélange complexes. Those identified include the Jefferson, Potomac, Smith River, Inner Piedmont, Falls Lake, Juliette, and Sussex terranes . The Bel Air-Rising Sun terrane (Baltimore Complex) in Maryland and Pennsylvania is the only terrane named separately as a possible oceanic crustal remnant. Similar mafic and ultramafic complexes are present in all of the disrupted terranes , but are too small to consider as separate terranes . Volcanic-arc terranes include the Chopawamsic, Carolina, Spring Hope, Roanoke Rapids, and Charleston terranes . The only terrane recognized as a continental terrane of Gondwanaland affinity is the Suwannee terrane , which contains rocks believed to correlate with those now exposed in west Africa. Metamorphic complexes of undetermined affinity are terranes that could not be clearly classified on the basis of available data. These include the Milton, Gaffney, Uchee, Crabtree, Goochland, Wilmington, and Hatteras terranes . The Penobscottian, Taconian, Acadian, and Alleghanian Paleozoic compressional events collectively assembled the various terranes into what is now the Appalachian orogen. Only the central and southern parts of the U.S. Appalachians are considered here. The Penobscottian orogeny, about 550 to 490 Ma, amalgamated the Potomac, the Chopawamsic, probably the Bel Air-Rising Sun, and possibly other exotic terranes at some unknown distance from Laurentia. This was followed by the Taconian orogeny, about 470 to 440 Ma, which accreted the previously amalgamated terranes and probably other terranes such as the Carolina terrane to Laurentia. The younger age limit for the Taconian event is partly constrained by Middle and Late Ordovician faunal assemblages in successor basin deposits of the Arvonia Slate and Quantico Formation. The significance of the Acadian orogeny, dated about 400 to 380 Ma in New England, is unclear in the central and southern Appalachians . In the Talladega block of Alabama and Georgia, an Early to Middle Devonian dynamothermal event is firmly bracketed between Early Devonian fossils and K-Ar ages that indicate a thermal peak no later than Middle Devonian time. A regional tectonothermal event and faulting of approximately this age are also suggested by isotopic studies in terranes to the east. The late Paleozoic (Alleghanian) continental collision between Laurentia and Gondwanaland, which formed the supercontinent Pangea, marks the final stage of accretionary history in the Appalachian -Caledonide orogen. Effects evident in the central and southern Appalachian region include: (1) the accretion of the Suwannee terrane and perhaps the Charleston terrane to what is now North America, (2) slicing and shifting of terranes along dextral strike-slip faults, particularly in the eastern Piedmont, (3) westward transport of native and previously accreted terranes in the western Piedmont and Blue Ridge as part of a composite crystalline thrust sheet, (4) deposition of clastic wedges in the Appalachian foreland, and (5) imbricate thrusting and folding of the resultant strata in the Valley and Ridge Province.

Appalachians↗

Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths

This study investigates connectivity in a small portion of the extremely heterogeneous aquifer at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi. A total of 19 fully penetrating soil cores were collected from a rectangular grid of 4 m by 4 m. Detailed grain size analysis was performed on 5 cm segments of each core, yielding 1740 hydraulic conductivity ( K ) estimates. Three different geostatistical simulation methods were used to generate 3‐D conditional realizations of the K field for the sampled block. Particle tracking calculations showed that the fastest particles, as represented by the first 5% to arrive, converge along preferential flow paths and exit the model domain within preferred areas. These 5% fastest flow paths accounted for about 40% of the flow. The distribution of preferential flow paths and particle exit locations is clearly influenced by the occurrence of clusters formed by interconnected cells with K equal to or greater than the 0.9 decile of the data distribution (10% of the volume). The fraction of particle paths within the high‐ K clusters ranges from 43% to 69%. In variogram‐based K fields, some of the fastest paths are through media with lower K values, suggesting that transport connectivity may not require fully connected zones of relatively homogenous K . The high degree of flow and transport connectivity was confirmed by the values of two groups of connectivity indicators. In particular, the ratio between effective and geometric mean K (on average, about 2) and the ratio between the average arrival time and the arrival time of the fastest particles (on average, about 9) are consistent with flow and advective transport behavior characterized by channeling along preferential flow paths.

Water Resources Research↗

Surface faults in the gulf coastal plain between Victoria and Beaumont, Texas

Displacement of the land surface by faulting is widespread in the Houston-Galveston region, an area which has undergone moderate to severe land subsidence associated with fluid withdrawal (principally water, and to a lesser extent, oil and gas). A causative link between subsidence and fluid extraction has been convincingly reported in the published literature. However, the degree to which fluid withdrawal affects fault movement in the Texas Gulf Coast, and the mechanism(s) by which this occurs are as yet unclear. Faults that offset the ground surface are not confined to the large (>6000-km 2 ) subsidence “bowl” centered on Houston, but rather are common and characteristic features of Gulf Coast geology. Current observations and conclusions concerning surface faults mapped in a 35,000-km 2 area between Victoria and Beaumont, Texas (which area includes the Houston subsidence bowl) may be summarized as follows: (1) Hundreds of faults cutting the Pleistocene and Holocene sediments exposed in the coastal plain have been mapped. Many faults lie well outside the Houston-Galveston region; of these, more than 10% are active, as shown by such features as displaced, fractured, and patched road surfaces, structural failure of buildings astride faults, and deformed railroad tracks. (2) Complex patterns of surface faults are common above salt domes. Both radial patterns (for example, in High Island, Blue Ridge, Clam Lake, and Clinton domes) and crestal grabens (for example, in the South Houston and Friendswood-Webster domes) have been recognized. Elongate grabens connecting several known and suspected salt domes, such as the fault zone connecting Mykawa, Friendswood-Webster, and Clear Lake domes, suggest fault development above rising salt ridges. (3) Surface faults associated with salt domes tend to be short (<5 km in length), numerous, curved in map view, and of diverse trend. Intersecting faults are common. In contrast, surface faults in areas unaffected by salt diapirism are frequently mappable for appreciable distances (>10 km), occur singly or in simple grabens, have gently sinuous traces, and tend to lie roughly parallel to the ENE-NE “coastwise” trend common to regional growth faults identified in subsurface Tertiary sediments. (4) Evidence to support the thesis that surface scarps are the shallow expression of faults extending downward into the Tertiary section is mostly indirect, but nonetheless reasonably convincing. Certainly the patterns of crestal grabens and radiating faults mapped on the surface above salt domes are more than happenstance; analogous fault patterns have been documented around these structures at depth. Similarly, some of the long surface faults not associated with salt domes seem to have subsurface counterparts among known regional growth faults documented through well logs and seismic data. Correlations between surface scarps and faults offsetting subsurface data are not conclusive because of the large vertical distances (1900- 3800 m) involved in making the most of the inferred connections. Nevertheless, the large number of successful correlations - in trend, movement sense, and position - suggests that many surface scarps represent merely the most recent displacements on faults formed during the Tertiary. (5) Upstream-facing fault scarps in this region of low relief can be significant impediments to streams. Locally, both abandoned, mud-filled Pleistocene distributary channels and, more commonly, Holocene drainage lines still occupied by perennial streams reflect the influence of faulting on their development. Some bend sharply near faults and have tended to flow along or pond against the base of scarps; others meander within topographically expressed grabens. Such evidence for Quaternary displacement of the ground surface is widespread in the Texas Gulf coast. In the general, however, streams in areas now offset by faulting show no disruption of their courses where they cross fault scarps. Such scarps are probably very young, and where they can be demonstrated to partly or wholly predate fluid withdrawal, very recent natural fault activity is indicated. (6) Early aerial photographs (1930) of the entire region and topographic maps (1915-16 surveys) of Harris County (Houston and vicinity) show that many faults had already displaced the land surface at a time when appreciable pressure declines in subjacent strata were localized to relatively few areas of large-scale pumping. Prehistoric faulting of the land surface, as noted above, appears to have affected much of the Texas Gulf Coast. (7) A relation between groundwater extraction and current motion on active faults is suspected because of the increased incidence of ground failure in the Houston-Galveston subsidence bowl. This argument is weakened somewhat by recognition of numerous surface faults, some of them active today, far beyond the periphery of the strongly subsiding area. Moreover, tilt beam records from two monitored faults in northwest Houston and accounts of fault damage from local residents demonstrate a complex, episodic nature of fault creep which can only partially be correlated with groundwater production. Nevertheless, although specific mechanisms are in doubt, the extraction of groundwater from shallow (<800-m) sands is probably a major factor in contributing to current displacement of the ground surface in the Houston-Galveston region. Within this large area, the number of faults recognizable from aerial photographs has increased at least tenfold between 1930 and 1970. Elsewhere in the Texas Gulf Coast only a moderate increase has been noted, some of which is possibly attributable to oil and gas production. Surface fault density in the Houston-Galveston region is far greater than in any other area of the Texas Gulf Coast investigated to date. A plausible explanation for these differences is that large overdrafts of groundwater over an extended period of time in the Houston-Galveston region have stimulated fault activity there. Throughout the Texas Gulf Coast, however, a natural contribution to fault motion remains a distinct possibility.

Texas↗

Estimating flood magnitude and frequency on streams and rivers in Connecticut, based on data through water year 2015

The U.S. Geological Survey, in cooperation with the Connecticut Department of Transportation, updated flood-frequency estimates with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities (2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively) for 141 streamgages in Connecticut and 11 streamgages in adjacent States using annual peak-flow data through water year 2015. Peak-flow regression equations were derived for estimating flows at ungaged stream sites with annual exceedance probabilities from 50 to 0.2 percent. Methods for estimating prediction intervals for the peak-flow regression equations are presented. The regression equations are applicable for basins in Connecticut with drainage areas ranging from 0.69 to 325 square miles that are not affected by flood-control regulation or flow diversions. The flood discharges for select annual exceedance probabilities were estimated following new (2018) national guidelines for flood-frequency analyses. New guidelines have improved statistical methods for flood-frequency analysis including (1) the expected moments algorithm to help describe uncertainty in annual peak flows and to better represent missing and historical record and (2) the generalized multiple Grubbs-Beck test to screen out potentially influential low outliers and to better fit the upper end of the peak-flow distribution. Additionally, a new regional skew (0.37) derived for New England was used in the flood-frequency analysis for the streamgages. Annual peak flows were analyzed for trends for four time periods (30, 50, 70, and 90 years) through 2015. Trend results show some statistical evidence of increasing peak flows in each of the time periods analyzed; however, multidecadal climate cycles may be influencing the number and magnitude of the trends. Historical peak-flow trends in and near Connecticut do not offer clear and convincing evidence for incorporating trends into flood-frequency analyses. For this study, the traditional assumption of stationarity is used with no adjustment for trends. Generalized least squares regression techniques were used to develop the final set of multivariable regression equations for estimating flood discharges with 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. The standard error of prediction for the regional regression equations ranged from 26.3 to 45.0 percent. The standard error of prediction was slightly smaller in the current study compared to the 2004 study, indicating an improvement in the predictive ability of the equations (6 percent smaller at the 50-percent annual exceedance probability to about 1 percent smaller at the 1-percent annual exceedance probability). Generalized least squares regression techniques also were used to develop a one-variable (drainage-area-only) equation. Drainage-area-only equations can be used as an alternative to the multiexplanatory variable statewide regression equations if decreased accuracy is acceptable. The revised statistical procedures and additional streamgage data applied in the current study result in a more accurate representation of peak-flow conditions in Connecticut than was previously available. The regional regression equations will be integrated in the U.S. Geological Survey StreamStats program, which estimates basin and climatic characteristics and streamflow statistics at user-selected ungaged stream sites.

Connecticut↗

Assessment of undiscovered conventionally recoverable petroleum resources of the Northwest European region

The estimates of undiscovered conventionally recoverable petroleum resources in the northwest European region at probability levels of 95 percent, 5 percent, statistical mean, and mode are for oil (in billions of barrels): 9, 34, 20, and 15; and for gas (in trillions of cubic feet): 92, 258, 167, and 162. The occurrence of petroleum can be accounted for in two distinct geological plays located in the various subbasins of the region. Play I is associated with the distribution of mature source rocks of Late Jurassic age relative to four distinct trapping conditions. The play has been demonstrated productive mostly in the Viking and Central Grabens of the North Sea, where the shale has been buried to optimum depths for the generation of both oil and gas. To the north of 62 ? N. latitude up to the Barents Sea, source rocks become increasingly deeply buried and are interpreted to be dominantly gas prone; a narrow band of potentially oil-prone shales tracks most of the coast of Norway, but water depths in favorable localities commonly range from 600 to 1,200 feet. To the south of the Central Graben, the Jurassic source rocks are either immature or minimally productive because of a change in facies. Undrilled traps remain within the favorable source-rock area, and exploration will continue to challenge the boundaries of conventional wisdom, especially on the Norwegian side where little has .been reported on the geology of the adjoining Bergen High or Horda Basin, though, reportedly, the Jurassic source rocks are missing on the high and are immature in the southern part of the basin. Play II is associated with the distribution of a coal facies of Carboniferous age that is mature for the generation of gas and locally underlies favorable reservoir and sealing rocks. The play is limited largely by facies development to the present area of discovery and production but is limited as well to the southeast into onshore Netherlands and Germany by the unfavorable economics of an increasing nitrogen content in the gas. This increase is apparently caused by excessive temperatures associated with increasing depth of burial of the source rock. The history of discovery in the North Sea would appear to deny the commonly held maxim that large fields are found first and early in the exploration process. However, if the discovery data are examined from the perspective of the award date of each exploration license, then it is clear that the largest fields and most of the reserves have indeed been found early in the exploration process of a particular license. Discoveries made within 1 year of granting the license are on average large giants, and they account for slightly less than two-thirds of the original reserves. Discoveries made within 2 to 5 years of the granting of the license are on average less than giant size and smaller than increment-l-year discoveries by a factor of 4; these fields account for a little less than one-third of the reserves. Those fields found 6 or more years after the granting of the license are relatively small and account for 20 percent of all discoveries but only 4 percent of total original reserves. These data suggest that a measure of an area's exploration maturity is the length of time elapsed since the award of the concession.

Circular↗