Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geophysical Journal”

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 145 records · Page 8Linked to original sources
USGS2010

Reflection—refraction of general P- and type-I S-waves in elastic and anelastic solids

The reflection and refraction of general (homogeneous or inhomogeneous) plane P and type-I S ( SV ) body waves incident on plane boundaries are considered for general linear viscoelastic solids. Reflection—refraction laws, physical characteristics of the waves, and the nature of critical angles are examined in detail at welded boundaries and a free surface. General visco-elasticity with no low-loss approximations predicts that contrasts in intrinsic absorption at boundaries give rise to inhomogeneous reflected and refracted waves with elliptical particle motions, velocities and maximum attenuations that vary with frequency and angle of incidence, energy propagation at speeds and directions different from phase propagation, phase propagation that in general is parallel to the boundary for at most one angle of incidence, and reflection—transmission coefficients dependent on energy flow due to wave interaction. None of these physical characteristics are predicted for waves incident on boundaries that respond instantaneously.

Geophysical Journal International↗
USGS1982

International geomagnetic reference field 1980: a report by IAGA Division I working group.

Describes the recommendations of the working group, which suggested additions to IGRF because of the cumulative effect of the inevitable uncertainties in the secular variation models which had led to unacceptable inaccuracies in the IGRF by the late 1970's. The recommendations were accepted by the International Association of Geomagnetism and Aeronomy on August 15, 1981 at the 4th Scientific Assembly, Edinburgh. An extended table sets out spherical harmonic coefficients of the IGRF 1980.-R.House

Geophysical Journal of the Royal Astronomical Soci↗
USGS1983

A description of the external and internal quiet daily variation currents at North American locations for a quiet-Sun year.

An order 4, degree 12 spherical harmonic analysis of the smoothed quiet geomagnetic daily variations was used to separate the external and internal geomagnetic Sq field at North American locations for the quiet-Sun year, 1965. These fields were represented by a month-by-month display of equivalent current vortex systems with dominant, pre-noon foci. The focus reached 40o latitude near the June solstice and about 30o latitude near the December solstice. The daily range of Sq current amplitudes was largest in late July to early August and smallest in mid-December. Semi-annual variations of Sq currents dominated only the equatorial region. Daily maxima in mid-latitudes, occurred mostly near local noon in December to February and about 1 hr before noon in June to mid-October. -Author

Geophysical Journal of the Royal Astronomical Soci↗
USGS1984

Inverse problems for torsional modes.

Considers a spherically symmetric, non-rotating Earth consisting of an isotropic, perfect elastic material where the density and the S-wave velocity may have one or two discontinuities in the upper mantle. Shows that given the velocity throughout the mantle and the crust and given the density in the lower mantle, then the freqencies of the torsional oscillations of one angular order (one torsional spectrum), determine the density in the upper mantle and in the crust uniquely. If the velocity is known only in the lower mantle, then the frequencies of the torsional oscillations of two angular orders uniquely determine both the density and the velocity in the upper mantle and in the crust. In particular, the position and size of the discontinuities in the density and velocity are uniquely determined by two torsional spectra.-Author

Geophysical Journal of the Royal Astronomical Soci↗
USGS1985

Inversion of seismic refraction data in planar dipping structure

A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time—slowness (tau— p ) wavefield. The tau— p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model. This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.

Geophysical Journal International↗
USGS1987

Coincident seismic reflection/refraction studies of the continental lithosphere: A global review

Vertical-incidence reflection profiling has identified several characteristic features of the continental lithosphere including a generally transparent upper crust, a reflective lower crust, reflections from the crust-mantle boundary, and a commonly transparent upper mantle. The underlying physical causes of these characteristic features remain poorly understood. This review summarizes additional information brought to bear on the physical properties of these characteristic crustal structures through the use of coincident wide-angle refraction profiling.

Geophysical Journal International↗
USGS1987

Crustal structure beneath exposed accreted terranes of Southern Alaska

The crustal structure beneath the exposed terranes of southern Alaska has been explored using coincident seismic refraction and reflection profiling. A wide-angle reflector at 8-9 km depth, at the base of an inferred low-velocity zone, underlies the Peninsular and Chugach terranes, appears to truncate their boundary, and may represent a horizontal decollement beneath the terranes. The crust beneath the Chugach terrane is characterized by a series of north-dipping paired layers having low and high velocities that may represent subducted slices of oceanic crust and mantle. This layered series may continue northward under the Peninsular terrane. Earthquake locations in the Wrangell Benioff zone indicate that at least the upper two low-high velocity layer pairs are tectonically inactive and that they appear to have been accreted to the base of the continental crust. The refraction data suggest that the Contact fault between two similar terranes, the Chugach and Prince William terranes, is a deeply penetrating feature that separates lower crust (deeper than 10 km) with paired dipping reflectors, from crust without such reflectors.

Alaska↗
USGS1988

Waveform modelling using locked-mode synthetic and differential seismograms: application to determination of the structure of Mexico

We have developed algorithms for modelling seismic waveforms to constrain regional Earth structure. The seismogram is represented as a sum of locked-mode travelling waves in a layered medium. This representation is convenient as it allows us to model structures with slowly varying heterogeneity and to construct differential seismograms. Describes the techniques we have implemented that enable us to compute synthetic and differential seismograms in an efficient and stable manner. The computational methods are sufficiently rapid that many modes can be included and in some cases the entire seismogram may be modified. These algorithms are applied to model a set of seismograms of southern Mexican earthquakes recorded in northern Mexico. The frequency bandwidth of these data is centred at 0.067 Hz and we demonstrate that even at these relatively high frequencies, many features of the seismogram can be successfully modelled. Our results suggest that the structure within the recording array in northern Mexico is resolvably different from that to the south. We find that the average shear velocity of the lower lithosphere of southern Mexico is very low, approximately 4.3 km s-1. If the low-velocity region is confined to the Trans Mexican Volcanic Belt, the shear velocities between 20-80 km depth are approximately 3.3 km s-1. This may be correlated with partial melt and is consistent with the active volcanism and high heat flow found in the region. -Authors

Geophysical Journal International↗
USGS1990

Brunhes chron excursion/polarity episode recorded during the late pleistocene, Albuquerque Volcanoes, New Mexico, USA

All basaltic lava flows of the Albuquerque Volcanoes (lat.: 35.2°N, long.: 253.2°E), Albuquerque-Belen Basin, New Mexico, record a short excursion/polarity episode. K-Ar isotopic age determinations (weighted average: 155 ± 47 ka) and evaluation of soil profiles on flow surfaces suggest the late Pleistocene ( circa between 250 and 80 ka) as the time of extrusion. Results from 63 sites in a minimum of eight flows yield a mean direction of D = 101.1°, I = -36.1°, α 95 = 1.2° (α 1 95 = 0.7, α 2 95 = 1.2), k = 2219 ( N = 8 flows) and a corresponding virtual geomagnetic pole position (VGP) of 354.1°E, 20.2°S, dp = 0.8°, dm = 1.4° ( A 1 95 = 0.5, A 2 95 = 1.4), and VGP angular standard deviation (ASD) = 1.8°. The unusual magnetization in flows of the Albuquerque Volcanoes is carried by fine-grained [single-domain (SD) and/or pseudo-single-domain (PSD)] low-Ti magnetite and is apparently not a function of complex sub-solidus alteration of magnetic phases. Underlying baked soils contain a relatively low coercivity magnetization with directions comprising two groups. One is roughly antipodal to the lava remanence; the other intermediate between the lava remanence and present field direction. Because all flows yield statistically indistinguishable directions, a field directional path clearly cannot be defined. At the precision level of our age determinations, the Albuquerque feature may correlate with other short polarity episodes or excursions of late Pleistocene age (e.g., Blake or Jamaica). Though characteristically limited in directional morphology, well-dated polarity episodes (i.e. sub-chrons) and excursions recorded in volcanic rocks may provide information on the frequency of significant dynamo instabilities.

New Mexico↗
USGS1991

Single-station decomposition of seismograms for subevent time histories

We have adapted an iterative least-squares inversion with positivity constraints to the problem of deconvolving the recordings of a small earthquake from the recordings of a large composite earthquake. The inversion uses an F test to minimize the number of non-zero amplitudes or subevents in the solution set, by requiring each added subevent to reduce the variance significantly. We apply this inversion to the analysis of strong motion recordings of a set of moderate-sized (4.0 ≤ M L ≤ 6.4) aftershocks of the 1976 Friuli earthquake. The earthquake pairs are selected by comparing S trigger times and polarization diagrams. The deconvolved time series are interpreted as the time history of slip events or stress release in the composite earthquake, depending on whether the recordings of the small earthquake are unfiltered or filtered to simulate the waveforms radiated by subevents. Deconvolving with the filtered recordings reduces the variance more rapidly, as a function of the number of subevents, than deconvolving with the unfiltered recordings. The deconvolved time series for three of the four earthquakes analysed indicate that the rupture processes resemble the failure of asperities where the initial subevent has the largest stress drop and the later subevents have longer durations and weaker stress drops.

Geophysical Journal International↗
USGS1993

Large-scale deformation associated with ridge subduction

Continuum models are used to investigate the large-scale deformation associated with the subduction of aseismic ridges. Formulated in the horizontal plane using thin viscous sheet theory, these models measure the horizontal transmission of stress through the arc lithosphere accompanying ridge subduction. Modelling was used to compare the Tonga arc and Louisville ridge collision with the New Hebrides arc and d'Entrecasteaux ridge collision, which have disparate arc-ridge intersection speeds but otherwise similar characteristics. Models of both systems indicate that diffuse deformation (low values of the effective stress-strain exponent n ) are required to explain the observed deformation. Deformation is somewhat insensitive to the vertically integrated strength of the arc (inversely proportional to the Argand number Ar ), but indicates that the arc lithosphere is not extremely weak ( Ar < 100). Low values of both Ar and n suggest that the thermal structure is typical of ‘cold’ or ‘normal’ arcs and that deformation is dominated by flow in the lower crust and mantle. In addition, low values of n (approaching Newtonian flow) may indicate that specific deformation mechanisms dictate deformation of the arc lithosphere. Possible mechanisms include low-stress, grain-size dependent creep, pyroxenite-controlled rheology and mechanisms associated with water weakening. Changes in the boundary conditions greatly affect deformation within island arcs. High rates of arc-ridge intersection speed (Tonga-Louisville system) yield arc-parallel tension and crustal thickening in the wake of ridge subduction. In contrast, low rates of arc-ridge intersection speed (New Hebrides-d'Entrecasteaux system) yield compressional deformation directly arcward of the collision zone and transverse strike-slip faulting adjacent to the region of compressional deformation. Localized regions of extensional deformation along the frontal part of the arc adjacent to the collision zone may contribute to the formation of re-entrants.

Geophysical Journal International↗
USGS1995

Finite-fault analysis of the 1979 March 14 Petatlan, Mexico, earthquake using teleseismic P waveforms

Vertical, teleseismic P waves recorded for the 1979 March 14 Petatlan, Mexico, earthquake were used to derive the distribution of coseismic slip using a linear finite-fault inversion scheme that solves for the amount of slip in each of a series of consecutive time windows. Data recorded by six stations of the Global Digital Seismograph Network were inverted in addition to digitized analogue long-period recordings available from nine Worldwide Standardized Seismograph Network stations. The digital data include four broad-band and short-period velocity waveforms reconstructed from the short- and long-period components. The time-window approach allows for a variable rise time on the fault and accounts for the source multiplicity evident in the recorded P waveforms. Synthetic tests conducted using the inversion method on the limited data set, however, reveal that the data are insufficient to identify the exact dislocation duration on the fault. The method is thus implemented by prescribing the fault rise time using five consecutive 1 s time windows. The coseismic slip inferred from the P waves shows a small 70 cm peak near the earthquake hypocentre and a large zone of dislocation (1.2 m maximum) further south-east. The slip pattern covers depths from 3 to 25 km and is located south-east of other recent large interplate ruptures on the Michoacan segment of the Mexican subduction zone. This result indicates that the 1979 Petatlan earthquake broke an independent, adjacent portion of the Cocos-North America plate boundary. The seismic moment of 1.5 × 10 27 dyn cm inferred from the P waves is approximately one-half the long-period moment estimated by other investigators from the observed surface waves. Although the discrepancy is within the uncertainty of the seismic-moment estimates, it may suggest the presence of a component of slow interplate motion that did not radiate significant P -wave energy.

Geophysical Journal International↗
USGS1996

The cyclic and fractal seismic series preceding an mb 4.8 earthquake on 1980 February 14 near the Virgin Islands

Seismic activity in the 10 months preceding the 1980 February 14, m b 4.8 earthquake in the Virgin Islands, reported on by Frankel in 1982, consisted of four principal cycles. Each cycle began with a relatively large event or series of closely spaced events, and the duration of the cycles progressively shortened by a factor of about 3/4. Had this regular shortening of the cycles been recognized prior to the earthquake, the time of the next episode of seismicity (the main shock) might have been closely estimated 41 days in advance. That this event could be much larger than the previous events is indicated from time-to-failure analysis of the accelerating rise in released seismic energy, using a non-linear time- and slip-predictable foreshock model. Examination of the timing of all events in the sequence shows an even higher degree of order. Rates of seismicity, measured by consecutive interevent times, when plotted on an iteration diagram of a rate versus the succeeding rate, form a triangular circulating trajectory. The trajectory becomes an ascending helix if extended in a third dimension, time. This construction reveals additional and precise relations among the time intervals between times of relatively high or relatively low rates of seismic activity, including period halving and doubling. The set of 666 time intervals between all possible pairs of the 37 recorded events appears to be a fractal; the set of time points that define the intervals has a finite, non-integer correlation dimension of 0.70. In contrast, the average correlation dimension of 50 random sequences of 37 events is significantly higher, close to 1.0. In a similar analysis, the set of distances between pairs of epicentres has a fractal correlation dimension of 1.52. Well-defined cycles, numerous precise ratios among time intervals, and a non-random temporal fractal dimension suggest that the seismic series is not a random process, but rather the product of a deterministic dynamic system.

Geophysical Journal International↗
USGS1996

Inelastic models of lithospheric stress - I. Theory and application to outer-rise plate deformation

Outer-rise stress distributions determined in the manner that mechanical engineers evaluate inelastic stress distributions within conventional materials are contrasted with those predicted using simple elastic-plate models that are frequently encountered in studies of outer-rise seismicity. This comparison indicates that the latter are inherently inappropriate for studies of intraplate earthquakes, which are a direct manifestation of lithospheric inelasticity. We demonstrate that the common practice of truncating elastically superimposed stress profiles so that they are not permitted to exceed laboratory-based estimates of lithospheric yield strength will result in an accurate characterization of lithospheric stress only under relatively restrictive circumstances. In contrast to elastic-plate models, which predict that lithospheric stress distributions depend exclusively upon the current load, inelastic plate models predict that stress distributions are also significantly influenced by the plate-loading history, and, in many cases, this influence is the dominant factor in determining the style of potential seismicity (e.g. thrust versus normal faulting). Numerous ‘intuitive’ interpretations of outer-rise earthquakes have been founded upon the implicit assumption that a unique relationship exists between a specified combination of plate curvature and in-plane force, and the resulting lithospheric stress distribution. We demonstrate that the profound influence of deformation history often invalidates such interpretations. Finally, we examine the reliability of ‘yield envelope’ representations of lithospheric strength that are constructed on the basis of empirically determined frictional sliding relationships and silicate plastic-flow laws. Although representations of this nature underestimate the strength of some major interplate faults, such as the San Andreas, they appear to represent a reliable characterization of the strength of intraplate oceanic lithosphere.

Geophysical Journal International↗
USGS1996

Inelastic models of lithospheric stress - II. Implications for outer-rise seismicity and dynamics

Outer-rise seismicity and dynamics are examined using inelastic models of lithospheric deformation, which allow a more realistic characterization of stress distributions and failure behaviour. We conclude that thrust- and normal-faulting outer-rise earthquakes represent substantially different states of stress within the oceanic lithosphere. Specifically, the normal-faulting events occur in response to downward plate bending, which establishes the ‘standard’, bending-dominated state of outer-rise stress, and the thrust-faulting events occur in response to an elevated level of in-plane compression, which develops only in response to exceptional circumstances. This interpretation accounts for the observation that normal-faulting outer-rise earthquakes occur more frequently and are more widely distributed than their thrust-faulting counterparts, an observation for which the simple bending model offers no explanation. In addition, attributing both thrust- and normal-faulting outer-rise earthquakes to plate bending implies that both classes of events should occur within relatively close lateral proximity to one another because both are allegedly a manifestation of the same bendingdominated stress distribution, whereas, in reality, this is not observed. We propose that the tendency for thrust-faulting outer-rise earthquakes to exhibit greater source depths than their normal-faulting counterparts (an observation that is frequently cited in support of the bending interpretation of the former) is merely a consequence of the fact that bending-induced tension is confined to the upper lithosphere. Our model predicts that outer-rise in-plane-force variations may promote thrust-faulting outerrise activity prior to an underthrusting interplate subduction earthquake and normalfaulting outer-rise activity following such an earthquake, but that both forms of outerrise activity are unlikely to be associated with the same subduction earthquake. A corollary implication of our model is t...

Geophysical Journal International↗
USGS1997

Seismic source study of the Racha-Dzhava (Georgia) earthquake from aftershocks and broad-band teleseismic body-wave records: An example of active nappe tectonics

The Racha-Dzhava earthquake ( Ms =7.0) that occurred on 1991 April 29 at 09:12:48.1 GMT in the southern border of the Great Caucasus is the biggest event ever recorded in the region, stronger than the Spitak earthquake ( Ms =6.9) of 1988. A field expedition to the epicentral area was organised and a temporary seismic network of 37 stations was deployed to record the aftershock activity. A very precise image of the aftershock distribution is obtained, showing an elongated cloud oriented N105°, with one branch trending N310° in the western part. The southernmost part extends over 80 km, with the depth ranging from 0 to 15 km, and dips north. The northern branch, which is about 30 km long, shows activity that ranges in depth from 5 to 15 km. The complex thrust dips northwards. A stress-tensor inversion from P -wave first-motion polarities shows a state of triaxial compression, with the major principal axis oriented roughly N-S, the minor principal axis being vertical. Body-waveform inversion of teleseismic seismograms was performed for the main shock, which can be divided into four subevents with a total rupture-time duration of 22 s. The most important part of the seismic moment was released by a gentle northerly dipping thrust. The model is consistent with the compressive tectonics of the region and is in agreement with the aftershock distribution and the stress tensor deduced from the aftershocks. The focal mechanisms of the three largest aftershocks were also inverted from body-wave records. The April 29th ( Ms =6.1) and May 5th ( Ms =5.4) aftershocks have thrust mechanisms on roughly E-W-oriented planes, similar to the main shock. Surprisingly, the June 15th ( Ms =6.2) aftershock shows a thrust fault striking N-S. This mechanism is explained by the structural control of the rupture along the east-dipping geometry of the Dzirula Massif close to the Borzhomi-Kazbeg strike-slip fault. In fact, the orientation and shape of the stress tensor produce a thrust on a N-S oriented plane. Nappe tectonics has been identified as an important feature in the Caucasus, and the source mechanism is consistent with this observation. A hidden fault is present below the nappe, and no large surface breaks were observed due to the main shock. The epicentral region is characterized by sediments that are trapped between two crystalline basements: the Dzirula Massif, which crops out south of Chiatoura, and the Caucasus Main Range north of Oni. Most, if not all, of the rupture is controlled by the thrusting of overlapping, deformed and folded sediments over the Dzirula Massif. This event is another example of blind active faults, with the distinctive feature that the fault plane dips at a gentle angle. The Racha Range is one of the surface expressions of this blind thrust, and its growth is the consequence and evidence of similar earthquakes in the past.

Geophysical Journal International↗
USGS1997

On the origin of diverse aftershock mechanisms following the 1989 Loma Prieta earthquake

We test the hypothesis that the origin at the diverse suite of aftershock mechanisms Following the 1989 M 7.1 Loma Prieta, California, earthquake is related to the post-main-shock static stress field We use a 3-D boundary-element algorithm to calculate static stresses combined with a Coulomb failure criterion to calculate conjugate failure planes at aftershock locations. The post-main-shock static stress field is taken as the of a sum of a pre-existing stress field and changes in stress due to the heterogeneous slip across the Loma Prieta rupture plane. The background stress field is assumed to be either a simple shear parallel to the regional trend of the San Andreas fault or approximately fault-normal compression A suite of synthetic aftershock mechanisms from the conjugate failure planes is generated and quantitatively compared (allowing for uncertainties in both mechanism parameters and earthquake locations) to well-constrained mechanisms reported in the US Geological Survey Northern California Seismic Network catalogue. We also compare calculated rakes with those observed by resolving the calculated stress tensor onto observed focal mechanism nodal planes, assuming either plane to be a likely rupture plane.

Geophysical Journal International↗