Geochemistry of loess,continental crustal composition and crustal model ages

Isotopic, major and trace element studies of loess deposits from America, China, Europe and New Zealand show general uniformity of composition. Silica, Zr and Hf are enriched relative to estimates of bulk composition of the upper continental crust. The REE data are indistinguishable from those of average shales, confirming the concept that these REE patterns ($\text{La}{}_{\mathrm{N}}\text{/}\text{Yb}{}_{\mathrm{N}}\text{= 9.5}\text{Eu}\text{/}\text{Eu}{}^1\text{= 0.66}$) represent the upper crustal average. Sm-Nd model ages are variable but <1700 m.y. They reflect derivation from younger elevated erogenic areas subject to Pleistocene glaciation. Although Sm-Nd model ages vary by a factor of two, the REE patterns remain constant. This indicates that processes responsible for formation of the upper crust have produced no secular change in composition since the mid-Proterozoic.     

Earth''s earliest continental lithosphere, hydrothermal flux and crustal recycling

The Kaapvaal craton in southern Africa and the Pilbara craton of northwestern Australia are the largest regions on Earth to have retained relatively pristine mid-Archaean rocks (3.0–4.0 Ga).The Kaapvaal craton covers about 1.2×10⁶ km², and varies in lithospheric thickness between 170 and 350 km. At surface, the craton can be subdivided into a number of Archaean sub-domains; some of the subdomains are also well defined at depth, and local variations in tomography of the lithosphere correspond closely with subdomain boundaries at surface.The Archaean history of the Kaapvaal craton spans about 1 Gyr and can be conveniently subdivided into two periods, each of about the same length as the Phanerozoic. The first period, from circa 3.7-3.1 Ga, records the initial separation of the cratonic lithosphere from the asthenosphere, terminating with a major pulse of accretion tectonics between 3.2 and 3.1 Ga, which includes the formation of “paired metamorphic belts”. This period of continental growth can be compared to plate tectonic processes occurring in modern-day oceanic basins. However, the difference is that in the mid-Archaean, these oceanic processes appear to have occurred in shallower water depths than the modern ocean basins. The second period, from circa 3.1-2.6 Ga, records intra-continental and continental-edge processes: continental growth during this period occurred predominantly through a combination of tectonic accretion of crustal fragments and subduction-related igneous processes, in much the same way as has been documented along the margins of the Pacific and Tethys oceans since the Mesozoic.The intra-oceanic processes resulted in small, but deep-rooted continental nucleii; the first separation of this early continental lithosphere could only have occurred when the mean elevation of mid-oceanicridges sank below sea-level. Substantial recycling of continental lithosphere into the mantle must have occurred during this period of Earth history. During the second period, at least two large continental nucleii amalgamated during collisional processes which, together with internal chemical differentiation processes, created the first stable continental landmass. This landmass, which is known to have been substantially bigger than its present outline, may have been part of the Earth's first supercontinent.The oldest known subdomains of the craton include the oceanic-like rocks of the Barberton greenstone belt. The comagmatic mafic-ultramafic rocks (3.48–3.49 Ga) of this belt represent a remnant of very early oceanic-like lithosphere (known as the Jamestown Ophiolite Complex), which was obducted, approximately 45 Ma after its formation, onto a volcanic arc-like terrain by processes similar to those which have emplaced modern ophiolites at convergent margins of Phanerozoic continents. The early metamorphic history, metamorphic mineralogy, oxygen isotope profiles and degree of hydration of the 3.49 Ga Jamestown Ophiolite Complex are similar to present day subseafloor hydrothermal systems. The ratio of ΔMg to ΔSi for hydrothermally altered igneous rocks, both present day and Archaean, are remarkably uniform at −5(±0.9) and the same as that of hydrothermal fluids venting on the present-day East Pacific Rise. This observation suggests that the process of Mg exchange for Si in hydrothermal systems was commonplace throughout Earth's history.The chemistry of vent fluids and hydrothermally altered igneous rocks was combined with an inventory of ³He in the mantle to model Earth's total hydrothermal flux. An Archaean flux (at 3.5 Ga) of about 10 times present day was accompanied by a correspondingly greater abundance of Mg(OH), SiO₂, carbonate and Fe---Mn metasomatic rock types as well as massive sulphides. Assuming a constant column of seawater since the Archaean, the average residence time of seawater in the oceanic crust was 1.65−8.90×10⁵ years in the Archaean. Assuming that ³He and heat are transported from the mantle in silicate melts in uniform proportions, the model stipulates that accretion of oceanic crust decreased from about 3.43−6.5×10¹⁷ g/yr to a present-day rate of 0.52−0.8×10¹⁷ g/yr, with a drop in heat flow from 1.4−2.6×10²⁰ cal/yr to 2.1−3.2×10¹⁹ cal/year.The total amounts of SiO₂ and Fe mobilised in marine hydrothermal systems since 3.5 Ga is less than their masses in the present exosphere reservoirs (crust, hydrosphere, atmosphere). The total amounts of Mg, K, CO₂, Ca and Mn are greater than their respective masses in exosphere reservoirs; therefore, they must have been recycled into mantle. The total mass of recycled hydrothermal components is small compared to the mass of the mantle. The flux of volatiles in hydrothermal systems is large compared to their volume in the atmosphere suggesting that the CO₂ and O₂ budgets of the atmosphere have been influenced by hydrothermal processes, especially in the Archaean.     

Newania carbonatites, Western India: example of mantle derived magnesium carbonatites

The key mineralogical features of the Newania carbonatites, that illustrate their derivation from primary mantle melts (Gruau et al. Terra Nova, Abstract Suppl 1:336, 1995; Viladkar Petrology 6(3):272–283, 1998; Basu and Murty Abstracts of Goldschmidt Conference A40, 2006), are the presence of magnesite, graphite and Cr-rich magnetite. Magnesite is an early crystallizing phase. Cr-rich magnetite and graphite coexist with carbonatite minerals and precipitated from carbonate magma. Graphite, as well as gaseous CO₂ and carbonate minerals such as dolomite and magnesite, can be stable in peridotite mantle. Coexistence of these minerals is controlled by fO ₂ and PT-conditions. Mineral geothermometers for the Newania carbonatite give temperatures from 463 to 950°C. The parental source for Newania carbonatites was characterized by a relatively high log (fHF/fH₂O) level which increased during the crystallization history of Newania. The estimated oxygen fugacity (for ilmenite–magnetite pairs) varies from ?1.5 to +3.5 (log-bar unit deviation from FMQ buffer), which is supported by the presence of Fe-columbite, and the composition of phlogopite, amphibole and pyroxene that have an elevated concentration of Fe³⁺. However, the oxygen fugacity range represented by co-existing early-crystallized graphite and magnesite is below that of the FMQ buffer and lies on the CCO buffer.     

Geochemistry of the Rare Earth elements in spilites from the oceanic and continental crust

Analytical data including major elements, the Rare Earths, Ni, Cu, Zn, Rb, and Sr are presented for twenty-three spilites from the Mid-Atlantic Ridge, the Hercynian part of the Variscan geosyncline in Germany, and several localities in Switzerland. Low grade metamorphism (up to approximately 400° C) and spilitization of basaltic rocks apparently do not alter the original Rare Earth element (REE) distributions. This fact permits comparison of the spilites and unaltered tholeiitic basalts from the Mid-Atlantic Ridge. The relative REE distributions thus appear suitable for delineating the original basalt types of spilites formed by metamorphism. The spilites from the Hercynian part of the Variscan geosyncline have a REE distribution pattern which is characteristic of continental tholeiites. It is thus probable that during the formation of this geosyncline the principal magma extruded was of tholeiitic composition and that these rocks were later converted to spilites by metamorphism.     

Evidence for crustal components in the mantle and constraints on crustal recycling mechanisms: pyroxenite xenoliths from Hannuoba,North China

Spinel and garnet pyroxenite xenoliths in Cenozoic basalts from Hannuoba, North China show extremely heterogeneous chemical and isotopic compositions (ε_Nd=−27 to +34). Most of these pyroxenites are relatively young, probably late Mesozoic in age, although a few Al-pyroxenites could be very old (∼2 Ga). While their texture and major element compositions suggest an origin of high pressure cumulates, the trace element and isotopic compositions of the Hannuoba pyroxenites require multiple segregation processes from different parental magmas. Strong LREE enrichment, ubiquitous HFSE depletion and some Eu anomalies of the Al- and Cr-pyroxenites indicate the involvement of crust components in their source. Their Sr–Nd isotopic ratios are negatively correlated and plot below the MORB–OIB–IAB–sediment trend, suggesting that the parental melts of the Cr- and Al-pyroxenites may have been derived from a mixture of asthenospheric melts and a long-term evolved continental crust. The garnet pyroxenites significantly deviate from the isotopic array defined by the Al-pyroxenites, due to their relatively high ⁸⁷Sr/⁸⁶Sr at given ε_Nd. They thus more likely represent segregates from melts derived from partial melting of hydrothermally altered oceanic crust (basalts+marine sediment). If the crustal component involved in the Al-pyroxenites is subducted terrigenous sediments or other continental materials from the Archean Sino-Korean Craton, the Al-pyroxenites and garnet pyroxenites may have formed contemporaneously at a palaeo-convergent plate margin. This may be related to the subduction of the Mongol–Okhotsk plate beneath North China during the late Jurassic. Alternatively, if the delaminated lower crust was involved, it implies that most of the Al-pyroxenites are younger than the garnet pyroxenites, and their formation may be temporally correlated with lithospheric thinning during the Cretaceous. This model is attractive because the inferred tectonic evolution from a convergent setting to an extensional environment is consistent with the geologic record in the area.     

Dynamic Earth: crustal and mantle heterogeneity

The dynamic processes within the Earth leave their record in geophysical and geochemical variation about the general stratification with depth. A snapshot of current structure is provided by geophysical evidence, whereas geochemical information provides a perspective over the age of the Earth. The combination of information on the distribution of heterogeneity from geophysical and geochemical sources provides enhanced insight into likely geodynamic processes. A variety of techniques can be used to examine crustal structure, but the major source of information on seismic heterogeneity within the Earth comes from tomographic studies, exploiting surface waves for the lithosphere and body waves for the bulk Earth. A powerful tool for examining the character of mantle heterogeneity is the comparison of images of bulk-sound and shear-wave speed extracted in a single inversion, since this isolates the dependencies on the elastic moduli. Such studies are particularly effective when a common path coverage is achieved for P and S as, for example, when common source and receiver pairs are extracted for arrival times of the phases. The relative behaviour of the bulk-sound and shear-wave speeds allows the definition of heterogeneity regimes. For subduction zones, a large part of the imaged structure comes from S-wave speed variations. The narrow segments of fast wave speeds in the lower mantle, in the depth range 900 – 1500 km, are dominated by S variations, with very little bulk-sound contribution, so images of P-wave speed are controlled by shear-wave variability. Deep in the mantle, there are many features with high seismic-wave speed without an obvious association with subduction in the last 100 million years, which suggests long-lived preservation of components of the geodynamic cycle. The base of the Earth's mantle is a complex zone with widespread indications of heterogeneity on many scales, discontinuities of variable character, and shear-wave anisotropy. Discordance between P- and S-wave speed anomalies suggests the presence of chemical heterogeneity rather than just the effect of temperature.     

Cold plumes trigger contamination of oceanic mantle wedges with continental crust-derived sediments: Evidence from chromitite zircon grains of eastern Cuban ophiolites

The origin of zircon grains, and other exotic minerals of typical crustal origin, in mantle-hosted ophiolitic chromitites are hotly debated. We report a population of zircon grains with ages ranging from Cretaceous (99 Ma) to Neoarchean (2750 Ma), separated from massive chromitite bodies hosted in the mantle section of the supra-subduction (SSZ)-type Mayarí-Baracoa Ophiolitic Belt in eastern Cuba. Most analyzed zircon grains (n = 20, 287 ± 3 Ma to 2750 ± 60 Ma) are older than the early Cretaceous age of the ophiolite body, show negative ε_Hf(t) (−26 to −0.6) and occasional inclusions of quartz, K-feldspar, biotite, and apatite that indicate derivation from a granitic continental crust. In contrast, 5 mainly rounded zircon grains (297 ± 5 Ma to 2126 ± 27 Ma) show positive ε_Hf(t) (+0.7 to +13.5) and occasional apatite inclusions, suggesting their possible crystallization from melts derived from juvenile (mantle) sources. Interestingly, younger zircon grains are mainly euhedral to subhedral crystals, whereas older zircon grains are predominantly rounded grains. A comparison of the ages and Hf isotopic compositions of the zircon grains with those of nearby exposed crustal terranes suggest that chromitite zircon grains are similar to those reported from terranes of Mexico and northern South America. Hence, chromitite zircon grains are interpreted as sedimentary-derived xenocrystic grains that were delivered into the mantle wedge beneath the Greater Antilles intra-oceanic volcanic arc by metasomatic fluids/melts during subduction processes. Thus, continental crust recycling by subduction could explain all populations of old xenocrystic zircon in Cretaceous mantle-hosted chromitites from eastern Cuba ophiolite. We integrate the results of this study with petrological-thermomechanical modeling and existing geodynamic models to propose that ancient zircon xenocrysts, with a wide spectrum of ages and Hf isotopic compositions, can be transferred to the mantle wedge above subducting slabs by cold plumes.     

Asymmetry in oceanic crustal structure of the South China Sea basin and its implications on mantle geodynamics

ABSTRACT

We investigated the oceanic crustal structure and lithospheric dynamics of the South China Sea (SCS) basin through a comprehensive analysis of residual gravity anomaly and bathymetry combined with seismic constraints and interpretation from geodynamic modelling. We first calculated the residual mantle Bouguer anomaly (RMBA) of the oceanic crustal regions of the SCS by removing from free-air gravity anomaly the predicted gravitational attractions of water-sediment, sediment-crust, and crust-mantle interfaces, as well as the effects of lithospheric plate cooling, using the latest crustal age constraints including IODP Expedition 349 and recent deep-tow magnetic surveys. We then calculated models of the gravity-derived crustal thickness and calibrated them using the available seismic refraction profiles of the SCS. The gravity-derived crustal thickness models correlate positively with seismically determined crustal thickness values. Our analysis revealed that the isochron-averaged RMBA are consistently more negative over the northern flank of the SCS basin than the southern conjugate for magnetic anomaly chrons C8n (~25.18 Ma) to C5Dn (~17.38 Ma), implying warmer mantle and/or thicker crust over much of the northern flank. Computational geodynamic modelling yielded the following interpretations: (1) Models of asymmetric and variable spreading rates based on the relatively high-resolution deep-tow magnetic analysis would predict alternating thicker and thinner crust at the northern flank than the southern conjugate, which is inconsistent with the observed systematically thicker crust on the northern flank. (2) Models of episodic southward ridge jumps could reproduce the observed N-S asymmetry, but only for crustal age of 23.6–20 Ma. (3) Southward migration of the SCS ridge axis would predict slightly thinner crust at the northern flank, which is inconsistent with the observations. (4) Models of higher mantle temperatures of up to 25–50°C or >2% less depleted mantle sources on the northern flank could produce large enough anomalies to explain the observed N-S asymmetries.     

Komsomolskaya diamondiferous eclogites: evidence for oceanic crustal protoliths

The Komsomolskaya kimberlite is one of numerous (>1,000) kimberlite pipes that host eclogite xenoliths on the Siberian craton. Eclogite xenoliths from the adjacent Udachnaya kimberlite pipe have previously been geochemically well characterized; however, data from surrounding diamond-bearing kimberlite pipes from the center of the craton are relatively sparse. Here, we report major- and trace-element data, as well as oxygen isotope systematics, for mineral separates of diamondiferous eclogite xenoliths from the Komsomolskaya kimberlite, suggesting two distinct subgroups of a metamorphosed, subducted oceanic crustal protolith. Using almandine contents, this suite can be divided into two subgroups: group B1, with a high almandine component (>20 mol%) and group B2, with a low almandine component (<20 mol%). Reconstructed REE profiles for B1 eclogites overlap with typical oceanic basalts and lack distinct Eu anomalies. In addition, elevated oxygen isotope values, which are interpreted to reflect isotopic exchange with seawater at low temperatures (<350 °C), are consistent with an upper-oceanic crustal protolith. Reconstructed REE profiles for B2 eclogites are consistent with oceanic gabbros and display distinct Eu anomalies, suggesting a plagioclase-rich cumulate protolith. In contrast to B1, B2 eclogites do not display elevated oxygen isotope values, suggesting an origin deep within the crustal pile, where little-to-no interaction with hydrothermal fluids has occurred. Major-element systematics were reconstructed based on mineral modes; group B1 eclogites have higher MgO wt% and lower SiO₂ wt%, with respect to typical oceanic basalts, reflecting a partial melting event during slab subduction. Calculated residues from batch partial melt modeling of a range of Precambrian basalts overlap with group B1 trace-element chemistry. When taken together with the respective partial melt trajectories, these melting events are clearly linked to the formation of Tonalite–Trondhjemite–Granodiorite (TTG) complexes. As a result, we propose that many, if not all, diamondiferous eclogite xenoliths from Komsomolskaya represent mantle ‘restites’ that preserve chemical signatures of Precambrian oceanic crust.     

Mapping crustal structure beneath southern Tibet: Seismic evidence for continental crustal underthrusting

Receiver function imaging along a temporary seismic array (ANTILOPE-2) reveals detailed information of the underthrusting of the Indian crust in southern Tibet. The Moho dips northward from ~ 50 km to 80 km beneath the Himalaya terrane, and locally reaches ~ 85 km beneath the Indus–Yalung suture. It remains at ~ 80 km depth across the Lhasa terrane, and shallows to ~ 70 km depth under the Qiangtang terrane. An intra-crustal interface at ~ 60 km beneath the Lhasa terrane can be clearly followed southward through the Main Himalaya Thrust and connects the Main Boundary Thrust at the surface, which represents the border of the Indian crust that is underthrusting until south of the Bangong–Nujiang Suture. A mid-crustal low velocity zone is observed at depths of 14–30 km beneath the Lhasa and Himalaya terranes probably formed by partial melt and/or aqueous fluids.     

Geochemical constraints on depth of origin of oceanic carbonatites: The Cape Verde case

We present new Sr-Nd isotope compositions together with major- and trace element concentrations measured for whole rocks and mineral separate phases (apatite, biotite and calcite) from fifteen Cape Verde oceanic carbonatites (Atlantic Ocean). Trace element patterns of calcio- and magnesio-carbonatites present a strong depletion in K, Hf, Zr and Ti and an overall enrichment in Sr and REE relative to Cape Verde basalts, arguing for distinct source components between carbonatites and basalts. Sr and Nd isotopic ratios show small, but significant variations defining a binary mixing between a depleted end-member with unradiogenic Sr and radiogenic Nd values and a ‘‘enriched’’ end-member compatible with old marine carbonates. We interpret the depleted end-member as the Cape Verde oceanic lithosphere by comparison with previous studies on Cape Verde basalts. We thus propose that oceanic carbonatites are resulting from the interaction of a deep rooted mantle plume carrying a lower ⁴He/³He signature from the lower mantle and a carbonated metasomatized lithosphere, which by low degree melting produced carbonatite magmas. Sr-Nd compositions and trace element patterns of carbonatites argue in favor of a metasomatic agent originating from partial melting of recycled, carbonated oceanic crust. We have successfully reproduced the main geochemical features of this model using a Monte-Carlo-type simulation.     

Constraints on crustal and mantle structure of the oceanic plate south of Iceland from ocean bottom recorded Rayleigh waves

From April to July 2002 we carried out a deployment of 6 ocean bottom seismometers and 4 ocean bottom hydrophones in the North Atlantic south of Iceland. During the deployment period we recorded clear Rayleigh waves from 2 regional and 14 teleseismic earthquakes. This corresponds to a Rayleigh wave detection rate of nearly 92% for events with M_W ≥ 6.06.0 and epicentral distance less than 110°, close to detection rate estimates based on noise level variability. We measured Rayleigh wave event-station group dispersion and inter-station phase dispersion for one Mid-Atlantic Ridge event. The group dispersion curve is sensitive to the structure of the North-East Atlantic with an average age of  39 Myr. The phase dispersion curve is sensitive to the structure just south of Iceland (average plate age 33 Myr). Both dispersion curves indicate faster velocities than previously postulated for oceanic plate generated at the Reykjanes Ridge. A grid search approach was used to constrain the range of models fitting the data. The high velocity seismic lid just south of Iceland in the model for the phase dispersion path is slower or thinner than in the group dispersion model, which averages over a larger area and a somewhat older plate age, but the velocities in the low velocity half space are similar. We further consider the residual bathymetry in the experimental area. The residual anomaly decreases by 300–400 m from the Reykjanes Ridge to the  30 Myr old plate south of Iceland. This decrease can be explained by the disappearance of a mantle thermal anomaly associated with the Iceland plume. Both the residual bathymetry and the surface wave data are thus consistent with the notion that the southward spreading of the Icelandic plume is channelised underneath the Reykjanes Ridge and does not spread far outside this channel.Based on the experience from the pilot experiment, we estimate that a minimum recording time of 13–15 months in favourable weather conditions (April–September) is required to record enough data to map the spreading plume with surface waves, and to produce a tomographic image to a depth of 1000 km using body waves. This can be achieved by a continuous deployment of at least  20 months, or by two or three deployments during the spring and summer of consecutive years.     

Accretion of oceanic plateaus at continental margins: Numerical modeling

The accretion of oceanic plateaus has played a significant role in continental growth during Earth's history, which is evidenced by the presence of oceanic island basalts (OIB) and plume-type ophiolites in many modern orogens. However, oceanic plateaus can also be subducted into the deeper mantle, as revealed by seismic tomography. The controlling factors of accretion versus subduction of oceanic plateaus remain unclear. Here, we investigate the dynamics of oceanic plateau accretion at active continental margins using a thermo-mechanical numerical model. Three major factors for the accretion of oceanic plateaus are studied: (1) a thinned continental margin of the overriding plate, (2) “weak” layers in the oceanic lithosphere, and (3) a young oceanic plateau. For a large oceanic plateau, the modes of oceanic plateau accretion can be classified into one-sided and two-sided subduction–collisional regimes, which mainly depend on the geometry of the continental margin (normal or thinned). For smaller-sized seamounts, accretion occurs only if all three factors are satisfied, of which a thinned continental margin is the most critical. Possible geological analogues for the two-sided subduction–collisional mode include the Taiwan orogenic belt and subduction of the Ontong Java Plateau. The accretion model for small oceanic plateaus applies to the Nadanhada Terrane in Northeast China.     

Tectonics and crustal structure of the Campania continental margin: relationships with volcanism

Summary ¶The crustal structure of the Campania continental margin is synthesized from outcrop, seismic reflection and gravimetric data. Outcrop and subsurface geological data reveal the presence of NE–SW faults, E–W faults and NW–SE faults. An older extensional event occurred along NW–SE faults and was followed by the main extensional event linked to the activity of NE–SW normal faults. The latter were active between 700 and 400ka producing half-grabens filled by more than 5km of Quaternary deposits. The stratigraphic signature of these tectonic events corresponds to a Lower Pleistocene marine unconformity-bounded unit overlain by Middle Pleistocene rocks belonging to a transgressive-regressive cycle. A crustal section of the Campania margin displays an asymmetric linked fault system characterized by a 10–12km-deep main detachment level, listric normal faults and rollover anticlines. Structural and stratigraphic data document that the inception of volcanic activity at Vesuvius occurred at 400ka, just after the main extensional event, and the volcano is located at the margin of a rollover anticline.Received June 26, 2002; revised version accepted November 9, 2002     

滇西新生代超钾质煌斑岩的元素和Sr-Nd同位素特征及其对岩石圈地幔组成的制约

系统分析了滇西早第三纪晚期煌斑岩的主元素、微量元素和Sr-Nd同位素组成。绝大多数煌斑岩属超钾质系列岩石（K2O/Na2O=2.1-5.2)，少数为钾玄质系列（K2O/Na2O=1.6-1.7)。这些超钾质煌斑岩富集大离子亲石元素、轻稀土元素和Pb，亏损高场强元素，具有高的初始^87Sr/^86Sr比值（O.70624-0.70924)和负的εNd(t)值（-1.7～-4.6），类似于与俯冲环境有关的高K/Ti-低Ti钾质系列岩石。这些超钾质煌斑岩母岩浆来源于含金云母的交代地幔，经历了不同程度的结晶分异和地壳物质的同化混染。与藏北钾玄质系列岩石相比，滇西超钾质煌斑岩的Th/U比较低、Rb/Sr比较高，而且Nd模式年龄系统偏低0.1-0.4Ga，表明青藏高原北部和东南部岩石圈地幔组成存在区域性的差异。