Numerical Modeling of Basin-Range Tectonics Related to Continent-Continent Collision

Continent-continent collision is the most important driving mechanism for the occurrence of various geological processes in the continental lithosphere. How to recognize and determine continent-continent collision,especially its four-dimensional temporal-spatial evolution, is a subject that geological communities have long been concerned about and studied. Continent-continent collision is mainly manifested by strong underthrnsting (subduction) of the underlying block along an intracontinental subduction zone and continuous obduction (thrusting propagation) of the overlying block along the intracontinental subduction zone, the occurrence of a basin-range tectonic framework in a direction perpendicular to the subduction zone and the flexure and disruption of the Moho. On the basis of numerical modeling, the authors discuss in detail the couplings between various amounts and rates of displacement caused by basin subsidence, mountain uplift and Moho updoming and downflexure during obduction (thrusting propagation) and subduction and the migration pattern of basin centers. They are probably indications or criteria for judgment or determination of continent-continent collision.     

Flexural and gravity modelling of the Mérida Andes and Barinas–Apure Basin, Western Venezuela

The kinematic evolution of the Barinas–Apure Basin and the southern Mérida Andes from Lower Miocene to the Present is numerically modelled using flexural isostatic theory and geophysical and geological data. Two published regional transects are used to build up a reference section, which is then used to constrain important parameters (e.g. shortenings and sedimentary thicknesses) for the flexural modelling. To control the location of the main fault system in the flexural model earthquake information is also used. The estimated flexural elastic thickness of the South American lithosphere beneath the Barinas–Apure Basin and the Mérida Andes Range is 25 km. The value for the final total shortening is 60 km. The flexural isostatic model shows that the Andean uplift has caused the South American lithosphere subsidence and the development of the Barinas–Apure Basin.In addition, gravity modelling was used to understand deep crustal features that could not be predicted by flexural theory. Consequently, the best-fit flexural model is used to build a gravity model across the Mérida Andes and the Barinas–Apure Basin preserving the best-controlled structural features from the flexural modelling (e.g. basin wavelength and depth) and slightly changing the main bodies density values and deep crustal structures. The final gravity model is intended to be representative of the major features affecting the gravity field in the study area. The predicted morphology in the lower crustal level of the final gravity model favours the hypothesis of a present delamination or megathrust of the Maracaibo crust over the South American Shield. This process would use the Conrad discontinuity as a main detachment surface within an incipient NW dipping continental subduction.     

Post-Collisional Transition from Subduction- to Intraplate-type Magmatism in the Westernmost Mediterranean: Evidence for Continental-Edge Delamination of Subcontinental Lithosphere

Post-collisional magmatism in the southern Iberian and northwesternAfrican continental margins contains important clues for theunderstanding of a possible causal connection between movementsin the Earth's upper mantle, the uplift of continental lithosphereand the origin of circum-Mediterranean igneous activity. Systematicgeochemical and geochronological studies (major and trace element,Sr–Nd–Pb-isotope analysis and laser ⁴⁰Ar/³⁹Ar-agedating) on igneous rocks provide constraints for understandingthe post-collisional history of the southern Iberian and northwesternAfrican continental margins. Two groups of magmatic rocks canbe distinguished: (1) an Upper Miocene to Lower Pliocene (8·2–4·8Ma), Si–K-rich group including high-K (calc-alkaline)and shoshonitic series rocks; (2) an Upper Miocene to Pleistocene(6·3–0·65 Ma), Si-poor, Na-rich group includingbasanites and alkali basalts to hawaiites and tephrites. Maficsamples from the Si–K-rich group generally show geochemicalaffinities with volcanic rocks from active subduction zones(e.g. Izu–Bonin and Aeolian island arcs), whereas maficsamples from the Si-poor, Na-rich group are geochemically similarto lavas found in intraplate volcanic settings derived fromsub-lithospheric mantle sources (e.g. Canary Islands). The transitionfrom Si-rich (subduction-related) to Si-poor (intraplate-type)magmatism between 6·3 Ma (first alkali basalt) and 4·8Ma (latest shoshonite) can be observed both on a regional scaleand in individual volcanic systems. Si–K-rich and Si-poorigneous rocks from the continental margins of southern Iberiaand northwestern Africa are, respectively, proposed to havebeen derived from metasomatized subcontinental lithosphere andsub-lithospheric mantle that was contaminated with plume material.A three-dimensional geodynamic model for the westernmost Mediterraneanis presented in which subduction of oceanic lithosphere is inferredto have caused continental-edge delamination of subcontinentallithosphere associated with upwelling of plume-contaminatedsub-lithospheric mantle and lithospheric uplift. This processmay operate worldwide in areas where subduction-related andintraplate-type magmatism are spatially and temporally associated. KEY WORDS: post-collisional magmatism; Mediterranean-style back-arc basins; subduction; delamination; uplift of marine gateways     

An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

长江中下游庐枞火山岩盆地南侧钾质侵入岩带的成因

庐枞火山岩盆地南侧的钾质侵入岩带由正长岩-石英正长岩-正长花岗岩组成,以石英正长岩为主。它们的形成时间介于123"130 Ma之间,峰值约为126 Ma,其中正长岩和石英正长岩的形成时间稍早,而正长花岗岩的形成时间略晚。整个钾质侵入岩带的侵位时间晚于庐枞盆地内的橄榄玄粗质火山作用约4"7 Ma,也是长江中下游地区除最东段的宁镇地区外中生代最晚的岩浆活动产物之一。地球化学上,该钾质侵入岩带以高钾、富碱、富集Rb、Th、U、K等强不相容元素和轻稀土元素、亏损高场强元素Nb、Ta和Ti为特征。它们的母岩浆主要是由富集型上地幔部分熔融形成的,从正长岩经石英正长岩到正长花岗岩的演化主要受矿物的分离结晶作用控制,地壳物质同化所起的影响不大。但与同样来自富集型上地幔部分熔融的庐枞盆地内火山杂岩的母岩浆相比,前者的母岩浆来源深度可能更大些或其中包含了更多来自软流圈地幔的组分。两者的演化路径也完全不同,钾质侵入岩带的母岩浆除经历过高压下的分离结晶作用外,晚期在低压下还经历过长石为主,可能还有黑云母的分离结晶,甚至上地壳物质一定程度的混染作用;而盆地内火山杂岩的母岩浆低压下矿物的分离结晶作用及上地壳物质的混染都不明显。庐枞盆地南缘的富钾侵入岩与盆地内的火山杂岩一样,地球化学上都具有明显的大陆弧的特征,暗示它们的岩浆源区可能形成于俯冲带环境,意味着扬子地块北缘先前(推测为古元古代晚期)曾发生过俯冲作用,上地幔的交代富集可能就与这次的俯冲作用有关。     

早古生代古亚洲洋俯冲记录:来自东天山卡拉塔格高镁安山岩的年代学、地球化学证据

吐哈盆地南缘荒草坡群大柳沟组火山岩出露于大草滩断裂以北的彩霞山、土屋铜矿北和卡拉塔格地区。该组火山岩主要为基性、中酸性熔岩及火山碎屑岩组合。LA-ICP MS锆石U-Pb分析获得安山岩和英安岩的年龄分别为434.8±3.8Ma和438.4±4.9Ma,表明大柳沟组火山岩形成于早志留世。安山岩Si O2含量为54.35%~60.70%;Al2O3含量为11.8%~16.3%,Mg O含量为4.94%~8.27%,Ti O2含量为0.38%~0.52%,Na2O和K2O的含量分别为0.50%~3.83%和0.08%~1.26%,铝饱和指数(A/CNK)为0.77~1.37,低ΣREE(43×10-6~70×10-6),为具有亏损源区特征的高镁安山岩。英安岩较安山岩Si O2含量高(63.74%~75.35%),具有相似的Al2O3(12.3%~16.0%%)和Ti O2(0.32%~0.52%)含量;Na2O和K2O含量分别为2.01%~5.53%和0.16%~2.19%,ΣREE含量较安山岩略高(61×10-6~84×10-6),铝饱和指数(A/CNK)为0.77~5.74。安山岩和英安岩均轻重稀土分异明显,具有中等的Eu负异常;富集大离子亲石元素(LILE),亏损高场强元素(HFSE),具有明显的Nb-Ta-Ti负异常。相似的岩石地球化学特征表明安山岩和英安岩可能为同源岩浆,英安岩是安山质岩浆经过角闪石等矿物结晶分异作用形成的,安山岩的母岩岩浆很可能来自亏损大洋岩石圈地幔。基于该组岩石组合、地球化学特征和区域地质资料,推测研究区高镁安山岩可能形成于洋内弧环境,这表明东天山吐哈盆地南缘卡拉塔格地区在早古生代已经存在古亚洲洋的俯冲。