Locating the Nordstream explosions without a velocity model using polarization analysis

The seismic events that preceded the leaks in the Nordstream pipelines in the Baltic Sea have been interpreted as explosions on the seabed, most likely man-made. We use a polarization-based location method initially developed for marsquakes to locate the source region without a subsurface velocity model. We show that the 2 largest seismic events can be unambiguously attributed to the methane plumes observed on the sea surface. The two largest events can be located with this method, using 4 and 5 stations located around the source, with location uncertainties of 30km and 10x60km. We can further show that both events emitted seismic energy for at least ten minutes after the initial explosion, indicative of resonances in the water column or the depressurizing pipeline.Comment: 6 pages, 2 figures, submitted as fast report to Seismic

Expected seismicity and the seismic noise environment of Europa

Seismic data will be a vital geophysical constraint on internal structure of Europa if we land instruments on the surface. Quantifying expected seismic activity on Europa both in terms of large, recognizable signals and ambient background noise is important for understanding dynamics of the moon, as well as interpretation of potential future data. Seismic energy sources will likely include cracking in the ice shell and turbulent motion in the oceans. We define a range of models of seismic activity in Europa's ice shell by assuming each model follows a Gutenberg-Richter relationship with varying parameters. A range of cumulative seismic moment release between $10^{16}$ and $10^{18}$ Nm/yr is defined by scaling tidal dissipation energy to tectonic events on the Earth's moon. Random catalogs are generated and used to create synthetic continuous noise records through numerical wave propagation in thermodynamically self-consistent models of the interior structure of Europa. Spectral characteristics of the noise are calculated by determining probabilistic power spectral densities of the synthetic records. While the range of seismicity models predicts noise levels that vary by 80 dB, we show that most noise estimates are below the self-noise floor of high-frequency geophones, but may be recorded by more sensitive instruments. The largest expected signals exceed background noise by $\sim$50 dB. Noise records may allow for constraints on interior structure through autocorrelation. Models of seismic noise generated by pressure variations at the base of the ice shell due to turbulent motions in the subsurface ocean may also generate observable seismic noise.Comment: 24 pages, 11 figures, Added in supplementary information from revision submission, including 3 audio files with sonification of Europa noise records. To view attachments, please download and extract the gzipped tar source file listed under "Other formats

Bayesian parameter-estimation of Galactic binaries in LISA data with Gaussian Process Regression

The Laser Interferometer Space Antenna (LISA), which is currently under construction, is designed to measure gravitational wave signals in the milli-Hertz frequency band. It is expected that tens of millions of Galactic binaries will be the dominant sources of observed gravitational waves. The Galactic binaries producing signals at mHz frequency range emit quasi monochromatic gravitational waves, which will be constantly measured by LISA. To resolve as many Galactic binaries as possible is a central challenge of the upcoming LISA data set analysis. Although it is estimated that tens of thousands of these overlapping gravitational wave signals are resolvable, and the rest blurs into a galactic foreground noise; extracting tens of thousands of signals using Bayesian approaches is still computationally expensive. We developed a new end-to-end pipeline using Gaussian Process Regression to model the log-likelihood function in order to rapidly compute Bayesian posterior distributions. Using the pipeline we are able to solve the Lisa Data Challange (LDC) 1-3 consisting of noisy data as well as additional challenges with overlapping signals and particularly faint signals.Comment: 12 pages, 10 figure

Accelerating global parameter estimation of gravitational waves from Galactic binaries using a genetic algorithm and GPUs

The Laser Interferometer Space Antenna (LISA) is a planned space-based gravitational wave telescope with the goal of measuring gravitational waves in the milli-Hertz frequency band, which is dominated by millions of Galactic binaries. While some of these binaries produce signals that are loud enough to stand out and be extracted, most of them blur into a confusion foreground. Current methods for analyzing the full frequency band recorded by LISA to extract as many Galactic binaries as possible and to obtain Bayesian posterior distributions for each of the signals are computationally expensive. We introduce a new approach to accelerate the extraction of the best fitting solutions for Galactic binaries across the entire frequency band from data with multiple overlapping signals. Furthermore, we use these best fitting solutions to omit the burn-in stage of a Markov chain Monte Carlo method and to take full advantage of GPU-accelerated signal simulation, allowing us to compute posterior distributions in 2 seconds per signal on a laptop-grade GPU.Comment: 13 pages, 11 figure

First Focal Mechanisms of Marsquakes

Since February 2019, NASA's InSight lander is recording seismic signals on the planet Mars, which, for the first time, allows to observe ongoing tectonic processes with geophysical methods. A number of Marsquakes have been located in the Cerberus Fossae graben system in Elysium Planitia and further west, in the Orcus Patera depression. We present a first study of the focal mechanisms of three well-recorded events (S0173a, S0183a, S0235b) to determine the processes dominating in the source region. We infer for all three events a predominantly extensional setting. Our method is adapted to the case of a single, multicomponent receiver and based on fitting waveforms of P and S waves against synthetic seismograms computed for the initial crustal velocity model derived by the InSight team. We explore the uncertainty due to the single-station limitation and find that even data recorded by one station constrains the mechanisms (reasonably) well. For the events in the Cerberus Fossae region (S0173a, S0235b) normal faulting with a relatively steep dipping fault plane is inferred, suggesting an extensional regime mainly oriented E-W to NE-SW. The fault regime in the Orcus Patera region is not determined uniquely because only the P wave can be used for the source inversion. However, we find that the P and weak S waves of the S0183a event show similar polarities to the event S0173, which indicates similar fault regimes

The Polarization of Ambient Noise on Mars

Seismic noise recorded at the surface of Mars has been monitored since February 2019, using the InSight seismometers. This noise can reach −200 dB. It is 500 times lower than on Earth at night and it increases of 30 dB during the day. We analyze its polarization as a function of time and frequency in the band 0.03–1 Hz. We use the degree of polarization to extract signals with stable polarization independent of their amplitude and type of polarization. We detect polarized signals at all frequencies and all times. Glitches correspond to linear polarized signals which are more abundant during the night. For signals with elliptical polarization, the ellipse is in the horizontal plane below 0.3 Hz. In the 0.3-1Hz high frequency band (HF) and except in the evening, the ellipse is in the vertical plane and the major axis is tilted. While polarization azimuths are different in the two frequency bands, they both vary as a function of local hour and season. They are also correlated with wind direction, particularly during the daytime. We investigate possible aseismic and seismic origins of the polarized signals. Lander or tether noise can be discarded. Pressure fluctuations transported by wind may explain part of the HF polarization but not the tilt of the ellipse. This tilt can be obtained if the source is an acoustic emission coming from high altitude at critical angle. Finally, in the evening when the wind is low, the measured polarized signals may correspond to the seismic wavefield of the Mars background noise