USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science and Technology”

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,659 records · Page 85Linked to original sources
USGS2024

Automated, near real-time ground-motion processing at the U.S. Geological Survey

We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.

Seismological Research Letters↗
USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗
USGS2025

Drone-based radiometric surveys provide high-resolution mine waste characterization

Airborne radiometric surveys use passive geophysical techniques to characterize geochemical variations at or near earth’s surface. These methods have been used for a variety of mapping applications, including mineral resource evaluation. However, detailed characterization of smaller geologic targets, including mine waste features, requires flying at lower altitudes and with tighter line spacing than is feasible with traditional aircraft. Here, a small uncrewed aircraft system (sUAS) equipped with a radiometric sensor was used to acquire high-resolution gamma-spectrometry over small mine waste features and a low-grade stockpile in southwestern New Mexico. The sUAS radiometric system mapped local variability within each survey area and revealed ~2–10 m wide zones where radioelements K, Th, and U may be elevated 2–10× the surrounding material. Additionally, the sUAS radiometric data revealed radioelement variability across survey sites, which correlated reasonably well with variability seen in geochemical samples at each survey site, even though samples collected from individual sites showed high local variability. The sUAS data characterized local heterogeneity within mine waste and other small geologic targets at scales of a few meters to tens of meters, which is not possible with traditional crewed aircraft, and with continuity of coverage that is not possible with ground surveys, thus filling a key gap in geophysical survey spatial resolution.

New Mexico↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2025

Dating a medieval tsunami with uranium-series techniques on Caribbean corals

Uranium-series dates from coral boulders constrain the timing of a medieval tsunami from the Puerto Rico Trench. Previously reported evidence for this tsunami includes hundreds of coral boulders that came to rest hundreds of meters inland on Anegada, British Virgin Islands. New U-series dates on these coral boulders provide limiting dates for the tsunami. The narrowest limits were by dating interior bands of a coral that retains the hemispherical form of a living coral colony, and which include adjustments for the number of annual density band couplets between the dated samples and the boulder exteriors. By those limits, the tsunami dates between 1381 and 1391 CE, and likely occurred during summer or fall. The tsunami is important as the only reported sign that the eastern Puerto Rico Trench has produced a great earthquake. The dating may aid in defining the earthquake source and in communicating tsunami hazards.

British Virgin Islands, Puerto Rico↗
USGS1999

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota↗
USGS2026

Tracking baseflow supply dynamics using SWOT data from small groundwater-dominated lakes

In situ surface-water monitoring strategies are biased towards larger perennial streams and lakes and are generally not designed to track mechanisms of baseflow supply contributed by the dynamic storage of aquifers. Additionally, small (< 1 km 2 ) groundwater-influenced lakes and wetlands globally have little in situ monitoring infrastructure. We explored the utility of remotely sensed Surface Water Ocean Topography Satellite (SWOT) data, collected from 2023 onward, to characterise the seasonal and multi-year water-level trends of groundwater flow-through kettle lakes distributed across the permeable sediments of eastern Massachusetts, USA. This analysis indicated that water levels for kettle lakes with areas down to approximately 0.05 km 2 are resolvable in the study area. Our examination of 17 kettle lakes found that SWOT water-surface elevation data closely tracked groundwater levels in adjacent monitoring wells where available, including the timing of seasonal patterns (highest levels generally in late spring), although there was some variation between years and there was a substantial lag in the timing of high water levels for a lake located downgradient from a 30-m-thick vadose zone. Furthermore, SWOT-observed water-level increases in kettle lakes tracked with baseflow increases in two adjacent groundwater-dominated streams, as would be expected from increased hydraulic gradients. Unlike spectral remote sensing, SWOT data are generally not affected by cloud cover, resulting in a potential for groundwater-dominated lakes to be sentinels of dynamic storage patterns, including identification of baseflow drought lags, which are currently ill-defined hydrological processes. SWOT monitoring of groundwater-influenced surface waters shows potential for augmenting existing monitoring wells and streamgages as continuous monitors of groundwater levels and baseflow supply in permeable terrain.

Massachusetts↗
USGS2025

Accurate simulation of flow through dipping aquifers with MODFLOW 6 using enhanced cell connectivity

In simulations of groundwater flow through dipping aquifers, layers of model cells are often “deformed” to follow the top and bottom elevations of the aquifers. When this approach is used in MODFLOW, adjacent cells within the same model layer are vertically offset from one another, and the standard conductance-based (two-point) formulation for flow between cells does not rigorously account for these offsets. The XT3D multi-point flow formulation in MODFLOW 6 is designed to account for geometric irregularities in the grid, including vertical offsets, and to provide accurate results for both isotropic and anisotropic groundwater flow. A recent study evaluated the performance of the standard formulation and XT3D using a simple, synthetic benchmark model of a steeply dipping aquifer. Although XT3D generally improved the accuracy of flow simulations relative to the standard formulation as expected, neither formulation produced accurate flows in cases that involved large vertical offsets. In this paper, we explain that the inability of XT3D to produce accurate flows in the steeply dipping aquifer benchmark was not due to an inherent limitation of the flow formulation, but rather to the limited cell connectivity inherent in the most commonly used discretization packages in MODFLOW 6. Furthermore, we demonstrate that XT3D is able to produce the expected accuracy when adequate cell connectivity is introduced using MODFLOW's unstructured grid type and the aquifer is discretized vertically using at least two model layers.

Groundwater↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2026

From model to action: Identifying the gaps in coastal marsh models for decision making

Coastal salt marsh modeling projects are designed to answer questions about where salt marshes may migrate, how they can migrate, and how they will be impacted by changing water levels. While this information is vital for climate-resilient marsh conservation, translating these findings into actionable results remains a challenge. The authors examined the actionability of marsh modeling efforts in the United States by analyzing data collected from semi-structured interviews of marsh model users ( n = 24) across two ongoing projects along the U.S. East Coast and the Gulf of America (Gulf of Mexico). By qualitatively analyzing the interview data, the authors found that the tasks of practitioners who use coastal marsh model outputs fall into three major themes: (1) marsh restoration, (2) planning with uncertainty, and (3) conserving habitat for marsh-reliant species. For each of these themes, the authors identify unmet needs, including high spatial resolution information for local planning, accessible descriptions of uncertainty to increase user confidence, and the incorporation of human dimensions data (e.g., human alterations to the landscape such as impoundments and culverts) for a comprehensive understanding of the coastal salt marsh. Marsh modeling projects should strategize how to fulfil these unmet user needs so that marsh modeling outputs can better support decisions related to marsh conservation and restoration.

Gulf of America, United States East Coast↗
USGS2024

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

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper↗
USGS2025

Groundwater surveillance of swine pathogens from private wells supplying swine farms in Iowa

Biosecurity practices are essential for maintaining pig health and productivity. Despite these measures, pathogen spread still occurs. Water is one of the largest daily inputs on swine farms by volume and is not routinely tested or disinfected before it is consumed by the animals [1-3], making it a poorly understood biosecurity risk. Groundwater from privately-owned wells is a common water source for swine farms. Pathogens in the landscape, such as bacteria, viruses, and protozoa, can reach groundwater more rapidly through soil macropores, maintaining viability and facilitating transmission of pathogens into aquifers [3-13].

Iowa↗
USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗
USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗
USGS2025

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2024

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗