Reduced seismic velocities in the source zone of New Madrid earthquakes

1988 ◽  
Vol 78 (4) ◽  
pp. 1491-1509
Author(s):  
Haydar J. Al-Shukri ◽  
Brian J. Mitchell

Abstract A three-dimensional inversion of P-wave travel-time residuals from local earthquakes reveals a remarkable pattern of low seismic velocities in crustal rocks immediately adjacent to the active portions of the New Madrid fault system. Seismic velocities are lowest in regions of greatest concentration of earthquake activity near two intersections of linear trends in seismicity. The maximum reduction in compressional wave velocity is at least 7 per cent in the upper 5 km of the crust and at least 4 per cent in the depth range of 5 to 14 km. The reductions are consistent with a velocity decrease which would be expected if crustal rocks in the source zone contain fluid-filled cracks in which pore pressure is a substantial fraction of external pressure. The presence or absence of such fluids may explain why some portions of the faults in and surrounding the upper Mississippi Embayment are active while others are not.

Author(s):  
Xu Zhang ◽  
Li-Sheng Xu ◽  
Lei Yi ◽  
Wanpeng Feng

Abstract On 8 August 2017, an Ms 7.0 earthquake struck the Jiuzhaigou town, Sichuan Province, China, rupturing an unmapped fault, which is adjacent to the Maqu seismic gap in the Min Shan uplift zone in the easternmost part of the Bayan Har block. Having summarized the previous studies on the source of this earthquake, we confirmed the rupture model by jointly inverting the teleseismic P-wave and SH-wave data, Interferometric Synthetic Aperture Radar line-of-sight displacement data, and the near-field seismic and strong-motion data, a most complete dataset until now. The confirmation showed that a scalar seismic moment of 6.6×1018  N·m was released (corresponding to a moment magnitude of Mw 6.5), and 95% of the release occurred in the first 10 s. The slip area was composed of two asperities, with a horizontal extension of ∼20  km and a depth range of ∼2–15  km. A bilateral extending occurred at shallow depths, but the rupturing upward from deep depth dominated in the early time. The rupture process was found generally featuring the slip-pulse mode, which was related to the weak prestress condition. The aftershocks almost took place in gaps of the mainshock slip because of the coulomb stress change. Combining the aftershock relocations, aftershock focal mechanism solutions, and our confirmed rupture model, we suggest that the seismogenic fault was a northward extension of the mapped Huya fault. The occurrence of this earthquake made the Maqu seismic gap at a higher level of seismic risk, in addition to the moderate to high strain accumulation on the easternmost tip of the Kunlun fault system and the weak lower crust below.


1985 ◽  
Vol 75 (5) ◽  
pp. 1285-1300
Author(s):  
Alan L. Kafka ◽  
Ellyn A. Schlesinger-Miller ◽  
Noel L. Barstow

Abstract Earthquakes recorded by stations of the Lamont-Doherty seismic network in the greater New York City area are analyzed to determine magnitudes and the relationship between seismicity and geologic structures. Between 1974 and 1983, the configuration of stations in this region remained relatively stationary and the type of recording devices (visual drum recorders and 16-mm photographic recorders) did not change. This distribution of stations and recording devices allows for a uniform measurement of magnitudes and seismicity. Magnitudes of these earthquakes are determined by comparing amplitudes and signal duration measured from high-frequency (5 to 10 Hz) data recorded by the local network with mbLg and ML determined from data at frequencies near 1 Hz. During the period of time studied (nearly 10 yr), 61 earthquakes were located in this region, but none of these earthquakes exceeded 3.0 on the mbLg scale. The largest event (mbLg = 3.0) occurred in the Coastal Plain province of northern New Jersey. The magnitude threshold for uniform detection of events throughout this region during the period of time studied is estimated to be mbLg = 1.6. With events below this threshold removed from the catalog of network seismicity, we find that about half of the earthquakes studied occurred within 10 km of the Ramapo fault system. This fault system lies about 30 km northwest of New York City and has been interpreted by several investigators to be the most active fault system in the greater New York City area. However, earthquakes at least as large as those recorded near the Ramapo fault were located as far as 50 km from this fault (and within 20 km of New York City), in geologic structures that surround the Newark basin. While the Ramapo fault can by no means be ruled out as a possible source zone for earthquakes in the greater New York City area, the geologic structures associated with most (if not all) earthquakes in this region are still unknown. Thus, the cause of earthquakes in this region remains an enigma.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. E1-E13 ◽  
Author(s):  
Michael A. Zimmer ◽  
Manika Prasad ◽  
Gary Mavko ◽  
Amos Nur

Knowledge of the pressure dependences of seismic velocities in unconsolidated sands is necessary for the remote prediction of effective pressures and for the projection of velocities to unsampled locations within shallow sand layers. We have measured the compressional- and shear-wave velocities and bulk, shear, and P-wave moduli at pressures from [Formula: see text] in a series of unconsolidated granular samples including dry and water-saturated natural sands and dry synthetic sand and glass-bead samples. The shear-wave velocities in these samples demonstrate an average pressure dependence approximately proportional to the fourth root of the effective pressure [Formula: see text], as commonly observed at lower pressures. For the compressional-wave velocities, theexponent in the pressure dependence of individual dry samples is consistently less than the exponent for the shear-wave velocity of the same sample, averaging 0.23 for the dry sands and 0.20 for the glass-bead samples. These pressure dependences are generally consistent over the entire pressure range measured. A comparison of the empirical results to theoretical predictions based on Hertz-Mindlin effective-medium models demonstrates that the theoretical models vastly overpredict the shear moduli of the dry granular frame unless the contacts are assumed to have no tangential stiffness. The models also predict a lower pressure exponent for the moduli and velocities [Formula: see text] than is generally observed in the data. We attribute this discrepancy in part to the inability of the models to account for decreases in the amount of slip or grain rotation occurring at grain-to-grain contacts with increasing pressure.


Geophysics ◽  
1976 ◽  
Vol 41 (5) ◽  
pp. 837-849 ◽  
Author(s):  
Robert H. Tatham ◽  
Paul L. Stoffa

Theoretically and experimentally, the shear‐wave velocity of a porous rock has been shown to be less sensitive to fluid saturants than the compressional wave velocity. Thus, observation of the ratio of the seismic velocities for waves which traverse a changing or laterally varying zone of undersaturation or gas saturation could produce an observable anomaly which is independent of the regional variation in compressional wave velocity. One source of shear‐wave data in reflection seismic prospecting is mode conversion of P waves to shear waves in marine areas of high water bottom P-wave velocity. A relatively simple interpretative technique, based on amplitude variation as a function of the angle of incidence, is a possible discriminant between shear and multiple compressional arrivals, and data for a real case are shown. A normal moveout velocity analysis, carefully coupled with this offset discriminant, leads to the construction of a shear‐wave reflection section which can then be correlated with the usual compressional wave section. One such a section has been constructed, the variation in the ratio of the seismic velocities can be mapped, and potentially anomalous subsurface regions observed.


2020 ◽  
Vol 223 (3) ◽  
pp. 1780-1793
Author(s):  
Theresa Rein ◽  
Katrin Hannemann ◽  
Christine Thomas ◽  
Michael Korn

SUMMARY Receiver function approaches have proven to be valuable for the investigation of crustal and upper mantle discontinuities whose sharp changes in seismic velocities cause wave conversions. While the crustal and mantle transition zone discontinuities are largely understood, the X-discontinuity at 250–350 km depth is still an object of controversial debate. The origin and global distribution of this structure with a velocity jump of 1.5–4.8 % for compressional and shear waves is still unexplained. Although the crustal and mantle transition zone discontinuities beneath SW Morocco and surroundings have been investigated, only a few studies observed the X-discontinuity and place the depth at 260–370 km beneath the region of western Morocco. In order to better locate and characterize the X-discontinuity beneath southwest Morocco, we create P-wave receiver functions using data recorded by the Morocco–Münster array and detect the X-discontinuity at apparent depths of 285–350 km. In the western part of our study region we find apparent depths of ∼ 310–340 km. The eastern part of the study area appears more complex: we locate two velocity jumps at apparent depths of around 285–295 km and 330–350 km in the northeast, and in the southeast we find a discontinuity at apparent depths of 340–350 km. Due to the large depth range and the twofold appearance of the X-discontinuity, we suggest that two different phase transitions cause the X-discontinuity beneath SW Morocco. The velocity contrasts at larger depths likely point to the coesite–stishovite phase transition occurring in deep eclogitic pools. The shallower depths can be explained by the transition from orthoenstatite to high-pressure clinoenstatite which requires the reaction between eclogite and peridotite to form orthopyroxene-rich peridotite. This reaction is likely related to previously proposed small-scale mantle upwellings beneath SW Morocco. Since both phase transitions require eclogite occurrence, the location of the X-discontinuity in this region can be used to indicate the location of recycled oceanic crust.


Geophysics ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. E15-E25 ◽  
Author(s):  
Michael A. Zimmer ◽  
Manika Prasad ◽  
Gary Mavko ◽  
Amos Nur

Unaccounted-for porosity variation in unconsolidated sediments can cloud the interpretation of the sediment’s seismic velocities for factors such as fluid content and pressure. However, an understanding of the effects of porosity variation on the velocities can permit the remote characterization of porosity with seismic methods. We present the results of a series of measurements designed to isolate the effects of sorting- and compaction-induced porosity variation on the seismic velocities and their pressure dependences in clean, unconsolidated sands. We prepared a set of texturally similar sand and glass-bead samples with controlled grain-size distributions to cover an initial porosity range from 0.26 to 0.44. We measured the compressional- and shear-wave velocities and porosity of dry samples over a series of hydrostatic pressure cycles from [Formula: see text]. Over this rangeof porosities, the velocities of the dry samples at a given pressure vary by [Formula: see text]. However, the water-saturated compressional-wave velocities, modeled with Gassmann fluid substitution, demonstrate a consistent increase with decreasing porosity. In both the dry and water-saturated cases, the porosity trend at a given pressure is approximately described by the isostress (harmonic) average between the moduli of the highest-porosity sample at that pressure and the moduli of quartz, the predominant mineral component of the samples. Empirical power-law fit coefficients describing the pressure dependences of the dry bulk, shear, and constrained (P-wave) moduli from each sample also demonstrate no significant, systematic relationship with the porosity. The porosity dependence of the water-saturated bulk and constrained moduli is primarily contained in the empirical coefficient representing the modulus at zero pressure.


Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.


2010 ◽  
Vol 47 (4) ◽  
pp. 389-408 ◽  
Author(s):  
Claire Perry ◽  
Carmen Rosieanu ◽  
Jean-Claude Mareschal ◽  
Claude Jaupart

Geothermal studies were conducted within the framework of Lithoprobe to systematically document variations of heat flow and surface heat production in the major geological provinces of the Canadian Shield. One of the main conclusions is that in the Shield the variations in surface heat flow are dominated by the crustal heat generation. Horizontal variations in mantle heat flow are too small to be resolved by heat flow measurements. Different methods constrain the mantle heat flow to be in the range of 12–18 mW·m–2. Most of the heat flow anomalies (high and low) are due to variations in crustal composition and structure. The vertical distribution of radioelements is characterized by a differentiation index (DI) that measures the ratio of the surface to the average crustal heat generation in a province. Determination of mantle temperatures requires the knowledge of both the surface heat flow and DI. Mantle temperatures increase with an increase in surface heat flow but decrease with an increase in DI. Stabilization of the crust is achieved by crustal differentiation that results in decreasing temperatures in the lower crust. Present mantle temperatures inferred from xenolith studies and variations in mantle seismic P-wave velocity (Pn) from seismic refraction surveys are consistent with geotherms calculated from heat flow. These results emphasize that deep lithospheric temperatures do not always increase with an increase in the surface heat flow. The dense data coverage that has been achieved in the Canadian Shield allows some discrimination between temperature and composition effects on seismic velocities in the lithospheric mantle.


1990 ◽  
Vol 80 (2) ◽  
pp. 430-449 ◽  
Author(s):  
William Menke ◽  
Arthur L. Lerner-Lam ◽  
Bruce Dubendorff ◽  
Javier Pacheco

Abstract Except for its very onset, the P wave of earthquakes and chemical explosions observed at two narrow-aperture arrays on hard-rock sites in the Adirondack Mountains have a nearly random polarization. The amount of energy on the vertical, radial, and transverse components is about equal over the frequency range 5 to 30 Hz, for the entire seismogram. The spatial coherence of the seismograms is approximately exp(−cfΔx), where c is in the range 0.4 to 0.7 km−1Hz−1, f is frequency and Δx is the distance between array elements. Vertical, radial, and transverse components were quite coherent over the aperture of the array, indicating that the transverse motion of the compressional wave is a property of relatively large (106 m3) volumes of rock, and not just an anomaly caused by a malfunctioning instrument, poor instrument-rock coupling, or out-crop-scale effects. The spatial coherence is approximately independent of component, epicentral azimuth and range, and whether P- or S-wave coda is being considered, at least for propagation distances between 5 and 170 km. These results imply a strongly and three-dimensionally heterogeneous crust, with near-receiver scattering in the uppermost crust controlling the coherence properties of the waves.


1976 ◽  
Vol 66 (6) ◽  
pp. 1931-1952
Author(s):  
Donald J. Stierman ◽  
William L. Ellsworth

abstract The ML 6.0 Point Mugu, California earthquake of February 21, 1973 and its aftershocks occurred within the complex fault system that bounds the southern front of the Transverse Ranges province of southern California. P-wave fault plane solutions for 51 events include reverse, strike slip and normal faulting mechanisms, indicating complex deformation within the 10-km broad fault zone. Hypocenters of 141 aftershocks fail to delineate any single fault plane clearly associated with the main shock rupture. Most aftershocks cluster in a region 5 km in diameter centered 5 km from the main shock hypocenter and well beyond the extent of fault rupture estimated from analysis of body-wave radiation. Strain release within the imbricate fault zone was controlled by slip on preexisting planes of weakness under the influence of a NE-SW compressive stress.


Sign in / Sign up

Export Citation Format

Share Document