The relative contributions of tectonic and magmatic processes to continental rifting are highly variable. Magnetic, gravity, and seismic reflection data from Surprise Valley, California, in the northwest Basin and Range, reveal an intrabasin, fault-controlled, ∼10-m-thick dike at a depth of ∼150 m, providing an excellent example of the interplay between faulting and dike intrusion. The dike, likely a composite structure representing multiple successive intrusions, is inferred from modeling a positive magnetic anomaly that extends ∼35 km and parallels the basin-bounding Surprise Valley normal fault on the west side of the valley. A two-dimensional high-resolution seismic reflection profile acquired across the magnetic high images a normal fault dipping 56°E with ∼275 m of throw buried ∼60 m below the surface. Densely spaced gravity measurements reveal a <1 mGal gravity low consistent with the fault offset inferred from the seismic data. Collinearity of the magnetic high and gravity low for ∼6 km implies normal fault control of the dike along that length. The unusually shallow angle of the dike suggests that motion along the fault (perhaps aided by reduced friction along the dike) and associated block rotation resulted in post-intrusion tilting of the dike. The source of the dike is likely related to a shallow brittle-ductile transition zone that was elevated following rapid slip on the Surprise Valley fault after 3 Ma. Prior to our work, the Surprise Valley fault was assumed to accommodate the vast majority of extension across the region. Our results indicate that subsurface features, although no longer active, are significant contributors to the processes, timing, and total amount of extension observed in continental rift environments.

In continental rift zones such as the Basin and Range or East African Rift, tectonic extension in the seismogenic crust can be accommodated by normal faulting or magmatism. In active systems, satellite geodesy coupled with studies of seismicity can identify events in which rift-related extension is accommodated by diking (Wright et al., 2006; Pallister et al., 2010), seismic slip and aseismic deformation (Payne et al., 2008; Bell et al., 2012), or a combination (Calais et al., 2008; Biggs et al., 2009). When extension is accommodated by both magmatic and tectonics processes, the proportion of strain accommodated by each process as well as the spatial and temporal distribution of the two processes are highly variable.

In systems that are no longer active, geologic mapping can assess the spatial and temporal distribution of strain over a longer time period. For example, in Mono Basin (California) Bursik and Sieh (1989) mapped fault scarps and measured offset along them, and compared the timing of events with eruption of the Mono Craters; on the basis of these data sets, they hypothesized that extension was accommodated by normal faulting prior to 40 ka, and was supplanted by dike intrusion since then. Parsons and Thompson (1991) expanded this hypothesis into the axiom that normal faulting and magmatism work together to accommodate strain in proportion with the magma supply. Their model suggests that when there is sufficient magma supply, magmatic intrusion suppresses normal faulting. While this may be the case in a broad sense (e.g., Parsons et al., 1998), the relationship between magmatism and faulting appears to be more complicated at the scale of individual events or on shorter time scales. For example, Valentine and Krogh (2006) hypothesized that dikes intruding into preexisting faults may actually promote slip along those faults by reducing friction.

Determining the total amount and distribution of strain, as well as the relative timing and contributions of dike intrusion and normal fault slip, requires looking beyond surface exposures. The basins of the Basin and Range hide a significant portion of the deformational history of the region; several geophysical studies have shown that normal faults with significant offset are buried beneath alluvium (Langenheim et al., 2001; Grauch and Hudson, 2007; Blackwell et al., 2009). It can be particularly difficult to assess the role of magmatic intrusions such as dikes in the subsurface because their narrow, vertical form makes them essentially invisible to seismic reflection profiles. Potential field data can provide more insight in many areas where mafic dikes present a significant contrast in both density and magnetic properties with the surrounding rocks. In combination, seismic and potential field data are capable of producing detailed maps and models of the subsurface that facilitate a more complete strain analysis than can be determined from surface mapping.

The northwestern margin of the Basin and Range has been the subject of numerous geological studies in the past decade (Colgan et al., 2006; Meigs et al., 2009; Scarberry et al., 2010; Egger and Miller, 2011). In the Surprise Valley region in particular, 12%–15% extension over ∼50 km has been documented through surface mapping of exposed faults (Egger and Miller, 2011). The estimated extension, however, did not take into account strain buried in the basin. The potential field and seismic reflection work presented here suggests that significant structures are hidden beneath the sediments of the basin, indicating not only more extension than has been previously estimated, but that a portion of that extension is occurring through dike intrusion. In addition, these intrabasin structures suggest a complex interaction between developing normal faults and dikes that may be influenced by a shallow brittle-ductile transition zone.

Surprise Valley is an extensional basin located along the northwestern margin of the Basin and Range Province (Fig. 1). It is bound on the west by the Surprise Valley fault, which has accommodated ∼8 km of dip-slip motion since the middle Miocene to expose a sequence of Eocene and younger volcanic and volcaniclastic rocks in the Warner Range (Egger and Miller, 2011). The northern part of the valley, referred to here as the upper basin, consists of a half-graben bound on the west by the Surprise Valley fault and by the Larkspur Hills to the east (Fig. 1). The Larkspur Hills consist of late Miocene–Pliocene (8–3 Ma) low-potassium, high-alumina olivine tholeiites interbedded with tuffs and tuffaceous sediments (Tbl), a sequence that crops out extensively in northeastern California and southern Oregon (Hart et al., 1984; Carmichael et al., 2006). Arc-derived Oligocene volcanic rocks (Tv), exposed to the south in the Hays Canyon Range and in the Warner Range to the west, underlie the basalts (Colgan et al., 2011).

The interbedded basalt flows and lake sediments that fill the basin generate strong contrasts ideally suited for geophysical methods. A potential field model was developed (Egger et al., 2010) along a seismic reflection profile (Fig. 1) acquired by Lerch et al. (2010); both the seismic profile and potential field model suggest the presence of numerous intrabasin faults with offsets to tens of meters. These buried faults are likely analogous to faults exposed east of the valley in the Larkspur Hills (Fig. 1), where numerous small-offset, east-dipping normal faults cut ca. 8–3 Ma basalts (Tbl), with a total extension of 5%–7% over 10 km (Strickley, 2014). These faults are no longer active, however, based on laterally continuous paleoshorelines that formed ca. 0.02 Ma (Egger, 2014). Strickley (2014) also mapped several linear volcanic vents that parallel fault trends, but did not observe any dikes. Ritzinger (2014) used paleomagnetic data to distinguish six distinct flow groups that were spatially controlled by normal faults, suggesting that normal faults developed concurrently (perhaps episodically) with volcanism.

Within Surprise Valley, a significant positive magnetic anomaly was imaged (Glen et al., 2013) that roughly parallels the orientation of the Surprise Valley fault, but is straighter and longer than any segment of it or any individual fault within the Larkspur Hills (Fig. 1B). Here we use potential field modeling that integrates data from gravity and magnetic profiles (Athens, 2011) with a high-resolution seismic reflection profile (Fontiveros, 2010) to examine the source of the magnetic high (Fig. 1).

Magnetic Data

Ground-based magnetic measurements were collected (Athens, 2011; Glen et al., 2013) using a cesium-vapor magnetometer mounted on or towed behind all-terrain vehicles (Athens et al., 2011) as well as a backpack-mounted system. Processing steps included subtracting diurnal variations of the Earth’s magnetic field, removing aberrant data points (either due to sensor error or cultural artifacts), and removing the International Geomagnetic Reference Field to derive the residual magnetic anomaly field (Fig. 1B). The data density is greatest in the upper basin and the central portion of the middle basin due to easy access to the playa and better terrain for the all-terrain vehicles. Detailed magnetic profiles in the upper basin are shown in Figure 2A.

Gravity Data

A detailed gravity survey complemented the magnetic survey in the upper basin; 313 gravity stations were collected along 17 east-west transects across the magnetic anomaly, with station spacing ranging from 50 to 250 m (Fig. 2B). Gravity stations were tied to a primary base station in Alturas, California (Jablonski, 1974), and were reduced using standard gravity methods that include Earth-tide, instrument drift, latitude, free-air, simple Bouguer, curvature, terrain, and isostatic corrections (e.g., Blakely, 1995).

Seismic Data

We acquired a high-resolution seismic reflection profile located over the narrowest portion of the magnetic high (Fig. 2) where the magnetic source is expected to be shallowest based on a simple rule-of-thumb depth-to-source calculation (Peter’s method; e.g., Blakely, 1995). The reflection profile was shot using a Betsy Seisgun and collected with a 955 m cabled recording array with 8 linked, 24-channel Geometrics Geode Ultra-Light Exploration ( seismographs. Details of the seismometer array, shot spacing, and recording parameters were provided in Fontiveros (2010). The ideal conditions of fine-grained lake deposits saturated by water almost immediately below the surface enabled penetration depths to 400 m (Fontiveros, 2010).

The acquired seismic data were processed using ProMAX seismic software. To address the variable frequency response of trace sections in the shot gathers (due to the two different types of geophones used), a 60–120 Hz bandpass filter was applied. Reverberations were suppressed using a predictive-deconvolution operator (60 ms operator length, 15 ms prediction distance). Amplitudes at depth were enhanced by applying a 175 ms automatic gain control. Airwave frequencies in the data set were similar to significant reflections (∼60–80 Hz); therefore, instead of filtering out the airwave, we applied a bottom mute (all samples recorded after the onset of the ground-roll were set to zero on each trace) that removed both the ground roll and the airwave. In addition, refractions at the top of the record were muted so only reflection energy was processed and stacked. A detailed velocity analysis was performed to characterize the complex lateral velocity variations, and a normal moveout correction was applied. The data were migrated using Kirchhoff pre-stack depth migration (Fontiveros, 2010). Our interpretation utilized both the migrated (Fig. 3A) and the unmigrated (see Supplemental Fig. 11) sections.

The magnetic anomaly map reveals an isolated, ∼35-km-long, nearly continuous magnetic high (Fig. 1B). The majority of the high is <500 m wide, although the northern and southern extents broaden to a width of 1–2 km. Profiles across the high show that its amplitude, wavelength, and shape are highly variable despite its continuity (Fig. 2A). In the northern profiles (lines 5–8) where the high is broadest, the shape is asymmetric with a gentle gradient west of the peak and a steeper gradient to the east. In the southern profiles (lines 11–19) the sense of asymmetry is opposite. In lines 9 and 10, the peaks of the high are 50 and 100 m east of the trend from other profiles, and several profiles have multiple peaks (lines 6, 7, 8, 9, and 17). Given its isolation from other anomalies (Fig. 1B), the variability within the anomaly is likely a primary feature of the causative body rather than due to interactions with other magnetic sources.

Gravity profiles reveal an ∼0.5 mGal low within a long-wavelength gradient (Fig. 2B). The trace of the gravity low is colocated with the peak of the magnetic high for ∼6 km, but is east of the magnetic high where it broadens to the north (Fig. 2C).

In the seismic reflection profile, the highest amplitude reflection (B in Fig. 3A) appears between 200 and 400 m depth, dips westward, and is laterally discontinuous, offset ∼340 m (dip-slip down to the east) in the middle of the profile. Both west and east of the offset, reflection B has a synformal shape that likely represents “smiles,” artifacts that result from migrating seismic data that are imperfect (e.g., Warner, 1987). In our data the limited line length did not capture the full amplitude response of reflection B at either end of the profile, and attenuation prevented imaging the complete diffractions where B is offset (Supplemental Fig. 1). We therefore interpret reflection B as homoclinal (Fig. 3A) as the simplest interpretation consistent with the data. Above reflection B to ∼100 m depth is a zone of low-amplitude west-dipping reflections (A in Fig. 3). No reliable reflections are imaged below reflection B, likely due to seismic ringing on a thick basalt layer.

Modeling Along Seismic Profile

In order to differentiate between possible sources of the magnetic high, we developed two-dimensional potential field forward models along the seismic profile using a commercial two-dimensional forward modeling package (GM-SYS; Because modeling solutions of potential field anomalies are nonunique, the highest amplitude reflection in the seismic data (B in Fig. 3) served as a key constraint. Physical property data (density, magnetic susceptibility, and magnetization measurements) collected from surface samples also provided important constraints for the modeling process. The values employed in the models (Tables 1 and 2) are primarily based on published physical property measurements (Ponce et al., 2009) and paleomagnetic measurements (Ritzinger, 2014) of units exposed at the surface in closest proximity to the model. In some cases, however, modeled layers cannot be sampled, or the use of known physical properties was not sufficient to match the observations, so in these cases, physical property values are estimated (these values are described where relevant).

Most exposed basalt flows in the area are reversely magnetized, although normal and transitional flows are also present (Ritzinger, 2014). A 292-m-deep geothermal gradient hole located 5 km south of the seismic profile along the magnetic high (Fig. 1) recovered reversely magnetized basalt core from 88 to 89 m depth and 143 to 145 m depth (J. Glen, personal data), while the rest of the core comprised sediments. To account for sediments interbedded in the basalts and for rubbly vesiculated layers between individual flows, physical property values for the basalt layer (Tbl) were calculated assuming that 50% of the layer had the properties of alluvium (Qal), which is less dense and less magnetic (Table 1).

Given the region’s history of extension and associated volcanism, there are limited geologically consistent possibilities for the structures and associated features likely to appear in the subsurface. For that reason, we considered three possible end-member models for the source of the magnetic high: (1) the normal fault model, which relies on faulted magnetic stratigraphy, (2) the ponded basalt model, where basalt flows fill faulted topography, and (3) the dike model, which relies on a mafic intrusion. End-member models are useful to determine the primary influences on the source of the anomaly, even if particular models are known to be incorrect, incomplete, or problematic. As described in the following section, the first step in the modeling process was to determine the characteristic magnetic and gravity fields produced by these end-member structures (Fig. 4). While these models are simplistic, they show that the end-member sources produce very distinct anomalies.

Results of End-Member Models

Our normal fault model (Fig. 4A) depicts a west-dipping magnetic layer that is cut and offset by an east-dipping fault based on offset reflections in the seismic profile that we interpret to be the top of the late Miocene–Pliocene basalts (Tbl). The transition from Quaternary lacustrine deposits to basalt provides the acoustic contrast necessary to produce the high-magnitude reflection. The dip of the Tbl layer produces a gentle gravity gradient that reflects the shape of the basin (Fig. 4A). Within the gravity gradient, the offset of the Tbl layer produces a gravity low where less dense alluvium (Qal) fills in above the hanging wall. The depth and orientation of Oligocene volcanics (Tv) were estimated by projecting mapped units (from Egger and Miller, 2011) into the subsurface, although this layer has relatively little influence in the model. This model correctly produces the observed gravity profile; therefore, an interpretation of a fault from the seismic data is supported by the gravity data. The fit of the magnetic data, however, is poor. In the first calculation (red line in Fig. 4A), magnetic parameters used for the Tbl layer (Table 2) are based on our best estimate from magnetic susceptibility and remanence measurements (Ritzinger, 2014), taking into account less magnetic sediment that is interbedded in the basalt. In the second calculation (blue dashed line in Fig. 4A), maximum magnetic parameters are used (Table 2) based on the upper end of measured magnetic values (Ritzinger, 2014). In both cases, however, the calculated magnetic high is not of sufficient amplitude.

In our ponded basalt model (Fig. 4B), we added a highly magnetic layer that conceptually represents a basalt flow (Tpb) that pooled on the hanging wall of the faulted topography (Fig. 4B); otherwise, this model is equivalent to the normal fault model using the best estimate for the Tbl layer’s magnetic parameters. In the first calculation (red line in Fig. 4B), the resulting magnetic profile is virtually unchanged from the normal fault model (Fig. 4A), due to the Tpb layer’s planar horizontal orientation and magnetic parameters that are similar to those of the Tbl layer (Table 2). In the second calculation (blue dashed line in Fig. 4B), changing the Tpb layer from normal to reversed polarity has the effect of increasing the amplitude of the magnetic high, which is closer to fitting the observed high (Fig. 4B). This model, however, is inconsistent with the seismic profile. If a pooled basalt flow were present in the thickness indicated in the model, we would expect it to appear in the seismic reflection profile.

Our dike model (Fig. 4C) shows a 10-m-thick dike in two orientations, vertical and east dipping, intruding horizontal stratigraphy. Although this does not conform to the dipping, faulted reflections in the seismic data, our intent was to assess the gravity and magnetic contributions of the dike alone rather than the faulted stratigraphy that is already depicted in the first two models. The dike is assumed to be highly magnetic and dense, similar in rock properties to a subsurface flow modeled by Egger et al. (2010). For convenience we have modeled the dike as a single, rectangular block, while recognizing that dikes generally thin toward the tip (e.g., Gudmundsson, 2003) and that at 10 m thickness, the dike may be a composite structure representing multiple successive intrusions. The depth to the top of the dike was chosen in order to minimally contribute to the gravity field (i.e., there is no observed gravity high). In the first calculation (red line in Fig. 4C), the magnetic high produced by the vertical dike does not match the sense of asymmetry of the observed magnetic high, which has a steeper gradient on the western side of the high. In the second calculation (blue dashed line in Fig. 4C), the sense of asymmetry matches the observed high more closely.

Final Model

As expected, none of the end-member models reproduce the observed high-amplitude magnetic high aligned with a gravity low. However, by combining the normal fault model (Fig. 4A) and the dipping dike model (Fig. 4C) into a single model, both the asymmetric high-amplitude magnetic high and gravity low can be reproduced (Fig. 5). The gravity fit was further improved by the addition of two alluvium layers (Qal) that increase the density of alluvium with depth, consistent with typical basin sediment depth profiles (Brocher, 2008). A low-density tuffaceous sediment layer (Tc) was also added to the model, consistent with tuffs modeled (by Egger et al., 2010) along the nearby seismic line acquired by Lerch et al. (2010) (Fig. 1A).

Nature and Development of the Magnetic Anomaly

Colocation of a magnetic high with a gravity low is unusual because of the common association of high magnetizations with high-density mafic bodies. Therefore, no end-member potential field model (Fig. 4) accurately reproduces the anomalies in Surprise Valley. Even if a single causative body were able to produce the magnetic high and gravity low along the seismic profile, the broadening of the magnetic high to the north and its divergence from the gravity low (Figs. 1 and 2) would preclude such a model. Instead, only a model that combines a normal fault with a dike fits the observed gravity, magnetic, seismic, and geologic data (Fig. 5).

The ∼10 m thickness and ∼150 m depth of the dike, as well as its relationship to the fault, are only partially constrained by the modeling. In order to fit the observed gravity low, the sole constraint on the dike is that it must be sufficiently thin and deep that it does not produce a gravity anomaly (Fig. 4C). However, because the proportion of thickness to depth may vary, and because the dike cannot be sampled for physical properties, the dike’s thickness and depth may vary by tens of meters without affecting the gravity fit. Nevertheless, the ∼10 m thickness is consistent with displacement length scaling relations that predict a maximum opening of 13 m for a 35-km-long dike (Schultz et al., 2008, fig. 5 therein). Furthermore, the location of the top of the dike at the contact between the basalt (Tbl) and sediment (Qal) is reasonable given field observations and modeling of dike arrest in layered crust (Gudmundsson, 2002).

Given the relatively shallow dip of the dike and its interaction with the fault, a key question is whether this scenario is likely, based on the processes of diking and faulting. Valentine and Krogh (2006) observed fault-captured dikes in Paiute Ridge, Nevada, finding that the dikes, 400–5000 m long and 1.2–9 m thick, only occupied normal faults that were steeper than 60°. Their field observations are supported by analytic and numerical modeling (Gaffney et al., 2007), which find that in addition to steep fault angles, fault capture of propagating dikes becomes more favorable at shallow depths and with high fracture toughness in the hanging wall. Assuming that the basalts (Tbl) have a fracture toughness of 1 MPa m1/2, which is likely an underestimate (Gaffney et al., 2007, Table 1 therein), and the fault dip angle is 60°, fault capture of a rising dike is permitted at depths <300–800 m, depending on preexisting cracks in the hanging wall. Although the fault in our model dips 56° to the east (shallower than expected for dike capture), stratigraphy dips 15° to the west, indicating that the fault probably formed at an angle as steep as 71° and, as is typical of normal faults in the Basin and Range, rotated to lower dips as motion occurred (Chamberlin, 1983). The dike could have been captured at any point during rotation of the fault plane.

We also note that the dike is located ∼8–10 km east of the main trace of the Surprise Valley fault, and is roughly parallel to the fault along its length. This places the dike and magnetic anomaly directly above the location of a predicted shallowing of the brittle-ductile transition (inset, Fig. 1A) (Lerch et al., 2010). A rise of this transition zone by as much as 3 km may have occurred during a period of rapid slip on the Surprise Valley fault after 3 Ma (Colgan et al., 2008), a suggestion supported by flexure observed in the Warner Range (Egger and Miller, 2011) and high heat flow in the basin (Blackwell et al., 1991; Benoit et al., 2005). The location of the dike is consistent with where decompression melting would most likely occur.

On the basis of these geological constraints, we have developed a conceptual model for the sequence of events leading to the observed phenomena (Fig. 6). In our schematic diagram, a steeply dipping normal fault forms in the crust in response to basin extension ca. 3 Ma; motion along the fault results in minor offset of late Miocene–Pliocene basalts (Tbl), indicating that faulting initiated after 3.8 Ma (Fig. 6A). Around the same time, rotation of the Surprise Valley fault and rapid uplift of the Warner Range raised the brittle-ductile transition zone (Lerch et al., 2010), generating the magma supply for dike intrusion. A dike rose subvertically through the crust and encountered the mechanically strong Miocene–Pliocene basalts in the shallow subsurface; as the dike moved through this strong layer, it was captured and diverted by the normal fault (Fig. 6B). Upon rising to the contact between the basalt and soft sediment just below the paleosurface, the dike was arrested (Fig. 6B). The age of the dike is not well known. An older age, during the post–3 Ma episode of rapid offset of the Surprise Valley fault, implies that the dike would have encountered the fault prior to significant tilt, and would have been likely to be captured by the fault. A younger age would have allowed greater accumulation of soft sediment above the fault tip, increasing the probability of dike arrest in the subsurface rather than eruption onto the paleolakebed (although examples are known where dikes were arrested within 5 m of the surface; Gudmundsson, 2003). Continued motion along the fault (Fig. 6C) and ongoing sedimentation resulted in tilting of the dike within the fault block and fanning of sedimentary deposits (Fig. 6D).

Extension in the Basin

Utilizing the results from our modeling, the extent of the seismically imaged and modeled fault can be mapped using the small gravity low as a proxy for its location (Fig. 2C). By connecting the location of the gravity low along several transects, we interpret the buried fault to have an orientation similar to east-dipping faults located to the east in the Larkspur Hills (Fig. 1), which range in strike from 350° to 010° (Strickley, 2014). In Egger and Miller (2011), ∼7.3 km of east-west extension across 50 km (or 15% extension) was calculated based on geologic mapping and offset calculations along exposed faults, including the Surprise Valley fault and several regional faults. Strickley (2014) calculated extension along several profiles across the Larkspur Hills, finding ∼460 m of east-west extension across 10 km (or 5% extension) at the latitude of our model. In comparison, the fault in our model accommodates 200 m of horizontal slip (Fig. 3B). While this alone is not a significant amount of extension, it is possible that there are additional, as yet unidentified, intrabasin faults whose combined slip may accommodate a sizable amount of extension.

Geophysical work in other valleys throughout the Basin and Range suggest that multiple intrabasin faults are likely (e.g., Okaya and Thompson, 1985). Our densely spaced gravity data, which targeted the magnetic high, did not extend far enough into the basin to reveal additional intrabasin faults, and existing regional gravity coverage (Ponce et al., 2009) on a 1.6 km grid is not sufficient to reveal faults that cause <1 mGal anomalies. Furthermore, deeper, basinward faults would produce a significantly smaller gravity signal that may preclude their detection by gravity alone. One way to estimate the number and size of faults that may exist in the basin is to assume that the fault population follows a fractal size distribution (Marrett and Allmendinger, 1992). Based on the average 2 km lateral spacing of faults in the Larkspur Hills (Strickley, 2014) and the intrabasin fault identified in the reflection data, we estimate that there could be 5–10 more faults in the basin resulting in 1–2 km of additional extension hidden beneath the playa, and regional extension of 16%–19%. As a result, extension calculations in the Surprise Valley region likely underestimate the total extension that has occurred. These results have broader implications for similar extensional basins that remain concealed under basin sediments and have not been characterized geophysically.

Implications Along Strike of the Anomaly

The results from our modeling also provide a basis for interpreting complex features in the magnetic data throughout the basin. Between the middle and upper basins, the magnetic high crosses a region previously referred to as the Lake City fault zone or Lake City fault (Fig. 1A) (Hedel, 1980; U.S. Geological Survey and California Geological Survey, 2006) that was widely cited as an important throughgoing structural element controlling geothermal circulation in the valley. However, it was concluded (Egger et al., 2014) that the region lacks a throughgoing fault on the basis that there is no consistent gravity, magnetic, or resistivity signature coinciding with the fault, and therefore should not be mapped as such. Our results support this interpretation; we see clear continuity of the magnetic high across the zone with no evidence of offset within this region or elsewhere along the length of the anomaly. There is, however, evidence of structural complexity in the shallow subsurface where the magnetic high is more diffuse, consistent with the presence of a deformation zone, such as may develop where faults interact, as suggested in Egger et al. (2014).

The broadening of the magnetic anomaly at its northern and southern extent (Fig. 1B) may be the result of the dike rising into more complex structure in the shallow subsurface (e.g., Keating et al., 2008). In the upper basin, the magnetic high not only broadens in the northern profiles, but also subtly changes in the sense of asymmetry and, in some profiles, has multiple peaks (Fig. 2A), suggesting that the dike is no longer controlled by a single east-dipping fault. Preliminary modeling suggests that the dike may be rising into a small horst or broader fault zone (Athens, 2011), although the lack of seismic reflection or well-log data means that the modeling is poorly constrained.

The acquisition of multiple, complementary geophysical data sets provided insight into subsurface features in Surprise Valley as it allowed us to identify features that are common across data sets or unique to one data set. Our combined analysis of geophysical data and modeling identifies a fault-controlled ∼10-m-thick dike at a depth of ∼150 m. The location of the dike, ∼8–10 km east of the main trace of the Surprise Valley fault, corresponds to an area of a predicted elevated brittle-ductile transition zone (Lerch et al., 2010), which is precisely where decompression melting and diking would occur. The concealed fault, imaged by reflection data, is located 3 km west of the eastern edge of the basin and accounts for an additional ∼200 m of extension. A further 1–2 km extension along unsampled intrabasin faults is considered possible based on the expected fault population.

Our results indicate that structures and magmatic features in the subsurface are important contributors to a complete assessment of the processes, timing, and total amount of extension observed in continental environments. Only through combined geological and geophysical analyses can we get a complete picture of the extensional history of regions such as the northwestern Basin and Range.

We are indebted to all the volunteers who helped collect data, in particular Claire Bouligand, Bruce Chuchel, Kevin Denton, Nate Levine, Erin Looby, Nellie Olsen, and Alex Yu. We also thank Rick Blakely, Agust Gudmundsson, Francesco Mazzarini, and one anonymous reviewer for providing comments that greatly improved the manuscript. Magnetic and gravity acquisition was funded by the U.S. Geological Survey and Stanford University. Seismic data acquisition was funded by Stanford University, using recording equipment from the IRIS-PASSCAL (Incorporated Research Institutions for Seismology–Portable Array Seismic Studies of the Continental Lithosphere) equipment pool. Landmark Geophysical generously donated the ProMAX software license used for reflection processing.

1 Supplemental Figure 1. Unmigrated common midpoint stack. The high-amplitude reflection (B in Fig. 3A) is clearly homoclinal rather than synformal. Depth conversion is based on the velocity model of Figure 3A. Please visit or the full-text article on to view Supplemental Figure 1.