Fluid-rock interaction recorded in black fault rocks in the Kodiak accretionary complex, Alaska

Abstract

Ultrafine-grained black fault rocks (BFRs) in the Pasagshak Point Thrust of the Kodiak accretionary complex are examples of fault rocks that have recorded seismicity along an ancient subduction plate boundary. Trace element concentrations and ⁸⁷Sr/⁸⁶Sr ratios of BFRs and surrounding foliated/non-foliated cataclasites were measured to explore the nature of fluid-rock interactions along a subduction thrust. Foliated and non-foliated cataclasites do not show significant geochemical anomalies, suggesting that they were formed by slowly distributed shear. BFRs are characterized by Li and Sr enrichment, Rb and Cs depletion, and a low ⁸⁷Sr/⁸⁶Sr ratio. These geochemical signatures can be explained by fluid-rock interactions at >350°C, which result in preferential removal of Rb and Cs and formation of plagioclase under the presence of fluids with high Li and Sr concentrations and low ⁸⁷Sr/⁸⁶Sr ratios. Geochemical anomalies recorded by the BFRs indicate both frictional heating and external fluid influx into the subduction thrust.

Findings

Introduction

In subduction zones, fluids are abundant due to porosity reduction and dehydration of hydrous minerals. Most of the seismic moments released by global earthquakes occur along subduction zones, where fluids play important mechanical and chemical roles in seismic faulting. From a geological perspective, pseudotachylytes formed by frictional melting during high-velocity seismic slip have classically been considered as the only definitive evidence of earthquakes (Sibson 1975; Cowan 1999). Although occurrences of pseudotachylytes in crystalline continental rocks were known until approximately the 1990s (Swanson 1992), several pseudotachylyte-bearing fault zones in subduction settings have been discovered in recent years (Ikesawa et al. 2003; Kitamura et al. 2005; Rowe et al. 2005, 2011; Mukoyoshi et al. 2006; Okamoto et al. 2006; Meneghini et al. 2010; Hashimoto et al. 2012).

The ultrafine-grained black fault rock (BFR) in the Kodiak accretionary prism, Alaska, was first discovered by Rowe et al. (2005), who considered this rock as an example of subduction zone pseudotachylytes. Because BFR is unique in terms of its large thickness (up to 30 cm) and widespread occurrence, it has been further investigated on both microscopic (Meneghini et al. 2010) and macroscopic (Rowe et al. 2011) scales. Meneghini et al. (2010) also characterized the geochemistry of fault rocks and detected Sr anomalies. Here, we describe new geochemical observations that emphasize the role of fluids in the development of BFR.

Anomalies of trace element concentrations and isotope compositions within fault zones have recently been used to estimate maximum fluid temperatures (Ishikawa et al. 2008; Hirono et al. 2009; Hamada et al. 2011; Honda et al. 2011) because trace elements and isotopes can reflect fluid-rock interactions at temperatures as high as approximately 350°C. The mineralogy of fault zones can also affect element distribution during fluid-rock interactions (Yamaguchi et al. 2011b). Taking these works into consideration, we attempt to reconcile trace element concentrations and Sr isotope ratios in the ultrafine-grained BFR and surrounding cataclasites in Kodiak Island to gain a better understanding of the role of fluids in faulting. Based on the geochemical result, we interpret these geochemical results in terms of fluid-rock interactions recorded in plate boundary thrusts.

Geological setting and occurrences of cataclasites

The Kodiak archipelago in Alaska (Figure 1a) consists of a Jurassic to Eocene accretionary complex composed of NW-dipping, thrust-bounded units showing progressively younger ages toward the southeast (Moore et al. 1983; Plafker et al. 1994) that are intruded by 58- to 59-Ma near-trench granitic rocks reflecting ridge subduction (Farris et al. 2006; Ayuso et al. 2009). The Ghost Rocks Formation, consisting of the latest Cretaceous-earliest Paleocene ridge-flank deposits, is distributed over the southeastern part of Kodiak Island (Fisher and Byrne 1987). This formation is overthrust by the Late Cretaceous Kodiak Formation and is underthrust by Eocene marine sediments (Moore 1969). The structural base of the Ghost Rocks Formation consists of a tectonic mélange composed of blocks of coherent turbidites and sandstones in a matrix of disrupted mudstones and argillites (Byrne 1982, 1984). Pillow basalts and hyaloclastites showing prehnite-pumpellyite facies metamorphism are included in the matrix (Moore 1969; Byrne 1982).

Figure 1

Geological setting of the study area. (a) Ghost Rocks Formation in Kodiak Island. (b) Simplified geological map of the Pasagshak Point Thrust and outcrop localities.

The Pasagshak Point Thrust (Rowe et al. 2005, 2011; Meneghini et al. 2010) is a strike-parallel fossil thrust zone exposed along the rocky shoreline of Pasagshak Point (Figure 1). The thrust is located within the basal mélange zone of the Ghost Rocks Formation, implying shear localization along an ancient subduction thrust (Rowe et al. 2005). Both the hanging wall and footwall of the thrust consist of mélanges composed of terrigenous sandstones, mudstones, and argillites (Additional file 1: Figure S1a) with minor greenstones. Large sandstone blocks (thickness >10 m) occur in the hanging wall in places. S-C structures, slickenlines, and rotated fold axes within the fault zone suggest a top-to-southeast sense of shear (Rowe et al. 2005, 2011; Meneghini et al. 2010). The paleotemperature of the syn-tectonic fluid is estimated to be 230°C to 260°C from fluid inclusions in quartz veins within the mélange zone (Vrolijk et al. 1988), whereas the maximum burial temperature of the wall rock is estimated to be 250°C ± 10°C from vitrinite reflectance analysis of carbonaceous material within the argillite (Rowe et al. 2011). The similarity of the two temperature estimates suggests that deformations recorded in the thrust zone occurred at approximately 250°C at seismogenic depths of the plate boundary décollement.

Various types of cataclasites, 7 to 31 m in thickness, occur within the thrust zone (Additional file 1: Figure S1) (Rowe et al. 2011). Detailed structures of cataclastic fault rocks have been described in previous accounts (Rowe et al. 2005, 2011; Meneghini et al. 2010). Cataclastic fault rocks are classified into three categories: foliated cataclasites (both clast-rich and matrix-rich; Additional file 1: Figure S1b), non-foliated cataclasites (Additional file 1: Figure S1c), and ultrafine-grained black fault rocks (BFRs; Additional file 1: Figure S1d). The BFRs, locally thicker than 30 cm, make up the principal shear zones in the thrust zone and are further classified into three types: aphanitic, grain-supported, and clast-bearing. The structural features of foliated and non-foliated cataclasites in comparison with the textures formed by experimentation suggest that foliated cataclasites were formed by slow creep (Niemeijer and Spiers 2006), whereas non-foliated cataclasites were formed at faster shear rates, such as through afterslips or delocalized slips associated with earthquakes (Rowe et al. 2011). The aphanitic BFRs contain crystalline microlayers that are regarded as melt-origin microtextures (Meneghini et al. 2010) and are therefore interpreted as pseudotachylytes formed during seismic slip. Conversely, grain-supported BFRs and clast-bearing BFRs are considered to be reworked pseudotachylytes that have mixed with surrounding cataclasites, respectively (Meneghini et al. 2010). BFRs show the greatest biomarker thermal maturity, whereas that of the cataclasites decreases with distance from the BFRs (Savage et al. 2014). BFRs, foliated cataclasites, and non-foliated cataclasites exhibit mutually crosscutting relationships (Additional file 1: Figure S1e) and represent repeated deformations with different strain rates, which may be related to seismic cycles (Rowe et al. 2011).

Geochemical analysis and results

Samples of two host rocks, six foliated cataclasites, two non-foliated cataclasites, two clast-bearing BFRs, and one aphanitic BFR were selected for geochemical analysis (Additional file 2: Table S1). These samples were collected from three representative outcrops (wpt 009, wpt 014, and wpt 015 of Meneghini et al. 2010; Rowe et al. 2011) and were carefully cut into approximately 10 mm × 10 mm × 10 mm chips to avoid the influence of surface weathering. The chips were ultrasonicated in deionized water and were then pulverized in an agate ball mill.

Major element concentrations were determined by an X-ray fluorescence spectrometer (XRF; MagiX PRO, Spectris, Egham, Surrey, UK) installed at the Kochi Core Center (KCC), Kochi, Japan, using the conventional glass bead method. Powdered splits for trace elements and isotope analyses were dissolved by HNO₃-HF-HClO₄ digestion and were then evaporated. The residue was dissolved with HNO₃ and was diluted to approximately 1:1,000. Concentrations of Sc, V, Cr, Co, Ni, Cu, Zn, As, Rb, Sr, Y, Zr, Nb, Mo, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta, Pb, Th, and U were determined by an inductively coupled plasma mass spectrometer (ICP-MS; Agilent 7500c, Agilent Technologies, Santa Clara, CA, USA) installed at the Department of Systems Innovation, the University of Tokyo, Tokyo, Japan (Kato et al. 2005). Li, Be, Ti, W, and Tl concentrations were measured by ICP-MS (ELAN DRC II, PerkinElmer, Waltham, MA, USA) at the KCC. Sr was separated using ion-exchange resin columns (Sr Resin, Eichrom Technologies Inc., Lisle, IL, USA), and Sr isotope ratios were determined by thermal ionization mass spectrometry (TIMS; TRITON, Thermo Fisher Scientific, Waltham, MA, USA) at the KCC. The ⁸⁷Sr/⁸⁶Sr value obtained for the NIST SRM 987 standard was 0.7102570 ± 0.000010.

Analytical results are shown in Additional file 2: Table S1. Major element and some trace element compositions are consistent with previously reported values as determined by XRF (Meneghini et al. 2010). Element compositions normalized by averaged host rock values are presented in Figure 2. The clast-bearing BFRs and aphanitic BFRs show systematic positive or negative anomalies in some elements such as Li and Sr enrichments and Rb, Cs, and As depletions. Conversely, foliated and non-foliated cataclasites do not show such clear systematic patterns. Sr concentrations in the BFRs are generally high, ranging from 228 to 941 ppm, with the highest concentration exhibited in sample 015-10BS (aphanitic BFR). ⁸⁷Sr/⁸⁶Sr ratios are 0.70877 and 0.70922 in the host rocks, 0.70667 to 0.70875 in the foliated and non-foliated cataclasites, and 0.70593 to 0.70627 in the BFRs. The aphanitic BFR exhibits the lowest ⁸⁷Sr/⁸⁶Sr ratio (0.70593).

Figure 2

Trace element anomalies of fault rocks normalized by averaged host rock compositions.

Discussion

The trace element composition of the BFRs exhibits tendencies similar to those of principal shear zones of other fault zones reported to possess coseismic fluid-rock interactions (Ishikawa et al. 2008; Hamada et al. 2011; Honda et al. 2011). Enriched and depleted elements in the Kodiak BFR (Li, Rb, Cs, and Sr) are known to be fluidly mobile elements easily affected by high-temperature fluid-rock interactions. The trace element concentration in a solid chemically equilibrated with fluid was modeled by Ishikawa et al. (2008) as

$${\mathit{C}}_{\mathrm{s}}=\mathit{D}\left({\mathit{C}}_{\mathrm{s}0}+\mathit{R}{\mathit{C}}_{\mathrm{f}0}\right)/\left(\mathit{R}+\mathit{D}\right),$$

where C_s0, C_s, and C_f0 are the element concentrations in the initial solid, solid equilibrated with fluid, and initial fluid, respectively; R is the fluid/solid mass ratio; and D is the temperature-dependent solid/fluid bulk distribution coefficient of the element. In our model analysis, the values for D were estimated by hydrothermal experimental results at 250°C, 300°C, and 350°C using hemipelagic sediments (You et al. 1996; Ishikawa et al. 2008). The values for C_f0 and C_s0 are based on the deep pore fluid data from offshore Texas (Kharaka and Hanor 2003) reported by Ishikawa et al. (2008) and averaged host rock, respectively. Figure 3 compares actual BFR composition (red and orange squares and triangles) and calculated fault rock compositions (solid and dashed lines, depending on R values) of selected trace elements (Rb, Sr, Cs, Ba, La, Sm, and Pb) at temperatures of 250°C, 300°C, and 350°C. The calculations at T = 350°C and R = 100 are closest to the actual trace element compositions of the aphanitic BFR. For the clast-bearing BFR, conditions of T = 350°C and R = 10 or T = 300°C and R = 100 are the most appropriate for reproducing trace element compositions, except for Li. The estimated R values (10 to 100) are significantly higher than those estimated in other fault zones that range from 3 to 7 (Ishikawa et al. 2008; Hamada et al. 2011; Honda et al. 2011). The reasons for this incongruity may be the lack of experimental data for temperatures higher than 350°C and that R values of 10 to 100 and T values of 300°C to 350°C were the maximum and minimum estimations. Alternatively, these results can be explained by a longer fluid residence time and greater duration of high temperatures due to the depth of the Pasagshak Point Thrust. In both cases, our geochemical datasets suggest that the overall composition of the BFR is controlled by fluid-rock interactions at temperatures higher than 350°C. This temperature is at least 100°C higher than the approximately 250°C estimated background temperature of the Pasagshak Point Thrust determined by vitrinite reflectance (Rowe et al. 2011) and is only achievable under conditions of coseismic frictional heating. High temperatures recorded in the BFR are consistent with recent results of biomarker thermal maturity (Savage et al. 2014). In contrast to BFR, foliated and non-foliated cataclasites exhibit no high-temperature signals. This is also consistent with the structural interpretation that foliated and non-foliated cataclasites are formed by slow creep and intermediate strain rate deformation, respectively (Rowe et al. 2011).

Figure 3

Comparison of actual BFR composition and calculated fault rock compositions of selected trace elements. Trace element anomalies of black fault rocks (BFRs; red diamonds and orange squares and triangles) and calculated fault rock compositions (solid and dashed lines, depending on R values) of selected trace elements (Rb, Sr, Cs, Ba, La, Sm, and Pb) at temperatures of 250°C, 300°C, and 350°C. The vertical axis shows host rock-normalized relative deviation.

Sr enrichment has been reported in the principal gouge zones of the Taiwan Chelungpu Fault, where it has been regarded as an indicator of plagioclase enrichment (Ishikawa et al. 2008). Ishikawa et al. (2008) concluded that plagioclase was hydrothermally precipitated from high-temperature fluid heated during coseismic frictional slip. Meneghini et al. (2010) reported that the tabular, euhedral, and zoned plagioclase crystals contained in the Pasagshak Point Thrust are enriched in aphanitic BFR. Particularly, in the crystalline microlayers, Ca-rich plagioclase is found only in the BFR. The Sr enrichment observed in the BFR would reflect the enrichment of plagioclase, as suggested by Meneghini et al. (2010). Based on this interpretation, Sr enrichment is caused by plagioclase microlites crystallized from fluid-rich frictional melt. Such an open-system fluid circulation during the formation of pseudotachylytes is also described in an Rb-Sr geochronometric study of pseudotachylytes from the North Cascade Mountains (Magloughlin 2003). An alternative possibility is that the Sr enrichment of the BFR was caused by hydrothermal precipitation of plagioclase from high-temperature fluid, as suggested by Ishikawa et al. (2008). Further detailed layer-by-layer microchemical analyses or Sr chemical mapping is needed to clarify the origin and distribution of Sr within the BFR because aphanitic BFR is composed of a micrometer-millimeter-scale ultrafine layering of granular and crystalline microlayers, which is much finer than the 1-cm³ scale of our sampling.

The ⁸⁷Sr/⁸⁶Sr ratio decreases with an increase in Sr concentration (Figure 4a). This reverse trend implies that the additional Sr incorporated into rocks with high Sr concentrations had a lower ⁸⁷Sr/⁸⁶Sr ratio (<0.7059) than the host rock, which has ⁸⁷Sr/⁸⁶Sr ratios of 0.7088 to 0.7092. If we assume that the incorporated Sr originated from pore fluid that existed during the formation of fault rocks, the fluid would be characterized by a high Sr concentration and low, less radiogenic ⁸⁷Sr/⁸⁶Sr ratio. The observed substantially lower ⁸⁷Sr/⁸⁶Sr ratio (<0.7059) cannot be explained by seawater or terrigenous sediments. Possible candidates for the source of this lower ⁸⁷Sr/⁸⁶Sr ratio include hydrated oceanic crust such as the approximately 62-Ma syn-sedimentary mid-ocean ridge basalts (MORBs) in the Ghost Rocks Formation related to ridge subduction (Moore et al. 1983) and subducted oceanic crust. Yamaguchi et al. (2012) detected calcite veins which precipitated from fluids with relatively low ⁸⁷Sr/⁸⁶Sr and large δ¹⁸O from a paleosubduction thrust in the Shimanto accretionary complex and suggested that fluids dehydrated from the oceanic crust might have acted as the sources of fluids along subduction thrusts. Recently, Sano et al. (2014) reported mantle-derived fluid from seawater at the bottom of the Japan Trench approximately 1 month after the 2011 Tohoku earthquake. Such co- or post-seismic discharges of deep-origin fluid are presumed to have the capability of producing the observed compositional and isotopic anomalies in the BFR.

Figure 4

⁸⁷ Sr/ ⁸⁶ Sr ratios of host and fault rocks plotted against (a) 1/Sr and (b) Li concentrations.

We observed a negative correlation between the ⁸⁷Sr/⁸⁶Sr ratio and Li concentration (Figure 4b). The Li concentrations of the fault rocks are significantly higher than their simulated values (Figure 3). Li usually shows depletion in the principal shear zones of faults (Ishikawa et al. 2008; Hamada et al. 2011; Honda et al. 2011) because its D value is small at high temperatures. The high concentration of Li in the fault rocks therefore suggests extremely high concentrations of Li (approximately 10 to 100 ppm) in the pore fluid that are 10² to 10³ times more concentrated than that of modern seawater (0.18 ppm; Wunder et al. 2006). The value of D_Li at the background temperature of 250°C is low enough to preferentially mobilize Li relative to Rb and Cs; D_Li, D_Rb, and D_Cs values at 250°C are 11, 290, and 80, respectively. The pore fluid that reacted with the BFR could have been significantly enriched in Li through Li extraction from the host rocks along the fluid pathways. High Li concentration in fluid along a seismogenic plate boundary thrust would explain the Li-enriched pore water up to 73 ppm discovered at the décollement zone of the Ocean Drilling Program Site 808 near the deformation front of the Nankai Trough (You et al. 1995). Direct measurement of fluid chemistry such as that through fluid inclusion analysis is a possible approach for understanding the geochemical characteristics of fluids existing in the deep portions of plate boundary thrusts.

The Kodiak complex is comparable to the Shimanto accretionary complex in southwest Japan: both complexes were formed by subduction of a relatively young oceanic plate, and they have similar lithologies characterized by thick terrigenous sediments with rare pelagic sediments. However, the occurrences of fault rock types and fluid-rock interaction patterns differ significantly, particularly considering their similar background temperatures (approximately 250°C). In the Shimanto complex, abundant quartz-calcite-ankerite veins are commonly observed along pseudotachylyte- and ultracataclasite-bearing large thrusts (Okitsu mélange: Ikesawa et al. 2003; Mugi mélange: Yamaguchi et al. 2012; Nobeoka Thrust: Okamoto et al. 2006; Yamaguchi et al. 2011a), whereas in the Pasagshak Point Thrust, mineral veins are relatively rare, except for minor quartz shear veins. Instead of mineral veins, plagioclase enrichment may reflect coseismic high-temperature fluid-rock interactions as recognized within the BFR. Such differences in mineralogy and precipitates may reflect variations in regional-scale fluid flow patterns and physicochemical characteristics of fluids along subduction plate boundary interfaces.

Conclusion

Trace element concentrations and ⁸⁷Sr/⁸⁶Sr ratios of fault rocks in the Pasagshak Point Thrust reflect fluid-rock interactions related to coseismic slip along an ancient subduction thrust. Sampled BFRs are characterized by Li and Sr enrichment, Rb and Cs depletions, and low ⁸⁷Sr/⁸⁶Sr ratios. Such geochemical signatures can be explained by plagioclase formation at temperatures higher than 350°C or the presence of a fluid with high Li and Sr concentrations and low ⁸⁷Sr/⁸⁶Sr ratios. The geochemical anomalies of the BFRs therefore suggest both frictional heating and influx of external fluid in a subduction thrust. In contrast to the BFRs, foliated and non-foliated cataclasites do not exhibit significant geochemical anomalies, suggesting that they were formed by slow deformation without large frictional heat or large-scale fluid flow. Further geochemical and comparative investigations are necessary to explore the nature of fluids along seismogenic subduction plate boundaries.