The earthquake cycle in the dry lower continental crust: insights from two deeply exhumed terranes (Musgrave Ranges, Australia and Lofoten, Norway)

This paper discusses the results of field-based geological investigations of exhumed rocks exposed in the Musgrave Ranges (Central Australia) and in Nusfjord (Lofoten, Norway) that preserve evidence for lower continental crustal earthquakes with focal depths of approximately 25–40 km. These studies have established that deformation of the dry lower continental crust is characterized by a cyclic interplay between viscous creep (mylonitization) and brittle, seismic slip associated with the formation of pseudotachylytes (a solidified melt produced during seismic slip along a fault in silicate rocks). Seismic slip triggers rheological weakening and a transition to viscous creep, which may be already active during the immediate post-seismic deformation along faults initially characterized by frictional melting and wall-rock damage. The cyclical interplay between seismic slip and viscous creep implies transient oscillations in stress and strain rate, which are preserved in the shear zone microstructure. In both localities, the spatial distribution of pseudotachylytes is consistent with a local (deep) source for the transient high stresses required to generate earthquakes in the lower crust. This deep source is the result of localized stress amplification in dry and strong materials generated at the contacts with ductile shear zones, producing multiple generations of pseudotachylyte over geological time. This implies that both the short- and the long-term rheological evolution of the dry lower crust typical of continental interiors is controlled by earthquake cycle deformation. This article is part of a discussion meeting issue ‘Understanding earthquakes using the geological record’.</jats:p

Hydrogen and 40Ar/39Ar isotope evidence for multiple and protracted paleofluid flow events within the long‐lived North Anatolian Keirogen (Turkey)

We present a new approach to identifying the source and age of paleofluids associated with low‐temperature deformation in the brittle crust, using hydrogen isotopic compositions (δD) and 40Ar/39Ar geochronology of authigenic illite in clay gouge‐bearing fault zones. The procedure involves grain‐size separation, polytype modeling, and isotopic analysis, creating a mixing line that is used to extrapolate to δD and age of pure authigenic and detrital material. We use this method on samples collected along the surface trace of today's North Anatolian Fault (NAF). δD values of the authigenic illite population, obtained by extrapolation, are −89 ± 3‰, −90 ± 2‰, and −97 ± 2‰ (VSMOW) for samples KSL, RES4‐1, and G1G2, respectively. These correspond to δD fluid values of −62‰ to −85‰ for the temperature range of 125°C ± 25°, indistinguishable from present‐day precipitation values. δD values of the detrital illite population are −45 ± 13‰, −60 ± 6‰, and −64 ± 6‰ for samples KSL, G1G2, and RES4‐1, respectively. Corresponding δD fluid values at 300°C are −26‰ to −45‰ and match values from adjacent metamorphic terranes. Corresponding clay gouge ages are 41.4 ± 3.4 Ma (authigenic) and 95.8 ± 7.7 Ma (detrital) for sample G2 and 24.6 ± 1.6 Ma (authigenic) and 96.5 ± 3.8 Ma (detrital) for sample RES4‐1, demonstrating a long history of meteoric fluid infiltration in the area. We conclude that today's NAF incorporated preexisting, weak clay‐rich rocks that represent earlier mineralizing fluid events. The samples preserve at least three fluid flow pulses since the Eocene and indicate that meteoric fluid has been circulating in the upper crust in the North Anatolian Keirogen since that time.Key Points:Illite preserves the hydrogen isotopic signature and age of paleofluids in the earth's upper crustThree fluid events are pinpointed in the NAKThe NAF exploited zones of preexisting weak clay material during its formationPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/112210/1/ggge20754.pd

Spatial correlation bias in late-Cenozoic erosion histories derived from thermochronology

International audienceThe potential link between erosion rates at the Earth's surface and changes in global climate has intrigued geoscientists for decades1,2 because such a coupling has implications for the influence of silicate weathering3,4 and organic-carbon burial5 on climate and for the role of Quaternary glaciations in landscape evolution1,6. A global increase in late-Cenozoic erosion rates in response to a cooling, more variable climate has been proposed on the basis of worldwide sedimentation rates7. Other studies have indicated, however, that global erosion rates may have remained steady, suggesting that the reported increases in sediment-accumulation rates are due to preservation biases, depositional hiatuses and varying measurement intervals8-10. More recently, a global compilation of thermochronology data has been used to infer a nearly twofold increase in the erosion rate in mountainous landscapes over late-Cenozoic times6. It has been contended that this result is free of the biases that affect sedimentary records11, although others have argued that it contains biases related to how thermochronological data are averaged12 and to erosion hiatuses in glaciated landscapes13. Here we investigate the 30 locations with reported accelerated erosion during the late Cenozoic6. Our analysis shows that in 23 of these locations, the reported increases are a result of a spatial correlation bias—that is, combining data with disparate exhumation histories, thereby converting spatial erosion-rate variations into temporal increases. In four locations, the increases can be explained by changes in tectonic boundary conditions. In three cases, climatically induced accelerations are recorded, driven by localized glacial valley incision. Our findings suggest that thermochronology data currently have insufficient resolution to assess whether late-Cenozoic climate change affected erosion rates on a global scale. We suggest that a synthesis of local findings that include location-specific information may help to further investigate drivers of global erosion rates

Fluid flow during unbending: Implications for slab hydration, intermediate-depth earthquakes and deep fluid subduction

We calculate the tectonic stress profile and associated direction of fluid flow during unbending and dehydration of oceanic plates, for a range of critical parameters that affect their combined elasto-plastic and viscous behaviour, such as bending curvature, age, pore fluid pressure and viscous flow laws. In all models, negative pressure gradients are established at Moho depths, down to the base of the slab elastic core. Fluids released at these depths flow downward across the plate, "wetting" or further hydrating the underlying dry levels, and ultimately accumulate at the base of the elastic core, increasing the pore fluid pressure and triggering deep seismicity. The thickness of the "wet" layer increases for low bending curvatures, low pore fluid pressure, old slabs and dry viscous rheologies. "Wetting" of the upper 10-30 km of the slab has important implications for its rheological and anisotropic structure. Unbending favours the redistribution and trapping of significant amounts of fluids in subducting oceanic plates that will subsequently be released at the base of the upper mantle. (C) 2010 Elsevier B.V. All rights reserved

Initiation and growth of strike-slip faults within intact metagranitoid (Neves area, eastern Alps, Italy)

Exhumed strike-slip faults exposed in the Neves area of the Tauern Window (eastern Alps, Italy) formed in the lower brittle crust under hydrous conditions within intact metagranitoids, where any precursor fractures were already healed due to earlier amphibolite-facies metamorphism. Faults initiated as newly formed en-echelon fractures delineating shear bands, with segmentation occurring over scales of 0.001 to 100 m. Due to the initial en-echelon pattern, stepovers between fault segments were almost invariably contractional during subsequent slip accumulation. In initial low-slip (centimeter to decimeter) stages, synthetic slip on the segments was associated with the development of a set of antithetic faults in the contractional stepovers, oriented at an angle of 30 \u201345 to the bounding faults. Slip on this antithetic set did not fully compensate for the decrease in slip on the main faults, implying an additional deformation mechanism within the stepovers (e.g., block rotation and/or out-of-plane movement). In larger faults with slip on the order of a few meters, the antithetic faults within the stepover were crosscut by synthetic sigmoidal faults that connected the overstepped fault segments and accommodated most of the subsequent displacement transfer. This second stage of evolution involving fault linkage is well documented in the Mesule fault, which has a current maximum offset of ~10 m. In the Mesule fault, the total horizontal slip summed along all secondary faults within the stepover accounts almost entirely for the net slip decrease toward the tips of the overstepping faults. However, the boundary faults remain nearly straight in spite of the sigmoidal ramp-like geometry of the connecting faults. Since the adjacent blocks are little deformed and there is no evidence for appreciable volume loss or block rotation in the stepover, signifi cant out-of- plane movement is implied, although it is difficult to quantify. The Neves area provides unusually detailed field constraints on fault initiation, linkage, and displacement accumulation within nonbedded and relatively isotropic granitoid rocks at the base of the brittle crust, where neither a free upper surface nor substantial volume change (e.g., by veining and pressure solution) was a controlling factor in accommodating fault linkage and displacement transfer

Brittle precursors, fluid\u2010rock interaction and the localization or spreading of shear zones

Heterogeneous ductile shear zones are generally considered to develop by progressive strain localization, implying that shear zones become narrower during their development and that individual zones should affect an ever decreasing volume of rock. We support a diametrically different model: in granitoid plutons and other non-layered rock bodies, shear zones are strongly localized at their very initiation, on pre-existing planar rheological discontinuities, and tend to spread into the adjacent rock with increasing strain. Precursor discontinuities can be either compositional layers (e.g. dykes or veins) or fractures (Fig. 1), with enhanced fluid flow and fluid-rock interaction along these fractures leading to localized compositional and rheological change of the original host rock. Spreading of strain reflects the interplay between two factors: (1) diffusion of fluid away from the central fracture, which broadens the zone of alteration, and (2) development of new fractures both in previously intact rock and in already sheared domains. Cycles of fracturing are driven by local stress concentrations in rocks that remain close to the critical stress state for fracture. Stress concentration can be due to local mechanical instability (e.g. dyke boudinage) or more generally due to inherent problems of strain accommodation during deformation of more \u201crigid\u201d blocks surrounded by a network of relatively discrete shear zones. As is clear from analogue and numerical models, such blocks between bounding shear zones must also deform internally to maintain strain compatibility. A distributed, more homogeneous background strain may develop in the intervening blocks under higher grade metamorphic conditions, but our field observations demonstrate that more localized shearing of intact granitic protolith in general develops from a brittle precursor. The conference location is particularly appropriate for this topic, as one of the first studies that proposed a brittle precursor to heterogeneous ductile shear zone development was that of Segall and Simpson (1986), who used the Roses shear zones as relevant examples. There are many published studies of shear zone development in previously undeformed granitoid plutons, especially with regard to gradients in microstructure and chemical and/or isotopic changes during \u201cshear localization\u201d. This interest was in part based on the assumption that plutons are relatively homogeneous, allowing a direct comparison between the heterogeneous shear zone and the homogeneous background. However, such intrusive bodies are certainly not homogeneous in detail: they show common compositional boundaries due to enclaves, intrusive contacts and dykes and veins. Cooling plutons also invariably develop a pervasive set of joints due to thermal contraction, typically involving a volume decrease on the order of 15% or more. These joints form conduits for late magmatic or subsequent metamorphic fluids, with the development of veins (especially quartz veins) and localized new mineral growth (commonly biotite in higher temperature cooling joints). These precursor discontinuities act as the controlling loci for localizing shear zones either during pluton cooling (Pennacchioni 2005; Pennacchioni et al. 2010) or later deformation (Mancktelow and Pennacchioni 2005; Pennacchioni and Mancktelow 2007). However, the process is not necessarily limited to precursor magmatic structures (dykes, veins or cooling joints). Deformation under higher grade metamorphic conditions (e.g. upper amphibolite facies, as in Fig. 2) can produce repeated cycles of fracture, fluid rock interaction and ductile shear localized on these brittle precursors. Already developed broader shear zones can themselves be cut by discrete fractures, with fluid-rock interaction and new mineral growth once more changing the local rheology and localizing further shearing to produce new heterogeneous shear zones oblique to the earlier zone (Fig. 2). Brittle precursors localizing ductile shearing have also been proposed for greenschist facies shear zones developed in schists of the Cap de Creus (Fusseis et al. 2006). References Fusseis, F., Handy, M.R., Schrank, C., 2006. Networking of shear zones at the brittle-to-viscous transition (Cap de Creus, NE Spain). Journal of Structural Geology 28, 1228-1243. Mancktelow, N.S., Pennacchioni, G., 2005. The control of precursor brittle fracture and fluid-rock interaction on the development of single and paired ductile shear zones. Journal of Structural Geology 27, 645-661. Pennacchioni, G., 2005. Control of the geometry of precursor brittle structures on the type of ductile shear zone in the Adamello tonalites, Southern Alps (Italy). Journal of Structural Geology 27, 627-644. Pennacchioni, G., Mancktelow, N.S., 2007. Nucleation and initial growth of a shear zone network within compositionally and structurally heterogeneous granitoids under amphibolite facies conditions. Journal of Structural Geology 29, 1757-1780. Pennacchioni, G., Menegon, L., Leiss, B., Nestola, F., Bromiley, G., 2010. Development of crystallographic preferred orientation and microstructure during plastic deformation of natural coarsegrained quartz veins. Journal of Geophysical Research-Solid Earth 115. Segall, P., Simpson, C., 1986. Nucleation of ductile shear zones on dilatant fractures. Geology 14, 56-59