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.@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.@en

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.

@en

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.@en

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.@en

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.@en

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.@en

Infrasound has a long history of monitoring sudden stratospheric warmings. Several pioneering studies have focused on the various effects of a major warming on the propagation of infrasound. A clear transition has been made from observing anomalous signatures towards the use of these signals to study anomalies in upper atmospheric conditions. Typically, the infrasonic signature of a major warming corresponds to summer-like infrasound characteristics observed in midwinter. More subtile changes occur during a minor warming, recognisable by the presence of a bidirectional stratospheric duct or propagation through a warm stratosphere leading to small shadow zones. A combined analysis of all signal characteristics unravels the general stratospheric structure throughout the life cycle of the warming. A new methodology to evaluate the state of the atmosphere as represented by various weather and climate models is demonstrated. A case study comparing regional volcano infrasound with simulations using various forecast steps indicates significant differences in stratospheric forecast skill, associated with a data assimilation issue during the warming.@en

Gravity waves are an important part of the momentum budget of the atmosphere. Despite this, parameterizations of gravity wave spectra in atmospheric models are poorly constrained. Gravity waves are formed by jet streams, flow over topography, and convection, all of which produce pressure perturbations as they propagate over the Earth's surface, detectable by microbarometer arrays used for sensing infrasound. In this study, observations of gravity waves between 2007 and 2011 at an infrasound station in the Ivory Coast, West Africa, are combined with meteorological data to calculate parameters such as intrinsic phase speed and wavenumber. Through spectral analysis, the seasonal and daily variations in all gravity wave parameters are examined. The gravity wave back azimuth varies with the migration of the Intertropical Convergence Zone, a region of intense convection, supporting previous studies. Daily variations in gravity wave arrivals at the station can be linked to two distinct convective cycles over the land and ocean. This was achieved by combining the gravity wave parameters with lightning strikes detected by the Met Office's Arrival Time Difference lightning detection system. Noise generated by turbulence in the middle of the day was found to attenuate smaller pressure amplitude gravity waves, artificially amplifying the daily variations in some gravity wave parameters. Detection of daily and seasonal variations in gravity wave parameters has the potential be used to improve the representation of gravity wave spectra in atmospheric models.

@en

The ShakeMap is a key component in the initial relief efforts following an earthquake disaster. It depicts the distribution of shaking intensity in the epicentral region and is used to guide emergency responders to the region. In regions where seismic instrumentation is limited, such ShakeMaps are poorly constrained and can take days to generate. We show, that pseudo-ShakeMaps that indicate the relative shaking intensity, can be generated, within minutes, from enhanced processing and modeling of infrasound. Furthermore, the source mechanism can be retrieved. This is illustrated with infrasound from the 2010 Mw 7.0 Port-au-Prince, Haiti earthquake, detected in Bermuda, over 1700 km away from Haiti. The pseudo-ShakeMap and focal mechanism retrieved in this study are in good agreement with the USGS estimated ShakeMap and the Global CMT moment tensor. Such observations are made possible by: (1) An advanced array processing technique that enables the detection of coherent wavefronts, even when amplitudes are below the noise level, and (2) Backprojection of observed pressure perturbations to ground motions in the epicentral region while accounting for advection effects in the atmosphere. We support our observations with an example using the Rayleigh integral to generate synthetic waveforms from four quadrants of an earthquake focal mechanism. Synthetics are then processed to retrieve the relative sense of motion in each quadrant. The current infrasound networks routinely detect earthquakes and allow for an unprecedented global coverage. This makes infrasound as an earthquake mitigation technique feasible for the first time.@en

The Royal Netherlands Meteorological Institute (KNMI) operates a three-dimensional microbarometer array at the Cabauw Experimental Site for Atmospheric Research observatory. The array consists of five microbarometers on a meteorological tower up to an altitude of 200 m. Ten ground-based microbarometers surround the tower with an array aperture of 800 m. This unique setup allows for the study of infrasound propagation in three dimensions. The added value of the vertical dimension is the sensitivity to wind and temperature in the atmospheric boundary layer over multiple altitudes. In this study, we analyze infrasound generated by an accidental chemical explosion at the Moerdijk petrochemical plant on 3 June 2014. The recordings of the tower microbarometers show two sequential arrivals, whereas the recordings on the ground show one wavefront. This arrival structure is interpreted to be the upgoing and downgoing wavefronts. The observations are compared with propagation modeling results using global-scale and mesoscale atmospheric models. Independent temperature and wind measurements, which are available at the Cabauw Experimental Site for Atmospheric Research, are used for comparison with model output. The modeling results explain the signal arrival times; however, the tower wavefront arrivals are not explained. This study is important for understanding the influence of the atmospheric boundary layer on infrasound detections and propagation.@en

The receiver-to-source backazimuth of atmospheric infrasound signals is biased when cross-winds are present along the propagation path. Infrasound from 598 surface explosions from over 30 years in northern Finland is measured with high spatial resolution on an array 178 km almost due North. The array is situated in the classical shadow-zone distance from the explosions. However, strong infrasound is almost always observed, which is most plausibly due to partial reflections from stratospheric altitudes. The most probable propagation paths are subject to both tropospheric and stratospheric cross-winds, and the wave-propagation modelling in this study yields good correspondence between the observed backazimuth deviation and cross-winds from the European Centre for Medium-Range Weather Forecasts Reanalysis (ERA)-Interim reanalysis product. This study demonstrates that atmospheric cross-winds can be estimated directly from infrasound data using propagation time and backazimuth deviation observations. This study finds these cross-wind estimates to be in good agreement with the ERA-Interim reanalysis.

@en