Abstract

The San Andreas Fault Observatory at Depth (SAFOD) drill hole encountered indurated, high-seismic-velocity arkosic sedimentary rocks west of the active trace of the San Andreas fault in central California. The arkosic rocks are juxtaposed against granitic rocks of the Salinian block to the southwest and against fine-grained Great Valley Group and Jurassic Franciscan rocks to the northeast. We identify three distinct lithologic units using cuttings, core petrography, electrical resistivity image logs, zircon fission-track analyses, and borehole-based geophysical logs. The upper arkose occurs from 1920 to 2530 m measured depth (mmd) in the borehole and is composed of five structural blocks defined by bedding orientations, wireline log character, physical properties, and lithologic characteristics. A clay-rich zone between 2530 and 2680 mmd is characterized by low Vp and an enlarged borehole. The lower arkose lies between 2680 and 3150 mmd. Fission-track detrital zircon cooling ages are between 64 and 70 Ma, appear to belong to a single population, and indicate a latest Cretaceous to Paleogene maximum depositional age. We interpret these Paleocene–Eocene strata to have been deposited in a proximal submarine fan setting shed from a Salinian source block, and they correlate with units to the southeast, along the western and southern edge of the San Joaquin Basin, and with arkosic conglomerates to the northwest. The arkosic section constitutes a deformed fault-bounded block between the modern strand of the San Andreas fault to the northeast and the Buzzard Canyon fault to the southwest. Significant amounts of slip appear to have been accommodated on both strands of the fault at this latitude.

INTRODUCTION

The primary component of the San Andreas Fault Observatory at Depth (SAFOD) is a 3-km-deep deviated borehole that crosses the San Andreas fault zone in central California (Hickman et al., 2004, 2005) (Fig. 1). Numerous questions about fault behavior, structure, composition, physical properties, and fluid-rock interactions in fault zones are being addressed in this project (i.e., Boness and Zoback, 2006; Solum et al., 2006; Erzinger et al., 2006; Hole et al., 2006; Schleicher et al., 2006; Bradbury et al., 2007).

Figure 1.

Geologic map of the SAFOD site, based on a compilation and modification of the maps of Walrond and Gribi (1963), Dickinson (1966), Dibblee (1971), Durham (1974), Sims (1990), and Thayer and Arrowsmith (2005). Figure modified from Bradbury et al. (2007). Inset map shows the location of the area relative to California.

Figure 1.

Geologic map of the SAFOD site, based on a compilation and modification of the maps of Walrond and Gribi (1963), Dickinson (1966), Dibblee (1971), Durham (1974), Sims (1990), and Thayer and Arrowsmith (2005). Figure modified from Bradbury et al. (2007). Inset map shows the location of the area relative to California.

A comprehensive suite of geophysical data and geological studies (see Geophysical Research Letters special volume 31, numbers 12 and 15, 2004; Dibblee, 1971; Sims, 1990; Page et al., 1998; Thayer and Arrowsmith, 2005; M. Rymer, 2005, personal commun.) was used to develop a subsurface model that predicted the borehole would encounter Salinian granitic rocks on the southwestern side juxtaposed against a Franciscan complex northeast of a nearly vertical San Andreas fault (Hickman et al., 2004; McPhee et al., 2004). A vertical pilot hole drilled in 2002 went through a 1 km thick Cenozoic sedimentary section overlying Salinian granite, consistent with the geologic model. In 2004 and 2005, phases 1 and 2 of the main SAFOD hole were drilled, and at 1920 m measured depth (mmd—the distance measured along the borehole), the borehole encountered an abrupt change from Salinian granodiorite into arkosic sedimentary rocks where the borehole is deviated northeast (Bradbury et al., 2007). This section of arkose is present from 1920 mmd to 3157 mmd along the track of the borehole. At 3157 mmd the borehole encountered fine-grained sandstones and mudstones, and at 3360 mmd it encountered rocks that are now identified as Cretaceous Great Valley Group (Evans et al., 2005) based on their texture, composition, and the presence of Late Cretaceous Great Valley fossils (K. McDougall, 2006, written commun.). The rocks between 3157 and 3410 mmd appear to be a complex zone of low-seismic-velocity rocks that are the subject of other studies (Hickman et al., 2008).

The presence of arkosic sedimentary rocks on the southwestern side of the fault and the presence of Great Valley sedimentary rocks northeast of the fault, and the apparent juxtaposition of Salinian granitic rocks against the arkosic rocks, pose several geological and tectonic questions regarding fault zone structure and the tectonic evolution of the San Andreas fault in this area. These questions include the following: What is the origin of the arkosic sedimentary rocks—are they a part of the Salinian block? Have these rocks been transported a long distance along older strands of the San Andreas fault? The offset of the Salinian block is distributed over several different parallel strike-slip fault systems (Graham, 1978; Dickinson and Butler, 1998; Whidden et al., 1998; Dickinson et al., 2005), which accounts for the ~500 km offset from the Sierra Nevada plutonic complex. In the vicinity of the SAFOD site the Buzzard Canyon fault is subparallel to the San Andreas fault, having an unknown amount of slip (Rymer et al., 2003; Thayer and Arrowsmith, 2005). The Buzzard Canyon fault exposed at the surface has been correlated to the fault that separates the granodiorite from the arkosic rocks in the SAFOD borehole (Bradbury et al., 2007). Arrowsmith (2007) suggests that several faults lie southwest of the San Andreas fault, and these faults may have accommodated a significant amount of the total slip across the San Andreas fault plate boundary in central California.

This paper provides the first detailed lithologic characterization of the arkosic rocks as a first step toward understanding their geologic history and answering the questions posed above. The composition of the arkosic rocks is evaluated to interpret provenance, depositional environment, and diagenetic history. The observations of composition are then used to identify surface units that may be equivalent to the SAFOD arkose section. These possible correlations constrain the geometry and evolution of the San Andreas fault in the area of the SAFOD site.

Methods

We use five data sets to characterize and interpret the arkosic sedimentary rocks: drill cuttings, spot and sidewall core, borehole geophysical logs, electrical image logs, and zircon fission-track analyses. Drill cuttings are silt- and sand-sized rock fragments created by rotary drill bit action on the rock at the bottom of the hole that circulate to the surface in the drilling mud, where they emerge as a mix of engineered drilling mud and rock fragments. Cuttings composition was examined as a function of depth along the borehole to identify major changes in wall rock composition. The spot cores range from centimeters to meters long, are 10–15 cm in diameter, and were collected using a rotary coring tool bit at the end of the drilling string. Side-wall samples are 5–8 cm long, 2 cm diameter cores that were collected with a percussive side-wall core tool that shoots a hollow bullet-style core tube laterally into the borehole wall. Borehole geophysical logs measure physical properties including velocity, density, and porosity. Resistivity-based image logs reveal the nature of lithology, bedding, faults, and fractures by measuring the resistivity of the borehole wall. High-quality data enable us to determine the bedding orientation and character, formation thickness, sedimentary structures, and fault, fracture, or vein orientations and character, and general lithologic characteristics such as grain size, shape, and sorting (Thompson, 2000). Dynamically normalized electrical image logs were collected in the arkosic section from 1920 mmd to 3050 mmd with the Schlumberger Formation Microresistivity Sonde (FMS) tool (Ekstrom et al., 1986). The Geomechanics International (GMI) Imager™ software was used to examine the electrical image logs and to plot bedding and fracture orientations.

Drill cuttings provide a nearly continuous sampling of the rocks encountered in drilling. Care must be taken when using cuttings to infer lithology or structure in the borehole. Contamination can occur as drilling fluids pick up and transport cuttings. This can be especially problematic in the deviated portion of the hole, where a bed of cuttings may accumulate in the out-of-gauge sections of the hole. The stability of the SAFOD main hole varied through the deviated portion of the hole; breakouts and caving were common, as seen in the caliper log (Fig. 2). These problems with borehole integrity may cause cuttings from higher in the hole to fall and mix with the cuttings lower in the hole. At the drill site, cuttings were gently washed with water through small-mesh sieves to eliminate coatings of drilling mud. Some of the clay-sized grains from the cuttings may have been washed away, possibly skewing the compositional analysis by an unquantifiable amount.

Figure 2.

Compiled wireline logs (caliper, Vp, gamma ray, density, neutron porosity) plotted along the borehole trajectory, along with composition plot of clay grains versus iron oxide–rich grains and fracture density histograms in 10 m wide bins. Every grid line in the log corresponds to the approximate depth of a cuttings thin section analyzed petrographically for this study. Intervals in velocity log that drop to zero are intervals where the tool malfunctioned and no data were collected. The likely locations of faults are indicated by red arrows, and black arrows indicate the location of block boundaries. Shadings of the figure represent the structural blocks defined by lithology and orientation of bedding imaged in the borehole.

Figure 2.

Compiled wireline logs (caliper, Vp, gamma ray, density, neutron porosity) plotted along the borehole trajectory, along with composition plot of clay grains versus iron oxide–rich grains and fracture density histograms in 10 m wide bins. Every grid line in the log corresponds to the approximate depth of a cuttings thin section analyzed petrographically for this study. Intervals in velocity log that drop to zero are intervals where the tool malfunctioned and no data were collected. The likely locations of faults are indicated by red arrows, and black arrows indicate the location of block boundaries. Shadings of the figure represent the structural blocks defined by lithology and orientation of bedding imaged in the borehole.

Within these limitations, compositional analysis of cuttings is effective in characterizing subsurface lithologic types (Winter et al., 2002), especially when used in conjunction with other methods (Solum et al., 2006; Bradbury et al., 2007). We examined cuttings and data collected in 2004 and 2005. We used optical microscopic methods to systematically quantify the mineralogy, deformation textures, and the degree and nature of alteration in the rocks. We used a riffle splitter to split an ~1 cm3 subsample of cuttings from every 15.3 mmd of the archived cuttings samples. Magnetic separation of the dry sub-sample removed shards of drill bit material. The magnetic material was examined in a binocular microscope, and magnetic rock fragments were integrated back into the sample to maintain the representative nature of the subsample. After magnetic separation, samples were lightly washed and decanted in distilled water to remove drilling mud additives such as walnut shells. These grains were used to make thin sections of grain mounts. A total of 83 thin sections were analyzed for this study. We used a modified Gazzi-Dickinson petrographic point-count method (Dickinson, 1970) to quantify the composition of the cuttings and core. The grain classification was modified to account for the limitations of the sampling procedures and the nature of cuttings. A grain may be classified both for its composition and for its texture, resulting in point counts that may sum to greater than 100%. Previous work on these rocks used the Gazzi-Dickinson method to count the bulk composition of the entire borehole, quantifying the ratio of quartz to feldspar to lithics, and dividing the lithic grains into several categories (Bradbury et al., 2007). We tailored this method and applied it to the arkosic section in order to deduce the diagenetic and deformational history as well as provenance of these rocks. For more detail regarding the standard compositional analyses of these rocks, see Bradbury et al. (2007). The raw data are presented in the GSA Data Repository (Appendix A).1

We counted grains of quartz, unaltered feldspar, opaque minerals, and accessory detrital minerals such as muscovite, fuchsite (a chromium-bearing muscovite), amphiboles, and pyroxenes to quantify the composition of the cuttings. We observed and recorded diagenetic features that include cements or matrix (these are grouped as it is difficult to differentiate between cement and matrix in cuttings), secondary quartz overgrowth, vein fragments, calcite, clay minerals, zeolites, and iron oxides that can also be considered an alteration feature derived from the alteration of mafic minerals. Alteration and deformational features such as cataclasite, altered feldspar, and biotite replaced by chlorite were recorded. The only clay minerals we counted as diagenetic features are those that show textural evidence of being either matrix or cement. This division of clay minerals is different from that of Solum et al. (2006), who used X-ray diffraction methods to quantify clay minerals in the subsidiary fault zones and the country rock of the SAFOD main hole. Iron oxides include free hematite grains and fine-grained material that is red in reflected light and in plane-polarized and cross-polarized transmitted light. Altered feldspar grains fall into three categories. The most common is the sericitized and vacuolated phase of alteration. Within the scope of the cuttings petrography, we define the alteration product of the first category to be illitic material, a phrase Moore and Reynolds (1989) use to encompass illite, mixed-layer phases, and other minerals such as sericite. The second altered feldspar category includes feldspars that show evidence of being altered directly to muscovite. In some cases the alteration products are coarse-grained clays or mica. The third category of alteration consists of feldspars that have been replaced almost completely by calcite; the grains still display feldspar twinning but with characteristic calcite birefringence. Damage characteristics are confined to cataclasites in the arkosic rocks and are discussed in a previous paper (Bradbury et al., 2007).

Twelve meters of nonoriented 10.16 cm diameter core was recovered from 3055 mmd to 3067 mmd (~2.8 km vertical depth) in the arkosic section (Fig. 3). Two 1.9 cm diameter, 3 cm long percussive sidewall cores were collected at 3084.6 mmd and 3121.2 mmd in the arkosic section. Fourteen petrographic thin sections were prepared from representative intervals in the spot core, and two thin sections were taken from the sidewall cores. Using thin sections we also evaluated evidence for crosscutting relationships and diagenetic features such as cementation, compaction structures, and mineral authigenesis.

Figure 3.

Photographs of the phase 1 core, runs 4 and 5 (runs 1–3 are in the granodiorite farther up in the borehole). Core is indexed by box number, smaller box numbers correspond to shallower depths. Core is 10.16 cm in diameter and each box is ~1 m long. Core grades from coarse-grained sandstones and matrix-supported granule-pebble conglomerate (boxes 10–16) to fine-grained sandstone (box 17) to siltstone (boxes 18–20) to very fine-grained siltstone/mudstone (boxes 21 and 22) and into a clay-rich shear zone (box 23). The small piece of rock to the right labeled “Run 5” is from a failed fifth coring attempt. Sample comes from the “core catcher,” which is a basket attached to the coring bit. That sample is a fine-grained iron oxide–stained sandstone. Red stars indicate where a sample was made into a thin section for our study (the samples in boxes 12, 14, 17, 20, and 22 were made into multiple thin sections in order to capture interesting features of the sample). The nomenclature used for labeling samples is of the form MHP1-B11, as in “main hole phase 1 box 11”. Here the stars are labeled with a condensed form of the sample number, i.e., B11.

Figure 3.

Photographs of the phase 1 core, runs 4 and 5 (runs 1–3 are in the granodiorite farther up in the borehole). Core is indexed by box number, smaller box numbers correspond to shallower depths. Core is 10.16 cm in diameter and each box is ~1 m long. Core grades from coarse-grained sandstones and matrix-supported granule-pebble conglomerate (boxes 10–16) to fine-grained sandstone (box 17) to siltstone (boxes 18–20) to very fine-grained siltstone/mudstone (boxes 21 and 22) and into a clay-rich shear zone (box 23). The small piece of rock to the right labeled “Run 5” is from a failed fifth coring attempt. Sample comes from the “core catcher,” which is a basket attached to the coring bit. That sample is a fine-grained iron oxide–stained sandstone. Red stars indicate where a sample was made into a thin section for our study (the samples in boxes 12, 14, 17, 20, and 22 were made into multiple thin sections in order to capture interesting features of the sample). The nomenclature used for labeling samples is of the form MHP1-B11, as in “main hole phase 1 box 11”. Here the stars are labeled with a condensed form of the sample number, i.e., B11.

A comprehensive suite of borehole-based, industry-standard geophysical logs was collected at the SAFOD site during all phases of drilling (see http://www.icdp-online.org/contenido/icdp/front_content.php?idcat=712). The primary logs we used for this study are P-wave velocities (Vp) derived from the sonic monopole velocity log, dipole S-wave velocity (Vs), resistivity, neutron porosity, density, caliper, gamma ray, and spectral gamma ray (Fig. 2). These logs were used to constrain the composition of the rocks, to determine bedding orientations and the orientation of fractures, and to evaluate regions where alteration products may be present. Raw data for the image log for unit thicknesses, bedding orientations, and fracture orientations and density are in the GSA Data Repository (Appendices B–D). A Hostile Environment Gamma Ray Neutron Sonde (HNGS) tool resolves the gamma ray signature into potassium, uranium, and thorium concentrations. These data are used to determine which elements have the dominant influence in the overall gamma ray signature.

Zircon fission-track analyses were performed on six samples from depths of 3121–3375 mmd using the method outlined in Bernet and Garver (2005). Zircons were separated from the arkosic sediment cuttings and analyzed for their fission-track dates. These dates provide some constraints on the age of sediment deposition. Four samples are from the coarse-grained sandstones and pebble conglomerates west of the San Andreas fault zone, and two are from fine-grained rocks in the region termed the “low-velocity zone” (Hickman et al., 2008) within the fault zone. Fifteen zircon grains per sample were dated. The annealing temperature of zircon is ~240 °C (Hurford and Carter, 1991), and in the case of sedimentary rocks, we do not date the time at which the sedimentary rocks cooled past 240 °C but rather when the source rock cooled through the annealing temperature, and as such these data are related to the provenance of the strata and not the thermal history of this part of the SAFOD hole.

RESULTS

Three main lithologic units have been identified in our analysis of the arkosic section. These units are subdivided into 11 structural domains (Fig. 2) based on bedding orientations determined from the analysis of the electrical image log data. The three main lithologic zones are referred to as (1) the “upper arkose,” which extends from 1920 to 2530 mmd, (2) the “clay-rich zone” from 2530 to 2680 mmd, and (3) the “lower arkose” from 2680 to 3157 mmd. The lithology and structure of each of these zones is provided below.

Upper Arkose

The upper arkose was encountered between 1920 and 2530 mmd in the borehole. Its porosity ranges from 0.01 to 0.44, with an average of 0.10, or 10%. The density is 2.13–2.72 g/cm3 (average of 2.56 g/cm3), and values of Vp range from 3.2 to 5.76 km/s, with an average of 4.57 km/s (Table 1). Sonic velocities generally increase with depth, although several low-velocity zones were encountered through the middle portions of the upper arkose.

TABLE 1.

RANGE AND MEAN VALUES FOR THE WIRELINE LOG DATA FROM THE ARKOSIC ROCKS AT SAFOD

Point-count analysis reveals that the rocks of the upper arkose are quartz and feldspar rich and contain 1%–10% clay-sized minerals as cement or matrix (Fig. 4). The layer silicate is typically chlorite or kaolinite. The seemingly anomalous gamma ray signatures averaging 129.08 API are typically associated with fine-grained rocks; however, the arkosic nature of the rocks makes the base level of the gamma ray log signature inherently higher. The gamma ray signature in the upper arkose is reasonably consistent through the section, with shallow peaks and troughs. This gamma ray pattern is interpreted to reflect relative homogeneity of the upper arkose. Quartz grains are generally undeformed aside from brittle fracturing, some of which is interpreted to have occurred as a result of drilling (Fig. 5A). Undulatory extinction and other plastic deformation of the quartz crystals are rarely observed in cuttings. Discrete quartz grains are the largest component of the unit, ranging from 45% to 60% (Fig. 4). Feldspar concentrations are higher in the upper arkose than in the other zones, constituting 4%–10% of the sample, showing approximately equal proportions of plagioclase to K-feldspar. The K-feldspar component is primarily microcline that contains classic tartan twinning. Orthoclase is present in small amounts.

Figure 4.

Composition of arkosic section derived from petrographic point counts, expressed in modal percentages. Graph is divided into the three different lithologic units: upper arkose, clay-rich zone, and lower arkose.

Figure 4.

Composition of arkosic section derived from petrographic point counts, expressed in modal percentages. Graph is divided into the three different lithologic units: upper arkose, clay-rich zone, and lower arkose.

Figure 5.

Examples of features commonly observed in cuttings. Photomicrographs are all taken with cross-polarized light except photo D, which is taken with plane-polarized light. Photos B, D, and E taken with a blue filter. A: Conchoidally fractured quartz grains. Fractured texture is drilling-induced. B: Clay-rich grain from the lower arkose, with several compositions of clay mineral present in the grain. Greenish color is due to the presence of fine-grained chlorite. C: Fuchsite in cross-polarized light. Interference colors are higher-order than in chlorite. D: The same fuchsite grain in plane-polarized light. Notice distinct green color. E: Plastically deformed quartz grain. F: Devitrified volcanic grain from the lower arkose. Photo light levels have been altered digitally in order to see features in the grain. Petrographically, volcanic grains are typically very dark.

Figure 5.

Examples of features commonly observed in cuttings. Photomicrographs are all taken with cross-polarized light except photo D, which is taken with plane-polarized light. Photos B, D, and E taken with a blue filter. A: Conchoidally fractured quartz grains. Fractured texture is drilling-induced. B: Clay-rich grain from the lower arkose, with several compositions of clay mineral present in the grain. Greenish color is due to the presence of fine-grained chlorite. C: Fuchsite in cross-polarized light. Interference colors are higher-order than in chlorite. D: The same fuchsite grain in plane-polarized light. Notice distinct green color. E: Plastically deformed quartz grain. F: Devitrified volcanic grain from the lower arkose. Photo light levels have been altered digitally in order to see features in the grain. Petrographically, volcanic grains are typically very dark.

The amount of iron oxide cements and staining ranges from 5% to 30% (Figs. 4 and 5). Some iron oxides occur as fracture fill or grain coatings, but the majority of iron oxide present forms cements and matrix of the grains (Fig. 6). Detrital micas such as muscovite, fuchsite, and biotite and Fe-Mg-bearing minerals are never present in an amount more than 1% of the sample (Fig. 7). Rock fragments other than granitic clasts are rarely observed. Altered feldspars make up 12%–20% of the sample (Fig. 8). Feldspar alteration ranges from slight with small amounts of illite to high where the feldspar is almost indistinguishable. The majority of altered feldspars contain sericite and clay minerals.

Figure 6.

Photomicrographs showing characteristic features from cuttings samples. Depths of the samples (meters) are noted at the lower right-hand corner of the photographs. All photographs were taken in cross-polarized light. Differences in color are due to different light levels. Scale bar applies to all photomicrographs. A: High-birefringence clay minerals fill space between quartz crystals. B: Altered plagioclase in center of photo shows alteration typical in the upper arkose. The grain is still easily recognizable as plagioclase with small areas of sericite formation. C: Grains displaying iron oxide–rich matrix/cement. The grain in the upper right of the photo shows deposition-related iron oxides (matrix), while the iron oxides in the grain in the lower left look postdepositional, either staining or an iron oxide–rich cement. D: Altered grain of feldspar to muscovite, with Fe oxide. E: Calcite within feldspar. F: Example of a cataclasite grain.

Figure 6.

Photomicrographs showing characteristic features from cuttings samples. Depths of the samples (meters) are noted at the lower right-hand corner of the photographs. All photographs were taken in cross-polarized light. Differences in color are due to different light levels. Scale bar applies to all photomicrographs. A: High-birefringence clay minerals fill space between quartz crystals. B: Altered plagioclase in center of photo shows alteration typical in the upper arkose. The grain is still easily recognizable as plagioclase with small areas of sericite formation. C: Grains displaying iron oxide–rich matrix/cement. The grain in the upper right of the photo shows deposition-related iron oxides (matrix), while the iron oxides in the grain in the lower left look postdepositional, either staining or an iron oxide–rich cement. D: Altered grain of feldspar to muscovite, with Fe oxide. E: Calcite within feldspar. F: Example of a cataclasite grain.

Figure 7.

Minor compositional components that make up less than 10% of any given sample. If plotted against the other components of the system, these would look insignificant in scale; thus they are plotted separately as modal percentages. Graphs have been split into the three lithologic units: upper arkose, clay-rich zone, and lower arkose.

Figure 7.

Minor compositional components that make up less than 10% of any given sample. If plotted against the other components of the system, these would look insignificant in scale; thus they are plotted separately as modal percentages. Graphs have been split into the three lithologic units: upper arkose, clay-rich zone, and lower arkose.

Figure 8.

Abundance of altered and unaltered feldspars, cataclastically deformed feldspar, and cement abundances in the drilled section, expressed in modal percentages.

Figure 8.

Abundance of altered and unaltered feldspars, cataclastically deformed feldspar, and cement abundances in the drilled section, expressed in modal percentages.

Throughout the upper arkose, the wireline logs and cuttings data suggest that there is little variation in composition. One exception to the homogeneity is in the region of 2022–2045 mmd where the Vp drops to a range of 3.20–4.6 km/s and porosity increases to 18%. This zone provides an excellent example of how cuttings microscopy can be combined with wireline logs to more fully characterize the rocks. The two cuttings thin sections evaluated here (2026.92 and 2041.4 mmd) show a marked increase in iron oxide–rich grains (17%–24% of the sample) from the samples above and below. Only 4.00%–8.00% of the grains in the samples from this depth range contain cataclasite, less than most of the other upper arkose samples. We interpret this zone to be a high-porosity, high-permeability zone that acted as a conduit for oxidizing fluids and that exhibits little evidence for faulting.

We divide the upper arkosic section into five structural blocks. Each block has distinct sedimentary and structural characteristics as revealed in their wireline log responses (Table 2). We calculate the sedimentary thickness of each package using simple trigonometry:
formula
(1)
where x is the true thickness, y is the apparent thickness, and φ is the angle between the trend and plunge of the borehole and the pole to bedding of the block. While there were minor variations in the bearing of the borehole, these were too small to significantly alter the true thickness result. We used an average bearing of 35° toward 045° for the borehole orientation in our calculations. An average pole to bedding was determined using the statistical methods of a stereonet program. A complete catalogue of stereograms and calculations for each block can be found in Table 3.
TABLE 2.

WIRELINE DATA VALUES FOR INDIVIDUAL STRUCTURAL BLOCKS

TABLE 3.

TRUE THICKNESSES OF STRUCTURAL BLOCKS

Block 1 is composed of fine- to coarse-grained, poorly bedded sandstone. Conductive and resistive spots likely represent different clast lithologies in the coarse sandstone, leading to a speckled or mottled appearance. The rocks are highly fractured in some portions. We define a fracture density of >15 fractures per 10 m drill distance to be highly fractured, 5–15 fractures per 10 m to be moderately fractured, and 0–5 fractures per 10 m to be slightly fractured. Three different bedding orientations with dip magnitudes that range from shallow to steep occur over a 10 m distance. If these trends were observed over a larger area, it could be interpreted to be a zone of postlithification tectonic folding. However, we interpret it as a soft-sediment slump fold.

In block 2, the rocks are well-bedded conglomeratic sandstone that dips 0° and 30° to the northwest, striking south to southwest. Assuming the beds are upright and nearly flat-lying, we calculate a true thickness for block 2 of 44.9 m (Table 3). With little bedding observed in block 1, the boundary between the two blocks is based on a change in rock character from massive pebble conglomerate in block 1 to distinctly bedded sandstone at 2145 mmd.

The interface between blocks 2 and 3 is marked by an abrupt dip change from shallow northwest-dipping beds to variably southwest-dipping beds (Fig. 9). Block 3 has a stratigraphic thickness of 20.1 m and consists of fine-grained sandstone and siltstone beds that dip 30°–90° southwest. Block 4 extends from 2250 to 2290 mmd, with a thickness of 18.8 m, and lithologically is nearly identical to block 2. Like block 2, the bedding in block 4 is nearly horizontal to shallowly dipping to the northwest. Bedding in block 3 is coherent and exhibits clear bedding trends throughout. The transitions from block 2 to block 3 and from block 3 to block 4 are abrupt. No diagnostic features of discrete fault planes were observed in the image logs. Block 5 has a true thickness of 72.6 m and consists of coarse-grained, speckled, sandstone-bearing zones of laminated bedding dipping 60° to 90° to the northeast. Without considering the thickness of block 1, the total measured thickness of the upper arkose is 156.4 m. If block 1 is similarly oriented to block 2, the unit would be a fining-downward sequence, and the thickness of block 1 would be 134.78 m. Adding the total maximum stratigraphic thickness, the upper arkose is 291 m thick.

Figure 9.

Southwest to northeast cross section of the entire arkosic section in the SAFOD borehole, with individual structural blocks labeled. Blocks 1–5 are the upper arkose, blocks 6–8 are the clay-rich zone, and blocks 9–11 are the lower arkose. Depths are in measured dimension, along the deviation of the borehole, in meters. Faults are noted but are not necessarily oriented correctly. Bold faults are projected to predicted surface fault equivalents. Stereonets for each block are equal-area nets with bedding plotted as great circles, black point is the trend/plunge of the borehole. Interpretation of the lowermost section is modified on the basis of core retrieved in 2007. SAF—San Andreas fault.

Figure 9.

Southwest to northeast cross section of the entire arkosic section in the SAFOD borehole, with individual structural blocks labeled. Blocks 1–5 are the upper arkose, blocks 6–8 are the clay-rich zone, and blocks 9–11 are the lower arkose. Depths are in measured dimension, along the deviation of the borehole, in meters. Faults are noted but are not necessarily oriented correctly. Bold faults are projected to predicted surface fault equivalents. Stereonets for each block are equal-area nets with bedding plotted as great circles, black point is the trend/plunge of the borehole. Interpretation of the lowermost section is modified on the basis of core retrieved in 2007. SAF—San Andreas fault.

The Clay-Rich Zone

A clay-rich zone containing as much as 60% clay minerals and clay-sized particles and registering high gamma ray signature, high porosity log values, and low seismic velocities is observed from 2530 to 2680 mmd (with the exception of a small block from 2565 to 2595 mmd). Clays cannot be fully characterized optically due to the small grain size. The clay-rich zone was initially identified at the SAFOD site by the dramatic increase in gamma ray and porosity log signatures and a decrease in Vp and Vs (Fig. 2) (Boness and Zoback, 2006). The lowest velocity in the clay-rich zone is 2.69 km/s with an average velocity of 4.00–4.27 km/s (Table 1), considerably lower than in the upper and lower arkose units. The average neutron porosity is 0.18–0.27, higher than the average porosities of 0.10 and 0.08 in the upper and lower arkose units, respectively. The clay minerals vary in composition from illite to smectite to mixed-layer illite-smectite (Fig. 5B), and from chlorite to muscovite (see Solum et al., 2006, for details). Calcite grains and grains stained by Fe oxides are more abundant in the clay-rich zone (8.26%) than in the lower arkose but less than in the upper arkose (Fig. 8). The chromium-bearing muscovite mineral fuchsite is present in the clay-rich zone and the lower arkose (Figs. 5 and 7).

Between 2565 and 2595 mmd the borehole encountered 27.33% clay, 37.66% quartz, and 10% feldspar in the clay-rich zone, similar to the composition of the lower arkose (Fig. 8). The increase in quartz and feldspar grains in this zone indicates an increase in grain size from the rest of the clay-rich zone, and wireline logs also indicate coarser-grained sediment in this interval (Table 2). Sonic velocities increase to 4.48 km/s, porosity decreases to 0.12, and there is a corresponding increase in density from 2.48 to 2.88 g/cm3. This thin interval represents a different lithology than the rest of the clay-rich interval and is not considered in the physical property calculations of the rest of the clay-rich lithological unit.

The borehole in the clay-rich zone was often elliptical, making the caliper in this unit rarely in gauge. The short axis of the caliper stayed approximately uniform at 25–30 cm diameter, but the long axis often exceeded the 50.8 cm maximum reach of the caliper tool (Fig. 2). The loss of borehole integrity results in compromised image logs throughout much of the entire clay-rich zone. Thus few fractures and no bedding orientations were recorded in the clay-rich portions of this zone. Bedding observed in block 7 strikes southeast and dips shallowly to the southwest from 2565 mmd to 2595 mmd, with a true thickness of 28.7 m (Table 3). Gamma ray and porosity values decreased while density and velocity increased compared to the rest of the clay-rich zone (Fig. 2). Point counts show a marked decrease in the amount of clay-rich grains. These data indicate that block 7, which is significantly different from the rest of the clay-rich zone, is coarser grained than blocks 6 and 8.

Lower Arkose

The lower arkose, from 2680 to 3157 mmd, is characterized by a variable wireline log signature that exhibits a wider range of values for gamma ray, velocity, and density than the upper arkose (Table 1; Fig. 2). Porosity exhibits a similar range in values in the lower arkose as in the upper arkose. However, the lower arkose exhibits a larger number of excursions, likely caused by the presence of differing lithologies (Fig. 2; see also Boness and Zoback, 2006). The porosity averages 0.08, which is generally less than the porosity calculated for the upper arkose. The average density of the lower arkose is 2.53 g/cm3. The average gamma ray value is 112.48 API, which is less than the upper arkose, and the gamma ray curve is more irregular than the homogeneous curve of the upper arkose.

The lower arkose is generally finer grained than the upper arkose. Clay-sized particles and minerals compose on average 15.5% of the samples. Clay minerals include illite and smectite and almost no kaolinite. There are similar amounts of chlorite present in the lower arkose as in the upper arkose, although altered biotite is much less common in the lower arkose.

Distinct compositional differences exist between the upper and lower arkoses. Quartz is less common in the lower arkose than the upper arkosic section, ranging from 33% to 60% of the total grains, averaging 48.12% (Fig. 4). The majority of quartz in the lower arkose is plastically deformed (Fig. 5E), and in places full sub-grains have developed in the crystal structure.

There are fewer feldspar grains overall (both altered and unaltered) in the lower arkose than the upper arkose, although the ratio of plagioclase to K-feldspar remains the same (Fig. 8). The upper arkose contains more altered feldspar grains, but the degree of alteration in the lower arkose is generally greater. Feldspar grains have been altered to muscovite in the lower arkose. There are some instances of nearly complete replacement by muscovite (Fig. 6D).

Iron oxide–rich grains have concentrations that range from 2% to 15%, averaging 5% of the sample (Figs. 4 and 6). The characteristics of the iron oxides are similar to those in the upper arkose and clay-rich zone, in some cases appearing dark red in plane- and cross-polarized light, and in other sites opaque under both forms of light but red in reflected light. Iron oxide is present as fracture fill, grain coatings, and fine-grained matrix. Porphyritic volcanic clasts that contain very fine-grained, often devitrified groundmass (Fig. 5F) are present in every thin section examined in the lower arkose, constituting an average of 1.8% of the sample. Fuchsite concentrations range from 0.1% to 2.0% of the sample in the lower arkose (Fig. 7). Fuchsite occurs as small grains within a larger fine-grained, clay-rich grain and rarely as discrete grains, in contrast to the upper arkose.

The lower arkose contains a larger amount of laumontite than the upper arkose. Solum et al. (2006) also detected relatively large amounts of laumontite through X-ray diffraction (XRD) analysis in the lower arkose with little to no laumontite in the upper arkose. This evidence led them to define two distinct arkosic units similar to our designations. The geophysical signature of the rocks (Boness and Zoback, 2006; Bradbury et al., 2007) likewise indicates that the two arkosic units are lithologically distinct.

Image log data, where quality is good, provide further constraints on lithofacies and bedding orientation (Figs. 10 and 11). We divide the lower arkose into three blocks (Fig. 9) on the basis of the image log data. Block 9 has a true thickness of 158.9 m (Table 3), dips shallowly to the south and southwest, and consists of homogeneous fine-grained sandstones cut by fractures parallel or perpendicular to bedding. Block 10 is lithologically similar to block 9, consisting of well-bedded, laminated, and locally fine-grained sandstone and siltstones. An abrupt change in dip magnitude and direction is observed at 2880 mmd. Bedding dips steeply and uniformly to the northeast to a depth of 3010 mmd. (The borehole extended to 3050 mmd, but due to the tool configuration, image logs were not collected to the bottom of the borehole.) The apparent thickness of block 10 is 170 m if we assume that the block continues to the shear zone encountered in the phase 1 spot core. The calculated thickness is 65.1 m (Table 3). Representative image logs from block 10, along with the corresponding composition (Fig. 12), show the bedding to be 5–10 cm thick. Block 11 is the last block in the sequence drilled during phase 2 drilling, and image logs are not available.

Figure 10.

Examples of different lithologies or image log “facies” for the upper arkose. Each example is displayed in a 1:1 aspect ratio but at different scales to highlight features. The color scale represents the electrical resistivity of the rock. The lighter shading represents higher resistivities. The four strips are the images resulting from the four pads that are arranged on the tool 90° apart. The sinusoidal lines in the images are examples of lithologic boundaries or fractures picked from the image. A: Conglomeratic interval in the upper arkose. Notice the abundance of resistive clasts, likely granitic in composition. B: Chaotic interval in the upper arkose with good image clarity but no discernable bedding or fractures. This interval is interpreted to be a slump fold or similar soft-sediment deformation due to the lack of brittle-deformation features that would be present in postlithification deformation. C: An example of poor image quality in the clay-rich zone. Image may represent a brecciated or cataclastic zone, but the clast/matrix features are more likely a result of blurring from poor data. The majority of the clay-rich zone images have these characteristics. D: An uninterpreted image of granite in the borehole above the arkoses. Image quality is good; granite is highly fractured.

Figure 10.

Examples of different lithologies or image log “facies” for the upper arkose. Each example is displayed in a 1:1 aspect ratio but at different scales to highlight features. The color scale represents the electrical resistivity of the rock. The lighter shading represents higher resistivities. The four strips are the images resulting from the four pads that are arranged on the tool 90° apart. The sinusoidal lines in the images are examples of lithologic boundaries or fractures picked from the image. A: Conglomeratic interval in the upper arkose. Notice the abundance of resistive clasts, likely granitic in composition. B: Chaotic interval in the upper arkose with good image clarity but no discernable bedding or fractures. This interval is interpreted to be a slump fold or similar soft-sediment deformation due to the lack of brittle-deformation features that would be present in postlithification deformation. C: An example of poor image quality in the clay-rich zone. Image may represent a brecciated or cataclastic zone, but the clast/matrix features are more likely a result of blurring from poor data. The majority of the clay-rich zone images have these characteristics. D: An uninterpreted image of granite in the borehole above the arkoses. Image quality is good; granite is highly fractured.

Figure 11.

A composite image from a 2 m interval in the lower arkose. Here the rocks are laminated fine-grained sandstone/siltstones. Solid purple lines on image log are interpreted fractures, while dashed purple lines are interpreted bedding planes. Image quality in this interval is some of the best in the entire arkosic interval.

Figure 11.

A composite image from a 2 m interval in the lower arkose. Here the rocks are laminated fine-grained sandstone/siltstones. Solid purple lines on image log are interpreted fractures, while dashed purple lines are interpreted bedding planes. Image quality in this interval is some of the best in the entire arkosic interval.

Figure 12.

Photomicrographs of phase 1 spot core and phase 2 sidewall cores. All photos except photo B are taken in cross-polarized light; photo B is taken in plane-polarized light. Scale bar applies to all photomicrographs except E and H. Depths are based on the distance from the lower end of the core to the sample collected. A: Fractured or brittlely damaged quartz grain. B: Quartz grain broken apart with individual slivers remaining in place, the beginning of grain size reduction. C: Fine-grained, devitrified volcanic clast takes up left half of photo. Opaque minerals at center of photo and edge of volcanic clast are not identifiable in reflected light. D: Altered granitic clast with amphibole beginning to oxidize. From sidewall core 69. E: Veins filled with laumontite. F: Iron oxide cementation in sidewall core 71. G: Localized calcite cementation in lower left-hand corner of photo. Grains are extremely fractured in this portion of the core. H: Postdepositional muscovite vein formation in quartz grain. I: Siltstone portion of spot core; grains are texturally immature with little to no mineral alignment. Closely resembles fine-grained cuttings.

Figure 12.

Photomicrographs of phase 1 spot core and phase 2 sidewall cores. All photos except photo B are taken in cross-polarized light; photo B is taken in plane-polarized light. Scale bar applies to all photomicrographs except E and H. Depths are based on the distance from the lower end of the core to the sample collected. A: Fractured or brittlely damaged quartz grain. B: Quartz grain broken apart with individual slivers remaining in place, the beginning of grain size reduction. C: Fine-grained, devitrified volcanic clast takes up left half of photo. Opaque minerals at center of photo and edge of volcanic clast are not identifiable in reflected light. D: Altered granitic clast with amphibole beginning to oxidize. From sidewall core 69. E: Veins filled with laumontite. F: Iron oxide cementation in sidewall core 71. G: Localized calcite cementation in lower left-hand corner of photo. Grains are extremely fractured in this portion of the core. H: Postdepositional muscovite vein formation in quartz grain. I: Siltstone portion of spot core; grains are texturally immature with little to no mineral alignment. Closely resembles fine-grained cuttings.

Block 11 has an apparent thickness of 107 mmd. The extent of block 11 ends at the abrupt lithology change seen at 3157 mmd, where the cuttings character changes from arkosic to very fine-grained mudstones and siltstones (Evans et al., 2005; Hickman et al., 2005; Bradbury et al., 2007). No bedding or fracture data were available for this section of the hole. The thicknesses of blocks 9 and 10 give a minimum thickness of 224 m if we assume the apparent thickness of block 11 to be the maximum possible thickness of that block. The total maximum thickness of the lower arkose is 331 m.

Fractures were identified on the basis of their resistivity character and disruption of bedding (Figs. 10 and 11). We observe four zones of relatively high fracture density (>10 fractures per 10 m) in the upper arkosic rocks. These are found at 1920–1980, 2020–2085, 2190–2215, and 2275–2530 mmd. Below the clay-rich zone, a fracture zone occurs from 2515 to 2725 mmd, and much of the lower arkosic section is highly fractured (Fig. 2). Several of these zones of increased fracture density are associated with boundaries between structural blocks, and the uppermost region is associated with the fault that juxtaposes the Salinian rocks against the arkosic rocks. The increase in fracture density around the clay-rich zone may suggest that this region is a fault zone, as suggested by Boness and Zoback (2006).

Core Analysis

The spot core (Fig. 3) collected during phase 1 drilling intersected a section of the coarse-grained sandstone and conglomerate and a small clay-rich shear zone at 3067 mmd. Other studies have examined the composition and frictional properties of this shear zone (Tembe et al., 2005; Schleicher et al., 2006; Solum et al., 2006). The sedimentary rocks in the spot core are thoroughly damaged and deformed (Almeida et al., 2005) (Figs. 3 and 12), perhaps as a result of the faulting that has occurred in the cored shear zone. The core samples a sequence of quartzo-feldspathic, indurated, and coarse-grained arkose and conglomerate. Over the 12 m length of core, grain sizes grade abruptly from a coarse pebble/cobble conglomerate to coarse sandstone, and from fine sandstone to siltstone (Fig. 3). The pebbles and cobbles of the conglomerate are granitic and volcanic, consisting of significantly more granite than volcanic clasts.

Porphyritic volcanic clasts are similar to those observed in the lower arkose cuttings. They contain large plagioclase crystals in a fine-grained devitrified groundmass (Fig. 12C). In some cases, phenocrysts are oriented in a primary flow texture. Granitic clasts commonly are altered (Fig. 12D) and composed of quartz, feldspar, amphiboles, and chloritized biotite. The majority of laumontite observed in the core occurs in granitic clasts as veins or cement between grains (Fig. 12E).

Macroscopically and microscopically the core is stained by iron oxides (Fig. 12F). The arkose is highly damaged, riddled with micro-fractures, and shows little evidence of healed fractures (Fig. 12G). Grains are locally broken, and evidence for grain size reduction is common (Fig. 12B). Evidence for cementation is not commonly observed. Rare veins are filled with laumontite, calcite (Fig. 12G), secondary quartz, or muscovite (Fig. 12H). Feldspar grains are generally unaltered through the core. The lack of significant amounts of laumontite, calcite, and secondary quartz corresponds to what has been observed in the cuttings. Authigenic muscovite was observed as vein material in several locations in the core samples (Fig. 12H). These veins, no more than 0.10 mm wide, cut across several grains, which is evidence that they are postdepositional features.

The fine-grained portions of the spot core are similar to grains counted as clay-rich in the cuttings (Fig. 12I). The fine-grained rocks are composed of texturally immature quartz and feldspar in a matrix of high-birefringence clay minerals. In most areas of the core, clay minerals are a depositional feature rather than a secondary cementation. Small angular mica grains consisting of muscovite and fuchsite are similar to those seen in the finer-grained cuttings and those observed in the spot core samples. The clay minerals do not exhibit any preferred orientation in the core samples. No microfossils were observed in either cuttings or core, even in the finest-grained intervals.

Diagenetic Textures

The thin-section petrography of the core and cuttings samples indicates several types of cement and alteration products that may have been caused by diagenesis of the sedimentary sequence. Calcite, laumontite, secondary quartz, iron oxide, and clay cements were all observed to some degree in the cuttings as well as the core samples. Iron oxide cement is most commonly observed as grain coatings and as interstitial cement (Figs. 6C and 12F). The iron oxide cements were most likely produced from the oxidation of biotite and amphibole coming from a granitic source (Dapples, 1979). This reaction is common in arkosic sandstones because biotite is rarely stable enough for the reaction products to remain long after deposition (i.e., Scholle, 1979; Carozzi, 1993).

Calcite and silica cements are rare. The majority of calcite grains are either discrete calcite crystals or feldspar altering to calcite. Very few calcite veins or grains cemented with calcite were encountered. The spot core also showed very little evidence of these cementing mechanisms as a pervasive trend. Most calcite veins observed in the spot core thin sections occur near cataclasite zones. Quartz overgrowths are rare with quartz veins ranging from 0.21% to 0.44% of the sample counted in the cuttings (Fig. 8). The compaction features observed in the quartz and feldspar grains are not as well developed as stylolites, but there are instances of convex-concave contacts and quartz indentor grains, similar to those features seen in deeply buried clastic sedimentary rocks (i.e., Liu, 2002).

The ubiquitous alteration products in these rocks are illite, muscovite, and laumontite. Illite appears to replace K-feldspar in the upper arkose, whereas muscovite is the more common alteration phase in the lower arkose. Other workers have examined the morphology of SAFOD clay minerals using SEM and TEM methods and have found that mixed-layer illite-smectite forms as films and fibers along shear surfaces along with fibrous to columnar laumontite forming in the pores of rock chips collected during coring (Solum et al., 2006; Schleicher et al., 2006).

If the illite in the SAFOD arkoses formed at 70–150 °C (Huang, 1992; Worden and Morad, 2003), then these rocks may have been buried at depths slightly greater than those at which they were encountered. The measured borehole temperatures yield a modern geothermal gradient of 38 °C/km (C. Williams, 2006, written commun.), within the range of or slightly higher than the temperature gradient for this area (Sass et al., 1997; Page et al., 1998; Blythe et al., 2004). Using the measured geothermal gradient, 70 °C to 150 °C corresponds to depths of 1.3 km to 3.4 km. This overlaps with the current depth of the arkose units but may be as much as 0.8 km deeper than the deepest present extent of the arkose.

Zircon Fission-Track Analysis

There are no direct measures of the age of these rocks available from either micropaleontology or radiometric analyses. Zircon fission-track analysis was used on detrital zircons to constrain the thermal history of the arkosic rocks (Wagner and Van der Haute, 1992; Bernet and Garver, 2005), and we combine the fission-track analyses with other tectonic studies of the area to determine the cooling history of the source rocks to the SAFOD arkoses and to constrain the age of deposition (Figs. 9 and 13). Uranium concentrations in the samples range from 50 to 450 ppm U, with a mean between 200 and 250 ppm (Table 4), which is typical for detrital suites (see Garver and Kamp, 2002). Four of the samples come from depths of 3122–3160 mmd, at the base of the arkosic section in block 11 and possibly slightly below or within a fault (Bradbury et al., 2007) (Figs. 9 and 13). Two samples come from depths of 3338 and 3375 mmd, within the region interpreted to be the fault zone (Hickman et al., 2008).

Figure 13.

Location of the six zircon fission-track samples in relation to the inferred subsurface geology and location of earthquakes below the hole. Four of the samples come from the lower arkose, and two come from within the low-velocity zone. This is an enlarged view of the lower part of the borehole, and some of this interpretation is a result of coring in 2007. This view encompasses the lower portion of block 11. The Vp data are plotted along the borehole, and the cross indicates the ±50 m uncertainty associated with the earthquake locations (C. Thurber, 2009, personal commun.). Two of the fault zones interpreted to be at the edge of or within the low-velocity zone are shown.

Figure 13.

Location of the six zircon fission-track samples in relation to the inferred subsurface geology and location of earthquakes below the hole. Four of the samples come from the lower arkose, and two come from within the low-velocity zone. This is an enlarged view of the lower part of the borehole, and some of this interpretation is a result of coring in 2007. This view encompasses the lower portion of block 11. The Vp data are plotted along the borehole, and the cross indicates the ±50 m uncertainty associated with the earthquake locations (C. Thurber, 2009, personal commun.). Two of the fault zones interpreted to be at the edge of or within the low-velocity zone are shown.

TABLE 4.

SUMMARY OF ZIRCON FISSION-TRACK DATA FOR SAFOD MAIN HOLE

Most of the zircon grains have statistically similar Late Cretaceous to earliest Paleocene annealing ages of 70–62 Ma, passing the chi-square test evaluating the likelihood that the zircon ages cluster (Table 4). Several samples show some age dispersion, but overall the data are consistent with a common thermal history for all sampled grains. This homogeneity implies that the dated zircons from the arkoses are likely derived from a single source terrain or from source regions with a homogeneous thermal history. The grains are mainly clear and euhedral: The suites of zircons from these units lack rounded, abraded grains or grains with significant radiation damage. We observe no significant variation in zircon fission-track age as a function of depth, and this, along with the narrow range of colors and grain shapes, suggests that the zircons were derived from the same thermally homogeneous source, and it was likely that the source of the zircons was first-cycle detritus, otherwise rounding and pitting of grains would have been common.

The source rock of the arkose units cooled through the effective closure temperature of zircon (~240 °C) during the Late Cretaceous, between 64 and 70 Ma. The zircons do not show evidence of being reset after deposition, or buried to temperatures exceeding 200–240 °C (see Garver et al., 2005). This inference is supported by low thermal maturation of organic detritus from this interval that has %Ro values of 0.74–1.01. In the SAFOD hole, we predict that temperatures would equal or exceed the zircon annealing temperature of 240 °C at ~5.7 km depth. Assuming a constant geothermal gradient over time, the lack of evidence for thermal resetting of the grains implies that the arkosic rocks have not been buried to a depth greater than the present depth of 5.7 km, or 3.1 km more than zircon fission-track samples 37–61 and 37–62 (Table 4) since deposition.

The four fission-track ages from the lower arkose samples and the two ages from the lower fine-grained rocks show that the two groups have similar thermal histories. We are confident that the higher set of samples was derived from a relatively restricted part of the arkosic section, as the hole was cased to a depth of 3065 mmd at the end of phase 1 drilling in 2004. The lower two samples have consistent cooling ages and %Ro values, pointing to a common thermal history for the samples. This indicates that the fault zone may consist of slivers of rocks from a range of sources, including the arkosic and clay-rich rocks we ascribe to the Salinian terrain.

The Maastrichtian–Paleocene cooling ages of the zircons recovered from the four samples derived from the arkosic section serve as a maximum depositional age for the SAFOD arkoses. Based on the petrography, composition, and zircon fission-track analyses, it seems unlikely that the SAFOD arkosic rocks are part of the Great Valley Group. Vermeesch et al. (2006) report zircon fission-track ages of 93.3–106.5 Ma for Upper Cretaceous Great Valley samples from ~15 km east of the SAFOD site, indicating a very different thermal history for rocks northeast of the San Andreas fault. Degraaff-Surpless et al. (2002) report similar ages, and cite the work of Linn (1991) to indicate that Great Valley Group rocks northeast of SAFOD have depositional ages 10–15 million years older than the fission-track cooling ages that we report. A Late Cretaceous exhumation rate of >2–3 mm/a has been calculated for the Salinian block (Kidder et al., 2003), the likely source terrain for the detrital zircon. Using this rate of exhumation, we estimate that the granitic body that supplied zircons in the SAFOD arkoses could have reached the surface 1.9–7 million years after passing through the zircon annealing temperature. Following this reasoning, the maximum depositional age for these rocks is latest Cretaceous (younger than 70 Ma) to Paleocene.

INTERPRETATION

We synthesize our data with other data from the SAFOD rocks, as well as other studies in the area, to interpret the depositional setting of the arkosic rocks, and to constrain the structure of the fault zone and offset history of the San Andreas fault at this site. We begin with a discussion of the origin of the sedimentary rocks. Next we discuss how these rocks fit into the structure of the fault zone, and in particular the tectonic setting of the arkosic sections relative to the San Andreas fault.

Structural Interpretation

We have identified three rock units in the section—the upper arkose, clay-rich zone, and lower arkose—which differ in framework composition, alteration, and deformational features. Interpretation of stratigraphic and age relationships between the three units is complicated by the likely presence of faults or folds in the arkosic section. With a total of ten possible faults identified, the entire arkosic section has been offset and deformed to such a degree that it may be impossible to precisely determine original depositional or stratigraphic relationships.

The data presented here, along with Boness and Zoback (2006), Solum et al. (2006), and Bradbury et al. (2007), indicate that the arkosic section is bounded and cut by faults. The three main faults lie at 1920, 2560, and 3067 mmd. The highest fault juxtaposes the Salinian granitic rocks and the arkosic rocks; the middle fault juxtaposes the two arkosic units. The lower fault juxtaposes the arkosic section against a fine-grained sedimentary unit. In addition, we interpret the presence of as many as nine other faults or narrow fold hinge zones within the arkosic section. Several studies provide a variety of data that support a model that the clay-rich region is a fault zone (Boness and Zoback, 2006; Solum et al., 2006; Bradbury et al., 2007). The clay-rich zone is a boundary between lithologically distinct units that are from different sources. The difference in composition between the upper and lower arkoses could be the result of juxtaposition of two units from different basins and different ages.

The uppermost fault contact between the Salinian granitic rocks and the arkosic section is interpreted to be the downdip continuation of the Buzzard Canyon fault and its associated damage zone. The lower fault encountered at the end of phase 2 drilling juxtaposes the arkosic section against sheared mudstone. These data indicate that the fault zone structure at depth is complex, and that the seismicity associated with the current San Andreas fault lies below and eastward of the lowest fault observed here (Hickman et al., 2008). Geophysical investigations of the site suggest that the active fault zone is narrow, and the region between what we infer to be the Buzzard Canyon fault and the San Andreas fault is a zone of high-Vp rocks (Thurber et al., 2004) overlain by low-velocity sedimentary rocks. Other workers suggest that the fault zone at this depth is relatively wide and marked by low-velocity or seismically imageable rocks associated with damage zones (Hole et al., 2001; Bleibenhaus et al., 2007; Catchings et al., 2007; Li and Malin, 2008), and that some of these faults extend to the near-surface region (Catchings and Rymer, 2002). Based on the work presented here and in related papers (Bradbury et al., 2007; T.N. Jeppson, 2009, personal commun.), we infer that the faults within the section have relatively modest offsets and do not affect the overall velocity character of the section. Wireline- and image-log characteristics along with the lithologic analyses indicate that rocks on either side of these inferred faults are similar, and the boundaries therefore displace or fold a relatively uniform package of rocks within each main rock body.

Depositional Environments

We use the petrographic and compositional data to infer depositional environments. The presence of subangular to subrounded feldspar grains and the textural immaturity of the arkosic rocks indicate that they were deposited in a proximal or high-energy setting. We cannot determine whether the arkosic rocks were deposited in a marine or terrestrial setting. If they are marine sediments, they likely represent the most proximal portion of a marine fan sequence. Clasts in the arkoses are almost entirely granitic, including no more than 2% volcanic clasts in the lower arkose, which implies that the arkosic units were shed from a granitic highland that contained a minor source for volcanic clasts.

Based on the range of grain sizes, the somewhat regular bed spacing, the angularity of clasts observed in image logs and core, and the textural immaturity associated with unaltered feldspars, we interpret the arkoses to have been deposited in a subaerial or subaqueous fan setting weathered from a granitic source terrain. Overall, there are more sandstone and siltstone intervals than conglomerate in the arkosic section, indicating periods of low energy alternating with intervals of high-energy deposition. The alternation of the fine-grained sequences and the coarse-grained arkosic sandstones and conglomerates may represent the alternating high- and low-energy portions of submarine flows (e.g., Falk and Dorsey, 1998; Clifton, 2007) in which coarse-grained sequences are deposited in relatively deep water. Chaotic intervals in the image logs may represent slump or slope failure deposits, both of which could be present in a number of environments, including a fluvial or marine fan setting.

Stratigraphic Correlations

How might these rocks correlate with rocks observed in the field or in drill holes in western California? The Late Cretaceous to early Eocene age estimate of these rocks combined with their petrography and physical properties suggests several likely correlative units. Late Cretaceous to Eocene siltstone to conglomerate units occur in small but widely distributed exposures throughout western California (Grove, 1993; Seiders and Cox, 1992; Burnham, 2009) (Fig. 14). Given the petrography and general age constraints, we suggest the rocks are latest Cretaceous to Paleocene–Eocene deposits on or associated with the Salinian block. This correlation may help constrain the evolution of the San Andreas fault in central California, because other correlation studies have successfully established amounts of offset on different strands of the San Andreas fault (Nilsen and Clarke, 1975; Nilsen, 1984; Graham, 1978; Graham et al., 1989).

Figure 14.

Three slip distributions for the San Andreas fault system in the vicinity of the SAFOD site. A sedimentary rock unit correlative to the SAFOD arkoses would be found in a different location depending on which hypothesis is correct. Gray unit represents the SAFOD arkoses and the equivalent sedimentary rocks. A: Slip on the San Andreas fault, with very small amounts of slip on Buzzard Canyon fault, results in roughly 300 km of slip, displacing the arkosic section such that arkosic rocks would lie northwest of SAFOD. B: Partitioning of slip on the Buzzard Canyon fault and San Andreas fault satisfies the total slip (see Fig. 14D) and results in arkosic rocks spread >30 km northwest of SAFOD. C: If most of the slip were on the Buzzard Canyon fault, arkosic rocks would be observed in the subsurface at the western edge of the Great Valley. Note that the star that shows the location of SAFOD shifts to show the relative location of the sediments, and is not shifted in an absolute sense. D: For reference, offset units along the San Andreas fault zone. Location of SAFOD drilling site marked by star. Eagle Rest Peak–Logan offset based on Ross (1970); Pinnacles–Neenach from Matthews (1976). Location of La Honda and San Joaquin Basins from Clarke and Nilsen (1973) and Graham et al. (1989). Map modified after Ingersoll (1988) and Irwin (1990). Dark ellipse near bottom of map represents the area that is ~315 km southeast of the SAFOD site on the opposite side of the San Andreas fault.

Figure 14.

Three slip distributions for the San Andreas fault system in the vicinity of the SAFOD site. A sedimentary rock unit correlative to the SAFOD arkoses would be found in a different location depending on which hypothesis is correct. Gray unit represents the SAFOD arkoses and the equivalent sedimentary rocks. A: Slip on the San Andreas fault, with very small amounts of slip on Buzzard Canyon fault, results in roughly 300 km of slip, displacing the arkosic section such that arkosic rocks would lie northwest of SAFOD. B: Partitioning of slip on the Buzzard Canyon fault and San Andreas fault satisfies the total slip (see Fig. 14D) and results in arkosic rocks spread >30 km northwest of SAFOD. C: If most of the slip were on the Buzzard Canyon fault, arkosic rocks would be observed in the subsurface at the western edge of the Great Valley. Note that the star that shows the location of SAFOD shifts to show the relative location of the sediments, and is not shifted in an absolute sense. D: For reference, offset units along the San Andreas fault zone. Location of SAFOD drilling site marked by star. Eagle Rest Peak–Logan offset based on Ross (1970); Pinnacles–Neenach from Matthews (1976). Location of La Honda and San Joaquin Basins from Clarke and Nilsen (1973) and Graham et al. (1989). Map modified after Ingersoll (1988) and Irwin (1990). Dark ellipse near bottom of map represents the area that is ~315 km southeast of the SAFOD site on the opposite side of the San Andreas fault.

The age and history of the Salinian block matches the interpreted thermal history of the source terrain for the SAFOD arkoses. The fission-track cooling ages of detrital zircons in the arkoses are younger than U/Pb geochronology that dates Salinian magmatism at between 130 and 70 Ma (Kistler and Champion, 2001; Barth et al., 2003; Kidder et al., 2003; Barbeau et al., 2005). The 40Ar/39Ar cooling ages of biotite of 76–75 Ma in the Salinas Valley, northwest of SAFOD, are consistent with the U/Pb ages of the Salinian block (Barth et al., 2003). From ca. 80 to 65 Ma the Salinian block was rapidly exhumed from 25 km depth (Barbeau et al., 2005). The zircon fission-track ages represent when the zircons passed through the effective closing temperature. Our data and the data of others indicate that the Paleogene exposed a Salinian granitic source with 70–64 Ma cooling ages at the surface. The homogeneous zircon fission-track cooling ages of the SAFOD arkoses likely reflect their derivation from a granitic source with a homogeneous history, and rule out a more complex and lithologically heterogeneous source terrain.

Clay alteration reactions in the arkosic rocks limit the depth of burial. Illitization of K-feldspar occurs between 70 °C and 150 °C, which corresponds to depths of 1.30–3.4 km. Authigenic muscovite forms at >100 °C, corresponding to a depth of 2.1 km given the modern geothermal gradient. Authigenic muscovite is only observed in the lower arkose at 2.3 km vertical depth; the upper arkose has limited feldspar alteration with almost no evidence of authigenic muscovite. The upper arkose extends from 1.9 to 2.2 km vertical depth; the lack of authigenic muscovite may indicate a different fluid chemistry at that level than in the lower arkose. In that case, authigenic clay mineral development is more dependent on fluid chemistry than temperature or pressure (Montoya and Hemley, 1975). In light of these observations, we suggest these rocks were not buried much more deeply than their present depth.

There are few direct constraints on the youngest possible age for the SAFOD arkosic rocks. They are interpreted to be unconformably overlain by the Etchegoin and Santa Margarita Formations (Thayer and Arrowsmith, 2005), which are 10–4 Ma old (see Sims, 1993). Stratigraphic studies of strata on the Salinian block suggest that most deposition ceased in the late Eocene to early Oligocene, although Oligocene deposition is recorded in some places of the western Coast Ranges (see discussion in Barbeau et al., 2005). Rocks younger than the Neenach/Pinnacles volcanic rocks (early Miocene, ca. 23 Ma) exhibit a markedly different character from the SAFOD section due to the marine sedimentation and initiation of the transform plate margin. Thus, we propose that the SAFOD arkoses were deposited sometime in the very latest Cretaceous to the Paleocene or Eocene, during exhumation of the nearby Salinian granite.

If this age assignment is correct for the arkosic rocks, then it is worth briefly considering possible correlatives. In evaluating possible correlative strata, we are interested in finding similar lithologies allied with the Salinian basement. This regional assessment will allow us to form a testable hypothesis for future stratigraphic work. Paleocene and Eocene sedimentary rocks are common in central California (i.e., Dibblee, 1971, 1973, 1980; Nilsen and Clarke, 1975; Graham, 1978; Clifton, 1984; Bent, 1988; Graham et al., 1989; Seiders and Cox, 1992; Sims, 1990, 1993; Cronin and Kidd, 1998; Dibblee et al., 1999; Dickinson et al., 2005; Burnham, 2009) and are mostly marine fan deposits deposited in basins that formed on or next to the Salinian/Sierran block (Nilsen and Clarke, 1975; Nilsen, 1984; Ingersoll, 1988; Graham et al., 1989; Wood and Saleeby, 1997) (Fig. 14). Paleocene conglomerates that lie on the Salinian block (Bachman and Abbott, 1988; Clifton, 1984; Seiders and Cox, 1992; Cronin and Kidd, 1998) may also correlate with the SAFOD arkoses. The San Francisquito Formation (Upper Cretaceous to Paleocene) and Carmelo Formation (Paleocene to lower Eocene) are composed of granite-rich conglomerates interbedded with arkosic sandstones and mudstones. The Paleocene–Eocene San Francisquito Formation crops out in the Ridge Basin, the Liebre Mountain block, and the Pinyon Ridge block (Seiders and Cox, 1992) (Fig. 14), and the Late Cretaceous Carmelo Formation lies on granodiorite at Point Lobos (Clifton, 1981, 2007). Rocks older than Miocene are not mapped at the surface in the immediate vicinity of the SAFOD site (Dibblee, 1971; Sims, 1990; Thayer and Arrowsmith, 2005). The closest Paleocene–Eocene sedimentary rocks similar to the SAFOD arkoses on the southwestern side of the San Andreas fault are thick sequences of unnamed marine sedimentary rocks in the La Panza and Sierra Madre Ranges ~50 km southwest of the SAFOD site (Fig. 14) (Jennings and Strand, 1976; Bachman and Abbott, 1988). The rocks in this area to the southwest consist of interbedded arkosic sandstones, siltstones, clay shale, and granitic conglomerates ranging in age from Upper Cretaceous to middle Eocene (Chipping, 1972; Vedder and Brown, 1968; Dibblee, 1973; Bachman and Abbott, 1988; Nilsen, 1986). The unnamed unit has been estimated to be between 600 and 910 m thick, but this encompasses the entire age range. The unit is not subdivided into different age units due to the paucity of fossils encountered (Dibblee, 1973). A complex series of folds and faults could result in similar rocks in the subsurface adjacent to the San Andreas fault.

The Butano Sandstone (Eocene) of the Santa Cruz Mountains and the Point of Rocks Formation (middle to upper Eocene) of the southwestern San Joaquin Basin are similar in provenance, sedimentary characteristics, and depositional environment (Clarke, 1973; Clarke and Nilsen, 1973). The Butano Sandstone is part of the La Honda Basin, ~200 km to the northwest of the SAFOD site (Fig. 14), and the Point of Rocks Formation is part of the southern San Joaquin Basin, ~100 km to the southeast of the SAFOD site (Clarke and Nilsen, 1973; Graham et al., 1989). The Butano–Point of Rocks Formation is composed of moderately well-indurated arkosic conglomerate, sandstone, and siltstone sequences hundreds to thousands of meters thick (Critelli and Nilsen, 1996), and in some ways is similar to the SAFOD arkoses. The Butano Sandstone in the Santa Cruz Mountains has been divided into three members: a 9–80 m thick bedded upper sandstone, a 75–230 m thick silt-stone member, and a >460 m thick sandstone interbedded with pebble conglomerate (Clark, 1981). The three members of the Butano Sandstone are similar in lithology and thickness to the SAFOD arkoses, although at the surface the Butano Sandstone is significantly less well indurated than the SAFOD arkosic rocks, and we find few of the fine-grained sandstones common in the Butano–Point of Rocks sections in the SAFOD arkoses.

A third candidate for rocks equivalent to the SAFOD arkosic section are the sandstones and conglomerates of the Eocene Tejon Formation, which are exposed at the southern end of the Great Valley (Fig. 14) (Critelli and Nilsen, 2000). These rocks consist of cobble- to pebble-rich conglomerates of granitic and metamorphic rock source, fine-grained sandstones, and siltstones and mudstones. The Tejon Formation rocks are petrologically different from the Butano–Point of Rocks sequence, having been deposited in a proximal fan setting (Critelli and Nilsen, 2000). However, the Tejon Formation contains a small amount of gabbro clasts derived from mafic basement upon which the Tejon Formation lies (Critelli and Nilsen, 2000), and we found no gabbro grains in the SAFOD arkoses. It is possible that mafic grains did not survive the drilling and washing for this study and might be present at depth.

There are a number of possible correlative units that are worth considering in our stratigraphic comparison. One problem we encounter is that although many of these units are lithologically and compositionally similar, there are few comparable diagnostic criteria that we can use to separate or distinguish between them. For example, our zircon fission-track results might allow for a distinction between the specific exhumation evolution of the source rocks through time and in different areas, but comparable data from other units are not available. So, one very fruitful avenue of future research would be provenance studies of these units using varietal studies, such as the zircon fission-track work on the SAFOD arkoses presented here.

Tectonic Implications

The structural and tectonic issues posed by the presence of these rocks include deciphering the fault zone structure at depth, how these rocks correlate with other latest Cretaceous to early Tertiary arkosic rocks in the region, and the implications for how slip is distributed along the San Andreas fault.

The tectonic history of the San Andreas fault is commonly thought to be reasonably simple in central California, by having the majority of offset occurring along a single strand of the modern trace of the fault (Irwin, 1990; Wright, 2000). However, other studies have documented the presence of several parallel strands along this portion of the fault in this area, and slip along the strands has switched over time (Dickinson, 1966; Rymer, 1981; Ross, 1984; Sims, 1993), as documented along the southern San Andreas fault (Powell, 1993). One such strand appears to be the Buzzard Canyon fault (Arrowsmith, 2007; M. Rymer, 2005, personal commun.; Thayer and Arrowsmith, 2005). The presence of the SAFOD arkosic section and a set of unnamed sedimentary rocks in the fault zone from 10 to 80 km northwest of the SAFOD site (Gribi, 1963; Walrond and Gribi, 1963; Forrest and Payne, 1964a, 1964b; Walrond et al., 1967; Rogers and Church, 1967; Dibblee, 1971) adds a measure of complexity to the offset history of the San Andreas system near the SAFOD site. The 315 km of total Neogene offset on the San Andreas fault zone is well established (Ross, 1970; Clarke, 1973; Clarke and Nilsen, 1973; Matthews, 1976; Graham et al., 1989; Sims, 1993; Dickinson et al., 2005; Burnham, 2009), with key offset markers including the San Joaquin–La Honda Basin sequences, the Pinnacles–Neenach andesite-dacite-rhyolite complexes, and the gabbros and gabbroic conglomerates of the Eagle Rest Peak–Logan–Gold Hill–Gualala system. The presence of multiple fault strands and fault-bounded zones of rock from SAFOD northwest to the central Gabilan Range on the southwest side of the fault (~95 km strike distance along the fault) and the presence of exotic lithologies, including Cretaceous (?) sedimentary rocks, metamorphosed limestones, granitic rocks, and small blocks of Franciscan rocks indicate that the fault consists of multiple strands along this section (Jennings and Strand, 1958; Dibblee, 1971; Rymer, 1981; Ross, 1984).

The presence of both the Buzzard Canyon fault and the seismically active San Andreas fault bounding the SAFOD arkosic section suggests three models for slip on the San Andreas fault system in the vicinity of SAFOD (Fig. 14): (1) The majority of dextral offset occurred on the current trace of the San Andreas fault, and only a small component of offset on the Buzzard Canyon fault (Fig. 14A); (2) the 315 km offset could have been distributed more evenly across both the Buzzard Canyon and related faults and the modern San Andreas fault (Fig. 14B); or (3) the Buzzard Canyon fault could have been the major strand of the system until recently, and the majority of offset has occurred along the Buzzard Canyon fault rather than the modern San Andreas fault (Fig. 14C). Because of the lack of a piercing point in the SAFOD arkosic rocks, we can only consider ranges of offset. It is important to explore the implications of these three models because each model results in a specific prediction that can be tested by future stratigraphic studies.

In model 1, rocks equivalent to the SAFOD arkosic section would lie ~300 km to the south, on the northeast side of the fault. This model predicts a small amount of offset across the Buzzard Canyon fault (Fig. 14A). If most of the slip occurred along the San Andreas fault, the correlative sequence should be found ~300 km southeast of the SAFOD site, southeast of the Neenach volcanic complex (Fig. 14A), and the SAFOD arkoses would have correlatives nearby and west of the SAFOD site.

Model 1 is not supported by the available data. Correlative rocks to the arkosic section are the Paleocene and Eocene marine fan deposits of the La Honda Basin (Nilsen and Clarke, 1975; Nilsen, 1984; Graham et al., 1989) (Fig. 14D). Restoring 315 km of right-lateral slip in the San Andreas system does not place these rocks at the SAFOD locality. No Paleocene and Eocene sedimentary rocks are found northeast of the San Andreas fault, and one possible correlative, the San Francisquito Formation, lies within a fault-bounded block southwest of the San Andreas fault zone and provides no offset constraints. Overlap relationships interpreted near SAFOD rule out a Miocene age for the arkosic section in the SAFOD drill hole (Thayer and Arrowsmith, 2005). Eocene rocks in the Orocopia Mountains (Crowell and Susuki, 1959) are the only other known arkosic rocks on the northeast side of the fault, but they are either much coarser or much finer than the arkosic section of the SAFOD rocks, and lack significant thicknesses of sand-bearing rock. Because they lie too far to the southeast, ~450 km southeast of SAFOD, a correlation is very unlikely. At the surface, there are no rocks similar to the SAFOD section nearby and northwest of SAFOD on the southwest side of the San Andreas fault. Subsurface data are scant; Upper Cretaceous or Paleogene rocks are reported to lie on Salinian rocks in wells 35–60 km to the northwest (Gribi, 1963; Forrest and Payne, 1964a; Rymer et al., 2006) (Fig. 14).

Model 2, which requires significant amounts of slip on both the San Andreas fault and the Buzzard Canyon fault, predicts that there are equivalent units to the southeast of the San Andreas fault midway between maximum and minimum possible offset distances (i.e., between 0 and 315 km), and as such could accommodate a broad area. There are a number of Paleocene–Eocene units distributed across the region from the San Andreas fault to the coast that are likely the remnants of a set of widespread deposits that are cut by the San Andreas, Sur-Nacimiento, and other strike-slip faults (Nilsen, 1987; Graham et al., 1989; Dickinson et al., 2005; Burnham, 2009). Rocks within these deposits might correlate with the SAFOD arkosic section, and depending on the correlation, up to 100 km of slip would be recorded on the Buzzard Canyon fault. The remainder of slip on the San Andreas fault modern strand would result in a correlation with early Tertiary arkosic rocks of the southern San Joaquin Valley–San Emigdio Mountain area, a slip of up to 200 km on the modern trace of the San Andreas fault (Fig. 14B). Model 2 appears to best explain the available data (Fig. 14D).

Finally, model 3 predicts that there has been only a small amount of offset on the modern San Andreas fault, and therefore rocks similar to the SAFOD arkoses lie somewhere nearby northeast of the San Andreas fault. Correlative arkosic rocks would lie ~300 km to the northwest on the southwest side of the San Andreas fault (Fig. 14C). To the northwest at this offset distance, no early Tertiary sedimentary rocks that resemble the SAFOD arkoses are present at the surface. The Tertiary sedimentary rocks east of the fault zone are the Miocene Monterey Formation (Dibblee, 1971, 1980; Sims, 1990; Dibblee et al., 1999), a sequence of argillaceous Miocene shale and mudstone (i.e., Dibblee, 1973; Sims, 1990), which do not resemble the SAFOD arkoses. Exposures of the Miocene Temblor Formation ~20 km southeast of the SAFOD site on the northeast side of the San Andreas fault (Sims, 1990) are only tens of meters thick and are volcanic-rich in composition. Likewise there are no reports of rocks similar to the SAFOD arkoses in well logs or imaged in seismic reflection profiles (Forrest and Payne, 1964b; Wentworth and Zoback, 1990) east or south of SAFOD. On the western side of the fault, the large displacement on the Buzzard Canyon fault would result in the equivalent rocks being faulted well to the north. Because Paleogene sedimentary rocks are relatively common northwest of SAFOD, and absent nearby to the southeast, we reject this model for the faulting of the arkosic rocks and the offset along the San Andreas fault system.

Our data, combined with the work of others in the region (Nilsen, 1984, 1986; Graham et al., 1989; Dickinson et al., 2005; Burnham, 2009) suggest that a model for partitioning of slip on the modern San Andreas fault and the Buzzard Canyon fault (model 2) provides a viable testable hypothesis for further work. The San Andreas fault consists of several strands from Parkfield northwestward to near Hollister, with small blocks of exotic rocks within the zone. Further work on the nature of these rocks, and provenance and petrographic analyses of sedimentary rocks in the SAFOD hole and northwest of SAFOD, would provide insight into the partitioning of displacement of the San Andreas fault zone in central California.

Summary

We infer that the faulted sedimentary arkosic section encountered in the SAFOD borehole, between the trace of the modern San Andreas fault and the Buzzard Canyon fault to the southwest, are latest Cretaceous to early Eocene rocks deposited in broad fan complexes proposed by Clarke and Nilsen (1973) and Nilsen and Clarke (1975). Both the continental borderland and granitic-cored Salinian terrains shed arkosic sediment westward and northward in a series of fan complexes (Nilsen and McKee, 1979; Graham et al., 1989) (see Fig. 14). Subsequent slip along faults of the San Andreas fault system have dissected and transported slivers within the fault zone northwestward 100–200 km. Similar complex structure and fault slivering has been mapped previously at a variety of scales from the Parkfield area northwestward up to 120 km from SAFOD (see Jennings and Strand, 1958; Dibblee, 1971; Rymer, 1981; Sims, 1993; Thayer and Arrowsmith, 2005).

CONCLUSIONS

We use wireline logs, cuttings microscopy, core samples, and electrical image logs to fully characterize the SAFOD arkosic rocks. As a whole, the 1230 m of arkosic rock between the modern-day San Andreas fault to the northeast and the Buzzard Canyon fault to the southwest is heavily faulted with at least 12 faults and 11 fault-bounded structural blocks. The arkosic section consists of three distinct lithologic units: the upper arkose, the clay-rich zone, and the lower arkose. The upper arkose is 156.4–381.4 m thick, interrupted by intraformational faults that separate five distinct structural blocks with different bedding orientations and fracture densities. The five blocks are coarse-grained quartzo-feldspathic conglomerates and sandstones that contain rare clay-rich grains. Feldspar altered to illitic material and iron oxide alteration products are abundant. The clay-rich zone is a 28.7–121.3 m thick, low-velocity zone that exhibits anomalously high gamma ray signatures. It is characterized by abundant clay-rich grains composed of high-birefringence clay minerals such as illite and smectite. A 27 m thick block of well-bedded, coarser-grained rock is incorporated into the clay-rich zone, from 2565 to 2595 mmd. The lower arkose has more variable log signatures than the upper arkose, with generally higher velocities and lower porosities. It is 224–331 m thick and is finer-grained with a higher abundance of clay-rich grains than the upper arkose. The lower arkose contains more volcanic clasts than the upper arkose.

We interpret the arkosic section to have been deposited in either a subaerial or subaqueous fan setting, perhaps as part of a transtensional subbasin system common to the Salinian block. Fission-track cooling ages of 64–70 Ma constrain the maximum age of the SAFOD arkoses to be latest Cretaceous to earliest Paleocene. Diagenetic mineral assemblages of authigenic muscovite and laumontite indicate that at the modern-day geothermal gradient of 38.5 °C/km the arkosic rocks have not likely been buried more than 0.8 km deeper than the depth they are today. Lithologic comparisons show that the rocks are no older than Eocene to latest Cretaceous.

The arkosic rocks are divided into blocks that are defined on the basis of common lithologies, bedding characteristics, and orientations that are cut by numerous faults and/or folded. We document the presence of the seismically inactive Buzzard Canyon fault, a western strand of the San Andreas fault, and an intermediate fault within the section, along with as many as eight other faults that define the boundaries of structural blocks. The three largest faults are associated with clay-rich mineralogy and locally high degree of imaged fractures and cataclastically deformed grains. We test several models for the partitioning of displacement across strands of the San Andreas fault system in the vicinity of the SAFOD site that predict where sedimentary rock units equivalent to the SAFOD arkoses would have lain before faulting. Based on our data and the tentative correlations to Paleocene–Eocene rocks of the region, we suggest that the SAFOD arkosic rocks are a stranded, fault-bounded sliver with correlative units 30–150 km to the northwest and >100 km to the southeast, on the northeast side of the fault.

This work was supported by NSF-Earthscope grant 06-043027 to Evans and grant EAR-0346190 to Kirschner and Garver. Additional funding to Springer is gratefully acknowledged. This work is part of a M.S. thesis funded in part by the Chevron Corporation, a DOSECC summer intern grant, and Springer and Williams endowments grants (Department of Geology, Utah State University). We thank Bill Ellsworth, Steve Hickman, and Mark Zoback for bringing the SAFOD project to fruition, and for their many discussions at the drill site and after drilling. We thank John Solum, Naomi Boness, Amy Day-Lewis, and Steve Graham for their invaluable insights into many aspects of the geology, mineralogy, and geophysics of the San Andreas fault, and for their generosity of time, and we thank Earl Brabb, Russ Graymer, Kellen Springer, and Carol Dehler for reviews of or discussions of this work. Many thanks to Laird Thompson for teaching us how to interpret image logs. We thank Robert Powell for his thoughtful and insightful review, and an anonymous reviewer for suggestions that helped improve the paper.

REFERENCES CITED

1.
Almeida
R.
Chester
J.S.
Chester
F.M.
Kirschner
D.L.
Moore
D.E.
Waller
T.D.
2005
,
Mesoscale structure and lithology of the SAFOD phase I and phase II core samples
:
Eos (Transactions, American Geophysical Union
 ), v.
86
,
Fall Meeting Supplement, abs. T21A-0454
.
2.
Arrowsmith
R.
2007
,
Structural geology, tectonic geomorphology, and recent slip history of the south-central San Andreas fault
:
Geological Society of America Abstracts with Programs
 , v.
39
,
no. 6
,
p
.
422
.
3.
Bachman
W.R.
Abbott
P.L.
1988
,
Lower Paleogene conglomerates in the Salinian block
,
in
Filewicz
M.V.
Squires
R.L.
eds
.,
Paleogene stratigraphy, West Coast of North America: Society of Economic Paleontologists and Mineralogists, Pacific Section
 , v.
58
,
p
.
135
150
.
4.
Barbeau
D.L.
Ducea
M.N.
Gehrels
G.E.
Kidder
S.
Wetmore
P.H.
Saleeby
J.B.
2005
,
U-Pb detrital-zircon geochronology of northern Salinian basement and cover rocks
:
Geological Society of America Bulletin
 , v.
117
,
p
.
466
481
,
doi: 10.1130/B25496.1
.
5.
Barth
A.P.
Wooden
J.L.
Grove
M.
Jacobson
C.E.
Pedrick
J.N.
2003
,
U-Pb zircon geochronology of rocks in the Salinas Valley region of California: A reevaluation of the crustal structure and origin of the Salinian block
:
Geology
 , v.
31
,
p
.
517
520
,
doi: 10.1130/0091-7613(2003)031<0517:UZGORI>2.0.CO;2
.
6.
Bent
J.V.
1988
,
Paleotectonics and provenance of Tertiary sandstones of the San Joaquin Basin, California
,
in
Graham
S.A.
ed
.,
Studies of the geology of the San Joaquin Basin: Society of Economic Paleontologists and Mineralogists, Pacific Section
 , v.
60
,
p
.
109
127
.
7.
Bernet
M.
Garver
J.I.
2005
,
Fission-track analysis of detrital zircon
,
in
Reiners
P.W.
Ehlers
T.A.
eds
.,
Low-temperature thermochronology: Techniques, interpretations, and applications: Reviews in Mineralogy and Geochemistry Series
 , v.
58
,
p
.
205
237
.
8.
Bleibinhaus
F.
Hole
J.A.
Ryberg
T.
Fuis
G.S.
2007
,
Structure of the California Coast Ranges and San Andreas fault at SAFOD from seismic waveform inversion and reflection imaging
:
Journal of Geophysical Research
 , v.
112
,
B06315
,
doi: 10.1029/2006JB004611
.
9.
Blythe
A.E.
d'Alessio
M.A.
Burgmann
R.
2004
,
Constraining the exhumation and burial history of the SAFOD pilot hole with apatite fission track and (U-Th)/He thermochronometry
:
Geophysical Research Letters
 , v.
31
,
L15S16
,
doi: 10.1029/2003GL019407
.
10.
Boness
N.L.
Zoback
M.D.
2006
,
A multiscale study of the mechanisms controlling shear velocity anisotropy in the San Andreas Fault Observatory at Depth
:
Geophysics
 , v.
71
,
F131
,
doi: 10.1190/1.2231107
.
11.
Bradbury
K.K.
Barton
D.C.
Draper
S.D.
Solum
J.G.
Evans
J.P.
2007
,
Mineralogic and textural analyses of drill cuttings from the San Andreas Fault Observatory at Depth (SAFOD) boreholes: Initial interpretations of fault zone composition and constraints on geologic models
:
Geosphere
 , v.
3
,
p
.
299
318
,
doi: 10.1130/GES00076.1
.
12.
Brandon
M.T.
1996
,
Probability density plot for fission track grain-age samples
:
Radiation Measurements
 , v.
26
,
p
.
663
676
,
doi: 10.1016/S1350-4487(97)82880-6
.
13.
Burnham
K.
2009
,
Predictive model of San Andreas fault system paleogeography, Late Cretaceous to early Miocene, derived from detailed multidisciplinary conglomerate correlations
:
Tectonophysics
 , v.
464
,
p
.
195
258
,
doi: 10.1016/j.tecto.2007.11.056
.
14.
Carozzi
A.V.
1993
,
Sedimentary petrography
 :
Englewood Cliffs, New Jersey
,
Prentice Hall
,
263
p
.
15.
Catchings
R.D.
Rymer
M.J.
2002
,
High-resolution seismic velocities and shallow structure of the San Andreas fault zone at Middle Mountain, Parkfield, California
:
Bulletin of the Seismological Society of America
 , v.
92
,
p
.
2493
2503
,
doi: 10.1785/0120010263
.
16.
Catchings
R.D.
Rymer
M.J.
Goldman
M.R.
2007
,
Internal structure of the San Andreas fault from seismic imaging and implications for damage zone width
:
Seismological Research Letters
 , v.
78
,
no. 2
,
p
.
307
.
17.
Chipping
D.H.
1972
,
Early Tertiary paleogeography of central California
:
American Association of Petroleum Geologists Bulletin
 , v.
56
,
p
.
480
493
.
18.
Clark
J.C.
1981
,
Stratigraphy, paleontology, and geology of the central Santa Cruz Mountains, California Coast Ranges
:
U.S. Geological Survey Professional Paper
1168
,
51
p
.
19.
Clarke
S.H.
1973
,
The Eocene Point of Rocks sandstone: Provenance, mode of deposition, and implications for the history of offset along the San Andreas fault in central California [Ph.D. thesis]
:
Berkeley, California
,
University of California
,
302
p
.
20.
Clarke
S.H.
Nilsen
T.H.
1973
,
Displacement of Eocene strata and implication for the history of offset of the San Andreas fault, central and northern California
,
in
Kovach
R.L.
Nur
A.
eds
.,
Proceedings of the Conference on Tectonic Problems of the San Andreas Fault System: Stanford University Publications in Geological Sciences
 , v.
13
,
p
.
358
367
.
21.
Clifton
H.E.
1981
,
Submarine canyon deposits, Point Lobos, California
,
in
Frizzell
V.
ed
.,
Upper Cretaceous and Paleocene turbidite deposits, central California coast: Society of Economic Paleontologists and Mineralogists, Pacific Section, Field Trip Guidebook
 , v.
20
,
p
.
79
92
.
22.
Clifton
H.E.
1984
,
Sedimentation units in stratified deep-water conglomerate, Paleocene submarine canyon fill, Point Lobos, California
,
in
Koster
E.H.
Steele
R.J.
eds
.,
Sedimentology of gravels and conglomerates
 :
Canadian Society of Petroleum Geologists Memoir 10, p
.
429
441
.
23.
Clifton
H.E.
2007
,
Conglomeratic submarine fill, Point Lobos Reserve
:
SEPM (Society for Sedimentary Geology), Pacific Section, Field Trip Guidebook
 , v.
105
,
79
p
.
24.
Critelli
S.
Nilsen
T.H.
1996
,
Petrology and diagenesis of the Eocene Butano Sandstone, La Honda Basin, California
:
Journal of Geology
 , v.
104
,
p
.
295
315
.
25.
Critelli
S.
Nilsen
T.H.
2000
,
Provenance and stratigraphy of the Eocene Tejon Formation, western Tehachapi Mountains, San Emigdio Mountains, and southern San Joaquin Basin, California
:
Sedimentary Geology
 , v.
136
,
p
.
7
27
,
doi: 10.1016/S0037-0738(00)00080-4
.
26.
Cronin
B.T.
Kidd
R.B.
1998
,
Heterogeneity and lithotype distribution in ancient deep-sea canyons: Point Lobos deep-sea canyon as a reservoir analog
:
Sedimentary Geology
 , v.
115
,
p
.
315
349
,
doi: 10.1016/S0037-0738(97)00099-7
.
27.
Crowell
J.C.
Susuki
T.
1959
,
Eocene stratigraphy and paleontology, Orocopia Mountains, southeastern California
:
Geological Society of America Bulletin
 , v.
70
,
p
.
581
592
,
doi: 10.1130/0016-7606(1959)70[581:ESAPOM]2.0.CO;2
.
28.
Dapples
E.C.
1979
,
Diagenesis of sandstones
,
in
Larsen
G.
Chilinger
G.V.
eds
.,
Developments in sedimentology 8: Diagenesis in sediments
 :
Amsterdam
,
Elsevier
,
p
.
91
127
.
29.
DeGraaff-Surpless
K.
Graham
S.A.
Wooden
J.L.
McWilliams
M.O.
2002
,
Detrital zircon provenance analysis of the Great Valley Group, California: Evolution of an arc-forearc system
:
Geological Society of America Bulletin
 , v.
114
,
p
.
1564
1580
,
doi: 10.1130/0016-7606(2002)114<1564:DZPAOT>2.0.CO;2
.
30.
Dibblee
T.W.
Jr.
1971
,
Geologic maps of 17 15-minute quadrangles along the San Andreas fault in the vicinity of King City, Coalinga, Panoche Valley and Paso Robles, with index map
:
U.S. Geological Survey Open-File Report 71–87, scale 1:62,500
.
31.
Dibblee
T.W.
Jr.
1973
,
Stratigraphy of the southern Coast Ranges near the San Andreas fault from Cholame to Maricopa, California
:
U.S. Geological Survey Professional Paper
764
,
45
p
.
32.
Dibblee
T.W.
Jr.
1980
,
Geology along the San Andreas fault from Gilroy to Parkfield
,
in
Streitz
R.
Sherburne
R.
eds
.,
Studies of the San Andreas fault zone in northern California
 :
California Division of Mines and Geology Special Paper
140
,
p
.
3
18
.
33.
Dibblee
T.W.
Graham
S.E.
Mahony
T.M.
Blissenbach
J.L.
Mariant
J.J.
Wentworth
C.M.
1999
,
Regional geologic map of San Andreas and related faults in Carrizo Plain, Temblor, Caliente, and La Panza Ranges and vicinity, California: A digital database
:
U.S. Geological Survey Open-File Report 99–14
.
34.
Dickinson
W.R.
1966
,
Structural relationships of San Andreas fault system, Cholame Valley and Castle Mountain Range, California
:
Geological Society of America Bulletin
 , v.
77
,
p
.
707
725
,
doi: 10.1130/0016-7606(1966)77[707:SROSAF]2.0.CO;2
.
35.
Dickinson
W.R.
1970
,
Interpreting detrital modes of graywacke and arkose
:
Journal of Sedimentary Petrology
 , v.
40
,
p
.
695
707
.
36.
Dickinson
W.R.
Butler
R.F.
1998
,
Coastal and Baja California paleomagnetism reconsidered
:
Geological Society of America Bulletin
 , v.
110
,
p
.
1268
1280
,
doi: 10.1130/0016-7606(1998)110<1268:CABCPR>2.3.CO;2
.
37.
Dickinson
W.R.
Ducea
M.
Rosenberg
L.I.
Greene
H.G.
Graham
S.A.
Clark
J.C.
Weber
G.E.
Kidder
S.
Ernst
G.W.
Brabb
E.E.
2005
,
Net dextral slip, Neogene San Gregorio–Hosgri fault zone, coastal California: Geologic evidence and tectonic implications
:
Geological Society of America Special Paper
391
,
43
p
.
38.
Durham
D.L.
1974
,
Geology of the southern Salinas Valley area, California
:
U.S. Geological Survey Professional Paper
819
,
111
p
.
39.
Ekstrom
M.P.
Dahan
C.A.
Chen
M.Y.
Lloyd
P.M.
Rossi
D.J.
1986
,
Formation imaging with microelectrical scanning arrays
:
Society of Professional Well Log Analysts 27th Annual Logging Symposium, paper BB
.
40.
Erzinger
J.
Wiersberg
T.
Zimmer
M.
2006
,
Real-time mud gas logging and sampling during drilling
:
Geofluids
 , v.
6
,
no. 3
,
p
.
225
233
.
41.
Evans
J.P.
Moore
D.E.
Kirschner
D.
Solum
J.G.
2005
,
Lithologic characterization of the deep portion of the SAFOD drillhole
:
Eos (Transactions, American Geophysical Union
 ), v.
86
,
Fall Meeting Supplement, abs.T21A-0450
.
42.
Falk
P.D.
Dorsey
R.J.
1998
,
Rapid development of gravelly high-density turbidity currents in marine Gilbert-type fan deltas, Loreto Basin, Baja California Sur, Mexico
:
Sedimentology
 , v.
45
,
p
.
331
349
,
doi: 10.1046/j.1365-3091.1998.0153e.x
.
43.
Forrest
L.
Payne
M.B.
1964a
,
San Andreas fault cross sections parallel along the east and west sides of the fault from Parkfield to King City
,
in Pacific Section of the American Association of Petroleum Geologists Committee for the Study of Lateral Faulting in California, Maps and Cross Sections
.
44.
Forrest
L.
Payne
M.B.
1964b
,
San Andreas fault cross sections parallel along the east and west sides of the fault from Parkfield to Soda Lakes
,
in Pacific Section of the American Association of Petroleum Geologists Committee for the Study of Lateral Faulting in California, Maps and Cross Sections
.
45.
Garver
J.I.
Kamp
P.J.
2002
,
Integration of zircon color and zircon fission track zonation patterns in orogenic belts: Application of the Southern Alps, New Zealand
:
Tectonophysics
 , v.
349
,
p
.
203
219
,
doi: 10.1016/S0040-1951(02)00054-9
.
46.
Garver
J.I.
Reiners
P.R.
Walker
L.J.
Ramage
J.M.
Perry
S.E.
2005
,
Implications for timing of Andean uplift based on thermal resetting of radiation-damaged zircon in the Cordillera Huayhuash, northern Perú
:
Journal of Geology
 , v.
113
,
p
.
117
138
,
doi: 10.1086/427664
.
47.
Graham
S.A.
1978
,
Role of Salinian block in the evolution of the San Andreas fault system, California
:
American Association of Petroleum Geologists Bulletin
 , v.
62
,
p
.
2214
2231
.
48.
Graham
S.A.
Stanley
R.G.
Bent
J.V.
Carter
J.B.
1989
,
Oligocene and Miocene paleogeography of central California and displacement along the San Andreas fault
:
Geological Society of America Bulletin
 , v.
101
,
p
.
711
730
,
doi: 10.1130/0016-7606(1989)101<0711:OAM-POC>2.3.CO;2
.
49.
Gribi
E.A.
Jr.
1963
,
Bitterwater oil field, San Benito County, California: Guidebook to the geology of Salinas Valley and the San Andreas fault
:
American Association of Petroleum Geologists Annual Field Trip Guidebook, p
.
74
75
.
50.
Grove
K.
1993
,
Latest Cretaceous basin formation within the Salinian terrane of west-central California
:
Geological Society of America Bulletin
 , v.
105
,
p
.
447
463
,
doi: 10.1130/0016-7606(1993)105<0447:LCBFWT>2.3.CO;2
.
51.
Hickman
S.
Zoback
M.D.
Ellsworth
W.
2004
,
Introduction to special section: Preparing for the San Andreas Fault Observatory at Depth
:
Geophysical Research Letters
 , v.
31
,
L12S01
,
doi: 10.1029/2004GL020688
.
52.
Hickman
S.
Zoback
M.
Ellsworth
B.
2005
,
Structure and composition of the San Andreas fault zone at Parkfield: Initial results from SAFOD phases 1 and 2
:
Eos (Transactions, American Geophysical Union
 ), v.
86
,
Fall Meeting Supplement, abs. T23E-05
.
53.
Hickman
S.
Zoback
M.
Ellsworth
W.
Chester
J.
Chester
F.
Evans
J.
Moore
D.
Kirschner
D.
Schleicher
A.
van der Pluijm
B.
Solum
J.
2008
,
Structure and composition of the San Andreas fault in central California: Recent results from SAFOD sample analyses
:
Eos (Transactions, American Geophysical Union
 ), v.
89
,
Fall Meeting Supplement, abs. T53F-01
.
54.
Hole
J.A.
Catchings
R.D.
St. Clair
K.C.
Rymer
M.J.
Okaya
D.A.
Carney
B.J.
2001
,
Steep-dip seismic imaging of the shallow San Andreas fault near Parkfield
:
Science
 , v.
294
,
p
.
1513
1515
,
doi: 10.1126/science.1065100
.
55.
Hole
J.A.
Ryberg
T.
Fuis
G.S.
Bleibinhaus
F.
Sharma
A.K.
2006
,
Structure of the San Andreas fault zone at SAFOD from a seismic refraction survey
:
Geophysical Research Letters
 , v.
33
,
L07312
,
doi: 10.1029/2005GL025194
.
56.
Huang
W.L.
1992
,
Illitic-clay formation during experimental diagenesis of arkoses
,
in
Houseknecht
D.W.
Pittman
E.D.
eds
.,
Origin, diagenesis, and petrophysics of clay minerals in sandstones
 :
SEPM (Society for Sedimentary Geology) Special Publication
47
,
p
.
49
63
.
57.
Hurford
A.J.
Carter
A.
1991
,
The role of fission track dating in discrimination of provenance
,
in
Morton
A.D.
Todd
S.P.
Haughton
P.D.W.
eds
.,
Developments in sedimentary provenance studies
 :
Geological Society, London, Special Publication
57
,
p
.
67
78
.
58.
Ingersoll
R.V.
1988
,
Development of the Cretaceous forearc basin on central California
,
in
Graham
S.A.
ed
.,
Studies of the geology of the San Joaquin Basin: Society of Economic Paleontologists and Mineralogists, Pacific Section
 , v.
60
,
p
.
141
155
.
59.
Irwin
W.P.
1990
,
Geology and plate tectonic development
,
in
Wallace
R.E.
ed
.,
The San Andreas fault system
 :
U.S. Geological Survey Professional Paper
1515
,
p
.
61
80
.
60.
Jennings
C.W.
Strand
R.G.
1958
,
Geologic map of California, Santa Cruz sheet
:
California Division of Mines and Geology, scale 1:250,000
.
61.
Jennings
C.W.
Strand
R.G.
1976
,
Geologic map of California, San Luis Obispo sheet
:
California Division of Mines and Geology, scale 1:250,000
.
62.
Kidder
S.
Ducea
M.N.
Gehrels
G.E.
Patchett
P.J.
Vervoot
J.
2003
,
Tectonic and magmatic development of the Salinian Coast Ridge Belt, California
:
Tectonics
 , v.
22
,
1058
,
doi: 10.1029/2002TC001409
.
63.
Kistler
R.W.
Champion
D.E.
2001
,
Rb-Sr whole rock and mineral ages, K-Ar, 40Ar/39Ar, and U-Pb mineral ages, and strontium, lead, neodymium, and oxygen isotopic compositions for granitic rocks from the Salinian Composite Terrane, California
:
U.S. Geological Survey Open-File Report 01-453
,
84
p
.
64.
Li
Y.-G.
Malin
P.E.
2008
,
San Andreas fault damage at SAFOD viewed with fault-guided waves
:
Geophysical Research Letters
 , v.
35
,
L08304
,
doi: 10.1029/2007GL032924
.
65.
Linn
A.M.
1991
,
Isotopic and geochemical stratigraphy and paleotectonic analysis: Mesozoic Great Valley forearc basin of California [Ph.D. thesis]
:
Los Angeles
,
University of California
,
133
p
.
66.
Liu
K.W.
2002
,
Deep burial diagenesis of the siliciclastic Ordovician Natal Group, South Africa
:
Sedimentary Geology
 , v.
154
,
p
.
177
189
,
doi: 10.1016/S0037-0738(02)00161-6
.
67.
Matthews
V.
III
1976
,
Correlation of the Pinnacles and Neenach volcanic formations and their bearing on the San Andreas fault problem
:
American Association of Petroleum Geologists Bulletin
 , v.
60
,
p
.
2128
2141
.
68.
McPhee
D.K.
Jachens
R.C.
Wentworth
C.M.
2004
,
Crustal structure across the San Andreas fault at the SAFOD site from potential field and geologic studies
:
Geophysical Research Letters
 , v.
31
,
L12S03
,
doi: 10.1029/2003GL019363
.
69.
Montoya
J.W.
Hemley
J.J.
1975
,
Activity relations and stabilities in alkali feldspar and mica alteration reactions
:
Economic Geology and the Bulletin of the Society of Economic Geologists
 , v.
70
,
p
.
577
583
.
70.
Moore
D.M.
Reynolds
R.C.
1989
,
X-ray diffraction and the identification and analysis of clay minerals
 :
Oxford University Press
,
378
p
.
71.
Nilsen
T.H.
1984
,
Offset along the San Andreas fault of Eocene strata from the San Juan Bautista area and western San Emigdio Mountains, California
:
Geological Society of America Bulletin
 , v.
95
,
p
.
599
609
,
doi: 10.1130/0016-7606(1984)95<599:OATSAF>2.0.CO;2
.
72.
Nilsen
T.H.
1986
,
Cretaceous Paleogeography of western North America
,
in
Abbott
P.L.
ed
.,
Cretaceous stratigraphy, western North America
 :
Society of Economic Paleontologists and Mineralogists, Pacific Section, Field Trip Guidebook
, v.
46
,
p
.
1
39
.
73.
Nilsen
T.H.
1987
,
Stratigraphy and sedimentology of the Eocene Tejon Formation, western Tehachapi and San Emigdio Mountains, California
:
U.S. Geological Survey Professional Paper
1268
,
110
p
.
74.
Nilsen
T.H.
Clarke
S.H.
1975
,
Sedimentation and tectonics in the early Tertiary continental borderland of central California
:
U.S. Geological Survey Professional Paper
925
,
64
p
.
75.
Nilsen
T.H.
McKee
E.H.
1979
,
Cenozoic paleogeography of the western United States
,
in
Armentrout
J.M.
Cole
M.R.
TerBest
H.
Jr.
eds., Pacific Coast Paleogeography Symposium 3, Anaheim, California, United States, p
.
257
276
.
76.
Page
B.M.
Thompson
G.A.
Coleman
R.G.
1998
,
Late Cenozoic tectonics of the central and southern Coast Ranges of California
:
Geological Society of America Bulletin
 , v.
110
,
p
.
846
876
,
doi: 10.1130/0016-7606(1998)110<0846:OLCTOT>2.3.CO;2
.
77.
Powell
R.E.
1993
,
Balanced palinspastic reconstruction of pre-late Cenozoic paleogeography, southern California: Geologic and kinematic constraints on evolution of the San Andreas fault system
,
in
Powell
R.E.
Weldon
R.J.
II
Matti
J.C.
eds
.,
The San Andreas fault system: Displacement, palinspastic reconstruction, and geologic evolution
 :
Geological Society of America Memoir 178, p
.
1
106
.
78.
Rogers
T.H.
Church
C.C.
1967
,
Upper Cretaceous rocks southwest of Hollister
,
in Guidebook to the Gabilan Range and adjacent San Andreas fault: American Association of Petroleum Geologists and Society of Economic Paleontologists and Mineralogists, Pacific Sections, p
.
45
47
.
79.
Ross
D.C.
1970
,
Quartz gabbro and anorthositic gabbro: Markers of offset along the San Andreas fault in the California Coast Ranges
:
Geological Society of America Bulletin
 , v.
81
,
p
.
3647
3661
,
doi: 10.1130/0016-7606(1970)81[3647:QGAAGM]2.0.CO;2
.
80.
Ross
D.C.
1984
,
Possible correlations of basement rocks across the San Andreas, San Gregorio–Hosgri, and Rinconada–Reliz–King City faults, California
:
U.S. Geological Survey Professional Paper
1317
,
37
p
.
81.
Rymer
M.J.
1981
,
Geologic map along a 12-kilometer segment of the San Andreas fault zone, southern Diablo Range, California
:
U.S. Geological Survey Open-File Report 81–1173, scale 1:12,000
.
82.
Rymer
M.J.
Catchings
R.D.
Goldman
M.R.
2003
,
Structure of the San Andreas fault zone as revealed by surface geologic mapping and high-resolution seismic profiling near Parkfield, California
:
Eos (Transactions, American Geophysical Union
 ), v.
85
,
Fall Meeting Supplement, abs. T-11F-08
.
83.
Rymer
M.J.
Hickman
S.H.
Stoffer
P.W.
2006
,
A field guide to the central, creeping section of the San Andreas fault and the San Andreas Fault Observatory at Depth
,
in
Prentice
C.S.
Scotchmoor
J.G.
Moores
E.M.
Kiland
J.P.
eds
.,
1906 San Francisco earthquake centennial field guides
 :
Field trips associated with the 100th anniversary conference: Geological Society of America Field Guide 7, p
.
237
272
,
doi: 10.1130/2006.1906SF(16)
.
84.
Sass
J.H.
Williams
C.F.
Lachenbruch
A.H.
Galanis
S.P.
Jr.
Grubb
F.V.
1997
,
Thermal regime of the San Andreas fault near Parkfield, California
:
Journal of Geophysical Research
 , v.
102
,
p
.
27,575
27,585
.
85.
Schleicher
A.M.
van der Pluijm
B.A.
Solum
J.G.
Warr
L.N.
2006
,
The origin and significance of clay-coated fractures in mudrock fragments of the SAFOD borehole (Parkfield, California)
:
Geophysical Research Letters
 , v.
33
,
L16313
,
doi: 10.1029/2006GL026505
.
86.
Scholle
P.A.
1979
,
A color illustrated guide to constituents, textures, cements, and porosities of sandstones and associated rocks
:
American Association of Petroleum Geologists Memoir 28
,
201
p
.
87.
Seiders
V.M.
Cox
B.F.
1992
,
Place of origin of the Salinian block, California, as based on clast compositions of Upper Cretaceous and lower Tertiary conglomerates
:
U.S. Geological Survey Professional Paper
1526
,
80
p
.
88.
Sims
J.D.
1990
,
Geologic map of the San Andreas fault in the Parkfield 7.5-minute quadrangle, Monterey and Fresno Counties, California
:
U.S. Geological Survey Miscellaneous Field Studies Map MF-2115, scale 1:24000
.
89.
Sims
J.D.
1993
,
Chronology of displacement on the San Andreas fault in central California: Evidence from reversed positions of exotic rock bodies near Parkfield, California
,
in
Powell
R.E.
Weldon
R.J.
II
Matti
J.C.
eds
.,
The San Andreas fault system: Displacement, reconstruction, and geologic evolution
 :
Geological Society of America Memoir 178, p
.
231
256
.
90.
Solum
J.G.
Hickman
S.H.
Lockner
D.A.
Moore
D.E.
van der Pluijm
B.A.
Schleicher
A.M.
Evans
J.P.
2006
,
Mineralogical characterization of protolith and fault rocks from the SAFOD main hole
:
Geophysical Research Letters
 , v.
33
,
L21314
,
doi: 10.1029/2006GL027285
.
91.
Tembe
S.
Lockner
D.A.
Solum
J.G.
Morrow
C.A.
Wong
T.
Moore
D.E.
2005
,
Strength of the San Andreas fault zone: Insight from cuttings and core
:
Eos (Transactions, American Geophysical Union
 ), v.
86
,
Fall Meeting Supplement, abstract T24B-03
.
92.
Thayer
M.
Arrowsmith
R.
2005
,
Fault zone structure of Middle Mountain, central California
:
Eos (Transactions, American Geophysical Union
 ), v.
86
,
Fall Meeting Supplement, abs. T21A-0458
.
93.
Thompson
L.B.
ed
.,
2000
,
Atlas of borehole imagery
 :
Tulsa, Oklahoma
,
American Association of Petroleum Geologists Discovery Series
4
,
CD-ROM
.
94.
Thurber
C.
Roecker
S.
Zhang
H.
Baher
S.
Ellsworth
W.
2004
,
Fine-scale structure of the San Andreas fault zone and location of the SAFOD target earthquakes
:
Geophysical Research Letters
 , v.
31
,
L12S02
,
doi: 10.1029/2003GL019398
.
95.
Vedder
J.G.
Brown
R.D.
1968
,
Structural and stratigraphic relations along the Nacimiento fault in the southern Santa Lucia Range and San Rafael Mountains, California
,
in
Dickinson
W.R.
Grantz
A.
eds
.,
Proceedings of Conference on Geologic Problems of the San Andreas Fault System: Stanford University Publications of the Geological Sciences
 , v.
11
,
p
.
243
258
.
96.
Vermeesch
P.
Miller
D.D.
Graham
S.A.
de Grave
J.
McWilliams
M.O.
2006
,
Multimethod detrital thermochronology of the Great Valley Group near New Idria, California
:
Geological Society of America Bulletin
 , v.
118
,
p
.
210
218
,
doi: 10.1130/B25797.1
.
97.
Wagner
G.A.
Van der Haute
P.
1992
,
Fission track dating
 :
Dordrecht and Boston
,
Kluwer
,
285
p
.
98.
Walrond
H.
Gribi
E.A.
Jr.
1963
,
Geologic map to accompany the guidebook to the geology of Salinas Valley and the San Andreas fault
:
American Association of Petroleum Geologists and Society of Economic Paleontologists and Mineralogists, Pacific Sections, Joint Annual Field Trip Guidebook, scale 1:126,720
.
99.
Walrond
H.
Thourp
R.R.
Gribi
E.A.
Jr.
Rogers
T.H.
1967
,
Geologic map of the Gabilan uplift and adjacent areas, Monterey and San Benito Counties
,
in Gabilan Range and adjacent San Andreas fault: American Association of Petroleum Geologists and Society of Economic Paleontologists and Mineralogists, Pacific Sections, Joint Annual Field Trip Guidebook, scale 1:126,720
.
100.
Wentworth
C.M.
Zoback
M.D.
1990
,
Structure of the Coalinga area and thrust origin of the earthquake
,
in
Rymer
M.J.
Ellsworth
W.L.
eds
.,
The Coalinga California earthquake of May 2, 1983
 :
U.S. Geological Survey Professional Paper
1487
,
p
.
41
68
.
101.
Whidden
K.J.
Lund
S.P.
Bottjer
D.J.
Champion
D.
Howell
D.G.
1998
,
Paleomagnetic evidence that the central block of Salinia (California) is not a far-traveled terrane
:
Tectonics
 , v.
17
,
no. 3
,
p
.
329
344
,
doi: 10.1029/97TC03021
.
102.
Winter
H.
Adelhardt
S.
Jerak
A.
Kuchenhoff
H.
2002
,
Characterization of cataclastic shear-zones of the KTB deep drill hole by regression analysis of drill cuttings data
:
Geophysical Journal International
 , v.
150
,
p
.
1
9
,
doi: 10.1046/j.1365-246X.2002.01738.x
.
103.
Wood
D.J.
Saleeby
J.B.
1997
,
Late Cretaceous–Paleocene extensional collapse and disaggregation of the southernmost Sierra Nevada batholith
:
International Geology Review
 , v.
39
,
p
.
973
1009
.
104.
Worden
R.H.
Morad
S.
2003
,
Clay minerals in sandstones: Controls on formation, distribution and evolution
,
in
Worden
R.H.
Morad
S.
eds
.,
Clay mineral cements in sandstones
 :
International Association of Sedimentologists Special Publication
34
,
p
.
3
42
.
105.
Wright
T.A.
2000
,
The evolving San Andreas fault system: Nascent tectonic succession of California?
in
Bokelmann
G.
Kovach
R.L.
eds
.,
Proceedings of the 3rd Conference on Tectonic Problems of the San Andreas Fault System
 :
Palo Alto, California
,
Stanford University Publications
,
p
.
144
149
.
1GSA Data Repository item 2009160, Appendices A–D, is available at www.geosociety.org/pubs/ft2009.htm, or on request from editing@geosociety.org, Documents Secretary, GSA, P.O. Box 9140, Boulder, CO 80301-9140, USA.