USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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 1,333 records · Page 74Linked to original sources
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↗
USGS1996

Using nonlinear forecasting to learn the magnitude and phasing of time-varying sediment suspension in the surf zone

The time-dependent response of sediment suspension to flow velocity was explored by modeling field measurements collected in the surf zone during a large storm. Linear and nonlinear models were created and tested using flow velocity as input and suspended-sediment concentration as output. A sequence of past velocities (velocity history), as well as velocity from the same instant as the suspended-sediment concentration, was used as input; this velocity history length was allowed to vary. The models also allowed for a lag between input (instantaneous velocity or end of velocity sequence) and output (suspended-sediment concentration). Predictions of concentration from instantaneous velocity or instantaneous velocity raised to a power (up to 8) using linear models were poor (correlation coefficients between predicted and observed concentrations were less than 0.10). Allowing a lag between velocity and concentration improved linear models (correlation coefficient of 0.30), with optimum lag time increasing with elevation above the seabed (from 1.5 s at 13 cm to 8.5 s at 60 cm). These lags are largely due to the time for an observed flow event to effect the bed and mix sediment upward. Using a velocity history further improved linear models (correlation coefficient of 0.43). The best linear model used 12.5 s of velocity history (approximately one wave period) to predict concentration. Nonlinear models gave better predictions than linear models, and, as with linear models, nonlinear models using a velocity history performed better than models using only instantaneous velocity as input. Including a lag time between the velocity and concentration also improved the predictions. The best model (correlation coefficient of 0.58) used 3 s (approximately a quarter wave period) of the cross-shore velocity squared, starting at 4.5 s before the observed concentration, to predict concentration. Using a velocity history increases the performance of the models by specifying a more complete description of the dynamical forcing of the flow (including accelerations and wave phase and shape) responsible for sediment suspension. Incorporating such a velocity history and a lag time into the formulation of the forcing for time-dependent models for sediment suspension in the surf zone will greatly increase our ability to predict suspended-sediment transport.

Journal of Geophysical Research C: Oceans↗
USGS1996

Factors affecting suspended-solids concentrations in South San Francisco Bay, California

Measurements of suspended-solids concentration (SSC) were made at two depths at three sites in South San Francisco Bay (South Bay) to determine the factors that affect SSC. Twenty-eight segments of reliable and continuous SSC time series data longer than 14 days were collected from late 1991 or 1992 through September 1993. Spectral analysis and singular spectrum analysis were used to relate these data segments to time series of several potential forcing factors, including diurnal and semidiurnal tides, the spring-neap tidal cycle, wind shear, freshwater runoff, and longitudinal density differences. SSC is greatest during summer, when a landward wind shear is applied to South Bay by the afternoon sea breeze. About one half the variance of SSC is caused by the spring-neap cycle, and SSC lags the spring-neap cycle by about 2 days. Relatively short duration of slack water limits the duration of deposition of suspended solids and consolidation of newly deposited bed sediment during the tidal cycle, so suspended solids accumulate in the water column as a spring tide is approached and slowly deposit as a neap tide is approached. Perturbations in SSC caused by wind and local runoff from winter storms during the study period were usually much smaller than SSC variations caused by the spring-neap cycle. Variations of SSC at the study sites at tidal timescales are tidally forced, and nonlinear physical processes are significant. Advective transport dominates during spring tides, when water with higher SSC due to wind wave resuspension is advected to the main channel from shallow water, but, during neap tides, advective transport is less significant. The findings of this and other studies indicate that the tidally averaged transport of suspended solids responds to seasonal variations of wind shear in South Bay.

Journal of Geophysical Research C: Oceans↗
USGS1996

Subtidal circulation patterns in a shallow, highly stratified estuary: Mobile Bay, Alabama

Mobile Bay is a wide (25-50 km), shallow (3 m), highly stratified estuary on the Gulf coast of the United States. In May 1991 a series of instruments that measure near-surface and near-bed current, temperature, salinity, and middepth pressure were deployed for a year-long study of the bay. A full set of measurements were obtained at one site in the lower bay; all but current measurements were obtained at a midbay site. These observations show that the subtidal currents in the lower bay are highly sheared, despite the shallow depth of the estuary. The sheared flow patterns are partly caused by differential forcing from wind stress and river discharge. Two wind-driven flow patterns actually exist in lower Mobile Bay. A barotropic response develops when the difference between near-surface and near-bottom salinity is less than 5 parts per thousand. For stronger salinity gradients the wind-driven currents are larger and the response resembles a baroclinic flow pattern. Currents driven by river flows are sheared and also have a nonlinear response pattern. Only near-surface currents are driven seaward by discharges below 3000 m3/s. At higher discharge rates, surface current variability uncouples from the river flow and the increased discharge rates drive near-bed current seaward. This change in the river-forced flow pattern may be associated with a hydraulic jump in the mouth of the estuary. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS1996

Origin of Bermuda's clay-rich Quaternary paleosols and their paleoclimatic significance

Red clayey paleosols that are chiefly the product of aerosolic dust deposition are interbedded in the Quaternary carbonate formations of the Bermuda oceanic island system. These paleosols provide a basis for reconstructing Quaternary atmospheric circulation patterns in the northwestern Atlantic. Geochemical analyses were performed on representative paleosol samples to identify their parent dust source. Fine-grained fractions were analyzed by energy-dispersive X ray fluorescence to determine trace element (Zr, Y, La, Ti, and Nb) concentrations and to derive geochemical signatures based on immobile element ratios. These ratios were compared with geochemical signatures determined for three possible sources of airborne dust: (1) Great Plains loess, (2) Mississippi River Valley loess, and (3) Saharan dust. The Zr/Y and Zr/La ratios provided the clearest distinction between the hypothesized dust sources. The low ratios in the paleosol B horizons most closely resemble Saharan dust in the <2-μm size class fraction. Contributions from the two North American loessial source areas could not be clearly detected. Thus Bermuda paleosols have a predominantly Saharan aerosolic dust signature. Saharan dust deposition on Bermuda during successive Quaternary glacial periods is consistent with patterns of general circulation models, which indicate that during glacial maxima the northeast summer trade winds were stronger than at present and reached latitudes higher than 30°N despite lower-than-present sea surface temperatures in the North Atlantic.

Journal of Geophysical Research D: Atmospheres↗
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↗
USGS1997

Permian deposition in the north central Brooks Range, Alaska: Constraints for tectonic reconstructions

Two opposing tectonic models have been offered to explain the regional structural relations in the north central Brooks Range fold-thrust belt of northern Alaska. The first suggests that rocks of the northern Endicott Mountains were thrust from south to north over the area of the present Mount Doonerak high and are therefore highly allochthonous. The second implies that the rocks of the northern Endicott Mountains were deposited in a basin that lay north of the Mount Doonerak high and later were thrust a short distance southward onto the northern flank of the high and are thus parautochthonous. To provide stratigraphic constraints for these models, this study examines Permian facies of the north central Brooks Range. Permian rocks in the north central Brooks Range comprise a thin (40 to 160 m thick), fining-upward succession of clastic, storm-influenced shelf deposits. When the rocks of the northern Endicott Mountains are restored south of the Mount Doonerak area, a minimum distance of 80 km, the Permian deposits grade systematically from distal facies (Siksikpuk Formation) in the southwest to proximal facies (Echooka Formation) in the northeast. Facies trends in the reconstructed Permian basin include, from southwest to northeast, (1) an increase in carbonate content and corresponding decrease in silica content, (2) a general darkening and thickening of shaley intervals, (3) an increase in proximal features of storm beds, including coarser, thicker, more abundant, and more closely spaced beds, and (4) an increase in abundance and diversity of the faunal assemblage with a corresponding decrease in age. These stratigraphic relations imply that rocks of the northern Endicott Mountains are allochthonous and structurally overlie a proximal stratigraphie succession similar to that exposed in the Mount Doonerak area and northeastern Brooks Range.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Source and path effects in the wave fields of tremor and explosions at Stromboli Volcano, Italy

The wave fields generated by Strombolian activity are investigated using data from small-aperture seismic arrays deployed on the north flank of Stromboli and data from seismic and pressure transducers set up near the summit crater. Measurements of slowness and azimuth as a function of time clearly indicate that the sources of tremor and explosions are located beneath the summit crater at depths shallower than 200 m with occasional bursts of energy originating from sources extending to a depth of 3 km. Slowness, azimuth, and particle motion measurements reveal a complex composition of body and surface waves associated with topography, structure, and source properties. Body waves originating at depths shallower than 200 m dominate the wave field at frequencies of 0.5–2.5 Hz, and surface waves generated by the surficial part of the source and by scattering sources distributed around the island dominate at frequencies above 2.5 Hz. The records of tremor and explosions are both dominated by SH motion. Far-field records from explosions start with radial motion, and near-field records from those events show dominantly horizontal motion and often start with a low-frequency (1–2 Hz) precursor characterized by elliptical particle motion, followed within a few seconds by a high-frequency radial phase (1–10 Hz) accompanying the eruption of pyroclastics. The dominant component of the near- and far-field particle motions from explosions, and the timing of air and body wave phases observed in the near field, are consistent with a gas-piston mechanism operating on a shallow (<200 m deep), vertical crack-like conduit. Models of a degassing fluid column suggest that noise emissions originating in the collective oscillations of bubbles ascending in the magma conduit may provide an adequate self-excitation mechanism for sustained tremor generation at Stromboli.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Strengths of serpentinite gouges at elevated temperatures

Serpentinite has been proposed as a cause of both low strength and aseismic creep of fault zones. To test these hypotheses, we have measured the strength of chrysotile-, lizardite-, and antigorite-rich serpentinite gouges under hydrothermal conditions, with emphasis on chrysotile, which has thus far received little attention. At 25°C, the coefficient of friction, μ, of chrysotile gouge is roughly 0.2, whereas the lizardite- and antigorite-rich gouges are at least twice as strong. The very low room temperature strength of chrysotile is a consequence of its unusually high adsorbed water content. When the adsorbed water is removed, chrysotile is as strong as pure antigorite gouge at room temperature. Heating to ∼200°C causes the frictional strengths of all three gouges to increase. Limited data suggest that different polytypes of a given serpentine mineral have similar strengths; thus deformation-induced changes in polytype should not affect fault strength. At 25°C, the chrysotile gouge has a transition from velocity strengthening at low velocities to velocity weakening at high velocities, consistent with previous studies. At temperatures up to ∼200°C, however, chrysotile strength is essentially independent of velocity at low velocities. Overall, chrysotile has a restricted range of velocity-strengthening behavior that migrates to higher velocities with increasing temperature. Less information on velocity dependence is available for the lizardite and antigorite gouges, but their behavior is consistent with that outlined for chrysotile. The marked changes in velocity dependence and strength of chrysotile with heating underscore the hazards of using room temperature data to predict fault behavior at depth. The velocity behavior at elevated temperatures does not rule out serpentinite as a cause of aseismic slip, but in the presence of a hydrostatic fluid pressure gradient, all varieties of serpentine are too strong to explain the apparent weakness of faults such as the San Andreas.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

A dynamic balance between magma supply and eruption rate at Kilauea volcano, Hawaii

The dynamic balance between magma supply and vent output at Kilauea volcano is used to estimate both the volume of magma stored within Kilauea volcano and its magma supply rate. Throughout most of 1991 a linear decline in volume flux from the Kupaianaha vent on Kilauea's east rift zone was associated with a parabolic variation in the elevation of Kilauea's summit as vent output initially exceeded then lagged behind the magma supply to the volcano. The correspondence between summit elevation and tilt established with over 30 years of data provided daily estimates of summit elevation in terms of summit tilt. The minimum in the parabolic variation in summit tilt and elevation (or zero elevation change) occurs when the magma supply to the reservoir from below the volcano equals the magma output from the reservoir to the surface, so that the magma supply rate is given by vent flux on that day. The measurements of vent flux and tilt establish that the magma supply rate to Kilauea volcano on June 19, 1991, was 217,000±10,000 m 3 /d (or 0.079±0.004 km 3 /yr). This is close to the average eruptive rate of 0.08 km 3 /yr between 1958 and 1984. In addition, the predictable response of summit elevation and tilt to each east rift zone eruption near Puu Oo since 1983 shows that summit deformation is also a measure of magma reservoir pressure. Given this, the correlation between the elevation of the Puu Oo lava lake (4 km uprift of Kupaianaha and 18 km from the summit) and summit tilt provides an estimate for magma pressure changes corresponding to summit tilt changes. The ratio of the change in volume to the change in reservoir pressure ( dV / dP ) during vent activity may be determined by dividing the ratio of volume erupted to change in summit tilt ( dV / d tilt ) by the ratio of pressure change to change in summit tilt ( dP / d tilt ). This measure of dV / dP , when combined with laboratory measurements of the bulk modulus of tholeitic melt, provides an estimate of 240±50 km 3 for the volume of Kilauea's magma reservoir. This estimate is much larger than traditional estimates but consistent with seismic tomographic imaging and geophysical modeling of Kilauea's magma system.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Quantifying Anderson's fault types

Anderson [1905] explained three basic types of faulting (normal, strike-slip, and reverse) in terms of the shape of the causative stress tensor and its orientation relative to the Earth's surface. Quantitative parameters can be defined which contain information about both shape and orientation [ Célérier , 1995], thereby offering a way to distinguish fault-type domains on plots of regional stress fields and to quantify, for example, the degree of normal-faulting tendencies within strike-slip domains. This paper offers a geometrically motivated generalization of Angelier's [1979, 1984, 1990] shape parameters ϕ and ψ to new quantities named A ϕ and A ψ . In their simple forms, A ϕ varies from 0 to 1 for normal, 1 to 2 for strike-slip, and 2 to 3 for reverse faulting, and Aψ/ranges from 0° to 60°, 60° to 120°, and 120° to 180°, respectively. After scaling, A ϕ and A ψ agree to within 2% (or 1°), a difference of little practical significance, although A ψ has smoother analytical properties. A formulation distinguishing horizontal axes as well as the vertical axis is also possible, yielding an A ϕ ranging from −3 to +3 and A ψ from −180° to +180°. The geometrically motivated derivation in three-dimensional stress space presented here may aid intuition and offers a natural link with traditional ways of plotting yield and failure criteria. Examples are given, based on models of Bird [1996] and Bird and Kong [1994], of the use of Anderson fault parameters A ϕ and A ψ for visualizing tectonic regimes defined by regional stress fields.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Seismic-geodynamic constraints on three-dimensional structure, vertical flow, and heat transfer in the mantle

Joint inversions of seismic and geodynamic data are carried out in which we simultaneously constrain global-scale seismic heterogeneity in the mantle as well as the amplitude of vertical mantle flow across the 670 km seismic discontinuity. These inversions reveal the existence of a family of three-dimensional (3-D) mantle models that satisfy the data while at the same time yielding predictions of layered mantle flow. The new 3-D mantle models we obtain demonstrate that the buoyancy forces due to the undulations of the 670 km phase-change boundary strongly inhibit the vertical flow between the upper and lower mantle. The strong stabilizing effect of the 670 km topography also has an important impact on the predicted dynamic topography of the Earth's solid surface and on the surface gravity anomalies. The new 3-D models that predict strongly or partially layered mantle flow provide essentially identical fits to the global seismic data as previous models that have, until now, predicted only whole-mantle flow. The convective vertical transport of heat across the mantle predicted on the basis of the new 3-D models shows that the heat flow is a minimum at 1000 km depth. This suggests the presence at this depth of a globally defined horizon across which the pattern of lateral heterogeneity changes rapidly.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Depth dependence of earthquake frequency-magnitude distributions in California: Implications for rupture initiation

Statistics of earthquakes in California show linear frequency-magnitude relationships in the range of M 2.0 to M 5.5 for various data sets. Assuming Gutenberg-Richter distributions, there is a systematic decrease in b value with increasing depth of earthquakes. We find consistent results for various data sets from northern and southern California that both include and exclude the larger aftershock sequences. We suggest that at shallow depth (∼0 to 6 km) conditions with more heterogeneous material properties and lower lithospheric stress prevail. Rupture initiations are more likely to stop before growing into large earthquakes, producing relatively more smaller earthquakes and consequently higher b values. These ideas help to explain the depth-dependent observations of foreshocks in the western United States. The higher occurrence rate of foreshocks preceding shallow earthquakes can be interpreted in terms of rupture initiations that are stopped before growing into the mainshock. At greater depth (9–15 km), any rupture initiation is more likely to continue growing into a larger event, so there are fewer foreshocks. If one assumes that frequency-magnitude statistics can be used to estimate probabilities of a small rupture initiation growing into a larger earthquake, then a small ( M 2) rupture initiation at 9 to 12 km depth is 18 times more likely to grow into a M 5.5 or larger event, compared to the same small rupture initiation at 0 to 3 km.

Journal of Geophysical Research B: Solid Earth↗