Applying Automatic Mapping Processing By GMT to Bathymetric and Geophysical Data: Cascadia Subduction Zone, Pacific Ocean

The Cascadia Trench is stretching along the convergent plate boundaries of Pacific Plate, North America Plate and Juan De Fuca Plate. It is an important geomorphological structural feature in the north-east Pacific Ocean. The aim of the paper is to analyse the geomorphology of the Cascadia Trench west of Vancouver Island (Canada and USA) using the GMT cartographic scripting toolset. The unique geomorphological feature of the Cascadia Trench is that the thick sediment layer completely obscures the subduction zone and abyssal hills. This results in the asymmetric profile in the cross-section of the trench. Bathymetric data were extracted from the GEBCO 2019 dataset (15 arc-second grid), sediment thickness by the GlobSed dataset. Due to the dominance of high sedimentary rate and complexity of the tectonic processes and geologic settings, Cascadia Trench develops very specific asymmetric geomorphic shape comparing to the typical V-form. The results of the geomorphic modelling show that eastern side of the trench has a gentle curvature (slope: 35.12°), partially stepped, due to the tectonic movements and faults. The opposite, oceanward side is almost completely leveled. The trench is narrow with maximal depth at the selected segment -3489 m and for the whole dataset -6201 m. The most repetitive depth is in a range -2500 to -2400 m (267 samples) and -2500 to -2600 m (261 samples). The bottom is mostly flat due to the high sedimentation rates indicating the accumulative leveling processes. Marine free-air gravity anomalies along the Cascadia Subduction Zone are characterized by weakly positive values (20 mGal) increasing rapidly in the zone of the continental slope (>200 mGal), which is associated with a decrease in thickness of the Earth’s crust.

Amphibious surface-wave phase-velocity measurements of the Cascadia subduction zone

SUMMARY A new amphibious seismic data set from the Cascadia subduction zone is used to characterize the lithosphere structure from the Juan de Fuca ridge to the Cascades backarc. These seismic data are allowing the imaging of an entire tectonic plate from its creation at the ridge through the onset of the subduction to beyond the volcanic arc, along the entire strike of the Cascadia subduction zone. We develop a tilt and compliance correction procedure for ocean-bottom seismometers that employs automated quality control to calculate robust station noise properties. To elucidate crust and upper-mantle structure, we present shoreline-crossing Rayleigh-wave phase-velocity maps for the Cascadia subduction zone, calculated from earthquake data from 20 to 160 s period and from ambient-noise correlations from 9 to 20 s period. We interpret the phase-velocity maps in terms of the tectonics associated with the Juan de Fuca plate history and the Cascadia subduction system. We find that thermal oceanic plate cooling models cannot explain velocity anomalies observed beneath the Juan de Fuca plate. Instead, they may be explained by a ≤1 per cent partial melt region beneath the ridge and are spatially collocated with patches of hydration and increased faulting in the crust and upper mantle near the deformation front. In the forearc, slow velocities appear to be more prevalent in areas that experienced high slip in past Cascadia megathrust earthquakes and generally occur updip of the highest-density tremor regions and locations of intraplate earthquakes. Beneath the volcanic arc, the slowest phase velocities correlate with regions of highest magma production volume.