USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal Research”

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

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

1,689 records · Page 42Linked to original sources

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia

Chemical evolution and estimated flow velocity of water in the Trinity Aquifer, south-central Texas

Three permeable zones with varying lithology and water chemistry compose the Trinity aquifer, a principal source of water in the 5,500- square-mile study area in south-central Texas. The upper permeable zone locally yields small quantities of water to wells and was not included in this study. The middle permeable zone primarily is composed of limestone with minor amounts of dolostone. Terrigenous sand and marine limestone, with minor amounts of dolostone, are the principal lithologic units in the lower permeable zone. Dissolved solids concentrations range from 329 to 1,820 milligrams per liter in water samples from the middle permeable zone and from 518 to 3,030 milligrams per liter in water samples from the lower permeable zone. Principal hydrochemical facies in the middle permeable zone are calcium magnesium bicarbonate and calcium magnesium sulfate. Hydrochemical facies in ground-water samples from the lower permeable zone vary. Tritium concentrations as large as 5.3 tritium units in the southeastern part of the study area are indicative of relatively recent recharge. Results of a geochemical mass balance simulation along a flowpath in the middle permeable zone indicate a mass transfer of 4.25 millimoles per liter of dolomite dissolved, 5.74 millimoles per liter of gypsum dissolved, 0.46 millimole per liter of sodium chloride dissolved, 8.07 millimoles per liter of calcite precipitated, and 0.67 millimole per liter of calcium-for-sodium cation exchange between solid and aqueous phases. These results support dedolomitization as a principal chemical process in the middle permeable zone of the Trinity aquifer. Results of a simulation along a flowpath in the lower permeable zone indicate a mass transfer of 0.41 millimole per liter of dolomite dissolved, 0.001 millimole per liter of gypsum dissolved, 9.58 millimoles per liter of sodium chloride dissolved, 1.09 millimoles per liter of calcite precipitated, and 1.11 millimoles per liter of sodium-forcalcium cation exchange between solid and aqueous phases. Lower permeable zone processes indicate sodium chloride dissolution, dedolomitization, and cation exchange. Ground-water-flow velocities determined from adjusted carbon-14 ages, calculated using NETPATH, for selected flowpaths in the middle and lower permeable zones were about 1.7 feet per year and less than about 4.4 feet per year, respectively.

Texas

Comparative assessment of a restored and natural wetland using 13C-DNA SIP reveals a higher potential for methane production in the restored wetland

Wetlands are the largest natural source of methane (CH 4 ), a potent greenhouse gas produced by methanogens. Methanogenesis rates are controlled by environmental factors such as redox potential, temperature, and carbon and electron acceptor availability and are presumably dependent on the composition of the active methanogen community. We collected intact soil cores from a restored and natural freshwater depressional wetland on Maryland’s Delmarva Peninsula (USA) to assess the effects of wetland restoration and redox shifts on microbial processes. Intact soil cores were incubated under either saturated (anoxic) or unsaturated (oxic) conditions and amended with 13 C-acetate for quantitative stable isotope probing (qSIP) of the 16S rRNA gene. Restored wetland cores supported a distinct community of methanogens compared to natural cores, and acetoclastic methanogens putatively identified in the genus Methanosarcina were among the most abundant taxa in restored anoxic and oxic cores. The active microbial communities in the restored wetland cores were also distinguished by the unique presence of facultatively anaerobic bacteria belonging to the orders Firmicutes and Bacteroidetes . In natural wetland incubations, methanogen populations were not among the most abundant taxa, and these communities were instead distinguished by the unique presence of aerobic bacteria in the phyla Acidobacteria , Actinobacteria , and class Alphaproteobacteria . Iron-reducing bacteria, in the genus Geobacter , were active across all redox conditions in both the restored and the natural cores, except the natural oxic–anoxic condition. These findings suggest an overall higher potential for methanogenesis in the restored wetland site compared to the natural wetland site, even when there is evidence of Fe reduction.

Maryland

Spatiotemporal variability of algal biomass and nitrate in Owasco and Seneca Lakes in the Finger Lakes Region, New York, in 2019

Cyanobacterial harmful algal blooms (CyanoHABs) have become increasingly common, threatening the security of water resources globally. The U.S. Geological Survey conducted high-resolution nearshore mapping surveys using boat-mounted multiparameter sondes and nitrate sensors during the summer and fall of 2019 on Owasco Lake and Seneca Lake, two lakes with documented CyanoHABs in the Finger Lakes region of New York State. Discrete sensor measurements and water-quality samples were collected at fixed points along survey routes and continuous data were generated at open-water monitoring platforms. This investigation examined whether water-quality information from nearshore surveys was representative of open-water conditions and if nearshore surveys could be used to identify areas with localized nearshore CyanoHABs and potential sources of nutrients not captured by tributary sampling. In addition to comparisons across methods, nearshore concentrations of nitrate and chlorophyll were evaluated relative to tributary outlets, cyanobacterial abundance and biovolume at discrete locations, and the locations of near-surface CyanoHABs that were designated as “confirmed with high toxins” by the New York State Department of Environmental Conservation. Nitrate and chlorophyll concentrations were comparable across methods for each lake, although concentration ranges were typically higher for nearshore mapping datasets than for nearshore discrete datasets. Nearshore surveys indicated areas of nitrate enrichment that varied temporally in both lakes. Orthophosphate was not routinely detected. Across methods, median chlorophyll concentrations were higher for the summer survey than for the fall survey in Owasco Lake. Nearshore chlorophyll concentrations varied more temporally in Owasco Lake than in Seneca Lake. Phytoplankton and cyanobacterial abundance and biovolume were about five times higher in Owasco Lake than in Seneca Lake. Cyanobacteria dominated the phytoplankton community in most samples, and Microcystis comprised the bulk of the cyanobacterial biovolume. The most abundant potential cyanotoxin-producing (specifically microcystins) genera were Microcystis , Synechococcus , Aphanocapsa , and Pseudanabaena . The cyanobacterial community composition was comparable between open-water monitoring platforms and nearshore samples. Microcystins were detected in seven survey samples from Owasco Lake, in one survey sample from Seneca Lake, and in one sample each from the open-water monitoring platforms on Owasco and Seneca Lakes that were collected about 7 days before the fall surveys. Microcystin detections were not consistently associated with high cyanobacterial cell counts or cyanotoxin-producing genera. Results from nearshore surveys were comparable to open-water monitoring platforms and discrete nearshore observations in the absence of nearshore or open-water CyanoHABs in these systems during the study. Patterns of nearshore concentrations of nitrate and chlorophyll from nearshore surveys may aid in the identification of areas with localized nitrate loading and shifts in phytoplankton abundance and community composition.

New York

Small earthquake moment magnitude and implications for frequency–magnitude scaling of injection induced earthquakes of the Raton Basin

Accurate estimation of earthquake source parameters—such as moment magnitudes, corner frequencies, and stress drops—is essential for improving seismic hazard assessments and understanding earthquake physics. In this study, moment magnitudes ( M W ) are calculated for 31,581 earthquakes associated with wastewater injection in the Raton Basin (located along the border between northern New Mexico and southern Colorado) between 2016 and 2024 using radiative transfer theory to fit coda decay envelopes. Our results show that it is feasible to estimate moment magnitudes down to M W ~1 with coda envelopes from a small local monitoring network. Significant differences were found between M W and local magnitudes ( M L ) for small earthquakes ( M < 3.0). A linear relationship was optimized to convert M L to M W : M W = 0.7 M L + 0.96 and M W = 0.73 M L + 0.99 (for the events reported by the U.S. Geological Survey), which can be applied in future studies of Raton Basin seismicity. We find that b -values calculated employing different methods and using M L are approximately 1.0, while those using M W range from 1.2 to 1.4. A larger estimate of the b -value could influence interpretations of the statistical behavior of earthquakes associated with injection and consequently seismic hazard assessments based on a magnitude–frequency distribution. The potential differences between local versus moment magnitude-based earthquake statistics should be considered in other seismically active regions.

Colorado, New Mexico

Assessment of western Oregon debris-flow hazards in burned and unburned environments

In the steep and mountainous environment of western Oregon, debris flows pose a considerable threat to property, infrastructure and life. Wildfire is commonly known to increase the susceptibility of steep slopes to debris flows, but the extent of this process in the western Cascades is not well understood. The US Geological Survey (USGS) currently estimates postfire debris-flow likelihood and triggering rainfall thresholds using a model calibrated to a southern California inventory of debris flows generated by excess runoff within the first year after fire. Because of a lack of available data, this model has not been tested in western Oregon, or in locations where postfire debris flows initiate via other mechanisms (e.g., shallow landslides or in-channel failures). Using repeat field observations and aerial imagery, we developed two new debris-flow inventories within and adjacent to the perimeters of five 2020 wildfires in western Oregon: Archie Creek, Holiday Farm, Beachie Creek, Lionshead and Riverside. The first inventory focuses on postfire debris flows (2020–2022); the second focuses on debris flows prior to fires (1995–2020). Our inventories of prefire and postfire debris flows were used to document initiation mechanisms in Oregon's western Cascades and to evaluate the effects of wildfire. We found that wildfire changed the distribution of debris-flow initiation mechanisms in the western Cascades. After the wildfires, annual rates of runoff-generated debris flows increased by 22% and the number of shallow landslide-initiated debris flows decreased by 17% relative to before the wildfires. Despite this shift, shallow landsliding was the dominant debris-flow initiation mechanism in both unburned and burned environments. We found the performance of the current USGS debris-flow likelihood model was degraded relative to other previously tested locations across the intermountain western United States. Our results highlight the need for improved postfire hazard assessment in western Oregon based on regional model calibration that is tuned to the dominant debris-flow initiation mechanisms.

Oregon

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California

In-situ geochronology using LA-ICP-MS/MS: Application of the Lu-Hf system in carbonate, apatite and fluorite

Geochronology is a fundamental tool for reconstructing earth history and constraining the timing of ore deposit formation. Traditionally, U-rich minerals like zircon, titanite, monazite and apatite have been analysed to determine the timing of mineralization and alteration events. However, not all deposits or alteration assemblages contain these minerals, and the U-Pb system can be disturbed due to open system behaviour. We present a methodology for obtaining in-situ Lu-Hf dates from various minerals to demonstrate the potential for dating ore deposits that lack U-rich minerals. Carbonate, apatite and fluorite minerals from different tectonic settings and geological environments were analysed. We acquired Lu-Hf and REE (rare earth element) data from carbonate minerals from a carbonatite orebody and dike samples from the Mountain Pass REE carbonatite deposit in California, USA. We also tested a ‘double-dating’ method, where U-Th-Pb dates and Lu-Hf dates were obtained from the same analysis of apatite from rocks in the Michigan’s Great Lakes Tectonic Zone. Preliminary results for Lu-Hf dating of fluorite from carbonatite rocks in the Lemitar Mountains of New Mexico are also included. The analytical methodology shows promising results in each of these cases where there is sufficient parent Lu and U and measurable daughter Hf and Pb respectively.

Conference Paper

Geologic map of the northwest flank of Mauna Loa volcano, Island of Hawai‘i, Hawaii

Mauna Loa, the largest active volcano on Earth, has erupted 34 times since written descriptions became available in A.D. 1832. The most recent eruption of Mauna Loa occurred on November 27, 2022, after a 38 year hiatus; it lasted for 12 days. Some eruptions began with only brief seismic unrest, whereas others followed several months to a year of increased seismicity. Once underway, Mauna Loa’s eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows that erupted from the Southwest Rift Zone in 1950 advanced at an average rate of 9.3 kilometers per hour (5.8 miles per hour); all three lobes reached the ocean within ~24 hours. Near the eruptive vents, the flows likely traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa has great volcanic-hazard potential for the Island of Hawai‘i. Volcanic hazards on Mauna Loa can be anticipated, and risk substantially mitigated, by documenting its past activity to refine our knowledge of the hazards, and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. The map of the north and west flanks of Mauna Loa shows the distribution and relation of volcanic and surficial sedimentary deposits. It incorporates previously reported work published as generalized small-scale maps and a more detailed map. Within the mapped area, lava has flowed from three different source regions: the Northeast Rift Zone (22 percent), the summit (64 percent), and radial vents (14 percent). All three have different points of origin which, in turn, affect the flow characteristics and periodicity of activity. The map area includes the uppermost part of the NERZ and extends from the highest elevation––13,040 feet at the south end of the Kokoolau quadrangle, just below the summit caldera––to the sea northwest and west of the summit. Lava that erupts from the north and west flanks typically flows to the west, northwest, or north, depending on the vent location. Both morphologic lava flow types—‘a‘ā and pāhoehoe—are present. Pāhoehoe units tend to spread out or widen in low-slope regions, such as in the saddle regions between Mauna Loa and Mauna Kea or between Mauna Loa and Hualālai. In comparison, ʻaʻā flows generally produce narrower flow lobes that have higher relief. This map is the fifth in a series of five maps that will cover Mauna Loa volcano. NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii

Geophysical characterization of an alkaline‑carbonatite complex using gravity and magnetic methods at Magnet Cove, Arkansas, USA

The Magnet Cove alkaline‑carbonatite complex (MCC), located in the Ouachita Mountains of south-central Arkansas in the United States, hosts an extensive variety of rare rock types and critical mineral resources with physical properties (density and magnetization) that contrast significantly with the sedimentary rocks into which they have intruded. Newly acquired ground-based gravity and magnetic data were used to develop two-dimensional and three-dimensional geophysical models of the Cretaceous-aged Magnet Cove intrusive complex. The models reveal that the MCC: (1) widens out at middle crustal depths to as much 22 km across, and may reach a depth of 20 km; (2) has a total volume (exposed and subsurface) that may be over 800 km 3 ; (3) is likely connected at depth to other intrusions in the Arkansas alkaline province; and (4) has a geometry that is aligned with pre-existing structures such as the Reelfoot rift and the Ouachita orogenic belt, some of which were likely structurally controlled by the Precambrian crystalline basement and the continent-ocean transition zone buried beneath the Ouachita orogen. For the first time, the magnetic models of the MCC account for the presence of strong remanent magnetization. This results in a geophysical workflow necessary to accurately interpret magnetic anomalies over the much larger Arkansas alkaline province, its geologic and structural framework, and critical mineral potential.

Arkansas

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America

Hydrogeologic framework and conceptual model of the Red River alluvial aquifer east of Lake Texoma, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. At present (2025), the Oklahoma Water Resources Board has not established a maximum annual yield for the Red River alluvial aquifer east of Lake Texoma. To support the evaluation and determination of a maximum annual yield, a hydrogeologic framework and conceptual groundwater-flow model were developed to assess groundwater availability in the Red River alluvial aquifer east of Lake Texoma. The scope of this hydrologic investigation is the alluvium and terrace containing the Red River alluvial aquifer in Oklahoma between Lake Texoma, the Texas State line, and the Arkansas State line, an extent referred to in this report as “the eastern part of the Red River alluvial aquifer.” Parts of the alluvium and terrace extent in Arkansas and Texas are included in some analyses to address hydrologic influences from outside the aquifer’s boundaries in Oklahoma. The eastern part of the Red River alluvial aquifer in southeastern Oklahoma consists of approximately 401,280 acres of Quaternary alluvium and terrace deposits associated with the Red River and its major tributaries. Mean annual recharge to the aquifer for the 1980–2022 study period was estimated to be 8.62 inches per year, or 17.98 percent of the mean annual precipitation over the same period (47.94 inches). This mean annual recharge rate is equivalent to an inflow of approximately 288,250 acre-feet per year for the eastern part of the Red River alluvial aquifer. Recharge estimated using the Soil-Water-Balance code accounts for 98.7 percent of the conceptual-model inflows to the eastern part of the Red River alluvial aquifer. Saturated-zone evapotranspiration accounts for 11.9 percent and net streambed seepage accounts for 87.4 percent of the outflows in the conceptual model.

Arkansas, Oklahoma, Texas

Virginia and Landsat

From the shores of Jamestown and spreading north, south, and west, the lands that became the State of Virginia were some of the first in North America top experience rapid landscape change from European settlement. Imagery and data from the USGS Landsat series of satellites offer an unparalleled resource for the study, understanding, and preservation of Virginia’s land and water resources. From monitoring the health of water bodies to managing invasive species to planning for a range of climate change effects, the USGS National Land Imaging Program’s stewardship and public delivery of Landsat data have benefitted Virginians in myriad ways.

Virginia

Development of vegetation based soil quality indices for mineralized terrane in arid and semi-arid regions

Soil quality indices (SQIs) are often management driven and attempt to describe key relationships between above- and below-ground parameters. In terrestrial systems, indices that were initially developed and modified for agroecosystems have been applied to non-agricultural systems in increasing number. We develop an SQI in arid and semi-arid ecosystems of the Western US impacted by different types of geologic mineralization using the relationship between vegetation community parameters and soil abiotic and biotic properties. We analyze these relations in soils associated with three different mineralization types: podiform chromite, Cu/Mo porphyry, and acid-sulfate gold vein systems at four different sites in California and Nevada. Soil samples were collected from undisturbed soils in both mineralized and nearby unmineralized substrates as well as from waste rock and tailings. Aboveground net primary productivity (ANPP), canopy cover and shrub density were measured for the vegetative communities. Minimum data sets were developed based on correlations between the soil and vegetation parameters, refined using principal components analysis, scored using non-linear functions, and combined into an overall SQI. The indices are comprised of one or two microbial parameters and three to six abiotic parameters, the latter consisting of nutrients and metals. Given the preliminary development of this approach, the parameters and combinations to arrive at an SQI for a given site cannot at this time be correlated or compared with that of another site. This SQI approach provides a means of quantifying disturbed ecosystem recovery resulting from mining, and could be applied to other disturbances in a way that readily distills the information for potential use by land managers. However, severely disturbed areas with little to no aboveground biomass, such as unreclaimed tailings, have likely crossed an ecological threshold that precludes the use of this type of monitoring tool.

California

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular

Magnitude and frequency of floods in the Suwannee River Water Management District, Florida

Flood-frequency statistics for 2-, 5-, 10-, 25-, 50-. 100-, 200-, and 500-year recurrence intervals, based on three methods of analysis, are presented for 25 continuous-record and seven peak flow partial-record gaging stations in the Suwannee River Water Management District. The first method, for gaged stations, utilizes station records; the second method, for ungaged sites, utilizes regional regression analysis; and the third method uses a weighted combination of the station and regional values. Because the weighted values utilize two more or less independent estimates of the peak flow statistic, they are considered more accurate than the station estimates or the regression estimates alone. Also, the use of another weighting scheme to improve estimates of flood frequency statistics at ungaged sites is demonstrated. The karstic nature of much of the Suwannee River Water Management District significantly attenuates flood peaks in some streams by providing substantial subsurface storage when river stages are high. At such times, springs discharging into rivers may reverse flow temporarily and become sinks.

Florida

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon