Strain localization across the brittle-ductile transition is a fundamental process in accommodating tectonic movement in the mid-crust. The tectonically active Cordillera Blanca shear zone (CBSZ), a ∼200-km-long normal-sense shear zone situated within the footwall of a discrete syn-convergent extensional fault in the Peruvian Andes, is an excellent field laboratory to explore this transition. Field and microscopic observations indicate consistent top-down-to-the-southwest sense of shear and a sequence of tectonites ranging from undeformed granodiorite through mylonite and ultimately fault breccia along the detachment.

Using microstructural analysis, two-feldspar and Ti-in-quartz (TitaniQ) thermometry, recrystallized quartz paleopiezometry, and analysis of quartz crystallographic preferred orientations, we evaluate the deformation conditions and mechanisms in quartz and feldspar across the CBSZ. Deformation temperatures derived from asymmetric strain-induced myrmekite in a subset of tectonite samples are 410 ± 30 to 470 ± 36 °C, consistent with TitaniQ temperatures of 450 ± 60 to 490 ± 33 °C and temperatures >400 °C estimated from microstructural criteria. Brittle fabrics overprint ductile fabrics within ∼150 m of the detachment that indicate that deformation continued to lower-temperature (∼280–400 °C) and/or higher-strain-rate conditions prior to the onset of pervasive brittle deformation. Initial deformation occurred via high-temperature fracturing and dissolution-precipitation in feldspar. Continued subsolidus deformation resulted in either layering of mylonites into monophase quartz and fine-grained polyphase domains oriented subparallel to macroscopic foliation or the interconnection of recrystallized quartz networks oriented obliquely to macroscopic foliation. The transition to quartz-controlled rheology occurred at temperatures near ∼500 °C and at a differential stress of ∼16.5 MPa. Deformation within the CBSZ occurred predominantly above ∼400 °C and at stresses up to ∼71.4 MPa prior to the onset of brittle deformation.


One of the most important long-standing problems in tectonics is the localization of strain during the transition from ductile to brittle deformation in the lithosphere (Kirby, 1983; Carter and Tsenn, 1987; Rutter et al., 2001; Platt and Behr, 2011a; Sullivan et al., 2013; Spruzeniece and Piazolo, 2015). Deformation mechanisms and the relative roles of temperature, strain rate, and fluids in localizing strain can vary widely along the brittle-ductile transition (BDT). Crustal-scale shear zones are ideal settings for quantifying deformation conditions and discerning the mechanisms facilitating strain localization along the BDT. However, many such examples have protracted and complicated histories, making it difficult to discriminate the relative timing of ductile and brittle behavior (Singleton and Mosher, 2012; Cooper et al., 2017). As crustal strength along the BDT is contingent upon lithology and tectonic regime (Sibson, 1983; Carter and Tsenn, 1987), it is important to describe naturally occurring deformation in a range of geologic environments. Although previous work on crustal-scale shear zones near the BDT in granitoids is extensive (Simpson, 1985; Fitz Gerald and Stünitz, 1993; Stünitz and Fitz Gerald, 1993; Tsurumi et al., 2002; Ree et al., 2005; Pennacchioni et al., 2006; Behr and Platt, 2011; Zibra et al., 2012; Singleton and Mosher, 2012; Wintsch and Yeh, 2013; Sullivan et al., 2013; Columbu et al., 2015; Spruzeniece and Piazolo, 2015; Wehrens et al., 2016), studies of young and/or currently active extensional shear zones are less common (Bessiere et al., 2018).

Rock strength along the BDT can be affected by variables that influence brittle deformation, including fluid pressure, strain rate, effective pressure, and mineralization along failure surfaces (Byerlee, 1978; Brace and Kohlstedt, 1980; Carter and Tsenn, 1987), and those that influence ductile deformation, such as temperature, activity of water (Kirby, 1980; Tullis and Yund, 1980; Yund and Tullis, 1980), interconnectedness of minerals (Wintsch and Yeh, 2013), differential stress, lithology, deformation mechanism, and preferred mineral orientations (Rutter, 1986; Carter and Tsenn, 1987; Behr and Platt, 2011; Platt and Behr, 2011b; Cooper et al., 2017). Experimental studies have constrained many of these factors (Tullis and Yund, 1991, 1987, 1980; Hirth and Tullis, 1992, 1994; Heilbronner and Tullis, 2006; Stünitz and Tullis, 2001; Holyoke and Tullis, 2006; Holyoke and Kronenberg, 2010; Stipp and Tullis, 2003), providing a framework for studying “natural laboratories” (e.g., Lloyd et al., 1992; Hirth et al., 2001; Stipp et al., 2002a, 2002b; Wells et al., 2005; Behr and Platt, 2011, 2014; Gottardi and Teyssier, 2013; Sullivan et al., 2013; Spruzeniece and Piazolo, 2015; Viegas et al., 2016; Rahl and Skemer, 2016).

The tectonically active Cordillera Blanca shear zone (CBSZ), an ∼200-km-long normal-sense shear zone in the Peruvian Andes, occurs almost entirely in the granodioritic Cordillera Blanca batholith (Atherton and Sanderson, 1987; Petford and Atherton, 1992; McNulty et al., 1998), which was emplaced at 14–5 Ma (U-Pb zircon; Mukasa, 1984; McNulty et al., 1998; Giovanni, 2007). This granodioritic intrusion is the predominant lithology in the exhumed footwall of the Cordillera Blanca detachment normal fault, which is the expression of the upper, brittle portion of the CBSZ and has been active since ca. 5 Ma (Bonnot, 1984; Giovanni, 2007). Tectonic fabrics progress from undeformed granodiorite in the core of the batholith in the east, through mylonites and ultramylonites toward the west, culminating in cataclasite and breccia along the detachment surface. The CBSZ is an excellent case study of the BDT, as the relatively simple deformation history and consistent lithology facilitate comparison of fabrics, deformation mechanisms, and deformation conditions at different structural levels across the BDT.

We present new data characterizing mid-crustal deformation within the CBSZ including deformation temperatures of quartz and feldspar, quartz crystallographic preferred orientations (CPOs), flow stress, and microstructural relationships relating to changing deformation mechanisms. Previous work describes the transition from magmatic to subsolidus tectonic fabrics in the CBSZ (Petford and Atherton, 1992), but no constraints on ductile deformation conditions are presently available. An additional motivation for this study is to offer the first quantitative constraints on ductile to brittle deformation conditions within the CBSZ and evaluate the deformation mechanisms that facilitate the fabric transition described by Petford and Atherton (1992). In this study, we present data supporting early deformation by high-temperature fracturing in a feldspar-controlled framework followed by predominantly quartz-controlled rheology, with at least local polyphase-controlled rheology within layered mylonites. Our results indicate that deformation during exhumation produced overprinting fabrics reflecting changes in deformation mechanisms and an increase in strain rate. Strain localized from an ∼450 m thick plastic shear zone to an ∼150 m thick brittle fault zone along a discrete detachment fault.


Subduction along the Andean margin of South America began by at least the Jurassic and cycled between steep– and shallow–to–flat slab subduction until its present configuration (James, 1971; Jordan et al., 1983; Ramos, 1999). The Cordillera Blanca range is situated above the Peruvian flat-slab segment of the present margin (Gutscher et al., 1999; McNulty and Farber, 2002), within the central-eastern Cordillera Occidental along the western margin of the Marañón fold-and-thrust belt (Fig. 1). West of the Cordillera Blanca lies the Cordillera Negra, a 54–15 Ma volcanic arc related to normal subduction along the Peru-Chile trench (Petford and Atherton, 1992; Scherrenberg et al., 2014, 2016). Beginning at ca. 15 Ma, magmatism migrated eastward, culminating in the intrusion and crystallization of the Cordillera Blanca batholith from 13.7 ± 0.3 Ma to 5.3 ± 0.3 Ma (U-Pb zircon; Mukasa, 1984; Giovanni, 2007) and eruption of the associated Yungay and Fortaleza ignimbrites between 8.7 ± 1.6 and 4.5 ± 0.2 Ma (40Ar/39Ar biotite; Bonnot, 1984; Giovanni et al., 2010) and between 6.2 ± 0.2 and 4.9 ± 0.2 Ma, respectively (K-Ar biotite; Wilson, 1975; Cobbing et al., 1981). These events mark the last episode of arc magmatism in this region (Cobbing et al., 1981; Bonnot, 1984; Mukasa, 1984; Petford and Atherton, 1996; McNulty et al., 1998; McNulty and Farber, 2002; Giovanni, 2007; Giovanni et al., 2010). Crustal shortening associated with convergence of the South American and Nazca plates also migrated inboard, east of the Marañón massif, to the present Subandean fold-and-thrust belt (Fig. 1).

The Cordillera Blanca range is cored by the granodiorite-leucogranodiorite Cordillera Blanca batholith (Atherton and Sanderson, 1987; Petford and Atherton, 1992). Batholith emplacement is estimated at ∼300 MPa (McNulty and Farber, 2002) and >200 MPa (Petford and Atherton, 1992) from mineral assemblages in the pelitic contact aureole, and at 90 ± 10 MPa to 260 ± 30 MPa at 720–800 °C from recent amphibole thermobarometry of batholith rocks (Margirier et al., 2016). The batholith cooled rapidly at 200 °C/m.y. until ca. 4 Ma and subsequently cooled at 25 °C/m.y. from 4 to 0 Ma (Margirier et al., 2015). Extension along the Cordillera Blanca detachment initiated at ca. 5 Ma (Bonnot, 1984; Giovanni, 2007), accommodating 12–15 km of normal-sense displacement (Giovanni et al., 2010). A series of glacially and fluvially incised ravines (quebradas) cut across the Cordillera Blanca nearly orthogonal to strike of the range (Wise and Noble, 2003), providing exposure of the full ductile-to-brittle fabric gradient developed in the CBSZ (Fig. 2). Fault scarps along the detachment truncate glacial moraines, offsetting moraine crests by ∼2–70 m (Bonnot et al., 1988; Schwartz, 1988; Siame et al., 2006). The most recent faulting is bracketed by detrital charcoal in colluvium between 2440 ± 1060 and 750 ± 80 14C yr B.P, and by positioning of pre-Incan structures relative to fault scarps at 1500–2000 yr B.P. (Schwartz, 1988). The hanging wall comprises Jurassic to Cretaceous sediments, Eocene–Pliocene volcanics, and Miocene–Pliocene sediments of the Lloclla Formation (Cobbing et al., 1981; Petford and Atherton, 1992; Giovanni et al., 2010; Margirier et al., 2015). Hot springs issue along the main detachment system and along steeply dipping faults and inferred synthetic faults proximal to the shear zone (Newell et al., 2015).

We present data from samples collected in Quebrada Gatay (Fig. 2) and describe CBSZ deformation relative to structural depths below the exposed detachment surface. Quebrada Gatay samples (Table S1 in the Supplemental Files1) were collected across the ductile-to-brittle fabric transition and are mostly leucogranodiorite, with a small number from quartz veins oriented parallel to foliation and phyllites presumably from a lens of the Jurassic Chicama Formation host rock. Supplemental information is available in Tables S1–S2, Figures S1–S2, and Text S1–S3 (footnote 1).


Thirty-two (32) samples were collected from Quebrada Gatay in July 2013 and were cut into thin sections perpendicular to foliation and parallel to lineation (X-Z sections). We used optical petrography, electron backscatter diffraction (EBSD) analysis, and electron probe microanalyzer (EPMA) analysis to evaluate microstructures, quartz CPOs, deformation temperature, and flow stress across the fabric progression from undeformed granite to ultramylonite and cataclasite.

Deformation Temperatures

Two-Feldspar Geothermometry

Deformation temperatures were estimated for samples with asymmetric strain-induced myrmekite (ASIM) using two-feldspar thermometry of coexisting feldspars within the symplectite (Simpson and Wintsch, 1989; Langille et al., 2010). ASIM develops when aqueous Na+ is concentrated at high-stress sites along K-feldspar grain boundaries (Fig. 3A; Simpson and Wintsch, 1989). High strain facilitates the reaction (Simpson and Wintsch, 1989; Tsurumi et al., 2002; Menegon et al., 2008):

The resulting ASIM microstructure comprises a K-feldspar porphyroclast impinged by asymmetrically distributed albite-oligoclase lobes that grow parallel to the maximum incremental shortening axis and contain vermicular quartz. K-feldspar recrystallizes as tails in low-strain sites (Fig. 3A; Simpson and Wintsch, 1989).

We derive ASIM deformation temperatures using chemical compositions from along the interface of the K-feldspar porphyroclast and the Na-feldspar myrmekitic lobe (Fig. 3B; Langille et al., 2010). Mineral compositions were obtained using the Cameca SX-100 EPMA at the University of Tennessee–Knoxville (UTK; Knoxville, Tennessee, USA) by spot analysis at 15 keV, 10 nA, 5 µm spot size, and 30 s counting time. X-ray maps (Na, K, Si, and Fe) were made at 15 keV, 20 nA, 50 ms dwell time, and 1 µm spot size to evaluate compositional zoning in porphyroclasts and myrmekite lobes and to characterize matrix phases. Multiple lobes from individual ASIM-bearing porphyroclasts were analyzed to determine the average deformation temperature. Temperatures were calculated using the two-feldspar thermometer (Stormer, 1975; Whitney and Stormer, 1977):
where Xab,AF is normalized molar concentration of albite in alkali feldspar, P is pressure (MPa), and αab,PL is the activity of albite in plagioclase. Because shear sense from ASIM is consistent with shear sense from other indicators, temperatures derived from ASIM are assumed to be deformation temperatures (Simpson and Wintsch, 1989; Langille et al., 2010). Temperatures were calculated at 100 MPa and 300 MPa pressures and are presented here calculated at 300 MPa with errors of the 2σ standard deviation of feldspar compositions propagated through the thermometer. Reported temperatures and errors do not include the ±50 °C error from the thermometer. Mineral compositions used for thermometry calculations are reported in Table S2 (footnote 1).

Titanium-in-Quartz (TitaniQ) Geothermometry

Titanium-in-quartz (TitaniQ) thermometry utilizes the log-linear relationship between increasing Ti concentration in quartz with increasing temperature (Wark and Watson, 2006). Established as a reliable thermometer in igneous and volcanic systems when pressure and Ti activity (αTiO2) are constrained, TitaniQ has also been increasingly used in greenschist facies shear zones to determine deformation conditions (e.g., Kidder et al., 2013). Some debate occurs regarding the reliability of TitaniQ in different dynamic recrystallization regimes, though there is general agreement that TitaniQ is reliable under conditions required for grain-boundary migration (GBM) recrystallization (Thomas et al., 2010; Behr et al., 2011; Grujic et al., 2011; Thomas and Watson, 2012; Kidder et al., 2013; Nachlas et al., 2014; Cross et al., 2015). Recent experimental evidence suggests that Ti concentrations are efficiently re-equilibrated during recrystallization by grain boundary formation and migration (Nachlas et al., 2018). Because high-temperature and/or low-strain-rate conditions tend toward preserving Ti concentrations on the order of 10–100 ppm, concentrations can be determined on an EPMA for quartz exhibiting GBM (high-temperature and/or low-strain-rate) recrystallization microstructures (Behr et al., 2011; Grujic et al., 2011; Nachlas et al., 2014, 2018). Samples without GBM recrystallization microstructures were not considered for analysis in this study due to the likelihood that Ti concentrations would be below the detection limit of the EPMA at UTK, or otherwise may be preserving inherited Ti concentrations from higher-temperature deformation (i.e., Grujic et al., 2011). Although Nachlas et al. (2018) suggested that Ti should be re-equilibrated for recrystallized quartz even at low temperatures, resolving possible inherited concentrations versus recrystallized Ti concentrations would require higher-precision instrumentation.

Cathodoluminescence (CL) imaging was conducted prior to EPMA analyses. CL imaging details and EPMA conditions for TitaniQ analyses are provided in Text S1 (footnote 1). Temperatures were calculated using the thermobarometer of Thomas et al. (2010):
where P is pressure (MPa), R is the gas constant, forumla is the concentration of Ti in quartz (ppm), and αTiO2 is the TiO2 activity, and the thermobarometer of Huang and Audétat (2012):
where T is temperature (Kelvin) and P is pressure (MPa). TiO2 activity for silicic igneous rocks can be difficult to constrain (Huang and Audétat, 2012), although for most igneous and metamorphic rocks, αTiO2 is above ∼0.5, and in silicic melts, typically above ∼0.6 (Hayden and Watson, 2007). Activity can be further constrained by presence of ilmenite or titanite (αTiO2 0.7–0.8) (Ghent and Stout, 1984; Peterman and Grove, 2010).

Quartz Crystallographic Preferred Orientations (CPOs)

Quartz CPOs were obtained using EBSD on the field emission gun–scanning electron microscope at the Imaging and Chemical Analysis Laboratory (ICAL) at Montana State University (Bozeman, Montana, USA). Analyses were conducted at 20 Pa and 20–30 kV, and Kikuchi band diffraction patterns were collected for quartz for 75 reflectors using Oxford Instruments HKL Channel 5 Flamenco software (ver. 5.5). Step sizes ranged from 1.5 to 5 µm depending on minimum grain size observed under petrographic microscope. Pole figures of quartz c- and a-axis orientations were constructed using the MTEX toolbox (Bachmann et al., 2010) for MATLAB. Details on raw EBSD data processing are provided in Text S2 (footnote 1). EBSD was applied to eight representative samples in interconnected quartz-rich domains as far away from other phases as possible to minimize potential effects from pinning. CPO pole figures are lower-hemisphere projections of c- and a-axes using one point per grain, where grains are defined using at least eight pixels and a critical misorientation of 10°, contoured at 7.5° half-width (see Text S2 [footnote 1] for details). Fabric strength was evaluated using the M-index of Skemer et al. (2005) and the J-index of Bunge (1982).

Flow Stress

Differential stresses were estimated using recrystallized grain-size piezometry, where differential stress and recrystallized grain size are inversely related by the following expression (Twiss, 1977; Stipp and Tullis, 2003):

Differential stresses were calculated using the piezometer calibration of Stipp and Tullis (2003), modified after Holyoke and Kronenberg (2010), where d is the root mean square (RMS) of the grain-size distribution, β is an empirically derived constant (3640 μm/MPam), σ is differential stress, and m is an experimentally derived stress exponent (−1.26). These values were compared to the recent “sliding-resolution” piezometer calibration for EBSD-derived grain sizes, where β and m are 104.22 ± 0.51 and −1.59 ± 0.26, respectively (Cross et al., 2017b). Grain-size distributions (Fig. S1 [footnote 1]) were derived from EBSD analysis on representative quartz-rich domains from mylonitic granodiorite and mylonitic quartz veins, where grain size was defined as the diameter of an equivalent-area circle. Errors on differential stress estimates are the standard deviation relative to RMS for the grain-size distribution propagated through the piezometer.


Mesoscale Observations

Mean foliation from Quebrada Gatay strikes 167°± 5° and dips 23° ± 5°, and mean stretching lineation trends toward 239° ± 6°, plunging 24° ± 6° (Fig. 4A). The detachment surface exposed at Quebrada Gatay (Fig. 4A) strikes ∼150° and dips ∼26° based on a series of three-point problems solved using Google Earth. Data from structural positions between 150 and 300 m below the detachment are sparse due to steep terrain and vegetative cover. An injection complex of leucogranodiorite veins intrudes into a host rock (shale) lens 130–150 m below the detachment surface (Fig. 4B). Localized zones of mylonite to ultramylonite from 5 to 50 cm thick occur at structural positions as deep as 385 m below the detachment within otherwise undeformed to protomylonitic granite (Fig. 4C). Quartz veins 10–60 cm thick occur parallel and subparallel to foliation throughout the shear zone (Fig. 4D). A progression in fabric from undeformed granodiorite to ultramylonite with consistent SSW-dipping foliation is observed (Fig. 4E) and is the basis for interpretations of deformation mechanisms and conditions during shear zone evolution.

Microstructural Observations

Observed microstructures can be divided based on interconnected phases into feldspar, quartz, mica, and polyphase domains. These domains are spatially variable such that the rheological framework may be controlled by either quartz, feldspar, mica, or a combination thereof.


Within ∼150 m of the detachment, throughgoing features such as fractures, microfaults, and cataclasite cut across feldspar and adjacent phases (Fig. 5A, Table 1). At structurally high positions, feldspar porphyroclasts are commonly fractured and filled with quartz and/or K-feldspar. These fractures tend to be isolated to feldspar porphyroclasts and do not extend into surrounding phases. Some feldspar fragments are rotated with respect to the other fragment(s) within the original porphyroclast, forming V-shaped pull-apart structures (Fig. 5B; i.e., Fukuda et al., 2012; Hippertt, 1993; Samanta et al., 2002). Undulose extinction and deformation lamellae are common within plagioclase porphyroclasts.

Fine-grained feldspar aggregates occur in layered mylonitic to ultramylonitic samples at structural depths from 430 to 111 m. Finely recrystallized grains mantle porphyroclasts and form tails on σ- and δ-type feldspar porphyroclasts that deflect into quartz + feldspar shear bands (Fig. 5C). Within these samples, myrmekite commonly has formed along high-strain sites of K-feldspar porphyroclasts as well as within fine-grained plagioclase–K-feldspar aggregates (Fig. 5D).

Interconnected networks of feldspar porphyroclasts or phenocrysts occur below 350 m. At 350 m, mylonitic samples contain interconnected plagioclase phenocrysts that exhibit bulging recrystallization at feldspar-feldspar grain boundaries (Fig. 5E) and contain biotite at triple junctions (Fig. 5F). Biotite also occurs in low-strain sites of boudinaged plagioclase where interconnected quartz and interconnected feldspar coincide. Plagioclase boudins are relatively symmetrical, are tapered yet connected, and have concave-faced boundaries (Fig. 5G). Weakly deformed samples at the deepest structural positions are characterized by interconnected networks of K-feldspar and plagioclase phenocrysts that exhibit intragranular fractures (Figs. 5H, 5I). Recrystallization has occurred along plagioclase grain boundaries and along intragranular fractures within plagioclase and K-feldspar phenocrysts (Fig. 5I). Large (∼500 µm) myrmekite lobes have grown into K-feldspar phenocrysts (Fig. 5J).


Quartz domains occur primarily either as interconnected networks surrounding other phases or as monophase aggregates in layered mylonites. Quartz commonly displays grain-shape fabrics that define S-planes and exhibits dynamic recrystallization microstructures spanning the range of characteristic grain-boundary microstructures described by Stipp et al. (2002a, 2002b; Fig. 6). Distribution with respect to structural depth of quartz recrystallization microstructures, along with that of feldspar microstructures discussed above and microfabrics including interconnected phase, inferred quartz crystallographic slip systems, and deformation temperatures discussed in the following sections, are summarized in Figure 7.

In the upper ∼150 m of the shear zone, ductile fabrics including bulging (BLG) and subgrain rotation (SGR) recrystallization microstructures in quartz are overprinted by brittle fabrics (Figs. 6A–6D; Table 1). Serrated grain boundaries and fine recrystallized grains characteristic of BLG recrystallization (Stipp et al., 2002b) occur in samples that also contain polygonal subgrains and new rotated grains within quartz porphyroclasts, characteristic of SGR recrystallization (Stipp et al., 2002a; Fig. 6E). Samples within the upper ∼150 m of the shear zone that exhibit only ductile fabrics contain microstructures indicative of SGR recrystallization or a combination of BLG and SGR.

Samples from deeper than ∼150 m below the detachment preserve ductile deformation in quartz and lack evidence for throughgoing brittle deformation. Quartz in most samples at these positions exhibits dominant SGR recrystallization microstructures, large grains with amoeboidal grain boundaries and commonly undulose extinction indicative of grain boundary migration (GBM) recrystallization (Stipp et al., 2002a), or both SGR and GBM recrystallization microstructures (Figs. 6F–6J). Two samples from 384 to 386 m below the detachment record a combination of BLG and SGR recrystallization in quartz (Fig. 6I).

At structurally deep positions, quartz forms discontinuous lenses and/or fills fractures in feldspar phenocrysts (Figs. 5H, 6J, 6K). Structurally deepest positions preserve large amoeboid quartz grains characteristic of GBM recrystallization microstructures within isolated quartz pockets at the interstices between feldspar phenocrysts (430 m; Fig. 5H). Planar, orthogonal subgrain boundaries indicative of chessboard extinction in quartz occur at deepest structural positions (446 m; Figs. 5H, 6K).


Microfabrics described here refer to the interconnected phase(s) present as well as which phases define S-, C-, and C′-surfaces and which of these surfaces is prominent within the foliation.

S-C mylonites. S-C mylonites are broadly characterized by prominent interconnected quartz domains with S-planes defined primarily by quartz grain-shape fabrics and mica. Type I S-C fabrics as described by Berthé et al. (1979) and Lister and Snoke (1984), where S-surfaces anastomose into and out of C-surfaces and C′-surfaces, are relatively well developed (Figs. 8A, 8B), and type II S-C fabrics, which contain prominent C-surfaces defined by tails of sheared mica fish and S-surfaces defined by quartz grain-shape fabrics (Fig. 8C), are observed in the CBSZ. Quartz in S-C mylonites surrounds isolated feldspar porphyroclasts that contain ASIM or are fractured with quartz infill (Fig. 8C). Biotite or white mica shear bands in some cases have become interconnected, enveloping feldspar porphyroclasts and dynamically recrystallized quartz clasts (Fig. 8D).

Banded mylonites. Banded mylonites comprising interlayered monophase aggregates of quartz or feldspar and polyphase aggregates of feldspar + quartz ± mica up to ∼1500 µm thick (Figs. 8E–8I) are observed between 111 and 430 m structural depth (Fig. 7; Table 1). Boundaries between layers tend to align with mesoscopic foliation, which are typically the most prominent foliation surfaces (Fig. 8E). Fine-grained feldspar ± mica bands alternate with monomineralic quartz layers containing microstructures characteristic of SGR and/or GBM recrystallization (Figs. 8F, 5D). Quartz + feldspar domains exhibit relatively equant grain shapes and grain sizes smaller than those of either feldspar or quartz monophase domains, and occur as continuous bands up to ∼500 µm thick (Fig. 8G). Domains of finely recrystallized quartz + feldspar + mica appear to originate at plagioclase or K-feldspar porphyroclasts and extend into bands up to ∼1.5 mm thick that are continuous at the thin-section to hand-sample scale (Fig. 8H). Fine-grained aggregates contain plagioclase with more albitic composition than porphyroclasts intermixed with K-feldspar, quartz, and biotite (Fig. 8I). K-feldspar occupies dilatant sites and triple junctions within quartz aggregates and between plagioclase grains (Fig. 8I).

Deformation Temperatures

Two-Feldspar Geothermometry

Deformation temperatures derived from ASIM for six samples from structural depths of 134 m to 446 m below the detachment range from 410 ± 30 °C to 470 ± 36 °C (Table 2; Fig. 7). Temperatures do not vary systematically with structural depth and instead overlap within the 50 °C error of the thermometer.

TitaniQ Geothermometry

Deformation temperatures at structural depths from 134 to 446 m are between 451 ± 60 and 489 ± 33 °C, using the Thomas et al. (2010) TitaniQ geothermometer, assuming αTiO2 of 0.5 and 300 MPa pressure. This assumption is made based on the absence in samples of a Ti-bearing phase other than biotite that is in clear equilibrium with quartz. Ti-bearing phases such as rutile and/or ilmenite are present only in three analyzed samples and were either along cleavage planes in biotite or within fractures that cut across recrystallized grain boundaries. The assumed activity of 0.5 is in line with Ti activities determined by geochemical modeling for granitoids of similar composition at similar pressures to the those of the CBSZ (Huang and Audétat, 2012) as well as experimentally grown quartzites where rutile is present but sparse (>200 µm between grains; Nachlas and Hirth, 2015). Ti activities as low as 0.2 have been reported in some plutons as well as volcanic rocks with similar compositions to the Cordillera Blanca batholith (Huang and Audétat, 2012; Kularatne and Audétat, 2014). Temperatures are ∼70 °C hotter for αTiO2 of 0.2 and ∼30 °C cooler for an activity of 0.8 (Table 2). Temperatures estimated from the Huang and Audétat (2012) TitaniQ calibration range from 520 ± 66 °C to 574 ± 30 °C. These values are within 10 °C of temperatures estimated using the Thomas et al. (2010) calibration for αTiO2 of 0.2 (Fig. 7, Table 2). TitaniQ thermobarometry demonstrates a slight decrease in temperature at decreasing structural depths (Fig. 7). However, this systematic variation is not readily observed with ASIM thermometry (Fig. 7, Table 2). Conservatively, deformation temperatures from TitaniQ thermobarometry are between 451 ± 60 °C and 574 ± 30 °C. TitaniQ temperatures for the CBSZ should be received with caution because, in addition to the uncertainty regarding αTiO2 in the CBSZ, the measured Ti concentrations are at or near the detection limit of ∼11 ppm for all samples (see Text S3 [footnote 1] for additional details).

Quartz CPOs

Quartz CPO patterns indicate predominant prism <a> and rhomb <a> slip (Fig. 9). Most samples from 94 m to 380 m from the detachment show slightly elongated y-maxima characteristic of dominant prism <a> slip with subsidiary rhomb <a> slip. Samples at 384 m and 430 m below the detachment had weak CPOs with M-indices ≤0.07. CPO strength otherwise appears uncorrelated with structural depth or slip system. At 152 m and 369 m, a- and prism-normal-axes do not form the six maxima expected for <a> slip and are instead smeared into three dispersed maxima.

Flow Stress

Differential stress, hereafter referred to as “stress”, estimates (Table 3) range from 16.5 ± 13 MPa to 71.5 ± 36 MPa using the Holyoke and Kronenberg (2010) piezometer, measured at 350 m (sample CB13-59) and 149 m (sample CB13-76), respectively (Fig. 10). The structurally highest sample, CB13-77c (94 m), yields a stress of 40.0 ± 22 MPa. The highest stress of 71.5 ± 36 MPa is from the deepest sample that contains throughgoing brittle overprinting of ductile fabrics (sample CB13-76, 149 m). Intermediately positioned sample CB13-78 (152 m) records the next highest stress of 64.3 ± 26 MPa. Sample CB13-59 (350 m) yields the lowest stress estimate of observed samples, 16.5 ± 13 MPa. Further below the detachment, sample CB13-55c (369 m) yields a stress of 48.2 ± 23 MPa. Stresses calculated at lower structural positions are 56.3 ± 26 MPa and 59.3 ± 27 MPa, from samples CB13-57a (384 m) and CB13-57c (386 m), respectively. The lowest structural position analyzed (sample CB13-55a, 430 m) yields a stress of 32.1 ± 14 MPa. Stress estimates using the Stipp and Tullis (2003) and Cross et al. (2017b) calibrations range from 22.5 ± 19 to 97.6 ± 59 MPa and 30.7 ± 17 to 98.1 ± 38 MPa, respectively and are within error of those calculated using the Holyoke and Kronenberg (2010) calibration (Fig. 10, Table 3). In summary, a consistent gradient in stress as a function of structural depth is not observed; rather, the highest stresses are observed at intermediate structural positions near 150 m where the deepest throughgoing brittle overprinting is observed.


Quartz CPOs

It is important to note that CPOs may be inherited from earlier deformation. Strong preexisting growth CPOs in quartz veins have been shown to influence CPOs of quartz veins deformed at lower temperatures, where quartz initially grew in an orientation favorable for easy slip at high temperatures (e.g., prism <a> slip at >500 °C) (Pennacchioni et al., 2010). Grains were oriented unfavorably for basal <a> slip when temperatures decreased during deformation such that even at relatively low temperatures, strong CPOs indicative of prism <a> ± rhomb <a> slip were preserved (Toy et al., 2008; Pennacchioni et al., 2010). Previous work also suggests that oriented grain growth during submagmatic flow may influence later CPO development during subsolidus deformation (Zibra et al., 2012). We thus consider the possibility that strong CPOs at positions nearer to the detachment may be inherited from early deformation at relatively high temperatures where prism <a> slip was initially favored, or from oriented grain growth during submagmatic flow. The weak fabric in the deepest sample analyzed for EBSD (sample CB13-55a, 430 m) resembles an elongate y-axis maximum expected for prism <a> ± rhomb <a> slip that is seen in stronger CPOs at structurally higher positions (Fig. 9). Although this could suggest some degree of inheritance from a preexisting CPO in the weakly deformed granodiorite, ASIM and TitaniQ estimates suggest that deformation for most of the CBSZ occurred at temperatures above ∼400 °C, which is consistent with temperatures where prism <a> and rhomb <a> slip are expected to occur more readily than basal <a> slip (Law et al., 1990; Toy et al., 2008). CPO inheritance is a possibility in the CBSZ but is not required to explain our data.

Deformation Mechanisms


Fracturing. Feldspar exhibits fracturing throughout most of the CBSZ (Fig. 7). At structurally higher positions, fractures occur predominantly in plagioclase porphyroclasts that are surrounded by quartz and/or fine-grained feldspar + quartz ± mica domains. At structurally deep positions, K-feldspar and plagioclase phenocrysts form an interconnected framework that features inter- and intragranular fractures.

Previous studies have documented fracturing at submagmatic to subsolidus conditions in networks of load-bearing interconnected feldspar (Rosenberg, 2001; Zibra et al., 2012; Bessiere et al., 2018). Syn-magmatic fractures in granitoids are characterized by K-feldspar and plagioclase films along grain boundaries at high angles to foliation planes and by fractures filled with at least one phase that would crystallize from a residual melt (Rosenberg, 2001). These characteristics are consistent with those described by Bouchez et al. (1992) for submagmatic fractures at melt fractions between 0.10 and 0.30. Submagmatic fractures are described as those: (1) that transect single grains, (2) that are filled with a phase that is crystallographically and compositionally continuous inside and outside of the fracture, and (3) whose filling phase is expected to crystallize from the residual melt (Bouchez et al., 1992). Within the CBSZ, fractures satisfying these criteria occur at structurally deep positions, where quartz exhibiting the same orientation extends from a wedge-shaped pocket into an adjacent narrow fracture (Fig. 5H, 446 m). The amount of melt required to produce submagmatic fractures has been debated. Melt fractions of 0.3–0.5 had been argued to be required for the formation of throughgoing melt-filled fractures (Arzi, 1978; Paterson et al., 1989), though further investigation into the experimental work from whence these values were derived suggest that melt-filled fractures could be produced at lower melt fractions, potentially as low as 0.07 (Rosenberg and Handy, 2005).

Although quartz crystallization would be consistent with the third criterion above as set forth by Bouchez et al. (1992), other evidence for magmatic deformation or submagmatic deformation, such as imbricated feldspar phenocrysts, is not observed in the CBSZ. The wedge-shaped nature of the quartz exhibiting chessboard extinction is consistent with some degree of rotation of the fractured K-feldspar (Hippertt, 1993), which must have occurred at or above the temperature required for chessboard extinction in quartz (∼630 °C; Kruhl, 1996), placing it in the field of high-temperature subsolidus deformation. Alternatively, the quartz exhibiting chessboard extinction may have been late-stage melt filling a pocket between plagioclase and K-feldspar porphyroclasts, with submagmatic fracturing occurring at the melt-porphyroclast interface of the melt pocket at a low melt fraction, still near the boundary between submagmatic and subsolidus conditions (Paterson et al., 1989). In either scenario, evidence for the role of fracturing in feldspar deformation at high temperatures is observed at deepest structural positions within the CBSZ.

Dislocation creep versus cataclastic flow. Deformation experiments on plagioclase single crystals with no initial microstructure other than growth twins suggest that feldspar experiences complex semi-brittle cataclastic flow, where brittle processes (i.e., microfracturing) and plastic processes (i.e., deformation twinning via dislocation climb) are interdependent (McLaren and Pryer, 2001). The coincidence of deformation lamellae with fractured feldspar porphyroclasts and undulose extinction within the upper ∼180 m of the CBSZ may suggest the activity of this type of semi-brittle flow. These microstructures are also consistent with those described for cataclastic flow by Tullis and Yund (1987). In their albite aggregate deformation experiments, Tullis and Yund (1987) described characteristic microstructures of cataclastic flow that include undulatory and/or patchy extinction and flattening, grain-scale faults with small offset, and microcracks. The authors noted that patchy undulose extinction, which has also been interpreted to represent dislocation creep (Tullis and Yund, 1977), is typically caused by microcrack arrays and crushed zones visible at the submicron scale.

Evidence for the activity of dislocation creep in feldspar is also observed within the CBSZ. Dynamic recrystallization in feldspar can produce relatively small, rounded, and undeformed feldspar porphyroclasts (Tullis and Yund, 1985). During dislocation creep in feldspar, recrystallization is primarily accommodated by migrating grain boundaries (Tullis and Yund, 1987), producing grain-boundary bulges and small strain-free recrystallized grains (Stünitz et al., 2003) characteristic of BLG recrystalliztion microstructures. Additional related microstructures can include undulatory extinction in relatively equant grains and augen surrounded by a fine-grained matrix at relatively high strains. Undulose extinction due primarily to dislocation creep sweeps continuously across a grain, as opposed to the patchy undulose extinction that can result from dislocation entanglement (Hirth and Tullis, 1992) or cataclastic flow (Tullis and Yund, 1987). Fine, equant, relatively strain-free, recrystallized grains occur along the margins of porphyroclasts, grain boundaries, and fractures, commonly extending out into tails under high strains (Tullis and Yund, 1985).

Within the CBSZ, fine recrystallized feldspar occurs at a wide range of structural depths (92–446 m; Table 1; Fig. 7). As these microstructures may be formed by either cataclastic flow or dislocation creep, interpreting deformation mechanisms by singular microstructures alone is difficult. Here, interpretation of the activity of dislocation creep in feldspar is limited to samples where grain-boundary bulges can be observed (i.e., Fig. 5E) and where undulatory extinction, if present, is continuous across grains rather than patchy. If patchy undulose extinction occurs, especially if coinciding with kinked or bent twins or cleavage planes, or pull-apart–type fractures, the prominent deformation mechanism is interpreted to be cataclastic flow. Under these considerations, dislocation creep in feldspar likely occurred at relatively deep structural positions (350 m to 446 m), suggesting deformation temperatures above ∼450 °C (Yund and Tullis, 1980), consistent with ASIM-derived temperatures of 440 ± 46 to 470 ± 36 °C and TitaniQ temperatures of 484 ± 99 to 489 ± 39 °C (Thomas et al. [2010] calibration; αTiO2 = 0.5) for the same structural positions (Fig. 10). Alternatively, cataclastic flow likely occurred at structural positions above ∼350 m. This relationship may be oversimplified, as a cataclastic overprint (in feldspar) could reduce evidence for dislocation creep, but also more simply because of the similarity of microstructures between cataclastic flow and dislocation creep (i.e., Tullis and Yund, 1987).


Dislocation creep. Recrystallization microstructures and CPOs can be used to distinguish between prominent deformation mechanisms (Hirth and Tullis, 1992; Stipp et al., 2002b). Recrystallized quartz CPOs mainly show patterns indicative of dominant prism <a> slip with variable lesser contributions of rhomb <a> slip, which is consistent with quartz dislocation creep at 94–386 m below the detachment surface (Fig. 9). Previous studies suggest that temperature is the primary contributing factor to activating different slip systems in quartz aggregates, with the transition from basal slip to prism slip occurring at high temperatures (Tullis et al., 1973; Mainprice et al., 1986; Morales et al., 2011). Above ∼500 °C, prism <a> occurs as the dominant slip system, and at increasingly high temperatures of ∼600 °C, slip in the <a> direction gives way to dominant prism [c] slip (Mainprice et al., 1986; Schmid and Casey, 1986; Kruhl, 1996; Passchier and Trouw, 2005). However, transitions between dominant slip mechanisms inferred from quartz CPOs may also be affected by changes of strain rate and/or the activity of water within the shear zone (Fig. 9; Mainprice and Paterson, 1984; Stipp et al., 2002b; Mancktelow and Pennacchioni, 2004; Law, 2014; Morales et al., 2014). Excluding the weak fabrics of samples CB13-57a and CB13-55a (384 and 430 m, respectively), quartz CPOs indicate predominant activity of prism <a> slip and rhomb <a> slip, which are consistent with dislocation creep in quartz at >400 °C (Schmid and Casey, 1986; Stipp et al., 2002a, 2002b; Stünitz et al., 2003).

Grain-boundary sliding. Weak CPOs at 384 and 430 m depth, as well as dispersed a-axis maxima at 152 m and 369 m depth, may suggest the activity of other deformation mechanisms, such as diffusion creep and/or grain-boundary sliding (GBS) (Bestmann and Prior, 2003; Wightman et al., 2006; Cross et al., 2017a) or relatively low amounts of accumulated strain (Lister and Hobbs, 1980; Toy et al., 2008). Diffusion creep is unlikely in coarse-grained quartz aggregates (Rutter and Brodie, 2004), such as those in sample CB13-55a (430 m). The weak CPO in sample CB13-55a more plausibly reflects relatively low strain, which is within reason for its deep structural position (430 m) relative to the detachment fault. This relationship is consistent with weak CPOs at positions structurally removed from the main detachment and/or shear zone in other natural settings (Haertel et al., 2013; Haertel and Herwegh, 2014) or in low finite strains in numerical models (Keller and Stipp, 2011; Lister and Hobbs, 1980).

The weak CPO at 384 m (sample CB13-57a) and the dispersed a-axis patterns at 152 and 369 m depth (samples CB13-78 and CB13-55c, respectively) may reflect GBS. Dispersed CPO patterns are often attributed to GBS (Jiang et al., 2000; Bestmann and Prior, 2003; Storey and Prior, 2005; Warren and Hirth, 2006; Cross et al., 2017a). Microstructures such as four-grain junctions, square or rectangular grain shapes, alignment of multiple grain boundaries parallel to foliation, and recrystallized grains smaller than subgrains are thought to be indicative of GBS (Halfpenny et al., 2012; Miranda et al., 2016; Rahl and Skemer, 2016). CPOs of sub-populations at 384 m (sample CB13-57a) show that the weak c-axis maxima derive mainly from relict grains and that recrystallized grains have dispersed orientations (Fig. S2 [footnote 1]). In this sample, four-grain junctions associated with square grain shapes and alignment of edges across multiple grains also occur within the finer recrystallized population (Fig. S2). These observations, in addition to the dispersion of the CPOs in the fine-grained recrystallized quartz population, are consistent with the activity of GBS in quartz (Bestmann and Prior, 2003; Halfpenny et al., 2006; Wightman et al., 2006; Gan et al., 2007; Kilian et al., 2011; Cross et al., 2017a).

Polyphase Domains

Microstructural observations suggest that polyphase quartz + feldspar ± mica domains may be promoted by phase mixing via GBS, dynamic recrystallization, and pressure-solution creep. Precipitation of soluble phases such as K-feldspar or quartz in dilatant openings, such as in recrystallized porphyroclast tails or pull-apart zones in fine-grained aggregates, has previously been interpreted to represent GBS (Behrmann and Mainprice, 1987; Kilian et al., 2011; Platt, 2015). At least local GBS is supported by K-feldspar in dilatant sites and at triple junctions within quartz aggregates and between plagioclase grains (Fig. 8I; e.g., Stünitz and Tullis, 2001; Fusseis et al., 2009; Oliot et al., 2014; Menegon et al., 2015; Spruzeniece and Piazolo, 2015; Czaplińska et al., 2015; Viegas et al., 2016). Disaggregation of dynamically recrystallized feldspar porphyroclast tails into polyphase quartz + plagioclase + K-feldspar ± biotite domains (Fig. 8H) is consistent with microstructures inferred to form from precipitation in dilatant openings following dynamic recrystallization in porphyroclast tails, where rotation of grains opens up pore space that is subsequently filled by K-feldspar (Kilian et al., 2011; Platt, 2015).

Strain Rates

Strain rates derived from the wet quartz flow law of Hirth et al. (2001) are 2.6 × 10‒14 to 2.6 × 10‒12 s‒1 using stresses derived from the Holyoke and Kronenberg (2010) piezometer and temperatures from TitaniQ (αTiO2 = 0.5) and ASIM thermometry. Similar strain rates of 3.5 × 10–14 to 3.6 × 10–12 s–1 are obtained using the new quartz flow law for prism <a> slip (Tokle et al., 2019) (Table 3). Strain rates increase to 3.1 × 10‒13 to 9.4 × 10‒12 using the Cross et al. (2017b) piezometer. Highest calculated strain rates (2–9 × 10‒12 s–1) are from ∼150 m below the detachment, near the position of the deepest throughgoing brittle overprint. Displacement of 12–15 km along the detachment since ca. 5 Ma (Bonnot, 1984; Giovanni, 2007; Giovanni et al., 2010) yields a time-averaged slip rate of 2.4–3 mm/yr. If it is assumed that the width of the shear zone between a given sample and the detachment was active during shear zone evolution, shear zone widths prior to onset of brittle deformation were between ∼450 m and ∼150 m. These widths correspond to time-averaged strain rates from 1.7 to 6.3 × 10‒13 s‒1 for 450 and 150 m widths, respectively, under the assumption that strain rate = velocity / width (Platt and Behr, 2011b). Despite uncertainties in temperature and stress values, strain rates estimated from individual samples generally agree with the time-averaged strain rate for the CBSZ.

Use of ASIM temperatures in the quartz flow law requires a simplified assumption that myrmekite formation and quartz deformation are coeval. TitaniQ temperatures at αTiO2 = 0.5 overlap with ASIM temperatures within error (Fig. 7; Table 2), supporting the use of ASIM temperatures as an approximation of deformation temperatures in quartz for strain rate calculations when quartz-based temperature data are unavailable. ASIM temperatures are broadly consistent with temperatures >400 °C inferred from quartz CPOs as well as temperatures expected for SGR recrystallization (400–500 °C; Stipp et al., 2002a), and are within error of the transition to GBM recrystallization at strain rates of ∼10‒12 s‒1 (480–530 °C; Stipp et al., 2002a). TitaniQ analyses were selectively conducted on samples with GBM recrystallization microstructures, which at strain rates of ∼10‒12 s‒1 occur at temperatures >500 °C (Stipp et al., 2002a). ASIM temperatures and TitaniQ temperatures for these samples, while in agreement with each other at αTiO2 = 0.5, are below the expected range from Stipp et al. (2002a). Deformation experiments suggest that one-order-of-magnitude variation in strain rate has a similar effect on recrystallization microstructure to a 200 °C shift in temperature, and that an increase in water content changes recrystallization microstructures in a similar matter to an increase in temperature (Stipp et al., 2006).

Shear Zone Evolution

Based on overprinting relationships and structural position relative to the detachment fault, we propose a conceptual model for the evolution of the CBSZ during post-magmatic cooling. In this model, initial deformation in the CBSZ was accommodated by high-temperature feldspar fracturing. Continuing subsolidus deformation occurred at temperatures above ∼400 °C through varying contributions of dislocation creep in feldspar, dislocation creep in quartz, cataclasis in feldspar, and at least locally by GBS in polyphase aggregates composed of quartz + feldspar ± mica. The evolution of these mechanisms in the CBSZ does not follow a single progression. Rather, our data support a conceptual model where banded mylonites and/or interconnected quartz domains formed as a result of multiple deformation mechanisms that likely reflect changes in temperature, strain rate, and presence of fluids prior to the onset of throughgoing brittle deformation (Fig. 11).

Initial Deformation

Earliest deformation occurred via fracturing and neocrystallization of K-feldspar and plagioclase and dislocation creep in quartz (Fig. 11A). Chessboard extinction in quartz within submagmatic, intragranular K-feldspar fractures (Figs. 5H, 11A) represents the onset of subsolidus deformation at >630 °C (Kruhl, 1996). Quartz with chessboard extinction and large grains with amoeboid grain boundaries and undulose extinction indicative of GBM recrystallization are common within isolated pockets at the interstices between feldspar phenocrysts. As such, these quartz domains are interpreted to be passive markers of deformation conditions, and rheology is inferred to be controlled by the interconnected network of feldspar phenocrysts.

Grain-Size Reduction

Grain-size reduction of the feldspar network occurred by dislocation creep in plagioclase, recrystallization to fine-grained feldspar or feldspar + quartz along intragranular fractures in K-feldspar phenocrysts, and myrmekite growth (Figs. 5I, 5J, 11A). BLG recrystallization occurs along fractures and the rims of magmatic plagioclase and K-feldspar grains (Fig. 5E, 5I), suggesting that deformation was partially accommodated by dislocation creep in feldspar during early stages of shear zone evolution. Dissolution along plagioclase and K-feldspar fractures facilitated by aqueous fluids is inferred from the crystallization of finer-grained, more sodic feldspar grains along these sites (Fitz Gerald and Stünitz, 1993; Fusseis and Handy, 2008; Viegas et al., 2016). At structurally deep positions within the CBSZ, prominent myrmekite growth along K-feldspar boundaries that is associated with white mica and quartz (Figs. 5J, 11A) also supports the role of fluids in early grain-size reduction of the load-bearing feldspar network (Menegon et al., 2006; Ree et al., 2005; Tsurumi et al., 2002). The presence of white mica and quartz along the inferred original K-feldspar edge supports myrmekite growth facilitated by a local addition of water, where excess K+ is released and myrmekite and muscovite are formed at the expense of K-feldspar following the reaction (Phillips et al., 1972):

Collectively, these observations suggest a role of fluids in promoting grain-size reduction via fracturing, recrystallization, and myrmekite formation within an initially interconnected network of feldspar phenocrysts at >400 °C within the CBSZ.

Interconnection of Quartz Domains

Quartz deformed by dislocation creep in interconnected networks around other isolated phases at a wide range of structural depths in the CBSZ (Fig. 7). The transition from a coarse-grained feldspar network to interconnected quartz network (Fig. 11B) is preserved at 350 m structural depth (i.e., sample CB13-59, Figs. 5E–5G, 6H), where quartz exhibits GBM recrystallization microstructures and plagioclase is boudinaged. The tapered, relatively symmetric nature of plagioclase boudins suggests that quartz was less competent than feldspar, and flowed into boudin necks (Goscombe et al., 2004). As plagioclase within this sample also exhibits BLG recrystallization (Fig. 5E), deformation at the transition from a coarse-feldspar network to a quartz network is consistent with the model for a two-phase system in which both rheological phases deform viscously (i.e., Handy et al., 1999, their figure 1b). In the CBSZ, this transition occurred under relatively low stress (16.5 ± 13 MPa in quartz) near ∼500 °C, at the lowest strain rate estimated for individual samples, ∼3 × 10‒14 s‒1. Similar differential stresses are reported for extensional mylonitic granitoids exhibiting microstructures indicative of GBM recrystallization in quartz and plastic deformation in feldspar in the core (10–16 MPa) and mylonitic front (∼16 MPa) of the Whipple Mountains metamorphic core complex (California, USA) in the North American Cordillera (Behr and Platt, 2011; Cooper et al., 2017).

Deformation within interconnected quartz domains proceeded predominantly by dislocation creep. GBM, SGR, or SGR + GBM recrystallization microstructures in quartz occur at structural positions deeper than ∼150 m below the detachment (Fig. 11C), whereas SGR ± BLG recrystallization microstructures primarily occur at structural positions shallower than ∼150 m (Fig. 11D), with the exception of a relatively narrow zone at ∼385 m depth (Fig. 7). These microstructures are accompanied by quartz CPOs that mainly show strong fabrics indicative of prism <a> slip and subsidiary rhomb <a> slip, consistent with quartz dislocation creep at >400 °C. The shift in GBM to BLG recrystallization microstructures in addition to SGR microstructures at ∼150 is likely due to increasing strain rate. This microstructural transition is often attributed to decreasing temperature, where GBM recrystallization occurs above ∼500 °C, and BLG recrystallization occurs below ∼400 °C (Stipp et al., 2002b). However, estimated deformation temperatures in the CBSZ are consistently >400 °C and within ∼70 °C of each other, suggesting that temperature is likely not the main contributing factor in microstructural evolution. Furthermore, a change in temperature alone would not explain the narrow zone of BLG recrystallization microstructures near ∼385 m depth, as temperature changes would be expected to be of relatively large scale due to pluton cooling and/or exhumation of the shear zone along the detachment. Although strain rate could only be derived for a limited number of samples, strain rates estimated for samples near the ∼150 m microstructural transition with differential stresses of ∼70 MPa (3 × 10‒12 s‒1; samples CB13-76, CB13-78) are an order of magnitude higher than those estimated in deeper positions, supporting the role of strain rate in shear zone evolution. Strain rate estimates near ∼150 m in the CBSZ are similar to the upper range of those calculated for the Ruby Mountains core complex (Nevada, USA; 10‒11 to 10‒13 s‒1) at similar differential stresses (∼64 MPa) (Hacker et al., 1990).

Layering of Banded Mylonites

Layering of banded mylonites (Fig. 11E) occurs at depths ranging from 430 m to 111 m (Fig. 7). Formation of layered mylonites in granitoids has previously been attributed to reaction weakening associated with myrmekite formation and recrystallization of myrmekite (Ishii et al., 2007; Menegon et al., 2008) or to cataclastic deformation in feldspar (Pryer, 1993; Hippertt, 1998; Viegas et al., 2016). Unlike with myrmekite formation, development of banded mylonites by cataclastic deformation can occur solely through brittle breakdown of feldspar with no reaction products (Viegas et al., 2016). The deepest banded mylonite (sample CB13-55a, 430 m) records a transition from a coarse-grained feldspar network (Figs. 11A, 11B) to banded mylonite (Fig. 11E) at ∼470 °C and moderately low stress (32.1 ± 14 MPa), corresponding to a strain rate of ∼1 × 10‒13 s‒1. In this sample, K-feldspar porphyroclasts with asymmetric strain-induced myrmekite have recrystallized tails that deflect into feldspar-rich layers. This transition, in addition to the abundance of myrmekite within the CBSZ along K-feldspar porphyroclasts and also within feldspar C-bands (Figs. 7, 8D), suggests that myrmekite reaction, rather than cataclastic deformation, was the main process that contributed to layering of banded mylonites within the CBSZ.

Development of polyphase domains within banded mylonites (Fig. 11F; i.e., sample CB13-55c, 369 m) may be promoted by dilation and precipitation accompanied by GBS. Dissolution and precipitation in polyphase domains within this sample are inferred from K-feldspar filling dilatant sites between fine recrystallized quartz as well as from biotite and K-feldspar at quartz and plagioclase triple junctions. K-feldspar dissolution is reasonable in this scenario considering the abundance of myrmekite within the CBSZ (Fig. 7) and the excess K2O produced in the myrmekite-forming reaction (Equation 6; Ishii et al., 2007). This type of behavior and microstructure is consistent with disaggregation of monophase domains by dynamic recrystallization, dissolution and precipitation in dilatant sites opened by rotation, and GBS (Kilian et al., 2011; Platt, 2015). Quartz domains also support a possible role of GBS, as the CPOs show smeared a-axes with three dispersed maxima, consistent with the activity of GBS (Cross et al., 2017a). Polygonal subgrains consistent with SGR recrystallization in quartz with a strong y-maximum c-axis orientation indicate that dislocation creep with preferential prism <a> slip was also active, perhaps suggesting that quartz deformation at least locally was accommodated by a combination of GBS and dislocation creep, or DisGBS (i.e., Miranda et al., 2016).


  • 1. The Cordillera Blanca shear zone (CBSZ) preserves an ∼450-m-thick extensional shear zone characterized by undeformed granite in the core of the Cordillera Blanca batholith and increasingly mylonitized fabrics toward the western margin of the batholith. Brittle deformation occurs in the uppermost ∼150 m, as evidenced by pseudotachylyte, fault breccia, and cataclasite, and culminates in a moderate- to low-angle normal detachment fault.

  • 2. Deformation along the CBSZ was initiated at near-solidus conditions above ∼630 °C in weakly deformed granodiorite at the structurally lowest position relative to the discrete detachment fault, and continued to subsolidus conditions by fracturing of feldspar phenocrysts, dislocation creep, and neocrystallization of fine-grained feldspar along with growth of myrmekite. Initial interconnection of quartz domains occurred at temperatures near 500 °C under low stress (16.5 ± 13 MPa). Layering of mylonites at structurally deep positions occurred at 470 ± 50 °C where quartz recrystallized via grain-boundary migration and subgrain rotation processes during deformation by dislocation creep. Myrmekite growth contributed to the formation of banded mylonites with quartz monophase and quartz + feldspar ± mica polyphase domains. Dissolution and precipitation in dilatant sites may also have played a role in formation of layered mylonites and at least locally contributed, along with grain-boundary sliding, to formation of fine-grained, mixed quartz + K-feldspar + Na-feldspar + biotite domains.

  • 3. Quartz deformed predominantly by dislocation creep between 400 and 500 °C, producing strong crystallographic preferred orientations indicative of prism <a> slip and rhomb <a> slip. Variations in differential stress estimates despite similar deformation temperatures across the shear zone suggest that strain rate, rather than temperature, may have been more important in shear zone evolution for the CBSZ.

  • 4. Throughgoing brittle fabrics overprint ductile quartz-rich fabrics within the upper ∼150 m. Differential stresses from the CBSZ are consistent with estimates from granitoids in other extensional settings such as the Ruby Mountains and Whipple Mountains core complexes in the North American Cordillera. The highest differential stress in the CBSZ was recorded by the deepest sample within the shear zone to preserve a throughgoing brittle overprint of ductile fabrics, indicating that within the CBSZ, rocks have a strength of at least 71.5 ± 36 MPa at the ductile to brittle transition.


We thank Andrew Cross and John Platt for thorough and constructive reviews. Quartz reference materials were graciously provided by W. Nachlas. We thank J. Mauch for sampling support in the field, along with L. Kellerman at the Montana State University Imaging and Chemical Analysis Laboratory for technical assistance in acquiring EBSD data, and A. Patchen at the University of Tennessee for assistance in acquiring EPMA data. J. Langille assisted with drafting of the modified Giovanni et al. (2010) map. A. Cafferata of Caraz, Peru, and J. Olaza of Huaraz, Peru, aided in field expeditions. This work would not have been possible without funding support from the U.S. National Science Foundation through grant EAR-1220237 to M. Jessup and grant EAR-1220216 to C. Shaw. TitaniQ geothermometry was funded by the Sigma Xi Grant-in-Aid of Research program, awarded to C. Hughes.

1Supplemental Files. Table S1: Sample locations. Table S2: Mineral chemistry used for two-feldspar thermometry. Text S1: Cathodoluminescence and Titanium-in-Quartz (TitaniQ) electron probe microanalyzer methods. Text S2: Raw data processing methods for electron backscatter diffraction (EBSD) analyses. Text S3: Discussion of TitaniQ reliability for this study. Figure S1: Recrystallized quartz grain size distributions. Figure S2: EBSD-derived evidence for grain boundary sliding. Please visit https://doi.org/10.1130/GES02040.S1 or access the full-text article on www.gsapubs.org to view the Supplemental Files.
Science Editors: Raymond M. Russo, Andrea Hampel.
Associate Editor: Jose M. Hurtado.
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.