Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “The Seismological Record”

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.

535 records · Page 30Linked to original sources

Uncertainties in intensity-based earthquake magnitude estimates

Estimating the magnitude of historical earthquakes is crucial for assessing seismic hazard. Magnitudes of early‐instrumental earthquakes can be inferred using a combination of instrumental records, field observations, and the observed distribution of shaking intensity determined from macroseismic observations. For earthquakes before 1900, shaking intensity distributions often provide the only information to constrain earthquake magnitude. Considerable effort has been made to develop methods to estimate the magnitude of moderate‐to‐large historical earthquakes using shaking intensities derived from macroseismic data. In this study, we consider earthquakes in California with known instrumental magnitudes to explore uncertainties in estimating the magnitude of historical earthquakes from intensity information alone. We use three California‐specific intensity prediction equations (IPEs) and an IPE based on a global ground‐motion model (GMM) to determine optimum intensity‐based magnitudes for 33 moderate‐to‐large California earthquakes between 1979 and 2021. Intensity‐based magnitudes are close to instrumental magnitudes on average. However, intensity‐based magnitudes for individual events differ by as much as 2.2 magnitude units from instrumental magnitudes. This result reflects the weak dependence of ground motions and shaking intensities on moment magnitude and their strong dependence on stress drop. Considering the intensity distributions of the 1906 San Francisco and 1989 Loma Prieta earthquakes, we show that information that could constrain rupture length is discarded when considering only the 2D decay of intensity with distance. We also show that ground‐motion intensity conversion equations used in a GMM‐based approach may cause a systematic overestimation of large historical earthquake magnitudes. This study underscores both the reducible and potentially irreducible uncertainties associated with using intensity data to estimate magnitudes of historical earthquakes using IPEs and highlights the value of using additional information to constrain rupture dimensions. Using intensity observations alone, moment magnitude uncertainties are typically on the order of a full unit.

California↗

Directivity of M 3.1 earthquake near Anza, California and the effect on peak ground motion

We show the effect of rupture directivity on peak ground‐motion values for a moderate magnitude event at Anza, California, and neighboring stations at the Imperial Valley. The event was located near Borrego Springs on the west side of the Salton Sea and was well recorded at broadband stations near Anza, California, and at stations on the west side of the Imperial Valley. After correcting for regional attenuation, an anomalously large residual in peak motion was observed at station ERR just to the southeast of the epicenter. Using the algorithm from Boatwright (2007) , peak motions from the regional seismic networks in southern California were inverted to determine directivity, which was to the southeast along the trend of the San Jacinto fault toward station ERR. This algorithm uses peak values compiled for the ShakeMap system mostly at regional distances. It does not capture the main features of the source time function (STF) predicted by directivity. Consequently, we determined the second‐degree moments for this earthquake, which confirmed that station ERR has a shorter and higher STF compared to stations to the northwest suggesting rupture propagated to the southeast. The azimuthal distribution of local stations is sparse, but nevertheless the largest amplitudes (such as at station ERR) correlate well with the maximum in the radiation pattern and smaller values with the minima, which is the radiation pattern for SH plus the effect of directivity. Using the data from the analysis of the second‐degree moments, the characteristic length of the fault is 0.58 km, assuming an idealized unilateral extended rupture with a rupture time of 0.09 s. This yields an apparent rupture velocity of 6.4    km / s "> 6.4 km / s for an idealized model, which is super shear. This value is model dependent and would change if, for example, the rupture was bilateral. Although this value is even greater than the P ‐wave velocity, it supports the idea that the rupture velocity is super shear and would enhance the correlation between the peak motions and the radiation pattern.

California↗

Ground motion in Anchorage, Alaska, from the 2002 Denali fault earthquake: Site response and displacement pulses

Data from the 2002 Denali fault earthquake recorded at 26 sites in and near Anchorage, Alaska, show a number of systematic features important in studies of site response and in constructing long-period spectra for use in earthquake engineering. The data demonstrate that National Earthquake Hazards Reduction Program (NEHRP) site classes are a useful way of grouping stations according to site amplification. In general, the sites underlain by lower shear-wave velocities have higher amplification. The amplification on NEHRP class D sites exceeds a factor of 2 relative to an average of motions on class C sites. The amplifications are period dependent. They are in rough agreement with those from previous studies, but the new data show that the amplifications extend to at least 10 sec, periods longer than considered in previous studies. At periods longer than about 14 sec, all sites have motion of similar amplitude, and the ground displacements are similar in shape, polarization, and amplitude for all stations. The displacement ground motion is dominated by a series of four pulses, which are associated with the three subevents identified in inversion studies (the first pulse is composed of P waves from the first subevent). Most of the high-frequency ground motion is associated with the S waves from subevent 1. The pulses from subevents 1 and 2, with moment releases corresponding to M 7.1 and 7.0, are similar to the pulse of displacement radiated by the M 7.1 Hector Mine earthquake. The signature from the largest subevent (M 7.6) is more subdued than those from the first two subevents. The two largest pulses produce response spectra with peaks at a period of about 15 sec. The spectral shape at long periods is in good agreement with the recent 2003 NEHRP code spectra but is in poor agreement with the shape obtained from Eurocode 8.

Alaska↗

Waveform signatures of earthquakes located close to the subducted Gorda Plate interface

Complex seismic velocity structure near the earthquake source can affect rupture dynamics and strongly modify the seismic waveforms recorded near the fault. Fault‐zone waves are commonly observed in continental crustal settings but are less clear in subduction zones due to the spatial separation between seismic stations and the plate boundary fault. We observed anomalously long duration S waves from earthquake clusters located near the interface of the subducted Gorda plate north of the Mendocino triple junction. In contrast, earthquakes located just a few kilometers below each cluster show impulsive S waves. A nodal array experiment was conducted around the Northern California Seismic Network station KCT for two months to investigate the origin of the complex S waves. Beamforming analysis shows that the S waves contain three arrivals that have different horizontal slownesses, which we term S 1, S 2, and S coda. Similar analysis on P waves also show two arrivals with different horizontal slownesses, which we term P 1 and P 2. P 1 and S 1 have larger horizontal slowness than P 2 and S 2, respectively, indicating that the phase pairs are body waves with different ray paths. Building upon a seismic refraction profile, we construct 1D velocity models and test different thicknesses and V P / V S ratios for the subducted oceanic crust. The arrival times and relative slownesses of P 1/ P 2 and S 1/ S 2 phases indicate that they are the direct and the Moho reflected phases, respectively. Their properties are consistent with a crustal thickness of ∼6 km and a moderate V P / V S ratio (∼1.8). The S coda is more difficult to characterize but has a clear dominant frequency that likely reflects the near‐source velocity and attenuation structure. Our study indicates that waveforms from earthquakes near the interface of the subducted slab can be used to infer detailed structural information about the plate‐boundary zone at seismogenic depths.

California↗

On the characteristics of local geology and their influence on ground motions generated by the Loma Prieta earthquake in the San Francisco Bay region, California

Strong ground motions recorded at 34 sites in the San Francisco Bay region from the Loma Prieta earthquake show marked variations in characteristics dependent on crustal structure and local geological conditions. Peak horizontal acceleration and velocity inferred for sites underlain by “rock” generally occur on the transverse component of motion. They are consistently greater with lower attenuation rates than the corresponding mean value predicted by empirical curves based on previous strong-motion data. Theoretical amplitude distributions and synthetic seismograms calculated for 10-layer models suggest that “bedrock” motions were elevated due in part to the wide-angle reflection of S energy from the base of a relatively thin (25 km) continental crust in the region. Characteristics of geologic and geotechnical units as currently mapped for the San Francisco Bay region show that average ratios of peak horizontal acceleration, velocity and displacement increase with decreasing mean shear-wave velocity. Ratios of peak acceleration for sites on “soil” (alluvium, fill/Bay mud) are statistically larger than those for sites on “hard rock” (sandstone, shale, Franciscan Complex). Spectral ratios establish the existence of predominant site periods with peak amplifications near 15 for potentially damaging levels of ground motion at some sites underlain by alluvium and fill/bay mud. Average spectral amplifications inferred for vertical and the mean horizontal motion are, respectively, (1,1) for sites on the Franciscan Complex (KJf), (1.4, 1.5) for sites on Mesozoic and Tertiary rocks (TMzs), (2.1, 2.0) for sites on the Santa Clara Formation (QTs), (2.3, 2.9) for sites on alluvium (Qal), and (2.1, 4.0) for sites on fill/Bay mud (Qaf/Qhbm). These mean values are not statistically different at the 5% significance level from those inferred from previous low-strain data. Analyses suggest that soil amplification and reflected crustal shear energy were major contributors to levels of ground motion sufficient to cause damage to vulnerable structures at distances near 100 km in the cities of San Francisco and Oakland.

California↗

Short baseline variations in site response and wave-propagation effects and their structural causes: Four examples in and around the santa clara valley, California

Ground motion records of local and regional events from a portable array are used to investigate the structural causes of variations in ground motion over distances of a few hundred meters to a few kilometers in the sedimentary basin environment of the Santa Clara Valley, California, and its margins. Arrays of portable seismic stations are used to target four study areas with different ground motion patterns: (1) an edge of the alluvial basin extending up onto a marginal ridge (Blossom Hill), (2) a Cenozoic basin with a nearly flat bottom (Cupertino Basin), (3) a long, narrow Cenozoic basin with a steep V profile (Evergreen Basin), and (4) a line perpendicular to the trace of the Hayward fault. Average peak velocities on Blossom Hill from local earthquakes are a factor of 2.5 times higher than nearby valley sites. Three-dimensional (3D) modeling is used to conclude that the majority of the amplification is due to lower shear-wave velocities along a local fault zone (Shannon–Berrocal). Site amplification over the Cupertino Basin in the frequency band 0.5–4 Hz is generally low (less than 2.0 relative to a Mesozoic rock site) and spatially uniform. This response is attributed to the shallow, flat-bottomed shape of the basin and the uniform, flat-laying sedimentary fill. In contrast, site amplification in the Evergreen Basin generally exceeds 3.0 and is attributed to the deep, V-shaped geometry of the basin and younger sedimentary fill. 3D waveform modeling shows the elongated shape of the Evergreen Basin causes more efficient trapping of long-period waves for sources along the long axis of the basin. A low-velocity zone is postulated along the Hayward fault with a width between 100 and 200 m, based on elevated site response along the fault trace and 4.5-Hz fault zone guided waves on the horizontal components of stations near the fault.

California↗

Site response for urban Los Angeles using aftershocks of the Northridge earthquake

Ground-motion records from aftershocks of the 1994 Northridge earthquake are used to estimate site response in the urban Los Angeles area. Over 1300 shear-wave records from 61 sources and 90 sites are used in a linear inversion for source and site-response spectra. The methodology makes no assumptions about the shape of the source spectrum. To obtain a stable unique inverse, a Q model and geometrical spreading factor are assumed. In addition, the site response at a hard-rock site is constrained to be approximately 1.0 with a kappa of 0.02. The site-response spectra compare favorably with the results of previous and on-going investigations in Los Angeles. A couple of first-order effects are lower site response in the surrounding mountains, dominated by Mesozoic and Tertiary rocks, and higher values in the San Fernando and Los Angeles Basins, containing surficial Pleistocene and Holocene alluvial deposits. Results show good correlation of high site-response spectral values with localized areas of severe damage (Interstate 10 collapse, Sherman Oaks, Northridge, Interstate 5/14 collapse). However, widespread trends in site response across the sedimentary basins are not obvious. The data suggest that site responses are lower near the southern margin of the San Fernando Valley for sources to the north, due to north-dipping sedimentary structures. But the general pattern of site response is characterized by high variability on length scales less than a kilometer. Variations of a factor of 2 in site response are observed over the length scale of 200 m and for the same surficial geologic unit. For some of the alluvial basin sites, surface-wave generation is a significant contributor to elevated site response at lower frequencies, below 2 Hz. The total damage pattern for the Northridge earthquake is influenced by strong source directivity to the north and strong local site effects. The correlation of weak-motion site-response estimates with areas of significant damage demonstrates the value of these field measurements in future urban planning and in the reduction of seismic risk in urban areas.

California↗

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗

The Alaska earthquake, March 27, 1964: Lessons and conclusions

One of the greatest earthquakes of all time struck south-central Alaska on March 27, 1964. Strong motion lasted longer than for most recorded earthquakes, and more land surface was dislocated, vertically and horizontally, than by any known previous temblor. Never before were so many effects on earth processes and on the works of man available for study by scientists and engineers over so great an area. The seismic vibrations, which directly or indirectly caused most of the damage, were but surface manifestations of a great geologic event-the dislocation of a huge segment of the crust along a deeply buried fault whose nature and even exact location are still subjects for speculation. Not only was the land surface tilted by the great tectonic event beneath it, with resultant seismic sea waves that traversed the entire Pacific, but an enormous mass of land and sea floor moved several tens of feet horizontally toward the Gulf of Alaska. Downslope mass movements of rock, earth, and snow were initiated. Subaqueous slides along lake shores and seacoasts, near-horizontal movements of mobilized soil (“landspreading”), and giant translatory slides in sensitive clay did the most damage and provided the most new knowledge as to the origin, mechanics, and possible means of control or avoidance of such movements. The slopes of most of the deltas that slid in 1964, and that produced destructive local waves, are still as steep or steeper than they were before the earthquake and hence would be unstable or metastable in the event of another great earthquake. Rockslide avalanches provided new evidence that such masses may travel on cushions of compressed air, but a widely held theory that glaciers surge after an earthquake has not been substantiated. Innumerable ground fissures, many of them marked by copious emissions of water, caused much damage in towns and along transportation routes. Vibration also consolidated loose granular materials. In some coastal areas, local subsidence was superimposed on regional tectonic subsidence to heighten the flooding damage. Ground and surface waters were measurably affected by the earthquake, not only in Alaska but throughout the world. Expectably, local geologic conditions largely controlled the extent of structural damage, whether caused directly by seismic vibrations or by secondary effects such as those just described. Intensity was greatest in areas underlain by thick saturated unconsolidated deposits, least on indurated bedrock or permanently frozen ground, and intermediate on coarse well-drained gravel, on morainal deposits, or on moderately indurated sedimentary rocks. Local and even regional geology also controlled the distribution and extent of the earthquake's effects on hydrologic systems. In the conterminous United States, for example, seiches in wells and bodies of surface water were controlled by geologic structures of regional dimension. Devastating as the earthquake was, it had many long-term beneficial effects. Many of these were socioeconomic or engineering in nature; others were of scientific value. Much new and corroborative basic geologic and hydrologic information was accumulated in the course of the earthquake studies, and many new or improved investigative techniques were developed. Chief among these, perhaps, were the recognition that lakes can be used as giant tiltmeters, the refinement of methods for measuring land-level changes by observing displacements of barnacles and other sessile organisms, and the relating of hydrology to seismology by worldwide study of hydroseisms in surface-water bodies and in wells. The geologic and hydrologic lessons learned from studies of the Alaska earthquake also lead directly to better definition of the research needed to further our understanding of earthquakes and of how to avoid or lessen the effects of future ones. Research is needed on the origins and mechanisms of earthquakes, on crustal structure, and on the generation of tsunamis and local waves. Better earthquake-hazard maps, based on improved knowledge of regional geology, fault behavior, and earthquake mechanisms, are needed for the entire country. Their preparation will require the close collaboration of engineers, seismologists, and geologists. Geologic maps of all inhabited places in earthquake-prone parts of the country are also needed by city planners and others, because the direct relationship between local geology and potential earthquake damage is now well understood. Improved and enlarged nets of earthquake-sensing instruments, sited in relation to known geology, are needed, as are many more geodetic and hydrographic measurements. Every large earthquake, wherever located, should be regarded as a full-scale laboratory experiment whose study can give scientific and engineering information unobtainable from any other source. Plans must be made before the event to insure staffing, funding, and coordination of effort for the scientific and engineering study of future earthquakes. Advice of earth scientists and engineers should be used in the decision-making processes involved in reconstruction after any future disastrous earthquake, as was done after the Alaska earthquake. The volume closes with a selected bibliography and a comprehensive index to the entire series of U.S. Geological Survey Professional Papers 541-546. This is the last in a series of six reports that the U.S. Geological Survey published on the results of a comprehensive geologic study that began, as a reconnaissance survey, within 24 hours after the March 27, 1964, Magnitude 9.2 Great Alaska Earthquake and extended, as detailed investigations, through several field seasons. The 1964 Great Alaska earthquake was the largest earthquake in the U.S. since 1700. Professional Paper 546, in 1 part, describes Lessons and Conclusions.

Alaska↗

Shear-wave velocity site characterization in Oklahoma from joint inversion of multi-method surface seismic measurements: Implications for central U.S. Ground Motion Prediction

We analyze multimethod shear (SH)‐wave velocity ( ⁠ V S "> V S ) site characterization data acquired at three permanent and 25 temporary seismograph stations in Oklahoma that recorded M 4+ earthquakes within a 50 km hypocentral distance of at least one of the 2016 M 5.1 Fairview, M 5.8 Pawnee, or M 5.0 Cushing earthquakes to better constrain earthquake ground‐motion modeling in the region. We acquired active‐source seismic data for time‐averaged V S "> V S to 30 m depth ( V S "> V S V S 30 "> 30 ) at 28 seismograph stations near the Fairview, Pawnee, and Cushing epicentral areas. The SH‐wave refraction travel times coupled with Rayleigh‐ and Love‐wave phase velocity dispersion were extracted and modeled in a nonlinear least‐squares (L2) joint inversion to obtain a best‐fit 1D V S "> V S versus depth profile for each site. At a subset of sites where the preferred L2 inverse model did not optimally fit each of the Love, Rayleigh, and SH travel‐time datasets, we explore application of simulated annealing in a joint inversion to find a more global solution. V S "> V S V S 30 "> 30 values range from 262 to 807    m / s "> 807 m / s for the preferred measured (in situ) V S "> V S profiles, or National Earthquake Hazards Reduction Program (NEHRP) site class D to B, and are broadly comparable with estimates from previous data reports in the region. Site amplification estimates were calculated next from 1D SH transfer functions of the preferred V S "> V S VS profiles and then compared against observed horizontal‐to‐vertical spectral ratios (HVSRs) from nearby seismograph stations. We generally see good agreement between the predicted in situ model and the observed HVSR resonant frequencies, with nominal amplifications between 2 and 10 within the 2–15 Hz frequency band. Next, using 40 known in situ V S "> V S V S 30 "> 30 measurements in the region, we demonstrate that the in situ V S "> V S V S 30 "> 30 values improve the fit for selected suites of ground‐motion models (GMMs) for M 4+ earthquakes within a 50 km hypocentral distance when compared with proxy methods, arguing for future development of GMMs implementing in situ V S "> V S profiles.

Oklahoma↗

The 26 January 2001 M 7.6 Bhuj, India, earthquake: Observed and predicted ground motions

Although local and regional instrumental recordings of the devastating 26, January 2001, Bhuj earthquake are sparse, the distribution of macroseismic effects can provide important constraints on the mainshock ground motions. We compiled available news accounts describing damage and other effects and interpreted them to obtain modified Mercalli intensities (MMIs) at >200 locations throughout the Indian subcontinent. These values are then used to map the intensity distribution throughout the subcontinent using a simple mathematical interpolation method. Although preliminary, the maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated toward the western edge of the inferred fault, consistent with western directivity. Significant sediment-induced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region, intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium, such as mudflats and salt pans. In addition, we use fault-rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0-1000 km. We then convert the predicted hard-rock ground-motion parameters to MMI by using a relationship (derived from Internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1-3 units than those estimated from news accounts, although they do predict near-field ground motions of approximately 80%g and potentially damaging ground motions on hard-rock sites to distances of approximately 300 km. For the most part, this discrepancy is consistent with the expected effect of sediment response, but it could also reflect other factors, such as unusually high building vulnerability in the Bhuj region and a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. The discrepancy may also be partly attributable to the inadequacy of the empirical relationship between MMI and peak ground acceleration (PGA), when applied to India. The MMI-PGA relationship was developed using data from California earthquakes, which might have a systematically different stress drop and therefore, a different frequency content than intraplate events. When a relationship between response spectra and MMI is used, we obtain larger predicted MMI values, in better agreement with the observations.

Bulletin of the Seismological Society of America↗

Estimation of ground motion for Bhuj (26 January 2001; Mw 7.6) and for future earthquakes in India

Only five moderate and large earthquakes ( M w ≥5.7) in India—three in the Indian shield region and two in the Himalayan arc region—have given rise to multiple strong ground-motion recordings. Near-source data are available for only two of these events. The Bhuj earthquake ( M w 7.6), which occurred in the shield region, gave rise to useful recordings at distances exceeding 550 km. Because of the scarcity of the data, we use the stochastic method to estimate ground motions. We assume that (1) S waves dominate at R < 100 km and Lg waves at R ≥ 100 km, (2) Q = 508 f 0.48 is valid for the Indian shield as well as the Himalayan arc region, (3) the effective duration is given by fc -1 + 0.05R, where fc is the corner frequency, and R is the hypocentral distance in kilometer, and (4) the acceleration spectra are sharply cut off beyond 35 Hz. We use two finite-source stochastic models. One is an approximate model that reduces to the ω 2 -source model at distances greater that about twice the source dimension. This model has the advantage that the ground motion is controlled by the familiar stress parameter, Δ σ . In the other finite-source model, which is more reliable for near-source ground-motion estimation, the high-frequency radiation is controlled by the strength factor, sfact , a quantity that is physically related to the maximum slip rate on the fault. We estimate Δ σ needed to fit the observed Amax and Vmax data of each earthquake (which are mostly in the far field). The corresponding sfact is obtained by requiring that the predicted curves from the two models match each other in the far field up to a distance of about 500 km. The results show: (1) The Δ σ that explains Amax data for shield events may be a function of depth, increasing from ∼50 bars at 10 km to ∼400 bars at 36 km. The corresponding sfact values range from 1.0-2.0. The Δ σ values for the two Himalayan arc events are 75 and 150 bars ( sfact = 1.0 and 1.4). (2) The Δ σ required to explain Vmax data is, roughly, half the corresponding value for Amax, while the same sfact explains both sets of data. (3) The available far-field Amax and Vmax data for the Bhuj mainshock are well explained by Δ σ = 200 and 100 bars, respectively, or, equivalently, by sfact = 1.4. The predicted Amax and Vmax in the epicentral region of this earthquake are 0.80 to 0.95 g and 40 to 55 cm/sec, respectively.

Bulletin of the Seismological Society of America↗

Microearthquake spectra from the Anza, California, seismic network: Site response and source scaling

We analyzed spectra of local microearthquakes recorded by the Anza, California, seismic network to isolate the effects of site response and to investigate the scaling of source parameters for small earthquakes. Spectra of microearthquakes (M < 2; Mo < 1019 dyne-cm) at Anza have shapes characteristic of the receiver sites and are generally independent of the source region. Thus, the site response is a major conditioner of the observed spectral shape. To remove the effects of site response from the spectra of a M ∼ 3 event and isolate its source spectrum, we divided by the spectra of an adjacent aftershock used as an empirical Green's function event. The spectral ratios indicate that the apparent corner frequencies of small earthquakes (Mo < 1019dyne-cm) observed at even the high-fmax stations on hard rock are much lower than the source corner frequencies. The spectral ratios are consistent with stress drop remaining constant with decreasing seismic moment, for events with moments as small as 1018 dyne-cm. The spectral ratios display remarkable agreement between sites which showed vast differences in their original spectra, indicating that the spectral division effectively removed the site response. The source spectrum of the M ∼ 3 event has a high-frequency spectral fall-off of about ω−2. An apparent dependence of high-frequency fall-off with seismic moment in the original spectra can also be explained by the effects of site response. The difference between the P- and S-wave corner frequencies and high-frequency roll-offs in the observed spectra for these events is the result of the site response and is not a source property. The shapes of the spectra of microearthquakes at Anza can largely be explained by attenuation at shallow depth with a frequency-independent Q. For some sites, near-surface resonances are also apparent in the spectra of microearthquakes. It is indicated by t* values determined for each site that Qp ∼ Qs for the shallow low-Q layer. Further evidence of low near-surface Q is observed in the anomalous spectra of an unusually shallow earthquake (source depth ∼ 1 km) in the network. The spectra from this shallow event are depleted in high-frequency energy at most stations, relative to those of deeper events. This observation can be explained by a low-Q surficial zone. For stations of the network situated on alluvium, this low-Q layer has a maximum thickness of about 3 km and maximum P- and S-wave Q values of 30 to 50.

California↗