Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

Gravity anomalies and associated tectonic features over the Indian Peninsular Shield and adjoining ocean basins

A combined gravity map over the Indian Peninsular Shield (IPS) and adjoining oceans brings out well the inter-relationships between the older tectonic features of the continent and the adjoining younger oceanic features. The NW–SE, NE–SW and N–S Precambrian trends of the IPS are reflected in the structural trends of the Arabian Sea and the Bay of Bengal suggesting their probable reactivation. The Simple Bouguer anomaly map shows consistent increase in gravity value from the continent to the deep ocean basins, which is attributed to isostatic compensation due to variations in the crustal thickness. A crustal density model computed along a profile across this region suggests a thick crust of 35–40 km under the continent, which reduces to 22/20–24 km under the Bay of Bengal with thick sediments of 8–10 km underlain by crustal layers of density 2720 and 2900/2840 kg/m³. Large crustal thickness and trends of the gravity anomalies may suggest a transitional crust in the Bay of Bengal up to 150–200 km from the east coast. The crustal thickness under the Laxmi ridge and east of it in the Arabian Sea is 20 and 14 km, respectively, with 5–6 km thick Tertiary and Mesozoic sediments separated by a thin layer of Deccan Trap. Crustal layers of densities 2750 and 2950 kg/m³ underlie sediments. The crustal density model in this part of the Arabian Sea (east of Laxmi ridge) and the structural trends similar to the Indian Peninsular Shield suggest a continent–ocean transitional crust (COTC). The COTC may represent down dropped and submerged parts of the Indian crust evolved at the time of break-up along the west coast of India and passage of Reunion hotspot over India during late Cretaceous. The crustal model under this part also shows an underplated lower crust and a low density upper mantle, extending over the continent across the west coast of India, which appears to be related to the Deccan volcanism. The crustal thickness under the western Arabian Sea (west of the Laxmi ridge) reduces to 8–9 km with crustal layers of densities 2650 and 2870 kg/m³ representing an oceanic crust.     

A possible Caledonide arm through the Barents Sea imaged by OBS data

The assembly of the crystalline basement of the western Barents Sea is related to the Caledonian orogeny during the Silurian. However, the development southeast of Svalbard is not well understood, as conventional seismic reflection data does not provide reliable mapping below the Permian sequence. A wide-angle seismic survey from 1998, conducted with ocean bottom seismometers in the northwestern Barents Sea, provides data that enables the identification and mapping of the depths to crystalline basement and Moho by ray tracing and inversion. The four profiles modeled show pre-Permian basins and highs with a configuration distinct from later Mesozoic structural elements. Several strong reflections from within the crystalline crust indicate an inhomogeneous basement terrain. Refractions from the top of the basement together with reflections from the Moho constrain the basement velocity to increase from 6.3 km s⁻¹ at the top to 6.6 km s⁻¹ at the base of the crust. On two profiles, the Moho deepens locally into root structures, which are associated with high top mantle velocities of 8.5 km s⁻¹. Combined P- and S-wave data indicate a mixed sand/clay/carbonate lithology for the sedimentary section, and a predominantly felsic to intermediate crystalline crust. In general, the top basement and Moho surfaces exhibit poor correlation with the observed gravity field, and the gravity models required high-density bodies in the basement and upper mantle to account for the positive gravity anomalies in the area. Comparisons with the Ural suture zone suggest that the Barents Sea data may be interpreted in terms of a proto-Caledonian subduction zone dipping to the southeast, with a crustal root representing remnant of the continental collision, and high mantle velocities and densities representing eclogitized oceanic crust. High-density bodies within the crystalline crust may be accreted island arc or oceanic terrain. The mapped trend of the suture resembles a previously published model of the Caledonian orogeny. This model postulates a separate branch extending into central parts of the Barents Sea coupled with the northerly trending Svalbard Caledonides, and a microcontinent consisting of Svalbard and northern parts of the Barents Sea independent of Laurentia and Baltica at the time. Later, compressional faulting within the suture zone apparently formed the Sentralbanken High.     

Poisson's ratio in the oceanic crust — a review

Laboratory samples from the upper oceanic crust (tholeiitic basalt flows) that have not been significantly weathered, hydrothermally altered or fractured have a typical Poisson's ratio of 0.30 ( ) and a compressional velocity of 6.0 km s⁻¹; from the middle crust (dolerite sheeted dykes) a ratio of 0.28 ( ) and a velocity of 6.7 km s⁻¹; from the lower crust (gabbro) a ratio of 0.31 ( ) and a velocity of 7.1 km s⁻¹; and from the uppermost mantle a ratio of 0.24 ( ) and a velocity of 8.4 km s⁻¹. These sample values are representative of the large scale insitu values for the middle and lower crust and for the upper mantle. The upper crust is modified by several processes that decrease the velocity and generally increase Poisson's ratio: (1) the formation of an irregular layer of low temperature weathering generally less than 50 m thick; (2) large scale porosity in the form of drained pillows and lava tubes, of talus and rubble and of large open fractures; (3) where there was a high sedimentation rate over the ridge that formed the crust, hydrothermal alteration and intercalation of basalt and sediments. The Poisson's ratios of both high velocity sediments and of crystalline continental crustal rocks generally are significantly lower than the ratios of oceanic crustal rocks of similar compressional wave velocity. Thus, the use of shear wave velocities should permit the separation of these different formations which frequently cannot be distinguished on the basis of compressional wave seismic refraction data alone.     

Two-dimensional numerical modeling of tectonic and metamorphic histories at active continental margins

The evolution of an active continental margin is simulated in two dimensions, using a finite difference thermomechanical code with half-staggered grid and marker-in-cell technique. The effect of mechanical properties, changing as a function of P and T, assigned to different crustal layers and mantle materials in the simple starting structure is discussed for a set of numerical models. For each model, representative P–T paths are displayed for selected markers. Both the intensity of subduction erosion and the size of the frontal accretionary wedge are strongly dependent on the rheology chosen for the overriding continental crust. Tectonically eroded upper and lower continental crust is carried down to form a broad orogenic wedge, intermingling with detached oceanic crust and sediments from the subducted plate and hydrated mantle material from the overriding plate. A small portion of the continental crust and trench sediments is carried further down into a narrow subduction channel, intermingling with oceanic crust and hydrated mantle material, and to some extent extruded to the rear of the orogenic wedge underplating the overriding continental crust. The exhumation rates for (ultra)high pressure rocks can exceed subduction and burial rates by a factor of 1.5–3, when forced return flow in the hanging wall portion of the self-organizing subduction channel is focused. The simulations suggest that a minimum rate of subduction is required for the formation of a subduction channel, because buoyancy forces may outweigh drag forces for slow subduction. For a weak upper continental crust, simulated by a high pore pressure coefficient in the brittle regime, the orogenic wedge and megascale melange reach a mid- to upper-crustal position within 10–20 Myr (after 400–600 km of subduction). For a strong upper crust, a continental lid persists over the entire time span covered by the simulation. The structural pattern is similar in all cases, with four zones from trench toward arc: (a) an accretionary complex of low-grade metamorphic sedimentary material; (b) a wedge of mainly continental crust, with medium-grade HP metamorphic overprint, wound up and stretched in a marble cake fashion to appear as nappes with alternating upper and lower crustal provenance, and minor oceanic or hydrated mantle interleaved material; (c) a megascale melange composed of high-pressure and ultrahigh-pressure metamorphic oceanic and continental crust, and hydrated mantle, all extruded from the subduction channel; (d) zone represents the upward tilted frontal part of the remaining upper plate lid in the case of a weak upper crust. The shape of the P–T paths and the time scales correspond to those typically recorded in orogenic belts. Comparison of the numerical results with the European Alps reveals some similarities in their gross structural and metamorphic pattern exposed after collision. A similar structure may be developed at depth beneath the forearc of the Andes, where the importance of subduction erosion is well documented, and where a strong upper crust forms a stable lid.     

Seismicity of the Hellenic subduction zone in the area of western and central Crete observed by temporary local seismic networks

Temporary local seismic networks were installed in western Crete, in central Crete, and on the island Gavdos south of western Crete, respectively, in order to image shallow seismically active zones of the Hellenic subduction zone.More than 4000 events in the magnitude range between −0.5 and 4.8 were detected and localized. The resulting three-dimensional hypocenter distribution allows the localization of seismically active zones in the area of western and central Crete from the Mediterranean Ridge to the Cretan Sea. Furthermore, a three-dimensional structural model of the studied region was compiled based on results of wide-angle seismics, surface wave analysis and receiver function studies. The comparison of the hypocenter distribution and the structure has allowed intraplate and interplate seismicity to be distinguished.High interplate seismicity along the interface between the subducting African lithosphere and the Aegean lithosphere was found south of western Crete where the interface is located at about 20 to 40 km depth. An offset between the southern border of the Aegean lithosphere and the southern border of active interplate seismicity is observed. In the area of Crete, the offset varies laterally along the Hellenic arc between about 50 and 70 km.A southwards dipping zone of high seismicity within the Aegean lithosphere is found south of central Crete in the region of the Ptolemy trench. It reaches from the interface between the plates at about 30 km depth towards the surface. In comparison, the Aegean lithosphere south of western Crete is seismically much less active including the region of the Ionian trench. Intraplate seismicity within the Aegean plate beneath Crete and north of Crete is confined to the upper about 20 km. Between 20 and 40 km depth beneath Crete, the Aegean lithosphere appears to be seismically inactive. In western Crete, the southern and western borders of this aseismic zone correlate strongly with the coastline of Crete.     

The Levantine Basin—crustal structure and origin

The origin of the Levantine Basin in the Southeastern Mediterranean Sea is related to the opening of the Neo-Tethys. The nature of its crust has been debated for decades. Therefore, we conducted a geophysical experiment in the Levantine Basin. We recorded two refraction seismic lines with 19 and 20 ocean bottom hydrophones, respectively, and developed velocity models. Additional seismic reflection data yield structural information about the upper layers in the first few kilometers. The crystalline basement in the Levantine Basin consists of two layers with a P-wave velocity of 6.0–6.4 km/s in the upper and 6.5–6.9 km/s in the lower crust. Towards the center of the basin, the Moho depth decreases from 27 to 22 km. Local variations of the velocity gradient can be attributed to previously postulated shear zones like the Pelusium Line, the Damietta–Latakia Line and the Baltim–Hecateus Line. Both layers of the crystalline crust are continuous and no indication for a transition from continental to oceanic crust is observed. These results are confirmed by gravity data. Comparison with other seismic refraction studies in prolongation of our profiles under Israel and Jordan and in the Mediterranean Sea near Greece and Sardinia reveal similarities between the crust in the Levantine Basin and thinned continental crust, which is found in that region. The presence of thinned continental crust under the Levantine Basin is therefore suggested. A β-factor of 2.3–3 is estimated. Based on these findings, we conclude that sea-floor spreading in the Eastern Mediterranean Sea only occurred north of the Eratosthenes Seamount, and the oceanic crust was later subducted at the Cyprus Arc.     

Crustal Structure across the Northwestern Margin of South China Sea: Evidence for Magma‐poor Rifting from a Wide‐angle Seismic Profile

We present results from a 484 km wide-angle seismic profile acquired in the northwest part of the South China Sea (SCS) during OBS2006 cruise. The line that runs along a previously acquired multi-channel seismic line (SO49-18) crosses the continental slope of the northern margin, the Northwest Subbasin (NWSB) of the South China Sea, the Zhongsha Massif and partly the oceanic basin of the South China Sea. Seismic sections recorded on 13 ocean-bottom seismometers were used to identify refracted phases from the crustal layer and also reflected phases from the crust-mantle boundary (Moho). Inversion of the traveltimes using a simple start model reveals crustal images in the study area. The velocity model shows that crustal thickness below the continental slope is between 14 and 23 km. The continental part of the line is characterized by gentle landward mantle uplift and an abrupt oceanward one. The velocities in the lower crust do not exceed 6.9 km/s. With the new data we can exclude a high-velocity lower crustal body (velocities above 7.0 km/s) at the location of the line. We conclude that this part of the South China Sea margin developed by a magma-poor rifting. Both, the NWSB and the Southwest Sub-basin (SWSB) reveal velocities typical for oceanic crust with crustal thickness between 5 and 7 km. The Zhongsha Massif in between is extremely stretched with only 6–10 km continental crust left. Crustal velocity is below 6.5 km/s; possibly indicating the absence of the lower crust. Multi-channel seismic profile shows that the Yitongansha Uplift in the slope area and the Zhongsha Massif are only mildly deformed. We considered them as rigid continent blocks which acted as rift shoulders of the main rift subsequently resulting in the formation of the Northwest Sub-basin. The extension was mainly accommodated by a ductile lower crustal flows, which might have been extremely attenuated and flow into the oceanic basin during the spreading stage. We compared the crustal structures along the northern margin and found an east-west thicken trend of the crust below the continent slope. This might be contributed by the east-west sea-floor spreading along the continental margin.     

Crustal structure and transform margin development south of Svalbard based on ocean bottom seismometer data

The Barents Sea is located in the northwestern corner of the Eurasian continent, where the crustal terrain was assembled in the Caledonian orogeny during Late Ordovician and Silurian times. The western Barents Sea margin developed primarily as a transform margin during the early Tertiary. In the northwestern part south of Svalbard, multichannel reflection seismic lines have poor resolution below the Permian sequence, and the early post-orogenic development is not well known here. In 1998, an ocean bottom seismometer (OBS) survey was collected southwest to southeast of the Svalbard archipelago. One profile was shot across the continental transform margin south of Svalbard, which is presented here. P-wave modeling of the OBS profile indicates a Caledonian suture in the continental basement south of Svalbard, also proposed previously based on a deep seismic reflection line coincident with the OBS profile. The suture zone is associated with a small crustal root and westward dipping mantle reflectivity, and it marks a boundary between two different crystalline basement terrains. The western terrain has low (6.2–6.45 km s⁻¹) P-wave velocities, while the eastern has higher (6.3–6.9 km s⁻¹) velocities. Gravity modeling agrees with this, as an increased density is needed in the eastern block. The S-wave data predict a quartz-rich lithology compatible with felsic gneiss to granite within and west of the suture zone, and an intermediate lithological composition to the east. A geological model assuming westward dipping Caledonian subduction and collision can explain the missing lower crust in the western block by subduction erosion of the lower crust, as well as the observed structuring. Due to the transform margin setting, the tectonic thinning of the continental block during opening of the Norwegian-Greenland Sea is restricted to the outer 35 km of the continental block, and the continent–ocean boundary (COB) can be located to within 5 km in our data. Distinct from the outer high commonly observed on transform margins, the upper part of the continental crust at the margin is dominated by two large, rotated down-faulted blocks with throws of 2–3 km on each fault, apparently formed during the transform margin development. Analysis of the gravity field shows that these faults probably merge to one single fault to the south of our profile, and that the downfaulting dominates the whole margin segment from Spitsbergen to Bjørnøya. South of Bjørnøya, the faulting leaves the continental margin to terminate as a graben 75 km south of the island. Adjacent to the continental margin, there is no clear oceanic layer 2 seismic signature. However, the top basement velocity of 6.55 km s⁻¹ is significantly lower than the high (7 km s⁻¹) velocity reported earlier from expanding spread profiles (ESPs), and we interpret the velocity structure of the oceanic crust to be a result of a development induced by the 7–8-km-thick sedimentary overburden.     

Global 1° × 1° thermal model TC1 for the continental lithosphere: Implications for lithosphere secular evolution

This paper reports a new 1° × 1° global thermal model for the continental lithosphere (TC1). Geotherms for continental terranes of different ages (> 3.6 Ga to present) constrained by reliable data on borehole heat flow measurements (Artemieva, I.M., Mooney, W.D. 2001. Thermal structure and evolution of Precambrian lithosphere: a global study. J. Geophys. Res 106, 16387–16414.), are statistically analyzed as a function of age and are used to estimate lithospheric temperatures in continental regions with no or low-quality heat flow data (ca. 60% of the continents). These data are supplemented by cratonic geotherms based on electromagnetic and xenolith data; the latter indicate the existence of Archean cratons with two characteristic thicknesses, ca. 200 and > 250 km. A map of tectono-thermal ages of lithospheric terranes complied for the continents on a 1° × 1° grid and combined with the statistical age relationship of continental geotherms (z = 0.04  t + 93.6, where z is lithospheric thermal thickness in km and t is age in Ma) formed the basis for a new global thermal model of the continental lithosphere (TC1). The TC1 model is presented by a set of maps, which show significant thermal heterogeneity within continental upper mantle, with the strongest lateral temperature variations (as large as 800 °C) in the shallow mantle. A map of the depth to a 550 °C isotherm (Curie isotherm for magnetite) in continental upper mantle is presented as a proxy to the thickness of the magnetic crust; the same map provides a rough estimate of elastic thickness of old (> 200 Ma) continental lithosphere, in which flexural rigidity is dominated by olivine rheology of the mantle.Statistical analysis of continental geotherms reveals that thick (> 250 km) lithosphere is restricted solely to young Archean terranes (3.0–2.6 Ga), while in old Archean cratons (3.6–3.0 Ga) lithospheric roots do not extend deeper than 200–220 km. It is proposed that the former were formed by tectonic stacking and underplating during paleocollision of continental nuclei; it is likely that such exceptionally thick lithospheric roots have a limited lateral extent and are restricted to paleoterrane boundaries. This conclusion is supported by an analysis of the growth rate of the lithosphere since the Archean, which does not reveal a peak in lithospheric volume at 2.7–2.6 Ga as expected from growth curves for juvenile crust.A pronounced peak in the rate of lithospheric growth (10–18 km³/year) at 2.1–1.7 Ga (as compared to 5–8 km³/year in the Archean) well correlates with a peak in the growth of juvenile crust and with a consequent global extraction of massif-type anorthosites. It is proposed that large-scale variations in lithospheric thickness at cratonic margins and at paleoterrane boundaries controlled anorogenic magmatism. In particular, mid-Proterozoic anorogenic magmatism at the cratonic margins was caused by edge-driven convection triggered by a fast growth of the lithospheric mantle at 2.1–1.7 Ga. Belts of anorogenic magmatism within cratonic interiors can be caused by a deflection of mantle heat by a locally thickened lithosphere at paleosutures and, thus, can be surface manifestations of exceptionally thick lithospheric roots. The present volume of continental lithosphere as estimated from the new global map of lithospheric thermal thickness is 27.8 (± 7.0) × 10⁹ km³ (excluding submerged terranes with continental crust); preserved continental crust comprises ca. 7.7 × 10⁹ km³. About 50% of the present continental lithosphere existed by 1.8 Ga.     

An explosion seismology investigation of the continental margin west of the Hebrides, Scotland, at 58°N

During summer 1975, a line of large shots was fired across the continental margin between the Rockall Trough and the Hebridean shelf along 58°N. Arrivals were observed at temporary seismic stations set up across Scotland and in northwestern Ireland. A clear P₂ phase was observed to cross the margin and a converted phase P₁ also seen on the records is interpreted as travelling in the sub-sedimentary oceanic crust of Rockall Trough and in the upper continental crust beneath the shelf.The continental crust beneath the Hebridean shelf is estimated to be 27 ± 2 km thick, with P_g = 6.22 ± 0.03 km/s and P_n = 8.01 ± 0.04 km/s as determined by time-term analysis. P_g delays on the outer shelf are interpreted in terms of a seaward thickening wedge of Mesozoic sediments which pre-date the split. P_n beneath the Rockall Trough was poorly determined at 8.20 ± 0.17 km/s and the Moho is estimated to be 18 ± 2 km deep at 58°N. This and other seismic and gravity work indicates a northward thickening of the crust along the Rockall Trough, accounting for the northward decrease in the height of the slope.Our results, and those of gravity interpretations, indicate a relatively abrupt transition between continental and oceanic crust, possibly correlating with the lack of major shelf subsidence. This is attributed to a relatively cool origin for this margin. The main thinning of the continental crust beneath the slope is attributed to outslip of continental crustal material into and beneath the newly forming oceanic crust during the first few million years after the split, possibly enhanced by pre-split stretching.     

The role of subduction erosion in the generation of Andean and other convergent plate boundary arc magmas,the continental crust and mantle

Subduction erosion, which occurs at all convergent plate boundaries associated with magmatic arcs formed on crystalline forearc basement, is an important process for chemical recycling, responsible globally for the transport of ~1.7 Armstrong Units (1 AU = 1 km³/yr) of continental crust back into the mantle. Along the central Andean convergent plate margin, where there is very little terrigenous sediment being supplied to the trench as a result of the arid conditions, the occurrence of mantle-derived olivine basalts with distinctive crustal isotopic characteristics (⁸⁷Sr/⁸⁶Sr ≥ 0.7050; ε_Nd ≤ −2; ε_Hf ≤ +2) correlates spatially and/or temporally with regions and/or episodes of high rates of subduction erosion, and a strong case can be made for the formation of these basalts to be due to incorporation into the subarc mantle wedge of tectonically eroded and subducted forearc continental crust. In other convergent plate boundary magmatic arcs, such as the South Sandwich and Aleutian Islands intra-oceanic arcs and the Central American and Trans-Mexican continental margin volcanic arcs, similar correlations have been demonstrated between regions and/or episodes of relatively rapid subduction erosion and the genesis of mafic arc magmas containing enhanced proportions of tectonically eroded and subducted crustal components that are chemically distinct from pelagic and/or terrigenous trench sediments. It has also been suggested that larger amounts of melts derived from tectonically eroded and subducted continental crust, rising as diapirs of buoyant low density subduction mélanges, react with mantle peridotite to form pyroxenite metasomatites that than melt to form andesites. The process of subduction erosion and mantle source region contamination with crustal components, which is supported by both isotopic and U-Pb zircon age data implying a fast and efficient connectivity between subduction inputs and magmatic outputs, is a powerful alternative to intra-crustal assimilation in the generation of andesites, and it negates the need for large amounts of mafic cumulates to form within and then be delaminated from the lower crust, as required by the basalt-input model of continental crustal growth. However, overall, some significant amount of subducted crust and sediment is neither underplated below the forearc wedge nor incorporated into convergent plate boundary arc magmas, but instead transported deeper into the mantle where it plays a role in the formation of isotopically enriched mantle reservoirs. To ignore or underestimate the significance of the recycling of tectonically eroded and subducted continental crust in the genesis of convergent plate boundary arc magmas, including andesites, and for the evolution of both the continental crust and mantle, is to be on the wrong side of history in the understanding of these topics.     

A first coincident normal-incidence and wide-angle approach to studying the extending Aegean crust

The Aegean Sea is a broad area of submerged continental crust undergoing active extension to varying degrees. A combined near-normal incidence and wide-angle seismic recording programme was conducted in the western Aegean Sea in 1993, with the principal objective of testing the popular hypothesis that lower crustal deformation (particularily extension) is expressed as a seismically “layered lower crust” (LLC). Across the southern margin of the Cretan trough (i.e. North Cretan offshore margin), a LLC was indicated by wide-angle arrivals that was not apparent on either the coincident near-normal-in-cidence profile or on older low-frequency refraction records. North of the northern margin of the Cretan Trough, beneath the Cyclades, a domain of strong reflectivity is recorded from the middle to lower crust. Here, the near-normal incidence sections also show this typical LLC reflectivity. On the wide-angle sections, a distinct interface is suggested in addition, at a larger depth than that previously assumed for the Moho discontinuity. The structural images and interpretations derived from the new seismic data so far do not clearly support either a pure-shear crustal stretching or an asymmetric simple-shear extension model for the Aegean Sea. Our results appear to be consistent with a tectonic model, where middle crust mobilised by flow coincides spatially with upper crust that has been thinned by active extension of an orogenically thickened crust and expressed near the surface as an exhumed metamorphic core complex.     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

Magmatic evolution of the South Shetland Islands,Antarctica, and implications for continental crust formation

Lavas from the South Shetland Islands volcanic arc (northern Antarctic Peninsula) have been investigated in order to determine the age, petrogenesis and compositional evolution of a long-lived volcanic arc constructed on 32-km-thick crust, a thickness comparable with average continental crust. New ⁴⁰Ar–³⁹Ar ages for the volcanism range between 135 and 47 Ma and, together with published younger ages, confirm a broad geographical trend of decreasing ages for the volcanism from southwest to northeast. The migration pattern breaks down in Palaeogene time, with Eocene magmatism present on both Livingston and King George islands, which may be due to a change in both subduction direction and velocity after c. 60 Ma. The lavas range from tholeiitic to calc-alkaline, but there is no systematic change with age or geographic location. The compositions of lavas from the north-eastern islands indicate magma generation in a depleted mantle wedge with relatively low Sr and high Nd isotopic compositions and low U/Nb, Th/Nd and Ba/Nb ratios that was metasomatized by hydrous fluids from subducted basaltic oceanic crust. Lavas from the south-western islands show an additional sedimentary influence most likely due to fluid release from subducted sediments into the mantle wedge. Although magmatic activity in the South Shetland arc extended over c. 100 m.y., there is no evolution towards more enriched or evolved magmas with time. Few South Shetland arc lavas are sufficiently enriched with incompatible elements to provide a potential protolith for the generation of average continental crust. We conclude that even long-established subduction zones with magmatic systems founded on relatively thick crust do not necessarily form continental crustal building blocks. They probably represent only the juvenile stages of continental crust formation, and additional re-working, for example during subsequent arc-continental margin collision, is required before they can evolve into average continental crust.     

Exhumation of oceanic blueschists and eclogites in subduction zones: Timing and mechanisms

High-pressure low-temperature (HP–LT) metamorphic rocks provide invaluable constraints on the evolution of convergent zones. Based on a worldwide compilation of key information pertaining to fossil subduction zones (shape of exhumation P–T–t paths, exhumation velocities, timing of exhumation with respect to the convergence process, convergence velocities, volume of exhumed rocks,…), this contribution reappraises the burial and exhumation of oceanic blueschists and eclogites, which have received much less attention than continental ones during the last two decades.Whereas the buoyancy-driven exhumation of continental rocks proceeds at relatively fast rates at mantle depths (≥ cm/yr), oceanic exhumation velocities for HP–LT oceanic rocks, whether sedimentary or crustal, are usually on the order of the mm/yr. For the sediments, characterized by the continuity of the P–T conditions and the importance of accretionary processes, the driving exhumation mechanisms are underthrusting, detachment faulting and erosion. In contrast, blueschist and eclogite mafic bodies are systematically associated with serpentinites and/or a mechanically weak matrix and crop out in an internal position in the orogen.Oceanic crust rarely records P conditions > 2.0–2.3 GPa, which suggests the existence of maximum depths for the sampling of slab-derived oceanic crust. On the basis of natural observations and calculations of the net buoyancy of the oceanic crust, we conclude that beyond depths around 70 km there are either not enough serpentinites and/or they are not light enough to compensate the negative buoyancy of the crust.Most importantly, this survey demonstrates that short-lived (< ～ 15 My), discontinuous exhumation is the rule for the oceanic crust and associated mantle rocks: exhumation takes place either early (group 1: Franciscan, Chile), late (group 2: New Caledonia, W. Alps) or incidentally (group 3: SE Zagros, Himalayas, Andes, N. Cuba) during the subduction history. This discontinuous exhumation is likely permitted by the specific thermal regime following the onset of a young, warm subduction (group 1), by continental subduction (group 2) or by a major, geodynamic modification of convergence across the subduction zone (group 3; change of kinematics, subduction of asperities, etc).Understanding what controls this short-lived exhumation and the detachment and migration of oceanic crustal slices along the subduction channel will provide useful insights into the interplate mechanical coupling in subduction zones.     

A-type granites induced by a breaking-off and delamination of the subducted Junggar oceanic plate,West Junggar,Northwest China

The A-type granites with highly positive ε_Nd(t) values in the West Junggar, Central Asian Orogenic Belt (CAOB), have long been perceived as a group formed under the same tectonic and geodynamic setting, magmatic sourceq and petrogenetic model. Geological evidence shows that these granites occurred at two different tectonic units related to the southeastern subduction of Junggar oceanic plate: the Hongshan and Karamay granites emplaced in the southeast of West Junggar in the Baogutu continental arc; whereas the Akebasitao and Miaoergou granites formed in the accretionary prism. Here the authors present new bulk-rock geochemistry and Sr-Nd isotopes, zircon U-Pb ages and Hf-O isotopes data on these granites. The granites in the Baogutu continental arc and accretionary prism contain similar zircon ε_Hf(t) values (+10.9 to +16.2) and bulk-rock geochemical characteristics (high SiO₂ and K₂O contents, enriched LILEs (except Sr), depleted Sr, Ta and Ti, and negative anomalies in Ce and Eu). The Hongshan and Karamay granites in the Baogutu continental arc have older zircon U-Pb ages (315–305 Ma) and moderate ¹⁸O enrichments (δ¹⁸O_zircon=+6.41‰–+7.96‰); whereas the Akebasitao and Miaoergou granites in the accretionary prism have younger zircon U-Pb ages (305–301 Ma) with higher ¹⁸O enrichments (δ¹⁸O_zircon=+8.72‰–+9.89‰). The authors deduce that the elevated ¹⁸O enrichments of the Akebasitao and Miaoergou granites were probably inherited from low-temperature altered oceanic crusts. The Akebasitao and Miaoergou granites were originated from partial melting of low-temperature altered oceanic crusts with juvenile oceanic sediments below the accretionary prism. The Hongshan and Karamay granites were mainly derived from partial melting of basaltic juvenile lower crust with mixtures of potentially chemical weathered ancient crustal residues and mantle basaltic melt (induced by hot intruding mantle basaltic magma at the bottom of the Baogutu continental arc). On the other hand, the Miaoergou charnockite might be sourced from a deeper partial melting reservoir under the accretionary prism, consisting of the low-temperature altered oceanic crust, juvenile oceanic sediments, and mantle basaltic melt. These granites could be related to the asthenosphere’s counterflow and upwelling, caused by the break-off and delamination of the subducted oceanic plate beneath the accretionary prism Baogutu continental arc in a post-collisional tectonic setting.©2022 China Geology Editorial Office.     

The building blocks of continental crust: Evidence for a major change in the tectonic setting of continental growth at the end of the Archean

Oceanic arcs are commonly cited as primary building blocks of continents, yet modern oceanic arcs are mostly subducted. Also, lithosphere buoyancy considerations show that oceanic arcs (even those with a felsic component) should readily subduct. With the exception of the Arabian–Nubian orogen, terranes in post-Archean accretionary orogens comprise < 10% of accreted oceanic arcs, whereas continental arcs compose 40–80% of these orogens. Nd and Hf isotopic data suggest that accretionary orogens include 40–65% juvenile crustal components, with most of these (> 50%) produced in continental arcs.Felsic igneous rocks in oceanic arcs are depleted in incompatible elements compared to average continental crust and to felsic igneous rocks from continental arcs. They have lower Th/Yb, Nb/Yb, Sr/Y and La/Yb ratios, reflecting shallow mantle sources in which garnet did not exist in the restite during melting. The bottom line of these geochemical differences is that post-Archean continental crust does not begin life in oceanic arcs. On the other hand, the remarkable similarity of incompatible element distributions in granitoids and felsic volcanics from continental arcs is consistent with continental crust being produced in continental arcs.During the Archean, however, oceanic arcs may have been thicker due to higher degrees of melting in the mantle, and oceanic lithosphere would be more buoyant. These arcs may have accreted to each other and to oceanic plateaus, a process that eventually led to the production of Archean continental crust. After the Archean, oceanic crust was thinner due to cooling of the mantle and less melt production at ocean ridges, hence, oceanic lithosphere is more subductable. Widespread propagation of plate tectonics in the late Archean may have led not only to rapid production of continental crust, but to a change in the primary site of production of continental crust, from accreted oceanic arcs and oceanic plateaus in the Archean to primarily continental arcs thereafter.     

A general model for the intrusion and evolution of 'mantle' garnet peridotites in high-pressure and ultra-high-pressure metamorphic terranes

Garnet‐bearing peridotite lenses are minor but significant components of most metamorphic terranes characterized by high‐temperature eclogite facies assemblages. Most peridotite intrudes when slabs of continental crust are subducted deeply (60–120 km) into the mantle, usually by following oceanic lithosphere down an established subduction zone. Peridotite is transferred from the resulting mantle wedge into the crustal footwall through brittle and/or ductile mechanisms. These ‘mantle’ peridotites vary petrographically, chemically, isotopically, chronologically and thermobarometrically from orogen to orogen, within orogens and even within individual terranes. The variations reflect: (1) derivation from different mantle sources (oceanic or continental lithosphere, asthenosphere); (2) perturbations while the mantle wedges were above subducting oceanic lithosphere; and (3) changes within the host crustal slabs during intrusion, subduction and exhumation. Peridotite caught within mantle wedges above oceanic subduction zones will tend to recrystallize and be contaminated by fluids derived from the subducting oceanic crust. These ‘subduction zone peridotites’ intrude during the subsequent subduction of continental crust. Low‐pressure protoliths introduced at shallow (serpentinite, plagioclase peridotite) and intermediate (spinel peridotite) mantle depths (20–50 km) may be carried to deeper levels within the host slab and undergo high‐pressure metamorphism along with the enclosing rocks. If subducted deeply enough, the peridotites will develop garnet‐bearing assemblages that are isofacial with, and give the same recrystallization ages as, the eclogite facies country rocks. Peridotites introduced at deeper levels (50–120 km) may already contain garnet when they intrude and will not necessarily be isofacial or isochronous with the enclosing crustal rocks. Some garnet peridotites recrystallize from spinel peridotite precursors at very high temperatures (c. 1200 °C) and may derive ultimately from the asthenosphere. Other peridotites are from old (>1 Ga), cold (c. 850 °C), subcontinental mantle (‘relict peridotites’) and seem to require the development of major intra‐cratonic faults to effect their intrusion.     

Cenozoic geodynamic evolution of the Aegean

The Aegean region is a concentrate of the main geodynamic processes that shaped the Mediterranean region: oceanic and continental subduction, mountain building, high-pressure and low-temperature metamorphism, backarc extension, post-orogenic collapse, metamorphic core complexes, gneiss domes are the ingredients of a complex evolution that started at the end of the Cretaceous with the closure of the Tethyan ocean along the Vardar suture zone. Using available plate kinematic, geophysical, petrological and structural data, we present a synthetic tectonic map of the whole region encompassing the Balkans, Western Turkey, the Aegean Sea, the Hellenic Arc, the Mediterranean Ridge and continental Greece and we build a lithospheric-scale N-S cross-section from Crete to the Rhodope massif. We then describe the tectonic evolution of this cross-section with a series of reconstructions from ~70 Ma to the Present. We follow on the hypothesis that a single subduction has been active throughout most of the Mesozoic and the entire Cenozoic, and we show that the geological record is compatible with this hypothesis. The reconstructions show that continental subduction (Apulian and Pelagonian continental blocks) did not induce slab break-off in this case. Using this evolution, we discuss the mechanisms leading to the exhumation of metamorphic rocks and the subsequent formation of extensional metamorphic domes in the backarc region during slab retreat. The tectonic histories of the two regions showing large-scale extension, the Rhodope and the Cyclades are then compared. The respective contributions to slab retreat, post-orogenic extension and lower crust partial melting of changes in kinematic boundary conditions and in nature of subducting material, from continental to oceanic, are discussed.