A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Geochemistry of oceanic flood basalts from the Caribbean plate

The thick oceanic crust of the Caribbean plate appears to be the tectonized remnant of an eastern Pacific oceanic plateau that has been inserted between North and South America. The emplacement of the plateau into its present position has resulted in the obduction and exposure of its margins, providing an opportunity to study the age relations, internal structure and compositional features of the plateau. We present the results of 40Ar-39Ar radiometric dating, major-, trace-element, and isotopic compositions of basalts from some of the exposed sections as well as drill core basalt samples from Leg 15 of the Deep Sea Drilling Project. Five widely spaced, margin sections yielded ages ranging from 91 to 88 Ma. Less well-constrained radiometric ages from the drill cores, combined with the biostratigraphic age of surrounding sediments indicate a minimum crystallization age of ~90 Ma in the Venezuelan Basin. The synchroneity of ages across the region is consistent with a flood basalt origin for the bulk of the Caribbean plateau i.e., large volume, rapidly erupted, regionally extensive volcanism.. The ages and compositions are also consistent with plate reconstructions that place the Caribbean plateau in the vicinity of the Galápagos hotspot at its inception. The trace-element and isotopic compositions of the ~90 Ma rocks indicate a depleted mantle and an enriched, plume-like mantle were involved in melting to varying degrees across the plateau. Within the same region, a volumetrically secondary, but widespread magmatic event occurred at 76 Ma, as is evident in Curacao, western Colombia, Haiti, and at DSDP Site 152/ODP Site 1001 near the Hess Escarpment. Limited trace-element data indicate that this phase of magmatism was generally more depleted than the first. We speculate that magmatism may have resulted from upwelling of mantle, still hot from the 90 Ma event, during lithospheric extension attending gravitational collapse of the plateau, andror tectonic emplacement of the plateau between North and South America. Still younger volcanics are found in the Dominican Republic (69 Ma) and the Quepos Peninsula of Costa Rica (63 Ma). The latter occurrence conceivably formed over the Galápagos hotspot and subsequently accreted to the western edge of the plateau during subduction of the Farallon plate.

Geochemistry of oceanic crust beneath Costa Rica

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge (Feigenson et al., 2004, doi:10.1029/2003GC000621; Herrstom et al., 1995, doi:10.1130/0091-7613(1995)023<0617:VILCAW>2.3.CO;2; Abratis and Woerner, 2001) or eroded fore-arc complexes (Goss and Kay, 2006, doi:10.1029/2005GC001163) but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (~63-190 mm/yr) and are comparable to the magnitude of subducting Cocos plate motion (approx85 mm/yr). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.

Pb concentrations, Pb isotopic composition, and amount of metalliferous components in pelagic sediments from the Pacific Ocean

The amount of lead annually transferred from oceanic crust to metalliferous sediments was estimated in order to test the hypothesis that a non-magmatic flux of lead causes the Pb surplus in the continental crust. A Pb surplus has been inferred from global crust-mantle lead mass balances derived from lead concentration correlations with other trace elements and from lead isotope systematics in oceanic basalts. DSDP/ODP data on the amount of metalliferous sediments in the Pacific Ocean and along a South Atlantic traverse are used to calculate the mean worldwide thickness of 3 (+/-1) m for purely metalliferous sediment componens. Lead isotope ratios of 39 metalliferous sediments from the Pacific define mixing lines between continent-derived (seawater) and mantle-derived (basaltic) lead, with the most metal-rich sediments usually having the most mantle-like Pb isotope composition. We used this isotope correlation and the Pb content of the 39 metalliferous sediments to derive an estimate of 130 (+/-70) µg/g for the concentration of mantle-derived lead in the purely metalliferous end-member. Mass balance calculations show that at least 12 (+/-8)% of the lead, annually transferred from upper mantle to oceanic crust at the ocean ridges, is leached out by hydrothermal processes and re-deposited in marine sediments. If all of the metalliferous lead is ultimately transferred to the continental crust during subduction, the annual flux of this lead from mantle to continental crust is 2.6 (+/-2.0) * 10**6 kg. Assuming this transfer rate to be proportional to the rate of oceanic plate production, one can fit the lead transfer to models of plate production rate variations through time. Integrating over 4 Ga, hydrothermal lead transfer to the continental crust accounts for a significant portion of the Pb surplus in the continental crust. It therefore appears to be one of the main reasons for the anomalous behavior of lead in the global crust-mantle system.

Tephrostratigraphy and radiometric dating of authigenic carbonates

At convergent margins, fluids rise through the forearc in response to consolidation of the upper plate and dewatering of the subducting plate, and produce various cold-seep-related features on the seafloor (mud diapirs, mud mounds). At the Central American forearc, authigenic carbonates precipitated from rising fluids within such structures during active venting while typical mixed-mud sediments were ejected onto the surrounding seafloor where they became intercalated with normal pelagic background sediments, indicating that mud mounds evolved unsteadily through alternating active and inactive phases. Intercalated regional ash layers from Plinian eruptions at the Central American volcanic arc provide time marks that constrain the ages of mud ejection activity. U/Th dating of drill core samples of authigenic carbonate caps of mud mounds yields ages agreeing well with those constrained by ash layers and showing that carbonate caps grow inward rather than outward during active venting. Both dating approaches show that offshore Nicaragua and Costa Rica (1) active and inactive phases can occur simultaneously at neighboring mounds, (2) mounds along the forearc have individual histories of activity, but there are distinct time intervals when nearly all mounds have been active or inactive, (3) lifetimes of mounds reach several hundred thousand years, and (4) highly active periods last 10-50 k.y. with intervening periods of >10 k.y. of relative quiescence.