USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Observer”

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.

At least 721 records · Page 40Linked to original sources
USGS1974

Apollo 11 voice transcript pertaining to the geology of the landing site

On July 20, 1969, America's Eagle touched down in southwestern Mare Tranquillitatis beginning man's firsthand exploration of the moon. This document is an edited record of the conversations between astronauts Neil Armstrong and Edwin "Buzz" Aldrin, Jr., at Tranquility Base, and Bruce McCandless at Mission Control in Houston during the approximately 22 hours spent on the lunar surface. It includes additional commentary during their return to Earth. It is a condensation hopefully of all the verbal data having geological significance. All discussions and observations documenting the lunar landscape, its geologic characteristics, the rocks and soils collected, and the photographic record are retained along with supplementary remarks essential to the continuity of events during the mission. We have deleted the words of mechanical housekeeping and engineering data, attempting not to lose the personal and philosophical aspects of this intensely human experience. The sources of this verbal transcript are the complete audio tapes recorded during the mission and the Technical Air-to-Ground Voice Transcription published by NASA. The voice record is listed chronologically given in days, hours, minutes, and seconds. These are the Ground Elapsed Times (GET) after launch from Kennedy Space Center which was 9:32 a.m. EDT on July 16, 1969. Figure 1 shows the vicinity of the landing site that was described, sampled, and photographed by the Apollo 11 crewmen.

Final Report↗
USGS2024

Empirical ground-motion basin response in the California Great Valley, Reno, Nevada, and Portland, Oregon

We assess how well the Next-Generation Attenuation-West 2 (NGA-West2) ground-motion models (GMMs), which are used in the US Geological Survey’s (USGS) National Seismic Hazard Model (NSHM) for crustal faults in the western United States, predict the observed basin response in the Great Valley of California, the Reno basin in Nevada, and Portland and Tualatin basins in Oregon. These GMMs rely on site parameters such as the time-averaged shear-wave velocity ( V S ) in the upper 30 m of Earth’s crust ( V S30 ) and depths to 1.0 and 2.5 km/s shear-wave isosurfaces ( Z 1.0 and Z 2.5 ) to capture basin effects and were developed using observations and simulations primarily from the Los Angeles region in southern California. Using ground-motion records from mostly small-to-moderate earthquakes and mixed-effects regression analysis, we find that the GMMs perform well with our local basin-depth models for the California Great Valley. With our local basin-depth models for Reno, the GMMs do not perform as well for this relatively shallow basin and exhibit little sensitivity to the basin parameters used in the NGA-West2 GMMs. We also find good performance for the local Z 1.0 model across the Portland region, whereas the local Z 2.5 model provides little predictive power except at sites in the deepest part of the Tualatin basin. Additional work could improve the performance of the site and basin terms in the NGA-West2 GMMs for regions with geologic structure different than the deep basins in southern California and the Great Valley. In addition, we find significant discrepancies among the GMMs in how the uncertainty in the ground motion varies with basin depth and pseudospectral period. Our results can help guide seismic hazard analyses on whether to include these local basin-depth models.

Earthquake Spectra↗
USGS2015

Multiscale geophysical imaging of the critical zone

Details of Earth's shallow subsurface—a key component of the critical zone (CZ)—are largely obscured because making direct observations with sufficient density to capture natural characteristic spatial variability in physical properties is difficult. Yet this inaccessible region of the CZ is fundamental to processes that support ecosystems, society, and the environment. Geophysical methods provide a means for remotely examining CZ form and function over length scales that span centimeters to kilometers. Here we present a review highlighting the application of geophysical methods to CZ science research questions. In particular, we consider the application of geophysical methods to map the geometry of structural features such as regolith thickness, lithological boundaries, permafrost extent, snow thickness, or shallow root zones. Combined with knowledge of structure, we discuss how geophysical observations are used to understand CZ processes. Fluxes between snow, surface water, and groundwater affect weathering, groundwater resources, and chemical and nutrient exports to rivers. The exchange of gas between soil and the atmosphere have been studied using geophysical methods in wetland areas. Indirect geophysical methods are a natural and necessary complement to direct observations obtained by drilling or field mapping. Direct measurements should be used to calibrate geophysical estimates, which can then be used to extrapolate interpretations over larger areas or to monitor changing processes over time. Advances in geophysical instrumentation and computational approaches for integrating different types of data have great potential to fill gaps in our understanding of the shallow subsurface portion of the CZ and should be integrated where possible in future CZ research.

Reviews of Geophysics↗
USGS2013

U–Pb, Rb–Sr, and U-series isotope geochemistry of rocks and fracture minerals from the Chalk River Laboratories site, Grenville Province, Ontario, Canada

As part of the Geologic Waste Management Facility feasibility study, Atomic Energy of Canada Ltd. (AECL) is evaluating the suitability of the Chalk River Laboratories (CRL) site in Ontario, situated in crystalline rock of the southwestern Grenville Province, for the possible development of an underground repository for low- and intermediate-level nuclear waste. This paper presents petrographic and trace element analyses, U–Pb zircon dating results, and Rb–Sr, U–Pb and U-series isotopic analyses of gneissic drill core samples from the deep CRG-series characterization boreholes at the CRL site. The main rock types intersected in the boreholes include hornblende–biotite (±pyroxene) gneisses of granitic to granodioritic composition, leucocratic granitic gneisses with sparse mafic minerals, and garnet-bearing gneisses with variable amounts of biotite and/or hornblende. The trace element data for whole-rock samples plot in the fields of within-plate, syn-collision, and volcanic arc-type granites in discrimination diagrams used for the tectonic interpretation of granitic rocks. Zircons separated from biotite gneiss and metagranite samples yielded SHRIMP-RG U–Pb ages of 1472 ± 14 (2 σ ) and 1045 ± 6 Ma, respectively, in very good agreement with widespread Early Mesoproterozoic plutonic ages and Ottawan orogeny ages in the Central Gneiss Belt. The Rb–Sr, U–Pb, and Pb–Pb whole-rock errorchron apparent ages of most of the CRL gneiss samples are consistent with zircon U–Pb age and do not indicate substantial large-scale preferential element mobility during superimposed metamorphic and water/rock interaction processes. This may confirm the integrity of the rock mass, which is a positive attribute for a potential nuclear waste repository. Most 234 U/ 238 U activity ratios (AR) in whole rock samples are within errors of the secular equilibrium value of one, indicating that the rocks have not experienced any appreciable U loss or gain within the past 1 Ma. However, 234 U/ 238 U AR in fracture mineral samples collected down to borehole lengths of about 740 m deviate from the secular equilibrium value and 234 U/ 238 U model ages calculated for fracture mineral samples showing excess 234 U range from 593 to 1415 ka, thus providing evidence of fracture flow in the associated bedrock during the past 1.5 Ma. Rare earth element patterns are variable in fracture-filling calcites and Fe oxides/hydroxides but are similar to those observed in associated whole-rock samples. The observed Ce anomalies are very small ( Ce N / Ce N ∗ ≈ 1 "> CeN/CeN∗≈1 ), do not vary with depth, and, therefore, do not contain evidence that the studied fracture minerals precipitated from oxidizing waters at the conceptual depth of a repository.

Ontario↗
USGS2019

Stormquakes

Seismic signals from ocean-solid Earth interactions are ubiquitously recorded on our planet. However, these wavefields are typically incoherent in most frequency bands limiting their utilization for understanding ocean dynamics or solid Earth properties. In contrast, we find that during large storms such as hurricanes and Nor’easters the interaction of long-period ocean waves with shallow seafloor features located near the edge of continental shelfs, known as ocean banks, excites coherent transcontinental Rayleigh wave packets in the 20 to 50 s period band. These “stormquakes” migrate coincident with the storms, but are spatiotemporally focused seismic point sources with equivalent earthquake magnitudes that can be greater than 3.5. Stormquakes thus provide new coherent sources to investigate Earth structure in locations that typically lack both seismic instrumentation and earthquakes. Moreover, they provide a new geophysical observable with high spatial and temporal resolution with which to investigate ocean wave dynamics during large storms.

Geophysical Research Letters↗
USGS2006

Comment on “222Rn Premonitory signals for earthquakes?”

Radon and other terrestrial gases in groundwater and soil air have been studied in recent years in seismically active areas, especially in China, Japan, the former Soviet Union, and the United States, in search of changes that may be useful for earthquake prediction. Concentrations of radon, helim, hydrogen, mercury vapor, carbon dioxide, and several other gases have been found to be high along active faults. Faults may thus be paths of least resistance along which terrestrial gases generated or stored in the Earth can escape into the atmosphere. Gas-concentration changes with durations of hours to months have been observed before and after earthquakes at relatively few favorably situated sites, and their patterns have not always been repeated. Some of these sites are located as far as several hundred kilometers away from earthquake epicenters. These “sensitive” sites are generally situated along active faults, whether the seismogenic ones or not, especially at their intersections or bends, or at other structurally weak zones, possibly because of concentration of tectonic strain and abundance of pore fluids at such places (see review papers by King , 1986, 1990). Due to this selectivity and for two other important reasons, studies of terrestrial gases for earthquake prediction have not been viewed favorably by most geophysicists, especially those in the United States.

Eos Science News↗
USGS2017

A land data assimilation system for sub-Saharan Africa food and water security applications

Seasonal agricultural drought monitoring systems, which rely on satellite remote sensing and land surface models (LSMs), are important for disaster risk reduction and famine early warning. These systems require the best available weather inputs, as well as a long-term historical record to contextualize current observations. This article introduces the Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS), a custom instance of the NASA Land Information System (LIS) framework. The FLDAS is routinely used to produce multi-model and multi-forcing estimates of hydro-climate states and fluxes over semi-arid, food insecure regions of Africa. These modeled data and derived products, like soil moisture percentiles and water availability, were designed and are currently used to complement FEWS NET’s operational remotely sensed rainfall, evapotranspiration, and vegetation observations. The 30+ years of monthly outputs from the FLDAS simulations are publicly available from the NASA Goddard Earth Science Data and Information Services Center (GES DISC) and recommended for use in hydroclimate studies, early warning applications, and by agro-meteorological scientists in Eastern, Southern, and Western Africa.

Scientific Data↗
USGS2024

Polar paleoenvironmental perspectives on modern climate change

In today’s rapidly changing climate, society needs a better understanding of climate impacts on sea level, ice sheets and glaciers, sea ice, ocean circulation, ecosystems, biodiversity, and other aspects of planet Earth. Paleoenvironmental records provide a unique and invaluable source of insight into these complex issues, and place recent observations into a broader historical context. This essay discusses why paleoclimate reconstructions from polar regions provide critical information to help anticipate possible future climate impacts. By highlighting some key research examples, this essay explains the value of expanding proxy-based research in Arctic/Antarctic regions, and makes a case for paying greater attention to the lessons already distilled from it.

PLOS Climate↗
USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗
USGS1981

Rare-earth element geochemistry of the island-arc volcanic rocks of Rabaul and Talasea, New Britain

The island-arc volcanic rocks of Rabaul and Talasea, New Britain, range in composition from basalt through rhyolite. Rare-earth elements have been determined by mass-spectrometric isotope dilution in 16 samples. Chondrite-normalized rare-earth element patterns are distinct for each volcanic center, but all are relatively flat (Ce/Yb E.F. = 1.1 to 3.0). Within each center, rare-earth element concentrations increase from basalt to dacite, and Eu anomalies become progressively more negative from basalt to dacite (Eu/Eu* = 1.0 to 0.8). Lavas inferred to contain cumulate phenocrysts show positive Eu anomalies (Eu/Eu* of up to 1.2). The observed variations are consistent with fractional crystallization of basalt or basaltic andesite by removal or accumulation of the observed phenocrysts, including olivine, plagioclase, pyroxene, and opaque minerals. Quantitative trace-element models for Talasea lavas indicate precipitation of 50 wt percent of phenocrysts from basalt to produce andesite and an additional 22 wt percent to produce dacite. A total-equilibrium model provides the best approximation to observed concentrations, implying that the process took place at a depth sufficient to allow slow cooling of the magma as precipitation proceeded. The “least fractionated” basalts and basaltic andesites are similar in rare-earth element pattern to those of corresponding type suites in other oceanic-island arcs, implying a widespread uniformity of source and process in the production of each magma type.

Geological Society of America Bulletin↗
USGS2024

Constraints from GPS measurements on plate coupling within the Makran subduction zone and tsunami scenarios in the western Indian Ocean

Plate-coupling estimates and previous seismicity indicate that portions of the Makran megathrust of southern Pakistan and Iran are partially coupled and have the potential to produce future magnitude 7+ earthquakes. However, the GPS observations needed to constrain coupling models are sparse and lead to an incomplete understanding of regional earthquake and tsunami hazard. In this study, we assess GPS velocities for plate coupling of the Makran subduction zone with specific attention to model resolution and the accretionary prism rheology. We use finite element model-derived Green's functions to invert for the interseismic slip deficit under both elastic and viscoelastic Earth assumptions. We use the model resolution matrix to characterize plate-coupling scenarios that are consistent with the limited spatial resolution afforded by GPS observations. We then forward model the corresponding tsunami responses at major coastal cities within the western Indian Ocean basin. Our plate-coupling results show potential segmentation of the megathrust with varying coupling from west to east, but do not rule out a scenario where the entire length of the megathrust could rupture in a single earthquake. The full subduction zone rupture scenarios suggest that the Makran may be able to produce earthquakes up to M w 9.2. The corresponding tsunami model from the largest earthquake event ( M w 9.2) estimates maximum wave heights reaching 2–5 m at major port cities in the northern Arabian Sea region. Cities on the west coast of India are less affected (1–2 m). Coastlines bounding eastern Africa, and the Strait of Hormuz, are the least affected (<1 m).

Geophysical Journal International↗
USGS2003

Variability in geyser eruptive timing and its causes: Yellowstone National Park

Field data from Upper Geyser Basin, Yellowstone, indicate that geyser frequency is less sensitive to elastic deformation than might be surmised from a review of the literature. Earth-tide influences are not identifiable in any of the geysers we monitored. Though atmospheric-pressure influences are observed, only long-period variations on the order of 5 mBars or greater seem to influence geyser frequency. Long-distance interconnections between geysers are common and add to the difficulty of identifying strain influences. Additional variations in geyser periodicity may be governed by the internal dynamics of the geysers rather than external influences.

Geophysical Research Letters↗
USGS2012

Hotspot of accelerated sea-level rise on the Atlantic coast of North America

Climate warming does not force sea-level rise (SLR) at the same rate everywhere. Rather, there are spatial variations of SLR superimposed on a global average rise. These variations are forced by dynamic processes, arising from circulation and variations in temperature and/or salinity, and by static equilibrium processes, arising from mass redistributions changing gravity and the Earth's rotation and shape. These sea-level variations form unique spatial patterns, yet there are very few observations verifying predicted patterns or fingerprints. Here, we present evidence of recently accelerated SLR in a unique 1,000-km-long hotspot on the highly populated North American Atlantic coast north of Cape Hatteras and show that it is consistent with a modelled fingerprint of dynamic SLR. Between 1950&ndash;1979 and 1980&ndash;2009, SLR rate increases in this northeast hotspot were ~ 3&ndash;4 times higher than the global average. Modelled dynamic plus steric SLR by 2100 at New York City ranges with Intergovernmental Panel on Climate Change scenario from 36 to 51 cm (ref. 3); lower emission scenarios project 24&ndash;36 cm (ref. 7). Extrapolations from data herein range from 20 to 29 cm. SLR superimposed on storm surge, wave run-up and set-up will increase the vulnerability of coastal cities to flooding, and beaches and wetlands to deterioration.

Atlantic Coast, Cape Cod, Cape Hatteras, North Ame↗
USGS2011

Rapid cooling rates at an active mid-ocean ridge from zircon thermochronology

Oceanic spreading ridges are Earth's most productive crust generating environment, but mechanisms and rates of crustal accretion and heat loss are debated. Existing observations on cooling rates are ambiguous regarding the prevalence of conductive vs. convective cooling of lower oceanic crust. Here, we report the discovery and dating of zircon in mid-ocean ridge dacite lavas that constrain magmatic differentiation and cooling rates at an active spreading center. Dacitic lavas erupted on the southern Cleft segment of the Juan de Fuca ridge, an intermediate-rate spreading center, near the intersection with the Blanco transform fault. Their U–Th zircon crystallization ages (29.3 − 4.6 + 4.8 ka; 1σ standard error s.e.) overlap with the (U–Th)/He zircon eruption age (32.7 ± 1.6 ka) within uncertainty. Based on similar 238 U− 230 Th disequilibria between southern Cleft dacite glass separates and young mid-ocean ridge basalt (MORB) erupted nearby, differentiation must have occurred rapidly, within ~ 10–20 ka at most. Ti-in-zircon thermometry indicates crystallization at 850–900 °C and pressures > 70–150 MPa are calculated from H 2 O solubility models. These time-temperature constraints translate into a magma cooling rate of ~ 2 × 10 − 2 °C/a. This rate is at least one order-of-magnitude faster than those calculated for zircon-bearing plutonic rocks from slow spreading ridges. Such short intervals for differentiation and cooling can only be resolved through uranium-series ( 238 U– 230 Th) decay in young lavas, and are best explained by dissipating heat convectively at high crustal permeability.

Juan de Fuca ridge and the Blanco transform fault ↗
USGS2009

Mesoscale raised rim depressions (MRRDs) on Earth: A review of the characteristics, processes, and spatial distributions of analogs for Mars

Fields of mesoscale raised rim depressions (MRRDs) of various origins are found on Earth and Mars. Examples include rootless cones, mud volcanoes, collapsed pingos, rimmed kettle holes, and basaltic ring structures. Correct identification of MRRDs on Mars is valuable because different MRRD types have different geologic and/or climatic implications and are often associated with volcanism and/or water, which may provide locales for biotic or prebiotic activity. In order to facilitate correct identification of fields of MRRDs on Mars and their implications, this work provides a review of common terrestrial MRRD types that occur in fields. In this review, MRRDs by formation mechanism, including hydrovolcanic (phreatomagmatic cones, basaltic ring structures), sedimentological (mud volcanoes), and ice-related (pingos, volatile ice-block forms) mechanisms. For each broad mechanism, we present a comparative synopsis of (i) morphology and observations, (ii) physical formation processes, and (iii) published hypothesized locations on Mars. Because the morphology for MRRDs may be ambiguous, an additional tool is provided for distinguishing fields of MRRDs by origin on Mars, namely, spatial distribution analyses for MRRDs within fields on Earth. We find that MRRDs have both distinguishing and similar characteristics, and observation that applies both to their mesoscale morphology and to their spatial distribution statistics. Thus, this review provides tools for distinguishing between various MRRDs, while highlighting the utility of the multiple working hypotheses approach. ?? 2008 Elsevier Ltd.

Planetary and Space Science↗
USGS2002

Challenges in observational seismology

Earthquake seismology became a quantitative scientific discipline after instruments were developed to record seismic waves in the late 19th century ( Dewey and Byerly, 1969 ; Chapter 1 by Agnew). Earthquake seismology is essentially based on field observations. The great progress made in the past several decades was primarily due to increasingly plentiful and high-quality data that are readily distributed. Our ability to collect, process, and analyze earthquake data has been accelerated by advances in electronics, communications, computers, and software (see Chapter 85 edited by Snoke and Garcia-Fernandez). Instrumental observation of earthquakes has been carried out for a little over 100 years by seismic stations and networks of various sizes, from local to global scales (see Chapter 87 edited by Lahr and van Eck). The observed data have been used, for example, (1) to compute the source parameters of earthquakes, (2) to determine the physical properties of the Earth's interior, (3) to test the theory of plate tectonics , (4) to map active faults, (5) to infer the nature of damaging ground shaking, and (6) to carry out seismic hazard analysis. Construction of a satisfactory theory of the earthquake process has not yet been achieved within the context of physical laws. Good progress, however, has been made in building a physical foundation of the earthquake source process, partly as a result of research directed toward earthquake prediction. This chapter is intended for a general audience. Technical details are not given, but relevant references and chapters in this Handbook are referred to. The first part of this chapter presents a brief overview of the observational aspects of earthquake seismology, concentrating on instrumental observations of seismic waves generated by earthquakes (i.e., seismic monitoring), and readers are referred to Chapter 49 by Musson and Cecic for noninstrumental observations. A few key developments and practices are summarized by taking a general view, since many national and regional developments have been chronicled in national and institutional reports (see Chapter 79 edited by Kisslinger). In the latter part of this chapter, the nature of seismic monitoring and some challenges in observational seismology are discussed from a personal perspective. Comments of a technical or philosophical nature are given in the Notes at the end of the chapter, and they are referenced by superscript numbers in the text.

International Geophysics↗
USGS2021

Reconstructing the dynamics of the highly similar May 2016 and June 2019 Iliamna Volcano, Alaska ice–rock avalanches from seismoacoustic data

Surficial mass wasting events are a hazard worldwide. Seismic and acoustic signals from these often remote processes, combined with other geophysical observations, can provide key information for monitoring and rapid response efforts and enhance our understanding of event dynamics. Here, we present seismoacoustic data and analyses for two very large ice–rock avalanches occurring on Iliamna Volcano, Alaska (USA), on 22 May 2016 and 21 June 2019. Iliamna is a glacier-mantled stratovolcano located in the Cook Inlet, ∼200 km from Anchorage, Alaska. The volcano experiences massive, quasi-annual slope failures due to glacial instabilities and hydrothermal alteration of volcanic rocks near its summit. The May 2016 and June 2019 avalanches were particularly large and generated energetic seismic and infrasound signals which were recorded at numerous stations at ranges from ∼9 to over 600 km. Both avalanches initiated in the same location near the head of Iliamna's east-facing Red Glacier, and their ∼8 km long runout shapes are nearly identical. This repeatability – which is rare for large and rapid mass movements – provides an excellent opportunity for comparison and validation of seismoacoustic source characteristics. For both events, we invert long-period (15–80 s) seismic signals to obtain a force-time representation of the source. We model the avalanche as a sliding block which exerts a spatially static point force on the Earth. We use this force-time function to derive constraints on avalanche acceleration, velocity, and directionality, which are compatible with satellite imagery and observed terrain features. Our inversion results suggest that the avalanches reached speeds exceeding 70 m s−1, consistent with numerical modeling from previous Iliamna studies. We lack sufficient local infrasound data to test an acoustic source model for these processes. However, the acoustic data suggest that infrasound from these avalanches is produced after the mass movement regime transitions from cohesive block-type failure to granular and turbulent flow – little to no infrasound is generated by the initial failure. At Iliamna, synthesis of advanced numerical flow models and more detailed ground observations combined with increased geophysical station coverage could yield significant gains in our understanding of these events.

Alaska↗