Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

Chemical composition of altered volcanic rocks from ODP Sites 193-1188 and 193-1189

The conventional model of leaching volcanic rocks as a source of metals in a seafloor hydrothermal systems has been tested by examining the behavior of Pb and other trace elements during hydrothermal alteration. ODP Leg 193 drill sites 1188 (Snowcap) and 1189 (Roman Ruins) on Pual Ridge in the eastern Manus Basin offshore eastern Papua New Guinea provide a unique three-dimensional window into an active back-arc hydrothermal system. We investigate by means of a LA-ICP-MS microbeam technique the capacity of Pb to be leached from a host volcanic rock exposed to various types and intensities of alteration. Our results are in general agreement with previous studies that utilized bulk analytical techniques but provide a more detailed explanation of the processes. Fresh representative dacitic lavas from the Pual Ridge have an average whole rock Pb content of 5.2 ppm, an average interstitial glass Pb content of 5.6 ppm and an average plagioclase Pb content of 1.0 ppm. Altered matrix samples have highly variable Pb values ranging from 0 to 52.4 ppm. High Pb values in altered samples are associated with a low temperature chlorite and clay mineral assemblage, in some cases overprinted by a high temperature (up to 350°C) silica-rich "bleaching" alteration. Only the most highly altered matrix samples have REE patterns that differ from the fresh Pual Ridge dacite. This may represent either different lava histories or alteration characteristics that have affected normally immobile REEs. Altered samples with the highest Pb values have similar REE patterns to those of the local unaltered lavas. They are compositionally similar to typical Pual Ridge dacites indicating a genetic relationship between the main regional volcanic suite and the subseafloor hydrothermally altered, Pb-enriched material. Relative loss/gain for Pb between the analyzed altered samples and a calculated precursor show a maximum relative gain of 901%. Samples with relative Pb gain from both drill sites are associated with lower temperature alteration mineral assemblages characterized by pervasive chloritization. The related lower temperature (220-250°C) neutral to slightly acidic fluids have been ascribed by others to return circulation of hydrothermal fluids that did not interact with seawater. Because altered samples have a higher Pb content than the fresh precursor, leaching of fresh volcanic rocks cannot be the source of Pb in the hydrothermal systems.

Lithogeochemical and isotope (Rb-Sr, Sm-Nd, Pb-Pb whole rock and Lu-Hf zircon) composition of samples from the Shackleton Range, East Antarctica

Three distinct, spatially separated crustal terranes have been recognised in the Shackleton Range, East Antarctica: the Southern, Eastern and Northern Terranes. Mafic gneisses from the Southern Terrane provide geochemical evidence for a within-plate, probably back-arc origin of their protoliths. A plume-distal ridge origin in an incipient ocean basin is the favoured interpretation for the emplacement site of these rocks at c. 1850 Ma, which, together with a few ocean island basalts, were subsequently incorporated into an accretionary continental arc/supra-subduction zone tectonic setting. Magmatic underplating resulted in partial melting of the lower crust, which caused high-temperature granulite-facies metamorphism in the Southern Terrane at c. 1710-1680 Ma. Mafic and felsic gneisses there are characterised by isotopically depleted, positive Nd and Hf initials and model ages between 2100 and 2000 Ma. They may be explained as juvenile additions to the crust towards the end of the Palaeoproterozoic. These juvenile rocks occur in a narrow, c. 150 km long E-W trending belt, inferred to trace a suture that is associated with a large Palaeoproterozoic accretionary orogenic system. The Southern Terrane contains many features that are similar to the Australo-Antarctic Mawson Continent and may be its furthermost extension into East Antarctica. The Eastern Terrane is characterised by metagranitoids that formed in a continental volcanic arc setting during a late Mesoproterozoic orogeny at c. 1060 Ma. Subsequently, the rocks experienced high-temperature metamorphism during Pan-African collisional tectonics at 600 Ma. Isotopically depleted zircon grains yielded Hf model ages of 1600-1400 Ma, which are identical to Nd model ages obtained from juvenile metagranitoids. Most likely, these rocks trace the suture related to the amalgamation of the Indo-Antarctic and West Gondwana continental blocks at ~600 Ma. The Eastern Terrane is interpreted as the southernmost extension of the Pan-African Mozambique/Maud Belt in East Antarctica and, based on Hf isotope data, may also represent a link to the Ellsworth-Whitmore Mountains block in West Antarctica and the Namaqua-Natal Province of southern Africa. Geochemical evidence indicates that the majority of the protoliths of the mafic gneisses in the Northern Terrane formed as oceanic island basalts in a within-plate setting. Subsequently the rocks were incorporated into a subduction zone environment and, finally, accreted to a continental margin during Pan-African collisional tectonics. Felsic gneisses there provide evidence for a within-plate and volcanic arc/collisional origin. Emplacement of granitoids occurred at c. 530 Ma and high-temperature, high-pressure metamorphism took place at 510-500 Ma. Enriched Hf and Nd initials and Palaeoproterozoic model ages for most samples indicate that no juvenile material was added to the crust of the Northern Terrane during the Pan-African Orogeny but recycling of older crust or mixing of crustal components of different age must have occurred. Isotopically depleted mafic gneisses, which are spatially associated with eclogite-facies pyroxenites, yielded late Mesoproterozoic Nd model ages. These rocks occur in a narrow, at least 100 km long, E-W trending belt that separates alkaline ocean island metabasalts and within-plate metagranitoids from volcanic arc metabasalts and volcanic arc/syn-collisional metagranitoids in the Northern Terrane. This belt is interpreted to trace the late Neoproterozoic/early Cambrian Pan-African collisional suture between the Australo-Antarctic and the combined Indo-Antarctic/West Gondwana continental blocks that formed during the final amalgamation of Gondwana.

Geochemistry at DSDP Leg 59

Deep Sea Drilling Project Legs 59 and 60 drilled 15 sites along an east-west transect at 18°N from the West Philippine Basin to the Mariana Trench (Fig. 1) in order to study the nature and genesis of the back-arc, marginal basins and the remnant and active arcs of the region. Leg 59 drilled at five sites at the western end of the traverse: Site 447 in the West Philippine Basin; Site 448 on the Palau-Kyushu Ridge; Sites 449 and 450 in the Parece Vela Basin; and Site 451 on the West Mariana Ridge. Penetration into basaltic basement of these sites was 183.5 meters at 447 (8 basalt flows); 623 meters at 448 (46 basalt flows, sills, and dikes and volcaniclastic units); 40.5 meters at 449 (2 basalt flows); 7 meters at 450 (1 basalt intrusion); and 4 meters of basalt breccia at 451 overlain by 861 meters of volcaniclastic sedimentary rocks.

Rhenium and platinum-group element concentrations and Os isotope data for Philippine Sea Plate basalts from DSDP Legs 31, 58, and 59

Platinum-group elements (PGE), rhenium and osmium isotope data are reported for basalts from Deep Sea Drilling Project cores in the Philippine Sea Plate (PSP). Lithophile trace element and isotopic characteristics indicate a range of source components including DMM, EMII and subduction-enriched mantle. MORB-like basalts possess smooth, inclined chondrite-normalised PGE patterns with high palladium-PGE/iridium-PGE ratios, consistent with previously published data for MORB, and with the inferred compatibility of PGE. In contrast, while basalts with EMII-type lithophile element chemistry possess high Pt/Ir ratios, many have much lower Pd/Ir and unusually high Ru/Ir of >10. Similarly, back-arc samples from the Shikoku and Parece-Vela basins have very high Ru/Ir ratios (>30) and Pd/Ir as low as 1.1. Such extreme Pd/Ir and Ru/Ir ratios have not been previously reported in mafic volcanic suites and cannot be easily explained by variable degrees of melting, fractional crystallisation or by a shallow-level process such as alteration or degassing. The data appear most consistent with sampling of at least two mantle components with distinct PGE compositions. Peridotites with the required PGE characteristics (i.e. low Pd, but relatively high Ru and Re) have not been documented in oceanic mantle, but have been found in sub-continental mantle lithosphere and are the result of considerable melt depletion and selective metasomatic enrichment (mainly Re). The long-term presence of subduction zones surrounding the Philippine Sea Plate makes this a prime location for metasomatic enrichment of mantle, either through fluid enrichment or infiltration by small melt fractions. The Re-Os isotope data are difficult to interpret with confidence due to low Os concentrations in most samples and the uncertainty in sample age. Data for Site 444A (Shikoku Basin) give an age of 17.7+/-1.3 Ma (MSWD = 14), consistent with the proposed age of basement at the site and thus provides the first robust radiometric age for these samples. The initial 187Os/188Os of 0.1298+/-0.0069 is consistent with global MORB, and precludes significant metasomatic enrichment of Os by radiogenic slab fluids. Re-Os data for Sites 446A (two suites, Daito Basin) and 450 (Parece-Vela Basin) indicate ages of 73, 68 and 43 Ma, which are respectively, 30, 17 and >12 Ma older than previously proposed ages. The alkalic and tholeiitic suites from Site 446A define regression lines with different 187Os/188Osinitial (0.170+/-0.033 and 0.112+/-0.024, respectively) which could perhaps be explained by preferential sampling of interstitial, metasomatic sulphides (with higher time-integrated Re/Os ratios) by smaller percentage alkalic melts. One sample, with lithophile elements indistinguishable from MORB, is Os-rich (146 pg/g) and has an initial 187Os/188Os of 0.1594, which is at the upper limit of the accepted OIB range. Given the Os-rich nature of this sample and the lack of evidence for subduction or recycled crust inputs, this osmium isotope ratio likely reflects heterogeneity in the DMM. The dataset as a whole is a striking indication of the possible PGE and Os isotope variability within a region of mantle that has experienced a complex tectonic history.

Geochemistry of Central American Volcanic Arc basement samples

The Central American Volcanic Arc (CAVA) has been the subject of intensive research over the past few years, leading to a variety of distinct models for the origin of CAVA lavas with various source components. We present a new model for the NW Central American Volcanic Arc based on a comprehensive new geochemical data set (major and trace element and Sr-Nd-Pb-Hf-O isotope ratios) of mafic volcanic front (VF), behind the volcanic front (BVF) and back-arc (BA) lava and tephra samples from NW Nicaragua, Honduras, El Salvador and Guatemala. Additionally we present data on subducting Cocos Plate sediments (from DSDP Leg 67 Sites 495 and 499) and igneous oceanic crust (from DSDP Leg 67 Site 495), and Guatemalan (Chortis Block) granitic and metamorphic continental basement. We observe systematic variations in trace element and isotopic compositions both along and across the arc. The data require at least three different endmembers for the volcanism in NW Central America. (1) The NW Nicaragua VF lavas require an endmember with very high Ba/(La, Th) and U/Th, relatively radiogenic Sr, Nd and Hf but unradiogenic Pb and low d18O, reflecting a largely serpentinite-derived fluid/hydrous melt flux from the subducting slab into a depleted N-MORB type of mantle wedge. (2) The Guatemala VF and BVF mafic lavas require an enriched endmember with low Ba/(La, Th), U/Th, high d18O and radiogenic Sr and Pb but unradiogenic Nd and Hf isotope ratios. Correlations of Hf with both Nd and Pb isotopic compositions are not consistent with this endmember being subducted sediments. Granitic samples from the Chiquimula Plutonic Complex in Guatemala have the appropriate isotopic composition to serve as this endmember, but the large amounts of assimilation required to explain the isotope data are not consistent with the basaltic compositions of the volcanic rocks. In addition, mixing regressions on Nd vs. Hf and the Sr and O isotope plots do not go through the data. Therefore, we propose that this endmember could represent pyroxenites in the lithosphere (mantle and possibly lower crust), derived from parental magmas for the plutonic rocks. (3) The Honduras and Caribbean BA lavas define an isotopically depleted endmember (with unradiogenic Sr but radiogenic Nd, Hf and Pb isotope ratios), having OIB-like major and trace element compositions (e.g. low Ba/(La, Th) and U/Th, high La/Yb). This endmember is possibly derived from melting of young, recycled oceanic crust in the asthenosphere upwelling in the back-arc. Mixing between these three endmember types of magmas can explain the observed systematic geochemical variations along and across the NW Central American Arc.

(Tables 1-3) Pb and Sr isotopic data for sediments and rocks from DSDP holes in the Mariana island-arc system

Analyses of the isotopic composition of Pb in (1) western Pacific Ocean sediments [Jurassic(?) to Pleistocene in age, including clays and biogenic oozes], (2) Pacific Ocean basaltic rocks, (3) Mariana frontal arc volcanic rocks (Eocene to Miocene), and (4) Mariana active arc volcanic rocks [Pliocene (?) to Holocene] indicate that Pacific Ocean sediments could not have been a significant component of the source material for the Mariana arc volcanic rocks. Calculations involving the average concentrations and isotopic compositions of Pb in oceanic sediments, sea-floor basaltic rocks, and the Mariana arc volcanic rocks suggest that the sediment component must have been less than 1 percent of this source material. The Pb isotopic compositions of the Mariana arc volcanic rocks lie, within experimental error, along the trend of available Pacific Ocean basalt analyses in versus 207Pb/204Pb versus 206Pb/204Pb and 208Pb/204Pb versus 206Pb/204Pb diagrams. Isotopic analyses of Pb in Pacific Ocean sediments do not lie along this trend; they have higher 207Pb/204Pb and 208Pb/204Pb values for comparable 206Pb/204Pb ratios. Clayey sediments generally have higher 208Pb/204Pb and 207Pb/204Pb ratios than biogenic oozes regardless of the age of the sediment. Comparison of combined Sr and Pb isotopic analyses for (1) mantle-derived materials erupted through oceanic crust, (2) altered ocean-floor basaltic rocks, and (3) volcanic rocks from oceanic island arcs suggests that the Mariana arc volcanic rocks were derived, at least in part, from altered Pacific lithosphere subducted beneath the Mariana arc. Unaltered basalts from the Mariana inter-arc basin (Mariana Trough) have Pb and Sr isotopic compositions that are very similar to those reported for some Hawaiian volcanic rocks but distinct from Mariana active and frontal arc compositions. These observations, in addition to existing major-and trace-element data, support a mantle origin for the interarc basin volcanic rocks. Dacites dredged from the Mariana remnant arc (South Honshu Ridge) have Pb isotopic compositions that are within experimental error of the active-arc analyses, consistent with a genetic relation.

Pleistocene rhyolithic tephra of ODP Leg 181 sites

Taupo Volcanic Zone (TVZ), in the North Island, New Zealand, is arguably the most active Quaternary rhyolitic system in the world. Numerous and widespread rhyolitic tephra layers, sourced from the TVZ, form valuable chronostratigraphic markers in onshore and offshore sedimentary sequences. In deep-sea cores from Ocean Drilling Program (ODP) Leg 181 Sites 1125, 1124, 1123 and 1122, located east of New Zealand, ca 100 tephra beds are recognised post-dating the Plio-Pleistocene boundary at 1.81 Ma. These tephras have been dated by a combination of magnetostratigraphy, orbitally tuned stable-isotope data and isothermal plateau fission track ages. The widespread occurrence of ash offshore to the east of New Zealand is favoured by the small size of New Zealand, the explosivity of the mainly plinian and ignimbritic eruptions and the prevailing westerly wind field. Although some tephras can be directly attributed to known TVZ eruptions, there are many more tephras represented within ODP-cores that have yet to be recognised in near-source on-land sequences. This is due to proximal source area erosion and/or deep burial as well as the adverse effect of vapour phase alteration and devitrification within near-source welded ignimbrites. Despite these difficulties, a number of key deep-sea tephras can be reliably correlated to equivalent-aged tephra exposed in uplifted marine back-arc successions of Wanganui Basin where an excellent chronology has been developed based on magnetostratigraphy, orbitally calibrated sedimentary cycles and isothermal plateau fission track ages on tephra. Significant Pleistocene tephra markers include: the Kawakawa, Omataroa, Rangitawa/Onepuhi, Kaukatea, Kidnappers-B, Potaka, Unit D/Ahuroa, Ongatiti, Rewa, Sub-Rewa, Pakihikura, Ototoka and Table Flat Tephras. Six other tephra layers are correlated between ODP-core sites but have yet to be recognised within onshore records. The identification of Pleistocene TVZ-sourced tephras within the ODP record, and their correlation to Wanganui Basin and other onshore sites is a significant advance as it provides: (1) an even more detailed history of the TVZ than can be currently achieved from the near-source record, (2) a high-resolution tephrochronologic framework for future onshore-offshore paleoenvironmental reconstructions, and (3) well-dated tephra beds correlated from the offshore ODP sites with astronomically tuned timescales provide an opportunity to critically evaluate the chronostratigraphic framework for onshore Plio-Pleistocene sedimentary sequences (e.g. Wanganui Basin, cf. Naish et al. (1998, doi:10.1016/S0277-3791(97)00075-9).

Geochemistry and radiogenic isotope ratios of the volcanic sequences from Izu-Bonin Arc, Japan

New radiogenic isotope and trace element data are presented for the volcanic sequences along 600 km of the active Izu-Bonin arc, the Oligocene Izu arc, and their associated rift basins. As with many intra-oceanic island arcs, the Pliocene-Recent Izu-Bonin frontal-arc lavas are highly depleted in Zr, Nb and the rare-earth elements relative to typical mid-ocean ridge basalt (MORB), indicating that the mantle wedge source has undergone a previous episode of melting. Ratios between these elements (such as Nb/Zr and La/Sm), as well as 143Nd/144Nd, do not vary significantly along the length of the frontal-arc. These parameters suggest that each of the arc volcanoes is derived from similar melt fractions of the mantle wedge. However, Ba/Zr, Ba/Rb and 87Sr/86Sr increase along the frontal-arc to the north. This leads us to propose that a variable enrichment in Ba and radiogenic Sr is superimposed on the mantle wedge. Sr-Nd and Pb-Nd isotope variation indicate that both Sr and Pb become more radiogenic after fluid addition. However, Pb isotope ratios do not correlate with increases in Pb concentration or ratios such as Ba/Zr and Nb/Pb. In other words, the Pb isotopic composition of the arc lavas appears to be independent of the amount of Pb introduced by subduction fluids into the mantle source. This buffering of Pb isotopes along the frontal-arc means that the isotopic composition of the lavas is indistinguishable from that of the fluid. Isotopic mixing models presented for the arc are only illustrative of the many plausible combinations of components and quantities. Despite this, we are able to determine that the mantle wedge has isotopic characteristics similar to Indian Ocean MORB, and that the subduction-fluid solute is primarily derived from subducted oceanic basalt with a <2% contribution from subducted sediment. Lavas in the Oligocene Izu arc and fore-arc basin were derived from a mantle wedge of similar composition to the active arc. Despite levels of Pb enrichment comparable to those of the modern arc, the Pb isotopes of the Oligocene volcanics indicate a lower sediment input into the melting region.

Age and nature of 560–520 Ma calc-alkaline granitoids of Biarjmand, northeast Iran: insights into Cadomian arc magmatism in northern Gondwana

The Biarjmand granitoids and granitic gneisses in northeast Iran are part of the Torud–Biarjmand metamorphic complex, where previous zircon U–Pb geochronology show ages of ca. 554–530 Ma for orthogneissic rocks. Our new U–Pb zircon ages confirm a Cadomian age and show that the granitic gneiss is ~30 million years older (561.3 ± 4.7 Ma) than intruding granitoids(522.3 ± 4.2 Ma; 537.7 ± 4.7 Ma). Cadomian magmatism in Iran was part of an approximately 100-million-year-long episode of subduction-related arc and back-arc magmatism, which dominated the whole northern Gondwana margin, from Iberia to Turkey and Iran. Major REE and trace element data show that these granitoids have calc-alkaline signatures. Their zircon O (δ18O = 6.2–8.9‰) and Hf (–7.9 to +5.5; one point with εHf ~ –17.4) as well as bulk rock Nd isotopes (εNd(t)= –3 to –6.2) show that these magmas were generated via mixing of juvenile magmas with an older crust and/or melting of middle continental crust. Whole-rock Nd and zircon Hf model ages (1.3–1.6 Ga) suggest that this older continental crust was likely to have been Mesoproterozoic or even older. Our results, including variable zircon εHf(t) values, inheritance of old zircons and lack of evidence for juvenile Cadomian igneous rocks anywhere in Iran, suggest that the geotectonic setting during late Ediacaran and early Cambrian time was a continental magmatic arc rather than back-arc for the evolution of northeast Iran Cadomian igneous rocks.

Two possible hydrothermal vents in the northern Okinawa Trough

As the Okinawa Trough is a back-are basin in early spreading, modern submarine hydrothermal activity and minerallization have many characteristics which have aroused wide attention. Up to now three well-known hydrothermal venting areas are all located in the middle part of the trough, During two cruise investigations to map and sample the seafloor numbers of Calyptogena sp, shells were dredged at two sites in the northern trough with comparatively thicker crust and numerous submarine volcanoes. Based on the fact that Calyptogena sp, is only observed around the hydrothermal vents and lives on hydrothermal activities, it is predicted that there is the possibility of modern hydrothermal activities in the northern part of the trough. In this note, the shell is carefully characterized and the sample locations with possible hydrothermal activity are given. it Is pointed out that the research of biogenic fossils to trace hydrothermal activity changes in venting time, strength fluctuations, evolution In chemical compositions and so on should be stressed in the future in addition to the study of the ecological characteristics of hydrothermal organisms.

芝麻坊石榴子石二辉橄榄岩多期变质演化及其地幔过程

Many garnet peridotite bodies are enclosed in ultrahigh-pressure (UHP) gneisses and/or migmatites in worldwide UHP terranes formed by subduction of continental crust. On the basis of petrochemical data, a group of garnet peridotites have been derived from depleted mantle and were subsequently metasomatized by melts and/or fluids derived from the subducted continental crust. However, their depletion and enrichment processes and tectonic evolutions are still in conflicts. New evidences for metamorphism of garnet lherzolite from Zhimafang, Donghai County, Sulu UHP terrane are reported. The garnet lherzolite have experienced a prolonged multistage metamorphic history. At least seven stages of recrystallization have been identified based on detailed analysis of reaction textures and mineral compositions. Stage I was a high-pressure and high-temperature enriched garnet lherzolite stage, which is inferred from the presence of high Ca-Cr core of garnet porphyroclast and inclusions of high-Mg clinopyroxene, high-Al-Cr orthopyroxene and high-Mg olivine. Stage II is a high-temperature and low-pressure depleted spinel-hurzbergite or spinel-dunite stage, as indicated by the presence of relict Al-rich spinel, very high-Mg and low-Ni olivine and high-Mg orthopyroxene included in the low-Cr mantle of the porphyroclastic garnet and core of fine-grained neoblastic garnet, clinopyroxene is absent in this stage. Stage III is an hydrous amphibole spinel-lherzolite stage, which recorded events of cooling and metasomatic re-enrichment, this stage is manifested by metasomatic origin of amphibole and phlogopite-bearing porphyroblastic clinopyroxene, and porphyroblastic orthopyroxene. Stage IV is a high-pressure amphibole garnet-lherzolite stage, which is indicated by the formation of low-Cr mantle of the porphyroclastic garnet and amphibole-bearing low-Cr core of neoblastic garnet. Stage V is an UHP metamorphic garnet-lherzolite stage, which is characterized by the formation of high-Cr rim of both porphyroclastic and neoblastic garnet and recrystallization of olivine, clinopyroxene and orthopyroxene in the matrix. During UHP metamorphism, the garnet lherzolite is dehydrated, hornblende decomposed to clinopyroxene and olivine. Stage VI is a high-pressure decompression amphibole garnet-lherzolite stage, indicated by formation of later coarse-grained pargasitic hornblende and phlogopite in the garnet stability field. Stage VII is a low-pressure decompression amphibole-chlorite spinel-lherzolite stage, indicated by replacement of garnet by kelyphite of high-Al orthopyroxene + aluminous spinel + tremolitic amphibole + chlorite + talc. The metamorphic evolutions of Zhimafang garnet lherzolite suggest that it displays progressive mantle wedge convection during the subduction of previous oceanic and subsequent continental slab. We propose that the Zhimafang garnet lherzolite were originated from enriched deep mantle wedge above the previously subducted oceanic slab, subduction of oceanic slab resulted in their convection to shallower back arc and sub-arc setting, decompressional melting transformed the enriched garnet-lherzolite to depleted spinel-hurzbergite or spinel-dunite, the spinel-hurzbergite or spinel dunite was then convected to the hydrous mantle wedge corner driven by corner flow and was cooled and metasomatized by slab-derived melts/fluids, and was transformed to enriched lherzolite. The lherzolites formed a downward mantle wedge layer above successively subducted continental crust. The peridotite subducted together with the underlying continental crust and suffered UHP metamorphism. Finally, the garnet-lherzolite exhumed to the earth surface together with the UHP terrane. Detailed analyses of reaction textures and mineral compositions revealed several stages of metasomatism related to continental subduction and exhumation.