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:
https://data.marine.copernicus.eu/product/GLOBAL_ANALYSISFORECAST_WAV_001_027/services
or
!motuclient --motu https://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
https://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.