USGS ScienceSearch

place

United STates

United STates: explore 1092 source-linked works published from 1918 to 2026, with original documents and citations.

This collection is a preview while coverage and quality are evaluated.

Search within this collection

Coverage and selection

Includes records with this source-supplied label or an explicit phrase match in their metadata. Matches indicate a mention, not proof that a paper uses a method or tests a material. Source versions are consolidated by DOI.

Sources: usgs. Collection updated 2026-09-15. Counts describe this index, not the complete source archives.

Projected sea-level rise and high tide flooding at Big Cypress National Preserve, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at Big Cypress National Preserve, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

Assessment of long-term changes in surface-water extent within Klamath Marsh, south-central Oregon, 1985–2021

The annual maximum extent of surface water in Klamath Marsh has naturally fluctuated in response to periods of wet and dry conditions in the surrounding basin. Field observations during the 2010s indicate that the annual maximum extent of surface water has been declining and the marsh is not responding to hydrologic inputs as it had historically. This report describes the results of a hydrologic evaluation of Klamath Marsh to characterize and understand multi-year declines in the surface-water extent and increased intermittency of streamflow exiting the marsh. Landsat imagery collected during 1985–2021 was processed to create a time series of annual maximum surface-water extent to assess changes in surface-water inundation within the marsh. A 50-percent decrease in the mean surface area of annual total open-water extent (TOWE) during the latter half of the study period (2003–21) compared to the first half (1985–2003) was observed in this 37-year time-series dataset. The change in open-water extent was offset by a corresponding increase in dry land in the marsh. Time series of streamflow, groundwater level, total annual precipitation, annual mean temperature, and anthropogenic water use and water management were compiled and evaluated to improve understanding of the factors affecting TOWE. Statistically significant downward trends in the regional groundwater table and streamflow into and out of the marsh were identified as well as statistically significant upward trends in annual mean temperature. Statistically significant correlations among TOWE, streamflow, and groundwater level also were identified. The decreasing trends could not be attributed to changes in total annual precipitation or changing anthropogenic groundwater use within the study area. Declines in the open-water extent of Klamath Marsh since 2000 principally are due to a decoupling of the groundwater and surface-water system beneath the marsh because of regional declines in groundwater level. Regional increases in air temperature and the reestablishment of more than 55,000 acres of forested land within the study area have likely contributed to increasing evapotranspiration, leaving less water available for groundwater recharge and stream base flow and resulting in basin-wide declines in streamflow and groundwater levels.

Oregon

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska

Collision structures of the Prince William terrane and Chugach terrane docking along the Shumagin and Unimak convergent margins, Alaska, USA

Western Alaska’s convergent margins are composed of tectonostratigraphic terranes. On land, terrane assembly is recognized along boundaries or sutures between neighboring geologic elements with distinctly different origins. In marine areas where rock outcrops are covered by sediment, recognizing terrane sutures is problematic. A fault in seismic dip line 5 of the ALEUT project has been interpreted as a terrane suture. It is imaged intermittently down to the 30+-km-deep plate interface. Processing of ALEUT strike line 7 revealed the suture at ~18 km depths extending 300 km along the margin. Upper structures in line 5 are like the structures of adjacent seismic transects where imaging is only 8−10 km deep. They were previously not recognized as the upper reaches of terrane sutures and show structural details obscured at greater depths. The composite data are the basis for a simple tectonic model of terrane docking.

Alaska

Interactive effects of salinity and hydrology on radial growth of bald cypress (Taxodium distichum (L.) Rich.) in coastal Louisiana, USA

Tidal freshwater forests are usually located at or above the level of mean high water. Some Louisiana coastal forests are below mean high water, especially bald cypress ( Taxodium distichum (L.) Rich.) forests because flooding has increased due to the combined effects of global sea level rise and local subsidence. In addition, constructed channels from the coast inland act as conduits for saltwater. As a result, saltwater intrusion affects the productivity of Louisiana’s coastal bald cypress forests. To study the long-term effects of hydrology and salinity on the health of these systems, we fitted dendrometer bands on selected trees to record basal area increment as a measure of growth in permanent forest productivity plots established within six bald cypress stands. Three stands were in freshwater sites with low salinity rooting zone groundwater (0.1–1.3 ppt), while the other three had higher salinity rooting zone groundwater (0.2–4.9 ppt). Water level was logged continuously, and salinity was measured monthly to quarterly on the surface and in groundwater wells. Higher groundwater salinity levels were related to decreased bald cypress radial growth, while higher freshwater flooding increased radial growth. With these data, coastal managers can model rates of bald cypress forest change as a function of salinity and flooding.

Louisiana

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado

Characterization of the water resources of the Pamunkey River watershed in Virginia—A review of water science, management, and traditional ecological knowledge

In central Virginia, the Pamunkey Indian Tribe and Reservation are facing increasingly complex water resource issues related to quantity and quality. Documentation of surface-water, groundwater, water quality, land subsidence, sea-level rise, and river ecology issues in the Pamunkey River watershed and incorporation of traditional ecological knowledge into these research topics may improve understanding of the water resources broadly. This report summarizes the relevant traditional ecological knowledge and scientific literature to elucidate gaps in the total combined knowledge of a suite of water science topics concerning the Pamunkey River watershed. This suite of water science topics includes some of the issues that the Pamunkey Indian Tribe and Reservation are facing within the watershed: fragmentation of streams and management of streamflow, flooding, degrading water quality, groundwater extraction, relative sea-level rise, rapid changes in ecosystems, loss of ecological and biological diversity, and climate change. Gaps in total combined knowledge are pervasive throughout each of these water resource topics. Identifying these knowledge gaps can help inform future research and management strategies.

Virginia

Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States

The National Water Model (NWM) provides critical analyses and projections of streamflow that support water management decisions. However, the NWM performs poorly in lower-elevation rivers of the western United States (US). The accuracy of the NWM depends on the fidelity of the model inputs and the representation and calibration of model processes and water sources. To evaluate the NWM performance in the western US, we compared observations of river water isotope ratios ( 18 O / "> 16 O and 2 H / "> 1 H expressed in δ notation) to NWM-flux-estimated (model) river reach isotope ratios. The modeled estimates were calculated from long-term (2000–2019) mean summer (June, July, and August) NWM hydrologic fluxes and gridded isotope ratios using a mass balance approach. The observational dataset comprised 4503 in-stream water isotope observations in 877 reaches across 5 basins. A simple regression between observed and modeled isotope ratios explained 57.9 % ( δ 18 O) and 67.1 % ( δ 2 H) of variance, although observations were 0.5 ‰ ( δ 18 O) and 4.8 ‰ ( δ 2 H) higher, on average, than mass balance estimates. The unexplained variance suggest that the NWM does not include all relevant water fluxes to rivers. To infer possible missing water fluxes, we evaluated patterns in observation–model differences using δ 18 O diff ( δ 18 O obs − δ 18 O mod ) and d diff ( δ 2 H diff - 8 ⋅ δ 18 O diff "> ). We detected evidence of evaporation in observations but not model estimates (negative d diff and positive δ 18 O diff ) at lower-elevation, higher-stream-order, arid sites. The catchment actual-evaporation-to-precipitation ratio, the fraction of streamflow estimated to be derived from agricultural irrigation, and whether a site was reservoir-affected were all significant predictors of d diff in a linear mixed-effects model, with up to 15.2 % of variance explained by fixed effects. This finding is supported by seasonal patterns, groundwater levels, and isotope ratios, and it suggests the importance of including irrigation return flows to rivers, especially in lower-elevation, higher-stream-order, arid rivers of the western US.

Hydrology and Earth Systems Science

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

A comparative analysis of OpenET for evaluating evapotranspiration in California almond orchards

The almond industry in California faces water management challenges that are being exacerbated by droughts, climate change, and groundwater sustainability legislation. The Tree-crop Remote sensing of Evapotranspiration eXperiment (T-REX) aims to explore opportunities to improve precision irrigation management for woody perennial cropping systems. Almond orchards in the California Central Valley were equipped with eddy covariance flux measurements to evaluate satellite remote sensing-based evapotranspiration (RSET) models. OpenET provides high-resolution (30-m spatial and daily temporal) RSET data, synthesizing decades of research for practical water management. This study provides an evaluation of OpenET performance at six almond sites covering a large range in soils, age, and variety. It also compares OpenET ensemble evapotranspiration (ET) data with applied irrigation and precipitation records over an additional 148 almond orchards located in the Central Valley of California. Results show OpenET models, including the ensemble ET value, produced reasonable and actionable ET values, with overall coefficient of determination (R 2 ) and mean absolute error values of 0.73- and 0.95-mm d −1 at the daily time step, respectively. However, given the temporal sampling of Landsat (8-day revisit) and the interpolation methods used, the assessed ET models had difficulty in capturing short-term variability in almond ET; for example, the rapid decline in measured ET observed as a response to lack of irrigation preceding and during almond harvest. The study also drew attention to the spatial complexity in scenarios where irrigated orchards are surrounded by hot/dry areas, causing discrepancies between measured and modeled ET values. In comparison with irrigation records, OpenET ensemble ET was capable of quantifying water input (applied irrigation + precipitation) in almond orchards to within 13 % when evaluating monthly data. Initial results presented here reinforce the idea that RSET models, such as in OpenET, are powerful tools, yet their application requires nuanced understanding and careful consideration of local conditions.

California

New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States

Four new 10 Be- 26 Al isochron burial ages ranging from 4.4 to 2 Ma on ancestral Colorado River deposits in the lower Colorado River corridor (LCRC) help constrain the river’s evolution during the Pliocene and early Pleistocene. They help fill a gap between previous work that focused on older and younger deposits: Older dated deposits include the 5 Ma Bouse Formation, which records the integration of the Colorado River through a series of preexisting basins to the Gulf of California and the ca. 4.5 - 3.5 Ma Bullhead Alluvium, a 200 to 300 m thick aggregational package that immediately followed integration. The much younger, 100 - 70 ka, Chemehuevi Formation is another major aggradation package mapped throughout the LCRC. The new burial ages on the facies of Santa Fe Railway (4.37 ± 0.71 Ma), boulder conglomerate of Bat Cave Wash (2.12 ± 0.26 and 2.05 ± 0.31 Ma), and the Palo Verde alluvium (3.03 ± 0.26 Ma) partially fill in a 3.5 M.y. gap between the deposition of the Bullhead Alluvium and the Chemeheuvi Formation and document the timescales over which the Colorado River was able to remove the Bullhead aggradational package and initiate newer and smaller aggradational pulses.

Arizona, California, Nevada

The 17 January 1994 Northridge, California, earthquake: A retrospective analysis

The 17 January 1994 Northridge, California, earthquake was a watershed event, with far-reaching societal and scientific impacts. The earthquake, which occurred in the early days of both broadband seismic networks and the Internet, spurred advances in seismic monitoring, real-time systems, and development of data products. Motivated by the 30th anniversary of the earthquake, we present a brief retrospective of the earthquake and its impact, and reconsider both ground motions and the aftershock distribution using modern tools and best-available data. With improvements in instrumentation and analysis methodology, recent earthquakes continue to reveal increasing complexity of ground motions, fault systems, and earthquake ruptures. Even in the absence of data from state-of-the art instrumentation, a retrospective consideration of ground motion data from the Northridge earthquake reveals complexities beyond what could be characterized (and modeled) thirty years ago. Aftershock relocations for both the 1971 Sylmar and 1994 Northridge earthquakes also reveal an updated view of fault complexity. Our study does provide a cautionary tale regarding legacy data sets and research results that are not easily accessible, which can result in discrepancies between catalog data and products from best-available science. We also briefly describe outreach products produced as part of the anniversary commemoration.

California

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida

Projected sea-level rise and high tide flooding at De Soto National Memorial, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at De Soto National Memorial, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California

Arsenic toxicity is an environmental health problem. Levels of arsenic in surface waters at some locations in the Klamath River Basin in southern Oregon and northern California can exceed the U.S. Environmental Protection Agency (EPA) standard for drinking water. There are both anthropogenic and natural sources of arsenic. The Klamath River Basin consists primarily of volcanic deposits and contains an underground geothermal system with hot springs and warm water wells, all known natural sources of arsenic. Anthropogenic sources of arsenic are related to the agricultural use of herbicides, fungicides, and insecticides. Surface water arsenic levels can also be affected by fertilizer amendments, evaporative concentration, oxygen-level depletion, and various geochemical transformations that can increase arsenic mobilization. In this study by the U.S. Geological Survey and the Bureau of Reclamation, dissolved concentrations of arsenic, copper, and lead were measured in surface waters at 39 unique sites within the upper portions of the Klamath River Basin between 2018 and 2022. In every year, except 2022, sites were sampled four times between April and November. Surface-water arsenic concentrations varied up to four-orders of magnitude among sites. Median arsenic concentration was lowest at Cherry Creek (0.03 micrograms per liter [μg/L]) and highest at Wood Kimball Spring (36.7 μg/L), two sites located north of Upper Klamath Lake. The highest arsenic concentrations (17.4±4.9 μg/L, n =3) were found in drain sites (defined here as a waterbody returning used irrigation water) while the lowest arsenic concentrations were found in an artesian well (0.8 μg/L, n =1). The elevated arsenic concentrations of the drain sites suggest that arsenic might be concentrated or mobilized by agricultural activities, water re-use practices, and (or) by geochemical processes occurring around water stored in drains (that is, in the water column and across sediment water boundaries). A source of arsenic in drain water in the Klamath Strait Drain area includes water used for irrigation originating from Ady Canal. Other potential sources include groundwater, geothermal water, and local soils and sediments. Seasonal differences in surface-water arsenic concentrations were detected at 13 sites, 10 of which had higher arsenic concentrations in summer than in either spring or fall. The sites sampled around Upper Klamath Lake, the impounded rivers, one of the two canal sites, and 5 of the 14 river sites had higher surface-water arsenic concentrations in the summer than in either spring or fall. Surface-water arsenic concentrations from groundwater sources (that is, springs and in the artesian well) did not vary significantly among seasons (p-values greater than 0.1). Median surface-water concentrations of copper and lead ranged from 0.03 to 3.7 μg/L, and from 0.013 to 0.175 μg/L ( n =2–18), respectively. Dissolved concentrations of both metals were below acute toxicity endpoints reported by the EPA for freshwater animals. Surface-water arsenic concentrations varied independently from corresponding changes in surface-water lead or copper concentrations. However, arsenic concentrations measured in bed-sediment samples collected from a subset of sites located north of Upper Klamath Lake correlated strongly and significantly with the corresponding sedimentary lead concentrations ( p =0.015). Aqueous arsenic speciation measured in a subset of sites in 2019 and 2022 showed that all the arsenic existed as arsenic (V), the most oxidized arsenic species, and presumably, the least toxic. The highest proportions of arsenite (As(III)), the presumably most toxic arsenic species, relative to total arsenic concentrations were found at drain sites. Our assessment of dissolved arsenic concentrations in various surface-water bodies in the Upper Klamath River Basin reveals geographical areas of consistently low (below 2.1 μg/L), moderate (below 10 μg/L) and high (above 10 μg/L) surface-water arsenic concentrations. South of Upper Klamath Lake, surface-water arsenic concentrations were consistently higher than 20 μg/L at two drain sites located in an area of predominant agricultural land use with extensive water re-use practices. North of Upper Klamath Lake, surface-water arsenic concentrations greater than 20 μg/L were consistently measured at sites with limited nearby agricultural activities, suggesting a geogenic source. The consistently high arsenic levels from the Wood River at Jackson F. Kimball State Park, Fort Creek, and Crooked Creek, which are sites located at or near headwater spring sources, suggest a natural background source of arsenic. Water flowing downstream from this area could be a potential source of arsenic to Upper Klamath Lake and the Upper Klamath River.

California, Oregon
Compare source metadata on this page
WorkPublishedSource identifierSource
Projected sea-level rise and high tide flooding at Big Cypress National Preserve, Florida2024-08-05fs20243022usgs
Assessment of long-term changes in surface-water extent within Klamath Marsh, south-central Oregon, 1985–20212024-07-30sir20245033usgs
3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope2024-07-2970258229usgs
Collision structures of the Prince William terrane and Chugach terrane docking along the Shumagin and Unimak convergent margins, Alaska, USA2024-07-2570256998usgs
Interactive effects of salinity and hydrology on radial growth of bald cypress (Taxodium distichum (L.) Rich.) in coastal Louisiana, USA2024-07-1970256391usgs
Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?2024-07-1670263545usgs
Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall2024-07-1570256026usgs
Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa2024-07-12sir20245059usgs
On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem2024-07-1070256002usgs
Characterization of the water resources of the Pamunkey River watershed in Virginia—A review of water science, management, and traditional ecological knowledge2024-07-08sir20245024usgs
Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States2024-07-0470255845usgs
Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California2024-07-0370265851usgs
A comparative analysis of OpenET for evaluating evapotranspiration in California almond orchards2024-07-0370269686usgs
New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States2024-07-0170259712usgs
The 17 January 1994 Northridge, California, earthquake: A retrospective analysis2024-07-0170263614usgs
Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network2024-06-25sir20245041usgs
Projected sea-level rise and high tide flooding at De Soto National Memorial, Florida2024-06-21fs20243016usgs
Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California2024-06-21sir20245029usgs

These are bibliographic comparisons, not experimental rankings. Follow the original document for methods and conditions.