SEARCH · Search USGS

Results for “Water Information Series Report”

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

SourceResearch connectionFrom yearThrough year

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

792 records · Page 44Linked to original sources

USGS2022

Estimating stream temperature in the Willamette River Basin, northwestern Oregon—A regression-based approach

The alteration of thermal regimes, including increased temperatures and shifts in seasonality, is a key challenge to the health and survival of federally protected cold-water salmonids in streams of the Willamette River basin in northwestern Oregon. To better support threatened fish species, the U.S. Army Corps of Engineers (USACE) and other water managers seek to improve the thermal regime in the Willamette River and key tributaries downstream of USACE dams by utilizing strategically timed flow releases from USACE dams. To inform flow management decisions, regression relations were developed for 12 Willamette River basin locations below USACE dams relating stream temperature with streamflow and air temperature utilizing publicly available datasets spanning 2000–18. The resulting relations provide simple tools to investigate stream temperature responses to changes in streamflow and climatic conditions in the Willamette River system. Regression relations on the Willamette River and key tributaries show that, at locations sufficiently distant from the direct temperature influence of upstream dam releases, air temperature and streamflow are reasonable proxies to predict the 7-day average of the daily mean (7dADMean) and 7-day average of the daily maximum (7dADMax) water temperature with errors generally ≤1 degrees Celsius (°C). To account for seasonal variations in the relation between air temperature, streamflow, and stream temperature, a transition-smoothed, seasonal regression approach was used. Stream temperature is inversely correlated with streamflow in all seasons except “winter” (January–March), when it is relatively independent. Stream temperature is positively correlated with air temperature in all seasons, but the slope decreases at very low or very high air temperatures. Generally, fit is best for seasonal models “winter” (January–March), “spring” (April–May), “summer” (June–August), and “early autumn” (September–October). Error in “autumn” (November–December) is larger, probably due to variation in the onset timing of winter storms. Simulated results from a climatological analysis of predicted stream temperature suggest that, excluding extremes and accounting for some seasonal variability, the 7dADMean and 7dADMax stream temperature sensitivity to air temperature and streamflow varies by location on the river. To investigate the potential range of stream temperature variability based on historical air temperature and streamflow conditions, stream temperature predictions were calculated using synthetic time series comprised of daily temperature values representing the 0.10, 0.33, 0.50, 0.67, and 0.90 quantile of air temperature and streamflow from 1954 (the year meaningful streamflow augmentation began) to 2018. Results show that from a “very hot” (0.90 quantile) and “very dry” (0.10 quantile) year to a “very cool” (0.10 quantile) and “very wet” (0.90; all quantiles from 1954 to 2018) year, the stream temperature sensitivity to air temperature and streamflow is about 3 °C at Harrisburg (river mile 161.0) and increases to about 5 °C at Keizer (river mile 82.2). While the number of days exceeding regulatory criteria are fewer in cooler, wetter years than in warmer, dryer years, the models suggest that the Willamette River will likely continue to exceed the State of Oregon maximum water-temperature criterion of 18 °C for sustained periods from late spring to early autumn and that the flow management practices evaluated in this study, while effective at influencing stream temperature, likely cannot prevent many or all such exceedances. As modeled for 2018, a representative very hot year with normal to below-normal streamflow, stream temperature sensitivity to changes in streamflow of ±100 to ±1000 cubic feet per second produced mean monthly temperature changes from 0.0 to 1.4 °C at Keizer, Albany, and Harrisburg during summer. For a specified change in flow, temperature sensitivity is greater at upstream locations where streamflow is less than that at downstream locations because the change in streamflow is a greater percentage of total streamflow at upstream locations. Similarly, temperature response to a set change in flow is greater in the summer and early autumn low-flow season than in spring when flows are higher. The regression models developed in this study thus indicate that flow management is likely to have a greater effect on stream temperature at upstream locations (such as Harrisburg or Albany) and during the low-flow season than at downstream locations (such as Keizer) or during periods of higher streamflow.

Oregon↗

USGS2018

Analysis of groundwater response to tidal fluctuations, Operable Unit 2, Area 8, Naval Base Kitsap, Keyport, Washington

Operable Unit 2, Area 8, at Naval Base Kitsap, Keyport is the site of a former chrome-plating facility that released metals (primarily chromium and cadmium), chlorinated volatile organic compounds, and petroleum compounds into the local environment. To ensure long-term protectiveness, as stipulated in the Fourth Five-Year Review for the site, Naval Facilities Engineering Command Northwest collaborated with the U.S. Environmental Protection Agency, the Washington State Department of Ecology, and the Suquamish Tribe, to collect data to monitor the contamination left in place and to ensure the site does not pose a risk to human health or the environment. To support these efforts, refined information was needed on the interaction of fresh groundwater with seawater in response to the up-to 13-ft tidal fluctuations at this nearshore site adjacent to Port Orchard Bay. The information was analyzed to meet the primary objective of this investigation, which was to determine the optimal time during the semi-diurnal and the neap-spring tidal cycles to sample groundwater for freshwater contaminants in Area 8 monitoring wells. Groundwater levels and specific conductance in five monitoring wells, along with marine water-levels (tidal levels) in Port Orchard Bay, were monitored every 15 minutes during a 3-week duration to determine how nearshore groundwater responds to tidal forcing. Time series data were collected from October 24, 2017, to November 16, 2017, a period that included neap and spring tides. Vertical profiles of specific conductance were also measured once in the screened interval of each well prior to instrument deployment to determine if a freshwater/saltwater interface was present in the well during that particular time. The vertical profiles of specific conductance were measured only one time during an ebbing tide at approximately the top, middle, and bottom of the saturated thickness within the screened interval of each well. The landward-most well, MW8-8, was completely freshwater, while one of the most seaward wells, MW8-9, was completely saline. A distinct saltwater interface was measured in the three other shallow wells (MW8-11, MW8-12, and MW8-14), with the topmost groundwater occurring fresh underlain by higher conductivity water. Lag times between minimum spring-tide level and minimum groundwater levels in wells ranged from about 2 to 4.5 hours in the less-than 20-ft deep wells screened across the water table, and was about 7 hours for the single 48-ft deep well screened below the water table. Those lag times were surprisingly long considering the wells are all located within 200-ft of the shoreline and the local geology is largely coarse-grained glacial outwash deposits. Various manmade subsurface features, such as slurry walls and backfilled excavations, likely influence and confuse the connectivity between seawater and groundwater. The specific-conductance time-series data showed clear evidence of substantial saltwater intrusion into the screened intervals of most shallow wells. Unexpectedly, the intrusion was associated with the neap part of the tidal cycle around November 13–16, when relatively low barometric pressure and high southerly winds led to the highest high and low tides measured during the monitoring period. The data consistently indicated that the groundwater had the lowest specific conductance (was least mixed with seawater) during the prior neap tides around October 30, the same period when the shallow groundwater levels were lowest. Although the specific conductance response is somewhat different between wells, the data do suggest that it is the heights of the actual high-high and low-low tides, regardless of whether or not they occur during the neap or spring part of the cycle, that allows seawater intrusion into the nearshore aquifer at Area 8. With all the data taken into consideration, the optimal time for sampling the shallow monitoring wells at Area 8 would be centered on a 2–5-hour period following the predicted low-low tide during neap tide, with due consideration of local atmospheric pressure and wind conditions that have the potential to generate tides that can be substantially higher than those predicted from lunar-solar tidal forces. The optimal time for sampling the deeper monitoring wells at Area 8 would be during the 6–8-hour period following a predicted low-low tide, also during the neap tide part of the tidal cycle. The specific time window to sample each well following a low tide can be found in table 5. Those periods are when groundwater in the wells is most fresh and least diluted by seawater intrusion. In addition to timing, consideration should be given to collecting undisturbed samples from the top of the screened interval (or top of the water table if below the top of the interval) to best characterize contaminant concentrations in freshwater. A downhole conductivity probe could be used to identify the saltwater interface, above which would be the ideal depth for sampling.

Washington↗

USGS1914

The Transportation of Debris by Running Water

Scope.-The finer debris transported by a stream is borne in suspension. The coarser is swept along the channel bed. The suspended load is readily sampled and estimated, and much is known as to its quantity. The bed load is inaccessible and we are without definite information as to its amount. The primary purpose of the investigation was to learn the laws which control the movement of bed load, and especially to determine how the quantity of load is related to the stream's slope and discharge and to the degree of comminution of the debris. Method.-To this end a laboratory was equipped at Berkeley, Cal., and experiments were performed in which each of the three conditions mentioned was separately varied and the resulting variations of load were observed and measured. Sand and gravel were sorted by sieves into grades of uniform size. Determinate discharges were used. In each experiment a specific load was fed to a stream of specific width and discharge, and measurement was made of the slope to which the stream automatically adjusted its bed so as to enable the current to transport the load. The slope factor.-For each combination of discharge, width, and grade of debris there is a slope, called competent slope, which limits transportation. With lower slopes there is no load, or the stream has no capacity for load. With higher slopes capacity exists; and increase of slope gives increase of capacity. The value of capacity is approximately proportional to a power of the excess of slope above competent slope. If S equal the stream's slope and sigma equal competent slope, then the stream's capacity varies as (S - sigma)n. This is not a deductive, but an empiric law. The exponent n has not a fixed value, but an indefinite series of values depending on conditions. Its range of values in the experience of the laboratory is from 0.93 to 2.37, the values being greater as the discharges are smaller or the debris is coarser. The discharge factor.-For each combination of width, slope, and grade of debris there is a competent discharge, k. Calling the stream's discharge Q, the stream's capacity varies as (Q - k)o. The observed range of values for o is from 0.81 to 1.24, the values being greater as the slopes are smaller or the debris is coarser. Under like conditions o is less than n; or, in other words, capacity is less sensitive to change3 of discharge than to changes of slope. The fineness factor.-For each combination of width, slope, and discharge there is a limiting fineness of debris below which no transportation takes place. Calling fineness (or degree of comminution) F and competent fineness o, the stream's capacity varies with (F - o)p. The observed range of values for p is from 0.50 to 0.62, the values being greater as slopes and discharges are smaller. Capacity is less sensitive to changes in fineness of debris than to changes in discharge or slope. The form factor.-Most of the experiments were with straight channels. A few with crooked channels yielded nearly the same estimates of capacity. The ratio of depth to width is a more important factor. For any combination of slope, discharge, and fineness it is possible to reduce capacity to zero by making the stream very wide and shallow or very narrow and deep. Between these extremes is a particular ratio of depth to width, p, corresponding to a maximum capacity. The values of p range, under laboratory conditions, from 0.5 to 0.04, being greater as slope, discharge, and fineness are less. Velocity.-The velocity which determines capacity for bed load is that near the stream's bed, but attempts to measure bed velocity were not successful. Mean velocity was measured instead. To make a definite comparison between capacity and mean velocity it is necessary to postulate constancy in some accessory condition. If slope be the constant, in which case velocity changes with discharge, capacity varies on the average with the 3.2 power of velocity. If discharge be the constant, in w

Professional Paper↗

USGS2024

Physics to fish—Understanding the factors that create and sustain native fish habitat in the San Francisco Estuary

Executive Summary The Bureau of Reclamation (Reclamation) operates the Central Valley Project (CVP), one of the nation’s largest water projects. Reclamation has an ongoing need to improve the scientific basis for adaptive management of the CVP and, by extension, joint operations with California’s State Water Project. The U.S. Geological Survey (USGS) works cooperatively with the Bureau of Reclamation to provide scientific support for the management of Reclamation’s CVP project. Major habitat restoration efforts and a new water-diversion point are planned to benefit delta smelt ( Hypomesus transpacificus ) and other species of concern while ensuring the reliability of water supply. In addition, various flow actions and management activities have been identified as possible methods to increase populations of delta smelt and salmonid ( Oncorhynchus spp.) runs of concern. The overarching goal of this cooperative project was to provide Reclamation with the scientific information needed to evaluate the efficacy of ongoing and future adaptive management actions and to improve the scientific basis for more flexible CVP operations that would achieve water-supply reliability and fish protection. The research and monitoring described in this report comprises the period 2015–19 and focuses on management issues related to native fish species of concern, especially delta smelt. Conserving the delta smelt population while providing a reliable water supply is a primary management and policy issue in California. Our approach for this cooperative project is based on the “physics to fish” concept, the idea that high-quality habitat is generated and sustained by the interaction between physical processes and the landscape. These interactions create a template for chemical and biological processes that can change across a variety of spatial and temporal scales. Following this concept, this project (hereafter referred to as “the physics to fish project”) included monitoring and studies of water flows, sediments, water quality, and invertebrate and fish dynamics across a range of spatial and temporal scales and in regions relevant to resource managers tasked with managing water supplies and ecosystem health in the San Francisco Estuary. The intent of this approach was to document the habitat conditions, important processes, and interactions among them that create high-quality habitat for native fishes so that the likely effects of future management actions (for example, habitat restoration) can be objectively assessed at the local (site-specific), regional (within subregions of the estuary), and landscape (across the entire estuary and beyond) scales. Hydrodynamically, the upper estuary (landward of Carquinez Strait) is characterized by a fixed volume of tidally exchanged water (for example, tidal prism) that interacts with the existing channel network and bathymetry to create regions with differing hydrodynamics. Our results indicate that careful study of construction or reoperation of existing infrastructure to perform management actions can help (1) improve the accuracy of hydrodynamic models; (2) further understanding of ecological effects; and (3) enhance abilities to predict ecological outcomes. At the local scale, we developed a new concept called the Lagrangian to Eulerian (LE) ratio that can be used as a tool for understanding the importance of various hydrodynamic processes in specific channels or channel networks and for forecasting transport dynamics. Channels with LE ratios<1 in a channel network or in a dead-end slough are hydrodynamically able to develop an exchange zone between two parcels of water that may have different chemical and physical properties. In a dead-end channel, there is a landward region with long residence time (no-exchange zone) and a seaward region with short residence time (high-exchange zone) that are well mixed with seaward waters. At the transition (exchange zone) between the high and no-exchange regions, a gradient will form in water-quality constituents that differ in concentration between the landward and seaward waters. Turbidity affects fish habitat and has declined through time in the San Francisco Estuary. Average turbidity across the Sacramento–San Joaquin Delta (hereafter referred to as “the Delta”) is dependent on annual hydrology. In dry years, the region around Cache Slough (known regionally as the “Cache Slough Complex”) in the northern Delta is generally more turbid than Suisun Bay and the lower Sacramento River. When the Yolo By-Pass (known regionally as “Yolo Bypass”), a large flood bypass that runs parallel to the Sacramento River in the northern Delta, is not flooding and river flows are lower, sediment is usually transported into the Cache Slough Complex because flood tides dominate ebb tides, resulting in transport of suspended sediment from seaward areas of the upper estuary into the Cache Slough Complex. These hydrodynamic conditions also favor the formation of turbidity maximums (TMs) in the Cache Slough Complex. The TMs are areas of higher suspended-sediment concentration, providing higher-turbidity habitat favored by some fishes, including delta smelt, and they can also concentrate other constituents, including phytoplankton and organic carbon that can be important in food webs. Pelagic primary production by phytoplankton is the basis for Delta food webs supporting pelagic fishes such as delta smelt; however, phytoplankton abundance in the Delta has declined during recent decades. We examined how nutrients, hydrodynamics, and other factors affect phytoplankton blooms. Based on our results, we developed three new concepts of phytoplankton bloom formation in the Delta, each associated with a distinct set of hydrologic conditions. First, productivity cascades highlighted how local processes can contribute to phytoplankton blooms observed at the regional scale. Second, we observed phytoplankton blooms in the upper San Francisco Estuary that were associated with transport out of Yolo By-Pass (transport blooms). Third, we also documented a series of phytoplankton blooms that were in the confluence area at the landward edge of Suisun Bay. The conditions leading to creation of confluence phytoplankton blooms are not yet understood, but the confluence region connects the Cache Slough Complex with Suisun Marsh. Therefore, blooms in this area have the potential to spread to large areas of the Delta. At the landscape scale, the distribution of the invasive clams ( Potamocorbula amurensis and Corbicula fluminea , hereafter referred to as “ Corbicula ”) is driven by salinity. At smaller spatial scales, the distribution of either species is sensitive to multiple factors affecting survival and reproduction, complicating efforts to predict distribution and abundance without considering local-scale conditions across the area of interest. In the Cache Slough Complex, the area landward of the exchange zone in regions with LE ratio<1 were characterized by low abundances of Corbicula probably because recruits from seaward areas are not transported past the exchange zone and because there are no landward tributaries with adult Corbicula to provide an upstream source of recruits. Corbicula biomass was highest near or downstream from the exchange zone consistent with Corbicula grazing on phytoplankton produced in the exchange zone or transported from the no-exchange zone. The severity of Corbicula grazing could be reduced by manipulating the hydrodynamic characteristics of waterways; however, the beneficial and harmful effects on the organisms meant to benefit from increased phytoplankton production, including zooplankton and fish species of concern, should be thoroughly examined before manipulating hydrodynamic characteristics. The distribution of fishes at the landscape scale is generally driven by the position of the salinity field in the estuary. The physics to fish project compared distributions of fishes at Ryer Island, a tidal wetland in Suisun Bay and a region of variable salinity, with fish distributions at the Cache Slough Complex, a freshwater region. At Ryer Island, there was an absence of freshwater invasive species and an abundance of native species, such as Sacramento splittail ( Pogonichthys macrolepidotus ), tule perch ( Hysterocarpus traskii ), and Sacramento pikeminnow ( Ptychocheilus grandis ). The native species were almost exclusively captured in wetland and nearshore shallow-water habitat regardless of water-quality conditions. In the Cache Slough Complex, our regional scale objective was to elucidate how hydrodynamic-physical habitat interactions drive fish-community structure. Our studies showed that dendritic channel systems were better able to support native species, while intertidal habitats supported those species best able to exploit the transient character of the habitat. Habitats upstream from the exchange zone were especially important in supporting high numbers of native fishes relative to within or downstream from the exchange zone. Many of the native species were associated with tidal marsh in the no-exchange zone. More pelagic-oriented, mobile species, such as Striped Bass ( Morone saxatilis ), threadfin shad ( Dorosoma petenense ), and Sacramento pikeminnow, were more affected by water-quality conditions, such as turbidity. The physics to fish concept developed in this project provides a framework for designing individual projects and for considering the cumulative effects of multiple projects in a region, using the LE ratio as a guiding metric. The physics to fish concept may also provide a suitable framework for coordinating management actions. Tidal wetlands can function in several ways in the hydrodynamic framework. Relatively small tidal wetlands with short channel networks and with LE ratios>1 are not able to maintain a landward no-exchange zone or an exchange zone. This likely means that any contributions to pelagic food webs would be limited to resources derived from wetland vegetation, which can include dissolved and particulate organic matter (detritus) and populations of consumers that can increase in abundance based on those resources. The fate of the contributed production from these channels depends on the characteristics of the receiving waters seaward of the tidal wetland. If these channels join a large system such as Suisun Bay, then any contribution is likely to be rapidly dispersed in the larger volume; however, the channel junction might provide a focal point for consumers, such as fishes, to congregate and feed on material leaving the wetland on ebb tides before it is dispersed in the larger volume. Fishes might also access these resources by entering the wetland. The physics to fish project has established a foundation and several new concepts for understanding how habitat restoration can benefit native fish populations at the local and regional levels. Many of the ideas regarding habitat restoration and channel modifications outlined in this report could help guide management actions that could improve conditions for native fishes at little or no water cost beyond water already dedicated to other management actions. A complete list of products originating from this work is provided in appendix 1.

California↗

USGS2010

Glaciers of Asia

This chapter is the ninth to be released in U.S. Geological Survey Professional Paper 1386, Satellite Image Atlas of Glaciers of the World, a series of 11 chapters. In each of the geographic area chapters, remotely sensed images, primarily from the Landsat 1, 2, and 3 series of spacecraft, are used to analyze the specific glacierized region of our planet under consideration and to monitor glacier changes. Landsat images, acquired primarily during the middle to late 1970s and early 1980s, were used by an international team of glaciologists and other scientists to study various geographic regions and (or) to discuss related glaciological topics. In each glacierized geographic region, the present areal distribution of glaciers is compared, wherever possible, with historical information about their past extent. The atlas provides an accurate regional inventory of the areal extent of glacier ice on our planet during the 1970s as part of a growing international scientific effort to measure global environmental change on the Earth’s surface. The chapter is divided into seven geographic parts and one topical part: Glaciers of the Former Soviet Union (F–1), Glaciers of China (F–2), Glaciers of Afghanistan (F–3), Glaciers of Pakistan (F–4), Glaciers of India (F–5), Glaciers of Nepal (F–6), Glaciers of Bhutan (F–7), and the Paleoenvironmental Record Preserved in Middle-Latitude, High-Mountain Glaciers (F–8). Each geographic section describes the glacier extent during the 1970s and 1980s, the benchmark time period (1972–1981) of this volume, but has been updated to include more recent information. Glaciers of the Former Soviet Union are located in the Russian Arctic and various mountain ranges of Russia and the Republics of Georgia, Kyrgyzstan, Tajikistan, and Kazakstun. The Glacier Inventory of the USSR and the World Atlas of Ice and Snow Resources recorded a total of 28,881 glaciers covering an area of 78,938 square kilometers (km 2 ). China includes many of the mountain-glacier systems of the world including the Himalaya, Karakorum, Tien Shan and Altay mountain ranges. The glaciers are widely scattered and cover an area of about 59,425 km 2 . The mountain glaciers may be classified as maritime, subcontinental or extreme continental. In Afghanistan, more than 3,000 small glaciers occur in the Hindu Kush and Pamir mountains. Most glaciers occur on north-facing slopes shaded by mountain peaks and on east and southeast slopes that are shaded by monsoon clouds. The glaciers provide vital water resources to the region and cover an area of about 2,700 km 2 . Glaciers of northern Pakistan are some of the largest and longest mid-latitude glaciers on Earth. They are located in the Hindu Kush, Himalaya, and Karakoram mountains and cover an area of about 15,000 km 2 . Glaciers here are important for their role in providing water resources and their hazard potential. The glaciers in India are located in the Himalaya and cover about 8,500 km 2 . The Himalaya contains one of the largest reservoirs of snow and ice outside the polar regions. The glaciers are a major source of fresh water and supply meltwater to all the rivers in northern India, thereby affecting the quality of life of millions of people. In Nepal, the glaciers are located in the Himalaya as individual glaciers; the glacierized area covers about 5,324 km 2 . The region is the highest mountainous region on Earth and includes the Mt. Everest region. Glaciers in the Bhutan Himalaya have a total area of about 1,317 km 2 . Many recent glacier studies are focused on glacier lakes that have the potential of generating dangerous glacier lake outburst floods. Research on the glaciers of the middle-latitude, high-mountain glaciers of Asia has also focused on the information contained in the ice cores from the glaciers. This information helps in the reconstruction of paleoclimatic records, and the computer modeling of global climate change.

Professional Paper↗

USGS2020

Smallmouth buffalo (Ictiobus bubalus) growth across a 1200km human use and ecological disturbance gradient in the Upper Mississippi River System

Smallmouth buffalo (Ictiobus bubalus) is a common and widely distributed large-bodied species of the family Catostomidae. It inhabits large rivers and reservoirs of the eastern continental United States (east of the continental Divide) and is most abundant and common in the large rivers of the Midwest and Central Plains, though it does occur as far north and east as the Hudson Bay drainage and as far south and west as Arizona (Edwards and Twoney 1982). Historically, smallmouth buffalo were an important component of commercial fisheries on both the Mississippi and Illinois Rivers. However, following the introduction of common carp (Cyprinus carpio) in the mid-1800s (Carlander 1954), the construction of a system of navigation dams on Upper Mississippi and Illinois River in the 1930s (USGS 1999), and water quality/pollution issues through the 1980s (Weiner 2010), the role of smallmouth buffalo in the overall UMRS fish community and commercial fishery has generally diminished relative to historical standards. Still, smallmouth buffalo remains an important and valued component of the UMRS commercial fishery. The study area is represented by three study reaches on the Illinois River and three study reaches on the Upper Mississippi River (Figure 1). Collectively, these study reaches represent nearly 1200 river km and exist across strong and pronounced ecological and disturbance gradients. For example, habitat composition, water quality, commercial navigation intensity, aquatic plant prominence, and the number and abundance of nonnative fish species vary strongly across the study domain, with northern Mississippi River reaches exhibiting less navigation traffic, better water quality, markedly greater aquatic plant prominence, more diverse habitat composition, and comparably much smaller numbers of nonnative species than the lower Mississippi River study reach and those on the Illinois River (USGS 1999; Johnson and Hagerty [eds] 2008; Irons et al. 2009). Long term monitoring efforts conducted under the auspices of the Upper Mississippi River Restoration program over the past 27 years have provided tremendous insights into shifts and changes of the overall UMRS fish community (Ickes et al. 2005; Garvey et al. 2010; Schramm and Ickes 2016). However, these monitoring efforts observe only the most basic aspects of the UMRS fish community (i.e., catch, length, weight, distribution, and occurrence). To gain a greater understanding of forces driving community level shifts and changes, more directed study is needed on the functional attributes of fish populations (i.e., growth, mortality, recruitment). Collectively, these functional attributes of populations are termed population dynamics and/or vital rates. It is important to note, the population dynamics of fishes in large rivers is generally poorly understood, especially for non-game species (Ickes 2018). The prevailing view is that abiotic factors largely govern inter-annual population dynamics, typically based upon rather short-term observations and correlations with assorted abiotic river attributes that vary on a seasonal or annual basis (for example, Risotto and Turner 1985). However, the role that longer-term abiotic factors play in regulating population abundance, or that biotic factors internal to the population (e.g., spawner-recruit dynamics, growth dynamics) or external to the population (e.g., predator-prey dynamics, sympatric competitors, disease) remain poorly understood. Achieving a greater understanding of these dynamics is important for stock, game, and invasive species management. In 2017, as part of a larger study designed to gain vital population rate information for smallmouth buffalo in the Upper Mississippi and Illinois Rivers (“Smallmouth Buffalo population demographics of the Upper Mississippi River System”; UMRR LTRM 2018SOW project items 2018MMBF1-2018MMBF6) annual growth patterns in smallmouth buffalo were determined and evaluated. This was accomplished by measuring growth histories recorded in annual growth increments on hard bony parts (here otoliths), a method known generically as biochronology, and somewhat analogous to dendrochronology practiced by foresters. These methods allow one to generate time-series of annual growth histories that depend upon age, year class (i.e., cohort), and annual environmental conditions experienced by the population over time (Weisberg, 1993). Biochronology methods were used to develop a 36-year time series of smallmouth buffalo growth in the Upper Mississippi and Illinois Rivers across a 1200 km ecological and human use disturbance gradient. Annual growth intervals were identified and measured from otoliths to determine fish age and growth history. A mixed model that parses the growth increment into age and year effects was fit to these data. Given the pronounced ecological and disturbance gradients inherent to the UMRS and the study domain, an a priori expectation of differing patterns in growth is accepted as a null hypothesis to test. The goal of this study was to model smallmouth buffalo growth as a function of the age of the fish and the growth year in which the growth was gained. The primary modeling objective was to parse growth observed on each annulus into a portion attributable to the age of the fish and the portion attributable to the year in which the growth was gained. In effect, this modeling approach removes the somewhat trivial age effects on growth so that a non-confounded growth year effect can be gained. Results attributable to growth year provide a time series of growth information that is of the same duration as the oldest fish observed and solely reflects environmental influences on growth. These model responses can then be investigated relative to environmental covariate time-series suspected of influencing growth of smallmouth buffalo in the Upper Mississippi and Illinois Rivers (e.g., temperature, discharge, population density, population mortality, forage availability, sympatric competition, habitat composition, navigation intensity, nonnative fish prominence, etc.). Thus, the primary scientific objective was to investigate if and how smallmouth buffalo growth varies in accordance with innate ecological and disturbance gradients across the study domain.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

USGS2026

Small creatures, big challenge: A review of current knowledge and assessment of research needs to inform water withdrawal management for the Duck River's freshwater mussel fauna

The Duck River is a remarkable national resource, boasting outstanding biodiversity (Ahlstedt et al. 2017). It contains 151 fish species, over 60 mussel species, and 22 snail species, many of which are federally listed (Etnier and Starnes 1993; Johnson et al. 2013; Ahlstedt et al. 2017; Womble and Rosenberger 2025a). The Duck River watershed is at a pivotal time in its management and ecological history due to significant growth in human population and economic activity. The consequent increase in surface water demand presents challenges for sustainable water resource management. The Tennessee Department of Environment and Conservation (TDEC) has permitted or proposed to permit eight water utilities to withdraw over 77 million gallons per day (mgd) from the Duck River, representing > 40% increase from previously permitted amounts (USFWS 2025). On November 20, 2024, Tennessee Governor Bill Lee signed Executive Order No 108 to ensure long-term stability of the Duck River watershed through balancing economic growth, water resource management, and environmental conservation. Some primary components of this plan are the creation of a Comprehensive Regional Drought Management Plan, the development of a Habitat Conservation Plan, and the establishment of the Duck River Watershed Planning Partnership (DRWPP). Although the planning partnership is not focused on a sole aspect of the Duck River, but the watershed as a whole, the mussel assemblages of the Duck River have been a particular focus due to their importance in providing ecosystem services, sensitivity to water quality, conservation status, and sensitivity to periodic drought conditions in this system. With multiple state and federal agencies involved in managing the Duck River’s unique water resources, a focused review on how future water management could affect the freshwater mussel communities in the Duck River will contribute to this process. Our objectives of this review are to: 1. Review the effects of changing flow regimes on mussels, with emphasis on the effect of droughts/low flows, highlighting key processes/factors affecting mussels. 2. Highlight critical data gaps that could allow future, evidence-based water management decisions in this system. 3. Propose research topics to address critical data gaps and enable evidence-based water management decisions that consider the needs of freshwater mussels. Freshwater mussels are a critical component of riverine ecosystems, providing ecosystem services including water filtration, nutrient cycling, and habitat provision (Vaughn et al. 2004; Haag 2012) that benefit human communities. Environmental flows, also referred to as ecological or e-flows, refer to the quantity, timing, and quality of water flows necessary to maintain the ecological health of rivers, wetlands, and other aquatic ecosystems (Acreman and Dunbar 2004; Acreman 2016). For freshwater mussels, environmental flows include maintaining a flow regime and associated water quality regimes that are favorable for both juveniles and adults while also seasonally supporting host fish and the survival of settling glochidia (i.e., larval mussels) in suitable environments (Gates et al. 2015). Freshwater mussels are adapted to the natural flow regime of their native streams, including seasonal variability in flows. As sedentary filter feeders, mussels require adequate water flow to provide oxygen and food (Allen and Vaughn 2009); however, they also require periods of low flow to increase access to fish hosts and to ensure juveniles are not displaced (Gates et al 2015). Mussels are sensitive to water quality conditions, including pollutants, sedimentation, and artificial alterations of temperature and pH (Strayer and Malcom 2012). Mussel community characteristics, including presence, abundance, and species richness, generally increase with stream size (e.g., stream order); however, responses to changes in flow conditions are often species-specific, and the magnitude and timing of extreme flow conditions have more predictive power for mussel community response than general changes in flow (Lopez and Vaughn 2021). Complex hydraulic variables, such as relative shear stress and Reynolds number, can be strong predictors of mussel community responses because they combine simple hydraulic variables (e.g., velocity, depth) with stream channel characteristics, such as substrate composition and dynamics. Droughts directly and indirectly affect mussels by influencing their physiology, behavior, reproduction, and ecosystem roles. Mass mussel die-offs during droughts and low water conditions are widely documented (Vaughn et al. 2015; Sousa et al. 2018; Dubose et al. 2019). Droughts can also lead to a shift in mussel community composition, with drought-tolerant species becoming more dominant while more sensitive species experience declines (Galbraith et al, 2010; Lopez et al. 2022; Vaughn and Atkinson 2025). Droughts affect mussels through direct mortality from stranding as well as physiological stress from altered water temperature and quality, reduced habitat and food availability and quality, increased predation and competition, disruption of reproductive cycles, and increased susceptibility to disease and pathogens. Although the effects of droughts, water withdrawals, and low flow conditions on mussels are an increasing area of study, there is a general paucity of data-driven consensus regarding optimal mitigation approaches and strategies (Cushway et al. 2024). Long-term management decisions regarding water withdrawals, particularly during extended or multi-year drought conditions, may consider factors beyond direct mortality such as the sublethal effects of drought and low flow conditions on freshwater mussels for a comprehensive understanding of the effects of prolonged low flows on freshwater mussels and aquatic communities.

Tennessee↗

USGS2011

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↗

USGS2020

Using advanced population genomics to better understand the relationship between offshore and spawning habitat use for Atlantic Sturgeon

Atlantic Sturgeon ( Acipenser oxyrinchus oxyrinchus ) are a large-bodied anadromous fish that historically supported important fisheries along the east coast of the United States. Following years of overharvest and habitat degradation, populations experienced severe declines. In 2012, the National Marine Fisheries Service listed Atlantic Sturgeon under the Endangered Species Act (ESA; 61 FR 4722). Their listing named five Distinct Population Segments (DPSs), predicated on genetic groups composed of geographically proximate populations. Federal management of Atlantic Sturgeon presents challenges, as sturgeon from each of the five DPSs mix extensively in coastal and marine habitats yet take and recovery progress must be evaluated separately for each unit. Genetic assignment testing based on mitochondrial and microsatellite markers allows individuals to be assigned back to their natal river and DPS. However, this approach is not perfect and some individuals may be incorrectly assigned. Recent advances in genomics offer the potential of a higher resolution approach to genetic assignment testing, and thus may reduce uncertainty associated with assignment testing. In addition, genomics allows a greater number of markers to be examined from across a broader portion of the sturgeon genome, thus may provide an enhanced perspective of population structure for the species, and potentially allow other previously intractable questions to be addressed (Bernatchez et al. 2017, Supple and Shapiro 2018). We used next-generation sequencing to develop a draft genome for Atlantic Sturgeon and identify single nucleotide polymorphisms (SNPs) that could be used to resolve the natal river and DPS of individual Atlantic Sturgeon. We identified 1,210 candidate SNPs within the nuclear genome as well as 49 SNPs within the mitochondrial genome. After filtering and review, we selected 161 nuclear SNPs and 39 mitochondrial SNPs for further testing and evaluation. We used genotyping-in-thousands by sequencing (GT-seq) to simultaneously sequence nuclear SNP loci, mitochondrial SNP loci, and the existing panel of twelve microsatellite loci. This effort required a pilot sequencing run on a single sturgeon sample to test marker amplification and refine primer strengths, followed by a series of sequencing runs to generate baseline data for 288 individuals representing nine populations of Atlantic Sturgeon in four DPSs. Using baseline data from the nine populations, we ran a series of genomic analyses to characterize diversity within and among populations, providing a benchmark for this species using the new SNP markers. Allelic richness was similar for all populations, although there was a general trend of more northern population containing greater levels of allelic richness. Interestingly, we observed linkage disequilibrium among many pairs of loci within many populations. This might be the result of physical linkage but could also suggest these populations are recovering from genetic bottlenecks and/or are effectively small, leading to specific haplotypes to be favored by chance. Pairwise differentiation among populations varied among the populations ( F ST range: 0.010-0.098) and was significantly correlated ( r = 0.771; P < 0.001) to pairwise F ST observed using microsatellite markers). Population clustering and ordination techniques using the new genomic data both support an overall population structure that is similar to the current DPS management units (which were developed primarily based on microsatellite genetic data). Overall, this suggests that existing microsatellite markers and the panel of SNP markers developed in this study provide similar information about the populations structure and ecology of Atlantic Sturgeon. Given the observed differences in allele frequencies among populations, our genomic baseline supports previous assertations that Atlantic Sturgeon show natal homing, despite mixing extensively in marine waters during non-breeding periods. Lower levels of differentiation between populations in the South Atlantic DPS suggest that populations in this region may have greater levels of gene flow relative to their more northerly conspecifics, which has also previously been suggested based on microsatellite data. The observed differentiation among populations provides the necessary foundation for determining the natal river and DPS of Atlantic Sturgeon using assignment testing. We tested the utility of our new genomic baseline for resolving the population and DPS of Atlantic Sturgeon. Our nuclear SNP markers showed utility for identifying the origin of unknown Atlantic Sturgeon samples, as 86.5% were assigned to the correct DPS and 66.3% were assigned to the correct natal river. However, since this study was funded the Conservation Genetics and Genomics Laboratory at Leetown Science Center has made significant improvements to their microsatellite genetic baseline, which now performs more effectively than our new genomic approach (the genetic baseline includes 12 populations and 5 DPSs, and correctly assigns 95.8% of individuals to DPS and 84.9% of individuals to their natal population using 12 microsatellite loci). We conducted an ad hoc exploration of how additional microsatellite or nuclear SNP loci may further improve the accuracy of assignment testing. We found that additional microsatellite markers are likely to result in greater improvements in assignment efficiency than additional nuclear SNPs. However, a much larger number of SNP loci (which if identified could be sequenced using other methods that are now available; e.g., the RAD-capture approach published by Ali et al. 2016) could produce assignment efficiencies that are greater than what is currently feasible using microsatellites. In the absence of further research and development of additional SNP markers for Atlantic Sturgeon (possibly using an approach other than GT-seq), the existing microsatellite loci are the most effective means available to determine the natal river and DPS of Atlantic Sturgeon encountered in offshore waters. Because our new genomic markers were less effective than the existing panel of 12 microsatellite markers, we chose to use the existing microsatellite markers to assign Atlantic Sturgeon captured in another BOEM-funded study (cooperative agreement M16AC00003; Monitoring endangered Atlantic Sturgeon and commercial finfish habitat use offshore New York) following consultation with our project officer. Using this approach, we genotyped and assigned 186 Atlantic Sturgeon captured in coastal waters off the Rockaway Peninsula, New York. The vast majority of these sturgeon were assigned to the New York Bight DPS (94.62%), and most appear to belong to the Hudson River population (87.10%) with smaller contributions from the Delaware River population (7.53%). Smaller contributions (2.15%) were observed from six other populations, including those from the James, York, Kennebec, Ogeechee, and Edisto rivers. Although most of the fish we assigned were assigned to the nearest spawning rivers (Hudson and Delaware), the contributions from distant rivers is consistent with the propensity of this species to move long distances and form mixed stock aggregations along the continental shelf. This finding indicates that spawning populations (and their corresponding DPS) from distant locations may potentially be impacted by offshore activities. In fact, activities in this region of the New York Bight could negatively impact Atlantic Sturgeon population from at least four different DPSs. Genetic or genomic assignment testing remains an essential tool to characterize potential impacts to Atlantic Sturgeon populations and should be applied more broadly to better characterize potential impacts of activities in other locations.

Atlantic Coast↗

USGS2015

In-reservoir behavior, dam passage, and downstream migration of juvenile Chinook salmon and juvenile steelhead from Detroit Reservoir and Dam to Portland, Oregon, February 2013-February 2014

In the second year of 2 years of study, the movements of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) and juvenile summer steelhead ( Oncorhynchus mykiss ) through Detroit Reservoir, passing Detroit Dam, and migrating downstream to Portland, Oregon, were studied during a 1-year-long period beginning in February 2013. The primary purpose of the study was to provide empirical data to inform decisions about future alternatives for improving downstream passage of salmonids at Detroit Dam. A secondary purpose was to design and assess the performance of a system to detect juvenile salmonids implanted with acoustic transmitters migrating in the Willamette River. Inferences about fish migration were made from detections of juvenile fish of hatchery origin at least 95 millimeters in fork length surgically implanted with an acoustic transmitter and released during the spring (March–May) and fall (September–November) of 2013. Detection sites were placed throughout the reservoir, near the dam, and at two sites in the North Santiam River and at three sites in the Willamette River culminating at Portland, Oregon. We based most inferences on an analysis period up to the 90th percentile of tag life (68–78 days after release, depending on species and season), although a small number of fish passed after that period as late as April 8, 2014. Chinook salmon migrated from the tributaries of release to the reservoir in greater proportion than steelhead, particularly in the fall. The in-reservoir migration behaviors and dam passage of the two species were similar during the spring study, but during the fall study, few steelhead reached the reservoir and none passed the dam within the analysis period. Migrations in the reservoir were directed and non-random, except in the forebay. Depths of fish within 25 meters of the dam were deeper in the day than at night for Chinook salmon and similar in the day and night for steelhead; steelhead generally were at shallower depths than Chinook salmon. The primary factors affecting dam passage rates were seasonal dam operating conditions and diel period. Fish passage rates were much greater during the spring and summer than in the fall and winter, and the difference was attributed to the availability and use of the spillway near the top of the dam during the spring and summer. The flood-control purpose of the reservoir prevented spillway use during much of the fall and winter because of the low forebay elevation. Passage rates at night were greater than in the day during spring and summer (4.2 times) and during the fall and winter (14.9 times). Fish length, dam discharge, and forebay elevation also affected dam passage rates. Travel times from Detroit Dam passage to the downstream sites were shorter during the fall and winter than during the spring and summer, and were less than a median of 8.68 days to Portland. The estimated survival in the 11 kilometers (km) between Detroit Dam and the Minto Dam forebay was lower than in the remaining 241 km to the Portland site. Estimated survival per 100 km in the free-flowing reach from Minto Dam to Portland was 0.675–0.836, depending on species and season, and was similar to other free-flowing rivers in the Western United States. The high probability of fish in the reservoir reaching the dam, the chance for repeated presence near the dam, the fish depths, and the factors known to affect passage rates suggest that a properly designed surface passage route could be a viable downstream passage alternative for juvenile Chinook salmon and steelhead at Detroit Dam. As part of the evaluations conducted at Detroit Dam, we continued to refine and improve methods for monitoring fish movements in the Willamette River. The goal was to develop stable, cost-effective, long-term monitoring arrays suitable for detection of any Juvenile Salmon Acoustic Telemetry System (JSATS)-tagged fish in the Willamette River. These data then could be used to estimate timing, migration rates, and survival of JSATS-tagged fish from various studies in the Willamette River Basin. The challenge, however, is that acoustic telemetry generally performs poorly in shallow, turbulent water, like that found in the Willamette River. We successfully designed, deployed, and maintained a series of monitoring sites near the Oregon cities of Salem, Wilsonville, and Portland. In the spring, detection probabilities at these sites ranged from 0.900 to 1.000. In the fall, the detection probabilities decreased and ranged from 0.526 to 1.000. The lower detection probabilities, particularly at the Salem site (0.526), were owing to loss of data caused by abnormally high flows as well as the 2013 Federal government shutdown, which prevented us from servicing the equipment. The monitoring sites that we installed seem to be robust and enable the efficient use of acoustic-tagged fish for studies of migration or survival in the Willamette River and similar environments.

Oregon↗

USGS2020

Recent planform changes in the Upper Mississippi River

Geomorphic changes in the Upper Mississippi River (UMR) have long been a concern of river agencies charged with maintaining and restoring river habitat (GREAT 1980; Jackson et al. 1981; USFWS 1992). Large meandering alluvial rivers like the UMR are expected to constantly change and adjust their fluvial landforms within their riparian corridors as a result of the natural interaction of hydrologic processes, sediment movement, and vegetation over time. However, present geomorphic changes in the UMR reflect altered hydrologic, hydraulic, and sediment conditions caused by regulated flows, constructed agricultural levees and navigation dams, altered land use in the watershed, and climate change. Levees reduce lateral hydrologic and sediment connectivity between channels and floodplains on many tributaries and on the Mississippi River downstream of Pool 13. Between each of the dams are a repeating series of landforms associated with tailwater, intermediate, and impounded conditions. The dams maintain a minimum water level, thus creating many off-channel areas that act as sediment traps. Whereas high-head dams cut off sedimentological connectivity longitudinally through the river corridor (Skalak et al., 2013), low head dams on the UMR only slightly altered transport longitudinally. Deltaic-like sedimentation can be common in the impounded sections of dammed rivers. Erosion of relict land surfaces that remained above the raised impounded water levels has been the dominant change in UMR impounded sections due to increased wind fetch leading to increased wave action. Even though upland sources of sediment from tributaries have decreased over the middle to late 20th century, increased annual precipitation, the interplay of increased variability in flood magnitudes from year to year, and more fall and winter flooding have likely changed erosion and sedimentation patterns in the UMR (Belby, et al., 2019). Paradoxically, monitoring and research indicates that the concentration of some water column constituents like total suspended solids and phosphorous has decreased during the 1991 to 2014 time period (Kreiling and Houser, 2016). In areas prone to increased sedimentation, bed elevations rise and thereby water depths are reduced at a given discharge, resulting in loss of fish habitat. Sediment deposition or erosion further influences water exchange rates between main channel and off-channel areas in the river by increasing resistance in connecting channels or enlarging existing connecting channels. Water depth and water exchange rates are the most prominent features describing habitat quality in the UMR (De Jager et al. 2018), and in some cases, the trajectory of planform change from heightened deposition promises to threaten deep backwater habitats particularly important for overwintering fish. Although information on the rate of vertical change in bed elevation is needed for a complete assessment of geomorphic change associated with the loss of deep backwater habitats, mapping planform changes over time (i.e., lateral changes between the land-water boundary) provide needed information on the location, potential cause, and progressive direction of deposition, especially in the mid sections between dams where deltaic processes are the most pronounced. Several types of planform changes have been observed and identified as concerns. For example, island loss in the large impounded areas of the upper part of the UMR was one of the concerns identified by river managers in the 1980s and 90s, and subsequently island construction became a common form of restoration implemented by the Upper Mississippi River Restoration (UMRR) Program (USACE 2012). Other subtler planform changes, such as channel bank erosion and delta formation in backwaters, are perceived to be important, but have largely gone unquantified. A systemwide reconnaissance of the UMR and IWW conducted in 1998 concluded that 14-percent of the river banks were eroding (Nakato and Anderson 1998). However, stabilization of existing river banks has never been widely pursued as a restoration measure, due to the high cost and uncertain benefits. Delta formation reduces the amount of backwater habitat; however, the deltas maintain and create a mix of riparian and aquatic habitats, and that is generally considered to be beneficial for wildlife and fish. If recent hydrologic trends of more frequent and longer duration flood events continue, a better understanding of planform changes can help in describing past changes, and then be used to forecast potential future trajectories of change. If UMR resource managers determine that past and forecasted conditions are undesirable, then UMRR projects could be identified and prioritized to address those concerns. Vegetative cover associations with landform changes have been used to detect and quantify planform changes in many rivers (Johnson 1985; Hiatt 2015; Volte et al. 2015). Freyer and Jefferson (2013) completed such a study in Pool 6 of the UMR using the landcover data from 12 dates over a 115-yr period, including the 1989, 2000, and 2010/2011 landcover/use (LCU) data from the UMRR Program. Planform change detected over the last 20 years represented by the UMRR Program data best reflect present-day geomorphic patterns, rates and processes. Changes occurring prior to dam construction and changes occurring soon after dam construction are likely not the same as those happening now, 50-70 years after dam construction and creation of the impoundments (McHenry et al., 1984; Bhowmik and Adams, 1986; WEST Consultants, 2000). The LCU data from each of the 1989, 2000, and 2010/2011 imagery was developed using similar methods and is available in a Geographical Information System (GIS) for the entire UMR and therefore provides the opportunity for a more comprehensive planform change analysis. This study used GIS overlays of LCU classes to map and quantify changes in planform features over two periods, looking specifically for depositional areas where terrestrial and wetland vegetation expanded at the expense of open water. The land expansion was grouped into four possible process-based types common in large floodplain rivers, some following that used by Lewin et al. (2017). The four types include: crevasse deltas emanating from a breach from a main channel through a natural levee or narrow floodplain into backwaters (crevasse deltas), tributary deltas expanding into backwaters (tributary deltas), deltaic bars at the upstream end of impoundments (impounded deltas), and linear-like bars extending from the downstream ends of narrow levees and remnant floodplains (bar-tail limbs). The methods deployed for change detection addressed possible errors from a variety of sources.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

USGS2018

California State Waters Map Series — Offshore of Point Conception, California

Introduction In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within the 3-nautical-mile limit of California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow subsurface geology. The Offshore of Point Conception map area is in the westernmost part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and this region is presently undergoing north-south shortening. The offshore part of the map area lies south of the stee p south and west flanks of t he Santa Ynez Mountains . The crest of the range , which has a maximum elevation of about 340 m in the map area, lies about 5 km north and east of the arcuate shoreline. The onland part of the coastal zone is remote and sparsely populated. The road to Jalama Beach County Park provides the only public coastal access in the entire map area. North of this county park, the coastal zone is part of Vandenberg Air Force Base. South of Jalama Beach County Park, most of the coastal zone is part of the Cojo-Jalama Ranch, purchased by the Nature Conservancy in December 2017. A relatively small part of the coastal zone in the eastern part of the map area lies within the privately owned Hollister Ranch. The nearest significant commercial centers are Lompoc (population, about 42,000), about 10 km north of the map area, and Goleta (population, about 30,000), about 50 km east of the map area. The Union Pacific railroad tracks run west and northwest along the coast through the entire map area, within a few hundred meters of the shoreline. The map area has a long history of petroleum exploration, and the seafloor notably includes large asphalt mounds and pockmarks that result from petroleum seepage. Several offshore gas and oil fields were discovered, and some were developed, in and on the margin of California’s State Waters. Much of the shoreline in the Offshore of Point Conception map area is characterized by narrow beaches that have thin sediment cover above bedrock platforms, backed by low (10- to 20-m-high) cliffs that are capped by a coastal terrace. Beaches are subject to wave erosion during winter storms, followed by gradual sediment recovery or accretion in the late spring, summer, and fall months during the gentler wave climate. The map area lies in the west-central part of the Santa Barbara littoral cell, which is characterized by west-to-east transport of sediment from Point Arguello on the northwest to Hueneme and Mugu Canyons on the southeast. Sediment supply to the map area is mainly from relatively small coastal watersheds, including the Jalama Creek–Espada Creek drainage basin (about 63 km 2 ), as well as Cañada del Jolloru, Black Canyon, Wood Canyon, Cañada del Cojo, and Barranca Honda. Coastal-watershed discharge and sediment load are highly variable, characterized by brief large events during major winter storms and long periods of low (or no) flow and minimal sediment load between storms. In recent (recorded) history, the majority of high-discharge, high-sediment-flux events have been associated with El Niño phases of the El Niño–Southern Oscillation climatic pattern. Following the coastline, the shelf bends to the north and northwest around Point Conception, and the trend of the shelf break changes from about 298° to 241° azimuth. Shelf width ranges from about 5 km south of Point Conception to about 11 km northwest of it; the slope ranges from about 1.0° to 1.2° to about 0.7° south and northwest of Point Conception, respectively. Southwest of Point Conception, the shelf break and upper slope are incised by a 600-m-wide, 20- to 30-m-deep, south-facing trough, one of five heads of the informally named Arguello submarine canyon. The map area is located at a major biogeographic transition zone between the east-west-trending Santa Barbara Channel region of the Southern California Bight and the northwest-trending central California coast. North of Point Conception, the coast is subjected to high wave exposure from the north, west, and south, as well as consistently strong upwelling that brings cold, nutrient-rich waters to the surface. Southeast of Point Conception, the Santa Barbara Channel is largely protected from strong north swells by Point Conception and from south swells by the Channel Islands; surface waters are warmer, and upwelling is weak and seasonal. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft, unconsolidated sediment interspersed with isolated areas of rocky habitat that support kelp-forest communities in the nearshore and rocky-reef communities in deeper water. The potential marine benthic habitat types mapped in the Offshore of Point Conception map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats. The fairly homogeneous seafloor of sediment and low-relief bedrock provides characteristic habitat for rockfish, groundfish, crabs, shrimp, and other marine benthic organisms. Several areas of smooth sediment form nearshore terraces that have relatively steep, smooth fronts, which are attractive to groundfish. Below the steep shelf break, soft, unconsolidated sediment is interrupted by the heads of several submarine canyons, gullies, and rills, also good potential habitat for rockfish. The map area includes the large (58.3 km 2 ) Point Conception State Marine Reserve.

California↗

USGS2025

The Grouse & Grazing Project: Effects of cattle grazing on demographic traits of greater sage-grouse

Greater sage-grouse (Centrocercus urophasianus) were once widespread within sagebrush -grassland ecosystems of western North America, but populations have declined since the mid-1960s. Though sage-grouse were not listed as threatened or endangered under the Endangered Species Act (ESA), when examined in 2015, they remain a species of interest and concern. Roughly half of the sage-grouse’s remaining habitat is on federal land, most of it managed by the Bureau of Land Management (BLM) and the U.S. Forest Service (USFS). Livestock grazing is the most extensive land use within sage-grouse habitat and the effects of livestock grazing on sage-grouse are often debated. The extensive decade-long research project summarized in this report was initiated to provide rigorous experimental research to inform the debate regarding the relationship between livestock grazing and sage-grouse. In 2012, the Idaho Grouse & Grazing Project was started with several partners including the University of Idaho, BLM, Idaho Department of Fish and Game (IDFG), and other partners to evaluate the effects of cattle grazing on sage-grouse vital rates. Many additional supporters have provided resources to this research effort including the Public Lands Council, Idaho Cattle Association, Idaho Governor’s Office of Species Conservation, Western Association of Fish & Wildlife Agencies, U.S. Fish and Wildlife Service, USFS, and numerous grazing associations and ranchers in Idaho. This 10-year research project was a scientifically rigorous and replicated experiment, occurring across five study sites in Idaho. This document is intended to provide a summary of the findings of this unprecedented study. Annual reports are available on the project’s website: https://idahogrousegrazing.org and scientific papers are being prepared and submitted to journals. The project focused on the influence of spring cattle grazing on sage-grouse vital rates across five study sites in Idaho including 21 BLM grazing pastures. From 2014-2023, we captured 1,343 grouse, documented the fate of 1,285 nests, and tracked 399 broods. Vegetation was characterized at 4,777 plots and grazing utilization levels were recorded at >30,000 locations. Because insects are an important food source for sage-grouse hens and their chicks, insect biomass and diversity were also examined in this study. We collected arthropods in 12,151 pitfall samples and 6,217 sweep-net samples across 786 plots within our five study sites. At each study site, three or four grazing treatments were implemented after two years of pre-treatment field investigations. These controlled cattle grazing treatments included spring-grazing in even years, spring-grazing in odd years, spring-and-fall grazing in alternating years, and a no grazing (or rested) control. Once grazing treatments were implemented at a study site, we measured sage-grouse demographic traits for 4-8 years post-treatment. Stocking rate (grazing intensity) was assessed across pastures each year and was influenced by vegetation communities, topography, and water sources. Grazed pastures exhibited lower grass cover and height compared to the no grazing pastures, and the extent of this difference varied based on annual precipitation levels. Rested pastures maintained higher grass cover and grass height, but the differences in habitat structure did not consistently translate to differences in sage-grouse demographic traits. Apparent nesting success varied annually and by site, ranging from 24% to 44% over the study period. Like some other studies, results from this research show that successful (i.e., hatched) sage-grouse nests have taller grass heights than failed nests. The average grass height surrounding successful nests in grazed pastures was shorter than that surrounding successful nests in non-grazed (i.e., rested) pastures. It is well documented that grazing reduces grass height, and these observations have led to widely held assumptions that livestock grazing reduces grass height which negatively affects sage-grouse nesting habitat. At the pasture scale, this study has found that sage-grouse nesting success is no greater in pastures that were rested for 4-8 years than those currently or recently grazed. This study gives no indication that removing cattle from pastures affected nesting success. We found some evidence that nest density varied among the grazing treatments, but we did not see compelling evidence of increases in density of nesting hens following cessation of grazing in the no grazing treatments. Brood survival varied by site and year but showed no strong effect of grazing treatment. Climatic conditions, particularly drought in 2021, had a greater effect on brood survival than grazing metrics. We also found no differences in hen survival among the grazing treatments. Results of this study suggest that hens nesting in spring and fall grazed pastures had similar or even slightly higher brood survival than hens in the rested pastures or the spring grazed pastures. Arthropod biomass and species diversity varied among our study sites and the differences between grazed and rested pastures also varied among study sites. Average biomass and diversity of arthropods was higher in the spring grazed pastures on two of three sites examined but higher in the rested pastures on the other site examined. Some taxa of arthropods were more abundant in grazed pastures while other taxa were more abundant in rested pastures. For example, Carabidae (Ground Beetles) and Formicidae (Ants) had higher biomass in grazed pastures, while Tenebrionidae (Darkling Beetles) and Acrididae (Grasshoppers) had higher biomass in non-grazed pastures. Results indicate that grazing effects on arthropod biomass and arthropod diversity are study site-dependent, suggesting a need to better quantify the most important prey taxa for sage-grouse chicks and to better control for other factors that influence arthropod abundance. Based on results of this research, livestock grazing, when properly managed, does not appear to negatively impact sage-grouse nest survival or brood success. This study provides critical insights for land managers balancing livestock production with sage-grouse conservation, supporting adaptive grazing strategies that maintain both economic and ecological objectives.

Idaho↗

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

USGS2021

Geologic and hydrogeologic characteristics of the White River Formation, Lance Formation, and Fox Hills Sandstone, northern greater Denver Basin, southeastern Laramie County, Wyoming

In cooperation with the Wyoming State Engineer’s Office, the U.S. Geological Survey studied the geologic and hydrogeologic characteristics of Cenozoic and Upper Cretaceous strata at a location in southeastern Laramie County within the Wyoming part of the Cheyenne Basin, the northern subbasin of the greater Denver Basin. The study aimed to improve understanding of the aquifers/aquifer systems in these strata, motivated in part by declining groundwater levels and interest in exploring future groundwater supplies. Based on detailed geologic characterization using information obtained by drilling and coring a 960-foot-(ft) deep exploratory borehole, and comparisons with previously published descriptions, identified Cenozoic lithostratigraphic units included 40 ft of Quaternary older alluvial fan deposits consisting of an unconsolidated mixture of sand and gravel with lesser quantities of silt and clay in varying proportions and the underlying 407.3-ft-thick White River Formation of late Eocene-Oligocene age consisting largely of mudrocks with sparse thin beds of sandstone, muddy gravel, and conglomeratic mudrocks. Identified Upper Cretaceous lithostratigraphic units included the 351.6-ft-thick Lance Formation, consisting of terrestrial sedimentary rocks including mudrocks (muddy shale and silty and sandy shale, siltstone, claystone, and mudstone) interbedded with much smaller quantities of very fine- to medium-grained muddy and silty sandstone and coal; the 79.6-ft-thick Fox Hills Sandstone, consisting of a transitional marine sequence of muddy or silty sandstone present in five individual beds; and 86.7 ft of the upper transition member of the Pierre Shale, consisting largely of marine sedimentary rocks such as muddy shale. Beds of the upper and lower Fox Hills Sandstone were separated by tongues of the Lance Formation and upper transition member of the Pierre Shale, respectively. The White River hydrogeologic unit, consisting of the entire White River Formation or Group at the study site, did not contain any substantial secondary permeability features in the mudrocks that composed almost all the unit. A monitoring well (BR–1) was completed in the White River aquifer with the well screen open to the only coarse-grained unit (muddy sandstone) that had sufficient thickness and permeability to be considered as an aquifer. Sampling of the well for a broad suite of constituents indicated groundwater generally was of excellent quality except dissolved arsenic was detected at a concentration greater than the U.S. Environmental Protection Agency (EPA) Maximum Contaminant Level, and dissolved sodium was measured at a concentration greater than several EPA Drinking Water Advisory Levels (DWAs) for the constituent. Well development, well purging for groundwater sampling, and calculated aquifer properties indicated the sandstone aquifer screened by monitoring well BR–1 was not very productive. Analysis of the well water-level responses in BR–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response with wellbore-storage effects. Hydraulic properties estimated based on these responses yielded values of hydraulic conductivity ( K , 0.057 foot per day [ft/d]), specific storage ( Ss , 1.6×10 −6 per foot [ft −1 ]) and porosity ( n , 0.43). Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated an upward trend (+1.13 foot per year [ft/yr]) during the period analyzed, September 5, 2014, to September 30, 2017. Lithologic characteristics of the Lance hydrogeologic unit, consisting of the entire Lance Formation at the study site, indicated a potential aquifer in a “sandy” interval in the upper part of the unit. Most of the Lance hydrogeologic unit below the “sandy” interval consisted of various low-permeability lithologies unlikely to yield substantial quantities of water. This lower part of the hydrogeologic unit likely functions as a confining unit separating the underlying Lance-Fox Hills aquifer. A geologic cross section constructed for this study indicated fine-grained sediments composed most of the Lance Formation/hydrogeologic unit not only at the study location, but also throughout southern Laramie County along the line of section and throughout the Wyoming and Colorado parts of the Cheyenne Basin. A monitoring well (LN–1) completed in a sandstone bed in the “sandy” interval of the Lance hydrogeologic unit produced a mean of about 23 gallons per minute (gal/min) during well development, indicating sandstone beds can form moderately productive confined subaquifers in this part of the hydrogeologic unit. Analysis of the well water-level responses in well LN–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response. Hydraulic properties estimated based on these responses yielded values for a lower bounding K of 0.60 ft/d, Ss of 1.6×10 −6 ft −1 , and n of 0.38. Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated a downward trend (−0.86 ft/yr) during the period analyzed (November 8, 2014, to September 30, 2017). Analyses for a broad suite of constituents in samples from well LN–1 indicated groundwater quality generally was excellent, although dissolved sodium was measured at a concentration greater than two EPA DWA levels for the constituent. Because of the absence of any overlying or intertonguing sandstone beds belonging to the lower/basal part of the Lance Formation, the Lance-Fox Hills aquifer at the study site consisted only of the five sandstone beds of the Fox Hills Sandstone. The cross section constructed for this study illustrated how the Fox Hills Sandstone, and thus, most of the Lance-Fox Hills aquifer, consists of a series of sandstone bodies that overlap (shingle) upward to the east across southern Laramie County. These bodies collectively form a fairly continuous body of sandstone, thus potentially forming an areally extensive aquifer across southern Laramie County, and by extension, throughout most of the formation’s extent in the Wyoming part of the Cheyenne Basin, as is the case in the Colorado part of the basin. A monitoring well (FH–1) completed in part of the thickest sandstone bed of the Lance-Fox Hills aquifer was moderately to highly productive and easily produced 25 to 30 gal/min after development. Substantially larger water production rates likely could be obtained by penetrating the full thickness of this bed and by completing a well open to the other overlying and underlying sandstone beds of the aquifer. Analysis of the water-level responses in well FH–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response. Hydraulic properties computed based on these responses yielded values for a lower bounding estimate for K of 0.26 ft/d, for Ss of 1.0×10 −6 ft −1 , and for n of 0.41. Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated a downward trend (−1.74 ft/yr) during the period analyzed, December 19, 2014, to September 30, 2017. Sampling of monitoring well FH–1 and two production wells completed in the Fox Hills Sandstone in other parts of Laramie County indicated groundwater quality generally is excellent, although pH exceeded a recommended EPA aesthetic drinking-water standard (Secondary Maximum Contaminant Level) in two of three sampled wells, total dissolved solids concentrations exceeded the Secondary Maximum Contaminant Level in one of the two sampled production wells, and dissolved sodium was measured in all three sampled wells at a concentration greater than two EPA DWA levels for the constituent. The Wyoming Class II agricultural (irrigation) sodium adsorption ratio standard of 8 was exceeded in all three sampled wells, indicating these waters are not suitable for irrigation use. Computed vertical hydraulic gradients indicated a strong potential for downward flow throughout the groundwater system at the study site, including from the low-yielding aquifer in the upper White River Formation/hydrogeologic unit (monitoring well BR–1) to the sandstone subaquifer in the Lance Formation/hydrogeologic unit (monitoring well LN–1), and from the Lance subaquifer (monitoring well LN–1) to the sandstone bed/aquifer that composes much of the Lance-Fox Hills aquifer thickness at the study site (monitoring well FH–1). However, large hydraulic-head differences between wells indicated high resistance to vertical flow attributable to the low vertical hydraulic conductivity of intervening strata, which consisted almost entirely of low-permeability mudrocks. The confined nature of the sandstone aquifers monitored by the various wells coupled with dissimilarities between groundwater-level fluctuations and trends in groundwater levels indicated downward flow through the intervening strata (primarily mudrocks in the various lithostratigraphic/hydrogeologic units) between the examined sets of wells likely was small.

Wyoming↗

USGS2015

Southern Salish Sea Habitat Map Series: Admiralty Inlet

In 2010 the Environmental Protection Agency, Region 10 initiated the Puget Sound Scientific Studies and Technical Investigations Assistance Program, designed to support research in support of implementing the Puget Sound Action Agenda. The Action Agenda was created in response to Puget Sound having been designated as one of 28 estuaries of national significance under section 320 of the U.S. Clean Water Act, and its overall goal is to restore the Puget Sound Estuary's environment by 2020. The Southern Salish Sea Mapping Project was funded by the Assistance Program request for proposals process, which also supports a large number of coastal-zone- and ocean-management issues. The issues include the recommendations of the Marine Protected Areas Work Group to the Washington State Legislature (Van Cleve and others, 2009), which endorses a Puget Sound and coast-wide marine conservation needs assessment, gap analysis of existing Marine Protected Areas (MPA) and recommendations for action. This publication is the first of four U.S. Geological Survey Scientific Investigation Maps that make up the Southern Salish Sea Mapping Project. The remaining three map blocks to be published in the future, located south of Admiralty Inlet, are shown in figure 1. Puget Sound is a deep, fjord-type estuary covering an area of 2,330 km 2 in the Pacific Northwest region of the United States (fig. 1). It is connected to the ocean by the Strait of Juan de Fuca, a turbulent passage approximately 160 km in length and 22 km wide at its west end, expanding to over 40 km wide at its east end (Thomson, 1994). During the Pleistocene, the area was occupied several times by lobes of continental ice, resulting in a complex basin-fill of glacial and interglacial deposits that are locally as thick as 1100 m (Johnson and others, 2001). The last glaciation, called the Fraser glaciation, began after 28,800±740 14 C yr B.P. when ice started a slow expansion (Clague, 1981). At peak advance the westward Juan de Fuca lobe reached the edge of the continental shelf through the Juan de Fuca Strait shortly before 14,460±200 14 C yr B.P. (Herzer and Bornhold, 1982). The southward Puget lobe advanced to its terminal position in Puget Sound by around 14,150 14 C yr B.P. (Porter and Swanson, 1998). Ice retreated from its maximum to northern Whidbey Island by 13,650±350 14 C yr B.P. (Dethier and others, 1995). Retreating glaciers resulted in a thick sequence of ice-contact, glacial-marine sediment, and early post-glacial sediments (Linden and Schurrer, 1988). These deposits have experienced the effects of a marine transgression followed by regression, resulting in a sea-level several tens of meters lower than the present day (Linden and Schurrer, 1988). A second transgression brought sea level to about the present level by around 5,470±120 14 C yr B.P. (Clague and others, 1982) establishing the present oceanographic and geologic environment Puget Sound is separated into four interconnected basins; Whidbey, Central (Main), Hood Canal, and South (Thomson, 1994). The Whidbey, Central, and Hood Canal basins are the three main branches of the Puget Sound estuary and are separated from the Strait of Juan de Fuca by a double sill at Admiralty Inlet. The Admiralty Inlet map area includes the Inlet and a portion of the Whidbey Basin (fig. 1). The shallower South Basin is separated by a sill at Tacoma Narrows and is highly branched with numerous finger inlets. Flow within Puget Sound is dominated by tidal currents of as much as 1 m/s at Admiralty Inlet, reducing to approximately 0.5 m/s in the Central Basin (Lavelle and others, 1988). The lack of silt and clay-sized sediments in the Admiralty Inlet map area is likely a result of the strong currents (see Ground-Truth Studies for the Admiralty Inlet Map Area, sheet 3). The subtidal component of flow reaches approximately 0.1 m/s and is driven by density gradients arising from the contrast in salty ocean water at the entrance and freshwater inputs from stream flow (Lavelle and others, 1988). The total freshwater input to Puget Sound is approximately 3.4 x 10 6 m 3 /day, primarily from the Skagit River (Cannon, 1983). The subtidal circulation mostly consists of a two-layered flow in the basins with fresher water exiting at the surface and saltier water entering at depth (Ebbesmeyer and Cannon, 2001). In general, surface waters flow north and deeper waters flow south; variations arise from wind effects that can drive a surface current in the same direction as the wind, and a baroclinic response in the lower layer to about 100-m depth (Matsuura and Cannon, 1997). Oceanographic properties are influenced by temporal forcing parameters such as reduced stream flow during the 2000-01 drought that increased surface salinity and decreased differences between surface and bottom waters (Newton and others, 2003). On offshore seismic-reflection profiles, Pleistocene strata (excluding latest Pleistocene glacial and post-glacial deposits) form a distinct seismic unit, bounded below by pre-Tertiary or Tertiary basement and above by typically flat-lying latest Pleistocene to Holocene deposits that fill in erosional or depositional relief (Johnson and others, 2001). Cores from central Puget Sound have accumulation rates that range from 85 to 1200 mg/cm 2 /yr, or 0.12 to 2.4 cm/yr; the highest accumulation rates are near the southern end of central Puget Sound (Carpenter and others, 1985). Carpenter and others (1985) un-weighted arithmetic mean of accumulation rates for central Puget Sound deeper stations is 480±340 (± one standard deviation) mg/cm 2 /yr. Lavelle and others (1985) also found rates as high as 1200 mg/cm 2 /yr over the past approximately 70 years in cores in the Central Basin off of and north and south of Elliott Bay. Puget Sound basin rates are comparable to rates in midshelf silt deposits on the Washington coast north of the Columbia River (Nittrouer and others, 1979). The deep subtidal (in other words, below SCUBA depths) habitats of Puget Sound are relatively poorly known. A few subtidal surveys exist for several habitat types from the 1960s and 1970s (reviewed in Dethier, 1990), using grab and box core data. The Dethier (1990) review divides habitat up into Coast and Marine Ecological Classification Standard (CMECS) substrate, water column energy, and depth zones but does not attempt to map these habitats, rather it is an inventory of habitats found in the area and the flora and fauna associated with each habitat. The approach of the Southern Salish Sea Mapping project is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, and bottom-sediment sampling data. This approach is based in part on methods presented and data collection and product needs identified at the Washington State Seafloor Mapping Workshop (Washington State Seafloor Mapping Workshop Steering Committee, 2008), attended by coastal and marine managers and scientists. The map products display seafloor geomorphology and substrate, and identify potential marine benthic habitats. It is emphasized that the more interpretive habitat and geology maps rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. Oceanographic current and wave data is not included in this analysis, however, the accompanying geographic information system (GIS) data set is designed and intended to be combined with oceanographic and biologic data sets assembled by others in the future and some of the GIS data has already been incorporated in the unpublished Nature Conservancy Benthic Habitats of Puget Sound database. This publication includes four map sheets, explanatory text, and a descriptive pamphlet. Each map sheet is published as a portable document format (PDF) file. ESRI ArcGIS compatible geotiffs (for example, bathymetry) and shapefiles (for example video observation points) will be available for download in the data catalog associated with this publication (Cochrane, 2015). An ArcGIS Project File with the symbology used to generate the map sheets is also provided. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html .

Washington↗

USGS1989

Geology of the Arabian Peninsula; shield area of western Saudi Arabia

Western Arabia lies within the low-latitude desert of north Africa and the Middle East, the core being the Arabian segment of the African Shield. The core of complex basement rocks accounts for about 670,000 km 2 , or one-third of the Arabian Peninsula. Reconnaissance mapping of these crystalline rocks, together with bordering sedimentary rocks and volcanic flows, begun in 1950, resulted during the next 13 years in a series of geologic and geographic maps without extensive texts. The maps served as general guides for development of natural resources, including water supplies, ore deposits, and building materials. An intensive exploration program that began in 1963 and involved numerous geologists has vastly increased geologic information. Rainfall in Arabia is meager and episodic, and vegetation is sparse except in isolated copses on the crest of the Hejaz Range. Comparison of flora with similar species in the Sudan, where records of rainfall have long been kept, allows evaluation of mean annual precipitation. Wandering bedouin following fodder created a delicate balance between population and water supply-now disturbed by wells drilled in alluvium and lava fields. A trapezoidal region of Precambrian crystalline rocks lies along the northeast flank of the Red Sea, with two long prongs extending northwest and southeast for a total of 1,800 km. These basement rocks of the Arabian Shield are well exposed on the uplands, scarp mountains, and coastal pediments where the Phanerozoic cover rocks have been stripped as a result of Paleozoic epeirogeny and Tertiary ramping. The shield outcrops are divided into three tectonic provinces by N. 45°W.- trending shear zones of the Najd fault system of latest Proterozoic and possibly earliest Paleozoic time. The southwestern province, the 'Asir1 upland, was sharply uplifted and tilted to the northeast during the Neogene. The northwestern province, consisting of the Ash Shifa'- Hisma upland as well as Jabal Shammar farther east, similarly was uplifted and tilted. These two provinces are separated by the flat-lying median N ajd province, which is chiefly bounded by the principal Najd faults. The outcrops of the shield rocks are of the Late Proterozoic Eonupper Riphean to Vendian or Infracambrian epochs, including the Ediacarian System. The most reliable isotopic ages range from about 900 to 560 m.y., but some Middle Proterozoic rocks may be present in the easternmost shield. The rocks are divided into six lithostratigraphic sequences, two plutonic suites, and an ophiolitic suite. The mafic and ultramafic volcanic and plutonic rocks of the ophiolitic suite everywhere were emplaced tectonically and are probably of different ages in different places. Some ophiolite occurs as obducted blocks, but most is highly deformed and altered to serpentinite in fault zones that mostly define sutures between different tectonic blocks or terranes within the shield. Three of the lithostratigraphic sequences consist of mafic to silicic volcanic rocks and volcanic-derived clastic rocks which, with their subvolcanic plutonic rocks of a dioritic suite, probably formed in oceanic island arcs during convergent plate tectonism. These rocks make up the primary, or first-formed, crust of the shield. Chemical analyses show that the primary shield rocks, regardless of age, are principally calc-alkalic with some associated tholeiitic varieties. Most of the layered rocks are andesitic, but they range from basalt to dacite and in places contain intercalated pillow basalt, marble, chert, and carbonaceous or graphitic schist. Most of the plutonic rocks of the dioritic suite are dioritic, but they range from gabbro to trondhjemite and rarely contain potassium feldspar. The sequences and an associated dioritic suite become younger toward the eastern shield, that is, the primary crust of the shield youngs toward the east. Two western sequences consist of the Jiddah (Samran) and BaishBahah Groups and range in radiometric age from about 900 to 800 m.y.; the eastern sequence consists of the Halaban (Hulayfah) Group and ranges from 800 to about 700 m.y. During subsequent orogeny, most of the rocks were intensely deformed and mostly metamorphosed to upper greenschist facies, but rising in places to the almandine-amphibolite facies. Two other lithostratigraphic sequences with an associated plutonic granitic suite are the products of two mountain-building episodes during which the primary crust was greatly thickened and converted into craton. The two sequences, including largely the Ablah (Al Ays) and Murdama (Shammar) Groups, consist of abundant sedimentary rocks, commonly arkosic, that are the erosional products of the orogenic mountains. They are several thousand meters thick. Less abundant calc-alkalic to alkalic volcanic rocks, commonly dacitic and rhyolitic, are intercalated with the sedimentary rocks. The plutonic rocks of the granitic suite in association with both sequences have syntectonic and posttectonic phases, are products of the orogenies, and are the principal new ingredients making up the craton. Gneiss domes were a significant part of these cratonization orogenies. In association with orogenic crustal heating, some of the low-density, more silicic tonalitic and trondhjemitic rocks of the primary crust rose as gneiss domes. Partial melting in the middle or lower crust below the gneiss domes produced large volumes of granitic magma that intruded the gneiss domes as granodioritic batholiths. The Ablah Group and the older part of the granitic suite are about 775 to 740 Ma old and are associated with the Ablah orogeny and early cratonization in the western and earlier formed half of the shield. The Murdama (Shammar) Group and the younger part of the granitic suite are about 660 to 580 Ma old and are associated with the culminant orogeny and late cratonization that was shieldwide. The granitic suite during both orogenies consists of early, syntectonic granodiorite batholiths associated with the gneiss domes and late, posttectonic monzogranite plutons. Only during the culminant orogeny, late magmatic evolution produced syenogranite and alkali-feldspar granite commonly in circular and ring-structured plutons and with associated explosive volcanic deposits (Shammar Group); final products, some of which have economic potential, were peralkalic and peraluminous. The late plutonism of the culminant orogeny was distinctly bimodal in that subordinate gabbroic rocks are associated with the granites. Various building blocks or terranes of the andesitic and dioritic primary crust were collisionally agglomerated during the Ablah orogeny, early cratonization, whereas the entire shield as currently exposed was further collisionally accreted and compressionally consolidated during the culminant orogeny, final cratonization. Thousands of kilometers of oceanic crust had to be subducted in about 300 m.y. to form the large primary crust of the Arabian Shield. The inevitable collisional events during consumption of such a large volume of oceanic crust invariably led to numerous collisional orogenies that collectively encompass the widely known Pan African tectonic episode. The youngest lithostratigraphic sequence, the Jubaylah Group, is essentially postcratonic, although it is the end product of the collisional culminant orogeny. Final east-west compression of the entire shield from about 580 to 560 m.y. caused the craton to fracture along the large northwest-trending, left-lateral faults and elsewhere along lesser, northeast-trending, right-lateral, conjugate faults of the N ajd fault system. Erosional products of this more localized deformation were the sedimentary rocks of the Jubaylah Group, which also includes intercalated andesitic to basaltic volcanic rocks of a mafic alkalic compositional trend. The collisional edge of an old continental plate (or tectonic fragments thereof), suspected on the eastern edge of the Arabian Shield, has not been shown with certainty to be exposed. Presumably, widespread contamination from such an old continental crust affects U/Pb, Sm/Nd, Rb/Sr, and common lead ratios in the young plutonic rocks of the easternmost shield. One mass of anorthosite near Jabal Khida' on the central eastern edge of the shield may be a fragment of this old continental plate in that associated granodiorite may be as old as 1,600 to 1,800 Ma. Epeirogenic uplift, erosion, and cooling of the uppermost shield during Early and Middle Cambrian time is indicated by an average fission track age of 510±52 m.y., on sphene from diorite (hornblende K-Ar age of 615±12 m.y.) in the southwestern part of the shield. The hiatus was followed by extensive deposition of the Cambro-Ordovician Saq Sandstone in the north and northeast and the Wajid Sandstone in the southeast and south of the shield. The Cambrian Siq Sandstone had already been deposited in the northern part. During the middle and late Paleozoic, broad epeirogeny caused further erosion of the shield until marine transgression deposited the Upper Permian Khuff Formation at least in the eastern part of the shield. In the southwestern shield, the nonmarine Upper Triassic Khums Sandstone was deposited variably on Wajid or Precambrian rocks and is overlain by limestone of the middle Upper Jurassic Amran Formation. Except for shallow marine sandstone of problematic Cretaceous age deposited on the Amran Formation in the south.western shield and on Precambrian rocks in the northwestern shield, the younger beds on the shield are Paleocene and younger, with the possibility that the lowermost are upper Maestrichtian. The early Tertiary beds contain vertebrate fossils of coastal marine or estuarine environment 250 km east of the Red Sea in the central shield. Marginal marine sediments were deposited in a western tongue of the latest Tethys Sea as late as Eocene on the western shield and at least as far south as Jiddah. The great harrats of flood basalt erupted on th~ western shield during late Oligocene and early Miocene at the same time a 2,000-kmlong continental rift valley developed along the future Red Sea axis. Within this rift valley, Baid freshwater tuffaceous lakebeds were deposited between mafic and silicic volcanoes. During late early Miocene time, the Red Sea opened at a rate of 4.4 cm/yr in a firststage movement while continental dikes and swarms of oceanic tholeiitic dikes, gabbro, and granophyre plutonic rocks were intruded into the rift sedimentary and volcanic rocks at the newly formed continental margin. The continental margin was deformed and greatly extended at this time. About 14 or 15 m.y., as the first-stage spreading stopped, the Red Sea Escarpment rose; its erosion caused deposition of coarse conglomerate of the Bathan Formation. About 3,000 m of evaporite was deposited on the young Red Sea oceanic crust during the late Miocene desiccation crisis. A second stage of sea-floor spreading about 4-5 m.y. produced the Red Sea axial trough, consisting of oceanic crust, as well as renewed uplift and tilting of the three tectonic provinces in response to compression from counterclockwise rotation against the Dead Sea Rift. This late movement caused widespread major stream capture, especially along the wadis that formerly drained southwesterly or northwesterly, the channels turning westward through narrow gorges to the coastal plain and the Red Sea.

Professional Paper↗