Search USGS⌕ Search

SEARCH · Search USGS

Results for “Geophysics Journal International”

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
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↗
USGS1997

Progressive failure on the North Anatolian fault since 1939 by earthquake stress triggering

10 M ≥ 6.7 earthquakes ruptured 1000 km of the North Anatolian fault (Turkey) during 1939–1992, providing an unsurpassed opportunity to study how one large shock sets up the next. We use the mapped surface slip and fault geometry to infer the transfer of stress throughout the sequence. Calculations of the change in Coulomb failure stress reveal that nine out of 10 ruptures were brought closer to failure by the preceding shocks, typically by 1–10 bar, equivalent to 3–30 years of secular stressing. We translate the calculated stress changes into earthquake probability gains using an earthquake-nucleation constitutive relation, which includes both permanent and transient effects of the sudden stress changes. The transient effects of the stress changes dominate during the mean 10 yr period between triggering and subsequent rupturing shocks in the Anatolia sequence. The stress changes result in an average three-fold gain in the net earthquake probability during the decade after each event. Stress is calculated to be high today at several isolated sites along the fault. During the next 30 years, we estimate a 15 per cent probability of a M ≥ 6.7 earthquake east of the major eastern centre of Ercinzan, and a 12 per cent probability for a large event south of the major western port city of Izmit. Such stress-based probability calculations may thus be useful to assess and update earthquake hazards elsewhere.

Geophysical Journal International↗
USGS2003

Motion of the Scotia sea plates

Earthquake data from the Scotia Arc to early 2002 are reviewed in the light of satellite gravity and other data in order to derive a model for the motion of plates in the Scotia Sea region. Events with magnitude ???5, which occurred on or near the boundaries of the Scotia and Sandwich plates, and for which Centroid Moment Tensor (CMT) solutions are available, are examined. The newer data fill some of the previous sampling gaps along the boundaries of the Scotia and Sandwich plates, and provide tighter constraints on relative motions. Variations in the width of the Brunhes anomaly on evenly spaced marine magnetic profiles over the East Scotia Ridge provide new estimates of Scotia-Sandwich plate spreading rates. Since there are no stable fracture zones in the east Scotia Sea, the mean azimuth of sea floor fabric mapped by sidescan is used to constrain the direction of spreading. 18 new rate estimates and four azimuths from the East Scotia Ridge are combined with 68 selected earthquake slip vectors from the boundaries of the Scotia Sea in a least-squares inversion for the best-fitting set of Euler poles and angular rotation rates describing the 'present-day' motions of the Scotia and Sandwich plates relative to South America and Antarctica. Our preferred model (TLP2003) gives poles that are similar to previous estimates, except for Scotia Plate motion with respect to South America, which is significantly different from earlier estimates; predicted rates of motion also differ slightly. Our results are much more robust than earlier work. We examine the implications of the model for motion and deformation along the various plate boundaries, with particular reference to the North and South Scotia Ridges, where rates are obtained by closure.

Geophysical Journal International↗
USGS2003

Hydrological response to earthquakes in the Haibara well, central Japan - I. Groundwater level changes revealed using state space decomposition of atmospheric pressure, rainfall and tidal responses

For the groundwater level observed at the Haibara well, Shizuoka Prefecture, central Japan, time series analysis using state-space modelling is applied to extract hydrological anomalies related to earthquakes. This method can decompose observed groundwater level time series into five components: atmospheric pressure, tidal, and precipitation responses, observation noise, and residual water level. The decomposed responses to atmospheric pressure and precipitation are independently determined and are consistent with the expected response to surface loading. In the groundwater level at the Haibara well, 28 coseismic changes can be discerned during the period from 1981 April to 1997 December. There is a threshold in the relationship between earthquake magnitude and the well-hypocentre distance, above which earthquakes cause coseismic changes in the residual water level. All of the coseismic water level changes at the Haibara well are decreases, although 33 per cent of the estimated coseismic volumetric strain steps are contraction, which would be expected to cause water level increases. The coseismic changes in groundwater level are more closely proportional to the estimated ground motion than to coseismic volumetric strain steps, suggesting that ground motion due to earthquakes is the major cause of the coseismic water level drops and that the contribution from static strain is rather small. Possible pre- or inter-earthquake water level changes have occurred at the Haibara well and may have been caused by local aseismic crustal deformation.

Geophysical Journal International↗
USGS2003

Hydrological response to earthquakes in the Haibara well, central Japan - II. Possible mechanism inferred from time-varying hydraulic properties

28 coseismic groundwater level decreases have been observed at the Haibara well, Shizuoka prefecture, central Japan, from 1981 to 1997. These groundwater level changes cannot be explained as the poroelastic response to coseismic static strain. We use the atmospheric pressure and tidal responses of the well, rock properties measured on core samples from the same formation and pumping test results to characterize the hydraulic and mechanical properties of the aquifer. The responses of the Haibara well to the M2 Earth tide constituent and to atmospheric pressure have varied over time. In particular, increasing amplitude and decreasing phase lags were observed after the 1993 pumping test, as well as after earthquakes that caused coseismic water level changes. The tidal response, together with the surface load efficiency derived from the atmospheric pressure response, is used to estimate the mechanical properties of the aquifer. The largest amplitude of the M2 constituent, 2.2 mm, is small enough to imply that pore fluid in this system is approximately twice as compressible as water, possibly due to the presence of a small amount of exsolved gas. Diffusion of a coseismic pressure drop near the well could account for the observed time histories of the water level changes. The time histories of the water level drops are well matched by the decay of a coseismic pressure drop at least 80 m away from the well. Removal of a small amount of gas from the formation in that location might in turn explain the coseismic pressure drops.

Geophysical Journal International↗
USGS2003

Quantitative analysis of seismic fault zone waves in the rupture zone of the 1992 Landers, California, earthquake: Evidence for a shallow trapping structure

We analyse quantitatively a waveform data set of 238 earthquakes recorded by a dense seismic array across and along the rupture zone of the 1992 Landers earthquake. A grid-search method with station delay corrections is used to locate events that do not have catalogue locations. The quality of fault zone trapped waves generated by each event is determined from the ratios of seismic energy in time windows corresponding to trapped waves and direct S waves at stations close to and off the fault zone. Approximately 70 per cent of the events with S-P times of less than 2 s, including many clearly off the fault, produce considerable trapped wave energy. This distribution is in marked contrast with previous claims that trapped waves are generated only by sources close to or inside the Landers rupture zone. The time difference between the S arrival and trapped waves group does not grow systematically with increasing hypocentral distance and depth. The dispersion measured from the trapped waves is weak. These results imply that the seismic trapping structure at the Landers rupture zone is shallow and does not extend continuously along-strike by more than a few kilometres. Synthetic waveform modelling indicates that the fault zone waveguide has depth of approximately 2-4 km, a width of approximately 200 m, an S-wave velocity reduction relative to the host rock of approximately 30-40 per cent and an S-wave attenuation coefficient of approximately 20-30. The fault zone waveguide north of the array appears to be shallower and weaker than that south of the array. The waveform modelling also indicates that the seismic trapping structure below the array is centred approximately 100 m east of the surface break.

Geophysical Journal International↗
USGS2004

Systematic variation in the depths of slabs beneath arc volcanoes

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

Geophysical Journal International↗
USGS2004

Estimating tectonic history through basin simulation-enhanced seismic inversion: Geoinformatics for sedimentary basins

A data assimilation approach is demonstrated whereby seismic inversion is both automated and enhanced using a comprehensive numerical sedimentary basin simulator to study the physics and chemistry of sedimentary basin processes in response to geothermal gradient in much greater detail than previously attempted. The approach not only reduces costs by integrating the basin analysis and seismic inversion activities to understand the sedimentary basin evolution with respect to geodynamic parameters-but the technique also has the potential for serving as a geoinfomatics platform for understanding various physical and chemical processes operating at different scales within a sedimentary basin. Tectonic history has a first-order effect on the physical and chemical processes that govern the evolution of sedimentary basins. We demonstrate how such tectonic parameters may be estimated by minimizing the difference between observed seismic reflection data and synthetic ones constructed from the output of a reaction, transport, mechanical (RTM) basin model. We demonstrate the method by reconstructing the geothermal gradient. As thermal history strongly affects the rate of RTM processes operating in a sedimentary basin, variations in geothermal gradient history alter the present-day fluid pressure, effective stress, porosity, fracture statistics and hydrocarbon distribution. All these properties, in turn, affect the mechanical wave velocity and sediment density profiles for a sedimentary basin. The present-day state of the sedimentary basin is imaged by reflection seismology data to a high degree of resolution, but it does not give any indication of the processes that contributed to the evolution of the basin or causes for heterogeneities within the basin that are being imaged. Using texture and fluid properties predicted by our Basin RTM simulator, we generate synthetic seismograms. Linear correlation using power spectra as an error measure and an efficient quadratic optimization technique are found to be most effective in determining the optimal value of the tectonic parameters. Preliminary 1-D studies indicate that one can determine the geothermal gradient even in the presence of observation and numerical uncertainties. The algorithm succeeds even when the synthetic data has detailed information only in a limited depth interval and has a different dominant frequency in the synthetic and observed seismograms. The methodology presented here even works when the basin input data contains only 75 per cent of the stratigraphic layering information compared with the actual basin in a limited depth interval.

Geophysical Journal International↗
USGS2004

Shear wave velocity, seismic attenuation, and thermal structure of the continental upper mantle

Seismic velocity and attenuation anomalies in the mantle are commonly interpreted in terms of temperature variations on the basis of laboratory studies of elastic and anelastic properties of rocks. In order to evaluate the relative contributions of thermal and non-thermal effects on anomalies of attenuation of seismic shear waves, Q −1 s , and seismic velocity, V s , we compare global maps of the thermal structure of the continental upper mantle with global Q −1 s and V s maps as determined from Rayleigh waves at periods between 40 and 150 s. We limit the comparison to three continental mantle depths (50, 100 and 150 km), where model resolution is relatively high. The available data set does not indicate that, at a global scale, seismic anomalies in the upper mantle are controlled solely by temperature variations. Continental maps have correlation coefficients of <0.56 between V s and T and of <0.47 between Q s and T at any depth. Such low correlation coefficients can partially be attributed to modelling artefacts; however, they also suggest that not all of the V s and Q s anomalies in the continental upper mantle can be explained by T variations. Global maps show that, by the sign of the anomaly, V s and Q s usually inversely correlate with lithospheric temperatures: most cratonic regions show high V s and Q s and low T , while most active regions have seismic and thermal anomalies of the opposite sign. The strongest inverse correlation is found at a depth of 100 km, where the attenuation model is best resolved. Significantly, at this depth, the contours of near-zero Q s anomalies approximately correspond to the 1000 °C isotherm, in agreement with laboratory measurements that show a pronounced increase in seismic attenuation in upper mantle rocks at 1000–1100 °C. East–west profiles of V s , Q s and T where continental data coverage is best (50°N latitude for North America and 60°N latitude for Eurasia) further demonstrate that temperature plays a dominant, but non-unique, role in determining the value of lithospheric V s and Q s . At 100 km depth, where the resolution of seismic models is the highest, we compare observed seismic V s and Q s with theoretical V T s and Q T s values, respectively, that are calculated solely from temperature anomalies and constrained by experimental data on temperature dependencies of velocity and attenuation. This comparison shows that temperature variations alone are sufficient to explain seismic V s and Q s in ca 50 per cent of continental regions. We hypothesize that compositional anomalies resulting from Fe depletion can explain the misfit between seismic and theoretical V s in cratonic lithosphere. In regions of active tectonics, temperature effects alone cannot explain seismic V s and Q s in the lithosphere. It is likely that partial melts and/or fluids may affect seismic parameters in these regions. This study demonstrates that lithospheric temperature plays the dominant role in controlling V s and Q s anomalies, but other physical parameters, such as compositional variations, fluids, partial melting and scattering, may also play a significant role in determining V s and Q s variations in the continental mantle.

Geophysical Journal International↗
USGS2004

Flow path of the 1993 Hokkaido-Nansei-oki earthquake seismoturbidite, suthern margin of the Japan sea north basin, inferred from anisotropy of magnetic susceptibility

A magnetic fabric analysis has been carried out on standard cube samples from one gravity and three multiple cores extracted from the Shiribeshi trough and Okushiri basin in the southern margin of the Japan sea north basin. It is aimed at tracing the flow path of turbidites that are assumed to have deposited in response to the 1993 Hokkaido-Nansei-oki earthquake. Magnetic remanence was used for reorientation to the geographic coordinates. Magnetomineralogical investigations including low-temperature magnetometry, magnetic hysteresis loops and isothermal remanent magnetization (IRM) acquisition experiments indicate that pseudosingle domain to multidomain magnetite is the principal magnetic carrier and is, therefore, capable of providing reliable anisotropy of magnetic susceptibility (AMS) palaeocurrent direction estimates. A well-developed near-horizontal magnetic foliation and minimum susceptibility axes lying close to vertical are recorded at all sites reflecting an original depositional fabric. Clearly defined magnetic lineation was observed at all sites and is considered to reflect the palaeocurrent direction. Down-core changes of susceptibility and key AMS parameters show good correspondence to occurrences of turbidite layers marking the increase of input of influx materials. In agreement with results from recent marine surveys and IZANAGI side-scan sonar images, an NNE transportation trend has been estimated for sediments at sites from the Shiribeshi trough with a possible depositing path initiating from the slope bounding the south and southeastern margin down to the trough floor. Similarly, a SSE palaeocurrent direction has been estimated for sediments from the Okushiri basin with evidence for a relatively strong transporting current flowing through the canyons along the steep slope bounding the north and northeastern margins of the basin. The present results agree with the view that slope failure is the most probable mechanism for the down-slope transport of the sand from the shelves and upper slopes down to floors of basins and troughs in the southern margin of the Japan sea north basin. They further support the ongoing assumption that the 1993 Hokkaido-Nansei-oki and other strong historical earthquakes together with associated tsunamis are the principal triggering forces for the down-slope mass gravitational transport and formation of turbidites in this seismically active area. ?? 2004 RAS.

Geophysical Journal International↗
USGS2004

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International↗
USGS2004

Shear wave velocity variation across the Taupo Volcanic Zone, New Zealand, from receiver function inversion

The Taupo Volcanic Zone (TVZ), New Zealand is a region characterized by very high magma eruption rates and extremely high heat flow, which is manifest in high-temperature geothermal waters. The shear wave velocity structure across the region is inferred using non-linear inversion of receiver functions, which were derived from teleseismic earthquake data. Results from the non-linear inversion, and from forward synthetic modelling, indicate low S velocities at ∼6–16 km depth near the Rotorua and Reporoa calderas. We infer these low-velocity layers to represent the presence of high-level bodies of partial melt associated with the volcanism. Receiver functions at other stations are complicated by reverberations associated with near-surface sedimentary layers. The receiver function data also indicate that the Moho lies between 25 and 30 km, deeper than the 15 ± 2 km depth previously inferred for the crust–mantle boundary beneath the TVZ.

Taupo volcanic zone↗
USGS2004

Apparent stress, fault maturity and seismic hazard for normal-fault earthquakes at subduction zones

The behavior of apparent stress for normal-fault earthquakes at subduction zones is derived by examining the apparent stress (?? a = ??Es/Mo, where E s is radiated energy and Mo is seismic moment) of all globally distributed shallow (depth, ?? < 70 km) earthquakes with normal-fault mechanisms that occurred in or near subduction zones between 1987 and 2001 for which Es and M o are available. Accurately determined hypocentres from the Engdahl-Hilst-Buland (EHB) catalogue establish the fine detail of the Wadati-Benioff zone. In many cases, we can relate trends in apparent stress to specific features within the subduction zone and compare these with trends for interplate-thrust earthquakes in the same subduction zones. There are two depth ranges over which ??a, for normal-fault earthquakes shows maxima. The highest and most anomalous values of ??a are found in the deeper depth range from 35 to 70 km. The high-??a, events (up to 5 MPa) are characteristically intraslab and in proximity to zones of intense deformation such as a sharp slab bend or where opposing slabs collide. High-??a events in the region of the shallower maximum (hypocentres between 10-35 km and ??a > 1 MPa) are also generally intraslab, but occur where the lithosphere has just begun subduction beneath the overriding plate. They usually occur in cold slabs near trenches where the direction of plate motion across the trench is oblique to the trench axis, or where there are local contortions or geometrical complexities of the plate boundary. Lower ??a (< 1 MPa) is associated with events occurring at the outer rise (OR) complex (between the OR and the trench axis), as well as with intracrustal events occurring just landward of the trench. The average apparent stress of intraslab-normal-fault earthquakes is considerably higher than the average apparent stress of interplate-thrust-fault earthquakes. In turn, the average ?? a of strike-slip earthquakes in intraoceanic environments is considerably higher than that of intraslab-normal-fault earthquakes. The variation of average ??a with focal mechanism and tectonic regime suggests that the level of ?? a is related to fault maturity. Lower stress drops are needed to rupture mature faults such as those found at plate interfaces that have been smoothed by large cumulative displacements (from hundreds to thousands of kilometres). In contrast, immature faults, such as those on which intraslab-normal-fault earthquakes generally occur, are found in cold and intact lithosphere in which total fault displacement has been much less (from hundreds of metres to a few kilometres). Also, faults on which high ??a oceanic strike-slip earthquakes occur are predominantly intraplate or at evolving ends of transforms. At subduction zones, earthquakes occurring on immature faults are likely to be more hazardous as they tend to generate higher amounts of radiated energy per unit of moment than earthquakes occurring on mature faults. We have identified earthquake pairs in which an interplate-thrust and an intraslab-normal earthquake occurred remarkably close in space and time. The intraslab-normal member of each pair radiated anomalously high amounts of energy compared to its thrust-fault counterpart. These intraslab earthquakes probably ruptured intact slab mantle and are dramatic examples in which Mc (an energy magnitude) is shown to be a far better estimate of the potential for earthquake damage than Mw. This discovery may help explain why loss of life as a result of intraslab earthquakes was greater in the 20th century in Latin America than the fatalities associated with interplate-thrust events that represented much higher total moment release. ?? 2004 RAS.

Geophysical Journal International↗
USGS2004

A viscoelastic damage model with applications to stable and unstable fracturing

A viscoelastic damage rheology model is presented that provides a generalization of Maxwell viscoelasticity to a non-linear continuum mechanics framework incorporating material degradation and recovery, transition from stable to unstable fracturing and gradual accumulation of non-reversible deformation. The model is a further development of the damage rheology framework of Lyakhovsky et al. for evolving effective elasticity. The framework provides a quantitative treatment for macroscopic effects of evolving distributed cracking with local density represented by an intensive state variable. The formulation, based on thermodynamic principles, leads to a system of kinetic equations for the evolution of damage. An effective viscosity inversely proportional to the rate of damage increase is introduced to account for gradual accumulation of irreversible deformation due to dissipative processes. A power-law relation between the damage variable and elastic moduli leads to a non-linear coupling between the rate of damage evolution and the damage variable itself. This allows the model to reproduce a transition from stable to unstable fracturing of brittle rocks and the Kaiser effect. 3-D numerical simulations based on the model formulation for homogeneous and heterogeneous materials account for the main features of rock behaviour under large strain. The model coefficients are constrained, using triaxial laboratory experiments with low-porosity Westerly granite and high-porosity Berea sandstone samples.

Geophysical Journal International↗