M. R. Goldman

Relationships among seismic velocity, metamorphism, and seismic and aseismic fault slip in the Salton Sea Geothermal Field region

The Salton Sea Geothermal Field is one of the most geothermally and seismically active areas in California and presents an opportunity to study the effect of high-temperature metamorphism on the properties of seismogenic faults. The area includes numerous active tectonic faults that have recently been imaged with active source seismic reflection and refraction. We utilize the active source surveys, along with the abundant microseismicity data from a dense borehole seismic network, to image the 3-D variations in seismic velocity in the upper 5 km of the crust. There are strong velocity variations, up to ~30%, that correlate spatially with the distribution of shallow heat flow patterns. The combination of hydrothermal circulation and high-temperature contact metamorphism has significantly altered the shallow sandstone sedimentary layers within the geothermal field to denser, more feldspathic, rock with higher P wave velocity, as is seen in the numerous exploration wells within the field. This alteration appears to have a first-order effect on the frictional stability of shallow faults. In 2005, a large earthquake swarm and deformation event occurred. Analysis of interferometric synthetic aperture radar data and earthquake relocations indicates that the shallow aseismic fault creep that occurred in 2005 was localized on the Kalin fault system that lies just outside the region of high-temperature metamorphism. In contrast, the earthquake swarm, which includes all of the M > 4 earthquakes to have occurred within the Salton Sea Geothermal Field in the last 15 years, ruptured the Main Central Fault (MCF) system that is localized in the heart of the geothermal anomaly. The background microseismicity induced by the geothermal operations is also concentrated in the high-temperature regions in the vicinity of operational wells. However, while this microseismicity occurs over a few kilometer scale region, much of it is clustered in earthquake swarms that last from hours to a few days and are localized near the MCF system.

Structure of the Koyna‐Warna Seismic Zone, Maharashtra, India: A possible model for large induced earthquakes elsewhere

The Koyna-Warna area of India is one of the best worldwide examples of reservoir-induced seismicity, with the distinction of having generated the largest known induced earthquake (M6.3 on 10 December 1967) and persistent moderate-magnitude (>M5) events for nearly 50 years. Yet, the fault structure and tectonic setting that has accommodated the induced seismicity is poorly known, in part because the seismic events occur beneath a thick sequence of basalt layers. On the basis of the alignment of earthquake epicenters over an ~50 year period, lateral variations in focal mechanisms, upper-crustal tomographic velocity images, geophysical data (aeromagnetic, gravity, and magnetotelluric), geomorphic data, and correlation with similar structures elsewhere, we suggest that the Koyna-Warna area lies within a right step between northwest trending, right-lateral faults. The sub-basalt basement may form a local structural depression (pull-apart basin) caused by extension within the step-over zone between the right-lateral faults. Our postulated model accounts for the observed pattern of normal faulting in a region that is dominated by north-south directed compression. The right-lateral faults extend well beyond the immediate Koyna-Warna area, possibly suggesting a more extensive zone of seismic hazards for the central India area. Induced seismic events have been observed many places worldwide, but relatively large-magnitude induced events are less common because critically stressed, preexisting structures are a necessary component. We suggest that releasing bends and fault step-overs like those we postulate for the Koyna-Warna area may serve as an ideal tectonic environment for generating moderate- to large- magnitude induced (reservoir, injection, etc.) earthquakes.

Detailed P- and S-Wave Velocity Models along the LARSE II Transect, Southern California

Abstract Structural details of the crust determined from P -wave velocity models can be improved with S -wave velocity models, and S -wave velocities are needed for model-based predictions of strong ground motion in southern California. We picked P - and S -wave travel times for refracted phases from explosive-source shots of the Los Angeles Region Seismic Experiment, Phase II (LARSE II); we developed refraction velocity models from these picks using two different inversion algorithms. For each inversion technique, we calculated ratios of P - to S -wave velocities ( V P / V S ) where there is coincident P - and S -wave ray coverage. We compare the two V P inverse velocity models to each other and to results from forward modeling, and we compare the V S inverse models. The V S and V P / V S models differ in structural details from the V P models. In particular, dipping, tabular zones of low V S , or high V P / V S , appear to define two fault zones in the central Transverse Ranges that could be parts of a positive flower structure to the San Andreas fault. These two zones are marginally resolved, but their presence in two independent models lends them some credibility. A plot of V S versus V P differs from recently published plots that are based on direct laboratory or down-hole sonic measurements. The difference in plots is most prominent in the range of V P =3 to 5 km/s (or V S ∼1.25 to 2.9 km/s), where our refraction V S is lower by a few tenths of a kilometer per second from V S based on direct measurements. Our new V S - V P curve may be useful for modeling the lower limit of V S from a V P model in calculating strong motions from scenario earthquakes.

Subsurface Geometry of the San Andreas Fault in Southern California: Results from the Salton Seismic Imaging Project (SSIP) and Strong Ground Motion Expectations

The San Andreas fault (SAF) is one of the most studied strike‐slip faults in the world; yet its subsurface geometry is still uncertain in most locations. The Salton Seismic Imaging Project (SSIP) was undertaken to image the structure surrounding the SAF and also its subsurface geometry. We present SSIP studies at two locations in the Coachella Valley of the northern Salton trough. On our line 4, a fault‐crossing profile just north of the Salton Sea, sedimentary basin depth reaches 4 km southwest of the SAF. On our line 6, a fault‐crossing profile at the north end of the Coachella Valley, sedimentary basin depth is ∼2–3  km and centered on the central, most active trace of the SAF. Subsurface geometry of the SAF and nearby faults along these two lines is determined using a new method of seismic‐reflection imaging, combined with potential‐field studies and earthquakes. Below a 6–9 km depth range, the SAF dips ∼50°–60° NE, and above this depth range it dips more steeply. Nearby faults are also imaged in the upper 10 km, many of which dip steeply and project to mapped surface fault traces. These secondary faults may join the SAF at depths below about 10 km to form a flower‐like structure. In Appendix D, we show that rupture on a northeast‐dipping SAF, using a single plane that approximates the two dips seen in our study, produces shaking that differs from shaking calculated for the Great California ShakeOut, for which the southern SAF was modeled as vertical in most places: shorter‐period (T<1  s) shaking is increased locally by up to a factor of 2 on the hanging wall and is decreased locally by up to a factor of 2 on the footwall, compared to shaking calculated for a vertical fault.

VS30 at three strong-motion recording stations in Napa and Napa County, California—Main Street in downtown Napa, Napa fire station number 3, and Kreuzer Lane—Calculations determined from s-wave refraction tomography and multichannel analysis of surface waves (Rayleigh and Love)

The August 24, 2014, moment magnitude (Mw) 6.0 South Napa earthquake caused an estimated

Faulting Apparently Related to the 1994 Northridge, California, Earthquake and Possible Co-seismic Origin of Surface Cracks in Potrero Canyon, Los Angeles County, California

400 million in structural damage to the City of Napa, California. In 2015, we acquired high-resolution Pand S-wave seismic data near three strong-motion recording stations in Napa County where high peak ground accelerations (PGAs) were recorded during the South Napa earthquake. In this report, we present results from three sites—Main Street in Downtown Napa (Northern California Seismic Network station, NCSN N016), Napa Fire Station Number 3 (National Strong Motion Project station, NSMP 1765), and Kreuzer Lane (station KRE, temporary deployment). To characterize the recording sites in terms of shallow-depth shear-wave velocities (VS), we used both surface waves (Rayleigh and Love) and body waves (S-wave) to evaluate the time-averaged VS in the upper 30 meters of the subsurface (VS30). We used two-dimensional multichannel analysis of surface waves (MASW) to evaluate VS from the surface waves, and a refraction tomography inversion algorithm, developed by Hole in 1992, to evaluate VS from the body waves. As determined by the various methods, we found VS30 near the strongmotion recording stations on Main Street in Downtown Napa, Napa Fire Station Number 3, and on Kreuzer Lane to be from 281 meters per second (m/s) to 286 m/s, 297 to 371 m/s, and 885 to 916 m/s, respectively. The VS30 calculated from Love waves were slightly lower (10 m/s) than those calculated from Rayleigh waves at the Downtown Napa location and at Napa Fire Station Number 3 (4 m/s); however, VS30 calculated from Love waves was higher (190 m/s) than those calculated from Rayleigh waves at Kreuzer Lane. We also found that VS30 determined from MASW for both Love and Rayleigh waves varied depending on the number of shots along the profiles, and VS30 was not systematic based on the number of shots used in the analysis. Furthermore, VS30 calculated from S-wave refraction tomography are closer to those determined from MASW calculated from Love waves than from using Rayleigh waves.