Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Stratigraphy, facies architecture and tectonic implications of the Upper Devonian to Lower Carboniferous Campwyn Volcanics of the northern New England Fold Belt

Upper Devonian to Lower Carboniferous strata of the Campwyn Volcanics of east central Queensland preserve a substantial sequence of first-cycle volcaniclastic sedimentary and coeval volcanic rocks that record prolonged volcanic activity along the northern New England Fold Belt. The style and scale of volcanism varied with time, producing an Upper Devonian sequence of mafic volcano-sedimentary rocks overlain by a rhyolitic ignimbrite-dominated sequence that passes upward into a Lower Carboniferous limestone-bearing sedimentary sequence. We define two facies associations for the Campwyn Volcanics. A lower facies association is dominated by mafic volcanic-derived sedimentary breccias with subordinate primary mafic volcanic rocks comprising predominantly hyaloclastite and peperite. Sedimentary breccias record episodic and high energy, subaqueous depositional events with clastic material sourced from a mafic lava-dominated terrain. Some breccias contain a high proportion of attenuated dense, glassy mafic juvenile clasts, suggesting a syn-eruptive origin. The lower facies association coarsens upwards from a lithic sand-dominated sequence through a thick interval of pebble- to boulder-grade polymict volcaniclastic breccias, culminating in facies that demonstrate subaerial exposure. The silicic upper facies association marks a significant change in eruptive style, magma composition and the nature of eruptive sources, as well as the widespread development of subaerial depositional conditions. Crystal-rich, high-grade, low- to high-silica rhyolite ignimbrites dominate the base of this facies association. Biostratigraphic age controls indicate that the ignimbrite-bearing sequences are Famennian to lower-mid Tournaisian in age. The ignimbrites represent extra-caldera facies with individual units up to 40 m thick and mostly lacking coarse lithic breccias. Thick deposits of pyroclastic material interbedded with fine-grained siliceous sandstone and mudstone (locally radiolarian-bearing) were deposited from pyroclastic flows that crossed palaeoshorelines or represent syn-eruptive, resedimented pyroclastic material. Some block-bearing lithic-pumice-crystal breccias may also reflect more proximal subaqueous silicic explosive eruptions. Crystal-lithic sandstones interbedded with, and overlying the ignimbrites, contain abundant detrital volcanic quartz and feldspar derived from the pyroclastic deposits. Limestone is common in the upper part of the upper facies association, and several beds are oolitic (cf. Rockhampton Group of the Yarrol terrane). Overall, the upper facies association fines upward and is transgressive, recording a return to shallow-marine conditions. Palaeocurrent data from all stratigraphic levels in the Campwyn Volcanics indicate that the regional sediment-dispersal direction was to the northwest, and opposed to the generally accepted notion of easterly sediment dispersal from a volcanic arc source. The silicic upper facies association correlates in age and lithology to Early Carboniferous silicic volcanism in the Drummond (Cycle 1) and Burdekin Basins, Connors Arch, and in the Yarrol terranes of eastern Queensland. The widespread development of silicic volcanism in the Early Carboniferous indicates that silicic (rift-related) magmatism was not restricted to the Drummond Basin, but was part of a more substantial silicic igneous province.

Microstructural evolution of the Serido Belt, NE-Brazil: the effect of two tectonic events on development of c-axis preferred orientation in quartz

The polyphase evolution of the Serido Belt (NE-Brazil) includes D, crust formation at 2.3-2.1 Ga, D-2 thrust tectonics at 1.9 Ga and crustal reworking by D-3 strike-slip shear zones at 600 Ma. Microstructural investigations within mylonites associated with D-2 and D-3 events were used to constrain the tectono-thermal evolution of the belt. D-2 shear zones commenced at deeper crustal levels and high amphibolite facies conditions (600-650 degreesC) through grain boundary migration, subgrain rotation and operation of quartz Q-prism slip. Continued shearing and exhumation of the terrain forced the re-equilibration of high-T fabrics and the switching of slip systems from (c)-prism to positive and negative (a)-rhombs. During D-3, enhancement of ductility by dissipation of heat that came from syn-D-3 granites developed wide belts of amphibolite facies mylonites. Continued shearing, uplift and cooling of the region induced D-3 shear zones to act in ductile-brittle regimes, marked by fracturing and development of thinner belts of greenschist facies mylonites. During this event, switching from (a)-prism to a basal slip indicates a thermal path from 600 to 350 degreesC. Therefore, microstructures and quartz c-axis fabrics in polydeformed rocks from the Serido Belt preserve the record of two major events, which includes contrasting deformation mechanisms and thermal paths. (C) 2003 Elsevier Ltd. All rights reserved.

Reconstructing the architectural evolution of volcanic islands from combined K/AR, morphologic, tectonic, and magnetic data. The Faial Island exemple (Azores)

The morpho-structural evolution of oceanic islands results from competition between volcano growth and partial destruction by mass-wasting processes. We present here a multi-disciplinary study of the successive stages of development of Faial (Azores) during the last 1 Myr. Using high-resolution digital elevation model (DEM), and new K/Ar, tectonic, and magnetic data, we reconstruct the rapidly evolving topography at successive stages, in response to complex interactions between volcanic construction and mass wasting, including the development of a graben. We show that: (1) sub-aerial evolution of the island first involved the rapid growth of a large elongated volcano at ca. 0.85 Ma, followed by its partial destruction over half a million years; (2) beginning about 360 ka a new small edifice grew on the NE of the island, and was subsequently cut by normal faults responsible for initiation of the graben; (3) after an apparent pause of ca. 250 kyr, the large Central Volcano (CV) developed on the western side of the island at ca 120 ka, accumulating a thick pile of lava flows in less than 20 kyr, which were partly channelized within the graben; (4) the period between 120 ka and 40 ka is marked by widespread deformation at the island scale, including westward propagation of faulting and associated erosion of the graben walls, which produced sedimentary deposits; subsequent growth of the CV at 40 ka was then constrained within the graben, with lava flowing onto the sediments up to the eastern shore; (5) the island evolution during the Holocene involves basaltic volcanic activity along the main southern faults and pyroclastic eruptions associated with the formation of a caldera volcano-tectonic depression. We conclude that the whole evolution of Faial Island has been characterized by successive short volcanic pulses probably controlled by brief episodes of regional deformation. Each pulse has been separated by considerable periods of volcanic inactivity during which the Faial graben gradually developed. We propose that the volume loss associated with sudden magma extraction from a shallow reservoir in different episodes triggered incremental downward graben movement, as observed historically, when immediate vertical collapse of up to 2 m was observed along the western segments of the graben at the end of the Capelinhos eruptive crises (1957-58).

Trypanites ichnofacies: palaeoenvironmental and tectonic implications. A case study from the miocene disconformity at Foz da Fonte (lower Tagus basin, Portugal)

Palaeogeography, Palaeoclimatology, Palaeoecology 292, 35–43

Probabilistic tsunami hazard in the Northeast Atlantic from near- and far-field tectonic sources

In this article, we present the first study on probabilistic tsunami hazard assessment for the Northeast (NE) Atlantic region related to earthquake sources. The methodology combines the probabilistic seismic hazard assessment, tsunami numerical modeling, and statistical approaches. We consider three main tsunamigenic areas, namely the Southwest Iberian Margin, the Gloria, and the Caribbean. For each tsunamigenic zone, we derive the annual recurrence rate for each magnitude range, from Mw 8.0 up to Mw 9.0, with a regular interval, using the Bayesian method, which incorporates seismic information from historical and instrumental catalogs. A numerical code, solving the shallow water equations, is employed to simulate the tsunami propagation and compute near shore wave heights. The probability of exceeding a specific tsunami hazard level during a given time period is calculated using the Poisson distribution. The results are presented in terms of the probability of exceedance of a given tsunami amplitude for 100- and 500-year return periods. The hazard level varies along the NE Atlantic coast, being maximum along the northern segment of the Morocco Atlantic coast, the southern Portuguese coast, and the Spanish coast of the Gulf of Cadiz. We find that the probability that a maximum wave height exceeds 1 m somewhere in the NE Atlantic region reaches 60 and 100 % for 100- and 500-year return periods, respectively. These probability values decrease, respectively, to about 15 and 50 % when considering the exceedance threshold of 5 m for the same return periods of 100 and 500 years.

On the applicability of joint inversion of gravity and resistivity data to the study of a tectonic sedimentary basin in Northern Portugal

The Chaves basin is a pull-apart tectonic depression implanted on granites, schists, and graywackes, and filled with a sedimentary sequence of variable thickness. It is a rather complex structure, as it includes an intricate network of faults and hydrogeological systems. The topography of the basement of the Chaves basin still remains unclear, as no drill hole has ever intersected the bottom of the sediments, and resistivity surveys suffer from severe equivalence issues resulting from the geological setting. In this work, a joint inversion approach of 1D resistivity and gravity data designed for layered environments is used to combine the consistent spatial distribution of the gravity data with the depth sensitivity of the resistivity data. A comparison between the results from the inversion of each data set individually and the results from the joint inversion show that although the joint inversion has more difficulty adjusting to the observed data, it provides more realistic and geologically meaningful models than the ones calculated by the inversion of each data set individually. This work provides a contribution for a better understanding of the Chaves basin, while using the opportunity to study further both the advantages and difficulties comprising the application of the method of joint inversion of gravity and resistivity data.

Depleted arc volcanism in the Alboran Sea and shoshonitic volcanism in Morocco: geochemical and isotopic constraints on Neogene tectonic processes

Samples of volcanic rocks from Alboran Island, the Alboran Sea floor and from the Gourougou volcanic centre in northern Morocco have been analyzed for major and trace elements and Sr-Nd isotopes to test current theories on the tectonic geodynamic evolution of the Alboran Sea. The Alboran Island samples are low-K tholeiitic basaltic andesites whose depleted contents of HFS elements (similar to0.5xN-MORB), especially Nb (similar to0.2xN-MORB), show marked geochemical parallels with volcanics from immature intra-oceanic arcs and back-arc basins. Several of the submarine samples have similar compositions, one showing low-Ca boninite affinity. Nd-143/Nd-144 ratios fall in the same range as many island-arc and back-arc basin samples, whereas Sr-87/Sr-86 ratios (on leached samples) are somewhat more radiogenic. Our data point to active subduction taking place beneath the Alboran region in Miocene times, and imply the presence of an associated back-arc spreading centre. Our sea floor suite includes a few more evolved dacite and rhyolite samples with (Sr-87/Sr-86)(0) up to 0.717 that probably represent varying degrees of crustal melting. The shoshonite and high-K basaltic andesite lavas from Gourougou have comparable normalized incompatible-element enrichment diagrams and Ce/Y ratios to shoshonitic volcanics from oceanic island arcs, though they have less pronounced Nb deficits. They are much less LIL- and LREE-enriched than continental arc analogues and post-collisional shoshonites from Tibet. The magmas probably originated by melting in subcontinental lithospheric mantle that had experienced negligible subduction input. Sr-Nd isotope compositions point to significant crustal contamination which appears to account for the small Nb anomalies. The unmistakable supra-subduction zone (SSZ) signature shown by our Alboran basalts and basaltic andesite samples refutes geodynamic models that attribute all Neogene volcanism in the Alboran domain to decompression melting of upwelling asthenosphere arising from convective thinning of over-thickened lithosphere. Our data support recent models in which subsidence is caused by westward rollback of an eastward-dipping subduction zone beneath the westemmost Mediterranean. Moreover, severance of the lithosphere at the edges of the rolling-back slab provides opportunities for locally melting lithospheric mantle, providing a possible explanation for the shoshonitic volcanism seen in northern Morocco and more sporadically in SE Spain. (C) 2004 Elsevier B.V. All rights reserved.

Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

40Ar/39Ar dating of Penninic Front tectonic displacement (W Alps) during the Lower Oligocene (31-34 Ma)

Direct absolute dating of the Penninic Frontal Thrust tectonic motion is achieved using the Ar-40/Ar-39 technique in the Pelvoux Crystalline Massif (Western Alps). The dated phengites were formed syn-kinematically in shear zones. They underline the brittle-ductile stretching lineation, pressure-shadow fibres and slickensides consistent with underthrusting of the European continental slab below the propagating Penninic Thrust. Chlorite-phengite thermobarometry yields 10-15 km and T similar to 280 degrees C, while Ar-40/Ar-39 phengite ages mainly range between 34 and 30 Ma, with one younger age at 27 Ma. This Early Oligocene age range matches a major tectonic rearrangement of the Alpine chain. Preservation of prograde Ar-40/Ar-39 ages is ascribed to passive exhumation of the Pelvoux shear zone network, sandwiched between more external thrusts and the Penninic Front reactivated as an E-dipping detachment fault. Partial resetting in the Low Temperature part of argon spectra below 24 Ma is ascribed to brittle deformation and alteration of phengites.

An alternative plate tectonic model for the Palaeozoic-Early Mesozoic Palaeotethyan evolution of Southeast Asia (northern Thailand-Myanmar)

An alternative model for the geodynamic evolution of Southeast Asia is proposed and inserted in a modern plate tectonic model. The reconstruction methodology is based on dynamic plate boundaries, constrained by data such as spreading rates and subduction velocities; in this way it differs from classical continental drift models proposed so far. The different interpretations about the location of the Palaeotethys suture in Thailand are revised, the Tertiary Mae Yuam fault is seen as the emplacement of the suture. East of the suture we identify an Indochina derived terrane for which we keep the name Shan-Thai, formerly used to identify the Cimmerian block present in Southeast Asia, now called Sibumasu. This nomenclatural choice was made on the basis of the geographic location of the terrane (Eastern Shan States in Burma and Central Thailand) and in order not to introduce new confusing terminology. The closure of the Eastern Palaeotethys is related to a southward subduction of the ocean, that triggered the Eastern Neotethys to open as a back-arc, due to the presence of Late Carboniferous-Early Permian arc magmatism in Mergui (Burma) and in the Lhasa block (South Tibet), and to the absence of arc magmatism of the same age East of the suture. In order to explain the presence of Carboniferous-Early Permian and Permo-Triassic volcanic arcs in Cambodia, Upper Triassic magmatism in Eastern Vietnam and Lower Permian-Middle Permian arc volcanites in Western Sumatra, we introduce the Orang Laut terranes concept. These terranes were detached from Indochina and South China during back-arc opening of the Poko-Song Ma system, due to the westward subduction of the Palaeopacific. This also explains the location of the Cathaysian West Sumatra block to the West of the Cimmerian Sibumasu block.

Tectonic evolution of the High Himalaya in upper Lahul (NW Himalaya, India)

The Upper Lahul region in the NW Himalaya is located in the transition zone between the High Himalayan Crystalline (HHC) to the SW and the Tethyan Zone sedimentary series to the NE. The tectonic evolution of these domains during the Himalayan Orogeny is the consequence of a succession of five deformation events. An early D1 phase corresponds to synmetamorphic, NE verging folding. This deformation created the Tandi Syncline, which consists of Permian to Jurassic Tethyan metasediments cropping out in the core of a large-scale synformal fold within the HHC paragneiss. This tectonic event is interpreted as related to a NE directed nappe stacking (Shikar Beh Nappe), probably during the late Eocene to the early Oligocene. A subsequent D2a phase caused SW verging folding in the HHC. This deformation is interpreted as contemporaneous with late Oligocene to early Miocene SW directed thrusting along the Main Central Thrust. In the Tethyan Zone, a D2b phase is marked by a decollement thrust, a system of reverse faults, and gentle folds, associated with SW directed tectonic movements. This deformation is related to an imbricate structure, characteristic of a shallow structural level, and developed in the frontal part of a nappe affecting the Tethyan Zone units of SE Zanskar (Nyimaling-Tsarap Nappe). A later D3 phase generated the Chandra Dextral Shear Zone (CDSZ), a large-scale, ductile, dextral strike-slip shear zone, located in the transition zone between the HHC and the Tethyan Himalaya. The CDSZ most likely represents a part of a system of early Miocene extensional and/or dextral, strike-slip shear zones-observed at the HHC-Tethyan Zone contact along the entire Himalaya. A final D4 phase induced large-scale doming and NE:verging back folding.