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

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

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

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

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

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

Instrumental timing and phase errors are a notorious problem in seismic data acquisition and processing. These can be frequency independent, for example due to clock drift, but may also be frequency dependent, for example due to imperfectly known instrument responses. A technique is presented that allows both types of errors to be recovered in a systematic fashion. The methodology relies on the time-symmetry usually inherent in time-averaged cross-correlations of ambient seismic noise: the difference between the arrival time of the direct surface-wave at positive time and the arrival time of the direct surface-wave at negative time is quantified. Doing this for all eligible receiver-receiver pairs of a large-N seismic array, including one or more receivers devoid of instrumental timing errors, the instrumental timing errors of all incorrectly timed receivers can be determined uniquely. Most notably, this is accomplished by means of a weighted least-squares inversion. The weights are based on the receiver-receiver distances and decrease the adverse effect of inhomogeneities in the noise illumination pattern on the recovered instrumental timing errors. Inversion results are furthermore optimized by limiting the inversion to receiver couples that (i) exceed a specific receiver-receiver distance threshold and (ii) whose time-averaged cross-correlations exceed a specific signal-to-noise ratio threshold. Potential frequency dependence of the timing errors is incorporated by means of an iterative, frequency-dependent approach. The proposed methodology is validated using synthetic recordings of ambient seismic surface-wave noise due to an arbitrary non-uniform illumination pattern. The methodology is successfully applied to time-averaged cross-correlations of field recordings of ambient seismic noise on and around the Reykjanes peninsula, SW Iceland.

@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

Continuous long-term sound sources are recorded at hydroacoustic station H03S, a threeelement hydrophone array south of Robinson Crusoe Island between 2014 April 23 and 2017 August 20. The origin of the signal between 3 and 20 Hz is investigated by using cross-correlation, array processing using plane wave beamforming and spectral analysis. Onebit normalization is successfully applied as a cross-correlation pre-processing step in order to suppress undesired earthquake events in the data. Traveltimes retrieved from averaged cross-correlations do not yield a coherent array direction of arrival. Averaged envelopes of the cross-correlations, however, indicate a coherent signal approaching H03S from a south. southwest direction. Beamforming indicates two dominant backazimuth directions: 172°-224° (Antarctica) and 242° (Monowai Volcanic Seamount). This continuous source field creates possibilities to investigate the applicability of acoustic thermometry at hydrophones H03 S1. S2. Cross-correlation and array processing indicate significant directional variation in local propagation, most likely related to the steep slope in the bathymetry near H03S. In addition, it is demonstrated that the ambient noise field is not sufficiently equipartitioned. It is shown that this causes a large error in the estimated temperature, primarily due to the short receiver spacing. These large errors have not been addressed in previous studies on deep-ocean acoustic thermometry. Hence, it is shown that acoustic thermometry does not perform well on small arrays such as H03S. The power spectral density yields a strong broadband signal in January. March, most likely related to iceberg noise. A narrow banded signal around 17 Hz between April and September corresponds to whale calls. The best-beam sound pressure levels towards Antarctica are compared to ERA5 climatologies for sea ice cover and normalized stress into the ocean, supporting the hypothesis of iceberg noise.

@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

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

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

Supersonic rockets generate low-frequency acoustic waves, that is, infrasound, during the launch and re-entry. Infrasound is routinely observed at infrasound arrays from the International Monitoring System, in place for the verification of the Comprehensive Nuclear-Test-Ban Treaty. Association and source identification are key elements of the verification system. The moving nature of a rocket is a defining criterion in order to distinguish it from an isolated explosion. Here, it is shown how infrasound recordings can be associated, which leads to identification of the rocket. Propagation modelling is included to further constrain the source identification. Four rocket launches by the Democratic People's Republic of Korea in 2009 and 2017 are analysed in which multiple arrays detected the infrasound. Source identification in this region is important for verification purposes. It is concluded that with a passive monitoring technique such as infrasound, characteristics can be remotely obtained on sources of interest, that is, infrasonic intelligence, over 4500+ km.@en

Low-frequency acoustic, i.e., infrasound, waves are measured by sparse arrays of microbarometers. Recorded data are processed by automatic detection algorithms based on array-processing techniques such as time-domain beam forming and f-k analysis. These algorithms use a signal-to-noise ratio (S/N) value as a detection criterion. In the case of high background noise or in the presence of multiple coinciding signals, the event's S/N decreases and can be missed by automatic processing. In seismology, detecting low-S/N events with geophone arrays is a well-known problem. Whether it is in global earthquake monitoring or reservoir microseismic activity characterization, detecting low-S/N events is needed to better understand the sources or the medium of propagation. We use an image-processing technique as a postprocessing step in the automatic detection of low S/N events. In particular, we consider the use of the Hough transform (HT) technique to detect straight lines in beam-forming results, i.e., a back azimuth (BA) time series. The presence of such lines, due to similar BA values, can be indicative of a low-S/N event. A statistical framework is developed for the HT parameterization, which includes defining a threshold value for detection as well as evaluating the false alarm rate. The method is tested on synthetic data and five years of recorded infrasound from glaciers. It is shown that the automatic detection capability is increased by detecting low-S/N events while keeping a low false-alarm rate.@en

The 2017 North Korean nuclear test gave rise to seismic and low-frequency acoustic signals, that is, infrasound. The infrasonic signals are due to seismo-acoustic coupling and have been detected on microbarometer array I45RU in the Russian Federation at 401 km from the test site. I45RU is part of the International Monitoring System for the verification of the Comprehensive Nuclear-Test-Ban Treaty. We analyze the seismo-acoustic coupling by making use of array-processing and backprojection techniques. The backprojections show that infrasound radiation is not confined to the epicentral region. More distant regions are found to be consistent with locations of topography, sedimentary basins, and underwater evanescent sources. The backprojections can be used to estimate the average infrasonic propagation speed through the atmosphere. We discuss these findings in the context of infrasound propagation conditions during the sixth nuclear test. It is suggested that propagation from the test site to I45RU may have occurred along unexpected paths instead of typical stratospheric propagation. We present several scenarios that could be considered in the interpretation of the observations. Electronic Supplement: Details on signal characterization and infrasound propagation conditions.

@en