High-frequency passive surface wave tomography at the Marathon PGE-Cu
deposit (Ontario, Canada): Bridging the gap between natural and
anthropogenic noise sources
Abstract
Ambient noise surface wave tomography is an environmentally friendly and
cost-effective seismic technique for subsurface imaging. However, noise
sources acting from preferential azimuths may introduce bias in the
Green’s function reconstruction and in the resultant velocity models.
This study, focused at the deposit scale, investigates how to correctly
merge the different phase velocity measurements at various frequencies,
in order to fill the gap between natural and anthropogenic noise sources
while adjusting the bias caused by changes in the azimuth of the source.
The target is the Marathon PGE-Cu deposit (Ontario, Canada), an alkaline
intrusion containing gabbros and syenites (ø = 25 km). Mineralisation is
hosted by gabbros close to the inward-dipping footwall of the intrusion.
The country rocks are Archaean volcanic breccias. 1024
vertical-component receivers were deployed for 30 days in two
overlapping grids: a 200 m x 6040 m dense array with node spacing of 50
m, and a 4000 m x 2500 m sparse array with node spacing of 150 m.
Beamforming analysis of the recorded data indicates variations in the
distribution of noise. Below 5 Hz, the Lake Superior (SSW) is the
dominant source of noise, while above 12 Hz, noise from the Canadian
Pacific Railway and Trans-Canada highway (SW) is prominent. In the 5 -
12 Hz frequency band, surface-wave energy is dominant, and it comes from
the Lake Superior and vehicle traffic. Between 12 Hz and 20 Hz, the
signal is characterized by body-wave energy combined with less energetic
surface waves, while above 20 Hz the imprint of body waves is dominant.
We retrieved the fundamental mode of Rayleigh wave propagation from the
recorded data set. The signal was down-sampled to 50 Hz, divided into
segments of 30 minutes, cross-correlated and stacked. Surface wave
dispersion curves were extracted from 2-km-long arrays. Besides, various
phase velocity measurements were applied. Phase-velocities were inverted
to S-wave velocity structures using different probabilistic approaches.
The overall results show a high-velocity shallow anomaly, probably
related to the gabbro intrusion hosting the mineralization, as well as
other structures consistent with the geological model inferred from
surface mapping and drill logs.