Contribution of the 2010 Maule Megathrust Earthquake to the Heat Flow at the Peru-Chile Trench

The 2010 Maule earthquake was a megathrust event that occurred along the Peru–Chile Trench. The earthquake source can be modelled as a fault with two asperities with different areas and strengths. By employing a discrete fault model, where asperities are the basic elements, the event can be described as a sequence of three dynamic modes involving simultaneous asperity slip. Interaction between asperities by mutual stress transfer plays a crucial role during fault slip. With a careful choice of values for the model parameters, the mode durations, the slip distribution, the seismic moment rate and the final moment calculated from the model are found to be consistent with the observed values. An important amount of frictional heat is produced by an event of this size and is calculated by summing up the contributions of each asperity. The seismic event produces a heat pulse propagating through the Earth’s crust and contributing to the average heat flow in the region. The calculated heat production is equal to about 2×1017 J and the peak value of the heat pulse is equal to 6×10−3 mW m−2 or about 10−4 of the average surface heat flow density, with a characteristic diffusion time in the order of 106 a.

Structural and thermal evolution of the eastern Aar Massif: insights from structural field work and Raman thermometry

The thermo-kinematic evolution of the eastern Aar Massif, Swiss Alps, was investigated using peak temperature data estimated from Raman spectroscopy of carbonaceous material and detailed field analyses. New and compiled temperature-time constraints along the deformed and exhumed basement-cover contact allow us to (i) establish the timing of metamorphism and deformation, (ii) track long-term horizontal and vertical orogenic movements and (iii) assess the influence of temperature and structural inheritance on the kinematic evolution. We present a new shear zone map, structural cross sections and a step-wise retrodeformation. From ca.\;26\,Ma onwards, basement-involved deformation started with the formation of relatively discrete NNW-directed thrusts. Peak metamorphic isograds are weakly deformed by these thrusts, suggesting that they initiated before or during the metamorphic peak under ongoing burial in the footwall to the basal Helvetic roof thrust. Subsequent peak- to post-metamorphic deformation was dominated by steep, mostly NNW-vergent reverse faults ( ca.  22-14 Ma). Field investigations demonstrate that these shear zones were steeper than 50 ∘ already at inception. This produced the massif-internal structural relief and was associated with large vertical displacements (7 km shortening vs. up to 11 km exhumation). From 14 Ma onwards, the eastern Aar massif exhumed "en bloc" (i.e., without significant differential massif-internal exhumation) in the hanging wall of frontal thrusts, which is consistent with the transition to strike-slip dominated deformation observed within the massif. Our results indicate 13 km shortening and 9 km exhumation between 14 Ma and present. Inherited normal faults were not significantly reactivated. Instead, new thrusts/reverse faults developed in the basement below syn-rift basins, and can be traced into overturned fold limbs in the overlying sediment, producing tight synclines and broad anticlines along the basement-cover contact. The sediments were not detached from their crystalline substratum and formed disharmonic folds. Our results highlight decreasing rheological contrasts between (i) relatively strong basement and (ii) relatively weak cover units and inherited faults at higher temperature conditions. Both the timing of basement-involved deformation and the structural style (shear zone dip) appear to be controlled by evolving temperature conditions.

Mineralization Epochs of Granitic Rare-Metal Pegmatite Deposits in the Songpan–Ganzê Orogenic Belt and Their Implications for Orogeny

Granitic pegmatite deposits, which are usually products of orogenic processes during plate convergence, can be used to demonstrate regional tectonic evolution processes. In the eastern Tibetan Plateau in China, the Jiajika, Dahongliutan, Xuebaoding, Zhawulong, and Ke’eryin rare metal pegmatite deposits are located in the southern, western, northern, midwestern, and central areas of the Songpan–Ganzê orogenic belt, respectively. In this study, we dated two muscovite Ar–Ar ages of 189.4 ± 1.1 Ma and 187.0 ± 1.1 Ma from spodumene pegmatites of the Dahongliutan deposit. We also dated one zircon U-Pb age of 211.6 ± 5.2 Ma from muscovite granite, two muscovite Ar–Ar ages of 179.6 ± 1.0 Ma and 174.3 ± 0.9 Ma, and one columbite–tantalite U-Pb age of 204.5 ± 1.8 Ma from spodumene pegmatites of the Zhawulong deposit. In addition, we dated one muscovite Ar–Ar age of 159.0 ± 1.4 Ma from spodumene pegmatite of the Ke’eryin deposit. Combining these ages and previous studies in chronology, we concluded that the granitic magma in the Jiajika, Xuebaoding, Dahongliutan, Zhawulong, and Ke’eryin deposits intruded into Triassic metaturbidites at approximately 223, 221, 220–217, 212, and 207–205 Ma, respectively, and that the crystallization of the corresponding pegmatite ceased at approximately 199–196, 195–190, 189–187, 180–174, and 159 Ma, respectively. In this study, we demonstrated that the peak in magmatic activity and the final crystallization age of the pegmatite lagged behind one another from the outer areas of the orogeny belt to the inner areas. The pegmatite–parented granitic magmas were sourced from Triassic metaturbidites that were melted by shear heating along the large-scale decollement resulting from Indosinian collisions along the North China block, Qiangtang–Changdu block, and Yangtze block. As a result, the above temporal and spatial regularities indicated that the tectonic–thermal stress resulting from the collisions of three blocks was transferred from the outer areas of the orogenic belt to the inner areas. A large amount of heat and a slow cooling rate at the convergent center of thermal stress in two directions will lead to crystallization and differentiation of magma in the Songpan–Ganzê orogenic belt, forming additional rare metal deposits.

Shear heating reconciles thermal models with the metamorphic rock record of subduction

Some commonly referenced thermal-mechanical models of current subduction zones imply temperatures that are 100-500 °C colder at 30-80-km depth than pressure-temperature conditions determined thermobarometrically from exhumed metamorphic rocks. Accurately inferring subduction zone thermal structure, whether from models or rocks, is crucial for predicting metamorphic reactions and associated fluid release, subarc melting conditions, rheologies, and fault-slip phenomena. Here, we compile surface heat flow data from subduction zones worldwide and show that values are higher than can be explained for a frictionless subduction interface often assumed for modeling. An additional heat source--likely shear heating--is required to explain these forearc heat flow values. A friction coefficient of at least 0.03 and possibly as high as 0.1 in some cases explains these data, and we recommend a provisional average value of 0.05 ± 0.015 for modeling. Even small coefficients of friction can contribute several hundred degrees of heating at depths of 30-80 km. Adding such shear stresses to thermal models quantitatively reproduces the pressure-temperature conditions recorded by exhumed metamorphic rocks. Comparatively higher temperatures generally drive rock dehydration and densification, so, at a given depth, hotter rocks are denser than colder rocks, and harder to exhume through buoyancy mechanisms. Consequently--conversely to previous proposals--exhumed metamorphic rocks might overrepresent old-cold subduction where rocks at the slab interface are wetter and more buoyant than in young-hot subduction zones.

Kinematic variables and water transport control the formation and location of arc volcanoes

The processes that give rise to arc magmas at convergent plate margins have long been a subject of scientific research and debate. A consensus has developed that the mantle wedge overlying the subducting slab and fluids and/or melts from the subducting slab itself are involved in the melting process. However, the role of kinematic variables such as slab dip and convergence rate in the formation of arc magmas is still unclear. The depth to the top of the subducting slab beneath volcanic arcs, usually approximately 110 +/- 20 km, was previously thought to be constant among arcs. Recent studies revealed that the depth of intermediate-depth earthquakes underneath volcanic arcs, presumably marking the slab-wedge interface, varies systematically between approximately 60 and 173 km and correlates with slab dip and convergence rate. Water-rich magmas (over 4-6 wt% H(2)O) are found in subduction zones with very different subduction parameters, including those with a shallow-dipping slab (north Japan), or steeply dipping slab (Marianas). Here we propose a simple model to address how kinematic parameters of plate subduction relate to the location of mantle melting at subduction zones. We demonstrate that the location of arc volcanoes is controlled by a combination of conditions: melting in the wedge is induced at the overlap of regions in the wedge that are hotter than the melting curve (solidus) of vapour-saturated peridotite and regions where hydrous minerals both in the wedge and in the subducting slab break down. These two limits for melt generation, when combined with the kinematic parameters of slab dip and convergence rate, provide independent constraints on the thermal structure of the wedge and accurately predict the location of mantle wedge melting and the position of arc volcanoes.

Cenozoic climate change as a possible cause for the rise of the Andes

Causal links between the rise of a large mountain range and climate have often been considered to work in one direction, with significant uplift provoking climate change. Here we propose a mechanism by which Cenozoic climate change could have caused the rise of the Andes. Based on considerations of the force balance in the South American lithosphere, we suggest that the height of, and tectonics in, the Andes are strongly controlled both by shear stresses along the plate interface in the subduction zone and by buoyancy stress contrasts between the trench and highlands, and shear stresses in the subduction zone depend on the amount of subducted sediments. We propose that the dynamics of subduction and mountain-building in this region are controlled by the processes of erosion and sediment deposition, and ultimately climate. In central South America, climate-controlled sediment starvation would then cause high shear stress, focusing the plate boundary stresses that support the high Andes.

The Interdependence of Deformational and Thermal Processes in Mountain Belts

Crustal temperatures within collisional orogens are anomalously high compared with temperatures at comparable depths in stable continents, which is evidence of thermal processes that are fundamental to orogenesis. These temperatures can be explained by the redistribution of crust enriched in heat-producing elements through the accretion of crust from the down-going plate to the upper plate and surface erosion. With the use of geologically reasonable rates, the model results predict high temperatures (over 600°C) and inverted upper-plate geotherms (about 100°C over 20 kilometers) at shallow depths (20 to 40 kilometers) by 25 to 35 million years after collision. This study emphasizes the interdependence of deformational, surficial, and thermal processes.