USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Analytical and Applied Pyrolysis”

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,688 records · Page 67Linked to original sources
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2025

A crosswalk of the 2015 World Terrestrial Ecosystems to the International Union for the Conservation of Nature Global Ecosystem Typology Framework

To support ecosystem mapping and accounting applications, we aligned the 2015 U.S. Geological Survey/Esri/The Nature Conservancy-World Terrestrial Ecosystems (WTEs) with the International Union for Conservation of Nature Global Ecosystem Typology (GET) framework. This process, known as “crosswalking,” enabled the development of a global map of GET level 3 Ecosystem Functional Groups (EFGs) at a 250-meter spatial resolution. Crosswalking involved manually assigning 1,781 biogeographically stratified WTEs to their most probable EFG based on similarities in climate, terrain, vegetation, and geographic distribution. We compared attributes of the WTE dataset with summary characteristics of the EFGs. The resulting crosswalked global map of International Union for Conservation of Nature GET ecosystems is intended to be useful for standardizing ecosystem classification and reporting under frameworks such as the Kunming-Montreal Global Biodiversity Framework and the United Nations System of Environmental-Economic Accounting. We discuss key challenges in reconciling non-identical classifications, such as many-to-one relationships and variation in data quality.

Open-File Report↗
USGS2024

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

Albuquerque Seismological Laboratory strategic vision

This circular presents a strategic outlook for the U.S. Geological Survey (USGS) Albuquerque Seismological Laboratory (ASL) for the next 10 years (2026–36). The ASL is a USGS field office in the Geological Hazards Science Center that operates portions of the Advanced National Seismic System and the Global Seismographic Network and focuses on fundamental research for instrumentation testing and data quality. The strategic outlook is categorized into two types of tasks: “Foundational Tasks” and “Aspirational Tasks.” Foundational Tasks are those that maintain the laboratory’s basic operations and services, including regional and global seismic monitoring, improving data quality, and providing instrument testing and support. A suite of Aspirational Tasks is also articulated; these can be considered priority targets of ASL that could improve ASL’s seismic monitoring capabilities and mission. Such tasks include improvements to remote stations, testing capabilities of nonseismic geophysical instruments, detection threshold monitoring, rapid aftershock deployments, and expanding seismic monitoring networks. This report was written with input from the USGS Geological Hazards Science Center, the USGS Earthquake Hazards Program (EHP), and colleagues with an interest in the work done by the ASL. Although the details of these tasks may change, this document can provide guidance on the overarching tasks at the ASL from 2026 to 2036 and an overview of the various components of the ASL and how they fit into the EHP and the Global Seismographic Network Program.

Circular↗
USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗
USGS2025

Engineering and risk research—A strategic vision for extending USGS earthquake science to risk assessment and mitigation

As a part of the U.S. Geological Survey (USGS) Natural Hazards Mission Area, the Engineering and Risk Project (ERP) of the Geologic Hazards Science Center conducts civil engineering research to extend natural hazards science to risk assessments. The ERP includes predominantly, but not exclusively, staff funded by the USGS Earthquake Hazards Program. A major focus of the ERP is to deliver actionable information to decisionmakers that helps promote the safety, security, and economic well-being of the Nation. The report discusses briefly the ERP’s evolution and scope, its organization structure, describes key products and deliverables, their efforts to engage with external partners, and its strategy and vision in conducting range of engineering and risk research. The research topics outlined in this report, as developed by the ERP, may be used to facilitate project planning activities during the next 10 years and serve as a longer-term planning resource for the USGS Earthquake Hazards Program and Geologic Hazards Science Center.

Circular↗
USGS2025

The Hardscrabble Creek complex: A newly discovered, mostly buried, Mesoproterozoic mafic-ultramafic pluton in the Wet Mountains, Colorado, USA

The origin of prolific ca. 1.4 Ga ferroan magmatism between the southwestern USA and eastern Canada is enigmatic and has been explained by various models, including extensional, mantle plume, and convergent plate-margin models. Rare mafic plutons are associated with the ferroan plutons, which may help constrain their mantle source and tectonic setting. In the southwestern USA, only two such mafic plutons are known to exist. We present the first evidence for a third, mostly buried, potentially layered, mafic-ultramafic Mesoproterozoic pluton, informally referred to as the Hardscrabble Creek complex, in the central Wet Mountains of Colorado, USA. Recent geophysical data show an elliptical magnetic and gravity high spatially coincident with local gabbroic outcrops. New field and petrographic analyses of these exposed rocks reveal that they consist of ultramafic to mafic cumulates, including orthopyroxenite, olivine norite, norite, and anorthosite. High-precision U-Pb dating of zircon from orthopyroxenite and norite yield weighted mean 206 Pb/ 238 U dates of 1352.36 ± 1.60 Ma and 1352.37 ± 1.71 Ma, respectively. These dates indicate that the complex formed over a narrow timeframe, after the adjacent 1362 ± 7 Ma ferroan San Isabel Granite, and during the waning stages of the regional ca. 1.4 Ga ferroan magmatism. Whole-rock geochemistry and Nd-Sr-Pb isotope compositions of samples from the Hardscrabble Creek complex are similar to those of the San Isabel Granite, suggesting that they were derived from the same or a similar mantle source. The mineral chemistry of the samples is comparable to Proterozoic massif-type anorthosites and related mafic intrusions, indicating that the Hardscrabble Creek complex and San Isabel Granite together represent a rare anorthosite-mangerite-charnockite-granite (AMCG) suite in the southwestern USA. The Hardscrabble Creek complex is unique because it formed ~80 m.y. after the other few mafic plutons in the southwestern USA, and it contains an ultramafic section that is absent from these plutons and rare to the AMCG suite in general. A combination of arc-like whole-rock geochemistry, chondrite uniform reservoir-like Nd-Sr-Pb isotopes, and ocean island basalt (OIB)-like zircon trace element chemistry suggests that the complex was derived from a partial melt of OIB-like mantle and interacted with metasomatically enriched lithospheric mantle. The enriched lithospheric mantle signature, combined with the long ~160 m.y. duration of magmatism in the region, is consistent with a period of protracted convergent tectonism.

Colorado↗
USGS1976

Non-radiometric borehole geophysical detection of geochemical halos surrounding sedimentary uranium deposits

Roll-type uranium deposits are formed by the concentration of uranium by ground water in geochemical cells. Non-uranium minerals having different solubilities may be deposited ahead of or behind the uranium minerals, forming halos that surround the ore. In addition, oxidizing and reducing environmental conditions may cause zones of mineral alteration to develop beyond the limits of the uranium deposit. Certain physical property anomalies that are commonly associated with halos can be detected by relatively fast and inexpensive borehole geophysical measurements made either in individual holes, or between two adjacent holes. Borehole measurements that have been found to be useful include electrical resistivity, induced polarization, and magnetic susceptibility. Electrical resistivity is increased by the presence of calcite and other cementing minerals that sometimes create permeability barriers in the neighborhood of uranium deposits. Induced polarization (IP) response is increased by sulfide and clay minerals that are commonly found in anomalous concentrations near roll-type deposits. Magnetic susceptibility is usually decreased by the oxidation of magnetite to hematite or limonite in the zone of chemical alteration that is left as a trail behind roll fronts. Borehole measurements of electrical resistivity, induced polarization and magnetic susceptibility were made in the vicinity of a uranium roll-type deposit in south Texas. Results indicate that mineral halos can be detected by borehole measurements made. in wide-spaced drill holes, and that the total amount of drilling needed to find a deposit can be reduced substantially by this exploration approach.

Open-File Report↗
USGS2025

Mineralization processes at Escanaba Trough

Escanaba Trough off northern California is the southernmost and only sediment covered segment of the Gorda Ridge. Hydrothermal mineralization was first discovered in Escanaba Trough in 1985, and subsequent investigations included the delineation of the NESCA (northern Escanaba) and SESCA (southern Escanaba) sections of Escanaba Trough, as well as drilling by Ocean Drilling Program (ODP) Leg 169 in 1996. Preliminary research conducted in the 1980s distinguished two classes of sulfide mineral assemblages within the Escanaba Trough: one dominated by pyrrhotite (Fe 1-x S) and one containing ‘mixed sulfides’ including chalcopyrite (CuFeS 2 ) and sphalerite (ZnS). In 2022, we conducted a research expedition to the Escanaba Trough to more thoroughly characterize mineralization processes and consider preservation of precipitated minerals throughout hydrothermally active and inactive sections of Escanaba Trough, with the specific objective of characterizing the host minerals, and processes of enrichment for specific elements which are currently listed as critical elements and can also be toxic under certain conditions (Sb, As, Ba, Co, Ga, Ge, Te). The following will include initial results from that work.

California↗
USGS1976

Distribution functions for statistics derived from bivariate normal and bivariate two-parameter log-normal populations

The distribution functions for statistics that may be used to assess the significance of differences between sample means, standard deviations, coefficients of skewness, and coefficients of variation are obtained by Monte Carlo experiments. The distributions are expressed as functions of sample size, cross correlation, and skewness. In general, the distributions are more sensitive to cross correlation than to skewness. As sample size increases, however, the distributions tend to become more sensitive to skewness.

Open-File Report↗
USGS2025

A global database of soil microbial phospholipid fatty acids and enzyme activities

Soil microbes drive ecosystem function and play a critical role in how ecosystems respond to global change. Research surrounding soil microbial communities has rapidly increased in recent decades, and substantial data relating to phospholipid fatty acids (PLFAs) and potential enzyme activity have been collected and analysed. However, studies have mostly been restricted to local and regional scales, and their accuracy and usefulness are limited by the extent of accessible data. Here we aim to improve data availability by collating a global database of soil PLFA and potential enzyme activity measurements from 12,258 georeferenced samples located across all continents, 5.1% of which have not previously been published. The database contains data relating to 113 PLFAs and 26 enzyme activities, and includes metadata such as sampling date, sample depth, and soil pH, total carbon, and total nitrogen. This database will help researchers in conducting both global- and local-scale studies to better understand soil microbial biomass and function.

Scientific Data↗
USGS2021

Factors influencing distribution of Coccidioides immitis in soil, Washington State, 2016

Coccidioides immitis and Coccidioides posadasii are causative agents of Valley fever, a serious fungal disease endemic to regions with hot, arid climate in the United States, Mexico, and Central and South America. The environmental niche of Coccidioide s spp. is not well defined, and it remains unknown whether these fungi are primarily associated with rodents or grow as saprotrophs in soil. To better understand the environmental reservoir of these pathogens, we used a systematic soil sampling approach, quantitative PCR (qPCR), culture, whole-genome sequencing, and soil chemical analysis to identify factors associated with the presence of C. immitis at a known colonization site in Washington State linked to a human case in 2010. We found that the same strain colonized an area of over 46,000 m 2 and persisted in soil for over 6 years. No association with rodent burrows was observed, as C. immitis DNA was as likely to be detected inside rodent holes as it was in the surrounding soil. In addition, the presence of C. immitis DNA in soil was correlated with elevated levels of boron, calcium, magnesium, sodium, and silicon in soil leachates. We also observed differences in the microbial communities between C. immitis -positive and -negative soils. Our artificial soil inoculation experiments demonstrated that C. immitis can use soil as a sole source of nutrients. Taken together, these results suggest that soil parameters need to be considered when modeling the distribution of this fungus in the environment.

Washington↗
USGS1997

Hydrologic and water-quality conditions in the Horse Creek basin, west-central Florida, October 1992-February 1995

A baseline study of the 241-square-mile Horse Creek basin was undertaken from October 1992 to February 1995 to assess the hydrologic and water-quality conditions of one of the last remaining undeveloped basins in west-central Florida. During the period of the study, much of the basin remained in a natural state, except for limited areas of cattle and citrus production and phosphate mining. Rainfall in 1993 and 1994 in the Horse Creek basin was 8 and 31 percent, respectively, above the 30-year long-term average. The lowest and highest maximum instantaneous peak discharge of the six daily discharge stations occurred at the Buzzard Roost Branch and the Horse Creek near Arcadia stations with 185 to 4,180 cubic feet per second, respectively. The Horse Creek near Arcadia station had the lowest number of no-flow days with zero days and the Brushy Creek station had the highest number with 113 days. During the study, the West Fork Horse Creek subbasin had the highest daily mean discharge per square mile with 30.6 cubic feet per second per square mile, and the largest runoff coefficient of 43.7 percent. The Buzzard Roost Branch subbasin had the lowest daily mean discharge per square mile with 5.05 cubic feet per second per square mile, and Brushy Creek and Brandy Branch shared the lowest runoff coefficient of 0.6 percent. Brandy Branch had the highest monthly mean runoff in both 1993 and 1994 with 11.48 and 19.28 inches, respectively. During the high-baseflow seepage run, seepage gains were 8.87 cubic feet per second along the 43-mile Horse Creek channel. However, during the low-baseflow seepage run, seepage losses were 0.88 cubic foot per second. Three methods were used to estimate average annual ground-water recharge in the Horse Creek basin: (1) well hydrograph, (2) chloride mass balance, and (3) streamflow hydrograph. Estimated average annual recharge using these three methods ranged from 3.6 to 8.7 inches. The high percentage of carbonate plus bicarbonate analyzed at the Carlton surficial aquifer well could indicate an upward ground-water flow from the underlying intermediate aquifer system. Based on constituent concentrations in water samples from the six daily discharge stations, concentrations generally are lower in the upper three subbasins, West Fork Horse Creek, Upper Horse Creek, and Brushy Creek than in the lower three subbasins. Typically, concentrations were highest for major ions at Buzzard Roost Branch and nutrients at Brushy Creek.

Florida↗
USGS2025

From hydrated silica to quartz: Potential hydrothermal precipitates found in Jezero crater, Mars

On Earth, silica-rich phases from opal to quartz are important indicators and tracers of geological processes. Hydrated silica, such as opal, is a particularly good matrix for the preservation of molecular and macroscopic biosignatures. Cherts, a type of silica-dominated rocks, provide a unique archive of ancient terrestrial life while quartz is the emblematic mineral of the Earth's continental crust. On Mars, hydrated silica has been detected in several locations based on remote sensing and rover-based studies. In the present article we report on the detection of cobbles made of hydrated silica (opal or chalcedony), as well as well-crystallized quartz. These detections were made with the SuperCam instrument onboard Perseverance (Mars 2020 mission), using a combination of LIBS, infrared and Raman spectroscopy. Quartz-dominated stones are detected unambiguously for the first time on the Martian surface, and based on grain size and crystallinity are proposed to be of hydrothermal origin. Although these rocks were all found as float, we propose that these detections are part of a common hydrothermal system, and represent different depths / temperatures of precipitation. This attests that hydrothermal processes were active in and around Jezero crater, possibly triggered by the Jezero crater-forming impact. These silica-rich rocks, in particular opaline silica, are very promising targets for sampling and return to Earth given their high biosignature preservation potential.

Earth and Planetary Science Letters↗
USGS1997

An optimized network for phosphorus load monitoring for Lake Okeechobee, Florida

Phosphorus load data were evaluated for Lake Okeechobee, Florida, for water years 1982 through 1991. Standard errors for load estimates were computed from available phosphorus concentration and daily discharge data. Components of error were associated with uncertainty in concentration and discharge data and were calculated for existing conditions and for 6 alternative load-monitoring scenarios for each of 48 distinct inflows. Benefit-cost ratios were computed for each alternative monitoring scenario at each site by dividing estimated reductions in load uncertainty by the 5-year average costs of each scenario in 1992 dollars. Absolute and marginal benefit-cost ratios were compared in an iterative optimization scheme to determine the most cost-effective combination of discharge and concentration monitoring scenarios for the lake. If the current (1992) discharge-monitoring network around the lake is maintained, the water-quality sampling at each inflow site twice each year is continued, and the nature of loading remains the same, the standard error of computed mean-annual load is estimated at about 98 metric tons per year compared to an absolute loading rate (inflows and outflows) of 530 metric tons per year. This produces a relative uncertainty of nearly 20 percent. The standard error in load can be reduced to about 20 metric tons per year (4 percent) by adopting an optimized set of monitoring alternatives at a cost of an additional $200,000 per year. The final optimized network prescribes changes to improve both concentration and discharge monitoring. These changes include the addition of intensive sampling with automatic samplers at 11 sites, the initiation of event-based sampling by observers at another 5 sites, the continuation of periodic sampling 12 times per year at 1 site, the installation of acoustic velocity meters to improve discharge gaging at 9 sites, and the improvement of a discharge rating at 1 site.

Florida↗
USGS2024

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska↗