bluuebunny/crossref_metadata_2025_split huggingface dataset:
Long-period ground motions and dynamic strain field of Los Angeles basin during large earthquakes
Authors: Teng T. L. , Qu J. 
Published: September 2021

AbstractDuring a big earthquake along the San Andreas fault in southern California, high excitation and low attenuation of long-period (3 to 10 sec) strong ground motions will cause wave motions to propagate efficiently far from the epicentral area. These ground motions could potentially be destructive to large-dimension structures in the Los Angeles basin. We performed calculations using the surface-wave Gaussian beam method for a 3D southern California crustal structure. Displacement field as well as the associated dynamic strain field produced by large propagating ruptures along the San Andreas fault are evaluated. Results indicate that in the presence of lateral heterogeneity, focusing and multipathing interference contribute significantly to a complex pattern of the displacement field and the associated dynamic strain field. For a big event on the San Andreas fault with a seismic moment of 1.8 × 1028 dyne-cm, long-period displacement in the Los Angeles basin could reach a maximum amplitude of meters in places. Since this calculation is fast, we have evaluated the displacement field for a dense grid of points; a differentiation gives the corresponding effective horizontal dynamic strain field. At times, the maximum effective dynamic strains may reach mid-10−3 to even 10−3—high enough to be of engineering concern. This computational result probably gives the upper bound values due to the large source assumed. For events of smaller seismic moment release along less extensive ruptures, these results can easily be scaled down proportionally. Different scenarios are considered in this study with different slip distributions. It is found that with a given seismic moment, a more evenly distributed fault slip over the rupture surface will result in lower peak values on both displacements and dynamic strains. Our displacement results give similar values to those obtained by Kanamori using empirical Green's functions but substantially higher than Bouchon and Aki's results.
High-frequency analysis of seismic background noise as a function of wind speed and shallow depth
Authors: Withers Mitchell M. , Aster Richard C. , Young Christopher J. , Chael Eric P. 
Published: September 2021

Abstract
We used a deep (1500 m) cased borehole near the town of Datil in west-central New Mexico to study high-frequency (&gt;1 Hz) seismic noise characteristics. The remote site had very low levels of cultural noise, but strong winds (winter and spring) made the site an excellent candidate to study the effects of wind noise on seismograms. Along with a three-component set of surface sensors (Teledyne Geotech GS-13), a vertical borehole seismometer (GS-28) was deployed at a variety of depths (5, 43, and 85 m) to investigate signal and noise variations. Wind speed was measured with an anemometer. Event-triggered and time-triggered data streams were recorded on a RefTek 72-02 data acquisition system located at the site. Our data show little cultural noise and a strong correlation between wind speed and seismic background noise. The minimum wind speed at which the seismic background noise appears to be influenced varies with depth: 3 m/sec at the surface, 3.5 m/sec at 43 m in depth, and 4 m/sec at 85 m in depth. For wind speed below 3 to 4 m/sec, we observe omni-directional background noise that is coherent at frequencies below 15 Hz. This coherence is destroyed when wind speeds exceed 3 to 4 m/sec. We use a test event (Md ∼ 1.6) and superimposed noise to investigate signal-to-noise ratio (SNR) improvement with sensor depth. For the low Q valley fill of the Datil borehole (DBH) site, we have found that SNR can be improved by as much as 20 to 40 dB between 23 and 55 Hz and 10 to 20 dB between 10 and 20 Hz, by deploying at a 43-m depth rather than at the surface. At the surface, there is little signal above noise in the 23- to 55-Hz frequency band for wind speeds greater than 8 m/sec. Thus, high-frequency signal information that is lost at the surface can be recorded by deploying at the relatively shallow depth of 40 m. Because we observe only minor further reductions in seismic background noise (SBN) at deeper depths, 40 m is likely to be a reasonable deployment depth for other high-frequency-monitoring sites in similar environmental and geologic conditions.
A statistical technique for validating velocity models
Authors: Iversen Edwin  S. , Lees Jonathan  M. 
Published: September 2021

Abstract
This study investigates the use of a station influence statistic to identify velocity model shortcomings in the earthquake hypocenter location problem. Two groups of microearthquake events are examined. The first is a group of 81 events from the Mount St. Helens region that occurred between November 1987 and September 1991; the second, 110 well-located events from the 1992 Joshua Tree after-shock sequence. We describe a method for validating a postulated earth model. Let λ denote the hypocenter estimates that Geiger's method obtains. Systematically remove each station observation from the location problem, and recompute the location estimate. Call this estimate λ(i) when the ith station is removed. For a single event, define a station's influence (SI) as a weighted difference between λ and λ(i). Distributional summaries of SI statistics across events are used to identify model shortcomings: Given a specified velocity model, SI distributions that are not homogeneous across stations provide evidence of model inadequacies and/or failures in the weighting scheme. We show that velocity model shortcomings detected using SI statistics for the Mount St. Helens sequence under a one-dimensional model appear to correlate with known physical anomalies; while SI distributions evaluated under a three-dimensional model are more homogeneous and reflect a modest improvement over the one-dimensional model. SI distributions provide evidence of model failure for the Joshua Tree sequence under a one-dimensional model, but no evidence of failure under a three-dimensional model. Finally, the weighting scheme's validity is verified for the Joshua Tree sequence under the three-dimensional model.
Predictability of strong motions from the Northridge, California, earthquake
Authors: Anderson John G. , Yu Guang 
Published: September 2021

AbstractThe composite source model for generating synthetic strong ground motions is tested for its ability to predict the statistical characteristics of Northridge accelerograms recorded in or adjacent to the San Fernando valley. The general problem is prediction of strong motions at a site of engineering interest with sufficient realism to be useful for engineering applications. The strongest test of any proposed method is a blind prediction. For this study, a completely blind test was not possible. Our objective was to use only a preliminary description of fault geometry and magnitude and previously published velocity models and, without iteration to improve the quality of fit, to evaluate the differences between predicted and observed accelerograms.The parameters that we predict are peak acceleration, peak velocity, peak displacement, Fourier spectra at seven frequencies, and pseudorelative velocity response (5% damping) at seven periods. Our results are given for 14 stations. For the horizontal components, these parameters are all predicted with a maximum bias of under 50% and an average bias of observations exceeding predictions by 6%. For peak acceleration and some response spectral periods, the bias for this model is smaller than at least some regressions, when applied to this specific earthquake. On the vertical component, the maximum bias is a factor of 2, and the average gives predictions exceeding observations by 25%. Standard deviations of the common logarithm of the ratio of observed-to-predicted parameters are typically about 0.3, which is perhaps 50% greater than the standard deviations typical of regressions but comparable to standard deviations of observations from this earthquake compared to regressions. In the future, it is likely that, in some cases, traditional regressions will be replaced with synthetic calculations of some type, such as the method used here. Based on the results of this study, the amount of progress that has been made in obtaining that goal is very encouraging.
Simulation of long-period ground motion near a large earthquake
Authors: Takeo Minoru , Kanamori Hiroo 
Published: September 2021

AbstractWe estimated the possible range of long-period ground motion for sites located on a soft sedimentary basin in the immediate vicinity of a large earthquake. Since many large cities in the world (e.g., Los Angeles, San Francisco, and Tokyo) where many large structures have been recently constructed are located in this type of environment, a better understanding of long-period ground motion is becoming increasingly important. Our objective is to estimate the possible range of long-period ground motion, rather than ground motion for a specific fault model. We computed ground-motion time series and pseudo-velocity response spectra (PVS) for more than 5,000 models for the 1923 Kanto, Japan, earthquake (MW = 7.9) using 180 slip distributions, eight rupture geometry, and rupture velocities ranging from 1.5 to 3.0 km/sec. Two seismograms recorded in Tokyo during the 1923 Kanto earthquake are used for comparison. The response spectra computed using seismologically reasonable sets of source parameters for the 1923 Kanto earthquake vary by more than an order of magnitude. At periods of 10 to 13 sec, they range from 25 to 170 cm/sec in Tokyo. For some combinations of model parameters, the response spectra exhibit peaks in the range of 10 to 13 sec. Many of the computed response spectra have peaks at periods longer than 10 sec, which is considerably longer than the dominant period (6 to 8 sec) estimated from studies of small earthquakes and microtremor measurements. Thus, the dominant period of the subsurface structure determined locally may not be representative of the dominant period of ground motion from a nearby large earthquake, which is controlled by rupture directivity and source depth. We performed a similar simulation for a hypothetical MW = 7.5 earthquake located beneath the Los Angeles basin. For a site just above the center of the fault, the ground-motion spectral amplitude at a period of 10 sec can vary from 50 to 350 cm/sec. This range, though very large, is what is expected for a seismologically plausible range of source parameters.
An automated Rayleigh-wave detection algorithm
Authors: Chael Eric P. 
Published: September 2021

Abstract
The desire to operate denser networks in order to monitor seismic activity at lower thresholds leads to greater emphasis on automated data processing. An algorithm for detecting and characterizing long-period Rayleigh-wave arrivals has been developed and tested. The routine continuously monitors all directions of approach to a station, in a manner similar to beamforming. The detector is based on cross-powers between the Hilbert-transformed vertical and rotated horizontal signals, so it is sensitive to both the power and polarization properties of the three-component wave field. Elliptically polarized Rayleigh arrivals are enhanced, while linearly polarized Love waves and body phases are suppressed. A test using one month of data from station ANMO demonstrated that this technique can, with high reliability, detect Rayleigh arrivals that are visible on the records. The measured arrival times and azimuths are accurate enough to permit automated association of the detections to events in a bulletin.
Upper-mantle shear velocities beneath southern California determined from long-period surface waves
Authors: Polet J. , Kanamori H. 
Published: September 2021

Abstract
We used long-period surface waves from teleseismic earthquakes recorded by the TERRAscope network to determine phase velocity dispersion of Rayleigh waves up to periods of about 170 sec and of Love waves up to about 150 sec. This enabled us to investigate the upper-mantle velocity structure beneath southern California to a depth of about 250 km. Ten and five earthquakes were used for Rayleigh and Love waves, respectively. The observed surface-wave dispersion shows a clear Love/Rayleigh-wave discrepancy that cannot be accounted for by a simple isotropic velocity model with smooth variations of velocity with depth. Separate isotropic inversions for Love- and Rayleigh-wave data yield velocity models that show up to 10% anisotropy (transverse isotropy). However, tests with synthetic Love waves suggest that the relatively high Love-wave phase velocity could be at least partly due to interference of higher-mode Love waves with the fundamental mode. Even after this interference effect is removed, about 4% anisotropy remains in the top 250 km of the mantle. This anisotropy could be due to intrinsic anisotropy of olivine crystals or due to a laminated structure with alternating high- and low-velocity layers. Other possibilities include the following: upper-mantle heterogeneity in southern California (such as the Transverse Range anomaly) may affect Love- and Rayleigh-wave velocities differently so that it yields the apparent anisotropy; higher-mode Love-wave interference has a stronger effect than suggested by our numerical experiments using model 1066A. If the high Love-wave velocity is due to causes other than anisotropy, the Rayleigh-wave velocity model would represent the southern California upper-mantle velocity structure. The shear velocity in the upper mantle (Moho to 250 km) of this structure is, on average, 3 to 4% slower than that of the TNA model determined for western North America.
The Kozani-Grevena (Greece) earthquake of 13 May 1995 revisited from a detailed seismological study
Authors: Hatzfeld D. , Karakostas V. , Ziazia M. , Selvaggi G. , Leborgne S. , Berge C. , Guiguet R. , Paul A. , Voidomatis P. , Diagnourtas D. , Kassaras I. , Koutsikos I. , Makropoulos K. , Azzara R. , Di Bona M. , Baccheschi S. , Bernard P. , Papaioannou C. 
Published: September 2021

AbstractThe Kozani earthquake (Ms = 6.6) of 13 May 1995 is the strongest event of the decade in Greece and occurred in a region of low seismic activity. Using regional data and the strong-motion record at the Kozani station, we relocate the mainshock at 40.183° N and 21.660° E, beneath the Vourinos massif at a depth of 14.2 km. We also compute a focal mechanism by body-waveform modeling at teleseismic distance, which confirms a normal mechanism. The most likely plane strikes 240° ± 1° N and dips 40° ± 1° N with a centroid depth of 11 ± 1 km. Modeling of the strong-motion record at Kozani confirms that nucleation started at the eastern termination of the bottom of the fault.Six days after the mainshock, we installed a network of 40 portable seismological stations for one week around the epicentral region. Several thousand aftershocks were recorded, among which we locate 622 with a precision better than 1 km. We compute 181 focal mechanisms that mostly show normal faulting. The aftershock seismicity is restricted between 5 and 15 km depth and defines a plane dipping north at an angle of about 35°, consistent with the mainshock mechanism. Seismic activity with the same pattern of normal fault mechanisms is also seen on an antithetic fault connected to the main one at 12 km depth, which cuts the ground surface north of the Vourinos ophiolite massif in the Siatista valley. These results suggest two possibilities for the active fault plane; either it is the Deskati fault that is flat and dips with a constant angle, and therefore the surface breaks are secondary features, or, more likely, it is the Paleohori fault that is new, of listric shape, and located ahead of the Deskati fault, which was not active during the earthquake.
Toward a homogeneous moment-magnitude determination for earthquakes in Greece and the surrounding area
Authors: Papazachos B. C. , Kiratzi A. A. , Karacostas B. G. 
Published: September 2021

Abstract
A strategy is presented in order to produce a statistically homogeneous catalog of historical and present-century earthquakes of the Aegean area by expressing the size of all these earthquakes in the moment-magnitude scale. Records of the Mainka and Wiechert intermediate period (T0 ∼ 4 sec) seismographs in Athens have been used to calculate the magnitudes of earthquakes that occurred in Greece and its surrounding area (Albania, S. Yugoslavia, S. Bulgaria, and W. Turkey) during the whole instrumental period (since 1911). This magnitude, M, has also been used to derive empirical relations with macroseismic data. It is shown here that M is equivalent to moment magnitude, MW, for a wide range of earthquake sizes (5.0 ≦ M ≦ 8.0), and to the surface-wave magnitude, MS, for large earthquake sizes (6 ≦ M ≦ 8.0). It is also shown that the local magnitude, ML, for earthquakes in Greece is half a unit smaller than the moment magnitude in the low-magnitude range (4.5 ≦ ML ≦ 6.0). This could be attributed to a magnification smaller than the nominal for the Wood-Anderson seismographs in Greece (V ∼ 1000). Relations are also proposed to calculate moment magnitude when other magnitudes (M, ML, MS, and mb) are known for earthquakes in Greece.
Analysis of topographic amplification of inclined shear waves in a steep coastal bluff
Authors: Ashford Scott A. , Sitar Nicholas 
Published: September 2021

AbstractThe effect of inclined shear waves on the seismic response of a steep bluff is analyzed using generalized consistent transmitting boundaries. The results of the frequency-domain analysis of a stepped half-space subjected to incident shear waves inclined from 0° to 30° show that the motion at the crest of the slope is amplified for waves traveling into the slope and attenuated for waves traveling away from the slope, as compared to the motion in the free field behind the crest of the slope. This amplification can be as much as twice that observed for vertically propagating waves. A time-domain analysis of bluffs at Seacliff State Beach, California, is used to estimate the effect of topography using realistic conditions, taking into account wave inclination and site effects. The analysis of the site shows that although topographic amplification does in fact nearly double the amplitude of the motion in some cases, this amplification is offset by reduced site amplification and by wave splitting at material interfaces. Thus, the actual peak acceleration occurring at the crest of the slope changes little with incident angle as compared to the amplification of the free-field motion and actually decreases in many cases. Though a more general study is recommended, these results suggest that wave orientation and inclination substantially increase topographic amplification; however, it may be adequate to only account for vertically propagating waves for site response and slope stability analyses where only the magnitude of acceleration is considered.
The large aftershocks of the Northridge earthquake and their relationship to mainshock slip and fault-zone complexity
Authors: Dreger Douglas 
Published: September 2021

AbstractThe distribution of fault slip for the Northridge mainshock and two MW ≧ 6 aftershocks is estimated using an approach that inverts the far-field momentrate functions derived by empirical Green's function deconvolution. The slip for these events was found to have a strong negative correlation with the locations of small aftershocks, and there is a notable clustering of aftershocks near the edges of slip, suggesting that these events are triggered by a redistribution of stress following the mainshock rupture. The two largest aftershocks ruptured the updip regions adjacent to the leading edge of the mainshock rupture. One aftershock, 11 hr after the main event, initiated 8 km west of the mainshock slip terminus and ruptured east to the western edge of mainshock slip. This event was found to have ruptured the same south-dipping fault as the mainshock. The two events, however, are separated by a lateral offset in the depth of seismicity that serves as a geometrical barrier to through-going rupture. An MW 6.2 aftershock approximately 1 min after the main event is inferred to have ruptured an updip region of the mainshock fault that is relatively devoid of small aftershock seismicity. This event was located at the eastern edge of the sequence that overlaps with the north-dipping fault responsible for the 1971 San Fernando earthquake. Although the south-dipping plane is favored for this aftershock, it is not possible to rule out that it occurred on the north-dipping fault structure, or within the hanging block of the mainshock.
Surface strain accumulation and the seismic moment tensor
Authors: Savage J. C. , Simpson R. W. 
Published: September 2021

Abstract
Although the scalar moment accumulation rate within the seismogenic zone beneath a given area is sometimes deduced from the observed average surface strain accumulation rate over that same area (e.g., Working Group on California Earthquake Probabilities, 1995), the correspondence between the two is very uncertain. The equivalence between surface strain accumulation and scalar moment accumulation is based on Kostrov's (1974) relation between the average strain rate over a volume and the moment-rate tensor for that volume. The average strain rate over the volume is replaced by the average strain rate measured at the free surface to deduce an approximate moment-rate tensor. Only in exceptional circumstances will that moment-rate tensor correspond to a double-couple mechanism, a mechanism that can be represented by a scalar moment accumulation rate. More generally, the moment tensor must be resolved into the superposition of two or more double-couple mechanisms, and that resolution can be done in many ways, each with its own scalar moment rate. Thus the resolution is not unique. This is demonstrated by deducing scalar moment accumulation rates for a GPS network that covers most of California south of San Francisco. It is shown that resolutions into different double-couple mechanisms lead to scalar moment accumulation rates differing by factors of ∼2. We suggest that the minimum scalar moment rate equivalent to principal surface strain rates ɛ1 and ɛ2 acting over the area A is M0(min) = 2μHA Max (¦ɛ1¦, ¦ɛ2¦, ¦ɛ1 + ɛ2¦), where μ is the rigidity and H the depth of seismogenic zone, and the function Max is equal to the largest of its arguments. Within the uncertainites of measurement, the scalar moment accumulation rate in southern California based on that approximation is in balance with the average historic seismic moment release rate so that no current earthquake deficit need be accumulating.
Seismicity at Warna reservoir (near Koyna) through 1995
Authors: Rastogi B. K. , Chadha R. K. , Sarma C. S. P. , Mandal P. , Satyanarayana H. V. S. , Raju I. P. , Kumar Narendra , Satyamurthy C. , Nageswara Rao A. 
Published: September 2021

Abstract
The Koyna Dam site in western India has been seismically active since its impounding in 1962. The seismicity peaked in 1967 when an earthquake of magnitude 6.3 claimed 200 human lives and caused some damage to the dam itself. Subsequently, bursts of seismicity with M ≧ 5 earthquakes have occurred during 1973, 1980, and 1993 to 1995. The epicenters during the latest burst are concentrated around the newly impounded Warna reservoir, 25 km south of Koyna. Several shocks of smaller magnitude have occurred in this area since 1967 under the influence of the Koyna (Shivajisagar) reservoir, but the impoundment of the Warna reservoir to over 60 m has led to the enhanced seismic activity in the region. The space-time pattern of epicenters located for the period 28 August 1993 to 31 December 1995 clearly delineates two parallel trends in the NNE-SSW direction. The correlation between seismic activity and water levels in both the reservoirs suggests that the present activity was influenced more by the Warna than the Koyna reservoir. The discriminatory characteristics of reservoir-induced seismicity based on b-value, foreshock-aftershock pattern, and decay rate of aftershocks continue to be valid.