Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

Isotopic geochemistry of granitoids in the Jinshajiang suture zone, SW China

The Jinshajiang suture zone, located in the eastern part of the Tethyan tectonic domain, is noticeable for a large-scale distribution of Late Jurassic to Triassic granitoids. These granitoids were genetically related to the evolution of the Paleo-Tethys Ocean. The Beiwu, Linong and Lunong granitoids occur in the middle zone of the Jinshajiang Suture Zone, and possess similar geochemical features, indicating they share a common magma source. SIMS zircon U-Pb dating reveals the Beiwu, Linong and Lunong granitic intrusions were emplaced at 233.9±1.4 Ma (2 sigma), 233.1 ±1.4 Ma (2 sigma) and 231.0±1.6 Ma (2 sigma), respectively. All of these granitoids are enriched in abundances of Si (SiO2 =65.2-73.5 wt.%), and large-ion-lithophile-elements (LILEs), but depleted in high-field-strength-elements contents (HFSEs, e.g., Nb, Ta, Ti). In addition, they have low P2O5 contents (0.06-0.11 wt.%), A/CNK values ([molecular Al2O3/(CaO+Na2O+K2O)], mostly<1.1) and 10000Ga/Al ratios (1.7-2.2), consistent with the characteristics of I-type granites. In terms of isotopic compositions, these granitoids have high initial 87Sr/86Sr ratios (0.7078-0.7148), Pb isotopic compositions [(206Pb/204Pb)t=18.213-18.598, (207Pb/204Pb)t=15.637-15.730 and (208Pb/204Pb)t=38.323-38.791], zircon d18O values (7. per mil-9.3 per mil) and negative eNd(t) values (-5.1 to -6.7), suggesting they were predominantly derived from the continental crust. Their Nb/Ta ratios (average value=8.6) are consistent with those of the lower continental crust (LCC). However, variable ?Hf(t) values (-8.6 to +2.8) and the occurrences of mafic microgranular enclaves (MMEs) suggest that mantle-derived melts and lower crustal magmas were involved in the generation of these granitoids. Moreover, the high Pb isotopic ratios and elevated zircon d18O values of these rocks indicate a significant contribution of the upper crustal composition. We propose a model in which the Beiwu, Linong and Lunong granitoids were generated under a late collisional or post-collisional setting. It is possible that this collision was completed before Late Triassic. Decompression induced mantle-derived magmas underplated and provided the heat for the anatexis of the crust. Hybrid melts including mantle-derived and the lower crustal magmas were then generated. The hybrid melts thereafter ascended to a shallow depth and resulted in some degree of sedimentary rocks assimilation. Such three-component mixing magmas source and subsequent fractional crystallization could be responsible for the formation of the Beiwu, Linong and Lunong granitoids.

Geochemistry of subducted volcaniclastic sediments from Mariana arc magmas

The phase relations of natural volcaniclastic sediments from the west Pacific Ocean were investigated experimentally at conditions of 3-6 GPa and 800-900 °C with 10 wt.% added H2O (in addition to ~ 10 wt.% structurally-bound H2O) to induce hydrous melting. Volcaniclastic sediments are shown to produce a sub-solidus assemblage of garnet, clinopyroxene, biotite, quartz/coesite and the accessory phases rutile ± Fe-Ti oxide ± apatite ± monazite ± zircon. Hydrous melt appears at temperatures exceeding 800-850 °C, irrespective of pressure. The melt-producing reaction consumes clinopyroxene, biotite and quartz/coesite and produces orthopyroxene. These phase relations differ from those of pelagic clays and K-bearing mid ocean ridge basalts (e.g. altered oceanic crust) that contain phengite, rather than biotite, as a sub-solidus phase. Despite their relatively high melt productivity, the wet solidus for volcaniclastic sediments is found to be higher (825-850 °C) than other marine sediments (700-750 °C) at 3 GPa. This trend is reversed at high-pressure conditions (6 GPa) where the biotite melting reaction occurs at lower temperatures (800-850 °C) than the phengite melting reaction (900-1000 °C). Trace element data was obtained from the 3 GPa run products, showing that partial melts are depleted in heavy rare earth elements (REE) and high field strength elements (HFSE), due to the presence of residual garnet and rutile, and are enriched in large ion lithophile elements (LILE), except for Sr and Ba. This is in contrast to previous experimental studies on pelagic sediments at sub-arc depths, where Sr and Ba are among the most enriched trace elements in glasses. This behavior can be partly attributed to the presence of residual apatite, which also host some light REE in our supra-solidus residues. Our new experimental results account for a wide range of trace element and U-series geochemical features of the sedimentary component of the Mariana arc magmas, including imparting a substantial Nb anomaly to melts from an anomaly-free protolith.

Sediment composition and sedimentation rates of four Holes in ODP Leg 166

To date, work on the Great Bahama Bank's western, leeward margin has centred chiefly on seismic-scale expressions of carbonate sequences and systems tracts. However, periplatform, slope sediments also exhibit very well developed cyclicity on scales of decimetres to several metres. It is these small-scale, high-frequency cycles within the larger-scale facies successions of the Quaternary which form the main topic of this paper. Previous studies have shown that the small-scale cycles correlate to the orbitally forced, high-frequency sea-level changes. Therefore these cycles should indicate how sea level has affected the slope development and thus platform-margin evolution during this period. Through detailed, high-resolution sequence stratigraphy of the Great Bahama Bank's leeward margin, obtained via delta18O isotope and mineralogical (XRD) analyses, confined by U/Th dating and nannofossil bioevents, a greater understanding of the bedding geometries within the Pleistocene-Holocene seismic sequences and clues as to the nature of the slope development has been achieved. The high-resolution seismic profiles indicate that since the Plio-Pleistocene change in geometry, in which the Great Bahama Bank developed into a rimmed platform, continued steepening and subsequent progradation of the leeward margin has typified slope development during the Quaternary, which is described as an accretionary slope. However, on the basis of our observations we conclude that only the early to lower middle Pleistocene section (isotope stages 45-20) and the Holocene (isotope stage 1) of the leeward margin is accretionary. This indicates that a degree of erosion and/or by-passing has occurred on the leeward margin since the lower middle Pleistocene (isotope stage 19). During the first part of this period (isotope stages 19-12) erosion and/or by-passing occurred in the middle to lower slope regions and toe-of-slope. By the end of the upper middle to late Pleistocene phase (isotope stages 11-2) erosion also occurred on the upper slope. This erosion by currents at the toe-of-slope and oversteepening of the upper and middle slopes have led to back-cutting upslope and resulted in the progressive retreat of the toe-of-slope towards the platform to the east. However, the rise in sea level since the Last Glacial Maximum to its present-day level has allowed high productivity on the platform top during the Holocene and the deposition of a thick sediment wedge on the slope and sedimentation across the entire leeward flanks. This has led to the redevelopment of an accretionary slope and continued westward progradation of the Great Bahama Bank's western, leeward margin.

Dead Sea lake level for the last 140 ka

In this paper we describe the stratigraphy and sediments deposited in Lake Samra that occupied the Dead Sea basin between ~135 and 75 ka. This information is combined with U/Th dating of primary aragonites in order to estimate a relative lake-level curve that serves as a regional paleohydrological monitor. The lake stood at an elevation of ~340 m below mean sea level (MSL) during most of the last interglacial. This level is relatively higher than the average Holocene Dead Sea (~400±30 m below MSL). At ~120 and ~85 ka, Lake Samra rose to ~320 m below MSL while it dropped to levels lower than ?380 m below MSL at ~135 and ~75 ka, reflecting arid conditions in the drainage area. Lowstands are correlated with warm intervals in the Northern Hemisphere, while minor lake rises are probably related to cold episodes during MIS 5b and MIS 5d. Similar climate relationships are documented for the last glacial highstand Lake Lisan and the lowstand Holocene Dead Sea. Yet, the dominance of detrital calcites and precipitation of travertines in the Dead Sea basin during the last interglacial interval suggest intense pluvial conditions and possible contribution of southern sources of wetness to the region.

Geochemistry of Lesser Antilles lavas

We present new U-series disequilibrium and radiogenic isotope data for 7 mafic lavas from the Lesser Antilles arc. These are combined with published data in an internally consistent model that quantitatively estimates the amount of sediment and fluid added to the source of the Lesser Antilles arc system. Some lavas form an array consistent with bulk sediment addition (0.2-2%) whereas others appear to require addition of 0.4-2% sediment melt, particularly in the south of the arc. Evidence for both bulk sediment and sediment melt addition can be found within both the northern and central sections of the arc suggesting a thermal structure whereby the upper portions of the subducted sediment pile lie close to their solidus beneath much of the arc. Addition of up to 5% fluid derived from altered oceanic crust to these sediment enriched mantle wedge source regions can simulate the majority of the lavas on a plot of 207Pb/204Pb versus Ce/Pb. By taking into account the range in calculated wedge compositions and allowing for some mobility of Th in the fluid, the same model can also account for much of the observed range of U-Th-Ra disequilibria, especially if the eclogitic residue contains trace amounts of rutile. The implication of this more complex model is that the time scales for fluid addition and differentiation could be significantly shorter than those estimated in some previous studies.

A deep-sea coral record of the North Atlantic

Our record of Younger Dryas intermediate-depth seawater D14C from North Atlantic deep-sea corals supports a link between abrupt climate change and intermediate ocean variability. Our data show that northern source intermediate water (~1700 m) was partially replaced by 14C-depleted southern source water at the onset of the event, consistent with a reduction in the rate of North Atlantic Deep Water formation. This transition requires the existence of large, mobile gradients of D14C in the ocean during the Younger Dryas. The D14C water column profile from Keigwin (2004) provides direct evidence for the presence of one such gradient at the beginning of the Younger Dryas (~12.9 ka), with a 100 per mil offset between shallow (<~2400 m) and deep water. Our early Younger Dryas data are consistent with this profile and also show a D14C inversion, with 35 per mil more enriched water at ~2400 m than at ~1700 m. This feature is probably the result of mixing between relatively well 14C ventilated northern source water and more poorly 14C ventilated southern source intermediate water, which is slightly shallower. Over the rest of the Younger Dryas our intermediate water/deepwater coral D14C data gradually increase, while the atmosphere D14C drops. For a very brief interval at ~12.0 ka and at the end of the Younger Dryas (11.5 ka), intermediate water D14C (~1200 m) approached atmospheric D14C. These enriched D14C results suggest an enhanced initial D14C content of the water and demonstrate the presence of large lateral D14C gradients in the intermediate/deep ocean in addition to the sharp vertical shift at ~2500 m. The transient D14C enrichment at ~12.0 ka occurred in the middle of the Younger Dryas and demonstrates that there is at least one time when the intermediate/deep ocean underwent dramatic change but with much smaller effects in other paleoclimatic records.

Coral Sr/Ca records from Bonaire Holocene corals

Proxy reconstructions of tropical Atlantic sea surface temperature (SST) that extend beyond the period of instrumental observations have primarily focused on centennial to millennial variability rather than on seasonal to multidecadal variability. Here we present monthly-resolved records of Sr/Ca (a proxy of SST) from fossil annually-banded Diploria strigosa corals from Bonaire (southern Caribbean Sea). The individual corals provide time-windows of up to 68 years length, and the total number of 295 years of record allows for assessing the natural range of seasonal to multidecadal SST variability in the western tropical Atlantic during snapshots of the mid- to late Holocene. Comparable to modern climate, the coral Sr/Ca records reveal that mid- to late Holocene SST was characterised by clear seasonal cycles, persistent quasi-biennial and prominent interannual as well as inter- to multidecadal-scale variability. However, the magnitude of SST variations on these timescales has varied over the last 6.2 ka. The coral records show increased seasonality during the mid-Holocene consistent with climate model simulations indicating that southern Caribbean SST seasonality is induced by insolation changes on orbital timescales, whereas internal dynamics of the climate system play an important role on shorter timescales. Interannual SST variability is linked to ocean-atmosphere interactions of Atlantic and Pacific origin. Pronounced interannual variability in the western tropical Atlantic is indicated by a 2.35 ka coral, possibly related to a strengthening of the variability of the El Niño/Southern Oscillation throughout the Holocene. Prominent inter- to multidecadal SST variability is evident in the coral records and slightly more pronounced in the mid-Holocene. We finally argue that our coral data provide a target for studying Holocene climate variability on seasonal and interannual to multidecadal timescales, when using further numerical models and high-resolution proxy data.

Sr/Ca ratios and oxygen isotopes from sclerosponges

We investigate aragonitic skeletons of the Caribbean sclerosponge Ceratoporella nicholsoni from Jamaica, 20 m below sea level (mbsl), and Pedro Bank, 125 mbsl. We use d18O and Sr/Ca ratios as temperature proxies to reconstruct the Caribbean mixed layer and thermocline temperature history since 1400 A.D. with a decadal time resolution. Our age models are based on U/Th dating and locating of the radiocarbon bomb spike. The modern temperature difference between the two sites is used to tentatively calibrate the C. nicholsoni Sr/Ca thermometer. The resulting calibration points to a temperature sensitivity of Sr/Ca in C. nicholsoni aragonite of about -0.1 mmol/mol/K. Our Sr/Ca records reveal a pronounced warming from the early 19th to the late 20th century, both at 20 and 125 mbsl. Two temperature minima in the shallow water record during the late 17th and early 19th century correspond to the Maunder and Dalton sunspot minima, respectively. Another major cooling occurred in the late 16th century and is not correlatable with a sunspot minimum. The temperature contrast between the two sites decreased from the 14th century to a minimum in the late 17th century and subsequently increased to modern values in the early 19th century. This is interpreted as a long-term deepening and subsequent shoaling of the Caribbean thermocline. The major trends of the Sr/Ca records are reproduced in both specimens but hardly reflected in the d18O records.