Seasonal variability in source activity and atmospheric temperature were retrieved from 11 years of infrasonic ambient noise. Variable lag times between an array and single microbarometer were obtained from envelopes of cross‐correlation functions. Beamforming and one‐bit normalization significantly enhanced the stationary phase. Both microbaroms and surf appeared abundantly present, in the 0.1–2.0 Hz frequency band. Modeling revealed both tropospheric and stratospheric propagation of the infrasound, following traditional and more unconventional propagation mechanisms. Virtual source‐receiver refractions from stratospheric altitudes appeared a plausible explanation for the unusual short lag times, which allows for new ways to passively probe the stratosphere.

Plain Language Summary
Microbarometers measure infrasound from the background noise field, that is, ambient noise. This ambient noise is continuously present and generated by oceanic waves. The sound speed is a function of wind and temperature and determines how fast the ambient noise travels from one to another microbarometer. The so‐called travel time can be obtained by comparing the recordings of ambient noise of different microbarometers. The use of advanced processing techniques with multiple microbarometers was explored and enhanced the results. The observed travel times provided evidence for propagation through the upper atmosphere, that is, the stratosphere. New ways are proposed to apply surface‐based infrasound recordings in passively probing stratospheric wind and temperature conditions.

The oceans are filled with acoustic waves, which are trapped in a low-velocity layer at about 1 km water depth. The sound speed depends on the temperature. From hydroacoustic arrays, travel times can be obtained through cross-correlating transient signals. Hydroacoustic station H10 (Ascension Island) appeared capable of measuring deep ocean temperature change. A decrease in modeled and observed travel times of −0.002 s/yr was derived between two arrays, corresponding to a warming of 0.007 °C/yr, at about 900 m water depth over two decades. As such, acoustic waves provide an independent and passively acquired measure of the temperature change in the deep ocean.

Eruptions at continental basaltic volcanoes can take and combine various forms, including lava lakes, lava flows and fountaining, explosions or structural collapses. Aside from a few well-instrumented cases worldwide, accurately reconstructing eruptive scenarios is hampered by the lack of detailed visual observations. However, volcanoes also have their own acoustic signatures composed of low-pitched inaudible sounds, called infrasounds. Here we analyze infrasound records close to Nyiragongo volcano (D.R. Congo) (<20 km) up to Kenya (~800 km), which are converted using acoustic numerical modeling into time-lapse observations of the catastrophic drainage of the world’s largest lava lake on 22 May 2021. The emitted infrasounds also enable tracking of fissure openings and lava eruptions along the flank, occurring simultaneously with the lava lake drainage. This striking example supports the growing role of infrasound as a key component of volcano monitoring and early-warning systems, as it provides unique information inaccessible to other ground-based instruments.

Mitigation of the threat from airbursting asteroids requires an understanding of the potential risk they pose for the ground. How asteroids release their kinetic energy in the atmosphere is not well understood due to the rarity of large impacts. Here we present a comprehensive, space-to-laboratory characterization of an impact of an L chondrite, which represents a common type of Earth-impacting asteroid. Small asteroid 2023 CX1 was detected in space and predicted to impact over Normandy, France, on 13 February 2023. Observations from several independent sensors and reduction techniques revealed an unusual but potentially high-risk fragmentation behaviour. The nearly spherical 650 ± 160 kg (72 ± 6 cm diameter) asteroid catastrophically fragmented at a dynamic pressure of 4 MPa around 28 km altitude, releasing 98% of its total energy in a concentrated region of the atmosphere. The resulting shock wave was spherical, not cylindrical, and released more energy closer to the ground. This type of fragmentation increases the risk of substantial damage at ground level. These results warrant consideration for a planetary defence strategy for cases where a >3–4 MPa dynamic pressure is expected, including planning for evacuation of areas beneath anticipated disruption locations.

Climate change is increasingly predisposing polar regions to large landslides. Tsunamigenic landslides have occurred recently in Greenland (Kalaallit Nunaat), but none have been reported from the eastern fjords. In September 2023, we detected the start of a 9-day-long, global 10.88-millihertz (92-second) monochromatic very-long-period (VLP) seismic signal, originating from East Greenland. In this study, we demonstrate how this event started with a glacial thinning–induced rock-ice avalanche of 25 × 106 cubic meters plunging into Dickson Fjord, triggering a 200-meter-high tsunami. Simulations show that the tsunami stabilized into a 7-meter-high long-duration seiche with a frequency (11.45 millihertz) and slow amplitude decay that were nearly identical to the seismic signal. An oscillating, fjord-transverse single force with a maximum amplitude of 5 × 1011 newtons reproduced the seismic amplitudes and their radiation pattern relative to the fjord, demonstrating how a seiche directly caused the 9-day-long seismic signal. Our findings highlight how climate change is causing cascading, hazardous feedbacks between the cryosphere, hydrosphere, and lithosphere.

In August 2021, Askja volcano in Iceland returned to the spotlight after a sudden onset of rapid uplift followed decades of continuous subsidence. In this study the extended record of microgravity data from Askja between 1988 to 2017 is revisited, and new microgravity data from 2021 and 2022 are introduced, which were collected after the uplift had started. Askja caldera had been steadily subsiding since at least 1984 and was characterised by a net decrease in microgravity, potentially signalling the contraction of its magma chamber or eviction of magma either laterally or to deeper levels. The microgravity data indicate that despite ongoing subsidence between 2017 and early 2021, a significant gravity increase can be detected in the center of the caldera between 2017 and August 2021. This increase may be introduced during – or leading up to – the period of uplift. The new microgravity data also indicate that during the period of 40 cm uplift after August 2021 to fall 2022, gravity changes approach the free-air gradient, suggesting subsurface density decreases as a driving process. This process may relate to the vesiculation of magma previously emplaced in the volcano roots, a change in the hydrothermal system, or replacement of dense basaltic magma with less dense rhyolitic magma, or a combination of these processes. However, uncertainties for this period are elevated and may obscure a gravity signal expected from additional mass accumulation. The timing and high uncertainties of some campaigns make it challenging to be conclusive on the driving process behind the uplift, but future microgravity campaigns could help solve the ambiguity. The study also provides a description of potential pitfalls in microgravity campaigns and recommendations on how the reliability of microgravity data can be improved.

A period of 18 years of infrasonic recordings was analyzed from a microbarometer array (I18DK) in northwestern Greenland, near Qaanaaq. A huge number of infrasonic detections, over 700,000, have been made in I18DKs soundscape during the Arctic summers. Simultaneously identified were both calving events from marine-terminating glaciers and discharge related acoustics from a land-terminating glacier. This infrasonic activity is correlated to sea-surface and atmospheric temperature, respectively. Inter-yearly to daily variations were retrieved showing a strong variability in infrasonic detection rates and hence glacier activity. The highest number of infrasonic detections were found in recent years from the land-terminating glacier. The latter is supported by actual discharge measurements and partly by a discharge model. It is concluded that monitoring infrasound from glaciers can complement other techniques to remotely and passively get insights into glacier dynamics with high temporal and spatial resolution.

Results from nine microgravity campaigns from Kı̄lauea, Hawaiʻi, spanning most of the volcano's 2008–2018 summit eruption, indicate persistent mass accumulation at shallow levels. A weighted least squares approach is used to recover microgravity results from a network of benchmarks around Kı̄lauea's summit, eliminate instrumental drift, and restore suspected data tares. A total mass of 1.9 × 10¹¹ kg was determined from these microgravity campaigns to have accumulated below Kı̄lauea Caldera during 2009–2015 at an estimated depth of 1.3 km below sea level. Only a fraction of this mass is reflected in surface deformation, and this is consistent with previously reported discrepancies between subsurface mass accumulation and observed surface deformation. The discrepancy, amongst other independent evidence from gas emissions, seismicity, and continuous gravimetry, indicate densification of magma in the reservoirs below the volcano summit. This densification may have been driven by degassing through the summit vent. It is hypothesized that during the final years of the summit eruption, magma densification resulted in a buildup of pressure in the reservoirs that may have contributed to the lower East Rift Zone outbreak of 2018. The observed mass accumulation beneath Kı̄lauea could not have been detected through other techniques and illustrates the importance of microgravity measurements in volcano monitoring.

The ambient infrasonic noise field is complex due to the interference of spatially distributed infrasound sources. Microbaroms are one of the most dominant omnipresent infrasonic sources within this wavefield. These microbaroms are generated by nonlinear ocean surface wave interactions, and have a characteristic and continuous signature within the infrasound spectrum. Under noisy conditions, microbaroms can mask infrasonic signals of interest, such as infrasound from volcanoes or explosions, which limits detection and identification of such sources. This study performs an infrasonic climatology for infrasound array I23FR, using five years of data between 2015-2020. The array is located on the Kerguelen Islands, within the Southern Ocean, and is part of the International Monitoring System for the verification of the Comprehensive Nuclear-Test-Ban Treaty. The climatology analysis addresses the expected ambient noise levels, propagation paths and potential sources within the vicinity of an infrasound sensor. Time- and frequency-domain beamforming methods have been applied to analyse the infrasonic wavefield from the I23FR observations. A recently introduced method is applied to compute so-called soundscapes, to be compared with beamform results. Although the comparison indicates a disagreement in amplitude, there is a good agreement in directionality and frequency between both.

We report on the extraction of deep ocean travel time variations from time-lapse cross-correlations between a hydrophone station and a three-component broadband seismometer. The signals we cross-correlate in this study result from repeated activity by the Monowai seamount, one of the most active submarine volcanoes of the Tonga-Kermadec ridge. In particular, we introduce a specific workflow to exploit repetitive hydroacoustic underwater source activity, which we detail to such an extent that it serves as an example (or “cookbook”). For this reason, we have made the source code publicly available. The workflow proposed in this study (a) overcomes differences in instrument sensitivity and sample rates, (b) involves the selection of eligible cross-correlations based on a source activity criterium as well as slowness analysis, and (c) extracts the travel time variations in distinct frequency bands. In our case, the two frequency bands are 3–6 and 6–12 Hz. We find that the estimated travel time variations in both frequency bands consist of a complex periodic pattern superimposed on a robust linear trend. This linear trend is decreasing, which we attribute to increasing water temperatures along the propagation path of the hydroacoustic signals.

Geophysical studies and real-time monitoring of natural hazards, such as volcanic eruptions or severe weather events, benefit from the joint analysis of multiple geophysical parameters. However, typical geophysical measurement platforms still provide logging solutions for a single parameter, due to different community standards and the higher cost per added sensor. In this work, the Infrasound and Environmental Atmospheric data Recorder (INFRA-EAR) is presented, which has been designed as a low-cost mobile multidisciplinary measurement platform for geophysical monitoring. In particular, the platform monitors infrasound but concurrently measures barometric pressure, accelerations, and wind flow and uses the Global Positioning System (GPS) to position the platform. Due to its digital design, the sensor platform can be readily integrated with existing geophysical data infrastructures and be embedded in geophysical data analysis. The small dimensions and low cost per unit allow for unconventional, experimental designs, for example, high-density spatial sampling or deployment on moving measurement platforms. Moreover, such deployments can complement existing high-fidelity geophysical sensor networks. The platform is designed using digital micro-electromechanical system (MEMS) sensors embedded on a printed circuit board (PCB). The MEMS sensors on the PCB are a GPS, a three-component accelerometer, a barometric pressure sensor, an anemometer, and a differential pressure sensor. A programmable microcontroller unit controls the sampling frequency of the sensors and data storage. A waterproof casing is used to protect the mobile platform against the weather. The casing is created with a stereolithography (SLA) Formlabs 3D printer using durable resin. Thanks to low power consumption (9Wh over 25 d), the system can be powered by a battery or solar panel. Besides the description of the platform design, we discuss the calibration and performance of the individual sensors.

Seabirds are amongst the most mobile of all animal species and spend large amounts of their lives at sea. They cross vast areas of ocean that appear superficially featureless, and our understanding of the mechanisms that they use for navigation remains incomplete, especially in terms of available cues. In particular, several large-scale navigational tasks, such as homing across thousands of kilometers to breeding sites, are not fully explained by visual, olfactory or magnetic stimuli. Low-frequency inaudible sound, i.e., infrasound, is ubiquitous in the marine environment. The spatio-temporal consistency of some components of the infrasonic wavefield, and the sensitivity of certain bird species to infrasonic stimuli, suggests that infrasound may provide additional cues for seabirds to navigate, but this remains untested. Here, we propose a framework to explore the importance of infrasound for navigation. We present key concepts regarding the physics of infrasound and review the physiological mechanisms through which infrasound may be detected and used. Next, we propose three hypotheses detailing how seabirds could use information provided by different infrasound sources for navigation as an acoustic beacon, landmark, or gradient. Finally, we reflect on strengths and limitations of our proposed hypotheses, and discuss several directions for future work. In particular, we suggest that hypotheses may be best tested by combining conceptual models of navigation with empirical data on seabird movements and in-situ infrasound measurements.

In the days following the January 12, 2010 M_w 7 Haiti earthquake the shaking intensity near the epicenter was overestimated and the spatial extent of the potentially damaging shaking was underestimated. This was due to the lack of seismometers in the near-source region at the time of the earthquake. Besides seismic waves, earthquakes generate infrasound, i.e., inaudible acoustic waves in the atmosphere. Here we show that infrasound signals, detected at distant ground-based stations, can be used to generate a map of the acoustic intensity, which is proportional to the shaking intensity. This is demonstrated with infrasound from the 2010 Haiti earthquake detected in Bermuda, over 1700 km away. Wavefront parameters are retrieved in a beamforming process and are backprojected to map the measured acoustic intensity to the source region. The backprojection process accounts for horizontal advection effects due to winds and inherent uncertainties with regard to the time of detection and the back azimuth resolution. Furthermore, we resolve the ground motion polarity in the epicentral region and use synthetics generated by an extended infrasound source model to support this result. Infrasound measurements are conducted globally for the verification of the Comprehensive Nuclear-Test-Ban Treaty and although the network was designed to provide global coverage for nuclear explosions in the atmosphere, it is shown in this paper that there is also global coverage for the estimation of acoustic shaking intensity. In this study, we lay the groundwork that can potentially make infrasound-based ShakeMaps a useful tool alongside conventional ShakeMaps and a valuable tool for earthquake disaster mitigation in sparsely monitored regions.

A method is introduced to reconstruct microbarom soundscapes in absolute values. The soundscapes are compared to remote infrasound recordings from infrasound array I23FR (Kerguelen Island) and in situ recordings by the INFRA-EAR, a biologger deployed near the Crozet Islands. The reconstruction method accounts for all-acoustic contributions, divided into evanescent microbaroms (detectable directly above the source) and propagating microbaroms (detectable over long ranges). It is computed by integrating acoustic intensities over the ocean surface, convolved with the transfer function quantifying the propagation losses and propagation time. The reconstructed soundscapes are found within 2.7 dB for (Formula presented.) of the measurements in the microbarom band of 0.1–0.3 Hz. Infrasonic soundscapes are essential for understanding the ambient infrasonic noise field and are a basic need for applications, such as atmospheric remote sensing, natural hazard monitoring, and verification of the Comprehensive Nuclear-Test-Ban Treaty.

The June 2, 2018 impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. Twenty‐three meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as an HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of ~156 cm diameter with high bulk density ~2.85 g cm−3, a relatively low albedo pV ~ 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of ~0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth‐impacting orbit via the ν6 resonance. The impact that ejected 2018 LA in an orbit toward Earth occurred 22.8 ± 3.8 Ma ago. Zircons record a concordant U‐Pb age of 4563 ± 11 Ma and a consistent 207Pb/206Pb age of 4563 ± 6 Ma. A much younger Pb‐Pb phosphate resetting age of 4234 ± 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin.

In seismology, the depth of a near-surface source is hard to estimate in the absence of local stations. The depth-yield trade-off leads to significant uncertainties in the source's depth and strength estimations. Long-range infrasound propagation from an underwater or underground source is very sensitive to variations in the source's depth and strength. This characteristic is employed in an infrasound based inversion for the submerged source parameters. First, a Bayesian inversion scheme is tested under the variations of the number of stations, the signal's frequency band, and the signal-to-noise ratio (SNR). Second, an ensemble of realistic perturbed atmospheric profiles is used to investigate the effect of atmospheric uncertainties on the inversion results. Results show that long-range infrasound signals can be used to estimate the depth and strength of an underwater source. Using a broadband signal proved to be a fundamental element to obtain the real source parameters, whereas the SNR was secondary. Multiple station inversions perform better than one-station inversions; however, variations in their position can lead to source strength estimations with uncertainties up to 50%. Regardless of the number of stations, their positions, and SNRs, all of the estimated depths were within 10% from the real source depth.

Human activity causes vibrations that propagate into the ground as high-frequency seismic waves. Measures tomitigate the coronavirus disease 2019 (COVID-19) pandemic caused widespread changes in human activity,leading to a months-long reduction in seismic noise of up to 50%. The 2020 seismic noise quiet period is thelongest and most prominent global anthropogenic seismicnoise reduction on record. Although the reduction isstrongestatsurfaceseismometersinpopulatedareas, this seismic quiescence extends for many kilometersradially and hundreds of meters in depth. This quiet period provides an opportunity to detect subtle signalsfrom subsurface seismic sources that would have been concealed in noisier times and to benchmark sources ofanthropogenic noise. A strong correlation between seismic noise and independent measurements of humanmobility suggests that seismology provides an absolute, real-time estimate of human activities.

The detection and characterization of signals of interest in the presence of (in)coherent ambient noise is central to the analysis of infrasound array data. Microbaroms have an extended source region and a dynamical character. From the perspective of an infrasound array, these coherent noise sources appear as interfering signals which conventional beamform methods may not correctly resolve. This limits the ability of an infrasound array to dissect the incoming wavefield into individual components. In this paper, this problem will be addressed by proposing a high-resolution beamform technique in combination with the CLEAN algorithm. CLEAN iteratively selects the maximum of the f/k spectrum (i.e., following the Bartlett or Capon method) and removes a percentage of the corresponding signal from the cross-spectral density matrix. In this procedure, the array response is deconvolved from the f/k spectral density function. The spectral peaks are retained in a ’clean’ spectrum. A data-driven stopping criterion for CLEAN is proposed that relies on the framework of Fisher statistics. This allows the construction of an automated algorithm that continuously extracts coherent energy until the point is reached that only incoherent noise is left in the data. CLEAN is tested on a synthetic data-set and is applied to data from multiple IMS infrasound arrays. The results show that the proposed method allows for the identification of multiple microbarom source regions in the Northern Atlantic, that would have remained unidentified if conventional methods had been applied.

In seismology and ocean acoustics, the interface with the atmosphere is typically represented as a free surface. Similarly, these interfaces are considered as a rigid surface for infrasound propagation. This implies that seismic or acoustic waves are not transmitted into the atmosphere from subsurface sources, and vice versa. Nevertheless, infrasound generated by subsurface sources has been observed. In this work, seismo-acoustic modeling of infrasound propagation from underwater and underground sources will be presented. The fast field program (FFP) is used to model the seismo-acoustic coupling between the solid earth, the ocean, and the atmosphere under the variation of source and media parameters. The FFP model allows for a detailed analysis of the seismo-acoustic coupling mechanisms in frequency-wavenumber space. A thorough analysis of the coupling mechanisms reveals that evanescent wave coupling and leaky surface waves are the main energy contributors to long-range infrasound propagation. Moreover, it is found that source depth affects the relative amplitude of the tropospheric and stratospheric phases, which allows for source depth estimation in the future.

Measurements of seismo-acoustic events by collocated seismic and infrasound arrays allow for studying the two wavefields that were produced by the same event. However, some of the scientific and technical constraints on the building of the two technologies are different and may be contradicting. For the case of a new station, an optimal design that will satisfy the constraints of the two technologies can be found. However, in the case of upgrading an existing array by adding the complementing technology, the situation is different. The site location, the array configuration and physical constraints are fixed and may not be optimal for the complementing technology, which may lead to rejection of the upgrade. The International Monitoring System (IMS) for the verification of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) includes 37 seismic arrays and 51 infrasound arrays. Although the CTBT verification regime is fixed in the treaty, an upgrade of the existing arrays by adding more technologies is possible. The Mount Meron seismic array (MMAI), which is part of the IMS, is composed of 16 sites. Microbarometers were installed at five MMAI sites to form the Mount Meron infrasound array. Due to regulation and physical constraints, it was not possible to relocate the sites nor to install analogue noise reduction filters (i.e. a pipe array). In this study, it is demonstrated that the installation of the MMAI infrasound array is beneficial despite the non-optimal conditions. It is shown that the noise levels of the individual array sites are between the high and median global noise levels. However, we claim that the more indicative measures are the noise levels of the beams of interest, as demonstrated by analysing the microbaroms originated from the Mediterranean Sea. Moreover, the ability to detect events relevant to the CTBT is demonstrated by analysing man-made events during 2011 from the Libya region.