Bulletin of the Seismological Society of America current issue

[Articles] An Overview of the Damaging and Low Magnitude Mw 4.8 La Paca Earthquake on 29 January 2005: Context, Seismotectonics, and Seismic Risk Implications for Southeast Spain

2007-06-11

This article presents an overview of the La Paca earthquake of magnitude m_bLg 4.7, which occurred on 29 January 2005, with its epicenter located near the town of Avilés in the Murcia region in southeast Spain. Despite its low magnitude, the earthquake caused important damage in two towns of the epicentral area, La Paca and Zarcilla de Ramos. These areas recorded intensities of VI–VII (European Macroseismic Scale, 1998) and sustained estimated economic losses amounting to 10 million ;eu[|20ac|]. Aftershocks continued for more than 2 weeks, producing considerable alarm in the population and mobilizing emergency services from the whole region. The La Paca seismic series is the third registered in the region in the past 8 years, being preceded by the Mula (1999) and southwest Bullas (2002) seismic series. These main events had also low magnitudes (m_bLg 4.8) and caused damage levels similar to the 2005 earthquake. The case is an example of a moderate seismic zone where low-magnitude and frequent earthquakes have important implications on the seismic hazard and risk of the region. Although these are not the largest expected earthquakes, they have yielded important information for improving the knowledge of the seismic activity of the area. With this aim in mind, different topics have been analyzed from a multidisciplinary perspective, including seismicity, local tectonics and surface geology, focal mechanisms, macroseismic effects, and ground motion. Results indicate a local tectonic interpretation, consistent with a strike-slip focal mechanism, the confirmation of a triggering process between the 2002 and 2005 earthquakes, a geotechnical and ground-motion characterization for the damaged sites (supporting local amplification effects and estimated peak ground acceleration values of ~0.1g), and an understanding of damage patterns in relation to local building trends. The results may be used as guidelines for future revisions of the Spanish Building Code (Norma de la Construcción Sismorresistente Española [NCSE-02], 2002). The study results should contribute to risk mitigation in a region where strong-motion records from the maximum expected earthquakes are not available. This approach can be extended to other regions with similar seismic backgrounds and a lack of strong-motion records.

[Articles] The 2004 Mw 4.4 Rotenburg, Northern Germany, Earthquake and Its Possible Relationship with Gas Recovery

2007-06-11

We study the 20 October 2004 M_w 4.4 Rotenburg (Wümme)/Neuenkirchen earthquake, located in a previously aseismic region in the northern German sedimentary basin. We constrain the source parameter by using different techniques. A possible relationship between this event, the regional tectonic setting, and local gas recovery is investigated. Different waveform inversion and modeling approaches constrain the depth of the mainshock between 5 and 7 km. The source mechanism was oblique normal faulting on planes striking roughly north–south. An inversion for kinematic rupture parameters indicates a unilateral rupture propagation toward the north, consistent with the higher macroseismic intensities found toward the north in the region of Hamburg compared with those at a similar distance toward the south in the region of Hannover. Relocations of the mainshock and three of the largest aftershocks indicate that these events occurred within a few kilometers of three major gas fields and at depth close to gas production intervals. Comparison with seismicity triggered in the northern Netherlands by depletion of similar gas reservoirs in a similar tectonic environment suggests that the M_w 4.4 Rotenburg event may be related to gas recovery.

Online material: Focal mechanism and waveform fit.

[Articles] Tsunami Hazard Evaluation of the Eastern Mediterranean: Historical Analysis and Selected Modeling

Salamon, A., Rockwell, T., Ward, S. N., Guidoboni, E., Comastri, A. — 2007-06-11

Seismic sea waves in the eastern Mediterranean have been reported since written history first emerged several thousand years ago. We collected and investigated these ancient and modern reports to understand and model the typical tsunamigenic sources, with the ultimate purpose of characterizing tsunami hazard along the Levant coasts. Surprisingly, only 35% of the tsunami reports could be traced back to primary sources, with the balance remaining questionable. The tsunamis varied in size, from barely noticeable to greatly damaging, and their effects ranged from local to regional. Overall, we list 21 reliably reported tsunamis that occurred since the mid second century b.c. along the Levant coast, along with 57 significant historical earthquakes that originated from the "local" continental Dead Sea Transform (dst) system. An in-depth evaluation shows that 10 tsunamis are clearly associated with on-land dst earthquakes, and therefore, as formerly suggested, they probably originated from offshore, seismogenically induced slumps. Eight tsunamis arrived from the "remote" Hellenic and Cypriot Arcs, one from Italy, and two are left with as yet unrecognized sources. A major conclusion from this work is that onshore earthquakes commonly produce tsunamis along the Levant coastline, and that analogous situations are present elsewhere in the Mediterranean, as well as along the California coast and in another regions with active faults near the coast.

We modeled three typical scenarios, and in light of the Sumatra experience, we examined the more likely severe magnitudes. This of course leads us toward the upper range of expected run-ups. The models show that sooner than five minutes after a strong earthquake produces an offshore slump, which occurs after close to a third of the large dst earthquakes, a 4- to 6-m run-up may flood part of the Syrian, Lebanese, and Israeli coasts. Tsunamis from remote earthquakes, however, arrive later and produce only 1- to 3-m run-ups, but are more regional in extent.

Online material: Tsunami modeling and reports.

[Articles] Unveiling the Sources of the Catastrophic 1456 Multiple Earthquake: Hints to an Unexplored Tectonic Mechanism in Southern Italy

Fracassi, U., Valensise, G. — 2007-06-11

We revisited data related to the 1456 seismic crisis, the largest earthquake to have ever occurred in peninsular Italy, in search of its causative source(s). Data about this earthquake consist solely of historical reports and their intensity assessment.

Because of the age of this multiple earthquake, the scarcity and sparseness of the data, and the unusually large damage area, no previous studies have attempted to attribute the 1456 events to specific faults. Existing analytical methods to identify a likely source from intensity data also proved inappropriate for such a sparse dataset, since historical evidence suggests that the cumulative damage pattern contains at least three widely separated events.

We subdivided the 1456 damage pattern into three independent mesoseismal areas; each of these areas falls onto east–west tectonic trends previously identified and marked by deep (>10 km) right-lateral slip earthquakes. Based on this evidence we propose (1) that the 1456 events were generated by individual segments of regional east–west structures and are evidence of a seismogenic style that involves oblique dextral reactivation of east–west lower crustal faults; (2) that each event may have triggered subsequent but relatively distant events in a cascade fashion, as suggested by historical accounts; hence (3) that the 1456 sequence reveals a fundamental but unexplored mechanism of tectonic deformation and seismic release in southern Italy. This style dominates the region that lies between the northwest–southeast system of large extensional faults straddling the crest of the southern Apennines and the buried outer front of the chain.

Although the quality of the available information concerning the 1456 earthquake is naturally limited, we show that the overlap of the damage distribution, the orientation and characteristics of regional tectonic structures, the seismicity patterns, and the focal mechanisms all concur with our interpretations and would be difficult to justify otherwise.

[Articles] 12,000-Year-Long Record of 10 to 13 Paleoearthquakes on the Yammouneh Fault, Levant Fault System, Lebanon

Daeron, M., Klinger, Y., Tapponnier, P., Elias, A., Jacques, E., Sursock, A. — 2007-06-11

We present results of the first paleoseismic study of the Yammoûneh fault, the main on-land segment of the Levant fault system within the Lebanese restraining bend. A trench was excavated in the Yammoûneh paleolake, where the fault cuts through finely laminated sequences of marls and clays. First-order variations throughout this outstanding stratigraphic record appear to reflect climate change at centennial and millennial scales. The lake beds are offset and deformed in a 2-m- wide zone coinciding with the mapped fault trace. Ten to thirteen events are identified, extending back more than ~12 kyr. Reliable age bounds on seven of these events constrain the mean seismic return time to 1127 ± 135 yr between ~12 ka and ~6.4 ka, implying that this fault slips in infrequent but large (M ~7.5) earthquakes. Our results also provide conclusive evidence that the latest event at this site was the great a.d. 1202 historical earthquake, and suggest that the Yammoûneh fault might have been the source of a less well-known event circa a.d. 350. These findings, combined with previous paleoseismic data from the Zebadani valley, imply that the parallel faults bounding the Beqaa release strain in events with comparable recurrence intervals but significantly different magnitudes. Our results contribute to document the clustering of large events on the Levant fault into centennial episodes, such as that during the eleventh through twelfth centuries, separated by millennial periods of quiescence, and raise the possibility of a M >7 event occurring on the Yammoûneh fault in the coming century. Such a scenario should be taken into account in regional seismic-hazard assessments and planned for accordingly.

[Articles] Tectonic Processes in the Jan Mayen Fracture Zone Based on Earthquake Occurrence and Bathymetry

2007-06-11

Jan Mayen is an active volcanic island situated along the mid-Atlantic Ridge north of Iceland. It is closely connected with the geodynamic processes associated with the interaction between the Jan Mayen Fracture Zone (jmfz) and the slowly spreading Kolbeinsey and Mohns Ridges. Despite the significant tectonic activity expressed by the frequent occurrence of medium to large earthquakes, detailed correlation between individual events and the causative faults along the jmfz has been lacking. Recently acquired detailed bathymetric data in the vicinity of Jan Mayen has allowed us to document such correlation for the first time. The earthquake of 14 April 2004 (M_w 6), which occurred along the jmfz, was studied in detail and correlated with the bathymetry. Locations of aftershocks within the first 12 hours after the mainshock outline a 10-km-long fault plane. Interactions between various fault systems are demonstrated through locations of later aftershocks, which indicate that supposedly normal fault structures to the north of the ruptured fault, in the Jan Mayen Platform, have been reactivated. Correlation of the waveforms shows that events located on these structures are significantly different from activity at neighboring structures. Coulomb stress modeling gives an explanation to the locations of the aftershocks but cannot reveal any information about their mechanisms.

[Articles] Aftershock Detection Thresholds as a Function of Time: Results from the ANZA Seismic Network following the 31 October 2001 ML 5.1 Anza, California, Earthquake

Kilb, D., Martynov, V. G., Vernon, F. L. — 2007-06-11

We examine aftershock detectability thresholds for events in the initial part of the 31 October 2001, M_L 5.1 sequence in southern California. This sequence occurred directly below the broadband anza seismic network, which recorded continuous waveform data at 13 azimuthally well-distributed stations within the study region (seven had epicentral distances < 20 km). Of the 608 aftershocks (0 < M_L < ~2.8) in the initial 2 hr of this sequence, the first five aftershocks recorded were only identifiable at stations within 30 km after applying a high-pass filter. Using a cluster (radius ≤ 1.1 km) of 200 representative aftershocks, we track the maximum seismogram amplitude versus earthquake magnitude. This relationship helps us quantify the visibility of aftershocks within the mainshock coda and assess our detection capabilities. We estimate that detectable aftershocks within the mainshock coda include (1) those over magnitude ~3 that are within 15 km of the network centroid that occur 12 sec or more into the sequence, and (2) those over magnitude ~2 that are within 30 km of the centroid of the network that occur 60 sec or more into the sequence. We find a lack of large aftershocks in this sequence. The largest aftershock (M_L ~2.8) is substantially smaller than the mainshock (M_L 5.1). We suggest this relatively large-magnitude differential is dictated by a combination of factors that includes complexity of the San Jacinto fault system and the lack of large earthquakes in the region in the past ~20 years.

Online material: Quicktime movies juxtaposing a 3.2 aftershock in the coda of a 5.1 mainshock.

[Articles] Wavelet Transform Methods for Azimuth Estimation in Local Three-Component Seismograms

Galiana-Merino, J. J., Rosa-Herranz, J., Jauregui, P., Molina, S., Giner, J. — 2007-06-11

We have developed a new method for estimating azimuth in local three- component seismograms through wavelet analysis. The proposed process proceeds in three stages. First, the seismogram is filtered through a wavelet packet approach. Second, the first arrival is determined through a P picker similar to a short-term- average long-term-average scheme that is applied on the wavelet domain. Finally, an adaptive-length window around the pick is selected and used for determining the azimuth, using the property of linear polarization of the first arrival.

The proposed method has been applied to three-component short-period seismograms for local earthquakes recorded by the seismic network of Alicante province in southeastern Spain. The locations of these events were previously obtained through the software HYPO71PC (Lee and Valdes, 1989) and the seismograms recorded by four analog stations of vertical component distributed within the province of Alicante. The results obtained by the wavelet-based algorithm have been compared with the azimuth angles obtained through the results from the location software. The comparison indicates that the proposed algorithm can determine the azimuth of the analyzed events to within a mean bias of 4.5°.

[Articles] Regional Travel-Time Uncertainty and Seismic Location Improvement Using a Three-Dimensional a priori Velocity Model

Flanagan, M. P., Myers, S. C., Koper, K. D. — 2007-06-11

We demonstrate our ability to improve regional travel-time prediction and seismic event location accuracy using an a priori 3D velocity model of Western Eurasia and North Africa (WENA1.0). Travel-time residuals are assessed relative to the iasp91 model for approximately 6000 Pg, Pn, and P arrivals, from seismic events having 2 epicenter accuracy between 1 km and 25 km (ground truth 1 [GT1] and GT25, respectively), recorded at 39 stations throughout the model region. Ray paths range in length between 0° and 40° (local, regional, and near teleseismic) providing depth sounding that spans the crust and upper mantle. The dataset also provides representative geographic sampling across Eurasia and North Africa including aseismic areas. The WENA1.0 model markedly improves travel-time predictions for most stations with an average variance reduction of 29% for all ray paths from the GT25 events; when we consider GT5 and better events alone, the variance reduction is 49%. For location tests we use 196 geographically distributed GT5 and better events. In 134 cases (68% of the events), locations are improved, and average mislocation is reduced from 24.9 km to 17.7 km. We develop a travel-time uncertainty model that is used to calculate location coverage ellipses. The coverage ellipses for WENA1.0 are validated to be representative of epicenter error and are smaller than those for iasp91 by 37%. We conclude that a priori models are directly applicable where data coverage limits tomographic and empirical approaches, and the development of the uncertainty model enables merging of a priori and data-driven approaches using Bayesian techniques.

Online material: Correction surfaces and histograms of travel-time residuals for 40 stations.

[Articles] Seismic Velocity Structure and Seismotectonics of the Eastern San Francisco Bay Region, California

Hardebeck, J. L., Michael, A. J., Brocher, T. M. — 2007-06-11

The Hayward Fault System is considered the most likely fault system in the San Francisco Bay Area, California, to produce a major earthquake in the next 30 years. To better understand this fault system, we use microseismicity to study its structure and kinematics. We present a new 3D seismic-velocity model for the eastern San Francisco Bay region, using microseismicity and controlled sources, which reveals a ~10% velocity contrast across the Hayward fault in the upper 10 km, with higher velocity in the Franciscan Complex to the west relative to the Great Valley Sequence to the east. This contrast is imaged more sharply in our localized model than in previous regional-scale models. Thick Cenozoic sedimentary basins, such as the Livermore basin, which may experience particularly strong shaking during an earthquake, are imaged in the model.

The accurate earthquake locations and focal mechanisms obtained by using the 3D model allow us to study fault complexity and its implications for seismic hazard. The relocated hypocenters along the Hayward Fault in general are consistent with a near-vertical or steeply east-dipping fault zone. The southern Hayward fault merges smoothly with the Calaveras fault at depth, suggesting that large earthquakes may rupture across both faults. The use of the 3D velocity model reveals that most earthquakes along the Hayward fault have near-vertical strike-slip focal mechanisms, consistent with the large-scale orientation and sense of slip of the fault, with no evidence for zones of complex fracturing acting as barriers to earthquake rupture.

Online material: Velocity model validation experiments and additional seismicity plots.

[Articles] A Regional Ground-Motion Excitation/Attenuation Model for the San Francisco Region

Malagnini, L., Mayeda, K., Uhrhammer, R., Akinci, A., Herrmann, R. B. — 2007-06-11

By using small-to-moderate earthquakes located within ~200 km of San Francisco, we characterize the scaling of the ground motions for frequencies ranging between 0.25 and 20 Hz, obtaining results for geometric spreading, Q(f), and site parameters using the methods of Mayeda et al. (2005) and Malagnini et al. (2004). The results of the analysis show that, throughout the Bay Area, the average regional attenuation of the ground motion can be modeled with a bilinear geometric spreading function with a 30-km crossover distance, coupled to an anelastic function exp(–fr/Q(f), where: Q(f) = 180 f^0.42. A body-wave geometric spreading, g(r) = r^–1.0, is used at short hypocentral distances (r < 30 km), whereas g(r) = r^–0.6 fits the attenuation of the spectral amplitudes at hypocentral distances beyond the crossover.

The frequency-dependent site effects at twelve of the Berkeley Digital Seismic Network stations were evaluated in an absolute sense using coda-derived source spectra. Our results show the following. (1) The absolute site response for frequencies ranging between 0.3 Hz and 2.0 Hz correlate with independent estimates of the local magnitude residuals (M_L) for each of the stations. (2) Moment magnitudes (M_w) derived from our path and site-corrected spectra are in excellent agreement with those independently derived using full-waveform modeling as well as coda-derived source spectra. (3) We use our weak-motion-based relationships to predict motions regionwide for the Loma Prieta earthquake, well above the maximum magnitude spanned by our data set, on a completely different set of stations. Results compare well with measurements taken at specific National Earthquake Hazards Reduction Program site classes. (4) An empirical, magnitude-dependent scaling was necessary for the Brune stress parameter to match the large-magnitude spectral accelerations and peak ground velocities with our weak-motion-based model.

Online material: Tables of peak ground acceleration, peak ground velocity, and pseudo-spectral acceleration at 0.3 sec, 1.0 sec, and 3.0 sec.

[Articles] Anomalous Propagation of Long-Period Ground Motions Recorded in Tokyo during the 23 October 2004 Mw 6.6 Niigata-ken Chuetsu, Japan, Earthquake

Furumura, T., Hayakawa, T. — 2007-06-11

Unusually large (>5 cm) and prolonged shaking associated with long- period ground motions at periods of about 7 sec were observed in central Tokyo during the M_w 6.6 Niigata-ken Chuetsu earthquake of 23 October 2004. The long- period ground motions caused significant resonance in high-rise buildings of about 70 floors in height. Thus, it is an urgent matter to understand the development and amplification properties of long-period ground motions in Tokyo associated with large earthquakes.

In this study, we use numerous waveform records from 585 stations in a nationwide accelerometer network (K-NET, KiK-net) and 495 intensity meters in the area around Tokyo. The data reveal that the long-period ground motion is characterized in most part by a surface, Rayleigh wave generated at the northern edge of Kanto basin, and the surface wave is developed as propagating through a thick cover of sediments (>3000–4000 m) that overlies rigid bedrock.

To complement the observational data, we conducted a large-scale computer simulation of seismic-wave propagation by employing the Earth Simulator supercomputer with a detailed source-slip model and a high-resolution 3D sedimentary structural model of central Japan. The results of the computer simulation demonstrate that the anomalously prolonged ground shaking of the long-period signal recorded in the center of Tokyo occurred because of the stagnation of seismic energy resulting from the multipathing and focusing of Rayleigh waves toward the bottom of the Kanto basin from surrounding mountain regions with interaction to the 3D basin structure.

Online material: Animations of long-period ground motions in Tokyo.

[Articles] Sensitivity of Ground-Motion Simulations to Earthquake Source Parameters: A Case Study for Istanbul, Turkey

Sorensen, M. B., Pulido, N., Atakan, K. — 2007-06-11

Following the disastrous earthquakes in Izmit and Düzce along the North Anatolian fault in 1999, the earthquake hazard in the Istanbul area became a great concern. In this study we simulate strong ground motions caused by a scenario earthquake (M 7.5) in the Marmara Sea, and investigate the effect of varying the input parameters on the broadband frequency ground motion. Simulations are based on a multiasperity source model that involves the combined rupture of the North Anatolian fault segments beneath the Marmara Sea. We use a hybrid model combining a deterministic simulation of the low frequencies (0.1–1.0 Hz) with a semistochastic simulation of the high frequencies (1.0–10.0 Hz). Computation at each frequency range is performed separately and the total ground motion is combined in the time domain. Computations are linear and are performed at bedrock level, thereby not taking any effect of local geological conditions into account. We calculate a total of 17 earthquake scenarios corresponding to different source and attenuation parameters to study their effect on the ground motion. The most significant parameters in terms of ground-shaking level are the rise time, rupture velocity, rupture initiation point, and stress drop. The largest variability of strong ground motions is observed in regions adjacent to asperities and is associated with frequencies higher than 5 Hz. For lower frequencies our simulated velocity spectra within the Istanbul area are fairly stable among scenarios. The average standard deviations of all ground-motion measures are less than 35% of the mean.

Online material: Figures of peak ground acceleration and peak ground velocity and their differences to the reference scenario values.

[Articles] Assessment of a Nonlinear Dynamic Rupture Inversion Technique Applied to a Synthetic Earthquake

Corish, S. M., Bradley, C. R., Olsen, K. B. — 2007-06-11

Dynamic rupture inversion is a powerful tool for learning why and how faults fail, but much more work has been done in developing inversion methods than evaluating how well these methods work. This study examines how well a nonlinear rupture inversion method recovers a set of known dynamic rupture parameters on a synthetic fault based on the 2000 western Tottori, Japan earthquake (M_w 6.6). Rupture evolution on the fault is governed by a slip-weakening friction law. A direct-search method known as the neighborhood algorithm (Sambridge, 1999) is used to find optimal values of both the initial stress distribution and the slip-weakening distance on the fault, based on misfit values between known and predicted strong-motion displacement records. The yield stress and frictional sliding stress on the fault are held constant. A statistical assessment of the results shows that, for this test case, the inversion succeeds in locating all parameters to within ±14% of their true values. With the model configuration used in this study, the parameters located in the central rupture area are better resolved than the parameters located at the sides and bottom of the fault. In addition, a positive linear correlation between the mean initial stress and the slip-weakening distance is identified. The investigation confirms that dynamic rupture inversion is useful for determining rupture parameters on the fault, but that intrinsic trade-offs and poor resolution of some parameters limit the amount of information that can be unambiguously inferred from the results. In addition, this study demonstrates that using a statistical approach to assess nonlinear inversion results shows how sensitive the misfit measure is to the various parameters, and allows a level of confidence to be attached to the output parameter values.

[Articles] A Constitutive Criterion for the Fault: Modified Velocity-Weakening Law

Miguel, L. F. F., Riera, J. D. — 2007-06-11

Expanding on a previous article, constitutive laws for solid friction are examined jointly with available experimental results. The models are evaluated by means of numerical dynamic analysis of two sliding blocks simulating adjacent fault sections. Effective constitutive laws are determined as relations between mean values of the relevant variables on a selected area of the sliding surfaces. The material is initially assumed elastic and homogeneous and the influence of nonhomogeneity evaluated next by modeling mass density, Young’s modulus, and friction coefficient as correlated random fields. The effect of fractures in rock close to the fault is also numerically assessed. Finally, the influence of rupture of protrusions (microasperities) between the sliding surfaces is analyzed. The influence of size of the averaging interface area on the parameters of the effective constitutive law is then obtained by means of Monte Carlo simulation. When the rock regions adjacent to the fault are assumed to be linearly elastic and homogeneous or nonhomogeneous no size effect is observed. On the other hand, when the friction coefficient is characterized by a random field, a size effect is detected. Fracture occurrence in the region surrounding the fault does not cause significant alteration of the macroconstitutive law, producing only minor perturbations of the mean law determined without fracture, but it should also introduce high-frequency slave vibrations. Finally, a macroconstitutive law that takes into account the shear rupture of microasperities on the sliding surfaces is suggested. The proposed modified velocity-weakening law, constitutes a more general and flexible constitutive law.

[Articles] Wave Gradiometry in the Time Domain

Langston, C. A. — 2007-06-11

A time-domain approach for solving for the change in geometrical spreading and horizontal wave slowness in wave gradiometry is presented based on the use of the analytic signal. The horizontal displacement gradient of a wave is linearly related to the displacement and its time derivative. The coefficients of this relationship give the change of geometrical spreading, the change in radiation pattern, and horizontal slowness. The new time-domain technique incorporates estimates of the instantaneous amplitude and frequency of the three time series to solve uniquely for the wave-field coefficients. The analysis is simpler and more suited to fast array processing of displacement gradient data sets compared with a spectral ratio method.

[Articles] Damping Correction Factors for Horizontal Ground-Motion Response Spectra

Cameron, W. I., Green, R. A. — 2007-06-11

Damping correction factors (dcfs) are used to adjust response spectral values corresponding to damping 5% of critical to other damping levels. Trends in the analytical response of viscously damped, linear-elastic sdof oscillators subjected to finite-duration, sinusoidal base excitations show that dcfs depend on both the frequency and duration of the ground motion, where the latter becomes significantly less influential as damping increases. These analytical trends, in conjunction with correlations relating duration and frequency content to earthquake magnitude, site- to-source distance, site classification, and tectonic setting, are used to explain/study observations in dcfs computed from a large ground-motion database for the central- eastern and western United States. For ≥ 2%, the dcfs proposed by the authors depend on earthquake magnitude, site classification, and tectonic setting, all of which significantly influence the frequency content of ground motions. For = 1%, the dcfs proposed by the authors additionally depend on site-to-source distance, which significantly influences the duration of ground motion. In comparison with the dcfs proposed by the authors, commonly used and recently proposed dcfs were shown to be both too low and too high, depending on the relation, period range, damping ratio, earthquake magnitude, site classification, and tectonic setting. Additionally, the dcf relations proposed in the literature for ground motions exhibiting near-fault effects should not be used for < 5%, and will likely be significantly too high for periods close to that of the near-fault velocity pulse(s) for ≥ 5%.

[Articles] Liquefaction, Ground Oscillation, and Soil Deformation at the Wildlife Array, California

Holzer, T. L., Youd, T. L. — 2007-06-11

Excess pore-water pressure and liquefaction at the Wildlife Liquefaction Array in 1987 were caused by deformation associated with both high-frequency strong ground motion and 5.5-second-period Love waves. The Love waves produced large (~1.5%) cyclic shear strains well after the stronger high-frequency ground motion abated. These cyclic strains generated approximately from 13 to 35% of the excess pore-water pressure in the liquefied layer and caused excess pore-water pressures ultimately to reach effective overburden stress. The deformation associated with the Love waves explains the "postearthquake" increase of pore-water pressure that was recorded at the array. This explanation suggests that conventional methods for predicting liquefaction based on peak ground acceleration are incomplete and may need to consider cyclic strains associated with long-period surface waves. A postearthquake survey of an inclinometer casing indicated permanent shear strain associated with lateral spreading primarily occurred in the upper part of the liquefied layer. Comparison of cone penetration test soundings conducted after the earthquake with pre-earthquake soundings suggests sleeve friction increased. Natural lateral variability of the liquefied layer obscured changes in tip resistance despite a ~1% reduction in volume. The large oscillatory motion associated with surface waves explains ground oscillation that has been reported at some liquefaction sites during earthquakes.

[Articles] Corner Frequency Scaling of Regional Seismic Phases for Underground Nuclear Explosions at the Nevada Test Site

Fisk, M. D. — 2007-06-11

Model fits are used to estimate source spectral corner frequencies of regional seismic phases Pn, Pg, and Lg from underground nuclear explosions (unes) at the Nevada Test Site (nts), based on recordings by Lawrence Livermore National Laboratory (llnl) nts Network stations. The spectra are corrected for instrument response, distance, and station effects, and then network averaged. Explicit information regarding geophysical working-point properties, depth of burial, and some announced yields for nts explosions (Springer et al., 2002) is utilized for the spectral modeling. Scaling of Pn, Pg, and Lg corner frequencies with source yield is examined for sets of explosions at Yucca Flat, Pahute Mesa, and Rainier Mesa with similar material properties (medium type, density, velocities, and gas porosity). The unes are limited to those in media with gas porosity of 10% or less so that a source spectral model with a rolloff of f^–2 is applicable. Model comparisons are also examined for regional P/Lg discriminants. A key result is that Lg corner frequencies for nts unes exhibit similar scaling with source size as for P waves, but shifted lower, analogous to observations by Fisk (2006) for the Lop Nor, Semipalatinsk, and Novaya Zemlya test sites. This appears to corroborate an important effect, that is, that major contributions to S-wave generation by explosions occur near the source with a similar length scale, comparable to the elastic radius, as for P waves from explosions. Although the explicit physical mechanism is not yet understood, the implications are important regarding where and perhaps how S waves are predominantly generated by explosions. A related, key result is that the increasing separation of P/Lg at higher frequencies between nts explosions and earthquakes has a consistent model-based explanation, as at other nuclear test sites, in terms of the difference between explosion P and S corner frequencies.

[Articles] Identification of Delay-Fired Mining Explosions Using Seismic Arrays: Application to the PDAR Array in Wyoming, USA

Arrowsmith, S. J., Hedlin, M. A. H., Arrowsmith, M. D., Stump, B. W. — 2007-06-11

We extend a time-frequency discrimination algorithm, developed in an earlier article (Arrowsmith et al., 2006), for application to seismic-array data. Spectrograms evaluated at each component of an array are stacked and then converted into binary form for computation of discriminants. Because noise can bias the discriminants, we develop a procedure for removing the effect of noise on the discriminants. The binary spectrograms are randomized where the spectral amplitude of the signal is similar to the mean spectral amplitude of the pre-event noise at that frequency. The formulism of Arrowsmith et al. (2006) is further extended by modifying the objective function used to optimize the values of input parameters and by removing high-frequency and low-frequency spectral content. We apply the method to a dataset of regional recordings of earthquakes and delay-fired mine blasts recorded at the Pinedale seismic array in Wyoming. Our results show that the utilization of array data improves the success rate for source identification. Furthermore, we find that incorporating the noise-correction procedure increases the separation between earthquakes and cast overburden blasts (the largest type of delay-fired mine blasts). In total, the algorithm successfully identifies 97.4% of the events (74 of a total of 76 events, which comprise earthquakes and cast overburden blasts).

[Articles] Imaging the Three-Dimensional Crust of the Korean Peninsula by Joint Inversion of Surface-Wave Dispersion and Teleseismic Receiver Functions

Yoo, H. J., Herrmann, R. B., Cho, K. H., Lee, K. — 2007-06-11

A detailed study of the 3D variation of shear-wave velocities in the southern part of the Korean Peninsula is made by combining high-frequency surface- wave tomography results of Cho et al. (2006b) with teleseismic P-wave receiver functions at 80 locations on the peninsula. Receiver functions were derived from high-gain acceleration, short-period, and broadband digital data streams of the Korea Meteorological Administration (kma) and Korean Institute for Geosciences and Mineral Resources (kigam) networks. Vertical cross sections trace the lateral variation in the depth to the Moho, the variation of low velocities near the surface, and the variable thickness of the transition from surface velocities to midcrustal velocities. The derived crustal structure provides new insights on the evolution of the Korean crust.

[Articles] Spatial Distribution of Coda Q in South Korea

Yun, S., Lee, W. S., Lee, K., Noh, M. H. — 2007-06-11

Analyzing spectral attenuation of coda waves, we first mapped coda Q- values in the crust over the whole of South Korea. During the period from 1995 to 2004, we selected 574 north–south-component seismograms with epicentral distances less than 100 km from 328 earthquakes with magnitudes between 1.4 and 5.2. We estimated coda Q-values using the single isotropic scattering model at center frequencies of 1.5, 3, 6, 9, 12, 15, and 18 Hz, and found significant spatial variation over all frequency ranges and strong frequency dependence of coda Q in the region. The estimated coda Q-value at 1 Hz (Q₀) and the -value ranges are 80–300 and 0.4–1.1, respectively. The values have strong correlation with the regional geology in the Korean peninsula. Low Q₀-values are mainly obtained in the regions comprising sedimentary strata in southeastern South Korea, whereas granite regions in the northern part of South Korea show high Q₀. The Q₀-values in the study area agree well with those of the eastern China and Kyushu, western Japan. Furthermore, our -values are also in good agreement with those of Japan.

Online material: Locations of seismometers, list of event-station pairs, and linear regression of Q_C measurements.

[Short Notes] Holocene Slip Rate for the Western Segment of the Castle Mountain Fault, Alaska

Willis, J. B., Haeussler, P. J., Bruhn, R. L., Willis, G. C. — 2007-06-11

The western segment of the Castle Mountain fault poses a significant seismic hazard to the most populated region of south-central Alaska. We identify a previously unrecognized margin of a postglacial outwash channel that is offset right laterally 36 ± 4 m across the western segment of the Castle Mountain fault. This offset occurred after glaciers withdrew from the lowland 11,300–15,380 cal yr b.p. and after outwash channel margins were cut and stabilized 11,210–13,470 cal yr b.p. Using these ages and the measured separation, we obtain a maximum slip rate of 3.0 ± 0.6 mm yr^–1 and a minimum slip rate of 2.8 ± 0.7 mm yr^–1. These are the first lateral slip rates for the Castle Mountain fault established by a field measurement. Based on timing of the most recent earthquake, 670 ± 60 yr b.p., the Castle Mountain fault could have accumulated an average single-event slip of about 1.9 m (extremes range from 1.3 to 2.6 m). The fault consists of two segments; a surface-rupturing earthquake likely will be limited to the 62-km-long western segment. Area-magnitude regression calculations suggest that such an earthquake on the western Castle Mountain fault would have a moment magnitude of 6.9 to 7.3.

[Short Notes] Wave Propagation in Irregularly Layered Elastic Models: A Boundary Element Approach with a Global Reflection/Transmission Matrix Propagator

Ge, Z., Chen, X. — 2007-06-11

A direct boundary element method that uses the full-space Green’s function is proposed for calculating elastic wave propagation in two-dimensional irregularly stratified models. The global matrix equation becomes larger as the number of layers increases. These equations are usually solved by improved block Gaussian elimination, conjugate gradient algorithms, or other approaches based on different approximations. In this article, we adopt the global generalized reflection/transmission matrix method (Chen, 1990, 1995, 1996) to solve this problem. This method can prevent excessive requirement of both computer memory and CPU time. The method is validated by comparing its results with those obtained using the finite- difference method.

[Errata] Earthquake Ground-Motion Prediction Equations for Eastern North America

Atkinson, G. M., Boore, D. M. — 2007-06-11