Search

We have integrated waveform and arrival-onset data collected in Costa Rica as part of the National Science Foundation (NSF)-sponsored Costa Rica Seismogenic Zone Experiment (CRSEIZE) and along central Costa Rica and Nicaragua as part of the German SFB 574 program. The five arrays, composed of different sensor types (one-and three-component land and ocean bottom seismometers and hydrophones), were archived using different software packages (Antelope and SEISAN) and were automatically and manually picked using various quality criteria resulting in a disparate set of pick weights. We evaluate pick quality using automated arrival detection and picking algorithm based on the wavelet transform and Akaike information criterion picker. The consistency of the arrival information over various scales provides a basis for assigning a quality to the analyst pick. Approximately 31% of P arrival times and 26% of S times have been classified as high-quality picks (quality 0-1). An additional 21% of P times and 27% of S arrivals are good quality (quality 2-3). The revised quality picks are mapped directly into new pick weights for inversion studies. We explore the effect of new weighting and removal of poor data by relocating hypocenters through a minimum 1D velocity model and conducting double-difference local earthquake tomography (LET). Analysis of the hypocenter relocation and seismic velocity tomography results suggest that using the improved quality determinations have a greater effect on improving sharpness in the velocity images than on the magnitude of hypocentral movement.

We present a new tomographic model of the mantle in the area of the 2010 M8.8 Maule earthquake and surrounding regions. Increased ray coverage provided by the aftershock data allows us to image the detailed subducting slab structure in the mantle, from the region of flat slab subduction north of the Maule rupture to the area of overlapping rupture between the 1960 M9.5 and the 2010 M8.8 events to the south. We have combined teleseismic primary and depth phase arrivals with available local arrivals to better constrain the teleseismic earthquake locations in the region, which we use to conduct nested regional–global tomography. The new model reveals the detailed structure of the flat slab and its transition to a more moderately dipping slab in the Maule region. South of the Maule region, a steeply dipping relic slab is imaged from ∼200 to 1000 km depth that is distinct from the moderately dipping slab above it and from the more northerly slab at similar depths. We interpret the images as revealing both horizontal and vertical tearing of the slab at ∼38°S to explain the imaged pattern of slab anomalies in the southern portion of the model. In contrast, the transition from a horizontal to moderately subducting slab in the northern portion of the model is imaged as a continuous slab bend. We speculate that the tearing was most likely facilitated by a fracture zone in the downgoing plate or alternatively by a continental scale terrane boundary in the overriding plate.

The subduction plate interface along the Nicoya Peninsula, Costa Rica, generates damaging large (Mw > 7.5) earthquakes. We present hypocenters and 3-D seismic velocity models (VP and VP/VS) calculated using simultaneous inversion of P- and S-wave arrival time data recorded from small magnitude, local earthquakes to elucidate seismogenic zone structure. In this region, interseismic cycle microseismicity does not uniquely define the potential rupture extent of large earthquakes. Plate interface microseismicity extends from 12 to 26 and from 17 to 28 km below sea level beneath the southern and northern Nicoya Peninsula, respectively. Microseismicity offset across the plate suture of East Pacific Rise-derived and Cocos-Nazca Spreading Center-derived oceanic lithosphere is ∼5 km, revising earlier estimates suggesting ∼10 km of offset. Interplate seismicity begins downdip of increased locking along the plate interface imaged using GPS and a region of low VP along the plate interface. The downdip edge of plate interface microseismicity occurs updip of the oceanic slab and continental Moho intersection, possibly due to the onset of ductile behaviour. Slow forearc mantle wedge P-wave velocities suggest 20–30 per cent serpentinization across the Nicoya Peninsula region while calculated VP/VS values suggest 0–10 per cent serpentinization. Interpretation of VP/VS resolution at depth is complicated however due to ray path distribution. We posit that the forearc mantle wedge is regionally serpentinized but may still be able to sustain rupture during the largest seismogenic zone earthquakes.