Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Washington Geology”

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

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

1,599 records · Page 89Linked to original sources

Stratigraphy and structure of the tatum salt dome area, southeastern Mississippi and northeastern Washington Parish, Louisiana

In the 3000-square-mile area of southeastern Mississippi and northeastern Washington Parish , Louisiana , which has Tatum dome in its center, rocks of known Jurassic to Recent age are more than 20,000 feet thick. They are underlain by an unknown thickness of Louann Salt of Jurassic (?) age. The age, thickness, and nature of the sedimentary rocks between the salt and the basement, as well as the character of the basement, are unknown. The salt in the northern half of the area has moved upward to shallow depths as diapirs in overlying rocks. Outside of the area of shallow-piercement domes there are more deeply seated domes, ridges, and anticlines that may be the result of vertical movement of the salt . The Wiggins anticline in the south and, southwest of it, an area dipping more steeply to the southwest probably were not affected by salt movement. Structural maps of four horizons and isopach maps of rocks in the three intervals between them give information about the sedimentary and structural history of the region since the end of Early Cretaceous time. The horizons contoured mark either unconformities or the upper boundaries of zones of fairly uniform structural history. The top of the Lower Cretaceous rocks dips generally to the west. Upper Cretaceous rocks dip to the west-southwest and thicken to the northeast. Pre-Miocene Tertiary rocks dip to the southwest and thicken to the northwest. Miocene and later rocks dip and thicken to the south-southwest. The structure of the older formations is considerably more irregular than that of the younger ones; domes, ridges, and basins that were conspicuous in the Early Cretaceous became less pronounced as younger sediments progressively covered them. Of the eight shallow-piercement domes with which the author is familiar, three are known to have penetrated Oligocene rocks; two, middle Eocene; one, lower Eocene; one, Upper Cretaceous; and one, Lower Cretaceous. Some movement has taken place at Tatum dome since the sandy, gravelly Citronelle Formation was laid down.

Louisiana, Mississippi↗

Yakima basalt of the Tieton River area, south-central Washington

Up to 1700 feet of the upper Miocene-lower Pliocene Yakima Basalt of the Columbia River Group underlie much of the eastern flank of the Cascade Range in the Tieton River area, Yakima County, Washington. Local prebasalt relief was more than 1700 feet, so thicknesses of each of the 15 exposed flows vary widely. Single flows can be traced for many miles, and terminate only against local topographic highs. The flows show typical colonnade-entablature jointing, and commonly overlie thin pillow-palagonite complexes. The basalts have Yakima-type chemistry. Plagioclase varies chiefly from An 65 (microphenocrysts) to An 45 (microlites). Clinopyroxenes range from Ca-rich pigeonite to augite, with subcalcic augite most abundant. Complex continuous zoning with respect to optic angles occurs between all the clinopyroxene phases. Olivine is sparse, and phenocrysts of all minerals are rare. Plots of 40 modal analyses indicate that plagioclase and pyroxene began crystallizing at about the same time, and crystallized at the same rate until the flows were at least 75 per cent crystalline. Individual flows differ in their plagioclase/pyroxene ratios and, to a lesser degree, in other microscopic characteristics. Therefore they can be correlated between the measured sections. The stratigraphy thus defined indicates that the youngest flow in the area is slightly older than the Vantage Sandstone Member (about 13.5 m.y.), a prominent sedimentary interbed farther east on the Columbia River Plateau. Flow directions show that the basalts advanced into the area from the east and southeast. Sedimentary interbeds between some of the flows contain a heavy mineral suite indicative of a northern or northeastern provenance, and crossbedding measurements suggest westward and southwestward current directions. The regional paleoslope, therefore, sloped westward or southwestward; and the basalts probably extended somewhat beyond the present Cascade Crest before being dammed by the ancient Western Cascades. Floods of pyrogenic hornblende in the interbeds and overlapping K/A ages suggest sporadic explosive activity of Tatoosh-type plutons in the Cascades contemporaneous with Yakima extrusions. The basalts are warped into five nearly west-trending folds and an eastward-sloping homocline. The homocline is related directly to Cascade uplift, which may have begun at about the time that Yakima-type flows ceased flooding the area.

Geological Society of America Bulletin↗

Coarse sediment dynamics in a large glaciated river system: Holocene history and storage dynamics dictate contemporary climate sensitivity

The gravel-bedded White River drains a 1279 km 2 basin in Washington State, with lowlands sculpted by continental glaciation and headwaters on an actively glaciated stratovolcano. Chronic aggradation along an alluvial fan near the river’s mouth has progressively reduced flood conveyance. In order to better understand how forecasted climate change may influence coarse sediment delivery and aggradation rates in this lowland depositional setting, we assessed the contemporary delivery and routing of coarse sediment through the watershed; this assessment was based on a rich set of topographic, sedimentologic, and hydrologic data from the past century, with a focus on repeat high-resolution topographic surveys from the past decade. We found that most of the lower river’s contemporary bed-load flux originates from persistent erosion of alluvial deposits in the lower watershed. This erosion is a response to a drop in local base level caused by a major avulsion across the fan in 1906 and then augmented by subsequent dredging. The 1906 avulsion and modern disequilibrium valley profiles reflect landscape conditioning by continental glaciation and a massive mid-Holocene lahar. In the proglacial headwaters, infrequent large sediment pulses have accomplished most of the observed coarse sediment export, with exported material blanketing downstream valley floors; during typical floods, transported bed material is largely sourced from erosion of these valley floor deposits. Throughout the watershed, we observe decadal-scale coarse sediment dynamics strongly related to the filling or emptying of valley-scale sediment storage over 10 2 −10 4 yr time scales, often in response to major disturbances that either emplace large deposits or influence their redistribution. Paraglacial responses in large watersheds are suggested to be inherently complicated and punctuated as a result of internal landform interactions and stochastic/threshold-dependent events. We argue, in combination, that Holocene disturbance, storage dynamics, and human flow modification make coarse sediment fluxes in the lower White River relatively insensitive to decadal climate variability. Results highlight the degree to which river sensitivity to contemporary disturbance, climatic or otherwise, may be contingent on local and idiosyncratic watershed histories, underscoring the need to unpack those histories while demonstrating the utility of watershed-scale high-resolution topography toward that end.

Washington↗

Late Mesozoic and possible early Tertiary accretion in western Washington State: The Helena-Haystack mélange and the Darrington-Devils Mountain fault zone

The Helena-Haystack mélange (HH mélange) and coincident Darrington-Devils Mountain fault zone (DDMFZ) in northwestern Washington separate two terranes, the Northwest Cascade System (NWCS) and the western and eastern mélange belts (WEMB). The two terranes of Paleozoic and Mesozoic rocks superficially resemble each other but record considerable differences in structural and metamorphic history. The HH mélange is a serpentinite-matrix mélange containing blocks of adjacent terranes but also exotic blocks of schistose metavolcanic rocks and Jurassic tonalite and associated amphibolite. The HH mélange must have formed between early Cretaceous and late middle Eocene time, because it contains tectonic clasts of early Cretaceous Shuksan Greenschist and is overlain by late middle Eocene sedimentary and volcanic rocks. Less certain constraints on its age are a tectonic clast of metarhyolite that yields 90 Ma metamorphic ages and the presumption that the mélange was emplaced before the outboard Olympic terrane arrived at about 50 Ma. The apparent continuity of the HH mélange and the Decatur terrane of the San Juan Islands suggests that the mélange is the strongly tectonized equivalent of the Fidalgo ophiolite. The out-crop pattern suggests that the HH mélange overlies rocks of the NWCS and it may have formed when the WEMB terranes were thrust over rocks of the NWCS. Much of the exposed belt of the HH mélange is overlain by late middle Eocene feldspathic sandstone and volcanic rocks of the Barlow Pass Volcanics of Vance (1957a), which are cut by numerous faults of the DDMFZ paralleling the mélange. The Barlow Pass Volcanics appear to overlie the Straight Creek fault without large offset, but a displaced exotic block of amphibolite with attached early or early middle Eocene(?) sandstone in the mélange suggests that strike-slip movement along the DDMFZ was synchronous with movement on the Straight Creek fault, and stretched cobbles in the conglomerates of the Barlow Pass Volcanics suggest post-Straight Creek movement. The possible continuation of the DDMFZ to the northwest as the San Juan and the West Coast faults on Vancouver Island suggests That the structure has had a major role in the emplacement of all the westernmost terranes in the Pacific Northwest. This major suture is strongly bowed to the northeast opposite the great oroclinal bend of the Olympic terrane, suggesting that the emplacement of that terrane may have deformed a once straighter strike-slip zone.

Washington↗

Rb-Sr, K-Ar, and stable isotope evidence for the ages and sources of fluid components of gold-bearing quartz veins in the northern Sierra Nevada foothills metamorphic belt, California

Gold-bearing quartz veins occur in and near major fault zones in deformed oceanic and island-arc rocks west of the main outcrop of the Sierra Nevada composite batholith. Veins typically occupy minor reverse faults that crosscut blueschist to amphibolite-grade metamorphic rocks whose metamorphic ages range from early Paleozoic to Jurassic. Vein micas and carbonate-quartz-mica assemblages that formed by hydrothermal metasomatism of ultramafic wall rocks in the Alleghany, Grass Valley, Washington, and Mother Lode districts yield concordant K-Ar and Rb-Sr ages. The dated veins are significantly younger than prograde metamorphism, penetrative deformation, and accretion of their host rocks to the continental margin. New and previously published mineralization ages from 13 localities in the Sierra foothills range from about 140 to 110 m.y. ago, with mean and median between 120 and 115 m.y. The age relations suggest that mineralizing fluids were set in motion by deep magmatic activity related to the resumption of east-dipping subduction along the western margin of North America following the Late Jurassic Nevadan collision event.CO 2 -bearing fluids responsible for metasomatism and much of the vein mica, carbonate, albite, and quartz deposition in several northern mines were isotopically heavy (delta 18 O [asymp] 8-14ppm; delta D between about -10 and -50ppm) and do not resemble seawater, magmatic, or meteoric waters. Metasomatic and vein-filling mica, dolomite, magnesite, and quartz in altered ultramafic rocks generally formed from fluids with similar Sr and O isotope ratios at a given locality. Consistent quartz-mica delta 18 O fractionations (delta 18 O (sub Q-M) = 4.5-4.9ppm) from various localities imply uniform equilibration temperatures, probably between 300 degrees and 350 degrees C. On a local (mine) scale, fluids responsible for both carbonate alteration of mafic and ultramafic wall rocks and albitic alteration of felsic and pelitic rocks had similar Sr isotope ratios.Samples from three veins in the central Alleghany district fit a 115.7 + or - 3-m.y. Rb-Sr isochron with a ( 87 Sr/ 86 Sr) i value of approximately 0.7119. Inferred 87 Sr/ 86 Sr ratios of metasomatic fluids from mines in different parts of the foothills region vary considerably (0.704-0.718), suggesting that Sr was derived from sources ranging from "western assemblage" Mesozoic ophiolitic or arc volcanic rocks to early Paleozoic continent-derived clastic rocks of the Shoo Fly Complex. Systematic geographic variations in both Sr and O isotopes can be rationalized by assuming extensive fluid interaction with rocks similar to the ones that are exposed within a few kilometers of the veins, but the ultimate sources of the fluids, and of Au and other constituents, may be independent of these. Isotopically lighter (meteoric?) fluids deposited some late quartz overgrowths and occupied secondary fluid inclusions in earlier vein quartz.

Economic Geology↗

Estimated average annualized tsunami losses for the United States

Tsunami hazards are substantial threats to coastal communities across the United States (U.S.) and its territories. U.S. states and territories collaborate through the National Tsunami Hazard Mitigation Program (NTHMP) to develop their own tsunami-hazard information for outreach and evacuation planning. An effort to curate this tsunami-hazard information to support comprehensive risk analysis at the national level has not yet been completed. In support of this effort, the Federal Emergency Management Agency (FEMA) collaborated with the NTHMP, the National Oceanic and Atmospheric Administration (NOAA) and the U.S. Geological Survey (USGS) starting in 2023. This collaboration included the collection and analysis of existing tsunami hazard data and methods in the U.S. Tsunami subject matter experts identified and selected scientifically defensible methods for estimating the risks to buildings and populations in coastal communities. These efforts may support decision making regarding resilience policies, priorities, strategies and funding levels. Tsunamis can be triggered by earthquakes, subaerial or submarine landslides, volcanic eruptions, glacial calving, near-earth objects, weather or other events. These events can cause severe destruction, injuries, and loss of life due to powerful currents and flooding. Tsunamis pose a substantial threat to the western United States and all U.S. territories, as described below. ■ Hawaii is threatened by distant tsunamis due to its central location in the Pacific Ocean basin and has a history of local events. ■ Alaska, particularly the Aleutian Islands, faces local tsunami threats due to proximity to the Alaska-Aleutian Subduction Zone, as well as distant tsunamis from around the Pacific Ocean basin. ■ The western coast of the U.S. is threatened by distant tsunamis from around the Pacific Ocean basin and local source tsunamis from earthquakes generated within the Cascadia Subduction Zone in the Pacific Northwest. ■ American Samoa faces local tsunami threats from earthquakes generated in the nearby Tonga Trench, as well as distant tsunami threats. ■ Guam and the Commonwealth of the Northern Mariana Islands are threatened by local tsunamis from the nearby Mariana Subduction Zone, as well as distant sources from around the Pacific Ocean Basin. ■ Puerto Rico and the United States Virgin Islands are threatened by multiple local and distant tsunami sources, such as the Puerto Rico Trench (PRT), given their location in the complex seismic region of the Caribbean Sea. Several historical events stand out because of their catastrophic impacts. ■ In the Pacific Northwest, the 1700 Cascadia earthquake caused a tsunami that affected coastal Native American communities, though the extent of the damage is not fully documented (Ludwin, et al., 2005). ■ In Puerto Rico, the 1918 earthquake triggered a tsunami that caused $77 million in damage in 2022 dollars and 116 fatalities, primarily along the western coast (Coffman et al., 1982). ■ The 1946 Aleutian Islands earthquake triggered a massive tsunami that devastated Hilo, Hawaii, killing 158 people and resulting in approximately $375 million in damage (adjusted to 2022 dollars) (Fisher et al., 2023). ■ The 1964 Alaska earthquake (M 9.2) generated tsunamis that caused severe destruction in some communities across Alaska, Oregon, and California. This disaster led to a total of 124 fatalities and approximately $2.9 billion in property damage (adjusted to 2022 dollars) (Brocher et al., 2014) (Alaska Science Center, 2024). ■ In American Samoa, a tsunami generated by the 2009 Samoa earthquake (Mw 8.1) caused widespread devastation, resulting in 34 confirmed fatalities (Apatu et al., 2013) and economic losses exceeding $160 million (adjusted to 2022 dollars) (DHS, 2011). More recent events, including the 2010 Chile earthquake, the 2011 Japan earthquake, and the 2022 Tonga volcanic eruption, resulted in millions of dollars in damage to numerous ports and harbors in the U.S. South Pacific territories, Hawaii, and along the west coast of the U.S. (Lynett, et al., 2022) (Wilson, et al., 2013). Since these events, the expansion of the built environment in lowlying areas along the coast has increased the exposure of buildings and people, thereby further escalating community risk from tsunamis. This report provides a comprehensive national assessment of earthquake-generated tsunami risk. It does not include impacts from tsunamis generated by landslides, volcanic eruptions, glacial calving, near-earth objects, weather, or other events. This study is based on the best available hazard data from the U.S. Pacific Coast (California, Oregon and Washington), Alaska, Hawaii, U.S. Pacific Territories (American Samoa, Guam and Commonwealth of the Northern Mariana Islands) and Caribbean Territories (Puerto Rico and United States Virgin Islands). Tsunami risks associated with states along the East Coast, Gulf Coast, and Great Lakes are not included in this study because Hazus 6.1 software (FEMA 2024a) does not currently include the ability to analyze tsunami risk in those states. Once modeling capabilities and tsunami hazard data become available for additional states, FEMA may incorporate these data into future editions of this study.

Alaska, California, Hawaii Oregon, Washington↗

Cascadia Margin cold seeps: Subduction zone fluids, gas hydrates, and chemosynthetic habitats

Priority Geographic Area: The outer continental shelf and upper continental slope from Canada/U.S. border offshore Washington State to the Mendocino Fracture Zone (Northern California), entirely within the U.S. Exclusive Economic Zone (EEZ), from the outermost shelf to at least 2000 m water depth (Figure 1). Description of Priority Area: Since 2015, over a thousand water column gas plumes originating at seafloor gas seeps have been discovered landward of the Cascadia deformation front (e.g., Embley et al., 2016; Johnson et al., 2015, 2019; Merle and Embley, 2016; NA-95 Cruise Report, 2018; Riedel et al., 2018), adding to those that had long been known on Hydrate Ridge (e.g., Heeschen et al., 2003; Tréhu et al., 2004). The recently-discovered seeps stretch from offshore Vancouver Island to the Mendocino Fracture Zone and from the outer shelf to ~2000 m water depth, occurring both landward and seaward of the nominal limit for gas hydrate stability zone on the upper continental slope (Figure 1). Hundreds of seeps likely remain undiscovered. Water column imaging is incomplete both within the target geographic area and farther seaward, between the 2000 m isobath and the deformation front, which is the subject of an imaging study described in a white paper by Watt et al. The recently-discovered Cascadia Margin cold seeps partially overlap an important active margin gas hydrate province (Spence et al., 2001; Tréhu et al., 2003, 2004), as well as an area where sediments on the North American plate are folded and faulted and affected by fluids generated in the subduction complex beneath the Cascadia forearc (e.g., Saffer and Tobin, 2011). Several Ocean Drilling Program expeditions have focused on hydrate systems offshore Vancouver and Oregon (e.g., Riedel et al., 2009; Tréhu et al., 2004) and on the connection between the shallow and deep hydrogeologic systems. Cabled observatories now continuously monitor physical, chemical, and venting processes on south Hydrate Ridge (OOI; e.g., Philip et al., 2016a) and offshore Vancouver Island (NEPTUNE; e.g. Römer et al., 2016). Outside of these well-studied gas hydrate areas, a subset of the recently-discovered Cascadia seeps, including some that we visited with R/V Falkor in 2019 (e.g., https://schmidtocean.org/cruise/methane-seeps-at-edge-of-hydrate-stability/), also likely emit methane associated with shallow subseafloor gas hydrate systems. Other seeps are delivering not only methane, but also deep-derived gases (Baumberger et al., 2018, 2020) to the seafloor. Many Cascadia Margin seeps have also been recognized at water depths too shallow (e.g., 175 m) to be connected to gas hydrate dynamics. These seeps are postulated to be emitting gas and fluids that originated deep in accretionary wedge before migrating up normal faults generated during forearc extension associated with large earthquakes (Johnson et al., 2019). Only a small fraction of the recently discovered U.S. Cascadia Margin water column gas plumes has so far been verified by ROVs (Hercules from E/V Nautilus in 2016 and 2018; SuBastian from R/V Falkor in 2018 and 2019) to correspond to seafloor seeps. Careful scientific mapping, investigation, and sampling at the seeps have also been limited (e.g., Baumberger et al., 2018, 2020; Merle and Embley, 2016; Seabrook et al., 2018; Greinert et al. 2019). This white paper focuses on expanding exploration of already-identified U.S. Cascadia Margin cold seeps through a multipronged and multidisciplinary discovery program that could be accomplished with a variety of NOAA assets. The goals of the proposed exploration activities are to develop high-resolution maps of seep fields from deep ocean vehicles; to verify (and sample) seafloor gas emissions at the locations of water column plumes for compositional and isotopic studies; to map, sample, and conduct analyses on chemosynthetic communities and deep-sea coral habitats near seep sites to document species distributions and habitats as a function of depth and latitude along the margin; to collect seep geologic samples that can constrain the timing of methane emissions through geochronology; and to record environmental data (e.g., CTD) near the seafloor and in the water column above the seeps. Seafloor mapping using shipboard systems (multibeam/backscatter) would be needed to characterize seafloor features near seep sites. Water column imaging (EK60/80 and/or multibeam WCD data) conducted before and after seafloor explorations would capture active methane plumes and constrain temporal variations in seep emissions (e.g., Kannberg et al., 2013; Philip et al., 2016a, 2016b), which are known to vary on time scales as rapid as tidal cycles on this margin (e.g., Römer et al., 2016). What are the characterization and data needs in this area? Check all that apply: __x_ Biology, Geology, Physical Oceanography, Chemistry ___ Marine Archaeology ___ Other Provide a list or brief description of the data needed within this area, from your perspective: 1. Water column backscatter to image active gas plumes 2. High-resolution multibeam bathymetry, seafloor backscatter, and shallow sub-bottom imaging 3. Visual characterization and ground truthing of potential seeps, including high-resolution mapping and photography from near-seafloor vehicles; collection of seep-associated species, corals, sediments, authigenic carbonates, gases, and seawater Describe relevance to national security, conservation, and/or the economy: The Cascadia margin seeps provide significant ecosystem services, including habitat for commercially important fishes and support for diversity along the continental margin. Methane seeps are also biological hotspots for krill, plankton, and crustaceans, which in turn sustain higher trophic levels (e.g., whales). Methane-derived authigenic carbonates serve as a hard substrate for deep-sea corals and sponges on millennial time scales. The studies proposed here will elucidate the relationship among seep environments, deep-sea corals, sponges, fisheries, and other organisms and provide new insight into subduction zone and hydrate-associated fluids in this important seismogenic zone. The studies address fishery management concerns and inform future conservation of sensitive species (e.g., deep-sea corals) and benthic habitats. From your perspective, what makes this area unique? The Cascadia Margin seeps are a critical component of the leaky margin that stretches from Baja California to the Aleutian Arc along the Pacific coastline of North America. Cold seeps have been intensely studied on the Gulf of Mexico and U.S. Atlantic passive margins with a focus on chemosynthetic communities, deep-sea corals, and leakage of microbially-generated and/or thermogenic hydrocarbons; however, the recently-discovered Cascadia Margin seeps, as well as active margin seep systems in general, remain more poorly characterized. Such seeps not only contribute to the ocean carbon cycle (e.g., Pohlman et al., 2011), thereby fueling the base of the food chain in these settings, but also emit subduction zone fluids that provide clues about processes within the seismogenic zone and the accretionary complex. The Cascadia seeps area allows both biological (e.g., benthic habitats, coral distributions) and physical processes (e.g., generation of subduction zone fluids) to be studied along both depth (perpendicular to the deformation front) and latitudinal gradients.

California, Oregon, Washington↗

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↗

U-Pb isotope systematics and age of uranium mineralization, Midnite mine, Washington.

Uranium ores at the Midnite mine, near Spokane, Washington, occur in phyllites and calcsilicates of the Proterozoic Togo Formation, near the margins of an anomalously uraniferous, porphyritic quartz monzonite of Late Cretaceous age. The present geometry of the ore zones is tabular, with the thickest zones above depressions in the pluton-country rock contact. Analyses of high-grade ores from the mine define a 207 Pb/ 204 Pb- 235 U/ 204 Pb isochron indicating an age of mineralization of 51.0 + or - 0.5 m.y. This age coincides with a time of regional volcanic activity (Sanpoil Volcanics), shallow intrusive activity, erosion, and faulting. U-Th-Pb isotopic ages of zircons from the porphyritic quartz monzonite in the mine indicate an age of about 75 m.y., hence the present orebodies were formed about 24 m.y. after its intrusion. The 51-m.y. time of mineralization probably represents a period of mobilization and redeposition of uranium by supergene ground waters, perhaps aided by mild heating and ground preparation and preserved by a capping of newly accumulated, impermeable volcanic rocks. It seems most likely that the initial concentration of uranium occurred about 75 m.y. ago, probably from relatively mild hydrothermal fluids in the contact-metamorphic aureole of the U-rich porphyritic quartz monzonite.Pitchblende, coffinitc, pyrite, marcasite, and hisingerite are the most common minerals in the uranium-bearing veinlets, with minor sphalerite and chalcopyrite. Coffinitc with associated marcasite is paragenetically later than pitchblende, though textural and isotopic evidence suggests no large difference in the times of pitchblende and colfinite formation.The U-Pb isotope systematics of total ores and of pitchblende-coffinite and pyrite-marcasite separates show that whereas open system behavior for U and Pb is essentially negligible for large (200-500 g) ore samples, Pb migration has occurred on a scale of 1 to 10 mm (out of pitchblende and coffinite and into pyrite and marcasite). Also, long-term continuous leakage of radioactive daughters of 238 U (probably 222 Rn) has occurred on scales of from approximately 100 mu m approximately 10 cm. The isotopic composition of unsupported radiogenic Pb in pyrite-marcasite seems to depend on the mineralogical microenvironment of the grains, so that the radiogenic Pb in pyrite-marcasite intimately intermixed with pitchblende-coffinite tends to be deficient in 206 Pb, and the radiogenic Pb in pyrite-marcasite in gangue tends to have excess 206 Pb. These systematics probably reflect differences between the average distances of Pb and 222 Rn diffusion since the formation of the ores.

Washington↗

Surface-water quality and suspended-sediment quantity and quality within the Big River Basin, southeastern Missouri, 2011-13

Missouri was the leading producer of lead in the United States—as well as the world—for more than a century. One of the lead sources is known as the Old Lead Belt, located in southeast Missouri. The primary ore mineral in the region is galena, which can be found both in surface deposits and underground as deep as 200 feet. More than 8.5 million tons of lead were produced from the Old Lead Belt before operations ceased in 1972. Although active lead mining has ended, the effects of mining activities still remain in the form of large mine waste piles on the landscape typically near tributaries and the main stem of the Big River, which drains the Old Lead Belt. Six large mine waste piles encompassing more than 2,800 acres, exist within the Big River Basin. These six mine waste piles have been an available source of trace element-rich suspended sediments transported by natural erosional processes downstream into the Big River. A study was performed by the U.S. Geological Survey in cooperation with U.S. Environmental Protection Agency, Region 7, to calculate and characterize suspended-sediment quantity and quality within the Big River basin after reclamation of the mine waste piles ended in 2012. Streamflow and suspended sediments were quantified and sampled at two locations along a 68-mile reach of the Big River between Bonne Terre and Byrnes Mill, Missouri. The results will help regulatory agencies, such as the U.S. Environmental Protection Agency and U.S. Fish and Wildlife Service, determine impaired reaches and ecosystems for remedial and restoration efforts. Continuous stream stage, water temperature, and turbidity, and discrete suspended-sediment concentration data were collected at the two sites between October 2011 and September 2013. Suspended-sediment samples were collected during various hydrologic conditions to develop a regression model between discrete suspended-sediment concentration and continuous turbidity. Suspended sediments collected during stormflow events were analyzed for concentrations of trace elements such as barium, cadmium, lead, and zinc within two sediment size fractions. Event loads and annual loads of suspended sediment and select trace elements in suspended sediments also were calculated. Suspended-sediment loads computed by the regression model increased downstream from about 201,000 tons at the upstream site to about 355,000 tons at the downstream site during the study period. Stormflow-event-based (hereinafter referred to as “event-based”) suspended-sediment loads ranged from 180 to 32,000 tons at the upstream sampling site and 390 to 53,000 tons at the downstream site along the Big River. Although only seven stormflow events at the upstream site and six at the downstream site were sampled, the event-based suspended-sediment loads accounted for nearly 30 percent of the total suspended-sediment loads computed at both sites, indicating most of the suspended sediment transported through the Big River occurs during higher streamflows. Sediment quality guidelines, known as the threshold effect concentration and the probable effect concentration, used to assess toxicity of trace-element concentrations in sediments were compared to the cadmium, lead, and zinc concentrations in suspended sediment samples collected during stormflow events. All concentrations of cadmium, lead, and zinc in event-based suspended sediment samples exceeded the threshold and probable effect concentrations. Lead and zinc concentrations in the sediment size fraction less than 0.063 millimeters also exceeded the toxic effect threshold, above which sediment is considered to be heavily polluted causing adverse effects on sediment-dwelling organisms. Concentrations of cadmium and zinc in event-based suspended sediment samples were notably higher in samples from the upstream site than samples from the downstream site, indicating the sources of sediments enriched in these trace elements decrease in the downstream area of the watershed. The reduction in concentration of cadmium and zinc could be from dissolution of the constituents during transport or possibly a decrease in downstream source material. The lead concentration exceedance of the probable effects concentration as well as the threshold effects concentration indicates that lead-rich suspended sediments in the fraction less than 0.063 millimeters are readily available within the Big River Basin for transport. These sediments remain in the system from historical mining, and as the reclamation of mine waste piles in the upstream area of the watershed reduce additional sediment loadings, these fine sediments may be continually released as the river scours the streambed and erodes stream banks causing the lead-rich suspended sediment to remain in a state of equilibrium. Barium concentrations in suspended-sediments were nearly twice as high in stormflow event samples collected at the downstream site as compared to samples from the upstream site. The source of barium in the Big River could be from Mineral Fork and Mill Creek, which flow through the historical barite (barium sulfate, also known as tiff) mining district in Washington County, and discharge into the Big River between the two study sites. Total trace-element loads and yields in suspended sediments were computed from the sampled events for each year in the study. The total barium loads in suspended sediments were higher for sampled events collected at the downstream site than the upstream site during both study years. Cadmium and zinc loads in suspended sediments were lower at the downstream site than the upstream site, although the decrease in total load was not substantial during the study period. Lead loads in suspended sediments were lower at the downstream site during the first study year, with a slightly higher load downstream in the second year though the increase from upstream to downstream was small. Event-based yields were higher at the upstream site, indicating that readily available sediment sources are closer to the upstream site where more mining affected areas are located. The estimates determined during large precipitation events indicate that large sources of suspended sediments with large concentrations of trace elements are still available for transport within the Big River.

Missouri↗

Crustal study of a continental strip from the Atlantic Ocean to the Rocky Mountains

Twenty aeromagnetic profiles over a 100-mile-wide strip along the arc of a great circle passing through Denver, Colorado, and Washington, D. C, reveal large anomalies of major crustal significance. Contoured data disclose several areas of distinct magnetic patterns reflecting basement lithology and structure. The mafic rocks of the Blue Ridge and Piedmont and the Keweenawan mafic belt in Iowa and Nebraska give rise to strong linear trends. Areas with a more random pattern of closely spaced magnetic anomalies appear in central Ohio, eastern Iowa, and central and western Nebraska. Except in the Blue Ridge and Piedmont areas, crystalline basement rocks are covered by a thick blanket of virtually nonmagnetic sedimentary rocks, and lithology must be inferred from correlations of the magnetic data with scattered drill-hole data and regional gravity data. The area of highly magnetic rocks in central Ohio has a sharp western boundary that coincides with the western limit of metamorphism associated with the probable extension into Ohio of the Grenville province of Canada. This area and a similar one in eastern Iowa are linked by an arcuate, nearly continuous belt of positive gravity anomalies that extends north into Wisconsin and northern Michigan and then swings southeast across central Michigan. This horseshoe-shaped feature is associated with lithologically diverse but highly magnetic basement rocks. A group of linear magnetic anomalies in western Iowa and eastern Nebraska correlates with the well-known midcontinent gravity high. In Nebraska the magnetic data provide a basis for grouping the extremely complex drill-hole data into three over all lithologic terranes. An analysis of the long-wavelength (>40 miles) variations of the profiles shows that they form a number of large coherent anomalies, many of which show little relation to the major tectonic trends and lithologic patterns of the basement surface. The very broad and less numerous anomalies in the east, which have more or less north-south trends, are significantly different from the more numerous anomalies in the western part of the strip , which tend to trend east-west. One linear anomaly extends for nearly 500 miles across Nebraska and Iowa and may mark a zone of rifting. Heat-flow data show that rocks at the Curie point, which determines the depth below which rock magnetization cannot occur, may be deep enough, at least in shield and other stable parts of the continent, to include a part of the upper mantle. The concentration of the large magnetic features in Iowa and Nebraska may indicate that the thickness of magnetized rock is greater in this area and that perhaps some of these features originate in the upper mantle.

Geological Society of America Bulletin↗

Summary of November 2010 meeting to evaluate turbidite data for constraining the recurrence parameters of great Cascadia earthquakes for the update of national seismic hazard maps

This report summarizes a meeting of geologists, marine sedimentologists, geophysicists, and seismologists that was held on November 18–19, 2010 at Oregon State University in Corvallis, Oregon. The overall goal of the meeting was to evaluate observations of turbidite deposits to provide constraints on the recurrence time and rupture extent of great Cascadia subduction zone (CSZ) earthquakes for the next update of the U.S. national seismic hazard maps (NSHM). The meeting was convened at Oregon State University because this is the major center for collecting and evaluating turbidite evidence of great Cascadia earthquakes by Chris Goldfinger and his colleagues. We especially wanted the participants to see some of the numerous deep sea cores this group has collected that contain the turbidite deposits. Great earthquakes on the CSZ pose a major tsunami, ground-shaking, and ground-failure hazard to the Pacific Northwest. Figure 1 shows a map of the Pacific Northwest with a model for the rupture zone of a moment magnitude M w 9.0 earthquake on the CSZ and the ground shaking intensity (in ShakeMap format) expected from such an earthquake, based on empirical ground-motion prediction equations. The damaging effects of such an earthquake would occur over a wide swath of the Pacific Northwest and an accompanying tsunami would likely cause devastation along the Pacifc Northwest coast and possibly cause damage and loss of life in other areas of the Pacific. A magnitude 8 earthquake on the CSZ would cause damaging ground shaking and ground failure over a substantial area and could also generate a destructive tsunami. The recent tragic occurrence of the 2011 M w 9.0 Tohoku-Oki, Japan, earthquake highlights the importance of having accurate estimates of the recurrence times and magnitudes of great earthquakes on subduction zones. For the U.S. national seismic hazard maps, estimating the hazard from the Cascadia subduction zone has been based on coastal paleoseismic evidence of great earthquakes over the past 5,000 years. The instrumental catalog of earthquakes is of little use for constraining the hazard of the CSZ, because there are virtually no recorded earthquakes on most of the plate interface of the CSZ. There are no historical accounts in the past 150 years of large earthquakes on most of the CSZ. Until about 20 years ago, some interpreted this lack of recent and historical earthquakes as an indicator that the subduction zone was slipping aseismically and could not produce a great earthquake. The work of Brian Atwater and others, in the late 1980s and the 1990s (Atwater, 1987, 1992; Atwater and others, 1995; Nelson and others, 1996; Clague, 1997; Atwater and Hemphill-Haley, 1997; Atwater and others, 2004) demonstrated that submerged forests, buried soils, tsunami deposits, and liquefaction along and near the coast were compelling evidence of repeated great Cascadia earthquakes over at least the past 5,000 years. Atwater and Hemphill-Haley (1997) concluded from paleoseismic evidence at Willapa Bay, Washington, that great earthquakes ruptured the CSZ with an average recurrence time of about 500 years. The date of the last great CSZ earthquake, January 26, 1700, was established from historical records of the so-called orphan tsunami in Japan that is inferred to have been produced by this earthquake (Satake and others, 1996, 2003; Atwater and others, 2005) and is consistent with tree-ring data from drowned forests in Washington and Oregon. From modeling the observations of the tsunami, Satake and others (2003) estimated a moment magnitude of about 9.0 for this earthquake. Many other paleoseismic sites have been investigated along the Pacific Northwest coast from Vancouver Island to northern California and show evidence of great CSZ earthquakes. Nelson and others (2006) summarized the dates found from these studies and proposed correlations between sites indicating the extent of rupture for individual events. Dating of inferred tsunami deposits in Bradley Lake, Oregon by Kelsey and others (2005), as well as tsunami and subsidence evidence from Six Rivers, Oregon (Kelsey and others, 2002) and Coquille River (Witter and others, 2003), indicates that there were probably M w 8 ruptures in the southern portion of the CSZ in addition to the M w 9 events that rupture the whole length of the CSZ (Nelson and others, 2006). A parallel development over the past 20 years or more is the use of deep-sea turbidite deposits for identifying and dating great Cascadia earthquakes over the past 10,000 years (Adams, 1990; Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Turbidites are sediment deposits in the deep ocean from turbidity currents, which are energetic flows of sediment and water along the continental shelf and slope. Adams (1990), using the counts of turbidites in deep-sea cores off the coast of Oregon and Washington collected and analyzed by Griggs (1969) and Griggs and others (1969), proposed that these turbidites were caused by the shaking of great Cascadia earthquakes. Part of his reasoning was that the number (13) of turbidite deposits that occurred since deposition of the Mazama Ash 7,000 years ago gave a recurrence time of about 500 years, consistent with that derived from the coastal submergence data. Adams (1990) also proposed the “confluence test” which evaluates the number of turbidites for submarine channels that form a confluence. He reported that the number of turbidites in the single downstream channel equaled the number in each of the tributary channels. He reasoned that this indicated that the turbidites in each tributary were simultaneously triggered and were, therefore, caused by a common forcing agent. He concluded that shaking from extended ruptures of great Cascadia earthquakes was the most likely cause of these turbidites. Based on the paleoseismic evidence of past great earthquakes, the hazard from the Cascadia subduction zone was included in the 1996 U.S. NSHM (Frankel and others, 1996), which were the basis for seismic provisions in the 2000 International Building Code. These hazard maps used the paleoseismic studies to constrain the recurrence rate of great CSZ earthquakes. Goldfinger and his colleagues have since collected many more deep ocean cores and done extensive analysis on the turbidite deposits that they identified in the cores (Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Using their dating of the sediments and correlation of features in the logs of density and magnetic susceptibility between cores, they developed a detailed chronology of great earthquakes along the CSZ for the past 10,000 years (Goldfinger and others, in press). These correlations consist of attempting to match the peaks and valleys in logs of density and magnetic susceptibility between cores separated, in some cases, by hundreds of kilometers. Based on this work, Goldfinger and others (2003, 2008, in press) proposed that the turbidite evidence indicated the occurrence of great earthquakes (M w 8) that only ruptured the southern portion of the CSZ, as well as earthquakes with about M w 9 that ruptured the entire length of the CSZ. For the southernmost portion of the CSZ, Goldfinger and others (in press) proposed a recurrence time of M w 8 or larger earthquakes of about 230 years. This proposed recurrence time was shorter than the 500 year time that was incorporated in one scenario in the NSHM’s. It is important to note that the hazard maps of 1996 and later also included a scenario or set of scenarios with a shorter recurrence time for M w 8 earthquakes, using rupture zones that are distributed along the length of the CSZ (Frankel and others, 1996; Petersen and others, 2008). Originally, this scenario was meant to correspond to the idea that some of the 500-year averaged ruptures seen in the paleoseismic evidence could have been a series of M w 8 earthquakes that occurred over a short period of time (a few decades), rather than M w 9 earthquakes. Figure 2 shows the logic tree for the CSZ used in the 2008 NSHM’s (Petersen and others, 2008). This logic tree includes whole CSZ rupture earthquakes (M w 8.8–9.2) and partial CSZ rupture earthquakes (M w 8.0–8.7). In this latest version of the NSHM’s, the effective recurrence time of earthquakes on the CSZ with moment magnitudes greater than or equal to 8.0 over the various models is about 270 years (Petersen and others, 2008). This recurrence time applies to the entire CSZ, so that the hazard from great earthquakes was approximately equal along the whole zone, although the hazard estimates taper on the northern and southern ends of the CSZ, because of the way rupture zones of M w 8 earthquakes were distributed along the strike of the CSZ. The NSHM will be updated in 2013, as part of the standard update cycle that corresponds to the update cycle of the national model building codes that are based on the seismic hazard maps. A meeting was necessary to assemble a wide group of experts to hear Dr. Goldfinger explain his methodology for dating and correlating the turbidites and for developing the earthquake chronology. The overall goal of the workshop was to evaluate observations of turbidite deposits to provide constraints on the recurrence times and rupture extents of great Cascadia subduction zone earthquakes for the next update of the NSHM. Before the meeting, participants were supplied with the U.S. Geological Survey (USGS) Professional Paper of Goldfinger and others (in press), as well as material from Brian Atwater and Alan Nelson. The agenda of the meeting was developed by Art Frankel, with assistance from Chris Goldfinger, Brian Atwater, Alan Nelson, Mark Petersen, and Craig Weaver. The meeting was hosted by Chris Goldfinger of Oregon State University. We stress that it is difficult to evaluate in a two-day meeting the large amount of work that Goldfinger and his colleagues have done over the past 15 years or more. This meeting is the first step in a process that develops the inputs to the update of the national maps. The conclusions of this workshop will be discussed and possibly modified at the regional Pacific Northwest workshop for the hazard maps to be held in early 2012. Vetting new research results using informed expert opinion is an integral part of updating the national maps and does not reflect on the veracity of these results.

Cascadia↗

Succession of Laramide magmatic and magmatic-hydrothermal events in the Patagonia Mountains, Santa Cruz County, Arizona

This investigation of the space-time progression of magmatism and hydrothermal activity in the Patagonia Mountains of southern Arizona is based on field and paragenetic relationships, and on U-Pb and 40 Ar/ 39 Ar geochronology of igneous and hydrothermal minerals. The Patagonia Mountains consist of Precambrian, Paleozoic, and Mesozoic sedimentary, granitic, and volcanic rocks, Laramide volcanic rocks, and a core of Laramide intrusions that comprise the Patagonia Mountains batholith. Laramide igneous rocks and adjacent Paleozoic and Mesozoic rocks contain significant porphyry Cu-Mo deposits, Mo-Cu breccia pipes, Ag replacement deposits, and numerous other Cu-Pb-Zn-Ag replacement and vein deposits. Ages of igneous and hydrothermal minerals from 20 U-Pb and 52 40 Ar/ 39 Ar determinations define four magmatic and magmatic-hydrothermal events that formed the batholith and altered parts of it and adjacent rocks; cumulatively the events span at least 16 m.y., from ~74 to 58 Ma. The oldest event of this succession includes the 74 Ma Washington Camp stock and spatially associated Cu-Pb-Zn-Ag replacement deposits in Paleozoic carbonate rocks of the Washington Camp-Duquesne district. Eruption of 73 to 68 Ma volcanic rocks in the northern part of the range was the next youngest event, which coincides temporally with replacement and vein deposits in Paleozoic carbonate rocks at the Flux mine (~71 Ma). An event at 65 to 62 Ma is marked by emplacement of small-volume quartz monzonite, granodiorite, and diorite intrusions, formation of the Ventura breccia deposit in Jurassic granite at 65 to 64 Ma, and formation of other Pb-Zn-Ag-Cu replacement and vein deposits (~62 Ma; Blue Nose and Morning Glory). The Red Mountain porphyry Cu-Mo system is hosted by ~62 Ma granodiorite and Laramide volcanic rocks (73–68 Ma) at the northern end of the batholith. It includes a deep, chalcopyrite-bornite resource (~60 Ma) that is associated with potassic and sericitic alteration and a near-surface chalcocite-enargite resource (60 Ma) that is associated with advanced, supergene-enriched argillic alteration. The youngest event includes the Sunnyside porphyry Cu-Mo system and a Cu-Mo breccia deposit at Red Hill (Four Metals mine), both of which formed in large-volume quartz monzonite, granodiorite, quartz monzonite porphyry, and quartz feldspar porphyry (~61–59 Ma). Similar to the Red Mountain system, the Sunnyside system consists of a deep chalcopyrite resource that occurs in ~60 to 59 Ma quartz feldspar porphyry, and a near-surface, slightly younger (~59–58 Ma) enargite-chalcocite-tennantite resource that occurs in quartz feldspar porphyry, quartz monzonite porphyry, and Mesozoic rocks. The Red Hill Cu-Mo breccia deposit is hosted by large-volume quartz monzonite, granodiorite, and quartz monzonite porphyry (~61–59 Ma). Discrepancies between field and paragenetic relationships and some analytic ages at Sunnyside and Red Hill preclude precise dating of mineralization stages, and may reflect disturbance of isotope systems by multiple, co-spatial to juxtaposed intrusive and hydrothermal events, and/or by unrecognized intrusions. Numerous vein and replacement deposits at the northern end of the batholith, including the Hardshell Ag resource and the Three R supergene chalcocite resource, are distal deposits of the Sunnyside and Red Mountain systems. Small, ~61 to 59 Ma Cu-Mo deposits in large-volume intrusions in the southern part of the batholith consist of hydrothermal quartz, biotite, K-feldspar, muscovite, chalcopyrite, and molybdenite. The age span of magmatic and magmatic-hydrothermal events in the Patagonia Mountains, minimally 16 m.y., is comparable to that of certain other magmatic-hydrothermal successions that contain porphyry Cu-Mo systems. Magmatic-hydrothermal events of the Wasatch-Oquirrh igneous trend, Utah, and the Boulder batholith, Montana, both span ~17 m.y. and include the Bingham and Butte porphyry Cu-Mo, vein and replacement deposits, respectively. Plutons and mineral deposits in the Pima district, Arizona, which includes the porphyry Cu-Mo deposits at Sierrita-Esperanza, Mission-Pima-San Xavier North, and Twin Buttes, formed over an interval of ~14 m.y. The diversity of igneous and hydrothermal products likely reflects evolutionary processes occurring at multiple sites in the lithosphere and at different time scales from >10 m.y. to less than the geochronologic precision currently achievable.

Arizona↗

Application of the Stream Salmonid Simulator (S3) model to assess fall Chinook salmon (Oncorhynchus tshawytscha) production in the American River, California

Executive Summary Anadromous fish returning to the lower American River are restricted to 36 kilometers of free-flowing river between Nimbus Dam and American River’s confluence with the Sacramento River, California. Salmon in the American River provide an important freshwater recreational fishery. However, annual salmon production in the American River in recent years has been low relative to the mid-1990s (Surface Water Resources, Inc., 2001). To investigate the low production of fall-run Chinook salmon ( Oncorhynchus tshawytscha ), the Bureau of Reclamation requested that the U.S. Geological Survey apply the Stream Salmonid Simulator (S3) model to the population of fall-run Chinook salmon on the American River. The American River was chosen among seven candidate Sacramento Basin rivers for S3 application. The American River was selected because of its management and public interest, recently low anadromous fish production, and rich time series of key demographic data needed for S3 application. Data that were not available, however, were empirical estimates on juvenile salmon habitat suitability in the American River. Therefore, a large component of applying S3 to the American River was devoted to the estimation of juvenile salmon habitat suitability and capacity. This entailed snorkeling the lower American River for 3 weeks in March 2021 during the early out-migration period for juvenile Chinook salmon. These efforts were fruitful and showed that the typically small fish (<55 millimeters) in the American River preferred much shallower depths than predicted by habitat suitability criteria derived from the literature for this population. Having empirical estimates on juvenile salmon in the American River provided a solid foundation from which to simulate the population using the S3 model. The S3 model is a spatially explicit population model that runs on a daily time step to simulate redd superimposition, egg maturation, fry emergence and the subsequent growth, survival, and emigration of juvenile Chinook salmon from the river. The key features of this model relevant to this report include (1) a temperature-dependent bioenergetics model driving daily growth rates; (2) density-dependent dynamics that are influenced by the effect of flow on suitable habitat area; and (3) within-year habitat, river flow, and water temperature effects specific to spawning, egg incubation, and fry, parr, and smolt life stages. We used estimates of spawning escapement and geo-referenced redd locations to quantify the spatial and temporal distribution of female spawners for brood years 2014–19. These estimates of female spawners initiate the simulation of each year’s juvenile salmon emergence and emigration over a spatial domain extending from Nimbus Dam to the river’s confluence with the Sacramento River. Using weekly estimates of juvenile salmon abundance and size (fork length) that passed the Watt Avenue fish trap (river kilometer 14.7), we calibrated the S3 model by estimating three key demographic parameters for each year, y : (1) S y , the average daily survival probability, (2) M 0y , the intercept for density-dependence in movement, representing the average daily probability of remaining in a habitat at zero abundance, and (3) C y , the average daily proportion of maximum consumption. These parameters were obtained by minimizing the Mallow’s distance (Lupu and others, 2017) between distributions of weekly abundances and sizes of fish at the traps and weekly simulated abundances and sizes (by S3). Investigation of model fit showed excellent agreement between simulated annual abundances and the abundance of fish passing the fish trap. However, when we compared weekly abundances at the fish trap, S3 under-predicted peaks and over-predicted troughs in the time series of weekly abundances at the fish trap. Thus, some unknown within-year effects have yet to be identified and incorporated in the S3 model. Identifying these important effects and incorporating them in the S3 model would help explain the lack of fit between estimated and simulated weekly abundances. We estimated parameters for 6 years that included a wide range of female spawner abundances (3,057–10,753) and water year types (Critical–Wet). We contrast our estimated parameters to the corresponding number of female spawners and the water year type for the Sacramento Valley. By happenstance, years having higher annual spawner abundances concurred with Critical to Dry water year types. Estimates of survival trended lower with higher spawner abundances and Critical to Dry conditions. In contrast, the extremely wet water year of 2017 had the lowest M 0y , suggesting less density-dependence in fish movement, and the lowest C y , suggesting lower average consumption in this year. When this high-flow year was excluded, a trend towards higher probabilities of fish remaining in a habitat at low abundance and lower proportions of maximum consumption was apparent from Critical to Wet conditions, but only 5 years of data were included. Except for 2017, daily proportions of maximum consumption were relatively high ( C y > 0.83), suggesting that fish were feeding at reasonably high proportions relative to the expected maximum consumption as defined by the “Wisconsin” bioenergetics model (Stewart and Ibarra, 1991). Survival estimates from fry emergence to outmigration at the Sacramento River confluence were generally low when integrated over time. The highest daily survival probability was S y = 0.93 in 2019, or 50 percent total mortality after 10 days. In contrast, our lowest daily survival probability was S y = 0.74 in 2015, or 95 percent total mortality after 10 days. Consequently, even our highest estimated daily survival probability might be considered low. This is especially true given that S y was estimated over a relatively short distance (<14.7 kilometers) from emergence to the Watt Avenue fish trap. Several factors, including our assumed and relatively high daily egg survival rate of 0.9975, could influence juvenile survival estimates. For example, an egg survival rate of 0.9975 results in 3-percent total mortality after 10 days. Egg mortality estimates used in S3 calibration were approximated from egg survivorship studies in the Yakima River, Washington (Johnson and others, 2012), and remains one of the greater uncertainties in S3 when estimating survival across life stages. By including bona fide estimates of egg survival in S3 simulations, the validity of the S3’s current daily egg survival rate could be assessed specifically for the American River. Tagging studies also could provide S3 with direct estimates of juvenile survival and movement; survival during egg incubation then could be estimated indirectly via model fitting.

California↗

Ground Water Atlas of the United States: Segment 11, Delaware, Maryland, New Jersey, North Carolina, Pennsylvania, Virginia, West Virginia

Segment 11 consists of the States of Delaware, Maryland, New Jersey, North Carolina, West Virginia, and the Commonwealths of Pennsylvania and Virginia. All but West Virginia border on the Atlantic Ocean or tidewater. Pennsylvania also borders on Lake Erie. Small parts of northwestern and north-central Pennsylvania drain to Lake Erie and Lake Ontario; the rest of the segment drains either to the Atlantic Ocean or the Gulf of Mexico. Major rivers include the Hudson, the Delaware, the Susquehanna, the Potomac, the Rappahannock, the James, the Chowan, the Neuse, the Tar, the Cape Fear, and the Yadkin-Peedee, all of which drain into the Atlantic Ocean, and the Ohio and its tributaries, which drain to the Gulf of Mexico. Although rivers are important sources of water supply for many cities, such as Trenton, N.J.; Philadelphia and Pittsburgh, Pa.; Baltimore, Md.; Washington, D.C.; Richmond, Va.; and Raleigh, N.C., one-fourth of the population, particularly the people who live on the Coastal Plain, depends on ground water for supply. Such cities as Camden, N.J.; Dover, Del.; Salisbury and Annapolis, Md.; Parkersburg and Weirton, W.Va.; Norfolk, Va.; and New Bern and Kinston, N.C., use ground water as a source of public supply. All the water in Segment 11 originates as precipitation. Average annual precipitation ranges from less than 36 inches in parts of Pennsylvania, Maryland, Virginia, and West Virginia to more than 80 inches in parts of southwestern North Carolina (fig. 1). In general, precipitation is greatest in mountainous areas (because water tends to condense from moisture-laden air masses as the air passes over the higher altitudes) and near the coast, where water vapor that has been evaporated from the ocean is picked up by onshore winds and falls as precipitation when it reaches the shoreline. Some of the precipitation returns to the atmosphere by evapotranspiration (evaporation plus transpiration by plants), but much of it either flows overland into streams as direct runoff or enters streams as base flow (discharge from one or more aquifers). The distribution of average annual runoff (fig. 2) is similar to the distribution of precipitation; that is, runoff is generally greatest where precipitation is greatest. Runoff rates range from more than 50 inches per year in parts of western North Carolina to less than 12 inches in parts of North Carolina, Virginia, and West Virginia. Parts of the seven following physiographic provinces are in Segment 11: the Coastal Plain, the Piedmont, the Blue Ridge, the New England, the Valley and Ridge, the Appalachian Plateaus, and the Central Lowland. The provinces generally trend northeastward (fig. 3). The northeastern terminus of the Blue Ridge Province is in south-central Pennsylvania, and the southwestern part of the New England Province, the Reading Prong, ends in east-central Pennsylvania. The topography, lithology, and water-bearing characteristics of the rocks that underlie the Blue Ridge Province and the Reading Prong are similar. Accordingly, for purposes of this study, the hydrology of the Reading Prong is discussed with that of the Blue Ridge Province. The Coastal Plain Province is a lowland that borders the Atlantic Ocean. The Coastal Plain is as much as 140 miles wide in North Carolina but narrows northeastward to New Jersey where it terminates in Segment 11 at the south shore of Raritan Bay. Although it is generally a flat, seaward-sloping lowland, this province has areas of moderately steep local relief, and its surface locally reaches altitudes of 350 feet in the southwestern part of the North Carolina Coastal Plain. The Coastal Plain mostly is underlain by semiconsolidated to unconsolidated sediments that consist of silt, clay, and sand, with some gravel and lignite. Some consolidated beds of limestone and sandstone are present. The Coastal Plain sediments range in age from Jurassic to Holocene and dip gently toward the ocean. The boundary between the Coastal Plain and the Piedmont Provinces is called the Fall Line (fig. 3) because falls and rapids commonly form where streams cross the contact between the consolidated rocks of the Piedmont (fig. 4) and the soft, semiconsolidated to unconsolidated sediments of the Coastal Plain. The increase in stream gradient at the Fall Line provided favorable locations for mills and other installations that harnessed water power during the early years of the Industrial Revolution, and on most major rivers, the Fall Line coincides with the head of navigation. The Piedmont Province is an area of varied topography that ranges from lowlands to peaks and ridges of moderate altitude and relief. The metamorphic and igneous rocks of this province range in age from Precambrian to Paleozoic and have been sheared, fractured, and folded. Included in this province, however, are sedimentary basins that formed along rifts in the Earth's crust and contain shale, sandstone, and conglomerate of early Mesozoic age, interbedded locally with basaltic lava flows and minor coal beds. The sedimentary rocks and basalt flows are intruded in places by diabase dikes and sills. The mountain belt of the Blue Ridge Province forms the northwestern margin of the Piedmont in most of Segment 11. This belt consists mostly of igneous and high-rank metamorphic rocks but also includes low-rank metamorphic rocks of late Precambrian age and small areas of sedimentary rocks of Early Cambrian age along its western margin. In this report, the Reading Prong of the New England Province, which is an upland that extends from east of the Susquehanna River in Pennsylvania northeastward into New Jersey (fig. 3), is treated as part of the Blue Ridge Province. Part of the Reading Prong in Pennsylvania and New Jersey and a small part of the Piedmont Province in northeastern New Jersey have been glaciated. Glacial deposits completely or partly fill some of the valleys, and the eroding action of the glacial ice removed some of the rock from the ridges. Thus, the glaciated parts of the province have a smoother topography and less relief than other parts. The Valley and Ridge Province is characterized by layered sedimentary rock that has been complexly folded and locally thrust faulted. As the result of repeated cycles of uplift and erosion, resistant layers of well-cemented sandstone and conglomerate form elongate mountain ridges and less resistant, easily eroded layers of limestone, dolomite, and shale form valleys. The rocks of the province range in age from Cambrian to Pennsylvanian. Parts of this province from central Pennsylvania into New Jersey have been glaciated, and glacial deposits fill or partially fill some of the valleys. The Appalachian Plateaus Province is underlain by rocks that are continuous with those of the Valley and Ridge Province, but in the Appalachian Plateaus the layered rocks are nearly flat-lying or gently tilted and warped, rather than being intensively folded and faulted. The boundary between the two provinces is a prominent southeast-facing scarp called the Allegheny Front in most of the northern part of Segment 11 (fig_ 5) and the Cumberland Escarpment in the southern part. The scarp faces the Valley and Ridge Province, and throughout most of the segment, the eastern edge of the Appalachian Plateaus Province is higher than the ridges in the Valley and Ridge. Like parts of the Reading Prong and the Valley and Ridge Province, the northern part of the Appalachian Plateaus Province in Pennsylvania has been glaciated. In the glaciated section, the surface is mantled by glacial drift, and the valleys are partly filled with glacial deposits. The northwestern corner of Segment 11 contains a small part of the Central Lowland Province. This flat lowland is underlain by gently dipping sedimentary rocks, some of which are the same geologic formations as those of the Appalachian Plateaus Province. The two provinces are separated by a northwest- facing scarp. Because of the small area of the Central Lowland Province within the segment and the similarity of aquifer properties with those of the glaciated part of the Appalachian Plateaus Province, the two provinces are discussed together in this report.

Delaware, Maryland, New Jersey, North Carolina, Pe↗