James F. Howle

Water-level changes induced by local and distant earthquakes at Long Valley caldera, California

Distant as well as local earthquakes have induced groundwater-level changes persisting for days to weeks at Long Valley caldera, California. Four wells open to formations as deep as 300 m have responded to 16 earthquakes, and responses to two earthquakes in the 3-km-deep Long Valley Exploratory Well (LVEW) show that these changes are not limited to weathered or unconsolidated near-surface rocks. All five wells exhibit water-level variations in response to earth tides, indicating they can be used as low-resolution strainmeters. Earthquakes induce gradual water-level changes that increase in amplitude for as long as 30 days, then return more slowly to pre-earthquake levels. The gradual water-level changes are always drops at wells LKT, LVEW, and CH-10B, and always rises at well CW-3. At a dilatometer just outside the caldera, earthquake-induced strain responses consist of either a step followed by a contractional strain-rate increase, or a transient contractional signal that reaches a maximum in about seven days and then returns toward the pre-earthquake value. The sizes of the gradual water-level changes generally increase with earthquake magnitude and decrease with hypocentral distance. Local earthquakes in Long Valley produce coseismic water-level steps; otherwise the responses to local earthquakes and distant earthquakes are indistinguishable. In particular, water-level and strain changes in Long Valley following the 1992 M7.3 Landers earthquake, 450 km distant, closely resemble those initiated by a M4.9 local earthquake on November 22, 1997, during a seismic swarm with features indicative of fluid involvement. At the LKT well, many of the response time histories are identical for 20 days after each earthquake, and can be matched by a theoretical solution giving the pore pressure as a function of time due to diffusion of a nearby, instantaneous, pressure drop. Such pressure drops could be produced by accelerated inflation of the resurgent dome by amounts too small to be detected by the two-color electronic distance-measuring network. Opening-mode displacement in the south moat, inferred to have followed a M4.9 earthquake on November 22, 1997, could also create extensional strain on the dome and lead to water-level changes similar to those following dome inflation. Contractional strain that could account for earthquake-induced water-level rises at the CW-3 well is inconsistent with geodetic observations. We instead attribute these water-level rises to diffusion of elevated fluid pressure localized in the south moat thermal aquifer. For hydraulic diffusivities appropriate to the upper few hundred meters at Long Valley, an influx of material at temperatures of 300°C can thermally generate pressure of 6 m of water or more, an order of magnitude larger than needed to account for the CW-3 water-level rises. If magma or hot aqueous fluid rises to within 1 km of the surface in the eastern part of the south moat, then hydraulic diffusivities are high enough to allow fluid pressure to propagate to CW-3 on the time scale observed. The data indicate that seismic waves from large distant earthquakes can stimulate upward movement of fluid in the hydrothermal system at Long Valley.

LiDAR-Assisted Identification of an Active Fault near Truckee, California

Abstract We use high-resolution (1.5–2.4 points/m 2 ) bare-earth airborne Light Detection and Ranging (LiDAR) imagery to identify, map, constrain, and visualize fault-related geomorphology in densely vegetated terrain surrounding Martis Creek Dam near Truckee, California. Bare-earth LiDAR imagery reveals a previously unrecognized and apparently youthful right-lateral strike-slip fault that exhibits laterally continuous tectonic geomorphic features over a 35-km-long zone. If these interpretations are correct, the fault, herein named the Polaris fault, may represent a significant seismic hazard to the greater Truckee–Lake Tahoe and Reno–Carson City regions. Three-dimensional modeling of an offset late Quaternary terrace riser indicates a minimum tectonic slip rate of 0.4±0.1 mm/yr. Mapped fault patterns are fairly typical of regional patterns elsewhere in the northern Walker Lane and are in strong coherence with moderate magnitude historical seismicity of the immediate area, as well as the current regional stress regime. Based on a range of surface-rupture lengths and depths to the base of the seismogenic zone, we estimate a maximum earthquake magnitude ( M ) for the Polaris fault to be between 6.4 and 6.9.

Deformation near the Casa Diablo geothermal well field and related processes Long Valley caldera, Eastern California, 1993–2000

Abstract Regional first-order leveling lines, which extend from Lee Vining, CA, to Tom’s Place, CA, have been surveyed periodically since 1957 by the U.S. Geological Survey (USGS), the National Geodetic Survey (NGS), and Caltrans. Two of the regional survey lines, or leveling networks, intersect at the Casa Diablo geothermal well field. These leveling networks, referenced to a distant bench mark (C916) near Lee Vining, provide time-series vertical control data of land-surface deformation that began around 1980. These data are also useful for delineating localized subsidence at Casa Diablo related to reservoir pressure and temperature changes owing to geothermal development that began in 1985. A comparison of differences in bench-mark elevations for five time periods between 1983 and 1997 shows the development and expansion of a subsidence bowl at Casa Diablo. The subsidence coincides spatially with the geothermal well field and temporally with the increased production rates and the deepening of injection wells in 1991, which resulted in an increase in the rate of pressure decline. The subsidence, superimposed on a broad area of uplift, totaled about 310 mm by 1997. The USGS established orthogonal tilt arrays in 1983 to better monitor deformation across the caldera. One tilt array (DBR) was established near what would later become the Casa Diablo geothermal well field. This array responded to magmatic intrusions prior to geothermal development, tilting away from the well field. With the start of geothermal fluid extraction in 1985, tilt at the DBR array reversed direction and began tilting into the well field. In 1991, geothermal power production was increased by a factor of four, and reservoir pressures began a period of steep decline. These changes caused a temporary three-fold increase in the tilt rate. The tilt rate became stable in 1993 and was about 40% lower than that measured in 1991–1992, but still greater than the rates measured during 1985–1990. Data from the local leveling networks spanning the well field and the bounding graben were analyzed for several 2-year periods (1993–1995, 1995–1997, and 1997–1999). Annual rates of change across the normal faults bounding the graben have steadily decreased for each 2-year period between 1993 and 1999, reflecting the slowing decline in geothermal reservoir pressure. Horizontal control data from a two-color electronic distance meter (EDM) defined the lateral extent of subsidence at Casa Diablo. The EDM and leveling data elucidate the localized effect of the shallow source of subsidence and the broader effect of the deeper magmatic inflation source. Data from bench marks common to both the vertical and the horizontal control networks were used to assess the effect of subsidence on the EDM base station (CASA). Modeling of geodetic data collected during periods of little or no magmatic inflation indicated that the CASA two-color EDM station is being drawn toward the well field at a rate of 3–5 mm/yr.

Airborne LiDAR analysis and geochronology of faulted glacial moraines in the Tahoe-Sierra frontal fault zone reveal substantial seismic hazards in the Lake Tahoe region, California-Nevada, USA

We integrated high-resolution bare-earth airborne light detection and ranging (LiDAR) imagery with field observations and modern geochronology to characterize the Tahoe-Sierra frontal fault zone, which forms the neotectonic boundary between the Sierra Nevada and the Basin and Range Province west of Lake Tahoe. The LiDAR imagery clearly delineates active normal faults that have displaced late Pleistocene glacial moraines and Holocene alluvium along 30 km of linear, right-stepping range front of the Tahoe-Sierra frontal fault zone. Herein, we illustrate and describe the tectonic geomorphology of faulted lateral moraines. We have developed new, three-dimensional modeling techniques that utilize the high-resolution LiDAR data to determine tectonic displacements of moraine crests and alluvium. The statistically robust displacement models combined with new ages of the displaced Tioga (20.8 ± 1.4 ka) and Tahoe (69.2 ± 4.8 ka; 73.2 ± 8.7 ka) moraines are used to estimate the minimum vertical separation rate at 17 sites along the Tahoe-Sierra frontal fault zone. Near the northern end of the study area, the minimum vertical separation rate is 1.5 ± 0.4 mm/yr, which represents a two- to threefold increase in estimates of seismic moment for the Lake Tahoe basin. From this study, we conclude that potential earthquake moment magnitudes (M w ) range from 6.3 ± 0.25 to 6.9 ± 0.25. A close spatial association of landslides and active faults suggests that landslides have been seismically triggered. Our study underscores that the Tahoe-Sierra frontal fault zone poses substantial seismic and landslide hazards.

Inferences on the hydrothermal system beneath the resurgent dome in Long Valley Caldera, east-central California, USA, from recent pumping tests and geochemical sampling

Abstract Quaternary volcanic unrest has provided heat for episodic hydrothermal circulation in the Long Valley caldera, including the present-day hydrothermal system, which has been active over the past 40 kyr. The most recent period of crustal unrest in this region of east-central California began around 1980 and has included periods of intense seismicity and ground deformation. Uplift totaling more than 0.7 m has been centered on the caldera’s resurgent dome, and is best modeled by a near-vertical ellipsoidal source centered at depths of 6–7 km. Modeling of both deformation and microgravity data now suggests that (1) there are two inflation sources beneath the caldera, a shallower source 7–10 km beneath the resurgent dome and a deeper source ∼15 km beneath the caldera’s south moat and (2) the shallower source may contain components of magmatic brine and gas. The Long Valley Exploration Well (LVEW), completed in 1998 on the resurgent dome, penetrates to a depth of 3 km directly above this shallower source, but bottoms in a zone of 100°C fluid with zero vertical thermal gradient. Although these results preclude extrapolations of temperatures at depths below 3 km, other information obtained from flow tests and fluid sampling at this well indicates the presence of magmatic volatiles and fault-related permeability within the metamorphic basement rocks underlying the volcanic fill. In this paper, we present recently acquired data from LVEW and compare them with information from other drill holes and thermal springs in Long Valley to delineate the likely flow paths and fluid system properties under the resurgent dome. Additional information from mineralogical assemblages in core obtained from fracture zones in LVEW documents a previous period of more vigorous and energetic fluid circulation beneath the resurgent dome. Although this system apparently died off as a result of mineral deposition and cooling (and/or deepening) of magmatic heat sources, flow testing and tidal analyses of LVEW water level data show that relatively high permeability and strain sensitivity still exist in the steeply dipping principal fracture zone penetrated at a depth of 2.6 km. The hydraulic properties of this zone would allow a pressure change induced at distances of several kilometers below the well to be observable within a matter of days. This indicates that continuous fluid pressure monitoring in the well could provide direct evidence of future intrusions of magma or high-temperature fluids at depths of 5–7 km.

Electrical resistivity investigation of fluvial geomorphology to evaluate potential seepage conduits to agricultural lands along the San Joaquin River, Merced County, California, 2012–13

.........................................................................................................................................................