scholarly journals Introducing noisi: a Python tool for ambient noise cross-correlation modeling and noise source inversion

Solid Earth ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 1597-1615
Author(s):  
Laura Ermert ◽  
Jonas Igel ◽  
Korbinian Sager ◽  
Eléonore Stutzmann ◽  
Tarje Nissen-Meyer ◽  
...  

Abstract. We introduce the open-source tool noisi for the forward and inverse modeling of ambient seismic cross-correlations with spatially varying source spectra. It utilizes pre-computed databases of Green's functions to represent seismic wave propagation between ambient seismic sources and seismic receivers, which can be obtained from existing repositories or imported from the output of wave propagation solvers. The tool was built with the aim of studying ambient seismic sources while accounting for realistic wave propagation effects. Furthermore, it may be used to guide the interpretation of ambient seismic auto- and cross-correlations, which have become preeminent seismological observables, in light of nonuniform ambient seismic sources. Written in the Python language, it is accessible for both usage and further development and efficient enough to conduct ambient seismic source inversions for realistic scenarios. Here, we introduce the concept and implementation of the tool, compare its model output to cross-correlations computed with SPECFEM3D_globe, and demonstrate its capabilities on selected use cases: a comparison of observed cross-correlations of the Earth's hum to a forward model based on hum sources from oceanographic models and a synthetic noise source inversion using full waveforms and signal energy asymmetry.

2020 ◽  
Author(s):  
Laura Ermert ◽  
Jonas Igel ◽  
Korbinian Sager ◽  
Eléonore Stutzmann ◽  
Tarje Nissen-Meyer ◽  
...  

Abstract. We introduce open-source tool noisi for the forward and inverse modeling of ambient seismic cross-correlations with spatially varying source spectra. It utilizes pre-computed databases of Green’s functions to represent seismic wave propagation between ambient seismic sources and seismic receivers, which can be obtained from existing repositories or imported from the output of wave propagation solvers. The tool was built with the aim of studying ambient seismic sources while accounting for realistic wave propagation effects. Furthermore, it may be used to guide the interpretation of ambient seismic auto- and cross-correlations, which have become pre-eminent seismological observables, in light of non-uniform ambient seismic sources. Written in the Python language, it is both accessible for usage and further development, as well as efficient enough to conduct ambient seismic source inversions for realistic scenarios. Here, we introduce the concept and implementation of the tool, compare its model output to the output of cross-correlations computed with SPECFEM3D_globe, and demonstrate its capabilities on selected use cases: A comparison of observed cross-correlations of the Earth’s hum to a forward model based on hum sources from oceanographic models, and a synthetic noise source inversion using full waveforms and signal energy asymmetry.


2020 ◽  
Author(s):  
Jonas Igel ◽  
Laura Ermert ◽  
Andreas Fichtner

<p>Common assumptions in ambient noise seismology such as Green’s function retrieval and equipartitioned wavefields are often not met in the Earth. Full waveform ambient noise tomography methods are free of such assumptions, as they implement knowledge of the time- and space-dependent ambient noise source distribution, whilst also taking finite-frequency effects into account. Such methods would greatly simplify near real-time monitoring of the sub-surface. Additionally, the distribution of the secondary microseisms could act as a new observable of the ocean state since its mechanism is well understood (e.g. Ardhuin et al., 2011).</p><p>To efficiently forward-model global noise cross-correlations we implement (1) pre-computed high-frequency wavefields obtained using, for example, AxiSEM (Nissen-Meyer et al., 2014), and (2) spatially variable grids, both of which greatly reduce the computational cost. Global cross-correlations for any source distribution can be computed within a few seconds in the microseismic frequency range (up to 0.2 Hz). Similarly, we can compute the finite-frequency sensitivity kernels which are then used to perform a gradient-based iterative inversion of the power-spectral density of the noise source distribution. We take a windowed logarithmic energy ratio of the causal and acausal branches of the cross-correlations as measurement, which is largely insensitive to unknown 3D Earth structures.</p><p>Due to its parallelisation on a cluster, our inversion tool is able to rapidly invert for the global microseismic noise source distribution with minimal required user interaction. Synthetic and real data inversions show promising results for noise sources in the North Atlantic with the structure and spatial distribution resolved at scales of a few hundred kilometres. Finally, daily noise sources maps could be computed by combining our inversion tool with a daily data download and processing toolkit.</p>


2020 ◽  
Author(s):  
Foivos Karakostas ◽  
Nicholas Schmerr ◽  
Samuel Hop Bailey ◽  
Daniella Dellagiustina ◽  
Namrah Habib ◽  
...  

<p>On July 25, 2018, a meteoroid-associated airburst occurred near the Qaanaaq town, in Greenland, at approximately 22:00 UTC (20:00 local time). The event generated seismic waves that were recorded by two stations of the Danish Seismological Network (TULEG and NEEM) and the bolide trajectory was consequently calculated by the NASA Center for Near-Earth Object Studies (CNEOS). The total impact energy, calculated by CNEOS was 2.1 kT of TNT and the brightest point on its trajectory corresponds to an altitude of around 43 km, at a distance of about 50 km S of the Qaanaaq town and 50 km N of the TULEG station and the Thule Air Force Base [1].</p> <p>An airburst occurring over the icy surface of Greenland is a rare terrestrial analog for regions of the Solar System, where both an atmosphere and an icy surface exist. In the past, a variety of works had indicated the presence of ice on Titan, the biggest moon of Saturn (e.g. [2] and more recently [3]) and more precisely, the icy composition of mountains which are formed by tectonic activity [4]. Titan has a relatively thick atmosphere, compared to those of other moons in the Solar System, composed mainly (94%) of nitrogen [5]. The characterization of atmospheric meteoroid-associated seismic sources for Titan has a particular interest, as it is found that, contrary to other moons of the solar system, the presence of craters on its surface is extremely low (only about 0.4% according to [3]). The reason for this low cratering of the surface is the presence of the thick atmosphere, into which many of the meteoroids are entirely ablated into dust. Therefore, a methodology for the characterization of airbursts as seismic sources and the modeling of the associated generated seismic waves is necessary for a future seismic experiment, as any recorded signal will either be a direct atmospheric wave (nonlinear shock wave, or linear acoustic wave) or a seismic wave generated through the coupling of the atmospheric and solid/ice part.<br /> <br />In the present study, our aim is to perform a seismic investigation of the Greenland ice shell with the use of the airburst-associated seismic source. The performed tasks into which this effort has been divided, include the application of a technique which approaches the bolide as an atmospheric seismic source, the calculation of the distance of shock wave propagation in the atmosphere, the description of the mechanism of generation of the seismic waves in the atmosphere and the solid-icy part.</p> <p>When the bolides enter the atmosphere of the Earth or that of any other body, shock waves are generated along the trajectory of the meteoroid. These waves are characterized by the overpressure that they generate, which create a clear pressure discontinuity in the atmosphere, referred to as the nonlinear part of the shock wave propagation. The propagation distance of this nonlinear wave is associated to the ratio of the meteoroid speed to the ambient sound speed, also known as the Mach number, as well as the physical diameter of the meteoroid. In this work, we compute this distance for the Earth case and for the known trajectory of the detected and examined bolide [1][6].</p> <p>The methodology developed in this study can serve the seismic investigation of structures covered by ice on planets or planetary bodies with a relatively thick atmosphere, where airbursts can occur due to the friction of the meteoroid with the ambient atmospheric material. An ideal example of this case are the icy mountains of Titan, which are known to be formed by tectonic activity on the Saturn’s moon [4]. The future Dragonfly mission to Titan will carry a seismometer as part of the DraGMet (Dragonfly Geophysics and Meteorology Package) payload [7]. Even if the primary goal of the mission is the characterization of the regolith properties, an eventual airburst and collection of seismic data near these mountainous icy structures, will be a great opportunity to investigate, through the identification of the associated waves and thus the investigation of the coupled seismic waves, the properties of this icy cover, its depth and composition.</p> <p>References: [1] https://cneos.jpl.nasa.gov/fireballs/ [2] Sohl, F. et al. (1995) Icarus, 115, 278–294 [3] Lopes R.M.C. et al. (2019) Nat Astron, [4] Radebaugh J. et al. (2007) Icarus, 192, 77-91, [5] Niemann H.B. et al. (2005) Nature, 438, 779–784 [6] Schmerr, N. et al. (2018) Abstract P21E-3406, AGU Fall Meeting 2018, Washington DC [7] Lorenz R. et al. (2018) Johns Hop- kins APL Technical Digest, 34, 3</p>


2020 ◽  
Author(s):  
Laura Ermert ◽  
Jonas Igel ◽  
Korbinian Sager ◽  
Eléonore Stutzmann ◽  
Tarje Nissen-Meyer ◽  
...  

2020 ◽  
Vol 53 (11) ◽  
pp. 5177-5201
Author(s):  
Shengquan He ◽  
Tuo Chen ◽  
Isaac Vennes ◽  
Xueqiu He ◽  
Dazhao Song ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document