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

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

The interplay between fracture and flow in the localization of crustal deformation

The Earth\u2019s crust is generally considered to consist of distinct brittle and viscous (or \u201cductile\u201d) rheological layers, corresponding to Navier-Coulomb failure or viscous flow, with a \u201cbrittle-ductile transition\u201d occurring over a specific and relatively limited depth interval. Depending on the assumed geothermal gradient, a compositionally layered crust could have several such brittle-ductile transitions, but the model still implies that large regions of the crust deform exclusively by either brittle fracture or viscous crystal-plastic flow. However, it is becoming increasingly clear from field observation that, in reality, there is an intimate interplay in space and time between precursor heterogeneities (either structural or compositional), brittle fracture, fluid-rock interaction and more distributed \u201cductile flow\u201d. In particular, there are now several well-documented examples of brittle precursors localizing subsequent ductile deformation under high grade metamorphic conditions ranging from upper amphibolite (Mancktelow and Pennacchioni 2005, 2007) to even eclogite facies (Austrheim and Boundy 1994). This interplay between fracture and crystal-plastic creep and/or diffusion occurs over a wide range of scales, from 100\u2019s of kilometres down to individual grains. Localization of strain in the crust can lead to the development of zones of very large relative displacement (such as low-angle thrusts and detachments, and steep strike-slip faults). The mechanics of this localization on a narrow zone and its repeated reactivation can only be considered in terms of a cyclical interaction between fracture, flow, and variation in local pore-fluid pressure. These relatively planar and discrete faults and shear zones are commonly observed to cross-cut layering and foliation at a small angle. Small-scale examples from the field, as well as numerical models, show that viscous localization is strongly controlled by existing compositional and rheological heterogeneity (such as bedding, dykes, veins etc), whereas fractures may crosscut such compositional layering at small angles. This suggests that major crosscutting faults, which may now be dominated by mylonitic fabrics characteristic of crystal-plastic flow (e.g., the Periadriatic Fault in the European Alps), could also have had a large-scale, brittle precursor that controlled subsequent ductile localization. On a smaller scale, flanking structures (Passchier 2001; Grasemann and St\ufcwe 2001) developed around brittle fractures of limited length are particularly clear examples of interacting brittle-ductile deformation, because their geometry can only be explained if discrete slip occurred synchronously with more distributed surrounding ductile flow (Exner et al. 2004). Flanking structures that formed in calcite marbles under amphibolite facies conditions (e.g., on the island of Naxos, Greece) demonstrate that brittle fracturing can play an important role even in weak rocks at high temperatures \u2013 conditions generally taken to imply exclusively ductile or viscous behaviour. Such flanking structures are common in mylonitic shear zones (e.g., in mylonites in the footwall of the major Simplon low-angle normal fault in the central Alps) and demonstrate the delicate balance between fracture and flow in such high strain zones, with switches back and forth varying locally in space and through time. This behaviour is not totally unexpected. The reduction of bulk porosity and permeability in rocks with depth raises the local pore fluid pressure from hydrostatic to near lithostatic (as usually assumed in metamorphic petrology), with the result that rocks are generally critically stressed and close to failure. Only minor local changes in the controlling parameters (strain rate, pore fluid pressure, dynamic or \u201ctectonic\u201d pressure) can cause a switch between fracture and flow. In natural examples, the interplay between fracture and flow is observed in middle to lower crustal rocks irrespective of whether they are weak (\u201cwet\u201d) or strong (\u201cdry\u201d). Excellent examples of interacting fracture and flow from glacier-polished outcrops of granodiorite in the Neves area of the eastern Alps developed under wet conditions, with very common quartz vein development and marked fluid-rock interaction along fractures. The deviatoric stress during both flow and fracture was low (<10 MPa), as demonstrated by little deformed calcite porphyroclasts in quartz mylonites, which did not even significantly twin during crystal plastic flow of the matrix quartz under upper amphibolite facies conditions (Mancktelow and Pennacchioni 2010). In contrast, in dry lower crust, such as from the Mont Mary area of the western Alps, stresses were high (as indicated by very small recrystallized quartz grain sizes; Fitz Gerald et al. 2006) and seismic fracture was associated with pseudotachlyte development. Pseudotachylytes subsequently act as rheologically weak layers that strongly localize ductile shearing under dry upper amphibolite facies conditions (Pennacchioni and Cesare 1997). References: Austrheim, H., Boundy, T.M., 1994. Pseudotachylytes generated during seismic faulting and eclogitization of the deep crust. Science 265, 82-83. Exner, U., Mancktelow, N.S., Grasemann, B., 2004. Progressive development of s-type flanking folds in simple shear. Journal of Structural Geology 26, 2191-2201. Grasemann, B., St\ufcwe, K., 2001. The development of flanking folds during simple shear and their use as kinematic indicators. Journal of Structural Geology 23, 715-724. Fitz Gerald, J.D., Mancktelow, N.S., Pennacchioni, G., Kunze, K., 2006. Ultrafine-grained quartz mylonites from high-grade shear zones: Evidence for strong dry middle to lower crust. Geology 34, 369-372. 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. Mancktelow, N.S., Pennacchioni, G., 2010. Why calcite can be stronger than quartz. Journal of Geophysical Research 115. Passchier, C.W., 2001. Flanking structures. Journal of Structural Geology 23, 951-962. Pennacchioni, G., Cesare, B., 1997. Ductile-brittle transition in pre-Alpine amphibolite facies mylonites during evolution from water-present to water-deficient conditions (Mont Mary Nappe, Italian Western Alps). Journal of Metamorphic Geology 15, 777-791. 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

Small-scale ductile shear zones: Neither extending, nor thickening, nor narrowing

The length, thickness and strain gradients of small-scale (10 123 \u201310 121 m thick) ductile shear zones are pre-determined by the presence of surface precursors (e.g. fractures or compositional layers) and associated \ufb02uid-rock interaction. Fractures, with their surrounding host-rock damage zone and concurrent tectonic underpressure during dilation, provide e\ufb03cient \ufb02uid pathways and networks with di\ufb00usive \ufb02uid-rock interaction at their margins. The \ufb02uid-induced, gradational to zoned compositional haloes that symmetrically surround a fracture control the strain gradients of developing shear zones and result in a diversity of geometric types, including single homogeneous-to-heterogeneous shear zones and paired shear zones. As a consequence, geochemical di\ufb00erences between shear zone and host rock may re\ufb02ect \ufb02uid-rock interaction during the precursor brittle history rather than during slip, especially considering that shear zones with even a small component of stretch may be over-pressured and therefore unable to drain \ufb02uids from the surrounding rocks. Support from \ufb01eld observation for this interpretation is given by (1) the lack of correlation between shear zone thickness and accumulated displacement; and (2) the similar thickness of shear zones and locally unexploited alteration haloes surrounding fracture precursors. Most small-scale shear zones in massive rocks (granitoids) are initially neither thickening nor narrowing with increasing strain. This concept may also apply to more foliated rocks, but does not necessarily hold for larger-scale shear zones