USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Monitoring and Assessment”

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,662 records · Page 21Linked to original sources
USGS2000

Assessment of potential for natural attenuation of chlorinated ethenes and ethanes in ground water at a petrochemical reclamation site, Harris County, Texas

Redox conditions in the Numerous Sand Channels Zone beneath a petrochemical reclamation site in Harris County, Texas, range from sulfate reducing to methanogenic as indicated by the presence of methane in ground water and the range of molecular hydrogen concentrations. Assessment of the potential for reductive dechlorination using BIOCHLOR as a screening tool indicated conditions favoring anaerobic degradation of chlorinated organic compounds in the Numerous Sand Channels Zone. Evidence supporting reductive dechlorination includes apparently biogenic cis -1,2-dichloroethene; an increased ratio of 1,2-dichloroethane to 1,1,2-trichloroethane downgradient from the assumed contaminant source area; ethene and methane concentrations greater than background concentrations within the area of the contaminant plume; and a positive correlation of the ratio of ethene to vinyl chloride as a function of methane concentrations. The body of evidence presented in this report argues for hydrogenolysis of trichloroethene to cis -1,2-dichloroethene; of 1,1,2-trichloroethane to 1,2-dichloroethane; and of vinyl chloride to ethene within the Numerous Sand Channels Zone. Simulations using BIOCHLOR yielded apparent first-order decay constants for reductive dechlorination in the sequence tetrachloroethene trichloroethene cis -1,2-dichloroethene vinyl chloride ethene within the range of literature values reported for each compound and apparent first-order decay constants for reductive dechlorination in the sequence 1,1,2-trichloroethane 1,2-dichloroethane slightly greater than literature values reported for each compound along the upgradient segment of a simulated ground-water flowpath. Except for vinyl chloride, apparent rates of reductive dechlorination for all simulated species show a marked decrease along the downgradient segment of the simulated ground-water flowpath. Evidence for reductive dechlorination of chlorinated ethenes within the Numerous Sand Channels Zone indicates potential for natural attenuation of chlorinated ethenes. Reductive dechlorination of chlorinated ethanes apparently occurs to a lesser extent, indicating relatively less potential for natural attenuation of chlorinated ethanes. Additional data are needed on the concentrations and distribution of chlorinated ethenes and ethanes in individual fine sand intervals of the Numerous Sand Channels Zone. This information, combined with lower minimum reporting levels for future chloroethane analyses, might enable a more complete and quantitative assessment of the potential for natural attenuation at the site.

Texas↗
USGS2025

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

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

Oregon↗
USGS2026

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗
USGS2025

Meet the people where they are: Assessing user needs for aftershock forecast products in El Salvador, Mexico and the United States

Aftershock forecasts can help communities reduce their seismic risk by conveying how many aftershocks can be expected following a large earthquake, and how the expected number of aftershocks and their corresponding ground shaking evolves over time and space. Prior work finds that graphical forecast products may communicate such information better than only text or numbers. To identify which visual products can serve multiple user groups, we held workshops with members of several professions, including emergency managers, engineers, critical infrastructure operators, public health specialists, science communicators, and more. We conducted these workshops in El Salvador, Mexico and the United States to understand which forecast products may be effective across different countries. In these workshops, users performed small-group activities to elicit the types of aftershock forecast information that would support decisions in their respective roles and how this information would optimally be displayed. Maps of shaking hazards were frequently requested across all professions and countries, even for dissimilar forecast uses. The design of these maps, and other product needs, were differentiated by profession, country and other factors. Other forecast products, including those showing forecasts about the magnitudes and time periods of future aftershocks, also served a variety of users, but for different types of decisions. We found a greater variation in user needs by profession than by country, and that user needs also vary with time, communication channel and other contextual factors. We discuss practical implications for user-centered visual communication of operational aftershock forecasts.

International Journal for Disaster Risk Reduction↗
USGS2025

Migration of seismicity from the mantle to the upper crust beneath Harrat Lunayyir volcanic field, Saudi Arabia

Harrat Lunayyir is a volcanic field in Saudi Arabia that experienced a Mw~5.4 earthquake driven by an upper-crustal dike intrusion in May 2009. This volcanic field has exhibited numerous forms of volcanic seismicity both prior to and since the 2009 dike intrusion. Significantly, earthquakes within the lithospheric mantle and, rarely, the lower crust are present in the two-decade long seismicity catalog of Harrat Lunayyir. Here we analyze 24 years of volcanic seismicity at Harrat Lunayyir from 1998 to 2022. We find that: 1) precursory seismicity began at least eight years prior to the 2009 event, with a particularly notable seismic episode one year prior; 2) lithospheric mantle seismicity is highly localized in space and in time, largely occurring in discrete sequences lasting on the order of a few hours to a few days; 3) one seismic sequence clearly migrates upward from the lithospheric mantle to the upper crust, including seismicity within the nominally ductile lower crust; 4) crustal seismicity has been slowly declining over time; and 5) lithospheric-mantle seismicity does not show any apparent decline with time. From these observations we infer that the seismicity is driven by magmatic fluids or volatiles, and seismic monitoring of this volcanic field should continue into the future.

Harrat Lunayyir volcanic field↗
USGS2024

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances↗
USGS2025

Ordovician stratigraphy, structure, and karst of the Falling Spring Valley, Alleghany County, Virginia, USA

This one-day trip highlights new findings on a preliminary bedrock geologic map that shows results from ongoing geologic mapping in the Falling Spring Valley of Alleghany County, Virginia, USA, which is the southern end of the larger Warm Springs Valley, an elongated anticlinal valley rimmed by Ordovician and Silurian siliciclastic rocks, and which is famous for its thermal springs. This mapping includes stratigraphic, structural, and karst field and lab research focused on the Ordovician strata exposed in the area, the oldest of which is the dolomitic upper part of the Beekmantown Formation (Lower Ordovician, Darriwilian), and the youngest of which is the Juniata Formation (Upper Ordovician, Katian), a sequence of siliciclastic redbeds. Warm Springs Valley is the location of the only known caves in the eastern United States—three at present—with thermal waters flowing in some of their passages. Stops on the trip will highlight key details from mapping efforts, primarily within the structurally deformed Ordovician carbonate sequence that is exposed in the core and limbs of the anticline, as well as the associated karst features that are developed in those carbonate rocks, including results of recent dye traces and water temperature monitoring that have improved our understanding of the karst hydrogeologic systems developed in these strata.

Virginia↗
USGS2025

Expression of corticoid-regulatory genes in the gills of Atlantic salmon (Salmo salar) parr and smolt and during salinity acclimation

In teleost fishes, cortisol is the major corticoid and has both glucocorticoid and mineralocorticoid actions. However, how fish tissues discriminate between these distinct corticosteroid actions is unclear. In mammals, the major factors responsible for intracellular corticosteroid regulation are glucocorticoid receptors ( gr s) and the mineralocorticoid receptor ( mr ), but their role in osmoregulation of fish is unclear. 11β-hydroxysteroid dehydrogenases ( hsd11b s) control the levels of intracellular corticosteroids by converting from bioactive forms to inert forms. To investigate how Atlantic salmon ( Salmo salar ) respond to cortisol in different physiological or environmental conditions, we performed comparisons of parr and smolt, and osmotic challenge experiments to examine the physiological responses and gill transcript levels of genes underlying cortisol-signalling, including gr1 , gr2 , mr , hsd11b2 and hsd11b3 . Because cortisol may interact with growth hormone and prolactin during salinity changes, transcript levels encoding growth hormone receptors ( ghr1 , ghr2 ) and the prolactin receptor ( prlr ) were also examined. Hsd11b2 transcript levels in seawater-acclimated fish were consistently lower compared to fish acclimated to fresh water. After transfer to seawater, prlr transcript levels in fish significantly decreased and transcript levels of ghr1 , ghr2 and hsd11b3 showed no change or were slightly higher than those of freshwater control groups. Gr1 , gr2 and mr transcript levels were slightly but consistently higher in fish acclimated to fresh water relative to seawater. Our results indicate that changes in corticosteroid receptor and hsd11b2 transcript levels in the gills may be important mechanisms that regulate corticoid signals to achieve ion homeostasis in Atlantic salmon.

Journal of Fish Biology↗
USGS2024

Models no not provide proof: An example of model ambiguity and application of isotopic data in a mine pit lake

Geochemical and hydrologic models of pit lakes are commonly used in environmental regulatory decisions to predict future water quality and hydrologic conditions and to understand existing pit lakes. Models may be used to quantify sulfide oxidation, predict thermal/chemical stratification and mixing, and better understand connections between pit lakes and aquifers. One concern related to the hydrologic character of pit lakes is if they are terminal (a groundwater sink with no outflow) or flowthrough (both receiving groundwater inflow and discharging to groundwater). This question was pertinent to the Liberty pit lake, a small acidic pit lake formed in a former Cu deposit in south-central Nevada where potentiometric and geochemical data potentially indicate pit-lake outflow. Potential discharge to groundwater from the pit lake was evaluated using a water-balance model, but uncertainty in hydraulic parameters led to ambiguity in the hydrologic character. Stable isotopes of water were then sampled from the pit lake and adjacent groundwater wells, which unambiguously indicated the lack of an evaporative signature in downgradient groundwater because the groundwater did not plot on a hypothetical mixing line between evaporated pit lake water and observed meteoric recharge. This methodology provided a more effective and more data-driven approach for understanding pit-lake hydrology. Although predictive models are required to quantify reasonable bounds on future conditions, many models contain substantial uncertainty and are not well suited in some environments. Datasets that provide more clear lines of evidence could be collected from existing pit lakes whenever possible to inform water-rock interaction, limnological behavior, and connectivity to adjacent groundwater.

Nevada↗
USGS2025

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv↗
USGS2025

Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?

In a previous study, bighorn sheep ( Ovis canadensis ) were periodically photographed by trail cameras when they visited desert tortoise ( Gopherus agassizii ) burrows in the San Bernardino Mountains near Palm Springs, California, USA. The authors suggested that bighorn sheep may utilize the excavated spoil material from burrows as a mineral lick. To test that supposition, we collected soil samples from the spoil piles and undisturbed surface soil up to 1 m away at two desert tortoise study sites with bighorn sheep at the aforementioned site and another nearby site. We hypothesized that sodium concentrations would be greater in spoil than in nearby surface soil samples. Common major cations were quantified in soil samples using inductively coupled plasma-mass spectrometry following microwave-assisted digestion to solubilize the fraction of elements that could become biologically or environmentally available. Mean sodium concentrations were greater in burrow spoil samples than in undisturbed soil samples at both sites, likely due to aqueous leaching of sodium from surface to subsurface soil horizons. Although we did not collect evidence of geophagy by bighorn sheep, the behavior is well known. However, trail camera photographs of bighorn sheep with their head next to the entrance of tortoise burrows suggests that they may periodically utilize sodium and other elements brought to the surface by desert tortoises that are otherwise unavailable to bighorn sheep on the surface.

California↗
USGS1996

Pesticides in surface water from three agricultural basins in south-central Georgia, 1993-97

Twenty-two of 43 pesticides analyzed were detected in 128 water samples collected from the Tucsawhatchee Creek, the Little River, and the Withlacoochee River. These streams drain agricultural basins in south-central Georgia and were sampled from March 1993 through June 1995. Herbicides were detected more frequently than insecticides. The most frequently detected herbicides were atrazine and metolachlor and the most frequently detected insecticide was carbaryl. Pesticide concentrations in the three streams were low and did not exceed U.S. Environmental Protection Agency drinking water standards. The maximum pesticide concentration was 2.6 mg/L (micrograms per liter) for propargite, a miticide detected in only one sample. The maximum concentrations of the remaining 21 pesticides were less than 0.25 mg/L. The median concentrations were equal to the method detection limit for all pesticides except atrazine (0.008 mg/L) and metolachlor (0.012 mg/L). The ratio of herbicide detections to nondetections was largest in the planting season, smaller in the harvest season and smallest in the fallow season for the three basins. The same pattern existed for the insecticide ratios in the Little River and the Withlacoochee River. Pairwise correlations between concentrations of atrazine and metolachlor, and four parameters (discharge, and concentrations of dissolved organic carbon, suspended organic carbon, and suspended sediment) were evaluated for each stream. The strongest correlations existed between metolachlor and mean daily discharge, and metolachlor and sediment in the Withlacoochee River. The only significant correlation for the Little River was between atrazine and suspended sediment.

Georgia↗
USGS2024

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↗
USGS2024

A global assessment of SAOCOM-1 L-band stripmap data for InSAR characterization of volcanic, tectonic, cryospheric, and anthropogenic deformation

SAOCOM-1 is an L-band (23.5 cm) synthetic aperture radar (SAR) constellation made up of two satellites launched in 2018 and 2020 by Comisión Nacional de Actividades Espaciales (CONAE, Argentina). In this contribution, we present a global summary of interferometric SAR (InSAR) observations of ground deformation with SAOCOM-1 stripmap data for tracking volcanic, tectonic, glacier, and anthropogenic deformation. These examples include: 1) episodes of unrest at volcanoes in the Aleutian Islands, Southern Andes, and Italy, with line-of-sight (LOS) deformation from 4 cm/yr in InSAR time series to ~70 cm in interferograms; 2) dike intrusions in Hawai’i; 3) earthquakes in the Andean fold and thrust belt and the East Anatolian fault; 4) ice flow of the Southern Patagonia icefield; and 5) subsidence due to lithium brine extraction in the Salar de Atacama basin (northern Chile). Comparisons between SAOCOM-1, ALOS-2 SM3, Sentinel-1, and TerraSAR-X/ TanDEM-X/PAZ (TSX/TDX/PAZ) mean velocities from InSAR time series show a 1:1 ± 3% correlation in the LOS velocity, which highlights the high accuracy of SAOCOM-1 data. The minimum deformation that we measured in individual interferograms is 4 ± 0.6 cm. One limitation of SAOCOM-1 is the lack of a global acquisition program, which reduces its global and broader applications. Considering the repeat periods, background observation program, and lack of a controlled orbital tube, the best suited targets for SAOCOM-1 InSAR are two. First, volcanoes that deform with secular rates located in vegetated regions in mid- and high-latitudes, and/or that undergo transient episodes of fast deformation in which C-band coherence is lost quickly. Second, glaciers where coherence can be sustained during the repeat period of eight days.

IEEE Transactions on Geoscience and Remote Sensing↗
USGS2022

South Carolina and Landsat

South Carolina, the eighth State admitted to the union, transcends its size with its deep, rich history; striking beauty; vast natural resources; and extensive cultural diversity. Home to part of the Blue Ridge Mountains of the Central Appalachians, the Upstate is graced with more than 100 waterfalls, while the Lowcountry borders the Atlantic Ocean with 187 miles of coastline and 35 barrier islands. Forests cover two-thirds of the State, and forestry and agriculture together, as agribusiness, make up South Carolina’s leading industry. Two historic crops—cotton and tobacco—still rank in the top 10 commodities, though corn and soybeans now rank higher. Poultry, cattle, peanuts, and flowers also make the list. South Carolina’s population totals more than five million. Other residents include a variety of wildlife, bird, reptile, and fish species, including Ursus americanus (black bears), Alligator mississippiensis (American alligators), and Tursiops truncatus (bottlenose dolphins). More than 100 tree species also reside in South Carolina, which pays homage to one with its “The Palmetto State” nickname. South Carolina’s subtropical climate, long coastline, and lower elevations make it highly susceptible to tornado and hurricane activity and coastal flooding. Projected sea-level rise is a growing concern. A view from space can help monitor and manage natural resources on the land and in rivers, marshes, and the coast. Landsat reveals not just what an area looks like now, but also insights from decades ago.

South Carolina↗
USGS2024

Participatory engagement to reduce communication gaps

Underserved communities, especially those in coastal areas in Puerto Rico, face significant threats from natural hazards such as hurricanes and rising sea levels. Limited funding hinders the investment in costly mitigation measures, increasing exposure to natural disasters. Providing coastal resources and data products through effective communication mechanisms is fundamental to improving the well-being of these underserved coastal communities. The overall objectives of the pilot effort to engage and connect with underserved coastal communities in Puerto Rico were the following: (1) compile a comprehensive database of the projects and resources relevant to natural hazards in Puerto Rico; (2) foster connections with Puerto Rican interested parties to better understand their priorities regarding coastal hazards and provide them with pertinent U.S. Geological Survey (USGS) resources; and (3) identify knowledge gaps to guide future USGS projects in Puerto Rico. Here we outline our participatory engagement framework and process, along with two specific resources developed with the information collected from this effort. These resources are available in English and Spanish and consist of user-friendly, non-technical information products. Among them are: (1) a website where users can learn about USGS research on landslides, hurricanes, earthquakes, water resources, coastal hazards, tsunamis, and ecosystem hazards and environmental contaminants, and (2) a geonarrative highlighting shoreline changes in Puerto Rico with sections on historical shoreline trends, hurricane impacts, and potential solutions that could help protect communities and mitigate coastal hazards. Continuing participatory engagement in future projects could enhance the accessibility and usability of natural hazards resources within the community.

San Juan↗
USGS2024

Comparison of water quality in shallow groundwater near agricultural areas in the Delaware Coastal Plain, 2014 and 2019

The State of Delaware has encouraged agricultural conservation practices to improve nutrient uptake by crops and mitigate nutrient transport to groundwater in the surficial aquifer. To study recent changes in groundwater quality, the U.S. Geological Survey and the Delaware Department of Agriculture (DDA) developed a network of shallow wells near agricultural areas throughout the Delaware Coastal Plain. This network was designed to characterize water quality related to agricultural practices and to detect any recent changes in shallow groundwater quality, in particular groundwater nitrate concentrations. The shallow well network was first sampled in 2014 and resampled in 2019. In 2019, field parameters (including dissolved oxygen, pH, specific conductance, and temperature), major ions, nutrients, stable isotopes of water, and isotopes of nitrate were measured in groundwater samples collected between October and December. Wells were organized into three groups based on their geochemical characteristics measured in 2014: the Agricultural, Urban, and Mixed Groups. Results from the 2019 sampling showed little change in water quality from the 2014 sampling. Land-use factors continued to be the driving influence between groups. Groundwater moves slowly and changes in groundwater quality are likely to respond slowly to changes in conservation practices. Continued sampling of both groundwater quality in this network and monitoring land management practices can help detect groundwater quality trends in the future.

Delaware↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗