Middle to Late Miocene Contractional Deformation in Costa Rica Triggered by Plate Geodynamics

Contractional deformation in Costa Rica is usually attributed to the subduction of the aseismic Cocos Ridge. In this work, we review the evidences for contraction in the middle to late Miocene, prior to the arrival of the Cocos Ridge at the Middle America Trench. We find that the Miocene phase of contractional deformation is found in all of Costa Rica, probably extending to Nicaragua as well. The widespread distribution of this event requires a regional or plate geodynamic trigger. We analyze the possible mechanisms that could produce the onset of contractional deformation, using the better known case of subduction orogeny, the Andes, as an analog. We propose that a change in the direction of the Cocos plate since ∼19 Ma led to a change from oblique to orthogonal convergence, producing contractional deformation of the upper plate.Fil: Mescua, Jose Francisco. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; Argentina. Universidad Nacional de Cuyo. Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Porras, Hernan. Observatorio Vulcanológico y Sismológico de Costa Rica; Costa RicaFil: Duran, Patrick. Universidad de Costa Rica; Costa RicaFil: Giambiagi, Laura Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Provincia de Mendoza. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales. Universidad Nacional de Cuyo. Instituto Argentino de Nivología, Glaciología y Ciencias Ambientales; ArgentinaFil: de Moor, Maerten. Observatorio Vulcanológico y Sismológico de Costa Rica; Costa RicaFil: Cascante, Monserrat. Observatorio Vulcanológico y Sismológico de Costa Rica; Costa RicaFil: Salazar, Esteban. Servicio Nacional de Geología y Minería; ChileFil: Protti, Marino. Observatorio Vulcanológico y Sismológico de Costa Rica; Costa RicaFil: Poblete, Fernando. Universidad de O’Higgins; Chil

Interplate seismicity at the CRISP drilling site: The 2002 Mw=6.4 Osa Earthquake at the southeastern end of the Middle America Trench

We investigate potential relations between variations in seafloor relief and age of the incoming plate and interplate seismicity. Westward from Osa Peninsula in Costa Rica, a major change in the character of the incoming Cocos Plate is displayed by abrupt lateral variations in seafloor depth and thermal structure. Here a Mw 6.4 thrust earthquake was followed by three aftershock clusters in June 2002. Initial relocations indicate that the main shock occurred fairly trenchward of most large earthquakes along the Middle America Trench off central Costa Rica. The earthquake sequence occurred while a temporary network of OBH and land stations ∼80 km to the northwest were deployed. By adding readings from permanent local stations, we obtain uncommon P wave coverage of a large subduction zone earthquake. We relocate this catalog using a nonlinear probabilistic approach within both, a 1-D and a 3-D P wave velocity models. The main shock occurred ∼25 km from the trench and probably along the plate interface at 5–10 km depth. We analyze teleseismic data to further constrain the rupture process of the main shock. The best depth estimates indicate that most of the seismic energy was radiated at shallow depth below the continental slope, supporting the nucleation of the Osa earthquake at ∼6 km depth. The location and depth coincide with the plate boundary imaged in prestack depth-migrated reflection lines shot near the nucleation area. Aftershocks propagated downdip to the area of a 1999 Mw 6.9 sequence and partially overlapped it. The results indicate that underthrusting of the young and buoyant Cocos Ridge has created conditions for interplate seismogenesis shallower and closer to the trench axis than elsewhere along the central Costa Rica margin

Late Cenozoic tephrostratigraphy offshore the southern Central American Volcanic Arc: 2. Implications for magma production rates and subduction erosion

Pacific drill sites offshore Central America provide the unique opportunity to study the evolution of large explosive volcanism and the geotectonic evolution of the continental margin back into the Neogene. The temporal distribution of tephra layers established by tephrochonostratigraphy in Part 1 indicates a nearly continuous highly explosive eruption record for the Costa Rican and the Nicaraguan volcanic arc within the last 8 M.y. The widely distributed marine tephra layers comprise the major fraction of the respective erupted tephra volumes and masses thus providing insights into regional and temporal variations of large-magnitude explosive eruptions along the southern Central American Volcanic Arc (CAVA). We observe three pulses of enhanced explosive magmatism between 0-1 Ma at the Cordillera Central, between 1-2 Ma at the Guanacaste and at >3 Ma at the Western Nicaragua segments. Averaged over the long-term the minimum erupted magma flux (per unit arc length) is ∼0.017 g/ms. Tephra ages, constrained by Ar-Ar dating and by correlation with dated terrestrial tephras, yield time-variable accumulation rates of the intercalated pelagic sediments with four prominent phases of peak sedimentation rates that relate to tectonic processes of subduction erosion. The peak rate at >2.3 Ma near Osa particularly relates to initial Cocos Ridge subduction which began at 2.91±0.23 Ma as inferred by the 1.5 M.y. delayed appearance of the OIB geochemical signal in tephras from Barva volcano at 1.42 Ma. Subsequent tectonic re-arrangements probably involved crustal extension on the Guanacaste segment that favored the 2-1 Ma period of unusually massive rhyolite production

The steeply subducting edge of the Cocos Ridge : evidence from receiver functions beneath the northern Talamanca Range, south-central Costa Rica

The deep structure of the south-central Costa Rican subduction zone has not been studied in great detail so far because large parts of the area are virtually inaccessible. We present a receiver function study along a transect of broadband seismometers through the northern flank of the Cordillera de Talamanca (south Costa Rica). Below Moho depths, the receiver functions image a dipping positive conversion signal. This is interpreted as the subducting Cocos Plate slab, compatible with the conversions in the individual receiver functions. In finite difference modeling, a dipping signal such as the one imaged can only be reproduced by a steeply (80°) dipping structure present at least until a depth of about 70–100 km; below this depth, the length of the slab cannot be determined because of possible scattering effects. The proposed position of the slab agrees with previous results from local seismicity, local earthquake tomography, and active seismic studies, while extending the slab location to greater depths and steeper dip angle. Along the trench, no marked change is observed in the receiver functions, suggesting that the steeply dipping slab continues until the northern flank of the Cordillera de Talamanca, in the transition region between the incoming seamount segment and Cocos Ridge. Considering the time predicted for the establishment of shallow angle underthrusting after the onset of ridge collision, the southern Costa Rican subduction zone may at present be undergoing a reconfiguration of subduction style, where the transition to shallow underthrusting may be underway but still incomplete

Ridge subduction at an erosive margin - the collision zone of the Nazca Ridge in southern Peru

The 1.5-km-high, obliquely subducting Nazca Ridge and its collision zone with the Peruvian margin have been imaged by wide-angle and reflection seismic profiles, swath bathymetry, and gravity surveying. These data reveal that the crust of the ridge at its northeastern tip is 17 km thick and exhibits seismic velocities and densities similar to layers 2 and 3 of typical oceanic crust. The lowermost layer contributes 10–12 km to the total crustal thickness of the ridge. The sedimentary cover is 300–400 m thick on most parts of the ridge but less than 100 m thick on seamounts and small volcanic ridges. At the collision zone of ridge and margin, the following observations indicate intense tectonic erosion related to the passage of the ridge. The thin sediment layer on the ridge is completely subducted. The lower continental slope is steep, dipping at ∼9°, and the continental wedge has a high taper of 18°. Tentative correlation of model layers with stratigraphy derived from Ocean Drilling Program Leg 112 cores suggests the presence of Eocene shelf deposits near the trench. Continental basement is located <15 km landward of the trench. Normal faults on the upper slope and shelf indicate extension. A comparison with the Peruvian and northern Chilean forearc systems, currently not affected by ridge subduction, suggests that the passage of the Nazca Ridge along the continental margin induces a temporarily limited phase of enhanced tectonic erosion superposed on a long-term erosive regime

Geodetic and seismic constraints on some seismogenic zone processes in Costa Rica

New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costeña fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence

Fast rates of subduction erosion along the Costa Rica Pacific margin: implications for non-steady rates of crustal recycling at subduction zones

At least since the middle Miocene (∼16 Ma), subduction erosion has been the dominant process controlling the tectonic evolution of the Pacific margin of Costa Rica. Ocean Drilling Program Site 1042 recovered 16.5 Ma nearshore sediment at ∼3.9 km depth, ∼7 km landward of the trench axis. The overlying Miocene to Quaternary sediment contains benthic foraminifera documenting margin subsidence from upper bathyal (∼200 m) to abyssal (∼2000 m) depth. The rate of subsidence was low during the early to middle Miocene but increased sharply in the late Miocene-early Pliocene (5–6.5 Ma) and at the Pliocene-Pleistocene boundary (2.4 Ma). Foraminifera data, bedding dip, and the geometry of slope sediment indicate that tilting of the forearc occurred coincident with the onset of rapid late Miocene subsidence. Seismic images show that normal faulting is widespread across the continental slope; however, extension by faulting only accounts for a minor amount of the post-6.5 Ma subsidence. Basal tectonic erosion is invoked to explain the subsidence. The short-term rate of removal of rock from the forearc is about 107–123 km3 Myr−1 km−1. Mass removal is a nonsteady state process affecting the chemical balance of the arc: the ocean sediment input, with the short-term erosion rate, is a factor of 10 smaller than the eroded mass input. The low 10Be concentration in the volcanic arc of Costa Rica could be explained by dilution with eroded material. The late Miocene onset of rapid subsidence is coeval with the arrival of the Cocos Ridge at the subduction zone. The underthrusting of thick and thermally younger ocean crust decreased the subduction angle of the slab along a large segment of the margin and changed the dynamic equilibrium of the margin taper. This process may have induced the increase in the rate of subduction erosion and thus the recycling of crustal material to the mantle

Seismic Structure of Cocos and Malpelo Volcanic Ridges and Implications for Hotspot-Ridge Interaction

The Cocos and Malpelo Volcanic Ridges are blocks of thickened oceanic crust thought to be the result of the interaction between the Galapagos hot spot and the Cocos‐Nazca Spreading Center during the last 20 m.y. In this work we investigate the seismic structure of these two aseismic ridges along three wide‐angle transects acquired during the Panama basin and Galapagos plume—New Investigations of Intraplate magmatism (PAGANINI)‐1999 experiment. A two‐dimensional velocity field with the Moho geometry is obtained using joint refraction/reflection travel time tomography, and the uncertainty and robustness of the results are estimated by performing a Monte Carlo‐type analysis. Our results show that the maximum crustal thickness along these profiles ranges from ∼16.5 km (southern Cocos) to ∼19 km (northern Cocos and Malpelo). Oceanic layer 2 thickness is quite uniform regardless of total crustal thickness variations; crustal thickening is mainly accommodated by layer 3. These observations are shown to be consistent with gravity data. The variation of layer 3 velocities is similar along all profiles, being lower where crust is thicker. This leads to an overall anticorrelation between crustal thickness and bulk lower crustal velocity. Since this anticorrelation is contrary to crustal thickening resulting from passive upwelling of abnormally hot mantle, it is necessary to consider active upwelling components and/or some compositional heterogeneities in the mantle source. The NW limit of the Malpelo Ridge shows a dramatic crustal thinning and displays high lower crustal velocities and a poorly defined crust‐mantle boundary, suggesting that differential motion along the Coiba transform fault probably separated Regina and Malpelo Ridges

Three-dimensional P-wave velocity structure on the shallow part of the Central Costa Rican Pacific margin from Local Earthquake Tomography using off- and onshore networks

The Central Costa Rican Pacific margin is characterized by a high-seismicity rate, coincident with the subduction of rough-relief ocean floor and has generated earthquakes with magnitude up to seven in the past. We inverted selected P-wave traveltimes from earthquakes recorded by a combined on- and offshore seismological array deployed during 6 months in the area, simultaneously determining hypocentres and the 3-D tomographic velocity structure on the shallow part of the subduction zone (<70 km). The results reflect the complexity associated to subduction of ocean-floor morphology and the transition from normal to thickened subducting oceanic crust. The subducting slab is imaged as a high-velocity perturbation with a band of low velocities (LVB) on top encompassing the intraslab seismicity deeper than ∼30 km. The LVB is locally thickened by the presence of at least two subducted seamounts beneath the margin wedge. There is a general eastward widening of the LVB over a relatively short distance, closely coinciding with the onset of an inverted forearc basin onshore and the appearance of an aseismic low-velocity anomaly beneath the inner forearc. The latter coincides spatially with an area of the subaerial forearc where differential uplift of blocks has been described, suggesting tectonic underplating of eroded material against the base of the upper plate crust. Alternatively, the low velocities could be induced by an accumulation of upward migrating fluids. Other observed velocity perturbations are attributed to several processes taking place at different depths, such as slab hydration through outer rise faulting, tectonic erosion and slab dehydration