Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

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.     

Late Quaternary activity and dextral strike-slip movement on the Karakax Fault Zone, northwest Tibet

Field observations and interpretations of satellite images reveal that the westernmost segment of the Altyn Tagh Fault (called Karakax Fault Zone) striking WNW located in the northwestern margin of the Tibetan Plateau has distinctive geomorphic and tectonic features indicative of right-lateral strike-slip fault in the Late Quaternary. South-flowing gullies and N–S-trending ridges are systematically deflected and offset by up to ~ 1250 m, and Late Pleistocene–Holocene alluvial fans and small gullies that incise south-sloping fans record dextral offset up to ~ 150 m along the fault zone. Fault scarps developed on alluvial fans vary in height from 1 to 24 m. Riedel composite fabrics of foliated cataclastic rocks including cataclasite and fault gouge developed in the shear zone indicate a principal right-lateral shear sense with a thrust component. Based on offset Late Quaternary alluvial fans, ¹⁴C ages and composite fabrics of cataclastic fault rocks, it is inferred that the average right-lateral strike-slip rate along the Karakax Fault Zone is ~ 9 mm/a in the Late Quaternary, with a vertical component of ~ 2 mm/a, and that a M 7.5 morphogenic earthquake occurred along this fault in 1902. We suggest that right-lateral slip in the Late Quaternary along the WNW-trending Karakax Fault Zone is caused by escape tectonics that accommodate north–south shortening of the western Tibetan Plateau due to ongoing northward penetration of the Indian plate into the Eurasian plate.     

P-wave velocity structure under the Changbaishan volcanic region, NE China, from wide-angle reflection and refraction data

We present velocity models determined by inverting refracted and reflected arrivals along two active source lines in the Changbaishan volcanic region, NE China. We resolve a prominent low-velocity zone (LVZ) in the crust, with velocities as low as 5.4 km/s. Away from the LVZ, the velocity gradients in the crust are relatively smooth, with average P-wave velocities of about 6.0–6.5 km/s. The Moho is at about 35 km depth, thickening to about 40 km under the Tianchi volcano, and thinning to about 30 km under the LVZ. The LVZ is located about 30–60 km to the north of the summit of the Tianchi volcano (the most recently active volcano in the region), is about 30–75 km in north–south extent, is at most 35 km in east–west extent, and is in the depth range of about 10–25 km below the surface. We use these results to constrain receiver function inversions, and show that the receiver functions in the region are compatible with our findings. With these data alone, the significance of the LVZ in non-unique, although we do not see any evidence to support the presence of partial melt in the crust, and favor the interpretation that the LVZ indicates a residual crustal magma chamber.     

A new kinematic evolutionary model for the growth of a duplex — an example from the Rangit duplex, Sikkim Himalaya, India

The Lesser Himalayan duplex (LHD) is a prominent structure through much of the Lesser Himalayan fold–thrust belt. In the Darjeeling - Sikkim Himalaya a component of the LHD is exposed in the Rangit window as the Rangit duplex (RD). The RD consists of ten horses of the upper Lesser Himalayan Sequence (Gondwana, Buxa, Upper Daling). The duplex varies from hinterland-dipping in the north, through an antiformal stack in the middle to foreland-dipping in the south. The Ramgarh thrust (RT) is the roof thrust and, based on a balanced cross-section, the Main Himalayan Sole thrust is the floor thrust at a depth of ~ 10 km and with a dip of ~ 3.5° N.Retrodeformation suggests that the RD initiated as a foreland-dipping duplex with the Early Ramgarh thrust as the roof thrust and the RT as the floor thrust. The RT became the roof thrust during continued duplexing by a combination of footwall imbrication and concurrent RT reactivation. This kinematic history best explains the large translation of the overlying MCT sheets. The restoration suggests that RD shortening is ~ 125 km, and the original Gondwana basin extended ~ 142 km northward of its present northernmost exposures within the window.     

Crustal structure of mainland China from deep seismic sounding data

Since 1958, about ninety seismic refraction/wide angle reflection profiles, with a cumulative length of more than sixty thousand kilometers, have been completed in mainland China. We summarize the results in the form of (1) a new contour map of crustal thickness, (2) fourteen representative crustal seismic velocity–depth columns for various tectonic units, and, (3) a Pn velocity map. We found a north–south-trending belt with a strong lateral gradient in crustal thickness in central China. This belt divides China into an eastern region, with a crustal thickness of 30–45 km, and a western region, with a thickness of 45–75 km. The crust in these two regions has experienced different evolutionary processes, and currently lies within distinct tectonic stress fields. Our compilation finds that there is a high-velocity (7.1–7.4 km/s) layer in the lower crust of the stable Tarim basin and Ordos plateau. However, in young orogenic belts, including parts of eastern China, the Tianshan and the Tibetan plateau, this layer is often absent. One exception is southern Tibet, where the presence of a high-velocity layer is related to the northward injection of the cold Indian plate. This high-velocity layer is absent in northern Tibet. In orogenic belts, there usually is a low-velocity layer (LVL) in the crust, but in stable regions this layer seldom exists. The Pn velocities in eastern China generally range from 7.9 to 8.1 km/s and tend to be isotropic. Pn velocities in western China are more variable, ranging from 7.7 to 8.2 km/s, and may display azimuthal anisotropy.     

The post-collisional stabilization of thickened crust in the Southern Urals

The Urals are characterized by a depression of the Moho to a depth of 57 km. This structure is interpreted as a relic orogenic root, which has been conserved because no significant post-collisional processes occurred. However, there is evidence that voluminous post-collisional magmatism affected the lower crust. In this paper, we use thermal finite element models to quantify the influence of the post-collisional magmatism on the stabilization of the root. We show that at least 70% of the heat producing elements migrated in granitic melts from the lower crust to the upper crust. As a result the crustal heat flow reduced and the lithosphere could stabilize at a thickness of 180 km. Furthermore, we propose that a granulite metamorphic event during the thermal relaxation of the collision zone prevented the 57 km thick crust from delamination. These results strongly indicate that post-collisional processes were necessary for the stabilization of the Uralian crust and lithosphere.     

Ultrahigh-temperature metamorphism in the Achankovil Zone: Implications for the correlation of crustal blocks in southern India

The Achankovil Zone of southern India, a NW–SE trending lineament of 8–10 km in width and > 100 km length, is a kinematically debated crustal feature, considered to mark the boundary between the Madurai Granulite Block in the north and the Trivandrum Granulite Block in the south. Both these crustal blocks show evidence for ultrahigh-temperature metamorphism during the Pan-African orogeny, although the exhumation styles are markedly different. The Achankovil Zone is characterized by discontinuous strands of cordierite-bearing gneiss with an assemblage of cordierite + garnet + quartz + plagioclase + spinel + ilmenite + magnetite ± orthopyroxene ± biotite ± K-feldspar ± sillimanite. The lithology preserves several peak and post-peak metamorphic assemblages including: (1) orthopyroxene + garnet, (2) perthite and/or anti-perthite, (3) cordierite ± orthopyroxene corona around garnet, and (4) cordierite + quartz symplectite after garnet. We estimate the peak metamorphic conditions of these rocks using orthopyroxene-bearing geothermobarometers and feldspar solvus which yield 8.5–9.5 kbar and 940–1040 °C, the highest P–T conditions so far recorded from the Achankovil Zone. The retrograde conditions were obtained from cordierite-bearing geothermobarometers at 3.5–4.5 kbar and 720 ± 60 °C. From orthopyroxene chemistry, we record a multistage exhumation history for these rocks, which is closely comparable with those reported in recent studies from the Madurai Granulite Block, but different from those documented from the Trivandrum Granulite Block. An evaluation of the petrologic and geochronologic data, together with the nature of exhumation paths leads us to propose that the Achankovil Zone is probably the southern flank of the Madurai Granulite Block, and not a unit of the Trivandrum Granulite Block as presently believed. Post-tectonic alkali granites that form an array of “suturing plutons” along the margin of the Madurai Granulite Block and within the Achankovil Zone, but are absent in the Trivandrum Granulite Block, suggest that the boundary between the Madurai Granulite Block and the Trivandrum Granulite Block might lie along the Tenmalai shear zone at the southern extremity of the Achankovil Zone.     

P-wave velocity features of the lithosphere under the Eromanga Basin, Eastern Australia, including a prominent MID-crustal (Conrad?) discontinuity

The crustal structure of the central Eromanga Basin in the northern part of the Australian Tasman Geosyncline, revealed by coincident seismic reflection and refraction shooting, contrasts with some neighbouring regions of the continent. The depth to the crust-mantle boundary (Moho) of 36–41 km is much less than that under the North Australian Craton to the northwest (50–55 km) and the Lachlan Fold Belt to the southeast (43–51 km) but is similar to that under the Drummond and Bowen Basins to the east.The seismic velocity boundaries within the crust are sharp compared with the transitional nature of the boundaries under the North Australian and Lachlan provinces. In particular, there is a sharp velocity increase at mid-crustal depths (21–24 km) which has not been observed with such clarity elsewhere in Australia (the Conrad discontinuity?).In the lower crust, the many discontinuous sub-horizontal reflections are in marked contrast to lack of reflecting horizons in the upper crust, further emphasising the differences between the upper and lower crust. The crust-mantle boundary (Moho) is characterised by an increase in velocity from 7.1–7.7 km/s to a value of 8.15 + 0.04 km/s. The depth to the Moho under the Canaway Ridge, a prominent basement high, is shallower by about 5 km than the regional Moho depth; there is also no mid-crustal horizon under the Canaway Ridge but there is a very sharp velocity increase at the Moho depth of 34 km. The Ridge could be interpreted as a horst structure extending to at least Moho depths but it could also have a different intra-crustal structure from the surrounding area.The sub-crustal lithosphere has features which have been interpreted, from limited data, as being caused by a velocity gradient at 56–57 km depth with a low velocity zone above it.Because of the contrasting crustal thicknesses and velocity gradients, the lithosphere of the central Eromanga Basin cannot be considered as an extension of the exposed Lachlan Fold Belt or the North Australian Craton. The lack of seismic reflections from the upper crust indicates no coherent accoustic impedance pattern at wavelengths greater than 100 m, consistent with an upper crustal basement of tightly folded meta-sedimentary and meta-volcanic rocks. The crustal structure is consistent with a pericratonic or arc/back-arc basin being cratonised in an episode of convergent tectonics in the Early Palaeozoic. The seismic reflections from the lower crust indicate that it could have developed in a different tectonic environment.     

New evidence from the western Canadian Arctic Archipelago for the resubmergence of Bering Strait

Widespread molluscan samples were collected from raised marine sediments to date the last retreat of the NW Laurentide Ice Sheet from the western Canadian Arctic Archipelago. At the head of Mercy Bay, northern Banks Island, deglacial mud at the modern coast contains Hiatella arctica and Portlandia arctica bivalves, as well as Cyrtodaria kurriana, previously unreported for this area. Multiple H. arctica and C. kurriana valves from this site yield a mean age of 11.5 ¹⁴C ka BP (with 740 yr marine reservoir correction). The occurrence of C. kurriana, a low Arctic taxon, raises questions concerning its origin, because evidence is currently lacking for a molluscan refugium in the Arctic Ocean during the last glacial maximum. Elsewhere, the oldest late glacial age available on C. kurriana comes from the Laptev Sea where it is < 10.3 ¹⁴C ka BP and attributed to a North Atlantic source. This is 2000 cal yr younger than the Mercy Bay samples reported here, making the Laptev Sea, ~ 3000 km to the west, an unlikely source. An alternate route from the North Atlantic into the Canadian Arctic Archipelago was precluded by coalescent Laurentide, Innuitian and Greenland ice east of Banks Island until ~ 10 ¹⁴C ka BP. We conclude that the presence of C. kurriana on northern Banks Island records migration from the North Pacific. This requires the resubmergence of Bering Strait by 11.5 ¹⁴C ka BP, extending previous age determinations on the reconnection of the Pacific and Arctic oceans by up to 1000 yr. This renewed ingress of Pacific water likely played an important role in re-establishing Arctic Ocean surface currents, including the evacuation of thick multi-year sea ice into the North Atlantic prior to the Younger Dryas geochron.     

Mantle flow-induced crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland (E Turkey) deduced from the geological and geophysical data

Eastern Anatolia consisting of an amalgamation of fragments of oceanic and continental lithosphere is a current active intercontinental contractional zone that is still being squeezed and shortened between the Arabian and Eurasian plates. This collisional and contractional zone is being accompanied by the tectonic escape of most of the Anatolian plate to the west by major strike-slip faulting on the right-lateral North Anatolian Transform Fault Zone (NATFZ) and left-lateral East Anatolian Transform Fault Zone (EATFZ) which meet at Karlıova forming an east-pointing cusp. The present-day crust in the area between the easternmost part of the Anatolian plate and the Arabian Foreland gets thinner from north (ca 44 km) to south (ca 36 km) relative to its eastern (EAHP) and western sides (central Anatolian region). This thinner crustal area is characterized by shallow CPD (12–16 km), very low Pn velocities (< 7.8 km/s) and high Sn attenuation which indicate partially molten to eroded mantle lid or occurrence of asthenospheric mantle beneath the crust. Northernmost margin of the Arabian Foreland in the south of the Bitlis–Pötürge metamorphic gap area is represented by moderate CPD (16–18 km) relative to its eastern and western sides, and low Pn velocities (8 km/s). We infer from the geophysical data that the lithospheric mantle gets thinner towards the Bitlis–Pötürge metamorphic gap area in the northern margin of the Arabian Foreland which has been most probably caused by mechanical removal of the lithospheric mantle during mantle invasion to the north following the slab breakoff beneath the Bitlis–Pötürge Suture Zone. Mantle flow-driven rapid extrusion and counterclockwise rotation of the Anatolian plate gave rise to stretching and hence crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland which is currently dominated by wrench tectonics.     

Gravity signals from the lithosphere in the Central European Basin System

We study the gravity signals from different depth levels in the lithosphere of the Central European Basin System (CEBS). The major elements of the CEBS are the Northern and Southern Permian Basins which include the Norwegian–Danish Basin (NDB), the North-German Basin (NGB) and the Polish Trough (PT). An up to 10 km thick sedimentary cover of Mesozoic–Cenozoic sediments, hides the gravity signal from below the basin and masks the heterogeneous structure of the consolidated crust, which is assumed to be composed of domains that were accreted during the Paleozoic amalgamation of Europe. We performed a three-dimensional (3D) gravity backstripping to investigate the structure of the lithosphere below the CEBS.Residual anomalies are derived by removing the effect of sediments down to the base of Permian from the observed field. In order to correct for the influence of large salt structures, lateral density variations are incorporated. These sediment-free anomalies are interpreted to reflect Moho relief and density heterogeneities in the crystalline crust and uppermost mantle. The gravity effect of the Moho relief compensates to a large extent the effect of the sediments in the CEBS and in the North Sea. Removal of the effects of large-scale crustal inhomogeneities shows a clear expression of the Variscan arc system at the southern part of the study area and the old crust of Baltica further north–east. The remaining residual anomalies (after stripping off the effects of sediments, Moho topography and large-scale crustal heterogeneities) reveal long wavelength anomalies, which are caused mainly by density variations in the upper mantle, though gravity influence from the lower crust cannot be ruled out. They indicate that the three main subbasins of the CEBS originated on different lithospheric domains. The PT originated on a thick, strong and dense lithosphere of the Baltica type. The NDB was formed on a weakened Baltica low-density lithosphere formed during the Sveco-Norwegian orogeny. The major part of the NGB is characterized by high-density lithosphere, which includes a high-velocity lower crust (relict of Baltica passive margin) overthrusted by the Avalonian terrane. The short wavelength pattern of the final residuals shows several north–west trending gravity highs between the Tornquist Zone and the Elbe Fault System. The NDB is separated by a gravity low at the Ringkøbing–Fyn high from a chain of positive anomalies in the NGB and the PT. In the NGB these anomalies correspond to the Prignitz (Rheinsberg anomaly), the Glueckstadt and Horn Graben, and they continue further west into the Central Graben, to join with the gravity high of the Central North Sea.     

Morphological characteristics and emplacement mechanism of the seamounts in the Central Indian Ocean Basin

The morphotectonic features of the Central Indian Ocean Basin (CIOB) provide information regarding the development of the basin. Multibeam mapping of the CIOB reveals presence of abundant isolated seamounts and seamount chains sub-parallel to each other and major fracture zones along 73° E, 79° E and 75°45′ E. Morphological analyses were carried out for 200 seamounts that occur either as isolated edifies or along eight sub-parallel chains. The identified eight parallel seamount chains that trend almost north–south and reflecting the absolute motion of the Indian plate, probably originated from the ancient propagative fractures. Inspite of the differences in their height, the seamounts of these eight chains are morphologically correlatable. In the study area the seamounts are clustered north and south of 12° S latitude. Interestingly, in the area north of 12° S (area II: 9°–12° S) the seamounts are distinctly smaller (≤ 400 m height) whereas, the area south of 12° S (area I: 12°–15° S) has a mixed population of seamounts. The normalized abundance of the CIOB seamount is 976 seamounts/10⁶ km² but on a finer scale this value varies from 500 to 1600 seamounts/10⁶ km², which is less than the seamount concentrations of the Pacific and Atlantic oceans (9000 to 16,000 seamounts/10⁶ km²). Three categories of seamounts are present in the CIOB e.g. (1) single-peaked (2) multi-peaked and (3) composite. The study indicate that single-peaked seamounts are dominant (89%) while multi-peaked is less (8%) and composite ones are rare (3%) in the CIOB.The progressive northward movement of the Indian continent caused collision between India and Asia at around 62 Ma ago. A majority of the near-axis originated seamounts in the CIOB seemed to have formed as a consequence of the temporally widespread (Cretaceous  65 Ma to late Eocene < 49 Ma) collision between India and Eurasia. The regional stress patterns in the Indian plate vary N to NE in the continent and N to NW in Indian Ocean areas. The combined effect of the regional stress patterns maintained the orientation of the seamount chains and the local stress regime helped in the upwelling of magma and formation of seamounts. The low heat flow, morphological features and geochemical signature indicate that the morphotectonic structures formed contemporaneously with the oceanic crust.     

Dynamic recrystallization near the brittle-plastic transition in naturally and experimentally deformed quartz aggregates

The electron backscattering diffraction technique (EBSD) was used to analyze bulging recrystallization microstructures from naturally and experimentally deformed quartz aggregates, both of which are characterized by porphyroclasts with finely serrated grain boundaries and grain boundary bulges set in a matrix of very fine recrystallized grains. For the Tonale mylonites we investigated, a temperature range of 300–380 °C, 0.25 GPa confining pressure, a flow stress range of ~ 0.1–0.2 GPa, and a strain rate of ~ 10^− 13 s^− 1 were estimated. Experimental samples of Black Hills quartzite were analyzed, which had been deformed in axial compression at 700 °C, 1.2–1.5 GPa confining pressure, a flow stress of ~ 0.3–0.4 GPa, a strain rate of ~ 10^− 6 s^− 1, and to 44% to 73% axial shortening. Using orientation imaging we investigated the dynamic recrystallization microstructures and discuss which processes may contribute to their development. Our results suggest that several deformation processes are important for the dismantling of the porphyroclasts and the formation of recrystallized grains. Grain boundary bulges are not only formed by local grain boundary migration, but they also display a lattice misorientation indicative of subgrain rotation. Dynamic recrystallization affects especially the rims of host porphyroclasts with a hard orientation, i.e. with an orientation unsuitable for easy basal slip. In addition, Dauphiné twins within porphyroclasts are preferred sites for recrystallization. We interpret large misorientation angles in the experimental samples, which increase with increasing strain, as formed by the activity of fluid-assisted grain boundary sliding.     

Comparison of gravity anomaly between mature and immature intra-oceanic subduction zones in the western Pacific

The western Pacific hosts major subduction systems such as Izu–Bonin–Mariana and Tonga–Kermadec, but also less conspicuous systems such as Yap, Mussau and Hjort trenches which constitute the young, incomplete, or ultraslow-member in the evolutionary spectrum of subduction zones. We used satellite-derived gravity data to compare well-developed and immature subduction systems. It is shown that at spatial resolution > 10–20 km or so, the satellite data have accuracy comparable to ship-board gravity measurements over intra-oceanic subduction zones. In the isostatic residual gravity anomaly map, the width of non-isostatically-compensated region of the mature subduction zones is much wider than that of immature ones. More importantly, when the gravitational attraction due to seafloor is removed, a large difference exists between the mature and immature subduction zones in the overriding plate side. Mature subduction zones exhibit broad low gravity anomalies of ~ 200–250 mGal centered at distances of 150–200 km from the trench which are not found over immature subduction zones. The cause of the broad low gravity anomalies over mature subduction zones is debatable due to lack of information on the deep crust and upper mantle structure and property. We discuss the following four causes: (1) serpentinization of the upper mantle beneath the forearc; (2) presence of partial melt in the mantle wedge caused by release of volatiles from the slab, frictional heating and distributed by mantle circulation; (3) difference in density structure between the overriding and subducting plates caused by difference in age and thermal structures with and without compositional stratification between crust and mantle; and (4) anomalous thickness of the arc not explained by isostasy. Our analysis suggests that serpentinization cannot explain the observed gravity anomaly which appears ~ 150–200 km from the trench. Although the extent and distribution of partial melt within the mantle wedge remain in question, to our best estimate, partial melting contributes little (< 50 mGal) to the total negative gravity anomaly. The difference in density structure reflecting temperature difference can only explain less than half of the low gravity anomaly. The sinking of lighter crustal material produces a large negative anomaly in the forearc but its location does not match the observed gravity anomaly. It appears that one cannot explain the total difference in gravity anomaly without invoking anomalous thickness of the arc. Although we could not identify the sole or combination of factors that give rise to the low gravity anomaly in mature subduction zones, the comparison of gravity anomalies between mature and immature subduction zones is likely to provide an important constraint for understanding the evolution and structure of subduction zones as more complementary evidences become available.     

Variable microstructure of peridotite samples from the southern Mariana Trench: Evidence of a complex tectonic evolution

We retrieved samples of peridotite from a dredge haul (KH92-1-D2) collected during Cruise KH92-1 undertaken by the research vessel (R/V) Hakuho in 1992 at the landward trench slope of the southern Mariana Trench (11°41.16′N, 143°29.62′E; depth 6594–7431 m), which is the deepest ocean in the world. Ten of 30 retrieved samples possessed both a foliation and lineation, as assessed from 46 thin sections of various orientations and observations of hand samples. The samples showed marked variation in microstructure, ranging from coarse (> 5 mm) equigranular and intensely elongated textures to finer (< 1 mm) porphyroclastic and fine-grained equigranular textures. Olivine fabrics also varied among the different samples, with (010)[100] and (010)[001] patterns (termed A- and B-type, respectively) observed in samples with coarse textures and no clear patterns observed in samples with fine textures. Even though the peridotite samples were retrieved from a single dredge site, some contain primary tectonic microstructures and some contain secondary microstructures. Recent bathymetric and topographic analyses indicate that the lithosphere in this region is as thin as 20 km. Such a thin lithosphere may have been intensely deformed, even perhaps in the ductile regime, during fore-arc extension; consequently, the observed variations in microstructure within the peridotite samples probably reflect the complex tectonic evolution of the southern Mariana region.     

Folding of Europa's icy lithosphere: an analysis of viscous-plastic buckling and subsequent topographic relaxation

Potential contractional folds on Jupiter's icy moon Europa have been identified. The best example is at the extensional band Astypalaea Linea, where a series of subtle topographic undulations, 25 km in wavelength, possess parasitic tectonic features that support a folding origin. A scenario has been qualitatively proposed, whereby folds form via unstable contraction of the icy lithosphere, compensate for extension elsewhere on Europa, and then subsequently relax. Here, we quantitatively address this scenario, applying a model for viscous-plastic buckling of planetary lithospheres and finite element simulations of topographic relaxation. Our results suggest that the lithosphere of Europa could indeed be unstable, but the low required surface temperatures limit fold formation to higher latitudes, and the high required driving stresses (9–10 MPa) are difficult to achieve on the satellite. The depth to the brittle–ductile transition is well constrained, and high thermal gradients are indicated, implying heat flows near 100 mW m⁻². In addition, topographic relaxation progresses so slowly even at these heat flows that it is not a viable mechanism to eliminate such features over the age of Europa's surface. Given the paucity of identified folds, we conclude that the necessary conditions for their formation are rare and that lithospheric folding is a minor mechanism for compensating the large amounts of extension seen elsewhere on Europa.     

Porphyry Cu–Au and associated polymetallic Fe–Cu–Au deposits in the Beiya Area, western Yunnan Province, south China

The Alkaline porphyries in the Beiya area are located east of the Jinshajiang suture, as part of a Cenozoic alkali-rich porphyry belt in western Yunnan. The main rock types include quartz-albite porphyry, quartz-K-feldspar porphyry and biotite–K-feldspar porphyry. These porphyries are characterised by high alkalinity [(K₂O + Na₂O)% > 10%], high silica (SiO₂% > 65%), high Sr (> 400 ppm) and ⁸⁷Sr/⁸⁶Sr (> 0.706)] ratio and were intruded at 65.5 Ma, between 25.5 to 32.5 Ma, and about 3.8 Ma, respectively. There are five main types of mineral deposits in the Beiya area: (1) porphyry Cu–Au deposits, (2) magmatic Fe–Au deposits, (3) sedimentary polymetallic deposits, (4) polymetallic skarn deposits, and (5) palaeoplacers associated with karsts. The porphyry Cu–Au and polymetallic skarn deposits are associated with quartz–albite porphyry bodies. The Fe–Au and polymetallic sedimentary deposits are part of an ore-forming system that produced considerable Au in the Beiya area, and are characterised by low concentrations of La, Ti, and Co, and high concentrations of Y, Yb, and Sc.The Cenozoic porphyries in western Yunnan display increased alkalinity away from the Triassic Jinshajiang suture. Distribution of both the porphyries and sedimentary deposits in the Beiya area are interpreted to be related to partial melting in a disjointed region between upper mantle lithosphere of the Yangtze Plate and Gondwana continent, and lie within a shear zone between buried Palaeo-Tethyan oceanic lithosphere and upper mantle lithosphere, caused by the subduction and collision of India and Asia.     

Cenozoic uplift and subsidence in the North Atlantic region: Geological evidence revisited

The topographic evolution of the “passive” margins of the North Atlantic during the last 65 Myr is the subject of extensive debate due to inherent limitations of the geological, geomorphological and geophysical methods used for studies of uplift and subsidence. We have compiled a database of sign, time and amplitude (where possible) of topographic changes in the North Atlantic region during the Cenozoic (65–0 Ma). Our compilation is based on published results from reflection seismic studies, AFT (apatite fission track) studies, VR (vitrinite reflectance) trends, maximum burial, sediment supply studies, mass balance calculations and extrapolation of seismic profiles to onshore geomorphological features. The integration of about 200 published results reveal a clear pattern of topographic changes in the North Atlantic region during the Cenozoic: (1) The first major phase of Cenozoic regional uplift occurred in the late Palaeocene–early Eocene (ca 60–50 Ma), probably related to the break-up of the North Atlantic between Europe and Greenland, as indicated by the northward propagation of uplift. It was preceded by middle Palaeocene uplift and over-deepening of some basins of the North Sea and the surrounding areas. (2) A regional increase in subsidence in the offshore marginal areas of Norway, the northern North Sea, the northern British Isles and west Greenland took place in the Eocene (ca 57–35 Ma). (3) The Oligocene and Miocene (35–5 Ma) were characterized by regional tectonic quiescence, with only localised uplift, probably related to changes in plate dynamics. (4) The second major phase of regional uplift that affected all marginal areas of the North Atlantic occurred in the Plio-Pleistocene (5–0 Ma). Its amplitude was enhanced by erosion-driven glacio-isostatic compensation. Despite inconclusive evidence, this phase is likely to be ongoing at present.     

Poly-orogenic multi-stage metamorphic evolution inferred via P–T pseudosections: An example from Aspromonte Massif basement rocks (Southern Calabria, Italy)

The metamorphic evolution of a key sector of the western Mediterranean internal Alpine orogenic belt (southern Calabrian Peloritani Orogen) is identified and described by means of P–T pseudosections calculated for selected metapelite specimens, showing evidence of multi-stage metamorphism.Attention focused on the two lowermost basement nappes of the Aspromonte Massif (southern Calabria), which were differently affected by poly-orogenic multi-stage evolution. After a complete Variscan orogenic cycle, the upper unit (Aspromonte Peloritani Unit) was involved in a late-Alpine shearing event. In contrast, the several underlying metapelite slices, here grouped together as Lower Metapelite Group, show exclusive evidence of a complete Alpine orogenic cycle.In order to obtain reliable P–T constraints, an integrated approach was employed, based on: a) garnet isopleth thermobarometry; and b) theoretical predictions of the P–T stability fields of representative equilibrium assemblages. This approach, which takes into account the role of the local equilibrium volumes in controlling textural developments, yielded reliable information about P–T conditions from early to peak metamorphic stages, as well as estimates of the retrograde trajectory in the pseudosection P–T space.According to inferred detailed P–T paths, the evolution of the Aspromonte Peloritani Unit is characterised by a multi-stage Variscan cycle, subdivided into an early crustal thickening stage with P–T conditions ranging from 0.56 ± 0.05 GPa at 570 ± 10 °C to 0.63–0.93 GPa at 650–710 °C (peak conditions) and evolving to a later crustal thinning episode in lower P–T conditions (0.25 GPa at 540 °C), as documented by the retrograde trajectory.Conversely, the prograde evolution of the rocks of the Lower Metapelite Group shows evidence of a HP-LT early Alpine multi-stage cycle, with P–T evolving from 0.75–0.90 GPa at 510–530 °C towards peak conditions, with pressure increasing northwards from 1.12 ± 0.02 GPa to 1.24 ± 0.02 GPa, and temperatures of 540–570 °C.A late-Alpine mylonitic overprint affected the rocks of both the Aspromonte Peloritani Unit and the Lower Metapelite Group. This overprint was characterised by an initial retrograde decompression trajectory (0.75 ± 0.05 GPa at 570–600 °C), followed by a joint cooling history, ranging from 0.38 ± 0.14 at temperature from 450 to 520 °C.These inferred results were then used: a) to interpret the Lower Metapelite Group as a single crystalline basement unit exclusively affected by a complete Alpine orogenic cycle, according to the very similar features of P–T paths, comparable petrography and analogous structural characteristics; b) as a tool for more reliable correlations between the Aspromonte Massif, the other Calabrian terranes and the north African Orogenic Complexes. They may therefore consider a contribution to the geodynamic modelling of the western Mediterranean.