Geology topics

W. Thatcher

Publications and source records attributed to W. Thatcher.

At least 19 recordsLinked to original sources

USGS2012

Lithospheric loading by the 1896 Riku-u earthquake, northern Japan: Implications for plate flexure and asthenospheric rheology

Under favorable circumstances the time-dependent aseismic deformation resulting from the loading of the lithosphere by the stress drop of large dip slip earthquakes can be used to determine both the effective elastic plate thickness and the asthenospheric viscosity. The deformation has several similarities with the deflection of the lithosphere by surface loads and with movements due to postglacial rebound. Level changes obtained in the 80 years since the M = 7.5, 1896 Riku-u earthquake, an intraplate thrust event in northern Honshu, provide convincing evidence that asthenospheric readjustments are responsible for the observed movements. Leveling surveys crossing the zone of surface faulting have been repeated five times since 1900 and delineate a localized depression that has subsided at a continually decreasing rate. The depression is centered close to the 1896 faulting, and its shape and width, about 75 km, are matched by our model using a plate thickness of 30 km. The decaying subsidence rate constrains the viscosity of the uppermost asthenosphere to be 1×10 20 P. A linear viscous rheology matches the observed decay quite well, although measurements are sparse during the several decades following the earthquake.

Journal of Geophysical Research Solid Earth↗

USGS2012

Quaternary and geodetically measured crustal movements in the Tokai District, Central Honshu, Japan

Geodetic measurements made in this century demonstrate that modern movements correlate well with deformation that has occurred during the past 2 m.y. Both show effects of subduction and of the Quaternary collision of Izu Peninsula with Central Honshu. Westward from Suruga Bay, deformation is greatest near the coast, decreases markedly towards the uplifting Akaishi Mountains, and then increases again in the conjugate block-faulted terranes of Chubu district. Quaternary faults, historic earthquakes, and the distribution of horizontal shear straining since ∼1880 all reflect this deformation pattern. Since 1900 the west coast of Suruga Bay has undergone secular subsidence, horizontal contraction and tilting towards Suruga Trough, all typical of the strain buildup process that occurs landward of subduction zones. However, futher westward, active mountain building is occuring, and the Akaishi Ranges are rising 3 mm/a. Unlike the episodic aseismic uplift of the Transverse Ranges of southern California, movements here have been relatively uniform. It seems likely that an increase in the tilt rate of the Suruga coast recognized since ∼1973 is part of a long-term readjustment to the great 1944 Tonankai earthquake. Data show that a monotonic increase has occurred since at least 1949, and similar effects are seen adjacent to the 1944 rupture zone over 200 km southwest along the Nankai Trough.

Journal of Geophysical Research Solid Earth↗

USGS2012

Nonlinear strain buildup and the earthquake cycle on the San Andreas fault

Two contrasting models of the earthquake deformation cycle on strike slip faults predict significant temporal declines in shear strain rate near the fault, accompanied by a progressive broadening of the zone of deformation adjacent to it. In the thin lithosphere model, transient deformation results from flow in the asthenosphere due to stress relaxation following faulting through most or all of the lithosphere. For an earth model with a thick elastic lithosphere (plate thickness » depth of seismic slip), transient motions are due to postearthquake aseisrnic slip below the coseismic fault plane. Data from the San Andreas fault indicate a long-term temporal decrease in strain rate that persists for at least 30 years and may extend through the entire earthquake cycle. Observations support a cycle-long rate decrease and a temporal spreading of the deformation profile only if movement cycles on the northern and southern locked sections of the fault are basically similar. If so, the usually lower strain rates and broader deformation zone currently observed on the southern San Andreas represent a later evolutionary stage of the northern locked section, where a great earthquake is a more recent occurrence. Although the data allow some extreme models to be discarded, no sufficiently strong constraints exist to decide between the thin and thick lithosphere models. Regardless of the appropriate model the geodetic observations themselves indicate that strain buildup is sufficiently nonlinear to cause significant departures from recurrence estimates based on linear strain accumulation and the time-predictable model.

Journal of Geophysical Research Solid Earth↗

USGS2012

Order and diversity in the modes of circum-Pacific earthquake recurrence

Recurrence characteristics of great circum-Pacific earthquakes and determinations of spatial distribution of seismic moment release are surveyed in order to delineate their general features and relate earthquake slip distribution to models of recurrent rupture. As noted by others, the pattern of moment release is typically very irregular, with strong concentrations in a few isolated regions of a much larger aftershock zone. Despite this complexity, rupture nucleation is notably systematic, with mainshock epicenter showing a strong tendency to locate in or immediately adjacent to identified regions of high moment release. This generalization suggests that earthquake recurrence is more likely to be controlled by maximum rather than average fault slip. Well-documented case histories from twelve plate boundary regions demonstrate that seismic strain release tends to be temporally well-ordered, while source dimensions, slip, and cumulative moment release vary considerably from cycle to cycle. On none of the segments studied do earthquakes consistently recur in nearly identical events. Instead, individual great earthquakes differ significantly from cycle to cycle or rerupture takes place in a sequence of two or more smaller events. Despite these differences the duration of the seismic cycle is approximately uniform. Furthermore, when strain release occurs in a sequence of large earthquakes, these events take place towards the end of the cycle and occur with increasing rupture length and magnitude. These ordered and irregular features of earthquake recurrence argue for the existence of corresponding elements on plate boundary faults. The ordered characteristics are identified with zones of concentrated moment release that slip comparable amounts in each cycle and have high shear strength. The irregular features are associated with intervening weaker regions that move in response to the stress concentration of dynamic rupture and slip by differing amounts in each event. Although these mechanistic associations are indirect and tentative, the features of recurrent behavior documented here have implications for long-term earthquake hazard assessment that are not dependent on the models proposed to explain them. Chief among these are the cycle-to-cycle differences among gap-filling events and the absence of shocks that fill major slip-deficient regions early in the seismic cycle.

Journal of Geophysical Research Solid Earth↗

USGS2011

Diffuse Pacific-North American plate boundary: 1000 km of dextral shear inferred from modeling geodetic data

Geodetic measurements tell us that the eastern part of the Basin and Range Province expands in an east-west direction relative to stable North America, whereas the western part of the province moves to the northwest. We develop three-dimensional finite element representations of the western United States lithosphere in an effort to understand the global positioning system (GPS) signal. The models are constrained by known bounding-block velocities and topography, and Basin and Range Province deformation is represented by simple plastic (thermal creep) rheology. We show that active Basin and Range spreading by gravity collapse is expected to have a strong southward component that does not match the GPS signal. We can reconcile the gravitational component of displacement with observed velocity vectors if the Pacific plate applies northwest-directed shear stress to the Basin and Range via the Sierra Nevada block. This effect reaches at least 1000 km east of the San Andreas fault in our models.

Nevada↗

USGS2008

Temporal evolution of continental lithospheric strength in actively deforming regions

It has been agreed for nearly a century that a strong, load-bearing outer layer of earth is required to support mountain ranges, transmit stresses to deform active regions and store elastic strain to generate earthquakes. However the dept and extent of this strong layer remain controversial. Here we use a variety of observations to infer the distribution of lithospheric strength in the active western United States from seismic to steady-state time scales. We use evidence from post-seismic transient and earthquake cycle deformation reservoir loading glacio-isostatic adjustment, and lithosphere isostatic adjustment to large surface and subsurface loads. The nearly perfectly elastic behavior of Earth's crust and mantle at the time scale of seismic wave propagation evolves to that of a strong, elastic crust and weak, ductile upper mantle lithosphere at both earthquake cycle (EC, ???10?? to 103 yr) and glacio-isostatic adjustment (GIA, ???103 to 104 yr) time scales. Topography and gravity field correlations indicate that lithosphere isostatic adjustment (LIA) on ???106-107 yr time scales occurs with most lithospheric stress supported by an upper crust overlying a much weaker ductile subtrate. These comparisons suggest that the upper mantle lithosphere is weaker than the crust at all time scales longer than seismic. In contrast, the lower crust has a chameleon-like behavior, strong at EC and GIA time scales and weak for LIA and steady-state deformation processes. The lower crust might even take on a third identity in regions of rapid crustal extension or continental collision, where anomalously high temperatures may lead to large-scale ductile flow in a lower crustal layer that is locally weaker than the upper mantle. Modeling of lithospheric processes in active regions thus cannot use a one-size-fits-all prescription of rheological layering (relation between applied stress and deformation as a function of depth) but must be tailored to the time scale and tectonic setting of the process being investigated.

GSA Today↗

USGS2007

Microplate model for the present-day deformation of Tibet

[1] Site velocities from 349 Global Positioning System (GPS) stations are used to construct an 11-element quasi-rigid block model of the Tibetan Plateau and its surroundings. Rigid rotations of five major blocks are well determined, and average translation velocities of six smaller blocks can be constrained. Where data are well distributed the velocity field can be explained well by rigid block motion and fault slip across block boundaries. Residual misfits average 1.6 mm/yr compared to typical one standard deviation velocity uncertainties of 1.3 mm/yr. Any residual internal straining of the blocks is small and heterogeneous. However, residual substructure might well represent currently unresolved motions of smaller blocks. Although any smaller blocks must move at nearly the same rate as the larger blocks within which they lie, undetected relative motions between them could be significant, particularly where there are gaps in GPS coverage. Predicted relative motions between major blocks agree with the observed sense of slip and along-strike partitioning of motion across major faults. However, predicted slip rates across Tibet's major strike-slip faults are low, only 5–12 mm/yr, a factor of 2–3 smaller than most rates estimated from fault offset features dated by radiometric methods as ∼2000 to ∼100,000 year old. Previous work has suggested that both GPS data and low fault slip rates are incompatible with rigid block motions of Tibet. The results reported here overcome these objections.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

The 1923 Kanto earthquake reevaluated using a newly augmented geodetic data set

This study revisits the mechanism of the 1923 Ms = 7.9 Kanto earthquake in Japan. We derive a new source model and use it to assess quantitative and qualitative aspects of the accommodation of plate motion in the Kanto region. We use a new geodetic data set that consists of displacements from leveling and angle changes from triangulation measurements obtained in surveys between 1883 and 1927. Two unique aspects of our analysis are the inclusion of a large number of second-order triangulation measurements and the application of a correction to remove interseismic deformation. The geometry of the fault planes is adopted from a recent seismic reflection study of the Kanto region. We evaluate the minimum complexity necessary in the model to fit the data optimally. Our final uniform-slip elastic dislocation model consists of two adjacent ???20?? dipping low-angle planes accommodating reverse dextral slip of 6.0 in on the larger, eastern plane and 9.5 m on the smaller, western plane with azimuths of 163?? and 121??, respectively. The earthquake was located in the Sagami trough, where the Philippine Sea plate subducts under Honshu. Compared to the highly oblique angle of plate convergence, the coseismic slip on the large fault plane has a more orthogonal orientation to the strike of the plate boundary, suggesting that slip partitioning plays a role in accommodation of plate motion. What other structure is involved in the partitioning is unclear. Uplift records of marine coastal terraces in Sagami Bay document 7500 years of earthquake activity and predict average recurrence intervals of 400 years for events with vertical displacement profiles similar to those of the 1923 earthquake. This means that the average slip deficit per recurrence interval is ???50% of the relative plate convergence. These findings of plate motion partitioning and slip deficit lead us to suggest that instead of a simple recurrence model with characteristic earthquakes, additional mechanisms are necessary to describe the accommodation of deformation in the Kanto region. So far, obvious candidates for these alternative mechanisms have not been discovered. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Inference of postseismic deformation mechanisms of the 1923 Kanto earthquake

Coseismic slip associated with the M7.9, 1923 Kanto earthquake is fairly well understood, involving slip of up to 8 m along the Philippine Sea-Honshu interplate boundary under Sagami Bay and its onland extension. Postseismic deformation after the 1923 earthquake, however, is relatively poorly understood. We revisit the available deformation data in order to constrain possible mechanisms of postseismic deformation and to examine the consequences for associated stress changes in the surrounding crust. Data from two leveling lines and one tide gage station over the first 7-8 years postseismic period are of much greater amplitude than the corresponding expected interseismic deformation during the same period, making these data suitable for isolating the signal from postseismic deformation. We consider both viscoelastic models of asthenosphere relaxation and afterslip models. A distributed coseismic slip model presented by Pollitz et al. (2005), combined with prescribed parameters of a viscoelastic Earth model, yields predicted postseismic deformation that agrees with observed deformation on mainland Honshu from Tokyo to the Izu peninsula. Elsewhere (southern Miura peninsula; Boso peninsula), the considered viscoelastic models fail to predict observed deformation, and a model of ???1 in shallow afterslip in the offshore region south of the Boso peninsula, with equivalent moment magnitude Mw = 7.0, adequately accounts for the observed deformation. Using the distributed coseismic slip model, layered viscoelastic structure, and a model of interseismic strain accumulation, we evaluate the post-1923 stress evolution, including both the coseismic and accumulated postseismic stress changes and those stresses contributed by interseismic loading. We find that if account is made for the varying tectonic regime in the region, the occurrence of both immediate (first month) post-1923 crustal aftershocks as well as recent regional crustal seismicity is consistent with the predicted stress pattern. This suggests that the influence of the 1923 earthquake on regional seismicity is fairly predictable and has persisted for at least seven decades following the earthquake.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Northwest Basin and Range tectonic deformation observed with the Global Positioning System, 1999-2003

We use geodetic velocities obtained with the Global Positioning System (GPS) to quantify tectonic deformation of the northwest Basin and Range province of the western United States. The results are based on GPS data collected in 1999 and 2003 across five new quasi-linear networks in northern Nevada, northeast California, and southeast Oregon. The velocities show ???3 mm/yr westward movement of northern Nevada with respect to stable North America. West of longitude 119??W the velocities increase and turn northwest, parallel to Sierra Nevada/Great Valley microplate motion, and similar to velocities previously obtained to the south. The observations are explained by a kinematic model with three domains that rotate around Euler poles in eastern Oregon and western Idaho. Northeast California experiences internal dextral shear deformation (11.2 ?? 3.6 nstrain/yr) subparallel to Pacific/North America motion. Relative motions of the domains imply 2-5 mm/yr approximately east-west extension in northwest Nevada and 1-4 mm/yr approximately north-south contraction near the California/Oregon border. The northward decreasing approximately east-west extension in northwest Nevada is consistent with the northern termination of Basin and Range deformation, faulting and characteristic topography. No significant extension is detected in the Oregon Basin and Range. The Oregon Cascade arc moves north at ???3.5 mm/yr and is possibly influenced by the approximately eastward motion of the Juan de Fuca plate. These results disagree with secular northwest trenchward motion of the Oregon forearc inferred from paleomagnetic rotations. South of latitude 43??, however, trenchward motion exists and is consistent with block rotations, approximately east-west Basin and Range extension, and northwest Sierra Nevada translation. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

New constraints on the active tectonic deformation of the Aegean

Site velocities from six separate Global Positioning System (GPS) networks comprising 374 stations have been referred to a single common Eurasia-fixed reference frame to map the velocity distribution over the entire Aegean. We use the GPS velocity field to identify deforming regions, rigid elements, and potential microplate boundaries, and build upon previous work by others to initially specify rigid elements in central Greece, the South Aegean, Anatolia, and the Sea of Marmara. We apply an iterative approach, tentatively defining microplate boundaries, determining best fit rigid rotations, examining misfit patterns, and revising the boundaries to achieve a better match between model and data. Short-term seismic cycle effects are minor contaminants of the data that we remove when necessary to isolate the long-term kinematics. We find that present day Aegean deformation is due to the relative motions of four microplates and straining in several isolated zones internal to them. The RMS misfit of model to data is about 2-sigma, very good when compared to the typical match between coseismic fault models and GPS data. The simplicity of the microplate description of the deformation and its good fit to the GPS data are surprising and were not anticipated by previous work, which had suggested either many rigid elements or broad deforming zones that comprise much of the Aegean region. The isolated deforming zones are also unexpected and cannot be explained by the kinematics of the microplate motions. Strain rates within internally deforming zones are extensional and range from 30 to 50 nanostrain/year (nstrain/year, 10-9/year), 1 to 2 orders of magnitude lower than rates observed across the major microplate boundaries. Lower strain rates may exist elsewhere withi the microplates but are only resolved in Anatolia, where extension of 13 ?? 4 nstrain/ year is required by the data. Our results suggest that despite the detailed complexity of active continental deformation revealed by seismicity, active faulting, fault geomorphology, and earthquake fault plane solutions, continental tectonics, at least in the Aegean, is to first order very similar to global plate tectonics and obeys the same simple kinematic rules. Although the widespread distribution of Aegean seismicity and active faulting might suggest a rather spatially homogeneous seismic hazard, the focusing of deformation near microplate boundaries implies the highest hazard is comparably localized.

Journal of Geophysical Research B: Solid Earth↗

USGS2004

Contemporary tectonic deformation of the Basin and Range province, western United States: 10 years of observation with the Global Positioning System

[1] We have estimated patterns and rates of crustal movement across 800 km of the Basin and Range at ∼39° north latitude with Global Positioning System surveys in 1992, 1996, 1998, and 2002. The total rate of motion tangent to the small circle around the Pacific‐North America pole of rotation is 10.4 ± 1.0 mm/yr, and motion normal to this small circle is 3.9 ± 0.9 mm/yr compared to the east end of our network. On the Colorado Plateau the east end of our network moves by ∼1–2 mm/yr westerly with respect to North America. Transitions in strain rates delimit six major tectonic domains within the province. These deformation zones coincide with areas of modern seismicity and are, from east to west, (1) east‐west extension in the Wasatch Fault zone, (2) low rate east‐west extension centered near the Nevada‐Utah border, (3) low rate east‐west contraction between 114.7°W and 117.9°W, (4) extension normal to and strike‐slip motion across the N10°E striking Central Nevada Seismic Zone, (5) right lateral simple shear oriented N13°W inside the Walker Lane Belt, and (6) shear plus extension near the Sierra Nevada frontal faults. Concentration of shear and dilatational deformation across the three westernmost zones suggests that the Walker Lane Belt lithosphere is rheologically weak. However, we show that linear gradients in viscosity and gravitational potential energy can also effectively concentrate deformation. In the Basin and Range, gradients in gravitational potential are spatially anticorrelated with dilatational strain rates, consistent with the presence of horizontal variations in viscosity of the lithosphere.

Arizona, California, Nevada, Utah↗

USGS2004

Systematic variation in the depths of slabs beneath arc volcanoes

The depths to the tops of the zones of intermediate-depth seismicity beneath arc volcanoes are determined using the hypocentral locations of Engdahl et al. These depths are constant, to within a few kilometres, within individual arc segments, but differ by tens of kilometres from one arc segment to another. The range in depths is from 65 km to 130 km, inconsistent with the common belief that the volcanoes directly overlie the places where the slabs reach a critical depth that is roughly constant for all arcs. The depth to the top of the intermediate-depth seismicity beneath volcanoes correlates neither with age of the descending ocean floor nor with the thermal parameter of the slab. This depth does, however, exhibit an inverse correlation with the descent speed of the subducting plate, which is the controlling factor both for the thermal structure of the wedge of mantle above the slab and for the temperature at the top of the slab. We interpret this result as indicating that the location of arc volcanoes is controlled by a process that depends critically upon the temperature at the top of the slab, or in the wedge of mantle, immediately below the volcanic arc.

Geophysical Journal International↗

USGS2003

Rheology of the lithosphere inferred from postseismic uplift following the 1959 Hebgen Lake earthquake

We have modeled the broad postseismic uplift measured by geodetic leveling in the epicentral area of the 1959 Mw = 7.3 Hebgen Lake, Montana earthquake, a normal faulting event in the northern Basin and Range province. To fit the observed uplift we calculate synthetic postseismic deformation using the relaxation response of a gravitational viscoelastic Earth to the earthquake. For a model with an elastic plate overlying a viscoelastic half-space, we find that the elastic thickness is 38 ?? 8 km, which isclose to the local crustal thickness. The half-space viscosity is estimated at 4 ?? 1018??0.5 Pa s. The leveling data do not require a viscous lower crust but permit a lower bound viscosity of 1020 Pa s. The observed broad uplift cannot be explained by physically plausible afterslip on and below the coseismic fault. However, local deformation across the coseismic surface rupture requires shallow afterslip reaching the surface. The postseismic deformation induced by the estimated viscoelastic structure decays exponentially with a time constant of ???15 years. Because of coupling between the elastic layer and the viscoelastic substrate, this relaxation time is significantly longer than the 2 year Maxwell relaxation time of the viscous half-space itself. Our result suggests the importance of postseismic relaxation in interpreting high-precision global positioning system velocities. For example, our model results suggest that postseismic transient velocities from both the 1959 Hebgen Lake and the 1983 Mw = 6.9 Borah Peak earthquakes are currently as large as 1-2 mm/yr.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

GPS constraints on the kinematics of continental deformation

Recent GPS observations from the western United States, New Zealand, central Greece, and Japan indicate that present-day continental deformation is typically focused in narrow deforming zones whose extent is much smaller than the intervening largely inactive regions. However, these narrow zones are heterogeneously distributed, reflecting the inherent heterogeneity of continental lithospheric strength and internal buoyancy. Plate driving and resisting forces stress plate boundary zones and plate interiors and drive deformation. These forces change continuously and discontinuously, leading to continental deformation that typically evolves and migrates with time. Magmatic and tectonic processes alter lithospheric rheology and internal buoyancy and also contribute to the time-varying character of continental deformation.

International Geology Review↗

USGS2002

Interferometric synthetic aperture radar studies of Alaska volcanoes

Interferometric synthetic aperture radar (InSAR) imaging is a recently developed geodetic technique capable of measuring ground-surface deformation with centimeter to subcentimeter vertical precision and spatial resolution of tens-of-meter over a relatively large region (/spl sim/10/sup 4/ km/sup 2/). The spatial distribution of surface deformation data, derived from InSAR images, enables the construction of detailed mechanical models to enhance the study of magmatic and tectonic processes associated with volcanoes. This paper summarizes our recent InSAR studies of several Alaska volcanoes, which include Okmok, Akutan, Kiska, Augustine, Westdahl, and Peulik volcanoes.

Alaska↗

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↗