Subduction-zone structure and magmatic processes beneath Costa Rica constrained by local earthquake tomography and petrological modelling

A high-quality data set of 3790 earthquakes were simultaneously inverted for hypocentre locations and 3-D P-wave velocities in Costa Rica. Tests with synthetic data and resolution estimates derived from the resolution matrix indicate that the velocity model is well constrained in central Costa Rica to a depth of 70 km; northwestern and southeastern Costa Rica are less well resolved owing to a lack of seismic stations and seismicity. Maximum H2O content and seismic wave speeds of mid-ocean ridge basalt and harzburgite were calculated for metamorphic phase transformations relevant to subduction. Both the 3-D P-wave velocity structure and petrological modelling indicate the existence of low-velocity hydrous oceanic crust in the subducting Cocos Plate beneath central Costa Rica. Intermediate-depth seismicity correlates well with the predicted locations of hydrous metamorphic rocks, suggesting that dehydration plays a key role in generating intermediate-depth earthquakes beneath Costa Rica. Wadati-Benioff zone seismicity beneath central Costa Rica shows a remarkable decrease in maximum depth toward southeastern Costa Rica. The presence of asthenosphere beneath southeastern Costa Rica, which entered through a proposed slab window, may explain the shallowing of seismicity due to increased temperatures and associated shallowing of dehydration of the slab. Tomographic images further constrain the existence of deeply subducted seamounts beneath central Costa Rica. Large, low P-wave velocity areas within the lower crust are imaged beneath the southeasternmost volcanoes in central Costa Rica. These low velocities may represent anomalously hot material or even melt associated with active volcanism in central Costa Rica. Tomographic images and petrological modelling indicate the existence of a shallow, possibly hydrated mantle wedge beneath central Costa Rica.