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

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

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

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

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

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

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

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

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

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

This chapter examines the potential improvements in tropospheric weather forecasts that might arise from an enhanced representation of the upper stratospheric state. First, the chapter reviews current operational practice regarding observation of the atmosphere and the relative paucity of observations in the altitude range 40–70 km. Then, we describe some idealised model calculations to quantify the potential gain in skill available from improved monitoring in this region. The idealised model experiments use a relaxation technique with the Hadley Centre General Environment Model, to assess the potential gain in skill from observations both of the whole stratosphere and the upper stratosphere. At weather forecasting timescales (up to forecast day 30), better knowledge of the stratosphere, close to the onset of a sudden stratospheric warming, improves forecasts of the tropospheric northern annular mode. Whole-stratosphere information significantly improved average surface temperature anomalies over northern North America, whilst upper stratosphere information improved anomalies over Central Siberia. These results suggest any new observational technique which can contribute to monitoring of the 40–70 km region would likely benefit tropospheric forecast skill during wintertime.@en

Infrasound recordings can be used as input to inversion procedures to delineate the vertical structure of temperature and wind in a range of altitudes where ground-based or satellite measurements are rare and where fine-scale atmospheric structures are not resolved by the current atmospheric specifications. As infrasound is measured worldwide, this allows for a remote sensing technique that can be applied globally. This chapter provides an overview of recently developed infrasonic remote sensing methods. The methods range from linearized inversionsInversions to direct search methods as well as interferometric techniques for atmospheric infrasound. The evaluation of numerical weather prediction (NWP) products shows the added value of infrasound, e.g., during sudden stratospheric warming (SSW) and equinox periods. The potential transition toward assimilation of infrasound in numerical weather prediction models is discussed.@en

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