Program For Array Seismic Studies Of The Continental Lithosphere Research Articles

Data Products at the IRIS DMC: Stepping Stones for Research and Other Applications

In 2010, the Incorporated Research Institutions for Seismology (IRIS) Data Management Center (DMC) expanded the preexisting effort to generate, archive, and distribute data products derived from the extensive data archives of the center, which include the IRIS Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL), Global Seismographic Network (GSN), EarthScope, and many other data sources. This expansion included dedicating two new full‐time staff to data products and to the development of a new system to manage and distribute them. With community guidance provided by a data products working group (DPWG), DMC has developed a number of freely available unique and increasingly popular data products for the seismological community. In addition, DMC is also archiving and distributing other data products that are produced by community members. A complete list of data products managed by DMC is maintained here: http://www.iris.edu/dms/products. Many of these are available from the Searchable Product Depository (SPUD, http://www.iris.edu/spud), the DMC data product management system. Traditionally, data products result from specialized processing applied during research by the principal investigators and their teams. Products span the range of simple plots or animations derived from waveform data to highly refined models. Some waveform‐processing techniques (e.g., receiver functions) have gained wide community use and can be derived by widely accepted algorithms. Such first‐order data products make ideal candidates for routine generation to provide a standard for the user community. Data products are typically derived from raw data to serve as insightful summaries of the raw time series data, or as quality control metrics, or as foundations or stepping stones to more specialized products and research. The overall goal of the DMC data product effort is to produce and manage a variety of community‐approved and utilized data products that complement the raw data managed at the DMC and ultimately serve …

An Improved Method to Extract Very-Broadband Empirical Green's Functions from Ambient Seismic Noise

The extraction of Green’s functions by cross correlation of continuous seismic records at station pairs is best achieved in a diffuse wave field, where energies radiated by random sources have equal power. To partially satisfy the diffuse wave‐field condition, original seismic records must be normalized as they are highly nonstationary and may have amplitude variations of orders of magnitude within the cross‐correlation time window. We adapt a frequency–time normalization method to obtain seismograms with an even spectrum at all times within the data‐processing unit. Compared with the commonly used one‐bit normalization, the new method improves the signal‐to‐noise ratio of empirical Green’s functions by a factor of ∼2 and increases the effective data‐recording duration by a factor of 4, assuming random local and instrument noise. It yields useful Rayleigh waves at periods up to 300 s for Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL)–type deployments and 600 s for permanent stations with very‐broadband sensors. Thus, the new method makes it possible to extend surface‐wave signals to frequencies beyond those in traditional earthquake‐based surface‐wave tomography at both the high‐ and low‐frequency ends.

Rayleigh wave phase velocity maps of Tibet and the surrounding regions from ambient seismic noise tomography

Ambient noise tomography is applied to the significant data resources now available across Tibet and surrounding regions to produce Rayleigh wave phase speed maps at periods between 6 and 50 s. Data resources include the permanent Federation of Digital Seismographic Networks, five temporary U.S. Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) experiments in and around Tibet, and Chinese provincial networks surrounding Tibet from 2003 to 2009, totaling ∼600 stations and ∼150,000 interstation paths. With such a heterogeneous data set, data quality control is of utmost importance. We apply conservative data quality control criteria to accept between ∼5000 and ∼45,000 measurements as a function of period, which produce a lateral resolution between 100 and 200 km across most of the Tibetan Plateau and adjacent regions to the east. Misfits to the accepted measurements among PASSCAL stations and among Chinese stations are similar, with a standard deviation of ∼1.7 s, which indicates that the final dispersion measurements from Chinese and PASSCAL stations are of similar quality. Phase velocities across the Tibetan Plateau are lower, on average, than those in the surrounding nonbasin regions. Phase velocities in northern Tibet are lower than those in southern Tibet, perhaps implying different spatial and temporal variations in the way the high elevations of the plateau are created and maintained. At short periods (<20 s), very low phase velocities are imaged in the major basins, including the Tarim, Qaidam, Junggar, and Sichuan basins, and in the Ordos Block. At intermediate and long periods (>20 s), very high velocities are imaged in the Tarim Basin, the Ordos Block, and the Sichuan Basin. These phase velocity dispersion maps provide information needed to construct a 3‐D shear velocity model of the crust across the Tibetan Plateau and surrounding regions.

Field Geophysics Class at Volcán Tungurahua, Ecuador

Ecuador's erupting Volcán Tungurahua was the recent site of a 3‐week graduate‐level geophysical course on volcanoes, hosted by Ecuador's Instituto Geofisico of the Escuela Politecnica Nacional (IG‐EPN) and the Department of Earth Science at the New Mexico Institute of Mining and Technology (NMT). Sixteen students from 12 universities and four countries participated in the intensive course, which entailed broadband seismometer and infrasound sensor deployment followed by subsequent data processing, analysis, interpretation, and result synthesis. Hardware for the course was provided by the Incorporated Research Institutes for Seismology (IRIS) through the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) as well as the IG‐EPN and NMT geophysics programs.

Teleseismic Receiver Function and Surface-Wave Study of Velocity Structure beneath the Yanqing-Huailai Basin Northwest of Beijing

Abstract Shear-wave velocities beneath the Yanqing-Huailai Basin, 90–140 km northwest of Beijing, are estimated from the joint inversion of surface-wave phase velocities and teleseismic receiver functions. The data set is from a temporary broadband seismic network supported by the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) in the basin and includes 34 teleseismic events from 2003 to 2005. Receiver functions from the teleseismic events are similar for the stations around the Yanqing-Huailai Basin and exhibit little variation with azimuth. The velocity models constrained by receiver functions and surface-wave dispersion curves are also similar. The resulting models reflect the low-velocity basin sediments to 2 km followed by a positive velocity gradient to 15 km with shear-wave velocity increasing from 2.0 to 3.55 km/sec. Evidence of a midcrust low-velocity layer starts at 15 km with a shear velocity decrease to 3.3 km/sec that extends to approximately 25 km. The total crustal thickness is 38–42 km with a smooth Moho transition to an upper-mantle shear velocity of 4.3 km/sec. The low-velocity zone is consistent with recent extension, geothermal activity, and earthquake locations above this depth. The average shear velocity model for the basin has similarities to other regional and global models but provides more detailed structure in the uppermost and lower portions of the crust. The new model includes the effect of the sediments in the basin, the low-velocity layer, and the gradual Moho transition. Predicted P - and S -travel times are 1–3.5 sec slower than the previous models at regional distances.

Seismic Recording on Drifting Icebergs: Catching Seismic Waves, Tsunamis and Storms from Sumatra and Elsewhere

Under the Incorporated Research Institutions for Seismology (IRIS) Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) project SOUTHBERG, we have operated for the past three years a number of seismic stations on giant Antarctic icebergs parked or drifting slowly along the margins of the Ross Sea. The purpose of this deployment was to investigate in situ the characteristics and origin of high-frequency tremors emanating from inside the icebergs, tremors that previously had been detected in the far field as hydroacoustic T phases recorded in French Polynesia (Talandier et al. 2002, 2006). Similar observations were made later in the Indian Ocean by Chapp et al. (2005). In this general context, the purpose of this paper is to report on the ancillary functioning of iceberg-sited seismometers as teleseismic observatories, and in particular their recording of the two Sumatra mega-earthquakes of 26 December 2004 and 28 March 2005. These data are available upon request from the authors. We are in the third and final year of the SOUTHBERG deployment. In November–December 2003, we installed four seismic stations on iceberg C16 (including an STS-2 broadband instrument at the central station C16A) and left one of them to winter over during 2004. Because all stations work on solar power, the station shut down on 28 May 2004 after the arrival of permanent winter darkness at that latitude, but it woke up and started recording uneventfully on 25 September 2004 around the austral spring equinox, thus providing a total of 141 days of seismic recording, which went beyond our expectations. Bolstered by this performance, we left four stations to winter over in 2005: one at C16A (hereafter referred to as C16 for short), one on iceberg B15A, one on iceberg B15K (a fragment of B15A that detached in 2004), and one on …

Variations in Short-period Geophone Responses in Temporary Seismic Arrays

Geophones with resonant frequencies of 1–4 Hz are commonly used as short-period sensors for site-response studies ( e.g. , Bonilla et al. , 1997; Hartzell et al. , 1997; Frankel et al. , 1999; Hartzell et al. , 2000) and crustal-reflection, refraction, and tomography experiments ( e.g. , Brocher et al. , 2001; Fuis et al. , 2003). In many cases these studies require only arrival-time picks, but site-response and attenuation measurements require accurate ground-velocity data from these sensors. Moreover, these geophones are often used to record seismic phases at or below their nominal resonant frequencies ( e.g. , Bonilla et al. , 1997; Hartzell et al. , 2000; Dutta et al. , 2001; Pratt et al. , 2003b; Pratt et al. , 2003a). For studies that require accurate ground-velocity measurements or that use a variety of sensor types, the instrument response must be removed from the data using either the manufacturer's specifications or results from geophone calibration tests. Geophone placement, particularly any tilt to the sensor, can alter the instrument response significantly. The geophone response ideally is measured while the sensors are in the ground, but such in situ calibrations generally are not carried out because of time constraints. More commonly, it simply is assumed the geophones have the response characteristics specified by the manufacturer. In a recent site-response study in the Seattle area, named the Seattle Seismic Hazard Investigation of Puget Sound (Seattle SHIPS 2002; Pratt et al. , 2003b), we used 72 three-component geophones from the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) instrument pool to determine the amplitudes of both horizontal and vertical ground shaking at sites over the Seattle sedimentary basin (Pratt et al. , 2003b). During the experiment, we used the signal-coil calibration method (Rogers …

IRIS Seismology Program marks 20 years of discovery

The year 2004 marked the twentieth anniversary of the Incorporated Research Institutions for Seismology (IRIS) Consortium and its four core programs: the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL), the Global Seismic Network (GSN), Education and Outreach (E&O),and the Data Management System (DMS).IRIS activities (www.iris.edu) are supported principally through the U.S. National Science Foundation's (NSF) Earth Sciences Instrumentation and Facilities Program, and through the community efforts of IRIS's more than 150 member institutions in the United States and abroad. The consortium provides fundamental capabilities for the exploration of the Earth's interior through the collection and distribution of seismographic data. PASSCAL is the IRIS facility program dedicated to seismographic data collection associated with temporary deployments using highly portable instrumentation.

Lithospheric assembly and modification of the SE Canadian Shield: Abitibi‐Grenville teleseismic experiment

This paper presents the results of a joint Lithoprobe‐Incorporated Research Institutions for Seismology (IRIS)/Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) teleseismic experiment that investigates portions of the Grenville and Superior Provinces of the Canadian Shield along the Québec‐Ontario border. Data from a 600‐km‐long, N‐S array of 28 broadband seismographs deployed between May and October 1996 have been supplemented with additional recordings from an earlier 1994 deployment and from stations of the Canadian National Seismograph Network and the Southern Ontario Seismic Network. Relative delay times of P and S waves from 123 and 40 teleseismic events, respectively, have been inverted for velocity perturbations in the upper mantle and reveal a low‐velocity, NW‐SE striking corridor that crosses the southern portion of the line at latitude 46°N and lies between 50 and 300 km depth. Multievent S K S‐splitting results yield an average delay time of 0.57±0.22 s and a direction of fast polarization of N93°E±18°, which is consistent with an earlier interpretation as being due to fossil strain fields related to the last major regional tectonic event. Subtle variations in splitting parameters over the low‐velocity corridor may suggest an associated disruption in mantle fabric. Profiling of radial receiver functions reveals large and abrupt variations in Moho topography, specifically, a gradual thickening in crust from 40 to 45 km between latitudes 45°N and 46°N, which is followed by an abrupt thinning to 35 km at 46.6°N, some 65 km southeast of the Grenville Front. This structure is interpreted as a subduction suture extending the full length of the Front and punctuating a major pre‐Grenvillian (Archean‐Proterozoic) episode of lithospheric assembly in the southeast Canadian Shield. The low‐velocity mantle corridor, by contrast, is better explained as the extension of the Monteregian‐White Mountain‐New England seamount hotspot track below the craton and is here postulated to represent interaction of the Great Meteor plume with zones of weakness within the craton developed during earlier rifting episodes.

Seismic anisotropy and mantle flow from the Great Basin to the Great Plains, western United States

Shear wave splitting and P, SKS, and S travel time residuals are calculated for teleseismic arrivals recorded on the Colorado Plateau‐Great Basin Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) portable broadband seismic deployment and for permanent stations in the western United States. Little shear wave splitting is observed for broadband recordings in the northern Colorado Plateau, the Rocky Mountains, or the central Great Basin. The transition between the Colorado Plateau and the Great Basin is marked by moderate shear wave splitting (1.0 s) and unusually late teleseismic phase arrivals. This suggests material with a higher content of mantle melt or volatiles than regions to either side. Splitting in the transition between the Colorado Plateau and Great Basin is part of a pattern of fast polarizations that align in a semicircle, surrounding a central Great Basin region of null (no splitting) measurements. Away from the California plate boundary, splitting to the north and south of our study region aligns roughly parallel to the absolute plate motion of the North American plate. No simple spatial relation of splitting with geological and geophysical features such as mountain ranges, velocity anomalies, gravity, magnetics, or heat flow is evident in most of the western United States. However, splitting in the Great Basin is compatible with asthenospheric flow. The smallest shear wave splitting delay times coincide with the Eureka Low in heat flow, also having low S velocity at 300 km depth. We suggest that the circumferential pattern of fast polarization directions ringing a central region of nulls in the Great Basin is caused by mantle flow, by the interaction of mantle up welling and the absolute motion of the North American plate.

Mantle discontinuity structure beneath the Colorado Rocky Mountains and High Plains

Analysis of mantle discontinuity structure using converted P to S (PaS) phases beneath Colorado from the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) Rocky Mountain Front (RMF) experiment reveals significant topography at the 410 and 660 km depth discontinuities and corresponding transition zone thickness variations. A stack of all radial receiver functions resolves the 410 and 660 km discontinuities at average depths of 419 and 677 km, respectively. Imaging of lateral variations in mantle discontinuity structure is accomplished by geographically binning the Pds conversion points and then stacking the receiver functions in each bin to form spatial images, analogous to common depth point stacking. Corrections for lateral velocity heterogeneity are calculated using the local S wave tomographic model of Lee and Grand [1996] and a constant ∂lnVs/∂lnVP scaling of 1.3. This scaling value is determined from the relative scaling between teleseismic P and S wave travel time residuals measured from the Rocky Mountain Front deployment. Mantle discontinuity images using 150 km square bins show 20 km of 410 km discontinuity topography, 30 km of 660 km discontinuity topography, and up to 40 km of transition zone thickness variation. Features of the discontinuity structure include a 20 km depression of the 660 km discontinuity beneath western Colorado and a gradual 10 km dip of the 410 km discontinuity beneath the High Plains. The thickening of the transition zone beneath southwest Colorado is consistent with the presence of the subducted Farallon slab in this region as imaged by Van der Lee and Nolet [1997]. In general, our results show that the transition zone discontinuity structure is more complex than that predicted by the simple model of olivine phase boundaries modulated by vertically coherent thermal anomalies.

Estimation of crustal stochastic parameters from seismic exploration data

Stochastic models for the crystalline crust produce synthetic seismograms that compare well with recorded data for a variety of crustal ages and tectonic environments. In this paper, we explore the parameter space describing such stochastic models as a basis for formulating the inverse problem; that is, we wish to estimate the parameters which define a stochastic model from the recorded backscattered wave field. We base the estimation on approximate relations between the primary reflectivity series, which is the ideal wave field response of a medium, and various seismic gathers. A two‐dimensional lateral correlation method is used to investigate the sensitivity of synthetic wave fields to horizontal characteristic length. A derived empirical relationship relates the scale length to the half width of the correlation coefficient. A horizontal wavenumber‐time domain spectral analysis successfully estimates the horizontal characteristic length (ax) and the fractal dimension (D) of the stochastic medium from which the wave field was backscattered. We employ a Monte Carlo based optimization scheme to fit the observed spectral values to a one‐dimensional von Kármán function and simultaneously estimate values for ax and the Hurst number ν. The method yields consistent results when tested on synthetic data. Application of this technique to data sorted into different domains shows that zero‐offset sections yield the best estimate. We also present results of our estimates using this method on a Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) dynamite source shot gather and a vibroseis‐source Consortium for Continental Reflection Profiling (COCORP) data set recorded in the Basin and Range Province.

Crustal thickness variations across the Colorado Rocky Mountains from teleseismic receiver functions

Variations in crustal thickness from the Great Plains of Kansas, across the Colorado Rocky Mountains, and into the eastern Colorado Plateau are determined by receiver function analysis of broadband teleseismic P waveforms recorded during the 1992 Rocky Mountain Front Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) experiment. The receiver functions are calculated using a time domain deconvolution approach and are interpreted in terms of a single crustal layer, with thickness determined by a grid‐search comparison of observed receiver functions with synthetics. The average crustal thicknesses determined by these methods are Kansas Great Plains, 43.8±0.4 km; Colorado Great Plains, 49.9±1.2 km; Colorado Rocky Mountains, 50.1±1.3 km; and northeast Colorado Plateau, 43.1±0.9 at latitudes of 38°–40°N. The main variations in crustal thickness that we observe are between the Kansas Great Plains and the Colorado Great Plains and between the Rocky Mountains and the Colorado Plateau. There is not a significant crustal thickness difference between the Colorado Great Plains and the Colorado Rocky Mountains. Together with gravity data and mass balance calculations, these results are incompatible with the hypothesis that the compensation of the Rocky Mountains relative to the Great Plains is accommodated purely by an Airy‐type crustal root or any other mechanism that restricts compensation solely to the crust and requires significant support for the excess topography of the Rocky Mountains to come from the mantle. Models with a rigid elastic plate may match receiver function estimates of crustal thickness but underpredict the amplitude of the gravity low over the Rockies. Our favored model includes lateral variations in crustal velocities obtained from refraction studies and crustal thickness variations constrained by the receiver functions. These models indicate that there is a profound transition in mantle density structure near the eastern range front.

Active thrust front of the Greater Caucasus: The April 29, 1991, Racha earthquake sequence and its tectonic implications

Although fault‐bounded thrust sheets are common in the geological record, seismic evidence for their motion is sparse. The April 29, 1991, Racha earthquake (Ms = 7.0), the largest instrumentally recorded earthquake in the Greater Caucasus, is one of the largest recent earthquakes in continental thrust belts and provides evidence on mechanisms of thrust sheet motion. Using data from a deployment of Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) digital seismographs and various other instruments, we locate 1952 aftershocks occurring between May 7 and June 30, 1991. The aftershocks form a zone ∼70 km long and 10–25 km wide striking E‐W, following the Racha ridge at the southern boundary of the Greater Caucasus thrust system. Teleseismic body waves are inverted for source parameters of the mainshock and the two largest aftershocks. The solutions show thrust faulting with centroid depths of 3–10 km, comparable to depths of locally recorded aftershocks (∼2–12 km). The shallow‐dipping nodal plane, the aftershock distribution, and surface geology demonstrate that the main event was caused by faulting on a thrust system dipping NNE at 20°–31° bounding the southern slope of the Greater Caucasus. This fault system thrusts the Greater Caucasus structures south over the Dzhirula basement massif. The inferred fault geometry suggests that the active fault is either a detachment between sediments and Dzhirula basement or cuts through the basement at shallow depths. The 1500‐m‐high Racha ridge overlies the aftershock zone and is a likely consequence of repeated similar earthquakes. Hence the 1991 earthquake sequence shows that the western Greater Caucasus is accommodating plate convergence at a rate possibly comparable to the eastern Greater Caucasus (a few millimeters per year). Along‐strike geological discontinuities above and below the thrust surface correspond to the eastern end of the mainshock rupture area. No strong evidence for transfer structures could be found along strike, suggesting that differences in collisional style between the western and eastern Greater Caucasus may reflect differences in mechanical properties rather than differences in convergence rate. A June 15, 1991, event and its aftershocks, southeast of the primary aftershock zone along strike, show fault planes and slip vectors rotated ∼41° clockwise from the mainshock. This rotation is consistent with an along‐strike change in direction of the thrust front, near 44°E longitude, and demonstrates strong local structural or topographic control on slip direction. The rotation requires along‐strike shortening within the Greater Caucasus thrust system at a rate comparable to the rate of thrusting.

Mantle anisotropy beneath the Tibetan Plateau: evidence from long-period surface waves

We investigate long-period surface waves ( T > 100 s) propagating across the Tibetan Plateau, based on the seismic records from a PASSCAL portable array and nearby Chinese Digital Seismic Network (CDSN) stations. Significant quasi-Love waveform anomalies, associated with fundamental Rayleigh-Love coupling, are consistently observed, suggesting that strong lateral gradients in azimuthal anisotropy exist beneath the plateau. More intense surface-wave scattering is observed at intermediate periods ( T ≈ 50 s), which suggests that the cause of the long-period quasi-Love waves lies in the upper mantle, not the thickened crust of the plateau. A detailed analysis of the data indicates that the long-period waveform anomalies are generated beneath the central Tibetan Plateau, where the structural trend implies deep deformations induced by continental collision. The absence of quasi-Love anomalies at the westernmost station of the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) array suggests that the east-west extent of this mantle deformation is limited. Some, but not all, of the quasi-Love observations are consistent with SKS splitting observations. Both sets of observations predict a strong gradient in anisotropic properties in central Tibet. However, the quasi-Love waveforms are absent from records collected in the northern plateau for northerly propagation paths, which is not consistent with a gradient at the northern edge of the plateau, as suggested by SKS studies. This discrepancy indicates that the simple models used to interpret both body- and surface-wave data may be inadequate. Either significant P-wave anisotropy exists under the plateau, or the S-wave anisotropy does not possess a uniformly horizontal symmetry axis, or strong east-west gradients in anisotropy bias our data-synthetic comparison.

The FDSN Archive at the IRIS Data Management Center

In August of 1990, the IRIS Data Management Center (IRIS DMC) was designated as the first FDSN Data Center. Since that time the IRIS DMC has also come to be known as the FDSN Archive. As the FDSN Archive, all member networks, through their own local Network Data Centers, routinely send data to the IRIS DMC. At the present time we have received data from seven member networks. These include: CNSN (Canada); CDSN (China); GEOSCOPE (France); Graefenberg (Germany); IRIS GSN (U.S.A.); MEDNET (Italy); POSEIDON (Japan). In addition to acting as the FDSN Archive, the IRIS DMC is responsible for archiving and distributing data from several sources. The Global Seismographie Network (GSN) and the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) and the Joint Seismic Program (JSP) are three IRIS programs that presently send data to the IRIS DMS. The DMC also archives data from some regional networks in the United States and in the Former Soviet Union (FSU). The DMC is one component of a larger Data Management System (DMS) consisting of several components including two Data Collection Centers (DCCs), one Data Management Center (DMQ, a DMC Host at the University of Washington, a Wave form Quality Center (WQC), and the Moscow Data Analysis Center (MDC) in Russia.

Three‐dimensional attenuation tomography at Loma Prieta: Inversion of t* for Q

Three‐dimensional Q−1 variations in the aftershock region of Loma Prieta are derived by tomographic inversion. The data set consists of over 4000 aftershock recordings at 22 PASSCAL (Program for Array Seismic Studies of the Continental Lithosphere) stations deployed after the Loma Prieta mainshock of 1989. Estimates of attenuation are determined from nonlinear least squares best fits to the Fourier amplitude spectrum of P and S wave arrivals. The linear attenuation inversion is accomplished by using three‐dimensional velocity variations derived previously in nonlinear velocity inversions. Low Q is observed near the surface and Q generally increases with depth. The southwest side of the San Andreas fault exhibits lower Q than does the northeast side and this feature apparently extends to approximately 7 km depth. The fault zone, as determined by the dipping plane of aftershock activity, is characterized by slightly higher Qp and lower Qs, compared to regions immediately adjacent to the fault. These correlate with high‐velocity anomalies associated with seismicity at depth. The results are in agreement with earlier observations regarding the association of high‐velocity anomalies, seismicity, and fault zone asperities.

Azimuthal shear wave velocity anisotropy in the Basin and Range Province using moho Ps converted phases

An eight‐station large aperture seismometer array was deployed in the Basin and Range province of west‐central Nevada for the 1988–1989 Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) Passive‐source seismic experiment. During the 10 months of data collection a total of 100 teleseismic events were recorded (Δ > 30°). Source‐equalized P wave receiver functions show waveform variations in timing, amplitude, and polarity in the Moho Ps phase that is diagnostic of shear wave velocity anisotropy. Information regarding the crustal anisotropic component can be recovered through the Ps phase since it is generated at the Moho and is sensitive to shear wave splitting in the crust only. We utilize single‐event, stacked, and array‐beam‐formed receiver functions to minimize the influence of structural variations and scattering on the observed converted phases. Analysis of shear wave splitting of the Ps conversion from the crust‐mantle boundary indicates a fast azimuth of anisotropy oriented approximately NW‐SE with observed time delays of 0.20 s. This strongly suggests coherent crustal anisotropy with a fast direction perpendicular to the strike of fault block mountain ranges. This fabric is most likely due to the preferred crystallographic orientation of seismically anisotropic minerals in the middle to lower crust rather than the distribution of Basin and Range upper crustal fracture systems.

Three-dimensional Vp/Vs variations along the Loma Prieta rupture zone

An inversion of S - P seismic wave arrival time data primarily from the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) array deployed in the aftershock zone of the 1989 Loma Prieta, California, earthquake has yielded a 3-D model for Vp/Vs variations along the Loma Prieta rupture zone. Significant effort was devoted to identifying S -to- P converted waves, which were prevalent at a few of the stations used. The inversion method uses 3-D ray tracing in an existing laterally heterogeneous P -wave velocity model to determine the initial paths for S waves, assuming a constant Vp/Vs . The model represents spatial variations in Vp/Vs on a 3-D grid with node spacings of 3 to 15 km horizontally and 3 to 5 km vertically. The reduction in data variance, computed exactly by applying the Vp/Vs perturbations to the S -velocity structure, exceeds 65% for the single-step inversion. We interpret the Vp/Vs model in conjunction with other available geological and geophysical information. Our major findings are 1) high Vp/Vs along the fault zone in the upper crust, indicative of upthrust gabbroic rocks, 2) predominantly low Vp/Vs in the lower crust, indicative of an underlying quartz-rich basement, and 3) two areas of low Vp/Vs in the near-surface associated with known quartz-rich rock sequences, one a sandstone and the other granite and/or granodiorite. We also tentatively identify a boundary between areas of low and slightly high Vp/Vs at depth in the immediate vicinity of the Loma Prieta mainshock hypocenter. Our results indicate the diagnostic power of a 3-D Vp/Vs model when it is combined with other geological and geophysical information.

The crustal structure of the Northwestern Basin and Range Province, Nevada, from wide‐angle seismic data

I present a model of lithospheric P‐wave velocity structure in the northwestern Basin and Range province, Nevada, based on wide‐angle seismic data collected in 1986 under the auspices of the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL). The lithosphere in the Basin and Range consists of (1) 2‐ to 3‐km‐thick low‐velocity (2.5–4.0 km/s) sedimentary basins and high‐velocity (>4.5 km/s) ranges overlying a basement of velocity 5.5–6.1 km/s, (2) a mid‐crust (6.1–6.2 km/s) which consists of rocks of probable granitic to granodioritic composition in greenschist to amphibolite facies and which includes two local low‐velocity zones (5.5–5.7 km/s), (3) a laterally varying lower crust which consists of rocks of probable gabbroic composition in upper amphibolite to granulite fades and which comprises layers of 6.6 and 7.0–7.2 km/s average velocity, (4) a complex crust‐mantle transition best modeled as a velocity gradient zone 1–2 km thick, and (5) an upper mantle with velocities of 7.7–7.9 km/s which includes a local low‐velocity zone (7.5 km/s). Crustal thickness in the region varies laterally from 28 to 37 km, significantly thicker than several previous estimates. The upper mantle velocities and velocity gradients differ slightly on the two perpendicular profiles, suggesting possible azimuthal anisotropy. A wedge‐shaped, eastward‐thickening lower crust of high seismic velocity (7.0–7.2 km/s) on the eastern pan of the survey may represent the westward edge of “rift‐stage” Précambrien continental crust or, alternatively, mafic material underplated and intruded into the crust during Cenozoic extension. Comparison of the velocity model with coincident seismic reflection data shows that the reflection and refraction Mohos are coincident across the entire profile, thus indicating that deep reflections (4–10 sec two‐way travel time) seen on COCORP and PASSCAL‐piggyback reflection profiles originate entirely within the crust. The top of the reflective zone, however, lies consistently above the mid‐crustal increase in average velocity, indicating that the reflectivity is not restricted to the mafic part of the crust; on PASSCAL piggyback profiles, in fact, the middle (silicic) crust appears more reflective than the lower (mafic) crust This can be explained as the result of either mafic, sill‐like intrusions or ductile strain fabrics. Both of these processes are likely to have accompanied Cenozoic extension in the region.