SEARCH · Search USGS

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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39Linked to original sources

USGS1997

Coseismic deformation of the Wrights tunnel during the 1906 San Francisco earthquake: A key to understanding 1906 fault slip and 1989 surface ruptures in the southern Santa Cruz Mountains, California

The Wrights tunnel is an abandoned railroad tunnel that crosses the San Andreas fault in the southern Santa Cruz Mountains in the vicinity of the 1989 Loma Prieta earthquake. The tunnel was damaged and deformed during the 1906 San Francisco earthquake and a plot showing postearthquake measurements made in the tunnel is given by Lawson [1908]. The amount of offset shown on this plot (1.5 m) has been used in several studies as being representative of the amount of fault offset along this segment of the San Andreas fault in 1906. Our historical research shows that different observers reported different amounts of fault offset in the tunnel and that the 1.5 m given on the plot is not a surveyed measurement. In addition, the plot of the tunnel has been interpreted in several previous studies as evidence of a broad (1.5 km) zone of faulting beneath Summit Ridge. Our analysis shows that this plot need not indicate a broad zone of deformation. Our historical research and modeling of the tunnel measurements indicate that faulting was confined to a zone less than 400 m wide and that 60–85% of the coseismic slip occurred across a single fault plane. There is no evidence for offset across a second shear zone beneath Summit Ridge in 1906. This implies that surface fractures reported on Summit Ridge in 1906 were not associated with significant deformation of the tunnel, implying that they were shallow, surficial features. By analogy, the very similar fractures that occurred on Summit Ridge in 1989 were also probably the result of shallow gravitational, rather than deep-seated tectonic, processes. Our modeling also indicates that total coseismic, near-surface slip across the San Andreas fault zone in the Wrights tunnel in 1906 was at least 1.7–1.8 m.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Permeability and porosity of the Illinois UPH 3 drillhole granite and a comparison with other deep drillhole rocks

Permeability, porosity, and volumetric strain measurements were conducted on granite cores obtained at depths of 0.7 to 1.6 km from the Illinois UPH 3 drillhole at effective confining pressures from 5 to 100 MPa. Initial permeabilities were in the range of 10-17 to 10-19 m2 and dropped rapidly with applied pressure to values between 10-20 and 10-24 m2 at 100 MPa, typical of other deep granite core samples. These values are several decades lower than equivalent weathered surface granites at comparable effective confining pressures, where weathering products in cracks and pores inhibit crack closure with applied pressure. Permeabilities of the Illinois cores were inversely related to sample depth, suggesting that stress relief and thermal microfractures induced during core retrieval dominated the fluid flow. Thus these samples provide an upper bound on in situ matrix permeability values. A comparison of core permeability from UPH 3 and other deep drillholes shows that stress relief damage can often dominate laboratory permeability measurements. We conclude that it may be difficult to make meaningful estimates of in situ permeability based on either borehole samples (possible damage during retrieval) or surface-derived analogs (altered by weathering). Volumetric strain determined from porosity measurements was compared with differential strain analysis (DSA) data reported by other investigators on samples from the same depths in the drillhole. Our strain measurements (0.002 to 0.005 at 100 MPa) were nearly twice as large as the DSA values, probably because of the crack-enhancing effects of fluids present in our samples that are absent in the dry DSA cores, as well as other time-dependent deformation effects. This difference in observed strain magnitudes between the two measurement methods may be an important consideration if strain and/or porosity data from deep core samples are used in models of stress, fluid circulation, and excess fluid pressure generation in the midcrust.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Three-dimensional upper crustal velocity structure beneath San Francisco Peninsula, California

This paper presents new seismic data from, and crustal models of the San Francisco Peninsula. In much of central California the San Andreas fault juxtaposes the Cretaceous granitic Salinian terrane on its west and the Late Mesozoic/Early Tertiary Franciscan Complex on its east. On San Francisco Peninsula, however, the present-day San Andreas fault is completely within a Franciscan terrane, and the Pilarcitos fault, located southwest of the San Andreas, marks the Salinian-Franciscan boundary. This circumstance has evoked two different explanations: either the Pilarcitos is a thrust fault that has pushed Franciscan rocks over Salinian rocks or the Pilarcitos is a transform fault that has accommodated significant right-lateral slip. In an effort to better resolve the subsurface structure of the peninsula faults, we established a temporary network of 31 seismographs arrayed across the San Andreas fault and the subparallel Pilarcitos fault at ∼1–2 km spacings. These instruments were deployed during the first 6 months of 1995 and recorded local earthquakes, air gun sources set off in San Francisco Bay, and explosive sources. Travel times from these sources were used to augment earthquake arrival times recorded by the Northern California Seismic Network and were inverted for three-dimensional velocity structure. Results show lateral velocity changes at depth (∼0.5–7 km) that correlate with downward vertical projections of the surface traces of the San Andreas and Pilarcitos faults. We thus interpret the faults as high-angle to vertical features (constrained to a 70°–110° dip range). From this we conclude that the Pilarcitos fault is probably an important strike-slip fault that accommodated much of the right-lateral plate boundary strain on the peninsula prior to the initiation of the modern-day San Andreas fault in this region sometime after about 3.0 m.y. ago.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Resolution of fault slip along the 470-km-long rupture of the great 1906 San Francisco earthquake and its implications

Data from all available triangulation networks affected by the 1906 earthquake have been combined to assess the trade-off between slip resolution and its uncertainty and to construct a conservative image of coseismic slip along the rupture. Because of varying network aperture and station density, slip resolution is very uneven. Although slip is determined within uncertainties of ??1.0 m along 60% of the fault, constraints are poor on the remaining, mostly offshore portions of the rupture. Slip decreases from maxima of 8.6 and 7.5 m at Shelter Cove and Tomales Bay to 4.5 m near Mount Tamalpais and 2.7 m at Loma Prieta. The geodetically derived slip distribution is in poor agreement with estimates based on analysis of S wave seismograms, probably because these waves register only 20-30% of the total seismic moment obtained from longer-period surface waves. Consideration of a range of fault geometries for 1906 slip near Loma Prieta indicates right-lateral motions lie between 2.3 and 3.1 m. These values are considerably greater than the 1.5 m of measured surface slip on which several assessments of high earthquake hazard for this fault segment were based. This factor, along with the absence of 1989 slippage where 1906 surface slip was used to make the forecasts, casts doubt on some claims of success in predicting the 1989 M = 6.9 Loma Prieta earthquake.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Empirical Green's function analysis: Taking the next step

An extension of the empirical Green's function (EGF) method is presented that involves determination of source parameters using standard EGF deconvolution, followed by inversion for a common attenuation parameter for a set of colocated events. Recordings of three or more colocated events can thus be used to constrain a single path attenuation estimate. I apply this method to recordings from the 1995–1996 Ridgecrest, California, earthquake sequence; I analyze four clusters consisting of 13 total events with magnitudes between 2.6 and 4.9. I first obtain corner frequencies, which are used to infer Brune stress drop estimates. I obtain stress drop values of 0.3–53 MPa (with all but one between 0.3 and 11 MPa), with no resolved increase of stress drop with moment. With the corner frequencies constrained, the inferred attenuation parameters are very consistent; they imply an average shear wave quality factor of approximately 20–25 for alluvial sediments within the Indian Wells Valley. Although the resultant spectral fitting (using corner frequency and κ) is good, the residuals are consistent among the clusters analyzed. Their spectral shape is similar to the the theoretical one-dimensional response of a layered low-velocity structure in the valley (an absolute site response cannot be determined by this method, because of an ambiguity between absolute response and source spectral amplitudes). I show that even this subtle site response can significantly bias estimates of corner frequency and κ, if it is ignored in an inversion for only source and path effects. The multiple-EGF method presented in this paper is analogous to a joint inversion for source, path, and site effects; the use of colocated sets of earthquakes appears to offer significant advantages in improving resolution of all three estimates, especially if data are from a single site or sites with similar site response.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

An integrated model for the tectonic development of the frontal Brooks Range and Colville Basin 250 km west of the Trans-Alaska Crustal Transect

We present a kinematic model for the sequence of deformation and sedimentation in the frontal Brooks Range and adjacent Colville Basin in the Etivluk River region, 250 km west of the Trans-Alaska Crustal Transect (TACT). The model is based on a tectonic subsidence analysis of the foreland basin, combined with structural, stratigraphic, and thermal studies of the northern edge of the Brooks Range thrust belt. We interpret six discrete tectonic events that led to the present-day configuration of the thrust belt in this area: (1) emplacement of ophiolitic allochthons over the distal continental margin rocks in Valanginian time, hundreds of kilometers south of this study, (2) Hauterivian uplift of the Barrow Arch rift margin, affecting the northern part of the Colville Basin, (3) Barremian contraction involving emplacement of distal continental margin and ophiolitic allochthons onto the Endicott Mountains allochthon and creation of a southward dipping flexural basin on the North Slope autochthon, (4) mid-Cretaceous exhumation of imbricated rocks in the Brooks Range during northward propagation of the thrust front into the foreland, (5) minor thrusting in Late Cretaceous-Paleocene in the northern foreland to the northern limit of contractional structures, and (6) regional exhumation of the orogen and the foreland in Paleocene-Eocene time. This sequence of deformation agrees well with a simple model of a forward propagating thrust system.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Tectonic evolution of the central Brooks Range mountain front: Evidence from the Atigun Gorge region

Atigun Gorge, at the northern front of the eastern Endicott Mountains, contains well-exposed rocks of the upper part of the Endicott Mountains allochthon and rocks of the structurally higher Picnic Creek or Ipnavik River allochthon. These allochthons contain rocks as young as Early Cretaceous (Valanginian) and are separated by a nearly vertical fault zone that contains exotic blocks of Triassic and Jurassic chert and silicified mudstone. Siliceous rocks of this type are not present in the Endicott Mountains allochthon but are characteristic of the Picnic Creek, Ipnavik River, and some of the other allochthons that structurally overlie the Endicott Mountains allochthon in the central and western Brooks Range. These exotic blocks, therefore indicate that structurally higher rocks of either the Picnic Creek or Ipnavik River allochthon were emplaced during the Early Cretaceous and are preserved along the northern flank of the eastern Endicott Mountains. The deformed thickness of this higher allochthon in the subsurface north of the mountains is unknown but probably exceeds 2 kilometers. Similar relations are mapped east of Atigun Gorge in an area of structural transition from the eastern Endicott Mountains into the northern Philip Smith Mountains, which are formed by the parautochthonous North Slope stratigraphic assemblage. The allochthonous rocks at the mountain front are regionally unconformably overlain by proximal Lower Cretaceous (Albian) foredeep conglomerate at the southern flank of the Colville basin, but at Atigun Gorge, the base of these deposits is interpreted as a possible back thrust at a triangle zone. Conglomerate clasts in the foredeep deposits are dominantly chert, mafic igneous rock, and other lithologies characteristic of the Picnic Creek and Ipnavik River allochthons and scattered clasts from the Endicott Mountains allochthon. The conglomerates show that the chert-rich allochthonous rocks and the Endicott Mountains allochthon were emplaced in the north-central Brooks Range by large-scale crustal shortening (>300 km) between the Valanginian and Albian (∼135 to ∼112 Ma). This orogenic event significantly postdates early stages of Brooks Range orogeny but predates later stages of orogeny documented by stratigraphic and apatite fission-track data. These relations reduce the magnitude of shortening inferred at the triangle zone at the Brooks Range mountain front. The outcrop data suggest that some of the strata preserved at a structurally low level north of the mountain front and visible in the seismic data of the Trans-Alaska Crustal Transect (TACT) may consist of clastic sedimentary rocks of the structurally higher Picnic Creek or Ipnavik River allochthon.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Crustal implications of bedrock geology along the Trans-Alaska Crustal Transect (TACT) in the Brooks Range, northern Alaska

Geologic mapping of the Trans-Alaska Crustal Transect (TACT) project along the Dalton Highway in northern Alaska indicates that the Endicott Mountains allochthon and the Hammond terrane compose a combined allochthon that was thrust northward at least 90 km in the Early Cretaceous. The basal thrust of the combined allochthon climbs up section in the hanging wall from a ductile shear zone in the south through lower Paleozoic rocks of the Hammond terrane and into Upper Devonian rocks of the Endicott Mountains allochthon at the Mount Doonerak antiform, culminating in Early Cretaceous shale in the northern foothills of the Brooks Range. Footwall rocks north of the Mount Doonerak antiform are everywhere parautochthonous Permian and Triassic shale of the North Slope terrane rather than Jurassic and Lower Cretaceous strata of the Colville Basin as shown in most other tectonic models of the central Brooks Range. Stratigraphic and structural relations suggest that this thrust was the basal detachment for Early Cretaceous deformation. Younger structures, such as the Tertiary Mount Doonerak antiform, deform the Early Cretaceous structures and are cored by thrusts that root at a depth of about 10 to 30 km along a deeper detachment than the Early Cretaceous detachment. The Brooks Range, therefore, exposes (1) an Early Cretaceous thin-skinned deformational belt developed during arc-continent collision and (2) a mainly Tertiary thick-skinned orogen that is probabty the northward continuation of the Rocky Mountains orogenic belt. A down-to-the-south zone of both ductile and brittle normal faulting along the southern margin of the Brooks Range probably formed in the mid-Cretaceous by extensional exhumation of the Early Cretaceous contractional deformation.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Composition of the crust in the Grenville and Appalachian Provinces of North America inferred from VP/VS ratios

We use the ratios between P and S wave velocities (V P /V S ), derived from seismic refraction data, to infer the composition of the crust in the Grenville and the Appalachian Provinces of North America. The crust exhibits V P /V S increasing with depth from 1.64 to 1.84; there is a clear distinction between the Grenville Province (average V P /V S =1.81) and the Appalachian Province (average V P /V S =1.73) which persists at all depths. The boundary between these provinces is east dipping extending for 100 km east of the Champlain thrust. In the Appalachian Province the increase in V P /V S ratios with depth from 1.67 to 1.74±0.02 may reflect a normal decrease of silica content in the continental crust. In the Grenville Province beneath the Central Granulite Terrane, an anomalous V P /V S ratio of 1.82±0.02 is observed extending to a depth of 10 km; this correlates with the abundance of Ca‐plagioclase in the Marcy Anorthosite. At greater depth (15–20 km), where seismic lamination and high electrical conductivity is observed, V P /V S is 1.84±0.02 and correlates with the Tahawus Complex, a layered mafic intrusion. Within the 25‐km‐thick lower crust of the Grenville Province the V P /V S is 1.84±0.02 and P‐velocity is 7.0±0.1 km/s, which are typical for plagioclase‐bearing rocks (gabbro‐norite). The high V P /V S ratio in the Grenville Province has not been reported in crust of any other age. Since the Grenville Province contains 75% of the world's known anorthosites, high V P /V S ratio is related to high plagioclase. We suggest that the composition of the Grenville lower crust was significantly modified by the emplacement of the anorthosites in the mid‐Proterozoic.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS

Continuous Global Positioning System (GPS) measurements at Long Valley Caldera, an active volcanic region in east central California, have been made on the south side of the resurgent dome since early 1993. A site on the north side of the dome was added in late 1994. Special adaptations for autonomous operation in remote regions and enhanced vertical precision were made. The data record ongoing volcanic deformation consistent with uplift and expansion of the surface above a shallow magma chamber. Measurement precisions (1 standard error) for “absolute” position coordinates, i.e., relative to a global reference frame, are 3–4 mm (north), 5–6 mm (east), and 10–12 mm (vertical) using 24 hour solutions. Corresponding velocity uncertainties for a 12 month period are about 2 mm/yr in the horizontal components and 3–4 mm/yr in the vertical component. High precision can also be achieved for relative position coordinates on short (less than 10 km) baselines using broadcast ephemerides and observing times as short as 3 hours, even when data are processed rapidly on site. Comparison of baseline length changes across the resurgent dome between the two GPS sites and corresponding two-color electronic distance measurements indicates similar extension rates within error (∼2 mm/yr) once we account for a random walk noise component in both systems that may reflect spurious monument motion. Both data sets suggest a pause in deformation for a 3.5 month period in mid-1995, when the extension rate across the dome decreased essentially to zero. Three dimensional positioning data from the two GPS stations suggest a depth (5.8±1.6 km) and location (west side of the resurgent dome) of a major inflation center, in agreement with other geodetic techniques, near the top of a magma chamber inferred from seismic data. GPS systems similar to those installed at Long Valley can provide a practical method for near real-time monitoring and hazard assessment on many active volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Färoe-Iceland Ridge Experiment: 1. Crustal structure of northeastern Iceland

Results from the Färoe-Iceland Ridge Experiment (FIRE) constrain the crustal thickness as 19 km under the Northern Volcanic Zone of Iceland and 35 km under older Tertiary areas of northeastern Iceland. The Moho is defined by strong P wave and S wave reflections. Synthetic seismogram modeling of the Moho reflection indicates mantle velocities of at least 8.0 km/s beneath the Tertiary areas of northeastern Iceland and at least 7.9 km/s beneath the neovolcanic zone. Crustal diving rays resolve the structure of the upper and lower crust. Surface P wave velocities are 1.1–4.0 km/s in Quaternary rocks and are rather higher, 4.4–4.7 km/s, in the Tertiary basalts that outcrop elsewhere. The highest crustal P wave velocities observed directly from diving rays are 7.1 km/s, from rays that turn at 24 km depth. Velocities of 7.35 km/s at the base of the crust are inferred from extrapolation of the lower crustal velocity gradient (0.024 s −1 ). A Poisson's ratio of approximately 0.27, equivalent to an S wave to P wave travel time ratio of 1.78, is measured throughout the crust east of the neovolcanic zone. The Poisson's ratio and the steep Moho topography (in places up to 30° from the horizontal) indicate that the entire crust outside the neovolcanic zone is cool (<800°C). Gravity data are well matched by a velocity/density conversion of our seismic crustal model and indicate a region of low mantle density beneath the neovolcanic zone, believed to be due to elevated mantle temperatures. The crustal thickness in the neovolcanic zone is consistent with geochemical estimates of the melt generation, placing constraints on the flow within the Iceland mantle plume.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Multispectral thermal infrared mapping of sulfur dioxide plumes: A case study from the East Rift Zone of Kilauea Volcano, Hawaii

The synoptic perspective and rapid mode of data acquisition provided by remote sensing are well suited for the study of volcanic SO 2 plumes. In this paper we describe a plume-mapping procedure that is based on image data acquired with NASA's airborne thermal infrared multispectral scanner (TIMS) and apply the procedure to TIMS data collected over the East Rift Zone of Kilauea Volcano, Hawaii, on September 30, 1988. These image data covered the Pu‘u ‘O ‘o and Kupaianaha vents and a skylight in the lava tube that was draining the Kupaianaha lava pond. Our estimate of the SO 2 emission rate from Pu‘u ‘O ‘o (17–20 kg s −1 ) is roughly twice the average of estimates derived from correlation spectrometer (COSPEC) measurements collected 10 days prior to the TIMS overflight (10 kg s −1 ). The agreement between the TIMS and COSPEC results improves when we compare SO 2 burden estimates, which are relatively independent of wind speed. We demonstrate the feasibility of mapping Pu‘u ‘O ‘o - scale SO 2 plumes from space in anticipation of the 1998 launch of the advanced spaceborne thermal emission and reflectance radiometer (ASTER).

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Deep seismic structure and tectonics of northern Alaska: Crustal-scale duplexing with deformation extending into the upper mantle

Seismic reflection and refraction and laboratory velocity data collected along a transect of northern Alaska (including the east edge of the Koyukuk basin, the Brooks Range, and the North Slope) yield a composite picture of the crustal and upper mantle structure of this Mesozoic and Cenozoic compressional orogen. The following observations are made: (1) Northern Alaska is underlain by nested tectonic wedges, most with northward vergence (i.e., with their tips pointed north). (2) High reflectivity throughout the crust above a basal decollement, which deepens southward from about 10 km depth beneath the northern front of the Brooks Range to about 30 km depth beneath the southern Brooks Range, is interpreted as structural complexity due to the presence of these tectonic wedges, or duplexes. (3) Low reflectivity throughout the crust below the decollement is interpreted as minimal deformation, which appears to involve chiefly bending of a relatively rigid plate consisting of the parautochthonous North Slope crust and a 10- to 15-km-thick section of mantle material. (4) This plate is interpreted as a southward verging tectonic wedge, with its tip in the lower crust or at the Moho beneath the southern Brooks Range. In this interpretation the middle and upper crust, or all of the crust, is detached in the southern Brooks Range by the tectonic wedge, or indentor: as a result, crust is uplifted and deformed above the wedge, and mantle is depressed and underthrust beneath this wedge. (5) Underthrusting has juxtaposed mantle of two different origins (and seismic velocities), giving rise to a prominent sub-Moho reflector.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

A numerical investigation of choked flow dynamics and its application to the triggering mechanism of long-period events at Redoubt Volcano, Alaska

We use numerical simulations of transonic flow through a crack to study the dynamics of the formation of shock waves downstream from a nozzle-like constriction inside the crack. The model solves the full set of Navier-Stokes equations in two dimensions via an explicit multifield finite difference representation. The crack walls are assumed to be perfectly rigid, and elastic coupling to the solid is not considered. The simulations demonstrate how the behavior of unsteady shock waves near the walls can produce recurring step-like pressure transients in the flow, which in turn induce resonance of the fluid-filled crack. The motion of the shock waves is governed primarily by smooth, low-amplitude pressure fluctuations at the outlet of the crack. The force induced on the walls scales with the amplitude of the shock, which is a function of the magnitude of the inlet pressure, aperture of the constriction, and thickness of the boundary layer. The applied force also scales in proportion to the spatial extent of the shock excursion, which depends on the fluctuation rate of outlet pressure. Using the source parameters of long-period (LP) events at Redoubt Volcano, Alaska, as a guide for our simulations, we infer that coupling of the shock to the walls occurs for crack inlet to outlet pressure ratios p i /p o >2.31 and that the position of the shock front becomes most sensitive to outlet pressure fluctuations for flow regimes with p i /p o >2.48. For such regimes, fluctuations of outlet pressure of up to ±0.5 MPa at rates up to 3 MPa/s are sufficient to induce pressure transients with magnitudes up to 12.5 MPa over 0.1–2.5 m of the walls within ∼0.5 s. These flow parameters may be adequate for triggering the LP events in the precursory swarm to the December 14, 1989, eruption of Redoubt. According to the flow model the recurrence rate and amplitudes of L.P events are inferred to be a manifestation of the response of a shallow hydrothermal reservoir to the sustained injection of superheated steam from a magma column roofing below this reservoir.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Post seismic deformation associated with the 1992 Mω = 7.3 Landers earthquake, southern California

Following the 1992 M ω =7.3 Landers earthquake, a linear array of 10 geodetic monuments at roughly 5‐km spacing was established across the Emerson fault segment of the Landers rupture. The array trends perpendicular to the local strike of the fault segment and extends about 30 km on either side of it. The array was surveyed by Global Positioning System 0.034, 0.048, 0.381, 1.27, 1.88, 2.60, and 3.42 years after the Landers earthquake to measure both the spatial and temporal character of the postearthquake relaxation. The temporal behavior is described roughly by a short‐term (decay time 84±23 days) exponential relaxation superimposed upon an apparently linear trend. Because the linear trend represents motions much more rapid than the observed preseismic motions, we attribute that trend to a slower (decay time greater than 5 years) postseismic relaxation, the curvature of which cannot be resolved in the short run (3.4 years) of postseismic data. About 100 mm of right‐lateral displacement and 50 mm of fault‐normal displacement accumulated across the geodetic array in the 3.4‐year interval covered by the postseismic surveys. Those displacements are attributed to postseismic, right‐lateral slip in the depth interval 10 to 30 km on the downward extension of the rupture trace. The right‐lateral slip amounted to about 1 m directly beneath the geodetic array, and the fault‐normal displacement is apparently primarily a consequence of the curvature of the rupture. These conclusions are based upon dislocation models fit to the observed deformation. However, no dislocation model was found with rms residuals as small as the expected observational error.

California↗

USGS1997

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

A test of the double-shearing model of flow for granular materials

The double-shearing model of flow attributes plastic deformation in granular materials to cooperative slip on conjugate Coulomb shears (surfaces upon which the Coulomb yield condition is satisfied). The strict formulation of the double-shearing model then requires that the slip lines in the material coincide with the Coulomb shears. Three different experiments that approximate simple shear deformation in granular media appear to be inconsistent with this strict formulation. For example, the orientation of the principal stress axes in a layer of sand driven in steady, simple shear was measured subject to the assumption that the Coulomb failure criterion was satisfied on some surfaces (orientation unspecified) within the sand layer. The orientation of the inferred principal compressive axis was then compared with the orientations predicted by the double-shearing model. The strict formulation of the model [ Spencer , 1982] predicts that the principal stress axes should rotate in a sense opposite to that inferred from the experiments. A less restrictive formulation of the double-shearing model by de Josselin de Jong [1971] does not completely specify the solution but does prescribe limits on the possible orientations of the principal stress axes. The orientations of the principal compression axis inferred from the experiments are probably within those limits. An elastoplastic formulation of the double-shearing model [ de Josselin de Jong , 1988] is reasonably consistent with the experiments, although quantitative agreement was not attained. Thus we conclude that the double-shearing model may be a viable law to describe deformation of granular materials, but the macroscopic slip surfaces will not in general coincide with the Coulomb shears.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Crustal structure along the west flank of the Cascades, western Washington

Knowledge of the crustal structure of the Washington Cascades and adjacent Puget Lowland is important to both earthquake hazards studies and geologic studies of the evolution of this tectonically active region. We present a model for crustal velocity structure derived from analysis of seismic refraction/wide-angle reflection data collected in 1991 in western Washington. The 280-km-long north-south transect skirts the west flank of the Cascades as it crosses three tectonic provinces including the Northwest Cascades Thrust System (NWCS), the Puget Lowland, and the volcanic arc of the southern Cascades. Within the NWCS, upper crustal velocities range from 4.2 to 5.7 km s-1 and are consistent with the presence of a diverse suite of Mesozoic and Paleozoic metasediments and metavolcanics. In the upper 2-3 km of the Puget Lowland velocities drop to 1.7-3.5 km s-1 and reflect the occurrence of Oligocene to recent sediments within the basin. In the southern Washington Cascades, upper crustal velocities range from 4.0 to 5.5 km s-1 and are consistent with a large volume of Tertiary sediments and volcanics. A sharp change in velocity gradient at 5-10 km marks the division between the upper and middle crust. From approximately 10 to 35 km depth the velocity field is characterized by a velocity increase from ???6.0 to 7.2 km s-1. These high velocities do not support the presence of marine sedimentary rocks at depths of 10-20 km beneath the Cascades as previously proposed on the basis of magnetotelluric data. Crustal thickness ranges from 42 to 47 km along the profile. The lowermost crust consists of a 2 to 8-km-thick transitional layer with velocities of 7.3-7.4 km s-1. The upper mantle velocity appears to be an unusually low 7.6-7.8 km s-1. When compared to velocity models from other regions, this model most closely resembles those found in active continental arcs. Distinct seismicity patterns can be associated with individual tectonic provinces along the seismic transect. In the NWCS and Puget Lowland, most of the seismicity occurs below the base of the upper crust as defined by a seismic boundary at 5-10 km depth and continues to 20-30 km depth. The region of transition between the NWCS and the Puget Lowland appears as a gap in seismicity with notably less seismic activity north of the boundary between the two. Earthquakes within the Cascades are generally shallower (0-20 km) and are dominated by events associated with the Rainier Seismic Zone.

Washington↗