Crustal Structure and Relocated Earthquakes in the Puget Lowland, Washington, from High-Resolution Seismic Tomography T
Total Page:16
File Type:pdf, Size:1020Kb
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B12, 2381, doi:10.1029/2001JB000710, 2002 Crustal structure and relocated earthquakes in the Puget Lowland, Washington, from high-resolution seismic tomography T. M. Van Wagoner,1 R. S. Crosson, K. C. Creager, G. Medema, and L. Preston Department of Earth and Space Sciences, University of Washington, Seattle, Washington, USA N. P. Symons Sandia National Laboratories, Albuquerque, New Mexico, USA T. M. Brocher U.S. Geological Survey, Menlo Park, California, USA Received 13 June 2001; revised 19 July 2002; accepted 21 August 2002; published 31 December 2002. [1] The availability of regional earthquake data from the Pacific Northwest Seismograph Network (PNSN), together with active source data from the Seismic Hazards Investigation in Puget Sound (SHIPS) seismic experiments, has allowed us to construct a new high- resolution 3-D, P wave velocity model of the crust to a depth of about 30 km in the central Puget Lowland. In our method, earthquake hypocenters and velocity model are jointly coupled in a fully nonlinear tomographic inversion. Active source data constrain the upper 10–15 km of the model, and earthquakes constrain the deepest portion of the model. A number of sedimentary basins are imaged, including the previously unrecognized Muckleshoot basin, and the previously incompletely defined Possession and Sequim basins. Various features of the shallow crust are imaged in detail and their structural transitions to the mid and lower crust are revealed. These include the Tacoma basin and fault zone, the Seattle basin and fault zone, the Seattle and Port Ludlow velocity highs, the Port Townsend basin, the Kingston Arch, and the Crescent basement, which is arched beneath the Lowland from its surface exposure in the eastern Olympics. Strong lateral velocity gradients, consistent with the existence of previously inferred faults, are observed, bounding the southern Port Townsend basin, the western edge of the Seattle basin beneath Dabob Bay, and portions of the Port Ludlow velocity high and the Tacoma basin. Significant velocity gradients are not observed across the southern Whidbey Island fault, the Lofall fault, or along most of the inferred location of the Hood Canal fault. Using improved earthquake locations resulting from our inversion, we determined focal mechanisms for a number of the best recorded earthquakes in the data set, revealing a complex pattern of deformation dominated by general arc-parallel regional tectonic compression. Most earthquakes occur in the basement rocks inferred to be the lower Tertiary Crescent formation. The sedimentary basins and the eastern part of the Olympic subduction complex are largely devoid of earthquakes. Clear association of hypocenters and focal mechanisms with previously mapped or proposed faults is difficult; however, seismicity, structure, and focal mechanisms associated with the Seattle fault zone suggest a possible high-angle mode of deformation with the north side up. We suggest that this deformation may be driven by isostatic readjustment of the Seattle basin. INDEX TERMS: 7205 Seismology: Continental crust (1242); 7218 Seismology: Lithosphere and upper mantle; 7230 Seismology: Seismicity and seismotectonics; KEYWORDS: tomography, structure, Puget, earthquakes Citation: Van Wagoner, T. M., R. S. Crosson, K. C. Creager, G. Medema, L. Preston, N. P. Symons, and T. M. Brocher, Crustal structure and relocated earthquakes in the Puget Lowland, Washington, from high-resolution seismic tomography, J. Geophys. Res., 107(B12), 2381, doi:10.1029/2001JB000710, 2002. 1. Introduction 1Now at ChevronTexaco North America Upstream, New Orleans, [2] Although most damaging historical earthquakes in Louisiana, USA. the region have been within the Wadati-Benioff zone of Copyright 2002 by the American Geophysical Union. the Cascadia subduction zone (e.g., the 2001 magnitude 0148-0227/02/2001JB000710$09.00 6.8 Nisqually earthquake [Malone et al., 2001]), much ESE 22 - 1 ESE 22 - 2 VAN WAGONER ET AL.: CRUSTAL STRUCTURE AND EARTHQUAKES Figure 1. Schematic geologic map of northwestern Washington and southern British Columbia, showing the model region (heavy inner rectangle). Abbreviations for major cities: E—Everett; O— Olympia; S—Seattle; T—Tacoma; V—Victoria. Abbreviations for faults: CRBF—Coast Range boundary fault; DF—Doty fault; DDMF-Darrington—Devils Mountain fault; HCF—Hood Canal fault; HRF—Hurricane Ridge fault; LRF—Leech River fault; OF—Olympia fault; SCF—Straight Creek fault; SF—Seattle fault; SHZ—St. Helens zone; SJF—San Juan fault; SqF—Sequim fault; SWIF—Southern Whidbey Island fault; TF—Tacoma fault; WRF—White River fault. Abbreviations for oil test wells: DS—Dungeness Spit #1; MK—Mobil Kingston #1; P18—Pope-Talbot #18-1; SC—Silvana Community #12-1; SS—Socal Schroeder #1; SW—Socal Whidbey #1; WS—Washington State #1. Abbreviations for geologic, structural, and topographic features: BH—Blue Hills; BkH—Black Hills; DB—Dabob Bay; EB—Everett basin; HC—Hood Canal; KA—Kingston arch; LS—Lake Sammamish; LW—Lake Washington; M—Metchosin igneous complex; PTB—Port Townsend basin; PS—Puget Sound; SB— Seattle basin; SU—Seattle uplift; TB—Tacoma basin. Abbreviations for Cascade volcanoes: MSH—Mt. St. Helens; MR—Mt. Rainier; MB—Mt. Baker. Sources: Blakely et al. [2002]; Brocher et al. [2001]; Brocher and Ruebel [1998]; Gower et al. [1985]; Johnson et al. [1999]; Schruben et al. [1994]; Tabor and Cady [1978a]. recent geological and geophysical research in the Puget investigations [e.g., Brocher et al., 2001; Blakely et al., Lowland of western Washington State (Figure 1) has been 2002]. prompted by increasing concern about hazards from crustal [3] Lees and Crosson [1990] carried out the first 3-D earthquake sources. A postulated earthquake about 1100 tomography for the Puget Lowland region, but were limited years ago on the Seattle fault [Bucknam et al., 1992] and by relatively few earthquake observations and sparse station concern about other possible crustal faults [e.g., Johnson et spacing. Symons and Crosson [1997] developed new inver- al., 1996, 1999] in the Lowland region motivated the sion technology and applied this to an improved earthquake recent SHIPS active source experiments [Brocher et al., data set plus limited explosion data. In their study, shallow 1999, 2000], and analyses based on these and other resolution was again limited by lack of surface source to VAN WAGONER ET AL.: CRUSTAL STRUCTURE AND EARTHQUAKES ESE 22 - 3 Figure 2. Map of the Puget Lowland model area showing source/receiver locations used to construct the high resolution model. The map area is the heavy inner rectangle shown in Figure 1. Temporary Wet- SHIPS stations are shown as diamonds (à) and Wet-SHIPS air gun tracks are shown as thin line segments. Permanent PNSN seismograph stations are shown as black triangles. Station distributions for both the Dry-SHIPS and the 1991 experiment are shown as line segments and experimental sources are shown as white stars. Higher confidence 3-D relocated earthquake epicenters are shown as white circles (6), lower confidence events as white squares (5). Earthquake observations were recorded only by PNSN stations. surface receiver observations. Seismic reflection profiling our knowledge of the crustal structure in this region. In this has been interpreted for structural details in only the upper study, the simultaneous analysis of earthquake data and few kilometers at most of the crust [e.g., Johnson et al., active source observations has allowed us to develop a high 1994, 1996; Pratt et al., 1997]. More recently, a first arrival quality 3-D model to a depth of 30 km, considerably greater P wave tomography model for the upper 11 km of the crust than previously reported [e.g., Brocher et al., 2001], and to was constructed by Brocher et al. [2001] in the same region simultaneously relocate nearly 1000 regional earthquakes as our study using data from the initial SHIPS experiment. for tectonic interpretation. In these previous studies, the structural transition to mid- crustal depths has not been accurately imaged. Since the 1.1. Geologic Setting midcrust is likely to be the seat of important tectonic [5] The tectonics of the Pacific Northwest region are processes in this forearc region, potentially valuable struc- dominated by the oblique, northeastward subduction of tural information has been missing. the Juan de Fuca plate beneath western North America. [4] The availability of active source data from SHIPS as Subduction of oceanic lithosphere during the Tertiary has well as two earlier experiments [Miller et al., 1997; Parsons led to the formation of two major geologic provinces in the et al., 1999], together with high quality crustal earthquake region: the Coast Range province on the west, which data acquired for the last 20 years from the PNSN regional includes the Olympic subduction complex (OSC) plus the network (Figure 2), provides an ideal opportunity to refine Puget Lowland, and the Cascade province on the east ESE 22 - 4 VAN WAGONER ET AL.: CRUSTAL STRUCTURE AND EARTHQUAKES (Figure 1) [Brandon and Calderwood, 1990; Snavely and fication of faults with possible reverse offsets [e.g., Johnson Wells, 1996]. et al., 1994, 1999; Brocher et al., 2001] is consistent with [6] In the OSC, metamorphosed Tertiary sedimentary the focal mechanisms of crustal earthquakes [Ma et al., rocks from the accretionary prism have been uplifted since 1996]. Northward compression is most likely the result of the Miocene to form the Olympic mountains [Tabor and oblique subduction [Wang, 1996], possibly in conjunction Cady, 1978a; Brandon et al., 1998]. Within the eastern with northward migration of the Cascadia forearc [Brandon Olympics is the steeply dipping Crescent Formation—thick and Calderwood, 1990; Wells et al., 1998]. submarine and subaerial Paleocene to Eocene basalts which [10] Geologic identification of shallow faults is hampered may have formed by extrusion associated with extension by thick glacial sediments covering much of the Lowland along the margin of north America [Wells et al., 1984; and diffuse crustal seismicity does not reveal distinct Babcock et al., 1992].