Formation time of the big mantle wedge beneath eastern China and a new lithospheric thinning mechanism of the North China craton—Geodynamic effects of deep recycled carbon

High-resolution P wave tomography shows that the subducting Pacific slab is stagnant in the mantle transition zone and forms a big mantle wedge beneath eastern China. The Mg isotopic investigation of large numbers of mantle-derived volcanic rocks from eastern China has revealed that carbonates carried by the subducted slab have been recycled into the upper mantle and formed carbonated peridotite overlying the mantle transition zone, which becomes the sources of various basalts. These basalts display light Mg isotopic compositions(δ26 Mg = –0.60‰ to –0.30‰) and relatively low87 Sr/86 Sr ratios(0.70314–0.70564) with ages ranging from 106 Ma to Quaternary, suggesting that their mantle source had been hybridized by recycled magnesite with minor dolomite and their initial melting occurred at 300-360 km in depth. Therefore, the carbonate metasomatism of their mantle source should have occurred at the depth larger than 360 km, which means that the subducted slab should be stagnant in the mantle transition zone forming the big mantle wedge before 106 Ma. This timing supports the rollback model of subducting slab to form the big mantle wedge. Based on high P-T experiment results, when carbonated silicate melts produced by partial melting of carbonated peridotite was raising and reached the bottom(180–120 km in depth) of cratonic lithosphere in North China, the carbonated silicate melts should have 25–18 wt% CO2 contents, with lower Si O2 and Al2 O3 contents, and higher Ca O/Al2 O3 values, similar to those of nephelinites and basanites, and have higher εNdvalues(2 to 6). The carbonatited silicate melts migrated upward and metasomatized the overlying lithospheric mantle, resulting in carbonated peridotite in the bottom of continental lithosphere beneath eastern China. As the craton lithospheric geotherm intersects the solidus of carbonated peridotite at 130 km in depth, the carbonated peridotite in the bottom of cratonic lithosphere should be partially melted, thus its physical characters are similar to the asthenosphere and it could be easily replaced by convective mantle. The newly formed carbonated silicate melts will migrate upward and metasomatize the overlying lithospheric mantle. Similarly, such metasomatism and partial melting processes repeat, and as a result the cratonic lithosphere in North China would be thinning and the carbonated silicate partial melts will be transformed to high-Si O2 alkali basalts with lower εNdvalues(to-2). As the lithospheric thinning goes on,initial melting depth of carbonated peridotite must decrease from 130 km to close 70 km, because the craton geotherm changed to approach oceanic lithosphere geotherm along with lithospheric thinning of the North China craton. Consequently, the interaction between carbonated silicate melt and cratonic lithosphere is a possible mechanism for lithosphere thinning of the North China craton during the late Cretaceous and Cenozoic. Based on the age statistics of low δ26 Mg basalts in eastern China, the lithospheric thinning processes caused by carbonated metasomatism and partial melting in eastern China are limited in a timespan from 106 to25 Ma, but increased quickly after 25 Ma. Therefore, there are two peak times for the lithospheric thinning of the North China craton: the first peak in 135-115 Ma simultaneously with the cratonic destruction, and the second peak caused by interaction between carbonated silicate melt and lithosphere mainly after 25 Ma. The later decreased the lithospheric thickness to about70 km in the eastern part of North China craton.     

The effect of mantle composition on density in the extending lithosphere

The density distribution of the lithosphere is non-linear and discontinuous due to complex mineralogy and, most importantly, phase transitions. We evaluate the influence of changes in mantle composition on lithospheric density and its evolution during horizontal stretching, using thermodynamic calculations of the density as a function of pressure, temperature and composition. We also develop a simple parameterization based on end-member mineral reactions and geometric relationships between the geotherm and the phase boundary for comparison. The garnet–spinel peridotite transition leads to a moderate decrease in density of the mantle part of the lithospheric column at the initial stages of stretching. When the crust is sufficiently thinned and temperature is relatively high, plagioclase peridotite becomes stable in the upper part of the mantle. The density reduction due to the plagioclase-in reaction is controlled by bulk Al₂O₃ in the mantle and by the depth of the plagioclase-in reaction, which is mainly governed by the Na₂O/Al₂O₃ ratio. Since Na₂O and Al₂O₃ increase with the fertility of the mantle the phase transition effect is most pronounced for relatively fertile mantle (and strong extension) and can lead to 2.3% density reduction. This is equivalent to heating the entire lithosphere by 700 °C if only the effect of thermal expansion on density is taken into account. The formation of plagioclase peridotite can explain syn-rift uplift in sedimentary basins that experienced large mantle stretching without invoking an unrealistically strong increase in temperature. It might also be responsible for the break-up unconformity observed at continental margins.     

Tectonic modelling in the Bjørnøya West Basin of the Western Barents Sea1

The Bjøirnøya West Basin lies between latitudes 73° and 74°, longitudes 16°E and 18°E, contains at least 8 km of sediments deposited from the Late Jurassic, and is of considerable interest for hydrocarbon exploration. The Cenozoic extensional tectonics in the basin can be clearly seen from seismic data with normal faulting and from subsidence curves with rapid subsidence. The extension occurred during the Late Palaeocene with active extension lasting about 6 million years (m.y.) followed by thermal cooling. The tectonic subsidence within the study area shows a three-phase development: phase 1, synrift (58–52 Ma (million years before the present day)), is characterized by rapid subsidence; phase 2, postrift (52–5 Ma), by slow subsidence with occasional uplift; and phase 3 (5–0 Ma), by rapid subsidence. An adaptive finite-element model, with consideration of the radiogenic heat production in the lithosphere, has been used to model the subsidence and heat flow. The modelling of subsidence shows the β-factor distribution varying from 1.9 to 3.5 with an average of 2.4 for the uniform lithospheric extension. The heat-flow modelling predicts a rapid increase of heat flow during the Early Palaeocene. The maximum heat flow at about 52 Ma, which could be as much as 3.0 hfu (10^?6 cal/cm²/s), was followed by a decrease in heat flow. A plate-weakening model has been proposed to explain the rapid subsidence for the last 5 m.y. by flexure of the elastic lithosphere which is weakened by a decrease in elastic thickness caused by an increase of the temperature gradient in the lithosphere. The plate-weakening model predicts a heat-flow increase at 5 Ma of up to 2.0 hfu. Our study, using quantitative modelling of the tectonic subsidence, provides a partial (if not a full) understanding of the tectonic development and thermal evolution of the Bjønøya West Basin.     

The effective elastic lithosphere under the Cook-Austral and Society islands

We have determined the elastic thicknessT_e of the oceanic lithosphere along two volcanic chains of the South Central Pacific: Cook-Austral and Society islands. We used a three-dimensional spatial method to model the lithospheric flexure assuming a continuous elastic plate. The model was constrained by geoid height data from the SEASAT satellite.Along the Cook-Austral chain the elastic thickness increases westward, from 2–4 km at McDonald hot spot to 14 km at Rarotonga. At McDonald seamount, however, the data are better explained by a local compensation model. The observed trend shows an increase ofT_e with age of plate at loading time. However, the elastic layer under the Cook-Austral appears systematically thinner by several kilometers than expected for “normal” seafloor, suggesting that substantial thermal thinning has taken place in this region. Considering the apparent thermal age of the plate instead of crustal age improves noticeably the results. Along the Society chainT_e varies from 20 km under Tahiti to 13 km under Maupiti which is located 500 km westward. When plotting together the Society and Cook-AustralT_e results versus age of load, we notice that within the first five million years after loading,T_e decreases significantly while tending rapidly to an equilibrium value. This may be interpreted as the effect of initial stress relaxation which occurs just after loading inside the lower lithosphere and suggests that the presently measured elastic thickness under the very young Tahiti load ( 0.8 Ma) is not yet the equilibrium thickness.     

Thermal regime transition in eastern North China and its tectonic implication

Many evidences published in recent years reveal that the thickness,chemical composition and thermal state of lithosphere in eastern North China have ex-perienced dramatic transition during the Phanero-zoic.Comparative study[1—8]of early Paleozoic dia-mond-and xenoliths-bearing kimberlites to Cenozoic mantle peridotite xenoliths-bearing alkali basalts indi-cates that the early Phanerozoic lithosphere is thick and stable to depths within the diamond stability field(180–200km)with depleted ma…     

Rapid early Miocene acceleration of uplift in the Gangdese Belt, Xizang (southern Tibet), and its bearing on accommodation mechanisms of the India-Asia collision

Five samples from a biotite-hornblende granodiorite phase of the 42.5 Ma Quxu pluton, Gangdese batholith, southern Tibet, have been collected at 250 m vertical intervals. Biotite from these rocks yields monotonically decreasing⁴⁰Ar/³⁹Ar isochron ages with decreasing elevation of 26.8 ± 0.2, 23.3 ± 0.5, 19.7 ± 0.3, 18.4 ± 0.4,and17.8 ± 0.1Ma (T_c = 335°C). Coexisting K-feldspars have virtually identical minimum apparent⁴⁰Ar/³⁹Ar ages of 17.0 ± 0.4Ma (T_c = 285°C). These data indicate parts of southern Tibet experienced a pulse of uplift in the early Miocene with the rate of uplift rising from 0.07 to 4.4 mm/year in the interval 20 to 17 Ma. An apatite fission track age of 9.9 ± 0.9Ma from this locality constrains the average uplift rate at this site to about 0.81 mm/year between 17 and 9.9 Ma and 0.30 mm/year from 9.9 Ma to present. K-feldspar from the Dagze granite, 30 km to the east, near Lhasa, yields a minimum apparent⁴⁰Ar/³⁹Ar age of 35.9 ± 0.9Ma (T_c = 227°C) which indicates an average uplift rate there of 0.21 mm/year since then. The marked pulse of uplift of the Quxu granodiorite and the difference in uplift history between the Dagze and Quxu plutons suggests southern Tibet has experienced discrete pulses of uplift variable in both space and time. These data are not consistent with models which require a large proportion of uplift of the Tibetan plateau to have occurred in the last 2 Ma. The data support the suggestion that convergence between India and Asia was largely accommodated by tectonic escape during the opening of the South China Sea 32 to 17 Ma ago and permit distributed shortening as a mechanism for crustal thickening and uplift of this part of the Tibetan plateau subsequent to 20 Ma.     

Ancient Os isotope signatures from the Ontong Java Plateau lithosphere: Tracing lithospheric accretion history

In order to better understand the nature and formation of oceanic lithosphere beneath the Early Cretaceous Ontong Java Plateau, Re–Os isotopes have been analysed in a suite of peridotite xenoliths from Malaita, Solomon Islands. Geological, thermobarometric and petrological evidence from previous studies reveal that the xenoliths represent virtually the entire thickness of the southern part of subplateau lithospheric mantle (< 120 km). This study demonstrates that vertical Os isotopic variations correlate with compositional variations in a stratified lithosphere. The shallowest plateau lithosphere (< 85 km) is dominated by fertile lherzolites showing a restricted range of ¹⁸⁷Os/¹⁸⁸Os (0.1222 to 0.1288), consistent with an origin from ~ 160 Ma Pacific lithosphere. In contrast, the basal section of subplateau lithospheric mantle (~ 95–120 km) is enriched in refractory harzburgites with highly unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios ranging from 0.1152 to 0.1196, which yield Proterozoic model ages of 0.9–1.7 Ga. Although the whole range of Os isotope compositions of Malaita peridotites is within the variations seen in modern abyssal peridotites, the contrasting isotopic compositions of shallow and deep plateau lithosphere suggest their derivation from different mantle reservoirs. We propose that the subplateau lithosphere forms a genetically unrelated two-layered structure, comprising shallower, typical oceanic lithosphere underpinned by deeper impinged material, which included a component of recycled Proterozoic lithosphere. The impingement of residual but chemically heterogeneous mantle, mechanically coupled to the recently formed, thin lithosphere, may have a bearing on the anomalous initial uplift and late subsidence history of the seismically anomalous plateau root.     

Characteristics of subcontinental lithospheric mantle beneath Baegryeong Island, Korea: Spinel peridotite xenoliths

Abstract   Geochemical data on Baegryeong Island spinel peridotites found in Miocene alkali basalt provide the information for lithosphere composition, chemical processes, equilibrium pressure and temperature conditions. Spinel peridotite xenoliths, showing transitional textures between protogranular and porpyroclastic textures, were accidentally trapped by the ascending alkali basalt magma. The xenoliths originate at depths from 50 to 70 km with a temperature range from 800 to 1100°C. The variations of modal and mineral compositions of the spinel peridotite xenoliths indicate that the xenoliths have undergone 1–10% fractional melting. The spinel peridotites from Baegryeong Island have undergone cryptic mantle metasomatism subsequent to melt extraction. Metasomatic agent of enriched spinel peridotite xenoliths was carbonatite melt.     

Evolution of lithospheric mantle in the north of Nain‐Baft oceanic crust (Neo‐Tethyan ophiolite of Ashin,Central Iran)

This study is focused on a plagioclase‐bearing spinel lherzolite from Chah Loqeh area in the Neo‐Tethyan Ashin ophiolite. It is exposed along the west of left‐lateral strike‐slip Dorouneh Fault in the northwest of Central‐East Iranian Microcontinent. Mineral chemistry (Mg#^olivine < ~ 90, Cr#^{clinopyroxene} < ~ 0.2, Cr#^spinel < ~ 0.5, Al₂O₃^{orthopyroxene} > ~ 2.5 wt%, Al₂O₃^{clinopyroxene} > ~ 4.5 wt%, Al₂O₃^spinel > ~ 41.5 wt%, Na₂O^{clinopyroxene} > ~ 0.11 wt%, and TiO₂^{clinopyroxene} > ~ 0.04 wt%) confirms Ashin lherzolite was originally a mid‐oceanic ridge peridotite with low degrees of partial melting at spinel‐peridotite facies in a lithospheric mantle level. However, some Ashin lherzolites record mantle upwelling and tectonic exhumation at plagioclase‐peridotite facies during oceanic extension and diapiric motion of mantle along Nain‐Baft suture zone. This mantle upwelling is evidenced by some modifications in the modal composition (i.e. subsolidus recrystallization of plagioclase and olivine between pyroxene and spinel) and mineral chemistry (e.g. increase in TiO₂ and Na₂O of clinopyroxene, and TiO₂ and Cr# of spinel and decrease in Mg# of olivine), as a consequence of decompression during a progressive upwelling of mantle. Previous geochronological and geochemical data and increasing the depth of subsolidus plagioclase formation at plagioclase‐peridotite facies from Nain ophiolite (~ 16 km) to Ashin ophiolite (~ 35 km) suggest a south to north closure for the Nain‐Baft oceanic crust in the northwest of Central‐East Iranian Microcontinent.     

Diapirism of depleted peridotite — a model for the origin of hot spots

A model is proposed for the origin of hot spots that depends on the existence of major-element heterogeneities in the mantle. Generation of basaltic crust at spreading centers produces a layer of residual peridotite ～20–25 km thick directly beneath the crust which is depleted in Fe/Mg, TiO₂, CaO, Al₂O₃, Na₂O and K₂O, and which has a slightly lower density than undepleted peridotite beneath it. Upon recycling of this depleted peridotite back into the deep mantle at subduction zones, it becomes gravitationally unstable, and tends to rise as diapirs through undepleted peridotite. For a density contrast of 0.05 g cm^?3, a diapir 60 km in diameter would rise at roughly 8 cm y^?1, and could transport enough heat to the base of the lithosphere to cause melting and volcanism at the surface. Hot spots are thus viewed as a passive consequence of mantle convection and fractionation at spreading centers rather than a plate-driving force.It is suggested that depleted diapirs exist with varying amounts of depletion, diameters, upward velocities and source volumes. Such variations could explain the occurrence of hot spots with widely varying lifetimes and rates of lava production. For highly depleted diapirs with very low Fe/Mg, the diapir would act as a heat source and the asthenosphere and lower lithosphere drifting across the diapir would serve as the source region of magmas erupted at the surface. For mildly depleted diapirs with Fe/Mg only slightly less than in normal undepleted mantle, the diapir could provide not only the source of heat but also most or all of the source material for the erupted magmas. The model is consistent with isotopic data that require two separate and ancient source regions for mid-ocean ridge and oceanic island basalts. The source for mid-ocean ridge basalts is considered to be material upwelling at spreading centers from the deep mantle. This material forms the oceanic lithosphere. Oceanic island basalts are considered to be derived from varying mixtures of sublithospheric and lower lithospheric material and the rising diapir itself.     

Geodynamic evolution of a forearc rift in the southernmost Mariana Arc

The southernmost Mariana forearc stretched to accommodate opening of the Mariana Trough backarc basin in late Neogene time, erupting basalts at 3.7–2.7 Ma that are now exposed in the Southeast Mariana Forearc Rift (SEMFR). Today, SEMFR is a broad zone of extension that formed on hydrated, forearc lithosphere and overlies the shallow subducting slab (slab depth ≤ 30–50 km). It comprises NW–SE trending subparallel deeps, 3–16 km wide, that can be traced ≥ ∼30 km from the trench almost to the backarc spreading center, the Malaguana‐Gadao Ridge (MGR). While forearcs are usually underlain by serpentinized harzburgites too cold to melt, SEMFR crust is mostly composed of Pliocene, low‐K basaltic to basaltic andesite lavas that are compositionally similar to arc lavas and backarc basin (BAB) lavas, and thus defines a forearc region that recently witnessed abundant igneous activity in the form of seafloor spreading. SEMFR igneous rocks have low Na₈, Ti₈, and Fe₈, consistent with extensive melting, at ∼23 ± 6.6 km depth and 1239 ± 40°C, by adiabatic decompression of depleted asthenospheric mantle metasomatized by slab‐derived fluids. Stretching of pre‐existing forearc lithosphere allowed BAB‐like mantle to flow along the SEMFR and melt, forming new oceanic crust. Melts interacted with pre‐existing forearc lithosphere during ascent. The SEMFR is no longer magmatically active and post‐magmatic tectonic activity dominates the rift.     

Flexural uplift of a lithospheric slab near the Vema transform (Central Atlantic): Timing and mechanisms

The Vema Transverse Ridge (VTR) is a prominent, long and narrow topographic anomaly that runs for over 300 km along a sea floor spreading flow line south of the Vema transform at 11° N in the Atlantic. It rises abruptly about 140 km from the axis of the Mid-Atlantic Ridge (MAR) in 10 Myr old crust and runs continuously up to 25 Myr old crust. It reaches over 3 km above the predicted lithospheric thermal contraction level. It is absent in crust younger than 10 Myr; thus, the uplift of the VTR must have ended roughly 10 Ma. The VTR is interpreted as the exposed edge of a flexured and uplifted slab of oceanic lithosphere that was generated at an 80 km long MAR segment. Based on satellite gravimetry imagery this MAR segment was born roughly 50 Ma and increased its length at an average rate of 1.6 mm/yr. Multibeam data show that the MAR-parallel sea floor fabric south of the VTR shifts its orientation by 5° to 10° clockwise in 11–12 Myr old crust, indicating a change at that time of the orientation of the MAR axis and of the position of the Euler rotation pole. This change caused extension normal to the transform, followed between 12 and 10 Ma by flexure of the edge of the lithospheric slab, uplift of the VTR at a rate of 2 to 4 mm/yr, and exposure of a lithospheric section (Vema Lithospheric Section or VLS) at the northern edge of the slab, parallel to the Vema transform. Ages of pelagic carbonates encrusting ultramafic rocks sampled at the base of the VLS at different distances from the MAR axis suggest that the entire VTR rose vertically as a single block within the active transform offset. A 50 km long portion of the crest of the VTR rose above sea level, subsided, was truncated at sea level and covered by a carbonate platform. Subaerial and submarine erosion has gradually removed material from the top of the VTR and has modified its slopes. Spreading half rate of the crust south of the transform decreased from 17.2 mm/yr between 26 and 19 Ma to 16.9 mm/yr between 19 and 10 Ma, to 13.6 mm/yr from 10 Ma to present. The slowing down of spreading occurred close in time to the change in ridge/transform geometry, suggesting that the two events are related. A numerical model relates lithospheric flexure to extension normal to the transform, suggesting that the extent of the uplift depends on the thickness of the brittle layer, consistent with the observed greater uplift of the older lithosphere along the VTR.     

Heat flow and depth versus age for the Mesozoic northwest Atlantic Ocean: results from the Sohm abyssal plain and implications for the Bermuda Rise

Heat flow in the Sohm abyssal plain is measured to be 53 mW/m² at an age of 163 Ma. This is 25% higher than predicted by conductive cooling models, even though the sediment-corrected basement depth of 6.5 km at this location is normal for its age. An analysis of existing heat flow, depth and geoid anomalies in the northwest Atlantic shows that there is little correlation between heat flow and depth throughout the entire region. Depth and geoid are clearly related to the Bermuda swell while the associated heat flow anomaly, once adjusted for variations with age, is limited to 5 mW/m² and only decays to the south. This means that the Bermuda swell is probably not caused by extensive thermal reheating within the lithosphere, but instead by dynamic uplift at its lower boundary due to the convective upwelling of a mantle plume. The regionally high heat flow in the northwest Atlantic may be a thermal remanent of previous plumes which passed beneath this region early in its history. Therefore, depth and heat flow anomalies from this region cannot be used to provide constraints on steady-state parameters of the lithosphere, such as the presence or absence of a long-term boundary layer at its base.     

Quantitative relationships between uplift and relief parameters for the Southern Alps,New Zealand,as determined by fission track analysis

The Southern Alps mountain chain, New Zealand, has formed as a consequence of late Cenozoic collision of the continental parts of the Pacific and Australia plates. Fission track analysis has yielded estimates of the amount, age of initiation, and rate of late Cenozoic rock uplift for 82 surface samples taken from transects across the Southern Alps. The mean surface, summit and valley elevations in the vicinity of each of the rock sample sites have also been measured. Regression of the geomorphic variables on the uplift variables has been used to establish quantitative relationships between uplift and geomorphology. There are strong and consistent linear associations between uplift and the elevations of the mean surface, summits and valleys. The preferred regression models have uniform slope but varying elevation response between transects. Substitution of space for time has allowed the evolution of landforms to be studied. To the east of the Main Divide, elevation and relief are proportional to, and closely related to, the age of initiation of rock uplift (‘uplift age’) and the amount of rock uplift (r² > 0·8). Mean surface uplift was delayed for ～2 Ma after the start of rock uplift, a result of the stripping of a soft cover rock succession that, prior to rock uplift, overlaid the harder greywacke basement. Inter-transect variations in regression response and x-intercept are inferred, therefore, to reflect the variable preuplift thickness of cover rocks. However, the regular regression slope for the transects reflects the consistent nature of the interaction between uplift and the erodibility of greywacke basement. Uplift of the mean surface proceeded at 0·4 km/km and 0·4 km/Ma of rock uplift, while the rock uplift rate was 0·8 km/Ma. Summit elevations have increased at a rate of 0·6 km/Ma and valley elevations have increased at 0·2 km/Ma. Regression lines relating mean surface, summit and valley elevations to rock uplift and uplift age diverge from common intercepts; it is concluded, therefore, that the mountains east of the Main Divide have continued to increase in elevation and relief and change in form over time since the start of mean surface uplift. Mountain elevation has little relationship with late Cenozoic mean rock uplift rates of 0·8–1·0 km/Ma or inferred contemporary rock uplift rates (r² ～ 0·3). In contrast, to the west of the Main Divide, elevation is shown to be closely related to rock uplift rate (r² > 0·3). In contrast, to the west of the Main Divide, elevation is shown to be closely related to rock uplift rate (r² > 0·8). Transects with higher rock uplift rates support higher topography. Landforms are therefore in a stable equilibrium with rock uplift rate, and the landscape contains no residual evidence of the total amount of rock uplift, or the age of uplift. Lithological variation appears to have no relationship with elevation.     

Petrology of the Green Knobs diatreme and implications for the upper mantle below the Colorado Plateau

Peridotite inclusions, crystal fragments, and kimberlite breccia at Green Knobs, New Mexico, have been studied to evaluate compositions and processes in the upper mantle below the Colorado Plateau. Most peridotite inclusions are spinel lherzolites and harzburgites, or their partly hydrated equivalents, in the Cr-diopside group. Orthopyroxene-rich websterites and olivine websterites comprise 3% of the peridotites and formed as cumulates. Typical anhydrous or slightly hydrated peridotites contain aluminous, calcic diopside (5–7% Al₂O₃), aluminous orthopyroxene (3–6% Al₂O₃), spinel, and olivine (near Fa₉). Geothermometers based on different mineral pairs yield temperatures from above 1100°C to below 700°C in single rocks. High values, derived from pyroxenes with included exsolution lamellae, may approximate temperatures of primary crystallization. Low values, based on olivine-spinel and olivine-clinopyroxene pairs, approach upper mantle temperatures before eruption. In rare samples, some spinel grains are rimmed by garnet while others are not rimmed; garnet formation was controlled by nucleation kinetics. About one-third of the peridotites were deformed shortly before eruption, with effects ranging from mild cataclasis to the production of ultramylonites.Discrete crystals of garnet, olivine (near Fa₈), and Cr-diopside represent garnet peridotite. Eclogites were not found. The garnet peridotite is more depleted than overlying spinel peridotite, and it is not a likely source for the minettes associated with the kimberlites.The mantle below Green Knobs consists of spinel peridotite from 45 to perhaps 60 km depth immediately underlain by more-depleted garnet peridotite. The position of the spinel-garnet transition may be fixed by kinetics. The kimberlite may have been produced when heat from ascending minette magma released volatiles from otherwise depleted garnet peridotite. Resulting gas-solid mixtures erupted along zones of deformation associated with Colorado Plateau monoclines. Sheared lherzolites formed during renewed movement along these zones.     

Thermal cooling of the oceanic lithosphere: new constraints from geoid height data

Up to now, tests of thermal models of the oceanic lithosphere as it cools and moves away from the ridge crest have been based mainly on topography and heat flow data. However, large areas of the ocean floor deviate from the normal subsidence due to thermal contraction and heat flow data are not very sensitive to the form of the model.

Cooling of the lithosphere causes a short-wavelength step in the geoid across fracture zones that can also be used to constrain thermal models. We have analyzed geoid data at fracture zones from the SEASAT altimeter measurements in the entire Pacific Ocean and redetermined parameters of the cooling models. We find that the data reveal two distinct regimes of cooling; one for seafloor ages in the range 0–30 Ma, the other beyond 30 Ma; this does not appear to be correlated with particular fracture zones but rather it is representative of the whole area studied, i.e., the entire south Pacific and northeast Pacific Ocean. These two trends may be interpreted in terms of two different (asymptotic) thermal thicknesses of the plate model. The smaller thermal thickness ( 65 km) found for ages <30 Ma—compared to 90 km in the age range 30–50 Ma—calls for some kind of thermal perturbation in the vicinity of the ridge crest.

From the results obtained in this study, we conclude that the half-space cooling model is unable to explain the data, that beyond 30 Ma, a simple plate model gives a satisfactory fit to the data but in the younger plate portion (ages < 30 Ma) the cooling history of the oceanic lithosphere is much more complex than predicted by the usual cooling models. Furthermore, the depth-age relationship obtained from the geoid-derived thermal parameters departs significantly beyond 30 Ma from the widely used Parsons and Sclater's depth-age curve, predicting a lesser subsidence.     

Crustal structure and origin of the Cape Verde Rise

The Cape Verde Islands are located on a mid-plate topographic swell and are thought to have formed above a deep mantle plume. Wide-angle seismic data have been used to determine the crustal and uppermost mantle structure along a ~ 440 km long transect of the archipelago. Modelling shows that ‘normal’ oceanic crust, ~ 7 km in thickness, exists between the islands and is gently flexed due to volcano loading. There is no direct evidence for high density bodies in the lower crust or for an anomalously low density upper mantle. The observed flexure and free-air gravity anomaly can be explained by volcano loading of a plate with an effective elastic thickness of 30 km and a load and infill density of 2600 kg m^− 3. The origin of the Cape Verde swell is poorly understood. An elastic thickness of 30 km is expected for the ~ 125 Ma old oceanic lithosphere beneath the islands, suggesting that the observed height of the swell and the elevated heat flow cannot be attributed to thermal reheating of the lithosphere. The lack of evidence for high densities and velocities in the lower crust and low densities and velocities in the upper mantle, suggests that neither a crustal underplate or a depleted swell root are the cause of the shallower than expected bathymetry and that, instead, the swell is supported by dynamic uplift associated with the underlying plume.     

Subduction and back-arc activity at the Hikurangi convergent Margin,New Zealand

The Hikurangi Margin is a region of oblique subduction with northwest-dipping intermediate depth seismicity extending southwest from the Kermadec system to about 42°S. The current episode of subduction is at least 16–20 Ma old. The plate convergence rate varies along the margin from about 60 mm/a at the south end of the Kermadec Trench to about 45 mm/a at 42°S. The age of the Pacific lithosphere adjacent to the Hikurangi Trench is not known.The margin divides at about latitude 39°S into two quite dissimilar parts. The northern part has experienced andesitic volcanism for about 18 Ma, and back-arc extension in the last 4 Ma that has produced a back-arc basin onshore with high heaflow, thin crust and low upper-mantle seismic velocities. The extension appears to have arisen from a seawards migration of the Hikurangi Trench north of 39°S. Here the plate interface is thought to be currently uncoupled, as geodetic data indicate extension of the fore-arc basin, and historic earthquakes have not exceededM _s=7.South of 39°S there is no volcanism and a back-arc basin has been produced by downward flexure of the lithosphere due to strong coupling with the subducting plate. Heatflow in the basin is normal. Evidence for strong coupling comes from historic earthquakes of up to aboutM _s=8 and high rates of uplift on the southeast coast of the North Island.The reason for this division of the margin is not known but may be related to an inferred increase, from northeast to southwest, in the buoyancy of the Pacific lithosphere.     

Thermal-rheological structure of lithosphere beneath the northern flank of Tarim Basin, western China: Implications for geodynamics

Based on the data of geo-temperature and thermophysical parameters of rocks in the Kuqa Depression and the Tabei Uplift, northern flank of the Tarim Basin, in terms of the analytical solution of 1-D heat transfer equation, the thermal structure of the lithosphere under this region is determined. Our results show that the average surface heat flow of the northern flank of the Tarim Basin is 45 mW/m2, and the mantle heat flow is between 20 and 23 mW/m2; the temperature at crust-mantle boundary (Moho) ranges from 514℃ to 603℃ and the thermal lithosphere where the heat conduction dominates is 138-182 km thick. Furthermore, in combination with the P wave velocity structure resulting from the deep seismic sounding profile across this region and rheological modeling, we have studied the local composition of the lithosphere and its rheological profile, as well as the strength distribution. We find that the rheological stratification of the lithosphere in this region is apparent. The lowermost of the lower crust is ductile; however,the uppermost of the mantle and the upper and middle parts of the crust are both brittle layers,which is typically the so-called sandwich-like structure. Lithospheric strength is also characterized by the lateral variation, and the uplift region is stronger than the depression region. The lithospheric strength of the northem flank of the Tarim Basin decreases gradually from south to north; the Kuqa Depression has the lowest strength and the south of the Tabei Uplift is strongest.The total lithospheric strength of this region is 4.77× 1012-5.03 × 1013 N/m under extension, and 6.5 × 1012-9.4× 1013 N/m under compression. The lithospheric brittle-ductile transition depth is between 20 km and 33 km. In conclusion, the lithosphere of the northern flank of the Tarim Basin is relatively cold with higher strength, so it behaves rigidly and deforms as a whole, which is also supported by the seismic activity in this region. This rigidity of the Tarim lithosphere makes it little deform interior, but only into flexure under the sedimentation and tectonic loading associated with the rapid uplift of the Tianshan at its northern margin during the Indian-Eurasian continental collision following the Late Eocene. Finally, the influences of factors, such as heat flow, temperature,crustal thickness, and especially basin sediment thickness, on the lithospheric strength are discussed here.     

Rounded mineral inclusions in upper mantle peridotite

Mantle-derived peridotite xenoliths contain abundant olivine with rounded spinel and orthopyroxene inclusions, and orthopyroxene with rounded olivine and spinel inclusions. The shape-change of spinel, olivine and orthopyroxene inclusions from the primarily polyhedral outline to the spherical outline is governed by interfacial diffusion of oxygen in spinel (the most sluggish atom) due to reduction of total interfacial energy of the host-inclusion system.The critical radius of maximum rounded inclusions of spinel in olivine is a function of temperature and annealing time. Assuming that the activation energy for the interfacial diffusion is 40–70 kcal mol^?1 and that the annealing time for the spinel lherzolite from Salt Lake Crater in Hawaii is 100 Ma, the annealing times for perioditites under the island arcs of Japan are estimated to be 1–10 Ma.