Error propagation in equations for geochemical modeling of radiogenic isotopes in two-component mixing

This paper presents error propagation equations for modeling of radiogenic isotopes during mixing of two components or end-members. These equations can be used to estimate errors on an isotopic ratio in the mixture of two components, as a function of the analytical errors or the total errors of geological field sampling and analytical errors. Two typical cases (“Small errors” and “Large errors”) are illustrated for mixing of Sr isotopes. Similar examples can be formulated for the other radiogenic isotopic ratios. Actual isotopic data for sediment and basalt samples from the Cocos plate are also included to further illustrate the use of these equations. The isotopic compositions of the predicted mixtures can be used to constrain the origin of magmas in the central part of the Mexican Volcanic Belt. These examples show the need of high quality experimental data for them to be useful in geochemical modeling of magmatic processes.     

Geodynamic evolution of southern Costa Rica related to low-angle subduction of the Cocos Ridge: constraints from thermochronology

The Late Tertiary shallow subduction of the Cocos ridge under the Caribbean plate controlled the evolution of the Cordillera de Talamanca in southeast Costa Rica, which is a mountain range that consists mainly of granitoids formed in a volcanic arc setting. Fission track thermochronology using zircon and apatite, as well as ⁴⁰Ar–³⁹Ar and Rb–Sr age data of amphibole and biotite in granitoid rocks constrain the thermal history of the Cordillera de Talamanca and the age of onset of subduction of the Cocos ridge. Shallow intrusion of granitoid melts resulted in fast and isobaric cooling. A weighted mean zircon fission track age (13 analyses) and Rb–Sr biotite ages of about 10 Ma suggest rapid cooling and give minimum ages for granitoid emplacement. In some cases ⁴⁰Ar–³⁹Ar and Rb–Sr apparent ages of amphibole and biotite are younger than the zircon fission track ages, which can be attributed to partial resetting by hydrothermal alteration. Apatite fission track ages range from 4.8 to 1.7 Ma but show no correlation with the 3090-m elevation span over which they were sampled. The apatite ages seem to indicate rapid exhumation caused by tectonic and isostatic processes. The combination of the apatite fission track ages with subduction parameters of the Cocos plate such as subduction angle, plate convergence rate and distance of the Cordillera de Talamanca to the trench implies that the Cocos ridge entered the Middle America Trench between 5.5 and 3.5 Ma.     

科科斯脊玄武岩斜长石矿物化学及地质意义

综合大洋钻探计划(IODP) 334和344航次在U1381站位处的两个钻孔(A孔和C孔)获得了中美洲西海岸外科科斯脊基底拉斑玄武岩,对其岩浆过程开展研究可为理解其岩石成因提供重要依据。本文对科科斯脊玄武岩中斜长石斑晶和微晶进行了详细的原位主微量元素分析,结果表明,斜长石种属为培长石、拉长石及少量中长石。部分斜长石斑晶具有正环带结构;但多数斜长石斑晶不具有明显环带,仅从核部到边部存在微弱的成分变化。斜长石斑晶与微晶的微量元素差别较大:斜长石斑晶富集轻稀土和大离子亲石元素、亏损高场强元素,且具有明显的Eu正异常;斜长石微晶不相容元素含量通常高于斜长石斑晶。根据斜长石温度计计算获得斜长石斑晶结晶温度为1 050～1 253℃,斜长石微晶结晶温度为866～1 033℃。基于以上特征,推测斜长石斑晶核部是相对原始岩浆的产物,而斑晶边部以及微晶是演化岩浆的结晶产物。斜长石斑晶的成分变化及熔蚀麻点结构是由于岩浆补给及岩浆减压上升造成的。最后,本研究推测科科斯脊基底玄武岩来自于开放的岩浆房,且岩浆房内可能存在原始岩浆的不断注入及岩浆对流。     

Compilation of radiogenic isotope data in Mexico and their petrogenetic implications

Seven hundred and twenty-five Sr, two hundred and forty-three Nd and one hundred and fifty-one Pb isotopic ratios from seven different Mexican magmatic provinces were compiled in an extensive geochemical database. Data were arranged according to the Mexican geological provinces, indicating for each province total number of analyses, range and mean of values and two times standard deviation (2σ). Data from seven provinces were included in the database: Mexican Volcanic Belt (MVB), Sierra Madre Occidental (SMO), Baja California (BC), Pacific Ocean (PacOc), Altiplano (AP), Sierra Madre del Sur (SMS), and Sierra Madre Oriental (SMOr). Isotopic values from upper mantle and lower crustal xenoliths, basement outcrops and sediments from the Cocos Plate were also compiled. In the MVB the isotopic ratios range as follows:⁸⁷Sr/⁸⁶Sr 0.703003-0.70841;¹⁴³Nd/¹⁴⁴Nd 0.512496-0.513098;²⁰⁶Pb/²⁰⁴Pb 18.567-19.580;²⁰⁷Pb/²⁰⁴Pb 15.466-15.647;²⁰⁸Pb/²⁰⁴Pb 38.065-38.632. The SMO shows a large variation in⁸⁷Sr/⁸⁶Sr ranging from ∼0.7033 to 0.71387.¹⁴³Nd/¹⁴⁴Nd ratios are relatively less variable with values from 0.51191 to 0.51286. Pb isotope ratios in the SMO are as follows:²⁰⁶Pb/²⁰⁴Pb 18.060-18.860;²⁰⁷Pb/²⁰⁴Pb 15.558-15.636;²⁰⁸Pb/²⁰⁴Pb 37.945-38.625. PacOc rocks show the most depleted Sr and Nd isotopic ratios (0.70232-0.70567 for Sr and 0.512631-0.513261 for Nd). Pb isotopes for PacOc show the following range:²⁰⁶Pb/²⁰⁴Pb 18.049-19.910;²⁰⁷Pb/²⁰⁴⁷Pb 15.425-15.734;²⁰⁸Pb/²⁰⁴Pb 37.449-39.404. The isotopic ratios of the AP rocks seem to be within the range of those from the PacOc. Most samples with reported Sr and Nd isotopic data are spread within and around the “mantle array”. The SMO seems to have been formed by a mixing process between mantle derived magmas and continental crust. The MVB appears to have a larger mantle component, with AFC as the dominant petrogenetic process for the evolved rocks. There is still a need for Pb isotopic data in all Mexican magmatic provinces and of Nd isotopes in BC, AP, SMS, and SMOr.     

Measurements of upper mantle shear wave anisotropy from a permanent network in southern Mexico

Upper mantle shear wave anisotropy under stations in southern Mexico was measured using records of SKS phases. Fast polarization directions where the Cocos plate subducts subhorizontally are oriented in the direction of the relative motion between the Cocos and North American plates, and are trench-perpendicular. This pattern is interpreted as subslab entrained flow, and is similar to that observed at the Cascadia subduction zone. Earlier studies have pointed out that both regions have in common the young age of the subducting lithosphere. Changes in the orientation of the fast axes are observed where the subducting plates change dip and/or are torn, and are thus indicative of 3-D flow around the slab edges. They are consistent with slab rollback, as previously shown by other authors. Some stations located away from the plate boundaries have their fast directions controlled by the absolute motion of the North American plate. The fast axis for station ZAIG, located in the Mesa Central, is oriented WNW-ESE and is different from all the other measurements in this study.     

Middle Miocene to present plate tectonic history of the southern Central American Volcanic Arc

New mid Miocene to present plate tectonic reconstructions of the southern Central American Volcanic Arc (CAVA) reveal that the inception of Cocos Ridge subduction began no earlier than 3 Ma, and possibly as late as 2 Ma. The Cocos Ridge has been displaced from the Malpelo Ridge to the southeast since 9 Ma along the Panama Fracture Zone (PFZ) system. Ambiguous PFZ and Coiba Fracture Zone (CFZ) interaction since 9 Ma precludes conclusively establishing the age of initial Cocos Ridge subduction. Detailed reconstructions based on magnetic anomalies offshore reveal several other variations in subduction parameters beneath southern Central America that preceded subduction of the Cocos Ridge, including southeastward migration of the Nazca–Cocos–Caribbean triple junction along the Middle America Trench (MAT) from 12 Ma to present, and subduction of ≤2 km high scarps both parallel and perpendicular to the trench from 6 to 1 Ma.The timing of changes in subduction processes has commonly been determined by (and correlated with) geologic changes in the upper plate. However, reliable ⁴⁰Ar/³⁹Ar dating of these events has become available only recently [Abstr. Programs-Geol. Soc. Am. (2002)]. These new dates better constrain the magmatic and structural history of southern Costa Rica. Observations from this data set include: a gap in the volcanic record from 11 to 6 Ma, which coincides temporally with emplacement of most plutons in southern Costa Rica, normal arc volcanism ceased after 3.5 Ma in southern Costa Rica, and Pliocene (mostly 1.5 Ma) adakite volcanism was widely distributed from central Panama to southern Costa Rica (though volumetrically insignificant).This new data reveals that many geologic phenomena, commonly attributed to subduction and underplating of the buoyant Cocos Ridge, in fact precede inception of Cocos Ridge subduction and seem to correlate more favorably in time with earlier tectonic events. Adakite volcanic activity corresponds in space and time with the subduction of a large scarp associated with a tectonic boundary off southern Panama. Regional unconformities and an 11–6 Ma gap in arc volcanism match temporally with oblique subduction of the Nazca plate beneath central and southern Costa Rica. Cessation of volcanic activity, low-temperature cooling of plutons in the Cordillera de Talamanca (CT), and rapid increases in sedimentation in the fore-arc and back-arc basins coincide with passage of the Nazca–Cocos–Caribbean triple junction and initiation of subduction of “rough” crust associated with Cocos–Nazca rifting 3.5 Ma, closely followed by initial subduction of the Cocos Ridge 2–3 Ma. None of the aforementioned geologic events occurred at a time that would allow for underplating by the Cocos Ridge. Rather they are probably related to complex interactions with subduction of complicated plates offshore. All of the aforementioned events indicate that the southern Central American subduction system has been in flux since at least 12 Ma.     

A glaring omission in Australia’s marine conservation planning

The primary goal of Australia’s National Representative System of Marine Protected Areas (NRSMPA) is to establish a comprehensive, adequate and representative system of MPAs. This study identifies a glaring contradiction to this policy. The provinces of Christmas and Cocos Islands are among the most unique and threatened marine bioregions in Australia, yet receive no protection from the NRSMPA. The lack of protection appears to be due to difficulties with multiple governance arrangements and other political priorities. These issues have already caused biodiversity loss in the terrestrial environment of these bioregions. The Australian Government must include the Christmas and Cocos provinces in the NRSMPA otherwise it risks irreversible loss of marine biodiversity in these unique bioregions.     

Tectonic analysis and paleo-stress determination of the upper lava section at ODP/IODP site 1256 (East Pacific Ocean)

Research on the deep sea is of great importance for a better understanding of the mechanism of magma emplacement and the tectonic evolution of oceanic crust. However, details of the internal structure in the upper levels of the oceanic crust are much less complete than that of the more fully studied sub-aerial areas. For the first time, this study proposes a dynamic analysis using the inversion method on core data derived from the drilled basement of the present-day intact oceanic crust at ODP/IODP Site 1256 in the Cocos plate. The research is based on an innovative core reorientation process and combines different stress hypothesis approaches for the analysis of heterogeneous failure-slip data via exploitation of two distinct techniques. From the analysis of the failure-slip data, both techniques produce 5 distinct subsystem datasets. All calculated subsystems are mechanically and geometrically admissible. Interpretation of the results allows the researchers to note a complex local and regional tectonic evolution deriving from the interplay of (1) the ridge push and rotation of both the East Pacific Rise and the Cocos-Nazca Spreading Center, (2) the effect of the slab pull of the Middle America Trench, (3) the influence of lava emplacement mechanisms, and (4) intra-plate deformation.     

Seismic activity along the central america volcanic arc: is it related to subduction of the Cocos plate?

We determine seismic strain rate of tectonic earthquakes along the Central America Volcanic Arc. We then compare this result to those obtained from earthquakes related to the convergence of the Cocos and Caribbean plates and to earthquakes in the back-arc region of northern Central America.

The seismic strain-rate tensor for shallow-focus earthquakes along the Central America volcanic arc since 1700, has a compressive eigenvector with a magnitude of 0.7 × 10⁻⁸ year⁻¹, and oriented in a 357° azimuth. The extensive eigenvector is oriented in a 86° azimuth, with a magnitude of 0.82 × 10⁻⁸ year⁻¹. When only Centroid Moment-tensor solutions (CMT) are considered, the respective eigenvectors are 1.2 × 10⁻⁸ year⁻¹ and 1.0 × 10⁻⁸ year⁻¹.

The compressive eigenvector from the seismic strain-rate tensor for earthquakes along the Cocos-Caribbean convergent margin is 2.0 × 10⁻⁸ year⁻¹, plunging at 25°, and oriented in a 29° azimuth. Its magnitude and direction are similar to those of the compressive eigenvector for earthquakes along the volcanic arc. The extensive eigenvector along the convergent margin, on the other hand, has a large vertical component. The compressive and extensive eigevenvectors are 4.9 × 10⁻⁸ year⁻¹ and 4.6 × 10⁻⁸ year⁻¹, using only CMTs as the database.

Earthquakes along the grabens of northern Central America yield a seismic strain-rate tensor whose extensive eigenvector has a magnitude of 2.4 × 10⁻⁸ year⁻¹, oriented in a 109° azimuth. Magnitude and direction are similar to those of the extensive eigenvector for earthquakes along the volcanic arc. The compressive eigenvector along the grabens is practically vertical.

Similarities in magnitudes and directions for compressive and extensive eigenvectors suggest to us that the strain field along the Central America volcanic arc is the result of compression along the convergent Cocos-Caribbean margin, and extension in the back-arc region, along the grabens of northern Central America. This field is resolved as strike-slip faulting along the arc.