USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Coastal 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,660 records · Page 65Linked to original sources

Preliminary ground and airborne-based geophysical mapping and modelling of an active hydrothermal system at Mammoth Lakes, California

Mammoth Lakes, California hosts a productive hydrothermal system within the seismically active south moat of Long Valley Caldera. Surficial evidence of the shallow hydrothermal system includes discrete zones of tree-kill dispersed between Shady Rest Park and the Casa Diablo Geothermal Power Plant (40 MW), as well as east of the power plant. The tree-kill areas are associated with elevated diffuse CO₂ emissions, heated ground, hydrothermal alteration, diffuse soil H₂S emissions, and gas vents. Previous mapping delineates prominent north and northwest trending structures within the south moat along the southwestern edge of the resurgent dome that may accommodate gas and fluid flow at the Shady Rest Park and Basalt Canyon Tree Kill Areas (SRTKA and BCTKA, respectively). Both tree-kill areas are also located along contacts between resurgent rhyolite, mafic lavas, and surficial deposits which may provide additional pathways for gas and fluid migration in the shallow subsurface. Characterizing structure and lithology using geophysical anomalies is critical to determining primary structural controls on the hydrothermal system and the extent of subsurface alteration at these sites. We conducted ground and airborne-based potential field geophysical surveys to map gravity and magnetic anomalies. These anomalies are then used to model subsurface geology, structure, and hydrothermal alteration. Here we present our preliminary geophysical mapping and modelling results at both tree-kill locations. Gravity and magnetic data suggest complex structural intersections are coincident with heated ground and gas emissions at the SRTKA and BCTKA. Hydrothermal systems are often observed or interpreted to exploit fault intersections which can serve as highly permeable pathways for hydrothermal fluid and gas discharge, enabling economic geothermal energy production. Geophysical mapping and modelling are an effective means of investigating such structural complexity at Mammoth Lakes due to the presence of unidentified and concealed structures.

California

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

Reservoir evolution, downstream sediment transport, downstream channel change, and synthesis of geomorphic responses of Fall Creek and Middle Fork Willamette River to water years 2012–18 streambed drawdowns at Fall Creek Lake, Oregon

Executive Summary Chapter A. Introduction Fall Creek Dam impounds Fall Creek Lake, a 10-kilometer-long reservoir in western Oregon and is operated by the U.S. Army Corps of Engineers (USACE) primarily for flood-risk management (or flood control) in late autumn through early spring months, as well as for water quality, irrigation, recreation, and habitat in late spring through early autumn. Since 2011 (water year [WY] 2012), Fall Creek Lake has been temporarily drawn down each year to facilitate downstream passage of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) through the 55-meter (m) high dam. This annual dam operation is temporary, typically lasting about 1–2 weeks from WY 2012 through 2020 (drawdown operations in WY 2022–24 have increased to more than 6 weeks). Drawdown of the reservoir results in lake levels being lowered to the elevation near the historical, pre-dam streambed. The annual streambed drawdowns of WY 2012–18 have improved fish passage and led the USACE to formally adopt streambed drawdowns as part of annual operations at Fall Creek Dam. However, temporarily lowering the lake to streambed creates free-flowing conditions in the reservoir that result in the erosion and episodic export of predominantly sand and finer-grained sediments (less than 2 millimeters [mm]) to the lower gravel-bed reaches of Fall Creek and the Middle Fork Willamette River. The introduction of large volumes of sand and finer-grain sediment into the dam-regulated reaches downstream from Fall Creek Dam prompted questions about the geomorphic responses to annual streambed drawdowns within Fall Creek Lake and downstream reaches along Fall Creek and the Middle Fork Willamette River. The U.S. Geological Survey (USGS) in partnership with USACE initiated a comprehensive geomorphic and sediment transport investigation to assess the coupled processes of reservoir erosion, sediment evacuation from Fall Creek Lake, and patterns of sediment transport and deposition in reaches downstream from the Fall Creek Dam that have resulted from annual streambed drawdowns. The purpose of this report is to systematically describe the processes of sediment erosion, transport, and deposition at Fall Creek Lake and geomorphic interactions between reaches upstream and downstream from Fall Creek Dam that relate to dam operations. Specifically, this report focuses on evaluating geomorphic responses to streambed drawdowns from WY 2012 through 2018 and placing drawdown-induced geomorphic responses within the broader context of physiographic and historical conditions and dam operations of Fall Creek and Middle Fork Willamette Rivers. Key objectives for this study were to characterize changes in reservoir morphology and substrate at Fall Creek Lake, describe the character and temporal pattern of sediment transport downstream from Fall Creek Dam, characterize geomorphic changes in channel reaches downstream from the Fall Creek Dam, and relate these data to the annual streambed drawdowns of WY 2012–18. This study uses multiple independent monitoring and measurement approaches to assess site, reach, and river-scale geomorphic responses to drawdowns to inform dam and reservoir management. Patterns and processes of reservoir evolution were assessed with geomorphic mapping and volumetric analyses of topography through comparison of multiple digital surface models (DSMs). Just downstream from Fall Creek Dam, analyses of sediment export from the reservoir focused on suspended sediment but also incorporated bedload analyses to assess sediment sizes. Geomorphic assessments downstream from the dam used reach-scale and site-scale approaches to document changes in channel morphology and substrate, including site measurements of sand and finer-grained sediment deposition and in-channel bed-material, volumetric change analyses from comparison of digital elevation models (DEMs), and repeat geomorphic mapping. Findings from this study inform river management and dam operations by providing an understanding of (1) coupled upstream-downstream geomorphic responses to the Fall Creek Lake streambed drawdowns, (2) geomorphic responses of Fall Creek Lake streambed drawdowns in comparison to drawdowns at other large dams, (3) controls on reservoir erosion and downstream geomorphic responses, and (4) implications for future hydrogeomorphic changes that may result from continued drawdowns and monitoring activities to assess those changes. Chapter B. Reservoir Morphology and Evolution Related to Dam Operations at Fall Creek Lake To understand the volume and distribution of sediment accumulation in Fall Creek Lake since dam closure in 1965, decadal-scale sedimentation patterns (spanning approximately 1965–2016) are evaluated using a combination of storage curve analyses and geomorphic mapping. Short-term (drawdown event-scale) patterns of erosion, sedimentation, and sediment export downstream are evaluated using a combination of geomorphic mapping and change detection analyses that quantify the distribution and total volume of sediment erosion and deposition within Fall Creek Lake. Geomorphic mapping of reservoir topography and analyses of historical datasets reveals four categories of landforms and sediment processes within Fall Creek Lake related to lake level operations: lacustrine sedimentation expressed in the reservoir floor, fluvial erosion and deposition within historical stream channels during streambed drawdowns, channel-like features created by erosion within the reservoir floor during streambed drawdowns, and erosion on reservoir hillslopes. Where the reservoir floor is mapped for this study as pelagic (deep water), deposition up to 3 meters (m) thick by lacustrine processes and burial of pre-dam topography with deposits thinning toward the edges of the valley floor and upstream areas of reservoir are observed. Despite over 50 years of sediment accumulation since dam construction, the main stream channels of Fall and Winberry Creeks (or reservoir thalwegs) through the reservoir are well defined, though their distinct morphology is likely influenced by a long history of recurring historical drawdowns to or near streambed since dam construction. Unregulated streamflow and sediment transport through the reservoir primarily are confined to these channels during the streambed drawdown periods. Erosional channel-like features created by drawdowns are carved through underlying, unconsolidated reservoir floor sediments and are most prominent in the lower reservoir below minimum conservation pool (the low pool elevation during winter flood season); sediment generated from the formation of these drawdown channels is more likely to be transported through and out of the reservoir than sediment deposits along the reservoir hillslopes at the valley margins that are separated from main channels by areas of low-gradient reservoir floor. Morphologic changes in the lower reservoir topography between January 2012 and November 2016 indicate overall net erosion of about 129,500 cubic meters (m 3 ). The most prominent geomorphic changes occurred along the main channels of Fall and Winberry Creeks near the Fall Creek Dam where incision, lateral migration, and slumping banks resulted in vertical and lateral adjustments to channel position, whereas most changes fell below the detectable limit on higher-elevation reservoir floor surfaces except where erosion occurred along features mapped as drawdown channels. Chapter C. Sediment Delivery from Fall Creek Lake and Transport through Downstream Reaches USGS implemented a sediment monitoring program in WY 2013–18 to evaluate the quantity and character of reservoir sediment exported from Fall Creek Lake during streambed drawdowns. Turbidity and suspended sediments were monitored annually autumn through spring to span the WY 2013–18 streambed drawdowns; however, unequal monitoring timeframes each year reduced the ability to compare results and factors affecting sediment export from the reservoir difficult between years. These data were originally measured to develop regressions and compute suspended-sediment loads (SSL). Bedload sediment monitoring from a cableway at the Fall Creek streamgage was completed in the autumn-winter of WY 2013 and 2017. The limited number of samples and presumed variability in sediment supply from the reservoir precluded construction of streamflow and bedload discharge relations to compute more than instantaneous bedload. Sand and finer-grained silts and clays were transported from the reservoir in suspension, though some coarser grains (up to 32 mm) were also mobilized and transported downstream from the dam as bedload. Observations of increased sediment transport downstream from Fall Creek Dam coincided with lake levels approaching about 3 m (10 feet [ft] or elevation 690 ft) above the streambed regulating outlets. Suspended-sediment loads computed for the full monitoring periods WY 2013–18 at the Fall Creek streamgage, located 1.4 kilometers (km) downstream from Fall Creek Dam, range from 54,700 metric tons (t) in WY 2013 to 13,900 t in WY 2018. Although the total annual SSL varied from year to year, the overall seasonal patterns of suspended sediment transport throughout each year were similar during monitoring in WY 2013-18. Suspended-sediment loads were low prior to the drawdown, then increased rapidly as lake levels lowered and approached the streambed. In the weeks following the drawdown period, as pool levels were increased, SSL remained slightly elevated above pre-drawdown levels but generally declined through the following winter and spring except during streamflow-driven pulses of suspended-sediment transport. WY 2013 had the greatest total computed SSL for each streambed drawdown and partial-year monitoring period. SSL computed for the partial-year period have generally decreased since WY 2013 and have varied by about 6,800 t with the exception of WY 2014. WY 2014 SSL reflects anomalously low sediment export due to low streamflows and freezing conditions that stabilized reservoir floor deposits. Bedload measurements in the short 1.4-km reach between Fall Creek Dam and the Fall Creek streamgage showed an inverse correlation between bedload transport rates and discharge, which probably reflects diminishing supply of coarse-sized sediment. Sand was more abundant (60–100 percent) than gravel in bedload samples confirming sand and finer-grained sediment dominated sediment evacuated from the reservoir during streambed drawdowns at Fall Creek Lake. Chapter D. Geomorphic Responses to Fall Creek Lake Streambed Drawdowns Downstream from Fall Creek Dam In the days, weeks, and months following streambed drawdown operations at Fall Creek Dam through WY 2018, sites downstream from the dam displayed a variety of geomorphic responses to reservoir sediment delivery within the main channel and overbank areas. Evaluation of streambed elevations at two streamgages located 1.4 km downstream from the dam on Fall Creek and 16.3 km downstream from the dam on the Middle Fork Willamette River indicated the effects of drawdown sediment on bed elevations were modest and transient. Repeat particle size measurements (October 2015 and September 2016) at five sites along Fall Creek and the Middle Fork Willamette River showed similar grain-sized distributions that do not reveal substantial deposition of fine-grained sediment related to the WY 2016 streambed drawdown. Altogether, these findings indicate that transport capacity in the main, low-flow channels of Fall Creek and Middle Fork Willamette River during WY 2012–18 was sufficient to mobilize and evacuate reservoir sediments from streambed drawdowns or other bank material and tributary sources. However, other monitoring for this study indicate low-velocity zones in off-channel areas are prime locations for sand and finer-grain sediment deposition. Patterns of overbank sediment accumulation indicate that the magnitude and timing of overbank deposition on bars and low-elevation floodplain varies with proximity to the dam, geomorphic setting, streamflows, and other factors. Sand and finer-grained reservoir sediments carried as suspended-sediment load in the reaches downstream from Fall Creek Dam were deposited in overbank areas as observed with clay-horizon markers during WY 2016–17. Overbank deposition quantified with Geomorphic Change Detection (GCD) software evaluated landform-scale patterns of erosion and deposition using repeat light detection and ranging (lidar) surveys at two sites in the Upper Fall Creek reach and one site in the Jasper reach for 3 years (2012–15) and one site in the Clearwater reach for 6 years (2009–15). Deposition thickness and spatial patterns from the GCD analysis were variable; some sites had dispersed but measurable deposition while at others, deposition was highly localized and exceeded 1 m in depth. Patterns of overbank deposition illustrate interactions among bar morphology, local hydraulics, and suspended-sediment transport dynamics that can create patches of highly localized deposition. The measured deposition at the two Fall Creek GCD sites likely resulted from reservoir sediments released from Fall Creek Lake during streambed drawdowns in WY 2016 and 2017 because the limited sediment inputs from bank material (geomorphically laterally stable reach) or tributaries (no significant tributaries) provided few other sediment sources. On the Middle Fork Willamette River, observed patterns of overbank deposition could reflect sediment sourced from upstream tributaries, bank erosion, or Fall Creek Lake streambed drawdown operations. Despite the introduction of several thousand tons of reservoir sediment delivered from the Fall Creek Lake streambed drawdowns to below-dam river corridors, reach-scale mapping of channel features downstream from Fall Creek Dam shows minimal evidence of changes in channel planform or landforms that can be attributed to a drawdowns in WY 2012–16. On Upper Fall Creek reach, widespread increases in gravel bars or other in-channel sediment did not result from the five streambed drawdowns. The main changes attributable to sediment releases from Fall Creek Lake were localized increases in vegetated bar area, particularly on channel margin areas where sand and finer-grain sediment was deposited and rapidly colonized by vegetation. The area of mapped secondary water features decreased between 2005 and 2016, but that may be due to lower discharges depicted in the 2016 aerial photographs and less mapped area of inundation. Primary changes along the Lower Fall Creek reach include a 6.4 percent decrease in area of secondary water features between 2011 and 2016, and a nearly twofold increase in the area of unvegetated bars. Immediately downstream from the Fall Creek confluence, there were negligible changes in the location and areas of vegetated bars and the main wetted channel between 2005 and 2016, and local increases in bar area cannot be attributed solely to deposition of reservoir sediments from Fall Creek Lake because (1) areas along the Middle Fork Willamette River just upstream from the Fall Creek confluence display similar type and magnitude of changes and (2) some of the increases at the confluence area pre-date the drawdowns. The cumulative effect of sediment releases from Fall Creek Lake streambed drawdowns from WY 2012 to 2016 on downstream channel planform and landforms are modest compared to the river-scale transformations and planform changes that occurred in the decades following dam construction. Chapter E. Discussion of Geomorphic Responses of Fall Creek and Middle Fork Willamette River to Streambed Drawdowns at Fall Creek Lake Multiple aspects of Fall Creek Dam infrastructure and operations exert first-order controls on the magnitudes of reservoir erosion that occur during the streambed drawdowns and ultimately determine the sediment delivery to downstream reaches. Key aspects of the dam and its operations that are most relevant to assessing geomorphic responses to streambed drawdowns include the (1) dam infrastructure, including configuration and size of regulating outlets and their proximity to the streambed which dictates the capacity and competence of the river to deliver sediment to downstream reaches and mode of sediment transport as suspended-sediment load or bedload; (2) frequency of historical drawdowns and long-term, year-round dam operations and lake level management, which partly dictate reservoir morphology and locations and magnitudes of readily erodible materials; (3) dam operations and hydroclimatic conditions during the streambed drawdown (including length of the drawdown and streamflows entering the reservoir), which directly control the timing, duration and magnitude of reservoir erosion and sediment evacuation; and (4) dam operations following the streambed drawdown operation that regulate streamflows (and thereby sediment transport conditions) downstream of Fall Creek Dam which primarily reflect interactions between hydroclimatic conditions and flood control operations. Patterns of sediment erosion and evacuation observed in this study at Fall Creek Lake from WY 2012–18 suggest that reservoir erosion during annual streambed drawdowns may remain similar or decrease in future years assuming (1) annual streambed drawdown operations are implemented in similar manner as the WY 2012–18 drawdowns (in terms of duration, late autumn or early winter implementation, rate of pool-level lowering to reach streambed, and other factors), (2) streambed drawdowns coincide with similar conditions as were observed WY 2012–18 (similar sediment yield into reservoir, low reservoir inflows, limited precipitation, moderate air temperature), and (3) no major geomorphic changes in the main reservoir channels of Fall and Winberry Creeks occur (for example, channel avulsion). Under such conditions, it is hypothesized that the stream channel within the reservoir would achieve a quasi-equilibrium state with respect to annual influx and export of sediment and aided by the substantial amount of in-channel bedrock, will remain laterally stable without erosion across reservoir deposits. Patterns of sediment transport measured at the Fall Creek streamgage downstream from Fall Creek Dam provide insight into the potential effects of future streambed drawdowns at Fall Creek Lake. Analyses of suspended sediment measured in WY 2013–18 show a major reduction in suspended-sediment loads between WY 2013 and later years, indicating streamflows transporting sediment through the reservoir to downstream reaches during streambed drawdowns have become supply limited. The 6-year suspended-sediment monitoring and sampling program is insufficient to make predictions about future sediment transport conditions because of uneven monitoring periods and varying controls on reservoir sediment erosion. It is likely that future suspended-sediment loads will be variable but similar to those observed in WY 2015–18 if operational, climatic, and geomorphological factors remain similar to those monitored WY 2015–18. Suspended-sediment loads downstream from Fall Creek Lake will likely remain highest when regulating outlets are fully open and Fall Creek is free flowing with the reservoir fully drained with little to no residual pool. Over time, it is possible that the suspended-sediment loads would reflect mobilization of reservoir sediment deposited in the previous year rather than erosion of sediment deposited years or decades earlier. Bedload is likely to remain a small fraction of the total sediment load evacuated from the reservoir and is relatively modest compared with pre-dam bedload transport rates because most coarse sediment remains trapped by the dam. If sediment releases from Fall Creek Lake and ensuing streamflow conditions follow a similar pattern in the future as was assessed in this study spanning WY 2012–18, near-term geomorphic adjustments downstream of the dam are expected to be modest. Barring major operational, climatic, and geomorphological changes, local site-scale deposition on bars, overbank areas, or off-channel features that persists several months after the streambed drawdown will likely continue to be highly variable, ranging from negligible to several centimeters of deposition. At the landform-scale, low velocity areas nearest to Fall Creek Dam will likely continue to undergo rapid deposition immediately during and after a streambed drawdown event, similar to patterns observed for WY 2012–18. Some of the sediment entering these off-channel features and margin areas may be temporarily stored, then later remobilized and dispersed farther downstream. But if newly deposited sediment persists through the following spring, there is a greater likelihood that local vegetation will establish, reinforce deposited material, and trap sediment during later drawdowns. The reach-scale geomorphic changes may become more apparent if (1) streambed drawdowns continued for several decades, and geomorphic changes were measured at decadal scales or (2) the amount of sediment introduced to downstream reaches substantially increased and (or) sediment transport capacity decreased. The continued streamflow regulation of Fall Creek Dam after sediment releases provides an opportunity to strategically manage streamflows during and after the streambed drawdowns to minimize downstream sediment impacts and ensure other operational thresholds are satisfied. This study provides a comprehensive foundation of datasets and geomorphic analyses to inform dam operations at Fall Creek Lake, monitor sediment transport downstream, and consider operational schemes for future drawdowns. The datasets from this study also provide baselines of sediment transport and geomorphic conditions to assess future changes in reservoir and downstream environments. Future monitoring could be tailored to address specific questions regarding the long-term geomorphic effects of streambed drawdowns on fluvial habitats, flood hazards, cultural resources, or downstream water quality. Future monitoring activities could focus on the relevant geomorphic processes and spatial domains within the three categories used for this study: (1) reservoir erosion and net sediment evacuation, (2) sediment delivery to downstream reaches, including magnitude and temporal pattern of sediment transport, and (3) geomorphic responses of downstream reaches to sediment delivery. Specifically, high priority future monitoring activities could include: Repeat topographic or photographic surveys in the reservoir to characterize changes occurring within individual drawdowns, to quantify sediment export, to determine temporal changes in reservoir storage, and to identify locations of erosion and deposition. Continuous, year-round turbidity monitoring supplemented with suspended-sediment measurements at a streamflow-gaging station immediately downstream from the dam to quantify sediment export. Repeat geomorphic monitoring, mapping, or modeling in downstream reaches to track changes in channel and over bank features using a combination of site- and reach-scale monitoring approaches. This could support assessments of sediment deposition and ensuing vegetation encroachment on flood hazards and habitats and examine how sediment transport and depositional processes may be affected by different sediment supply, streamflow, or dam management scenarios.

Oregon

Pesticides in surface waters of the Upper Colorado River Basin, Colorado, 1996-98

Forty-four river, stream, and drain sites in the Upper Colorado River Basin in Colorado were sampled during 1996-98 to determine the occurrence and distribution of pesticides in the basin. In a fixed-station study, 57 surface-water samples were collected from October 1996 through January 1998 at four sites. Each site was sampled approximately monthly for up to a year, with more frequent sampling during the spring and summer growing season. In a synoptic study, surface-water samples were collected at 43 sites in the agricultural areas of the Grand Valley and the Uncompahgre River Valley in May 1998. Each site was sampled once. Pesticide concentrations generally were low and varied seasonally and across land uses. Thirty-five pesticides were detected at least once in the two studies, and 11 pesticides accounted for almost 82 percent of the detections. Herbicides were more commonly detected than insecticides, and only the herbicides alachlor and atrazine were detected in more than 50 percent of the samples. Carbofuran was the most commonly detected insecticide and was detected in 19 percent of the samples. Pesticide detections increased and were measured at higher concentrations in the summer months and at the agriculture sites. All pesticide concentrations were less than drinking-water standards, and most complied with human-health advisories and criteria for the protection of freshwater aquatic life. Pesticides in the Upper Colorado River Basin in Colorado were detected slightly less frequently and generally at lower concentrations than in 20 National Water-Quality Assessment Program study-unit basins that collected water-quality data from 1992 through 1996. Results from surface-water sampling conducted during 1996?98 in the Upper Colorado River Basin in Colorado indicate that beneficial uses of water were not being impaired by the presence of pesticides in surface waters in the basin.

Colorado

Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in-situ distribution of fault damage remain limited and along-strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub-bottom profiles. The combination of high-resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi-trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

California

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California

Hydrogeologic framework of the Mountain Home area, southern Idaho

In the arid western Snake River Plain around the City of Mountain Home, Idaho, declining groundwater levels concern agricultural, municipal, and other water users who rely on groundwater for sustenance because surface-water resources are limited. The U.S. Geological Survey developed this hydrogeologic framework to provide an updated characterization of groundwater resources in the western Snake River Plain around the City of Mountain Home. The hydrogeologic framework comprises: (1) a conceptual description of hydrogeologic units, (2) a three-dimensional hydrogeologic model and borehole database, (3) a map of groundwater levels and change, and (4) a discussion of groundwater occurrence and movement within the study area. Hydrogeologic units were defined based on existing literature and the borehole database compiled for this study; the five hydrogeologic units are granite, rhyolite, basalt, fine-grained sediments, and coarse-grained sediments. Each unit can bear water, but the main regional aquifer in the study area occurs in the basalt and fine-grained sediment units with depth to water ranging from 150 to 765 feet. A perched groundwater zone near the City of Mountain Home is primarily hosted in basalt and used domestically with most depths to water ranging from 30 to 100 feet. Interflow zones, scoria, and vertical fractures create heterogeneity within the basalt hydrogeologic unit that exerts strong control on groundwater movement, creating horizontal perching conditions and zones of enhanced vertical conductivity that facilitate downward groundwater percolation. In the fine- and coarse-grained sediments and rhyolite units, inferred faults both impede and enhance groundwater movement. The borehole database was constructed by digitizing 540 well-driller reports and was used to build a three-dimensional hydrogeologic framework model which reasonably represents the spatial distribution of hydrogeologic units in the study area. Generally, fine-grained sediments underlie much of the study area, with basalt concentrated in the central and western study area and rhyolite and granite in the uplands to the north. Groundwater levels were measured in 180 wells in March and November 2023; these data were used to develop water-table contour maps and describe groundwater-level change over an irrigation season. Groundwater generally flows south-southwest to the Snake River and groundwater levels declined across most of the study area (from 0.03 to 22.01 feet) between spring and autumn 2023, which is consistent with long-term declines in the Cinder Cone Butte Critical Groundwater Area and Mountain Home Groundwater Management Area. Groundwater levels rose (0.6 to 15.44 feet) over the irrigation season in most wells in the perched groundwater zone near the City of Mountain Home and near the Snake River, indicating the importance of surface-water recharge to groundwater in areas where surface water irrigation occurs. In aggregate, this hydrogeologic framework provides an updated characterization of and new insights into groundwater resources in the study area to help inform water resources management.

Idaho

Widespread anhydrite saturation in Laramide-age arc magmas of southwestern USA

Anhydrite is considered a rare mineral phase in magmas, with only ∼33 documented occurrences worldwide. However, anhydrite readily decomposes in the near-surface environment, making it difficult to recognize its former presence in rocks collected at or near Earth’s surface. In such samples, only small anhydrite inclusions fully shielded within other minerals can have survived. During a recent field trip to the southwestern USA, we sampled 17 Laramide-age (ca. 40−80 Ma) magma systems, most of which are associated with porphyry copper deposits. A systematic search for anhydrite inclusions preserved within apatite, amphibole, plagioclase, and quartz phenocrysts in ∼100 rock samples by optical microscopy and Raman spectroscopy revealed that each of these 17 magma systems was at least temporarily anhydrite-saturated. Also, most previously identified magmatic anhydrite-bearing intrusions are associated with porphyry copper deposits, and both intrusive and volcanic rocks containing magmatic anhydrite show high Sr/Y ratios. These observations suggest that anhydrite saturation and porphyry copper formation are linked via magma fractionation at high pressure. Compared to average arc magmas, anhydrite-bearing magmas are unusually oxidized and sulfur-rich and seem to also be unusually water-rich. Hence, our preferred interpretation is that magma generation and/or fractionation at high pressure promotes the formation of superhydrous and oxidized magmas, which in turn promotes high sulfur contents and ultimately the precipitation of anhydrite. The high mineralization potential of these magmas does not need to result from their high sulfur content but could be caused by other properties of high-pressure magmas.

Arizona, New Mexico

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone

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

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota

Mafic alkaline magmatism and rare earth element mineralization in the Mojave Desert, California: The Bobcat Hills connection to Mountain Pass

Occurrences of alkaline and carbonatite rocks with high concentrations of rare earth elements (REE) are a defining feature of Precambrian geology in the Mojave Desert of southeastern California. The most economically important occurrence is the carbonatite stock at Mountain Pass, which constitutes the largest REE deposit in the United States. A central scientific goal is to understand the genesis of the carbonatite ore body in the context of widespread REE-rich igneous activity. A swarm of mafic alkaline (shonkinite) dikes has been mapped and sampled at Bobcat Hills, 65 km southeast of the Mountain Pass mine. Whole-rock geochemistry and zircon geochronology demonstrate a clear affinity to the ca. 1.4 Ga Mountain Pass intrusive system. Bobcat Hills dikes have comparably high REE concentrations (La ∼1,000× chondritic) and an error-weighted mean 207 Pb/ 206 Pb zircon crystallization age of 1,426 ± 2 Ma (2 σ ). Unlike the alkaline intrusions at Mountain Pass, which have abundant inherited zircon from Paleoproterozoic basement rocks and crustally influenced oxygen isotope compositions (δ 18 O zircon = 6.5–7.5‰), the Bobcat Hills dikes lack any evidence of crustal assimilation and have oxygen isotope values that overlap a mantle range (Bobcat Hills average δ 18 O zircon = 5.6 ± 0.3‰). The dikes were a high-temperature, early center of mafic alkaline magmatism in the Mojave Desert that serve as a snapshot of melt generation from a spatially extensive, metasomatized mantle source. We propose that modification of the crust over many tens of Myr at Mountain Pass created an environment that favored crustal assimilation and enabled ascent of late-stage, REE-rich carbonatite magmas.

California

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States

Geology and geochemistry of Jurassic plutonic rocks, Baboquivari Mountains, south-central Arizona

Among the plutons of the Jurassic magmatic arc segment in southern Arizona, northern Sonora, and southern California, two distinctive rock types dominate, occurring together in numerous mountain ranges: porphyritic, titanite-bearing hornblende-biotite granodiorite; and biotite leucogranite. These characteristic rocks are particularly well preserved and well exposed around Kitt Peak, in the Baboquivari Mountains, southern Arizona. The Middle to early Late Jurassic Kitt Peak Plutonic Suite (KPPS) comprises three units: Aguirre Peak Quartz Diorite (APQD;<170 Ma, >=165 Ma, U-Pb), chiefly hornblende mesodiorite; Kitt Peak Granodiorite (KPGD; 165 Ma), hornblende-biotite granodiorite and monzogranite; and Pavo Kug Granite (PKG; 159 Ma), equigranular biotite leucogranite. This range also hosts a fourth, regionally unusual, Jurassic plutonic unit: the Baboquivari Peak Perthite Granite, related to post-arc, Late Jurassic crustal extension. Binary variation diagrams, incompatible element patterns, and REE spectra indicate APQD and KPGD are consanguineous. APQD basaltic mesodiorite evidently approximates the parental magma for the APQD–KPGD series. More evolved rocks of this series evidently formed by combinations of fractional crystallization and mixing with or assimilation of crustal material. PKG is less directly related to APQD and KPGD than these two are to one another.

Arizona, Sonora

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive

Technique for simulating peak-flow hydrographs in Maryland

The efficient design and management of many bridges, culverts, embankments, and flood-protection structures may require the estimation of time-of-inundation and (or) storage of floodwater relating to such structures. These estimates can be made on the basis of information derived from the peak-flow hydrograph. Average peak-flow hydrographs corresponding to a peak discharge of specific recurrence interval can be simulated for drainage basins having drainage areas less than 500 square miles in Maryland, using a direct technique of known accuracy. The technique uses dimensionless hydrographs in conjunction with estimates of basin lagtime and instantaneous peak flow. Ordinary least-squares regression analysis was used to develop an equation for estimating basin lagtime in Maryland. Drainage area, main channel slope, forest cover, and impervious area were determined to be the significant explanatory variables necessary to estimate average basin lagtime at the 95-percent confidence interval. Qualitative variables included in the equation adequately correct for geographic bias across the State. The average standard error of prediction associated with the equation is approximated as plus or minus (+/-) 37.6 percent. Volume correction factors may be applied to the basin lagtime on the basis of a comparison between actual and estimated hydrograph volumes prior to hydrograph simulation. Three dimensionless hydrographs were developed and tested using data collected during 278 significant rainfall-runoff events at 81 stream-gaging stations distributed throughout Maryland and Delaware. The data represent a range of drainage area sizes and basin conditions. The technique was verified by applying it to the simulation of 20 peak-flow events and comparing actual and simulated hydrograph widths at 50 and 75 percent of the observed peak-flow levels. The events chosen are considered extreme in that the average recurrence interval of the selected peak flows is 130 years. The average standard errors of prediction were +/- 61 and +/- 56 percent at the 50 and 75 percent of peak-flow hydrograph widths, respectively.

Maryland

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land