Workflow for multi-array beamforming in seismology to detect ocean microseism

The code and data in this repository is linked to the research article: "High-Resolution Spatiotemporal Monitoring of Secondary Microseisms via Multi-Array Analysis, Gao, Y., Rietbrock, A., Tilmann, F., Dushi, E., Geophysical Journal International (2025, under revision) This data repository contains the following files and subdirectories: Beamforming: Hinet: Slowness_data_hourly: beampower data and backprojection distance_baz_slowness_*: distance and back azimuth from the global grid to the center of network assuming seismic phase, e.g. P or PKP Step1_beamforming_hourly.py: beamforming code for the cross-correlation data (data can be obtained on request) Step2_beamforming.ipynb: plot the beamforming from step1 Step3_backprojection.py: backprojection to generate distance_baz_slowness_* ANTICS (SCSN) Slowness_data_hourly: beampower data and backprojection distance_baz_slowness_*: distance and back azimuth from the global grid to the center of network assuming seismic phase, e.g. P or PKP Waveheight: data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_WAV_001_027/services or !motuclient --motu nrt.cmems-du.eu/motu-web/Motu --service-id GLOBAL_ANALYSISFORECAST_WAV_001_027-TDS --product-id cmems_mod_glo_wav_anfc_0.083deg_PT3H-i --date-min 2022-11-18 00:00:00 --date-max 2022-11-1 WWIII: surface spectral density F_5_10s_glob: Equivalent force calculated using WMSAN github.com/lystom/WMSAN Summary_dailyplot_hourly_withwwiii3.py: plot the multi-array beamforming and compare with significant waveheight, surface spectral density, and P wave equivalent force on the ocean floor CCF: Is an empty dir but should store the cross-correlation data in ASDF format, due to the large size of dataset, the data could be obtained on request.