William D. Stanley

Surface seismic and electrical methods to detect fluids related to faulting

In the absence of drilling, surface-based geophysical methods are necessary to observe fault zones and fault zone physical properties at seismogenic depths. These in situ physical properties can then be used to infer the presence and distribution of fluids along faults, although such observations are by nature indirect and become less exact with greater depth. Multiple observations of a range of such geophysical properties as compressional and shear seismic velocity (Vp and Vs), Vp/V5 ratio (related to Poissons ratio), resistivity and attenuation in and adjacent to fault zones offer the greatest hope of making inferences of the fault zone geometry, fluids in the fault zone, and fluid reservoirs in the surrounding crust. For simple geometries, fault zone guided waves can provide information on fault zone width and velocities for faults of the order of 200 m wide. To address the question of whether a narrow fault zone can be imaged well enough at depths of seismic rupture to infer the presence of anomalously high fluid/rock ratios, we present synthetic seismic tomography and magnetotelluric examples for an ideal case of a narrow fault zone with a simple geometry, large changes in material properties, and numerous earthquakes within the fault zone. A synthetic 0.5-km wide fault zone with 20% velocity reduction is well imaged using local earthquake tomography. When sequential velocity inversions are done, the true fault width is found, even to 9 km depth, although the calculated amplitude of the velocity reduction is lower than the actual amplitude. Vp/Vs is as well determined as Vp. Magnetotelluric imaging of a synthetic fault zone shows that a conductive fault zone can be well imaged within the upper 10 km. Further, a narrow (1 km) very low resistivity (3 ohm m) fault core can be imaged within a broad (5 km) low resistivity (10 ohm m) fault zone, illustrating that regions of a fault containing large quantities of interconnected fluids within a broader, conductive fault zone should be detectable. Thus variations in fluid content and fluid pressure can be inferred from electrical and seismic methods but there will always be uncertainty in these inferences due to the trade-off with other factors, such as intrinsic variations in porosity, mineralogy, and pore geometry. The best approach is combined modeling of varied seismic and electrical data.

Active tectonics of the Seattle fault and central Puget sound, Washington - Implications for earthquake hazards

We use an extensive network of marine high-resolution and conventional industry seismic-reflection data to constrain the location, shallow structure, and displacement rates of the Seattle fault zone and crosscutting high-angle faults in the Puget Lowland of western Washington. Analysis of seismic profiles extending 50 km across the Puget Lowland from Lake Washington to Hood Canal indicates that the west-trending Seattle fault comprises a broad (4–6 km) zone of three or more south-dipping reverse faults. Quaternary sediment has been folded and faulted along all faults in the zone but is clearly most pronounced along fault A, the northernmost fault, which forms the boundary between the Seattle uplift and Seattle basin. Analysis of growth strata deposited across fault A indicate minimum Quaternary slip rates of about 0.6 mm/yr. Slip rates across the entire zone are estimated to be 0.7–1.1 mm/yr. The Seattle fault is cut into two main segments by an active, north-trending, high-angle, strike-slip fault zone with cumulative dextral displacement of about 2.4 km. Faults in this zone truncate and warp reflections in Tertiary and Quaternary strata and locally coincide with bathymetric lineaments. Cumulative slip rates on these faults may exceed 0.2 mm/yr. Assuming no other crosscutting faults, this north-trending fault zone divides the Seattle fault into 30–40-km-long western and eastern segments. Although this geometry could limit the area ruptured in some Seattle fault earthquakes, a large event ca. A.D. 900 appears to have involved both segments. Regional seismic-hazard assessments must (1) incorporate new information on fault length, geometry, and displacement rates on the Seattle fault, and (2) consider the hazard presented by the previously unrecognized, north-trending fault zone.

The Denali fault system and Alaska Range of Alaska: Evidence for underplated Mesozoic flysch from magnetotelluric surveys

Regional magnetotelluric surveys recently completed across the central and eastern Alaska Range of Alaska provide evidence for large volumes of conductive rocks beneath the core of the range. These conductive rocks may represent a formerly extensive, but now collapsed, Mesozoic flysch basin formed on the leading edge of the Talkeetna superterrane (amalgamated Wrangellia, Peninsular, and Alexander terranes). The docking of the Talkeetna superterrane caused large-scale oblique thrusting, folding, and metamorphism in the flysch basin, and formation of a megasuture along which the Cenozoic strike-slip Denali fault system developed. The deep magnetotelluric soundings and seismic reflection data suggest the possibility that the highly conductive rocks were tectonically emplaced beneath the thin crystalline sheet constituting the southern Yukon-Tanana terrane over a broad region of the Alaska Range. The conductive rocks are locally correlated with surface outcrops of Mesozoic black shales that are part of Upper Jurassic and Cretaceous flysch but may be composed of Paleozoic carbonaceous shales as well. In either case, their extremely low resistivities make them a valuable marker horizon for tectonic studies. The conductive rocks are interpreted to extend to depths of greater than 20 km and were mapped north and northeast of the Denali fault for more than 50 km. The magnetotelluric surveys represent the first large-scale surveys done in Alaska, but the structures mapped are similar to those observed in large, compressed flysch basins in the eastern Alps and Carpathian Mountains of Europe. The results of these surveys bear on several key tectonic questions, including development of the ancestral Denali fault, and collapse and possible underplating of an extensive Mesozoic flysch system and associated igneous arc.

Comparison of geoelectrical/tectonic models for suture zones in the western U.S.A. and eastern Europe: are black shales a possible source of high conductivities?

Abstract Large-scale geoelectrical anomalies have been mapped with geomagnetic depth sounding (GDS) and magnetotelluric (MT) surveys in the Carpathian Mountains region. These anomalies are associated with the zone of closure between stable Europe and a complex of microplates in front of the converging African plate. The zone of closure, or suture zone, is largely occupied by an extensive deformed flysch belt. The models derived to fit the observed geoelectrical data are useful in the study of other suture zones, and Carpathian structures have been compared with areas currently being studied in the western Cordillera of the U.S.A. Models derived for a smaller-scale suture zone mapped in western Washington State have features that are similar to the Carpathian models. The geoelectrical models for both the Carpathian and Washington anomalies require dipping conductive slabs of 1–5 Ω m material that extends to depths > 20 km. In both instances there is evidence that these materials may merge with lower crustal-mantle conductors along the down-dip margins of the slab. The main conductive units are interpreted to be sedimentary rocks that have been partially subducted due to collisional processes. Heat flow is low in both regions and it is difficult to explain fully the deep conduction mechanisms; however, evidence suggests that the conduction at depth may include electronic conduction in sulfide mineral or carbon films as well as ionic conduction in fluids or partial melt.

Tectonic controls on magmatism in The Geysers–Clear Lake region: Evidence from new geophysical models

In order to study magmatism and geothermal systems in The Geysers–Clear Lake region, we developed a detailed three-dimensional tomographic velocity model based on local earthquakes. This high-resolution model resolves the velocity structure of the crust in the region to depths of approximately 12 km. The most significant velocity contrasts in The Geysers–Clear Lake region occur in the steam production area, where high velocities are associated with a Quaternary granitic pluton, and in the Mount Hannah region, where low velocities occur in a 5-km-thick section of Mesozoic argillites. In addition, a more regional tomographic model was developed using traveltimes from earthquakes covering most of northern California. This regional model sampled the whole crust, but at a lower resolution than the local model. The regional model outlines low velocities at depths of 8–12 km in The Geysers–Clear Lake area, which extend eastward to the Coast Range thrust. These low velocities are inferred to be related to unmetamorphosed Mesozoic sedimentary rocks. In addition, the regional velocity model indicates high velocities in the lower crust beneath the Clear Lake volcanic field, which we interpret to be associated with mafic underplating. No large silicic magma chamber is noted in either the local or regional tomographic models. A three-dimensional gravity model also has been developed in the area of the tomographic imaging. Our gravity model demonstrates that all density contrasts can be accounted for in the upper 5–7 km of the crust. Two-dimensional magnetotelluric models of data from a regional, east-west profile indicate high resistivities associated with the granitic pluton in The Geysers production area and low resistivities in the low-velocity section of Mesozoic argillites near Mount Hannah. No indication of midcrustal magma bodies is present in the magnetotelluric data. On the basis of heat flow and geologic evidence, Holocene intrusive activity is thought to have occurred near the Northwest Geysers, Mount Hannah, Sulphur Bank Mine, and perhaps other areas. The geophysical data provide no conclusive evidence for such activity, but the detailed velocity model is suggestive of intrusive activity near Mount Hannah similar to that in the “felsite” of The Geysers production area. The geophysical models, seismicity patterns, distribution of volcanic vents, heat flow, and other data indicate that small, young intrusive bodies that were injected along a northeast trend from The Geysers to Clear Lake probably control the thermal regime.

The Geysers-Clear Lake geothermal area, California - an updated geophysical perspective of heat sources

Abstract The Geysers-Clear Lake geothermal area encompasses a large dry-steam production area in The Geysers field and a documented high-temperature, high-pressure, water-dominated system in the area largely south of Clear Lake, which has not been developed. Both systems have been extensively studied with geophysical techniques, drilling, and geological mapping during the past 20 years. An updated view is presented of the geological/geophysical complexities of the crust in The Geysers-Clear Lake region in order to address key unanswered questions about the heat source and tectonics. Early geophysical interpretations used a gravity low centered in the area between Clear Lake and The Geysers to suggest that a large magma chamber existed at depths starting at about 7 km. This first-order assumption of a large magma chamber expressed in the gravity data was used as a guide in subsequent geophysical and geological interpretations. Drill-hole temperature evidence is strongly suggestive of a shallow, hot-intrusive body, but in this paper the complexities are documented of the geological and geophysical data sets that make it difficult to pinpoint the location of “magma” or hot, solidified intrusive material. Forward modeling, multidimensional inversions, and ideal body analysis of the gravity data, new electromagnetic sounding models, and arguments made from other geophysical data sets suggest that many of the geophysical anomalies have significant contributions from rock property and physical state variations in the upper 7 km and not from ”magma“ at greater depths. Regional tectonic and magmatic processes are analyzed to develop an updated scenario for pluton emplacement that differs substantially from earlier interpretations. In addition, a rationale is outlined for future exploration for geothermal resources in The Geysers-Clear Lake area.

Computer graphics in geophysics

It has become increasingly apparent that the large quantities of geophysical data being collected by modern instrumentation must be efficiently stored, processed, and displayed before meaningful interpretation can be undertaken. To meet these needs a man‐machine interactive computing system has been developed to provide the geophysicist with a direct communication link to the digital computer and methods for providing two‐ and three‐dimensional representations of geophysical data on cathode‐ray tubes (CRT). Gravity, magnetic, resistivity, and induced‐polarization responses have been modeled on the interactive system and displayed on a CRT. The interactive modeling processes allow the direct interchange of geophysical variables and data followed by the display of the recomputed models and corresponding theoretical responses. Model parameters can be monitored visually and continuously adjusted for a better solution. In addition, color‐shaded perspective views of three‐dimensional surfaces representing aerom...

Geophysical Study of Unconsolidated Sediments and Basin Structure in Cache Valley, Utah and Idaho

During the summer of 1969, an integrated geophysical study was completed in Cache Valley near the state line between Utah and Idaho. The survey covered about 30 sq mi and was designed to provide data on unconsolidated sediments and basin structure in a study of the ground-water resources of the area. Gravity, magnetic, resistivity, and seismic reflection methods were employed to obtain the desired information. Detailed gravity data revealed that the west side of Cache Valley is underlain by a trough which is about 5 mi wide and which is filled with as much as 7,000 ft of Quaternary and Tertiary sedimentary units. The magnetic data revealed a probable mafic dike intruded into Paleozoic and Precambrian bedrock in the Lewiston area and possibly a buried volcanic flow unit in the Tertiary sediments. Interpretations of 38 vertical electrical soundings supported by three well logs showed that an area of about 8 sq mi in the center of the surveyed area is underlain by conductive clays (resistivity less than 10 ohm-meters) that are as much as 400 ft thick. Most of the area surveyed is underlain at depths of 100 to 500 ft by more resistive materials (20 to 40 ohmmeters) interpreted to contain a higher percentage of coarse-grained sediments. These sediments may constitute a potential aquifer. Seismic reflection data further confirmed the existence of the massive clays in the central part of the surveyed area, revealed complex sedimentary geology at the valley margins, and defined subtle structures within the valley fill.

Reply by Robert B. Smith and W. D. Stanley to discussion by Nelson C. Steenland

The point of Mr. Steenland’s discussion appears that he does not consider that the computer‐generated perspective views of geophysical maps are as clear or distinct as contour maps. Further he asserts that they lack a quantitative use. On the contrary the perspective views were produced from the same data as the contoured map but in a quantitative digital form. We believe that perspective views are by no means a substitute for contour maps but provide a complementary form of three‐dimensional data.