Large-N cross-correlation functions of ambient seismic noise recorded in the Vienna basin, Austria.

This repository contains the C₁ cross-correlation functions of the ambient seismic field and station metadata required to reproduce the results described in the manuscript "Source effects in higher-order ambient seismic field correlations" by Schippkus et al., 2025.

We provide a complete set of Jupyter notebooks that implements all processing described in the manuscript for full reproducibility. They make use of the data provided in this repository. The notebooks are hosted on GitHub at https://github.com/schipp/higher_order_correlations_c2.

The files hosted here are:

• stations_receivers.csv: List of 1990 stations used as receiver stations. Station OMV.GDT is used as the master station in the manuscript. Station names are arbitrary.
• stations_auxiliary.csv: List of 304 stations used as auxiliary stations.
• correlations_for_c1_data.pt: C₁ cross-correlations of all 1990 receiver stations with the master station OMV.GDT in the center. Saved as a PyTorch torch.tensor with shape [1990, 3001]. Sampling rate is 5 Hz, correlation functions are limited to 300 seconds of anti-causal and causal lapse time (=> 3001 samples). First dimension (the receiver stations) is sorted alphabetically by station name.
• correlations_for_c2_data.pt: C₁ cross-correlations of all 1990 receiver stations with the 304 auxiliary stations surrounding them. Saved as a PyTorch torch.tensor with shape [1990, 304, 3001]. Sampling rate is 5 Hz, correlation functions are limited to 300 seconds of anti-causal and causal lapse time (=> 3001 samples). First dimension (the receiver stations) and second dimension (the auxiliary stations) are sorted alphabetically by station name.

The cross-correlations are computed as described in the manuscript: ~4 weeks of continuous recordings are cut into 1-hr windows and spectrally whitened. All windows are cross-correlated and linearly stacked. No additional processing is applied. For more details on these data, please see the manuscript.

We provide a minimal Python code snippet to load the C₁ correlation functions with the master station ("OMV.GDT"), read station locations, filter the correlation functions with a narrowband-filter around 0.3 Hz, and plot the focal spot of the C₁ correlation wavefield in space:

import torch import polars as pl import matplotlib.pyplot as plt from scipy.signal import butter, filtfilt # load correlation functions sampling_rate = 5 lapse_times = torch.arange(-300, 300 + 1 / sampling_rate, 1 / sampling_rate) corrs = torch.load("correlations_for_c1_data.pt", weights_only=False) # load metadata stations = pl.read_csv("stations_receivers.csv") names = stations.select("station").to_series().to_list() x, y = stations.select(["X", "Y"]).to_numpy().T master_station = "OMV.GDT" # apply acausal narrowband filter frequency_band = (0.29, 0.31) a, b = butter(4, frequency_band, btype="bandpass", fs=sampling_rate) filtered_corrs = torch.tensor(filtfilt(a, b, corrs, axis=-1).copy()) # extract focal spot and set master station (=auto-correlation) to zero focal_spot = filtered_corrs[:, torch.argmin(torch.abs(lapse_times))] focal_spot[names.index(master_station)] = 0 # plot focal spot fig, ax = plt.subplots() ax.scatter( x, y, c=focal_spot, cmap="RdBu_r", vmin=-focal_spot.abs().max(), vmax=focal_spot.abs().max(), ) ax.set( xlabel="Distance (km)", ylabel="Distance (km)", title="Focal spot", aspect="equal" )