USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29Linked to original sources
USGS2001

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington↗
USGS2002

Deformation associated with the 1997 eruption of Okmok volcano, Alaska

Okmok volcano, located on Umnak Island in the Aleutian chain, Alaska, is the most eruptive caldera system in North America in historic time. Its most recent eruption occurred in 1997. Synthetic aperture radar interferometry shows deflation of the caldera center of up to 140 cm during this time, preceded and followed by inflation of smaller magnitude. The main part of the observed deformation can be modeled using a pressure point source model. The inferred source is located between 2.5 and 5.0 km beneath the approximate center of the caldera and ???5 km from the eruptive vent. We interpret it as a central magma reservoir. The preeruptive period features inflation accompanied by shallow localized subsidence between the caldera center and the vent. We hypothesize that this is caused by hydrothermal activity or that magma moved away from the central chamber and toward the later vent. Since all historic eruptions at Okmok have originated from the same cone, this feature may be a precursor that indicates an upcoming eruption. The erupted magma volume is ???9 times the volume that can be accounted for by the observed preeruptive inflation. This indicates a much longer inflation interval than we were able to observe. The observation that reinflation started shortly after the eruption suggests that inflation spans the whole time interval between eruptions. Extrapolation of the average subsurface volume change rate is in good agreement with the long-term eruption frequency and eruption volumes of Okmok.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Response of the San Andreas fault to the 1983 Coalinga-Nuñez earthquakes: an application of interaction-based probabilities for Parkfield

The Parkfield-Cholame section of the San Andreas fault, site of an unfulfilled earthquake forecast in 1985, is the best monitored section of the world's most closely watched fault. In 1983, the M = 6.5 Coalinga and M = 6.0 Nuñez events struck 25 km northeast of Parkfield. Seismicity rates climbed for 18 months along the creeping section of the San Andreas north of Parkfield and dropped for 6 years along the locked section to the south. Right-lateral creep also slowed or reversed from Parkfield south. Here we calculate that the Coalinga sequence increased the shear and Coulomb stress on the creeping section, causing the rate of small shocks to rise until the added stress was shed by additional slip. However, the 1983 events decreased the shear and Coulomb stress on the Parkfield segment, causing surface creep and seismicity rates to drop. We use these observations to cast the likelihood of a Parkfield earthquake into an interaction-based probability, which includes both the renewal of stress following the 1966 Parkfield earthquake and the stress transfer from the 1983 Coalinga events. We calculate that the 1983 shocks dropped the 10-year probability of a M ∼ 6 Parkfield earthquake by 22% (from 54 ± 22% to 42 ± 23%) and that the probability did not recover until about 1991, when seismicity and creep resumed. Our analysis may thus explain why the Parkfield earthquake did not strike in the 1980s, but not why it was absent in the 1990s. We calculate a 58 ± 17% probability of a M ∼ 6 Parkfield earthquake during 2001–2011.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Constraints on dike propagation from continuous GPS measurements

The January 1997 East Rift Zone eruption on Kilauea volcano, Hawaii, occurred within a network of continuous Global Positioning System (GPS) receivers. The GPS measurements reveal the temporal history of deformation during dike intrusion, beginning ∼8 hours prior to the onset of the eruption. The dike volume as a function of time, estimated from the GPS data using elastic Green's functions for a homogeneous half-space, shows that only two thirds of the final dike volume accumulated prior to the eruption and the rate of volume change decreased with time. These observations are inconsistent with simple models of dike propagation, which predict accelerating dike volume up to the time of the eruption and little or no change thereafter. Deflationary tilt changes at Kilauea summit mirror the inferred dike volume history, suggesting that the rate of dike propagation is limited by flow of magma into the dike. A simple, lumped parameter model of a coupled dike magma chamber system shows that the tendency for a dike to end in an eruption (rather than intrusion) is favored by high initial dike pressures, compressional stress states, large, compressible magma reservoirs, and highly conductive conduits linking the dike and source reservoirs. Comparison of model predictions to the observed dike volume history, the ratio of erupted to intruded magma, and the deflationary history of the summit magma chamber suggest that most of the magma supplied to the growing dike came from sources near to the eruption through highly conductive conduits. Interpretation is complicated by the presence of multiple source reservoirs, magma vesiculation and cooling, as well as spatial variations in dike-normal stress. Reinflation of the summit magma chamber following the eruption was measured by GPS and accompanied a rise in the level of the Pu'u O'o lava lake. For a spheroidal chamber these data imply a summit magma chamber volume of ∼20 km 3 , consistent with recent estimates from seismic tomography. Continuous deformation measurements can be used to image the spatiotemporal evolution of propagating dikes and to reveal quantitative information about the volcanic plumbing systems.

Hawai'i↗
USGS2001

Inferring rate and state friction parameters from a rupture model of the 1995 Hyogo-ken Nanbu (Kobe) Japan earthquake

We consider the applicability of laboratory-derived rate-and state-variable friction laws to the dynamic rupture of the 1995 Kobe earthquake. We analyze the shear stress and slip evolution of Ide and Takeo's [1997] dislocation model, fitting the inferred stress change time histories by calculating the dynamic load and the instantaneous friction at a series of points within the rupture area. For points exhibiting a fast-weakening behavior, the Dieterich-Ruina friction law, with values of d c = 0.01–0.05 m for critical slip, fits the stress change time series well. This range of d c is 10–20 times smaller than the slip distance over which the stress is released, D c , which previous studies have equated with the slip-weakening distance. The limited resolution and low-pass character of the strong motion inversion degrades the resolution of the frictional parameters and suggests that the actual d c is less than this value. Stress time series at points characterized by a slow-weakening behavior are well fitted by the Dieterich-Ruina friction law with values of d c ≥0.01–0.05 m. The apparent fracture energy G c can be estimated from waveform inversions more stably than the other friction parameters. We obtain a G c ≈ 1.5×l0 6 J m −2 for the 1995 Kobe earthquake, in agreement with estimates for previous earthquakes. From this estimate and a plausible upper bound for the local rock strength we infer a lower bound for D c of about 0.008 m.

Kobe, Osaka↗
USGS2001

Steady state deformation of the Coso Range, east central California, inferred from satellite radar interferometry

Observations of deformation from 1992 to 1997 in the southern Coso Range using satellite radar interferometry show deformation rates of up to 35 mm yr −1 in an area ∼10 km by 15 km. The deformation is most likely the result of subsidence in an area around the Coso geothermal field. The deformation signal has a short-wavelength component, related to production in the field, and a long-wavelength component, deforming at a constant rate, that may represent a source of deformation deeper than the geothermal reservoir. We have modeled the long-wavelength component of deformation and inferred a deformation source at ∼4 km depth. The source depth is near the brittle-ductile transition depth (inferred from seismicity) and ∼1.5 km above the top of the rhyolite magma body that was a source for the most recent volcanic eruption in the Coso volcanic field [ Manley and Bacon , 2000]. From this evidence and results of other studies in the Coso Range, we interpret the source to be a leaking deep reservoir of magmatic fluids derived from a crystallizing rhyolite magma body.

California↗
USGS2001

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24-hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real-time or near real-time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground-truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Hawai'i↗
USGS2001

Spatial extent of a hydrothermal system at Kilauea Volcano, Hawaii, determined from array analyses of shallow long-period seismicity 2. Results

Array data from a seismic experiment carried out at Kilauea Volcano, Hawaii, in February 1997, are analyzed by the frequency-slowness method. The slowness vectors are determined at each of three small-aperture seismic antennas for the first arrivals of 1129 long-period (LP) events and 147 samples of volcanic tremor. The source locations are determined by using a probabilistic method which compares the event azimuths and slownesses with a slowness vector model. The results show that all the LP seismicity, including both discrete LP events and tremor, was generated in the same source region along the east flank of the Halemaumau pit crater, demonstrating the strong relation that exists between the two types of activities. The dimensions of the source region are approximately 0.6×1.0×0.5 km. For LP events we are able to resolve at least three different clusters of events. The most active cluster is centered ∼200 m northeast of Halemaumau at depths shallower than 200 m beneath the caldera floor. A second cluster is located beneath the northeast quadrant of Halemaumau at a depth of ∼400 m. The third cluster is <200 m deep and extends southeastward from the northeast quadrant of Halemaumau. Only one source zone is resolved for tremor. This zone is coincident with the most active source zone of LP events, northeast of Halemaumau. The location, depth, and size of the source region suggest a hydrothermal origin for all the analyzed LP seismicity.

Hawai'i↗
USGS2001

Spatial extent of a hydrothermal system at Kilauea Volcano, Hawaii, determined from array analyses of shallow long-period seismicity 1. Method

We present a probabilistic method to locate the source of seismic events using seismic antennas. The method is based on a comparison of the event azimuths and slownesses derived from frequency-slowness analyses of array data, with a slowness vector model. Several slowness vector models are considered including both homogeneous and horizontally layered half-spaces and also a more complex medium representing the actual topography and three-dimensional velocity structure of the region under study. In this latter model the slowness vector is obtained from frequency-slowness analyses of synthetic signals. These signals are generated using the finite difference method and include the effects of topography and velocity structure to reproduce as closely as possible the behavior of the observed wave fields. A comparison of these results with those obtained with a homogeneous half-space demonstrates the importance of structural and topographic effects, which, if ignored, lead to a bias in the source location. We use synthetic seismograms to test the accuracy and stability of the method and to investigate the effect of our choice of probability distributions. We conclude that this location method can provide the source position of shallow events within a complex volcanic structure such as Kilauea Volcano with an error of ±200 m.

Hawai'i↗
USGS2001

Vertical tectonics in northern Escanaba Trough as recorded by thick late Quaternary turbidites

Escanaba Trough, the southernmost segment of the Gorda Ridge, is filled by as much as 500 m of late Quaternary turbidite and hemipelagic sediment. Coring at Deep Sea Drilling Project Site 35 and Ocean Drilling Program (ODP) Sites 1037 and 1038 together with 4.5-kHz deep-tow and 3.5-kHz surface-ship seismic reflection profiles enable a distinct pattern of reflections to be mapped throughout Escanaba Trough in the upper part of this sediment fill. The uppermost 80 m of turbidite sediment, which includes at least 11 turbidity current events, were deposited in <7 kyr. Nine of these turbidites are found throughout Escanaba Trough at water depths >3200 m. The turbidity currents were trapped upon entering Escanaba Trough, resulting in all of the sediment in suspension in the flows being deposited. The thickness of the turbidite layers reflects both the flow thickness and the vertical grain concentration within the flow that deposited the layer. Variations in the turbidite thickness with respect to water depth can be used to estimate the degree of relative vertical movement within the floor of Escanaba Trough. In the area of hydrothermal activity near ODP Site 1038, uplift of as much as 140 m has occurred over the past 8 kyr.

California↗
USGS2001

Viscoelastic shear zone model of a strike-slip earthquake cycle

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

California↗
USGS2002

Lower crustal deformation beneath the central Transverse Ranges, southern California: Results from the Los Angeles Region Seismic Experiment

We present a P wave velocity model derived from active source seismic data collected during the 1994 Los Angeles Region Seismic Experiment. Our model extends previously published upper crustal velocity models to mantle depths. Our model was developed by both ray tracing through a layered model and calculating travel times through a gridded model. It includes an 8-km-thick crustal root centered beneath the surface trace of the San Andreas fault, north of the highest topography in the San Gabriel Mountains. A simple mass balance calculation suggests that ∼36 km of north-south shortening across the San Andreas fault in the central Transverse Ranges could have formed this root. If north-south compression began when the “Big Bend” in the San Andreas fault formed at ∼5 Ma, 36 km of shortening implies a north-south contraction rate of ∼7.1 mm/yr across the central Transverse Ranges. If, instead, north-south compression began when the Transverse Ranges formed at 3.4–3.9 Ma, 36 km of shortening implies a contraction rate of 9.2–10.6 mm/yr. North of the San Andreas fault, the Mojave Desert crust has a low-velocity (6.3 km/s) mid and lower crust and a 28-km-deep Moho. South of the San Andreas fault, beneath the Los Angeles and San Gabriel Valley basins, there is a fast (6.6–6.8 km/s), thick (10–12 km) lower crust with a 27-km-deep Moho. Farther south still, the lower crust of the Continental Borderland is fast (6.6–6.8 km/s) and thin (5 km) with a shallow (22 km deep) Moho.

California↗
USGS2002

Draped aeromagnetic survey in Transantarctic Mountains over the area of the Butcher Ridge igneous complex showing extent of underlying mafic intrusion

A draped aeromagnetic survey over the area surrounding the Butcher Ridge igneous complex (BRIC), Transantarctic Mountains, was acquired in 1997–1998 as part of a larger Transantarctic Mountains Aerogeophysical Research Activity survey. The BRIC is a sill-like hypoabyssal intrusion ranging in composition from tholeiitic basalt to rhyolite. An 40 Ar/ 39 Ar age of 174 Ma and the chemical character of the basaltic rocks show the BRIC to be part of the widespread Jurassic Ferrar suite of continental tholeiitic rocks, that extends for 3500 km across Antarctica. The aeromagnetic survey shows a horseshoe-shaped pattern of anomalies reaching amplitudes as great as 1900 nT generally associated with the bedrock topography where it is exposed. It is apparent that the high-amplitude anomaly pattern is more extensive than the 10-km-long exposed outcrop, first crossed by a single 1990 aeromagnetic profile. The highest-amplitude anomalies appear south of the profile acquired in 1990 and extend out of the survey area. The new aeromagnetic data allow determination of the extent of the interpreted Butcher mafic(?) intrusion beneath exposures of Beacon sedimentary rock and ice in the area covered, as well as beneath the small BRIC exposure. The magnetic anomalies show a minimum area of 3000 km 2 , a much greater extent than previously inferred. Magnetic models indicate a minimum thickness of ∼1–2 km for a horizontal intrusion. However, nonunique models with magnetic layers decreasing in apparent susceptibility with depth are consistent with of a 4- to 8-km-thick layered intrusion. These magnetic models indicate progressively deeper erosion of the interpreted mafic-layered body from the south to north. The erosion has removed more magnetic upper layers that mask the magnetic effects of the lower less magnetic layers. The probable minimum volume of the intrusion in the area of the survey is ∼6000 km 3 . An alternate, but less likely, interpretation of a series of dikes can also fit the observed magnetic anomalies.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Carbon dioxide emission rate of Kīlauea Volcano: Implications for primary magma and the summit reservoir

We report a CO 2 emission rate of 8500 metric tons per day (t d &minus;1 ) for the summit of Kīlauea Volcano, several times larger than previous estimates. It is based on three sets of measurements over 4 years of synchronous SO 2 emission rates and volcanic CO 2 /SO 2 concentration ratios for the summit correlation spectrometer (COSPEC) traverse. Volcanic CO 2 /SO 2 for the traverse is representative of the global ratio for summit emissions. The summit CO 2 emission rate is nearly constant, despite large temporal variations in summit CO 2 /SO 2 and SO 2 emission rates. Summit CO 2 emissions comprise most of Kīlauea's total CO 2 output (&sim;9000 t d &minus;1 ). The bulk CO 2 content of primary magma determined from CO 2 emission and magma supply rate data is &sim;0.70 wt %. Most of the CO 2 is present as exsolved vapor at summit reservoir depths, making the primary magma strongly buoyant. Turbulent mixing with resident reservoir magma, however, prevents frequent eruptions of buoyant primary magma in the summit region. CO 2 emissions confirm that the magma supply enters the edifice through the summit reservoir. A persistent several hundred parts per million CO 2 anomaly arises from the entry of magma into the summit reservoir beneath a square kilometer area east of Halemaumau pit crater. Since most of the CO 2 in primary magma is degassed in the summit, the summit CO 2 emission rate is an effective proxy for the magma supply rate. Both scrubbing of SO 2 and solubility controls on CO 2 and S in basaltic melt cause high CO 2 /SO 2 in summit emissions and spatially uncorrelated distributions of CO 2 and SO 2 in the summit plume.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Magmatic inflation at a dormant stratovolcano: 1996-1998 activity at Mount Peulik volcano, Alaska, revealed by satellite radar interferometry

A series of ERS radar interferograms that collectively span the time interval from July 1992 to August 2000 reveal that a presumed magma body located 6.6 ??? 0.5 km beneath the southwest flank of the Mount Peulik volcano inflated 0.051 ??? 0.005 km3 between October 1996 and September 1998. Peulik has been active only twice during historical time, in 1814 and 1852, and the volcano was otherwise quiescent during the 1990s. The inflation episode spanned at least several months because separate interferograms show that the associated ground deformation was progressive. The average inflation rate of the magma body was ???0.003 km3/month from October 1996 to September 1997, peaked at 0.005 km3/month from 26 June to 9 October 1997, and dropped to ???0.001 km3/month from October 1997 to September 1998. An intense earthquake swarm, including three ML 4.8 - 5.2 events, began on 8 May 1998 near Becharof Lake, ???30 km northwest of Peulik. More than 400 earthquakes with a cumulative moment of 7.15 ?? 1017 N m were recorded in the area through 19 October 1998. Although the inflation and earthquake swarm occurred at about the same time, the static stress changes that we calculated in the epicentral area due to inflation beneath Peulik appear too small to provide a causal link. The 1996-1998 inflation episode at Peulik confirms that satellite radar interferometry can be used to detect magma accumulation beneath dormant volcanoes at least several months before other signs of unrest are apparent. This application represents a first step toward understanding the eruption cycle at Peulik and other stratovolcanoes with characteristically long repose periods.

Alaska↗
USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗
USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 2. Hydrogeologic implications

As part of a regional groundwater investigation, geophysical logs were obtained in two municipal water wells located near the west Texas city of Alpine. These boreholes are 252 and 285 m deep and penetrate extrusive rocks of Tertiary age. The deeper well was drilled in the central valley and the other along the northern flank of an east-west trending valley-ridge setting. Analysis and interpretation of the logs reveal that the two wells are subjected to significantly different stress environments because of topographic effects and exhibit significantly different hydrogeologic properties. Water production is associated with two specific types of features common to both wells: (1) the upper and lower contacts of a dense trachyte unit located in the shallow part of the wells and (2) deeper zones of highly fractured rocks within the interior of a basalt formation. The transmissivity of the trachyte boundaries is twice as large in the central valley well as it is in the ridge flank well, whereas the transmissivity of the deeper basalts is an order of magnitude greater in the flank well than it is in the central well. This discrepancy is examined from the perspective of rock failure, fracture opening, and flow enhancement by computing values for a Drucker-Prager stability factor that is based on the magnitudes of the normal and deviatoric stress invariants as a function of depth. Thus the field measurements and subsequent stress analysis offer evidence of a coupled tectonic-hydrologic interaction at this site.

Texas↗