Segmentation in Episodic Tremor and Slip All Along Cascadia
|
|
- Pearl McDaniel
- 5 years ago
- Views:
Transcription
1 Segmentation in Episodic Tremor and Slip All Along Cascadia Michael R. Brudzinski 1, Richard M. Allen 2 1 Miami University, Oxford, OH (brudzimr@muohio.edu) 2 University of California, Berkeley, CA (rallen@berkeley.edu) Submitted to Geology, January 24, 2007; Revised May 2, 2007 The recent discovery of Episodic Tremor and Slip (ETS) in subduction zones is based on slow slip episodes visible in GPS observations correlated with non-volcanic tremor signals on seismometers. ETS occurs just inboard from the region capable of great megathrust earthquakes, however, whether there is any communication between these two processes remains unknown. In this study we use new singlestation methods to compile an ETS catalog for the entire Cascadia subduction zone and compare the patterns with a variety of along-strike trends for the subducting and overriding plates. Correlated ground vibrations and strain observations are found all along the subduction zone, demonstrating ETS is an inherent part of the subduction process. There are 3 broad ( km), coherent zones with different recurrence intervals (14±2, 19±4, 10±2 months), where the interval duration is inversely proportional to upper plate topography and the spatial extent correlates with geologic terranes. These zones are further divided into segments of ETS that occur at times typically offset from each other. The 7 largest ( km) segments appear to be located immediately landward from fore-arc basins interpreted as manifestations of megathrust asperities, implying there is a spatial link between ETS and earthquake behavior. It is not yet clear if any temporal link exists, but the regional hold time between ETS episodes could be controlled by strength variations due to composition of geologic terranes. As oceanic plates subduct down into the mantle, friction on the interface with the overriding plate causes stick-slip behavior in the megathrust zone, pulling the upper plate down until it pops back up during a potentially devastating earthquake (Ruff and Kanamori, 1983). Recent observations have also revealed slow slip episodes (SSE) that occur regularly on parts of the deeper plate interface with motion indicating release of accumulated strain (Dragert et al., 2001; Lowry et al., 2001; Ozawa et al., 2002). Their frequency and amount of slip (Mw ~ ) (Kostoglodov et al., 2003; Melbourne et al., 2005) imply they are a substantial portion of the interplate deformation budget. The duration of these episodes is much greater than earthquakes, yet they are accompanied by weeks of non-volcanic tremor (NVT) (Obara, 2002; Rogers and Dragert, 2003; Szeliga et al., 2004). As such, they represent another section of the strain rate continuum between earthquake and geologic time scales. Processes that govern ETS or potential relationships to major earthquakes and local geology remain unknown, although ETS has been proposed to impact the likelihood of megathrust earthquakes (Mazzotti and Adams, 2004). Previous observations of ETS in Cascadia have focused on southern Vancouver Island, northern Washington, and northern California given the density of geophysical observatories (Fig. 1). When observed, SSE representing weeks of transient 1
2 displacements are visible on GPS stations within a few hundred kilometers of the trench and show westward motions back toward the trench (Dragert et al., 2001; Melbourne et al., 2005). Corresponding NVT composed of relatively small, non-earthquake signals are most prominent in the 1-10 Hz frequency band and appear to occur in a region from the plate interface into the upper plate (Kao et al., 2005; McCausland et al., 2005; Obara, 2002). ETS was first established on southern Vancouver Island with regular recurrence at ~14 months (Fig. 2) (Rogers and Dragert, 2003). Prevalent instrumentation in Northern California has detected similar ETS signals at the other end of the subduction zone, but revealed a shorter recurrence interval of ~11 months (Fig. 2) (Szeliga et al., 2004). Previous observations from Oregon are sparse, with a few seismic stations in southern Oregon showing a period of NVT in 2005 (McCausland et al., 2005), and one GPS station in central Oregon showing SSE with a longer recurrence interval (Szeliga et al., 2004). This initial set of observations raises some important questions: 1) Is ETS a phenomenon present throughout the entire subduction zone and thus an inherent part of the subduction process? 2) What generates the spatial variability in recurrence rates of ETS despite their temporal regularity in a given region? 3) Is ETS linked in any way to properties of the megathrust seismogenic zone? Each of these questions requires a broad set of observations in space and time. We utilize a new set of ETS information generated by automated identification of NVT and SSE at individual GPS and seismic stations that circumvents the need for dense networks. These methods are not replacements for network solutions, they are simply used as surrogates to perform a uniform investigation of ETS over the entire subduction zone while network density is still quite heterogeneous. For GPS data, we identify SSE by applying a hyperbolic tangent fit (Larson et al., 2004; Lowry et al., 2001) over a scrolling window of the publicly available PANGA time series. An f-test confirms when the fit is better than a linear fit at 99% confidence, and a threshold value for the transient displacement marks events that are larger than background noise. For seismic data, we calculate a time series from the mean amplitude of filtered envelope seismograms for each nighttime hour recorded at individual stations, instead of using a station network to judge hours when NVT is present (Obara, 2002; Rogers and Dragert, 2003). After a moving average and normalization, large peaks rising above background noise at 99% confidence are identified as periods during which NVT dominates. See data repository for details on automated data processing. Figure 2 summarizes results from 7 different segments throughout Cascadia. For each region a typical seismic trace indicating the times of NVT, and a GPS time series showing the times of SSE are plotted. In regions with previous ETS estimates (McCausland et al., 2005; Rogers and Dragert, 2003; Szeliga et al., 2004), our results are in good agreement with timing of ETS (e.g., Fig. 2C and data repository). More important, NVT that correlate with SSE are apparent in several new locations along the subduction zone, particularly along central Cascadia despite more limited observatories (Fig. 2D F). Corresponding seismic and GPS data availability ranges from 1 to 8 years, with 30 stations reporting SSE and 55 stations reporting NVT (Fig. 1), of the over 300 stations that have been investigated with our automated techniques (Brudzinski and Allen, 2006; Holtkamp et al., 2006, ms. in prep.). It is clear that ETS occurs along the entire subduction zone, meaning that localized geological conditions special to a particular site are not controlling factors that prohibit ETS. This finding is supported by a 2
3 growing number of observations in other subduction zones showing NVT source locations and displacements from SSE over broad regions (Brudzinski et al., 2007; Hirose and Obara, 2006), ruling out hypotheses that require particular along-strike variations to produce ETS (Mitsui and Hirahara, 2006). While ETS is observed throughout the Cascadian subduction zone, the characteristics vary coherently along-strike revealing clear segmentation in the recurrence interval and relative timing of ETS events. First, there are 3 broad geographic zones with different recurrence intervals of ETS (Fig. 1). For example, in the Siletzia Zone in the central part of Cascadia, data from COR extend the history of NVT back by a factor of two to 1989 (Fig. 2E), giving a robust estimate of 19±3 months for the recurrence interval at this station. The average interval across the Siletzia Zone (19±4) is longer than those observed on Vancouver Island to the north (14±2) and is nearly twice as long as that from California to the south (10±2). The broader geographic extent of our ETS measurements relative to previous studies allows us to identify that a coherent Wrangellia Zone extends from northern Vancouver Island down to ~47.5 N, and that a Klamath Zone extends up from the southern end of the subduction zone to ~42.8 N (Fig. 1). This pattern of recurrence intervals is not tied to the overall rate of subduction which drives the earthquake cycle as a whole. Overall convergence velocities decrease slowly from the north to the south (Fig. 1) (DeMets et al., 1990), while the longest recurrence interval occurs in the middle of the subduction zone. Even if one considers the Oregon forearc as a separate rotating microplate (e.g., Wells et al., 1998), the impact on convergence rates would be gradual and not aid in matching the coherent patterns in recurrence intervals. We also find the 3 zones of relatively uniform recurrence intervals cannot be explained by age of the subducting plate, implying along strike variations in ETS are not due to temperature changes. Since age of oceanic lithosphere is proportional to distance from the mid-ocean ridge and the thermal state of oceanic plate is proportional to the square-root of its age (Parsons and Sclater, 1977), one would look for ETS recurrence to be related to distance from the ridge. However, observations along Vancouver Island and northernmost Washington run nearly perpendicular to the Explorer ridge, with station-toridge distances varying from km (c.f., ~500 km in central Oregon). This means the age of oceanic plate beneath the Wrangellia zone should vary rapidly, causing both warmer and colder plate than that below central Cascadia where the longest recurrence dominates. Also, boundaries between the 3 main ETS zones occur ~100 km north of expected age/temperature boundaries in the subducting plate if the existing fracture zones are extended beneath the continent. We suggest that the recurrence interval of ETS is related to properties of the overriding continental plate instead of the subducting oceanic plate. Age and temperature of the subducting plate likely has some impact on generating ETS, because initial work has shown SSE and/or NVT are prominent in other young, warm subduction zones like southwest Japan and Mexico (Hirose and Obara, 2006; Larson et al., 2004; Lowry et al., 2001; Obara, 2002). Yet the oceanic plates subducting beneath Cascadia are relatively uniform compared to the heterogeneity of the continental plate they dive beneath. In fact, the central Siletzia Zone with an ~18 month recurrence interval corresponds to the relatively low lying and young Coastal Range Block of central and northern Oregon and southern Washington (mostly thick Siletzia terrane) (Trehu et al., 1994). The shorter- 3
4 recurrence interval zones to the north (Wrangellia) and south (Klamath) correspond to older Pre-Tertiary blocks with higher topography consisting of a mélange of old oceanic material with later silicic intrusion in a continental environment (Harden, 1998; Jones et al., 1977). Figure 3A shows how ETS recurrence intervals are inversely proportional to onshore fore-arc topography. Correlation of these continental blocks with along-strike patterns of ETS is also consistent with the observation that NVT appears to occur throughout the continental crust at depths above the interface with the subducting oceanic crust (Kao et al., 2005). The zones of spatially coherent recurrence intervals (Fig. 1) are further divided into segments where individual events recur over roughly the same location. While the average recurrence intervals of ETS are similar within a given zone, the relative timing between ETS events shows variation with location (Fig. 2), a phenomenon that is particularly clear when comparing northern and southern Vancouver Island (Fig. 2A, 2C) (Dragert, 2004). The extent of these segments is now emerging from the increased number of ETS observations. Figure 3B illustrates the phase shift in time between different segments by displaying the timing of ETS observations all along Cascadia, with horizontal lines estimating the along-strike extent of a given episode. Since these are station locations instead of source locations, we would expect the grey lines to extend on the order of 50 km beyond the actual source locations. Dashed vertical lines are approximate boundaries defined by events on either side that are separated in time by over a month for greater than 50% of the episodes. We find 7 large segments with alongstrike widths of km (Fig. 3B). The long-term timing of events for each segment is also illustrated by representative GPS and seismic data pairs (Fig. 2 and data repository). Close inspection of NVT time series for the 2003 ETS event that appears to cross segment boundaries reveals distinct offsets in time of 1-2 weeks precisely at the proposed boundaries (Fig. 2 inset and data repository). The largest segments of ETS occur immediately landward from the proposed locations of asperities on the Cascadia megathrust (Wells et al., 2003). The asperity locations are based on large, low gravity, sedimentary basins in the forearc that have been interpreted to indicate potential seismogenic segmentation at depth. This inference is from global surveys finding most of a megathrust earthquake s seismic moment and area of high coseismic slip (asperities) occur beneath the deep-sea forearc basin features from gravity, bathymetry, and seismic profiling (Song and Simons, 2003; Wells et al., 2003). Figure 3B shows the along-strike pattern of prominent fore-arc basins for comparison with the spatial extent of ETS segments (We find a stronger correlation using locations of fore-arc basins based on bathymetry than on gravity). The apparent correlation between segmentation of the seismogenic zone and segmentation of the ETS zone suggests that effects of locking (or lack thereof) on the megathrust are transmitted to greater depths where slow slip is believed to occur (Dragert et al., 2001). This spatially links megathrust structure and anticipated seismogenic behavior with ETS characteristics. A remaining question is whether upper plate structure controls plate interface behavior or vice versa. Both models have been proposed for fore-arc basins, with either basins developing in response to locking on the subduction interface (Song and Simons, 2003; Wells et al., 2003) or thickness of the upper plate critical wedge controlling the frictional behavior on the plate interface (Fuller et al., 2006). For ETS recurrence, the accreted terranes comprising the upper plate above ETS generate inherently sizable 4
5 along-strike variations in structure, composition and age that are presumably more significant than long-term effects of ETS on upper plate structure. This supports an interpretation where variations in the Wrangellia, Siletzia, and Klamath blocks control behavior of the ETS source zone. A clue to how continental blocks could be responsible for differences in ETS recurrence is geochemical evidence that the different terranes have different fluid content (Schmidt and Grunder, 2006), which is a possible catalyst for ETS (Kodaira et al., 2004; Obara, 2002). An intriguing hypothesis is that different terrane composition affects rheology of the upper plate and hence the plate interface (Kohlstedt et al., 1995). For example, the Siletzia terrane would represent denser, stronger, more oceanic-like crust, while the Klamath terrane represents lighter, weaker, more continental-like crust. Such a scenario would suggest that the low-lying Siletzia region has a longer recurrence interval because the upper plate has the strength to accumulate strain for longer periods between SSE. Although it is not yet clear whether rate- and state-dependent friction processes are the best way to explain ETS (Chen and Brudzinski, 2006), initial laboratory fault-sliding experiments suggest variable fluid pressure and rock composition would both be expected to generate coherent variations in recurrence intervals of transient slip (Liu and Rice, 2005, 2006). While it is well established that properties of the subducting plate play key roles in determining plate interface behavior, our result adds to growing evidence that the overriding plate is equally important in megathrust and ETS characteristics. More characterization of the upper plate geologic framework will be essential to assess the nature of deformation at convergent margins. Acknowledgments We relied heavily on data from PANGA, IRIS, USGS, CNSN, PNSN, BDSN, and EarthScope. We benefited from discussions/feedback from R. Blakely, W.-P. Chen, B. Currie, K. Creager, C. DeMets, H. Deshon, W. Hart, S. Holtkamp, H. Kao, A. Lowry, W. McCausland, T. Melbourne, T. Niemi, W. Szeliga, B. Tikoff, C. Thurber, R. Wells, and an anonymous reviwer. NSF provided support for M.B. (EAR , EAR ) and R.A. (EAR , EAR ). UW-Madison also provided support for M.B. and deployment of the OATS array. 5
6 Figure 1. Map illustrating patterns in ETS along the entire Cascadia subduction zone. Colored basemap shows topography and bathymetry. Dashed line onshore marks 40 km depth contour of the subduction interface. Arrows and associated annotations show directions and speeds of subduction relative to North America. Locations of continuous GPS stations (squares) and broadband seismometers (triangles) which exhibit ETS are shown, with colors indicating the recurrence interval when multiple ETS events were observed. Recurrence intervals establish 3 zones that are labeled based on the continental terrane block they associate with. 6
7 Figure 2. Summary of SSE and NVT measurements that characterize seven segments along Cascadia. Colored circles show processed east displacement of GPS monument and black curves are normalized, mean seismic amplitudes, where high spikes are bursts of NVT activity. Red rectangles and black arrows mark the peak of NVT and center of SSE, respectively, with clear temporal correlation in any given segment. GPS data are colored by recurrence interval of SSE as in Figure 1. Inset shows time offsets in NVT activity for the 2003 ETS event from (A) and (B), plotting stations according to distance from the Oregon-Washington border (see data repository for details). 7
8 Figure 3. Plot of along strike patterns of ETS and upper plate features. (A) Top panel shows distinct variations in ETS recurrence with symbols as in Fig 1 and vertical bars show 1σ of observed intervals. Bottom panel shows topography above the 40 km depth contour of the subduction interface, in the middle of ETS observations. Topography is inversely correlated with ETS recurrence, roughly matching the primary continental terranes of different age and composition (Wrangellia, Siletzia, Klamath). (B) Top panel shows "phase" of ETS for 7 different segments along the subduction zone from ETS timing at individual stations. Horizontal grey lines connect stations that record ETS within a month of one another. Bottom panel shows color shaded gravity anomalies and locations of offshore fore-arc sedimentary basins (white lines), features which have been correlated with megathrust asperities on the subduction interface in recent global studies (Wells et al., 2003). Vertical dashed lines show apparent edges of ETS segmentation from currently available data that seem to correlate with megathrust segmentation from the 5 largest sedimentary basins. To deal with trench curvature, station latitudes are those when projected on to 40-km contour (black curve). 8
9 References Brudzinski, M., Cabral, E., DeMets, C., Márquez-Azúa, B., and Correa-Mora, F., 2007, Multiple slow slip transients along the Oaxaca subduction segment from : Geophysical Journal International, p. in review. Brudzinski, M. R., and Allen, R. M., 2006, Segmentation in Episodic Tremor and Slip All Along Cascadia: Eos trans. AGU, v. 87, p. Abstract T53G-05. Chen, W.-P., and Brudzinski, M. R., 2006, Repeating Earthquakes, Episodic Tremor and Slip: Emerging Patterns in Complex Earthquake Cycles?: Complexity, v. in press. DeMets, C., Gordon, R. G., Argus, D. F., and Stein, S., 1990, Current plate motions: Geophysical Journal International, v. 101, p Dragert, H., G.C. Rogers, J. Cassidy, H. H. Kao, and K. Wang, 2004, Episodic tremor and slip in northern Cascadia: USGS Progress Report, v. 04HQGR0047, p Dragert, H., Wang, K. L., and James, T. S., 2001, A silent slip event on the deeper Cascadia subduction interface: Science, v. 292, p Fuller, C. W., Willett, S. D., and Brandon, M. T., 2006, Formation of forearc basins and their influence on subduction zone earthquakes: Geology (Boulder), v. 34, p Harden, D., 1998, California Geology, Prentice Hall, 479 p. Hirose, H., and Obara, K., 2006, Short-term slow slip and correlated tremor episodes in the Tokai region, central Japan: Geophysical Research Letters, v. 33, p. 5. Holtkamp, S., Brudzinski, M. R., and DeMets, C., 2006, Determination of Slow Slip Episodes and Strain Accumulation along the Cascadia Margin: Eos trans. AGU, v. 87, p. Abstract T41A Jones, D. L., Silberling, N. J., and Hillhouse, J., 1977, Wrangellia; a displaced terrane in northwestern North America: Canadian Journal of Earth Sciences, v. 14, p Kao, H., Shan, S. J., Dragert, H., Rogers, G., Cassidy, J. F., and Ramachandran, K., 2005, A wide depth distribution of seismic tremors along the northern Cascadia margin: Nature, v. 436, p Kodaira, S., Iidaka, T., Kato, A., Park, J. O., Iwasaki, T., and Kaneda, Y., 2004, High pore fluid pressure may cause silent slip in the Nankai Trough: Science, v. 304, p Kohlstedt, D. L., Evans, B., and Mackwell, S. J., 1995, Strength of the Lithosphere - Constraints Imposed by Laboratory Experiments: Journal of Geophysical Research-Solid Earth, v. 100, p. 17,587-17,602. Kostoglodov, V., Singh, S. K., Santiago, J. A., Franco, S. I., Larson, K. M., Lowry, A. R., and Bilham, R., 2003, A large silent earthquake in the Guerrero seismic gap, Mexico: Geophysical Research Letters, v. 30. Larson, K. M., Lowry, A. R., Kostoglodov, V., Hutton, W., Sanchez, O., Hudnut, K., and Suarez, G., 2004, Crustal deformation measurements in Guerrero, Mexico: Journal of Geophysical Research-Solid Earth, v. 109, p. B04409, doi: /2003jb Liu, Y. J., and Rice, J. R., 2005, Aseismic slip transients emerge spontaneously in threedimensional rate and state modeling of subduction earthquake sequences: Journal of Geophysical Research-Solid Earth, v Liu, Y. J., and Rice, J. R., 2006, Recurrence interval and magnitude of aseismic deformation transients: an investigation using rate- and state-dependent friction: EOS Trans. Am. Geophys. Union, v. 87, p. T41A Lowry, A. R., Larson, K. M., Kostoglodov, V., and Bilham, R., 2001, Transient fault slip in Guerrero, southern Mexico: Geophysical Research Letters, v. 28, p
10 Mazzotti, S., and Adams, J., 2004, Variability of near-term probability for the next great earthquake on the Cascadia subduction zone: Bulletin of the Seismological Society of America, v. 94, p McCausland, W., Malone, S., and Johnson, D., 2005, Temporal and spatial occurrence of deep non-volcanic tremor: From Washington to northern California: Geophysical Research Letters, v. 32. Melbourne, T., Szeliga, W. M., Miller, M. M., and Santillan, V. M., 2005, Extent and duration of the 2003 Cascadia slow earthquake: Geophysical Research Letters, v. 32, p. L04301, doi: /2004gl Mitsui, N., and Hirahara, K., 2006, Slow slip events controlled by the slab dip and its lateral change along a trench: Earth and Planetary Science Letters, v. 245, p Obara, K., 2002, Nonvolcanic deep tremor associated with subduction in southwest Japan: Science, v. 296, p Ozawa, S., Murakami, M., Kaidzu, M., Tada, T., Sagiya, T., Hatanaka, Y., Yarai, H., and Nishimura, T., 2002, Detection and monitoring of ongoing aseismic slip in the Tokai region, central Japan: Science, v. 298, p Parsons, B., and Sclater, J. G., 1977, An analysis of the variation of ocean floor bathymetry and heat flow with age: J. Geophys. Res., v. 82, p Rogers, G., and Dragert, H., 2003, Episodic tremor and slip on the Cascadia subduction zone: The chatter of silent slip: Science, v. 300, p Ruff, L. J., and Kanamori, H., 1983, Seismic coupling and uncoupling subduction zones: Tectonophysics, v. 99, p Schmidt, M. E., and Grunder, A. L., 2006, Segmentation of the Cascade Arc Based on Compositional and Sr and Nd Isotopic Variations in Primitive Volcanic Rocks: EOS Trans. Am. Geophys. Union, v. 87, p. T53G-02. Song, T. R. A., and Simons, M., 2003, Large trench-parallel gravity variations predict seismogenic behavior in subduction zones: Science, v. 301, p Szeliga, W., Melbourne, T. I., Miller, M. M., and Santillan, V. M., 2004, Southern Cascadia episodic slow earthquakes: Geophysical Research Letters, v. 31, p. L16602, doi: /2004gl Trehu, A. M., Asudeh, I., Brocher, T. M., Luetgert, J. H., Mooney, W. D., Nabelek, J. L., and Nakamura, Y., 1994, Crustal Architecture of the Cascadia Fore-Arc: Science, v. 266, p Wells, R. E., Blakely, R. J., Sugiyama, Y., Scholl, D. W., and Dinterman, P. A., 2003, Basincentered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?: Journal of Geophysical Research Solid Earth, v. 108, p. 2507, doi: /2002jb Wells, R. E., Weaver, C. S., and Blakely, R. J., 1998, Fore-arc migration in Cascadia and its neotectonic significance: Geology (Boulder), v. 26, p
11 Data Repository: Methods for Automated Data Processing and Observations of ETS in Cascadia Michael R. Brudzinski 1, Richard M. Allen 2 1 Miami University, Oxford, OH (brudzimr@muohio.edu) 2 University of California, Berkeley, CA (rallen@berkeley.edu) A.1. Slow Slip Episodes In this section, we discuss the relevant information for identification of slow slip episodes, with additional details presented elsewhere (Holtkamp et al., 2006, ms. in prep.). The GPS data analyzed for slow slip episodes is from the PANGA network provided by the Central Washington University clearinghouse (2005). PANGA time series are network solutions with phase ambiguities resolved and typically consist of one sample per day. Provided time series have been postprocessed with GIPSY, and have linear trends, steps due to earthquakes or hardware upgrades, and annual and semi-annual sinusoids signals simultaneously estimated and removed (Szeliga et al., 2004). In this study, the signal to noise ratio of the GPS data is further enhanced by taking a 10 sample running average of the time series that reduces the scatter among neighboring data points. As shown in previous studies examining transient episodes (Brudzinski et al., 2007; Larson et al., 2004; Lowry et al., 2001), anomalous displacements during slow slip events can be estimated by fitting the GPS coordinate time series with combination of a linear and hyperbolic tangent function. By using a grid search over hyperbolic tangent function parameters (timing and duration of transient), the linear parameters of steady-state velocity and transient displacement are estimated via least squares minimization, weighted by the formal inverse variance of GPS coordinate estimates. This technique has been particularly useful for characterizing the precise timing and overall displacement when a slow slip episode has been identified (Brudzinski et al., 2007; Larson et al., 2004; Lowry et al., 2001). Slow slip episodes are automatically identified by 1) applying the hyperbolic tangent fit over a scrolling window of the time series (6 months before event time and 6 months after) incremented at 0.01 years, 2) using an f-test to confirm when the combined linear and hyperbolic tangent fit is significantly better than a linear fit at 99% confidence within the window, and 3) establishing a threshold value for the transient displacement (~2.0 mm) to identify potential events that are larger than background noise from examination of stations within the stable continental interior (DeMets et al., 2004). Slow slip episodes are verified and event times are selected with an algorithm that first identifies one of the components with an event magnitude that exceeds the displacement threshold and the chi squared value associated with that component must be low. Precise times are selected when the ftest value is near its maximum relative to neighboring times, and the chi squared value is minimized relative to neighboring times. If multiple events within 0.25 years are identified, we take the time with the lowest ratio of the chi squared value to the average chi squared value for all times in the time series. Of these potential events, each one goes through a series of tests (density of data samples, chi squared relative to all time series, and coherence with other neighboring stations) to determine whether it should be included or excluded. Lastly, we perform a manual closer inspection of events, including a handful of extra events that are clear in the time series but only met a smaller displacement threshold (1.5 mm) and excluding a small set 11
12 of events that appear as noise spikes in the time series. Further details of this analysis will be presented in separate work (Holtkamp et al., 2006, ms. in prep.). A.2. Non-volcanic Tremor (NVT) In this section, we discuss the relevant information for identification of episodes of NVT, with additional details presented elsewhere (Brudzinski and Allen, 2006, ms. in prep.). The seismic data for analysis of NVT comes from a variety of networks that spread across the Cascadia subduction zone including the PNSN, USGS-NCSN, CNSN, OATS, IRIS-GSN, BDSN, and the new EarthScope Transportable Array. These networks primarily consist of short period vertical component seismometers, but these are ample instruments to identify NVT which has characteristic frequencies between 1 and 10 Hz (Obara, 2002). Although signals appear to be more pronounced on horizontal components when they are available, vertical components have been sufficient to clearly identify NVT episodes. PNSN data comprise the majority of data analyzed in this study, but continuous data only goes back to 2002, so longer running stations of other networks are critical for identifying longer term trends. Typically a permanent broadband seismic station uses 24-bit digitizers to continuously record all three components of ground velocity at a rate of 20 samples per second. To fully analyze this voluminous dataset which lasts up to 15 years, we use an automated scheme which begins by removing the instrument response and applying a filter with a pass-band of 2-6 Hz. We then take absolute values of the time series and calculate the envelope. Up to this point, our analysis is the same as other researchers (Obara, 2002; Rogers and Dragert, 2003; Szeliga et al., 2004). Next, instead of judging whether tremors are present in each segment of data by the hour, we calculate the mean amplitude of the envelope (Figure A1a). Since various non-tectonic factors can influence the amplitude, we limit our analysis to nighttime hours when culture-generated noise is at a minimum and exclude hours during which signals from known, large earthquakes or other spikes, such as those used for instrument calibration, are present (Figure A1b). To further avoid seismic events, we exclude hours with a high max-to-mean ratio. To help reduce longer-term noise trends that exist at a few stations, we high-pass filter the envelope time series to signals with periods less than about a week (Figure A1c). Lastly, we average amplitudes over a four-day moving window to help accentuate sustained tremor activity, and we normalize to values between 0 and 1 (Figure A1d). In the final time series, large peaks are identified as periods during which tremors dominates. The first cut for identifying NVT periods is by marking peaks that rise above background noise at greater than the 3σ confidence interval. In some cases, analyst intervention is necessary to identify NVT periods due to data loss. To confirm NVT, individual seismograms are investigated as well as mean amplitude time series from neighboring stations. For station VGZ, there is data loss in early 2003 preventing large enough amplitudes to meet the threshold value (Figure A1 and A2), but an NVT period is confirmed by further analysis. We also find that in some cases there are peaks that do not meet the 3σ confidence interval threshold, but the timing corresponds to peaks that meet the threshold at neighboring stations. In these cases, we also find that it is typical for this smaller peak to be accompanied by another smaller peak that is close in time. For these situations, we have added a second criteria that can meet the threshold for an event: at least two successive peaks near in time which rise above background noise at greater than 2σ confidence interval. When calculating the recurrence, we take the mean of the times for the smaller peaks that comprise the event. Our rationale for this 12
13 second criteria is that these cases predominantly occur at stations near the boundary between regions, suggesting they are recording signals from events in neighboring segments. The technique used here has two distinct advantage over previous analyses: 1) it does not requiring the scientist to examine each record to evaluate whether NVT has occurred and 2) it does not require an array of stations to identify NVT through station coherence. The first advantage results in the ability to automate the process and analyze data over large ranges of time and location. Nevertheless, in each case of a large peak, individual seismograms have been examined to confirm the presence of NVT. The second advantage is important for investigating NVT in regions where station density is significantly reduced, like central Cascadia. A.3. Comparison with Previous Studies In order to confirm the validity of the automated seismic and GPS single-station methods, we compare the results of the techniques in this study with previous identifications of ETS. Figure A2 compares the results from a previous for southern Vancouver Island (Rogers and Dragert, 2003) with analysis from this study for stations ALBH and VGZ. Despite the different analysis approaches, there is one-to-one consistency among identified tremor peaks and slow slip episodes between the two studies. This confirms that for single stations with a low enough noise level, the automated techniques produce results very similar to that of previous studies relying on inspection of a network of observations. These methods are not replacements for network solutions, they are simply used as surrogates to perform a uniform investigation of ETS over the entire subduction zone while network density is still quite heterogeneous. We also demonstrate the ability to resolve accurate timing of ETS by illustrating our identification of NVT for the early 2003 event in northern Cascadia that has been characterized in previous papers (e.g., Melbourne, 2004). Figure A3 shows a zoomed in time frame for our seismic time series around the 2003 event, with stations plotted according to distance along strike. Our seismic time series show the same progression to the north and south away from the Puget Sound region seen in the previous GPS analyses, demonstrating that our single station data have the resolution to accurate timing of ETS. Furthermore, our data show that the propagation of ETS is somewhat discontinuous, with changes in timing of the maximum NVT activity by 1-2 weeks around 150 km and 310 km. These locations correspond to the northern and southern boundaries of segment C defined from offsets in timing of several months between other ETS episodes (47.5 and 48.4 in Figure 3 of the main text). These offsets in time during the 2003 event reinforce our proposed segmentation boundaries for ETS. A.4. Observations of Episodic Tremor and Slip in Cascadia To illustrate the new set of ETS observations in Cascadia, we plot time series for a total of 43 GPS and seismic stations all along the subduction zone (Figures A4-A5). These time series were selected to show an equal mixture of slow slip episode and non-volcanic tremor events over an along-strike station spacing of ~25 km. The coherence in timing between stations in each segment provide convincing evidence of common recurrence intervals and segmentation of ETS along strike. Recurrence intervals for each station are calculated from the mean of each interval between times of slow slip or NVT. Then for each segment, we take the mean of all observed intervals for both slow slip and NVT within that segment. Uncertainties on the recurrence shown in Figure 2 are the 1s standard deviation of observed intervals from the mean within that segment. 13
14 We have also constructed a movie of ETS occurrence at stations along Cascadia (Figure A6), to help readers to visualize the segmentation of ETS episodes over time and space. The images in this movie color the station symbols by the time since the last ETS event when data is available. The movie progresses from 1997 to 2006 by 5 day increments, and the time differences between different segments can be seen in when a group of stations reset to a dark red color. Also, one can begin to see the propagation of ETS onsets during a given event from our dataset. Figure A1. Processed time series for seismic activity at station VGZ on southern Vancouver Island illustrating the effects of different components of our processing. (a) Time series are created by taking the mean amplitude of each hour segment of bandpass filtered envelope seismograms (a). Time series are limited to nighttime hours and those without significant earthquake or instrument glitches (b). Time series are high-pass filtered to amplify signals with periods less than about a week (c). Amplitudes are averaged over a four-day moving window to accentuate sustained tremor activity, and normalized to values between 0 and 1 (d). 14
15 Figure A2. Comparison of slow slip and non-volcanic tremor measurements made in this study with previous analysis of episodic tremor and slip. (a) Plot taken from Rogers and Dragert (2003) showing GPS time series for station ALBH relative to North America including both the long-term trend of eastward motions and the slow slip episodes, and hours of tremor activity from a network of seismic stations. (b) Plot showing GPS time series provided by the PANGA network for station ALBH relative to North America with long-term trend removed and stations corrections made. We have further reduced the scatter by using a 10 sample running mean to help accentuate slow slip episodes. Tremor activity is identified for station VGZ by the singlestation method of this study. Despite the different analysis approaches, there is one-to-one consistency among identified tremor peaks and slow slip episodes between the two studies. 15
16 Figure A3. Zoomed view of NVT activity for the 2003 ETS event showing stations from northcentral Cascadia plotted according to their distance along strike from the Oregon-Washington border. Distance is estimated by projecting stations onto the 40 km contour of the slab interface. Vertical bars mark week-long increments. Note that the progression to the north and south from the central Puget Sound region is discontinuous, with jumps in timing around 150 km and 310 km. These locations correspond to the northern and southern boundaries of this segment defined from offsets in timing of several months between other ETS episodes (Figure 3, main text). 16
17 Figure A4. Summary map of seismic (triangles) and GPS (squares) stations in the Cascadia subduction zone whose processed time series and ETS event times are illustrated in Figure A5. The expected source zone is near the dashed 40-km contour of the plate interface. 17
18 18
19 19
20 Figure A5. Examples of slow slip and non-volcanic tremor measurements from stations all along the Cascadia margin. Stations are sorted by latitude from the north to the south, separated into the segments of common ETS timing determined in this study. Station locations are shown in Figure A4. (File: ets c.gif) Figure A6. Movie illustrating ETS occurrence at stations along Cascadia, to help visualize the segmentation of ETS episodes over time and space. Individual images color station symbols by time since the last ETS event, when data is available. The movie progresses from 1997 to 2006 by 5 day increments, and the time differences between different segments can be seen in when a group of stations reset to a dark red color when an ETS event occurs. Also, one can see some indications of the propagation of ETS along strike during a given event as illustrated in Figure A3. 20
21 References Brudzinski, M., Cabral, E., DeMets, C., Márquez-Azúa, B., and Correa-Mora, F., 2007, Multiple slow slip transients along the Oaxaca subduction segment from : Geophysical Journal International, p. in review. Brudzinski, M. R., and Allen, R. M., 2006, Segmentation in Episodic Tremor and Slip All Along Cascadia: Eos trans. AGU, v. 87, p. Abstract T53G-05. Central Washington University, 2005, Geodesy Laboratory & PANGA Data Analysis Facility, DeMets, C., Brudzinski, M., Cabral, E., Márquez-Azúa, B., and Correa-Mora, F., 2004, Largescale seismic and aseismic deformation patterns associated with subduction: Constraints from continuous GPS measurements in Mexico (Abstract): EOS Trans. Am. Geophys. Union, v. 85. Holtkamp, S., Brudzinski, M. R., and DeMets, C., 2006, Determination of Slow Slip Episodes and Strain Accumulation along the Cascadia Margin: Eos trans. AGU, v. 87, p. Abstract T41A Larson, K. M., Lowry, A. R., Kostoglodov, V., Hutton, W., Sanchez, O., Hudnut, K., and Suarez, G., 2004, Crustal deformation measurements in Guerrero, Mexico: Journal of Geophysical Research-Solid Earth, v. 109, p. B04409, doi: /2003jb Lowry, A. R., Larson, K. M., Kostoglodov, V., and Bilham, R., 2001, Transient fault slip in Guerrero, southern Mexico: Geophysical Research Letters, v. 28, p Obara, K., 2002, Nonvolcanic deep tremor associated with subduction in southwest Japan: Science, v. 296, p Rogers, G., and Dragert, H., 2003, Episodic tremor and slip on the Cascadia subduction zone: The chatter of silent slip: Science, v. 300, p Szeliga, W., Melbourne, T. I., Miller, M. M., and Santillan, V. M., 2004, Southern Cascadia episodic slow earthquakes: Geophysical Research Letters, v. 31, p. L16602, doi: /2004gl
Segmentation in episodic tremor and slip all along Cascadia
Segmentation in episodic tremor and slip all along Cascadia Michael R. Brudzinski and Richard M. Allen Geology 35 (10) 907-910, 2007, doi: 10.1130/G23740A.1 Data Repository: Methods for Automated Data
More informationThe Non-volcanic tremor observation in Northern Cascadia. Hsieh Hsin Sung 3/22
The Non-volcanic tremor observation in Northern Cascadia Hsieh Hsin Sung 3/22 Reference Kao, H., S. J. Shan, H. Dragert, and G. Rogers (2009), Northern Cascadia episodic tremor and slip: A decade of observations
More informationAseismic slip and low-frequency earthquakes in the Bungo channel, southwestern Japan
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L769, doi:1.19/3gl19381, Aseismic slip and low-frequency earthquakes in the Bungo channel, southwestern Japan Shinzaburo Ozawa, 1 Yuki Hatanaka, 1 Masaru Kaidzu,
More informationAutomated Detection and Location of Tectonic Tremor along the Entire Cascadia Margin From
*Manuscript Click here to view linked References 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 Automated Detection and Location of Tectonic Tremor along the Entire Cascadia
More informationABSTRACT SPATIAL AND TEMPORAL PATTERNS OF NON-VOLCANIC TREMOR ALONG THE SOUTHERN CASCADIA SUBDUCTION ZONE. by Devin C. Boyarko
ABSTRACT SPATIAL AND TEMPORAL PATTERNS OF NON-VOLCANIC TREMOR ALONG THE SOUTHERN CASCADIA SUBDUCTION ZONE by Devin C. Boyarko Episodic tremor and slip (ETS), the spatial and temporal correlation of slow
More informationDetection and location of non-volcanic tremor beneath the Central Range in Taiwan
Detection and location of non-volcanic tremor beneath the Central Range in Taiwan Aileen Zhang Advisor: Prof. Kazushige Obara Introduction Volcanic tremor is a long-duration episode of weak seismic motions
More informationSeismic and geodetic constraints on Cascadia slow slip
1 2 3 4 5 6 Seismic and geodetic constraints on Cascadia slow slip Aaron G. Wech 1, Kenneth C. Creager 1, & Timothy I. Melbourne 2 1 University of Washington, Department of Earth and Space Science, Box
More informationScaling relationship between the duration and the amplitude of non-volcanic deep low-frequency tremors
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L07305, doi:10.1029/2007gl029391, 2007 Scaling relationship between the duration and the amplitude of non-volcanic deep low-frequency tremors Tomoko Watanabe, 1 Yoshihiro
More informationArray-analysis of Tremors in Shikoku Triggered by the 2012 Sumatra Earthquake
Array-analysis of Tremors in Shikoku Triggered by the 2012 Sumatra Earthquake Tianyi Li 1 Instructor: Prof. Kazushige Obara 2 1. Department of Geophysics, Peking University 2. Earthquake Research Institute,
More informationoverlie the seismogenic zone offshore Costa Rica, making the margin particularly well suited for combined land and ocean geophysical studies (Figure
Chapter 1 Introduction Historically, highly destructive large magnitude (M w >7.0) underthrusting earthquakes nucleate along the shallow segment of subduction zone megathrust fault, and this region of
More informationModeling short and long term slow slip events in the seismic cycles of large subduction earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007566, 2010 Modeling short and long term slow slip events in the seismic cycles of large subduction earthquakes Takanori Matsuzawa, 1 Hitoshi
More informationAn autocorrelation method to detect low frequency earthquakes within tremor
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L16305, doi:10.1029/2008gl034560, 2008 An autocorrelation method to detect low frequency earthquakes within tremor Justin R. Brown, 1 Gregory C. Beroza, 1 and David
More informationNumerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University
More informationDepth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone
2010/10/11-14 Earthscope Workshop Depth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone Kazushige Obara, ERI, Univ. Tokyo Recurrence
More informationSlip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation
Slip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation Toshitaka Baba Research Program for Plate Dynamics, Institute for Frontier
More informationAlong strike variations in short term slow slip events in the southwest Japan subduction zone
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006059, 2010 Along strike variations in short term slow slip events in the southwest Japan subduction zone Shutaro Sekine, 1,2 Hitoshi Hirose,
More informationSendai Earthquake NE Japan March 11, Some explanatory slides Bob Stern, Dave Scholl, others updated March
Sendai Earthquake NE Japan March 11, 2011 Some explanatory slides Bob Stern, Dave Scholl, others updated March 14 2011 Earth has 11 large plates and many more smaller ones. Plates are 100-200 km thick
More informationANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017
ANOTHER MEXICAN EARTHQUAKE! Magnitude 7.1, Tuesday Sept. 19, 2017 Why is there no oceanic crust older than 200 million years? SUBDUCTION If new oceanic crust is being continuously created along the earth
More informationHitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction
Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose
More informationPlate Tectonics. entirely rock both and rock
Plate Tectonics I. Tectonics A. Tectonic Forces are forces generated from within Earth causing rock to become. B. 1. The study of the origin and arrangement of Earth surface including mountain belts, continents,
More informationTABLE S1. CATALOG OF DETECTED EVENTS. Depth (km) M D C? 07/26/11 01:02:
GSA Data Repository 2018361 Templates/ Catalogs* Seamount Events Supplementary material for Newly detected earthquakes in Cascadia subduction zone linked to seamount subduction and deformed upper plate
More informationGlobal Tectonics. Kearey, Philip. Table of Contents ISBN-13: Historical perspective. 2. The interior of the Earth.
Global Tectonics Kearey, Philip ISBN-13: 9781405107778 Table of Contents Preface. Acknowledgments. 1. Historical perspective. 1.1 Continental drift. 1.2 Sea floor spreading and the birth of plate tectonics.
More informationElectrical Conductivity Structures around Seismically Locked Regions
Electrical Conductivity Structures around Seismically Locked Regions Tada-nori Goto Program for Deep Sea Research, IFREE, Japan Agency for Marine-Earth Science and Technology 1. Introduction Existence
More informationAn earthquake like magnitude frequency distribution of slow slip in northern Cascadia
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl044881, 2010 An earthquake like magnitude frequency distribution of slow slip in northern Cascadia Aaron G. Wech, 1,2 Kenneth C. Creager, 1 Heidi
More informationEvolution of seismic signals and slip patterns along subduction zones: insights from a friction lab scale experiment.
Evolution of seismic signals and slip patterns along subduction zones: insights from a friction lab scale experiment. Christophe Voisin, Jean-Robert Grasso, Eric Larose, François Renard * Laboratoire de
More informationDynamic Crust Practice
1. Base your answer to the following question on the cross section below and on your knowledge of Earth science. The cross section represents the distance and age of ocean-floor bedrock found on both sides
More informationScientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures
Scientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures Art Frankel U.S. Geological Survey Seattle, WA GeoPrisms-Earthscope
More informationModulation of Tremor Amplitudes by Tidal Stresses in Cascadia Amanda Klaus Manuscript for Master s project. Abstract
Modulation of Tremor Amplitudes by Tidal Stresses in Cascadia Amanda Klaus Manuscript for Master s project Abstract We recorded the 2010 and 2011 episodic tremor and slip events in Northern Cascadia on
More informationStudy megathrust creep to understand megathrust earthquakes
1 Study megathrust creep to understand megathrust earthquakes Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada, kelin.wang@canada.ca Introduction Once upon a time, there was a belief that
More informationMAR110 Lecture #5 Plate Tectonics-Earthquakes
1 MAR110 Lecture #5 Plate Tectonics-Earthquakes Figure 5.0 Plate Formation & Subduction Destruction The formation of the ocean crust from magma that is upwelled into a pair of spreading centers. Pairs
More information4 Deforming the Earth s Crust
CHAPTER 7 4 Deforming the Earth s Crust SECTION Plate Tectonics BEFORE YOU READ After you read this section, you should be able to answer these questions: What happens when rock is placed under stress?
More informationGeneral Oceanography Geology 105 Expedition 8 Plate Boundaries Beneath the Sea
General Oceanography Geology 105 Expedition 8 Plate Boundaries Beneath the Sea Name Not attempting to answer questions on expeditions will result in point deductions on course workbook (two or more blank
More informationTremor patches in Cascadia revealed by seismic array analysis
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L17316, doi:10.1029/2009gl039080, 2009 Tremor patches in Cascadia revealed by seismic array analysis Abhijit Ghosh, 1 John E. Vidale, 1 Justin R. Sweet, 1 Kenneth
More informationUNIT 6 PLATE TECTONICS
UNIT 6 PLATE TECTONICS CONTINENTAL DRIFT Alfred Wegner proposed the theory that the crustal plates are moving over the mantle. He argued that today s continents once formed a single landmass, called Pangaea
More informationUSArray the first five years
www.earthscope.org USArray the first five years USArray A Continental-Scale Seismic By the Numbers (2003 2008) More than 600 Transportable Array sites have been occupied and more than 535 permits acquired
More informationEffect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological
More informationThe influence of short wavelength variations in viscosity on subduction dynamics
1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.
More informationActivity Pacific Northwest Tectonic Block Model
Activity Pacific Northwest Tectonic Block Model The Cascadia tectonic margin is caught between several tectonic forces, during the relentless motions of the giant Pacific Plate, the smaller subducting
More informationCoseismic slip distribution of the 1946 Nankai earthquake and aseismic slips caused by the earthquake
Earth Planets Space, 53, 235 241, 2001 Coseismic slip distribution of the 1946 Nankai earthquake and aseismic slips caused by the earthquake Yuichiro Tanioka 1 and Kenji Satake 2 1 Meteorological Research
More informationTopic 5: The Dynamic Crust (workbook p ) Evidence that Earth s crust has shifted and changed in both the past and the present is shown by:
Topic 5: The Dynamic Crust (workbook p. 65-85) Evidence that Earth s crust has shifted and changed in both the past and the present is shown by: --sedimentary horizontal rock layers (strata) are found
More informationA) B) C) D) 4. Which diagram below best represents the pattern of magnetic orientation in the seafloor on the west (left) side of the ocean ridge?
1. Crustal formation, which may cause the widening of an ocean, is most likely occurring at the boundary between the A) African Plate and the Eurasian Plate B) Pacific Plate and the Philippine Plate C)
More informationA Brownian walk model for slow earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L17301, doi:10.1029/2008gl034821, 2008 A Brownian walk model for slow earthquakes Satoshi Ide 1 Received 29 May 2008; revised 10 July 2008; accepted 16 July 2008;
More informationPhysics of the Earth and Planetary Interiors 146 (2004) Received 22 November 2003; received in revised form 20 May 2004; accepted 20 May 2004
Physics of the Earth and Planetary Interiors 146 (2004) 513 530 Interplate coupling and a recent aseismic slow slip event in the Guerrero seismic gap of the Mexican subduction zone, as deduced from GPS
More informationOngoing Oroclinal bending in the Cascadia forearc and its relation to concave-outboard plate margin geometry
GSA DATA REPOSITORY 2019058 Ongoing Oroclinal bending in the Cascadia forearc and its relation to concave-outboard plate margin geometry Theron Finley 1,2, Kristin Morell 1,3, Lucinda Leonard 1, Christine
More informationLong-term Crustal Deformation in and around Japan, Simulated by a 3-D Plate Subduction Model
Long-term Crustal Deformation in and around Japan, Simulated by a 3-D Plate Subduction Model Chihiro Hashimoto (1) and Mitsuhiro Matsu ura (2) (1) Institute of Frontier Research for Earth Evolution, Japan
More informationLab 1: Plate Tectonics April 2, 2009
Name: Lab 1: Plate Tectonics April 2, 2009 Objective: Students will be introduced to the theory of plate tectonics and different styles of plate margins and interactions. Introduction The planet can be
More informationPlate Tectonics Unit II: Plate Boundaries (3.5 pts)
T. James Noyes, El Camino College Plate Tectonics Unit II: The Plate Boundaries (Topic 11A-2) page 1 Name: Section: Plate Tectonics Unit II: Plate Boundaries (3.5 pts) Plate Boundaries We will now discuss
More informationGeneral Oceanography Geology 105 Expedition 8 Plate Boundaries Beneath the Sea Complete by Thursday at 11:00 PM
General Oceanography Geology 105 Expedition 8 Plate Boundaries Beneath the Sea Complete by Thursday at 11:00 PM Name Expedition Objectives Learn about the types of plate boundaries and their key characteristics
More informationTremor asperities in the transition zone control evolution of slow earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2012jb009249, 2012 Tremor asperities in the transition zone control evolution of slow earthquakes Abhijit Ghosh, 1 John E. Vidale, 2 and Kenneth
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/4/3/eaao4915/dc1 Supplementary Materials for Global variations of large megathrust earthquake rupture characteristics This PDF file includes: Lingling Ye, Hiroo
More informationLab 2: Plate tectonics
Geology 101 Name(s): Lab 2: Plate tectonics Plate tectonics is the theory that is used to explain geological phenomena worldwide. For this reason, most of the useful maps that illustrate plate tectonics
More informationDynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake
Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 016 Mw 7.0 Kumamoto Earthquake Heng-Yi Su 1 *, Aitaro Kato 1 Department of Earth Sciences, National Central University, Taoyuan
More informationMigration of low frequency tremors revealed from multiple array analyses in western
Migration of low frequency tremors revealed from multiple array analyses in western Shikoku, Japan Tomotake Ueno, Takuto Maeda*, Kazushige Obara, Youichi Asano, and Tetsuya Takeda National Research Institute
More informationSEISMIC CONSTRAINTS ON SLOW SLIP EVENTS WITHIN THE CASCADIA SUBDUCTION ZONE. A Thesis. Presented to. The Graduate Faculty
SEISMIC CONSTRAINTS ON SLOW SLIP EVENTS WITHIN THE CASCADIA SUBDUCTION ZONE A Thesis Presented to The Graduate Faculty Central Washington University In Partial Fulfillment of the Requirements for the Degree
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
More informationOccurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan
LETTER Earth Planets Space, 9, 11 1, Occurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan Shinzaburo Ozawa, Hisashi Suito, and Mikio Tobita Geographical Survey
More informationPresent-day deformation across the southwest Japan arc: Oblique subduction of the Philippine Sea plate and lateral slip of the Nankai forearc
LETTER Earth Planets Space, 55, 643 647, 2003 Present-day deformation across the southwest Japan arc: Oblique subduction of the Philippine Sea plate and lateral slip of the Nankai forearc Takao Tabei 1,
More informationPlate Tectonics: The New Paradigm
Earth s major plates Plate Tectonics: The New Paradigm Associated with Earth's strong, rigid outer layer: Known as the lithosphere Consists of uppermost mantle and overlying crust Overlies a weaker region
More informationSubduction zones are complex plate boundaries in which variable geometry and structure can be
1 Chapter 1 Introduction Subduction zones are complex plate boundaries in which variable geometry and structure can be seismically observed. The along-strike transition from flat to normal subduction is
More informationMoment release rate of Cascadia tremor constrained by GPS
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005909, 2009 Moment release rate of Cascadia tremor constrained by GPS Ana C. Aguiar, 1 Timothy I. Melbourne, 1
More informationPlate Tectonics. Structure of the Earth
Plate Tectonics Structure of the Earth The Earth can be considered as being made up of a series of concentric spheres, each made up of materials that differ in terms of composition and mechanical properties.
More information12. The diagram below shows the collision of an oceanic plate and a continental plate.
Review 1. Base your answer to the following question on the cross section below, which shows the boundary between two lithospheric plates. Point X is a location in the continental lithosphere. The depth
More informationNumerical modelling: The governing equation used in this study is: (K T ) c T H 0,
GSA DATA REPOSITORY 2012254 Cozzens and Spinelli Numerical modelling: The governing equation used in this study is: (K T ) c T H 0, where K is thermal conductivity, T is temperature, ρ is density, c is
More informationNews Release December 30, 2004 The Science behind the Aceh Earthquake
News Release December 30, 2004 The Science behind the Aceh Earthquake PASADENA, Calif. - Kerry Sieh, the Robert P. Sharp Professor of Geology at the California Institute of Technology and a member of Caltech's
More informationOBJECTIVE: For each boundary type, give an example of where they occur on Earth.
OBJECTIVE: Explain the theory of Plate Tectonics. COMPARE AND CONTRAST DIVERGENT, CONVERGENT AND TRANSFORM BOUNDARIES. ***very important. Describe what geologic features form at each of the three CONVERGENT
More informationLETTERS. Non-volcanic tremor driven by large transient shear stresses
Vol 448 2 August 27 doi:1.138/nature617 Non-volcanic tremor driven by large transient shear stresses LETTERS Justin L. Rubinstein 1, John E. Vidale 1, Joan Gomberg 2, Paul Bodin 1, Kenneth C. Creager 1
More informationElectronic supplement for Forearc motion and deformation between El Salvador and Nicaragua: GPS, seismic, structural, and paleomagnetic observations
DR2011053 Electronic supplement for Forearc motion and deformation between El Salvador and Nicaragua: GPS, seismic, structural, and paleomagnetic observations by D. Alvarado et al., Lithosphere, April,
More informationArray analysis and precise source location of deep tremor in Cascadia
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006041, 2010 Array analysis and precise source location of deep tremor in Cascadia Mario La Rocca, 1 Danilo Galluzzo,
More informationKnowledge of in-slab earthquakes needed to improve seismic hazard estimates for southwestern British Columbia
USGS OPEN FILE REPORT #: Intraslab Earthquakes 1 Knowledge of in-slab earthquakes needed to improve seismic hazard estimates for southwestern British Columbia John Adams and Stephen Halchuk Geological
More informationMulti-station Seismograph Network
Multi-station Seismograph Network Background page to accompany the animations on the website: IRIS Animations Introduction One seismic station can give information about how far away the earthquake occurred,
More informationWhole Earth Structure and Plate Tectonics
Whole Earth Structure and Plate Tectonics Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 4/5/2017 14:45 We Discuss Whole Earth Structure and Plate
More informationL.O: THE CRUST USE REFERENCE TABLE PAGE 10
USE REFERENCE TABLE PAGE 10 1. The oceanic crust is thought to be composed mainly of A) granite B) sandstone C) basalt D) rhyolite 2. To get sample material from the mantle, drilling will be done through
More informationPlates & Boundaries The earth's continents are constantly moving due to the motions of the tectonic plates.
Plates & Boundaries The earth's continents are constantly moving due to the motions of the tectonic plates. As you can see, some of the plates contain continents and others are mostly under the ocean.
More informationFull file at
Chapter 2 PLATE TECTONICS AND PHYSICAL HAZARDS MULTIPLE-CHOICE QUESTIONS 1. What direction is the Pacific Plate currently moving, based on the chain of Hawaiian Islands with only the easternmost island
More informationMarine Science and Oceanography
Marine Science and Oceanography Marine geology- study of the ocean floor Physical oceanography- study of waves, currents, and tides Marine biology study of nature and distribution of marine organisms Chemical
More informationEmpirical Green s Function Analysis of the Wells, Nevada, Earthquake Source
Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad
More informationSupplementary Material
1 Supplementary Material 2 3 4 Interseismic, megathrust earthquakes and seismic swarms along the Chilean subduction zone (38-18 S) 5 6 7 8 9 11 12 13 14 1 GPS data set We combined in a single data set
More informationPlate Tectonics. By Destiny, Jarrek, Kaidence, and Autumn
Plate Tectonics By Destiny, Jarrek, Kaidence, and Autumn .The Denali Fault and San Andreas Fault - The San Andreas Fault is a continental transform fault that extends roughly 1300 km (810 miles) through
More informationCreep Events Slip Less Than Ordinary Earthquakes. Emily E. Brodsky 1 and James Mori 2
Creep Events Slip Less Than Ordinary Earthquakes Emily E. Brodsky 1 and James Mori 2 1 Dept. of Earth and Planetary Sciences, UC Santa Cruz, CA, USA 2 Disaster Prevention Research Institute, Kyoto University,
More informationRELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE
RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT
More informationJournal of Geophysical Research Letters Supporting Information for
Journal of Geophysical Research Letters Supporting Information for InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan H. Fattahi 1, F. Amelung
More information6. In the diagram below, letters A and B represent locations near the edge of a continent.
1. Base your answer to the following question on the cross section below and on your knowledge of Earth science. The cross section represents the distance and age of ocean-floor bedrock found on both sides
More informationGeology 300, Physical Geology Spring 2019 Quiz Ch 19, Plate Tectonics Name
Geology 300, Physical Geology Spring 2019 Quiz Ch 19, Plate Tectonics Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) The portion of a fracture
More information10. Paleomagnetism and Polar Wandering Curves.
Map of ocean floor Evidence in Support of the Theory of Plate Tectonics 10. Paleomagnetism and Polar Wandering Curves. The Earth's magnetic field behaves as if there were a bar magnet in the center of
More informationIntroduction To Plate Tectonics Evolution. (Continents, Ocean Basins, Mountains and Continental Margins)
Introduction To Plate Tectonics Evolution (Continents, Ocean Basins, Mountains and Continental Margins) Geo 386 (Arabian Shield Course) Dr. Bassam A. A. Abuamarah Mohanna G E O 3 8 6 A R A B I A N G E
More informationChapter Review USING KEY TERMS. asthenosphere uplift continental drift. known as. tectonic plates move. object. UNDERSTANDING KEY IDEAS
Skills Worksheet Chapter Review USING KEY TERMS 1. Use the following terms in the same sentence: crust, mantle, and core. Complete each of the following sentences by choosing the correct term from the
More informationNorth America subducted under Rubia. Are there modern analogs for Hildebrand s model of North America subducting under Rubia?
North America subducted under Rubia Are there modern analogs for Hildebrand s model of North America subducting under Rubia? In the Geological Society of America Special Papers Did Westward Subduction
More informationPSc 201 Chapter 3 Homework. Critical Thinking Questions
PSc 201 Chapter 3 Homework Critical Thinking Questions 1. (adapted from text) Seawater is denser than fresh water. A ship moving from the Atlantic Ocean into the Great Lakes goes from seawater to fresh
More informationWhat Forces Drive Plate Tectonics?
What Forces Drive Plate Tectonics? The tectonic plates are moving, but with varying rates and directions. What hypotheses have been proposed to explain the plate motion? Convection Cells in the Mantle
More information3. PLATE TECTONICS LAST NAME (ALL IN CAPS): FIRST NAME: PLATES
LAST NAME (ALL IN CAPS): FIRST NAME: PLATES 3. PLATE TECTONICS The outer layers of the Earth are divided into the lithosphere and asthenosphere. The division is based on differences in mechanical properties
More informationPlate Tectonics. Essentials of Geology, 11 th edition Chapter 15
1 Plate Tectonics Essentials of Geology, 11 th edition Chapter 15 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Plate Tectonics: summary in haiku form Alfred Wegener gave us Continental Drift. Fifty years later...
More informationJapan Agency for Marine-Earth Science and Technology (JAMSTEC), 2-15 Natsushima-cho, Yokosuka JAPAN
Dense Ocean floor Network System for Mega Thrust Earthquakes & Tsunamis(DONET) -Towards Understanding Mega Thrust Earthquakes, the Geohazard & Disaster Mitigation- Yoshiyuki KANEDA, Katsuyoshi KAWAGUCHI,
More informationPlate Tectonics: A Scientific Revolution Unfolds
Chapter 2 Lecture Earth: An Introduction to Physical Geology Eleventh Edition Plate Tectonics: A Scientific Revolution Unfolds Tarbuck and Lutgens From Continental Drift to Plate Tectonics Prior to the
More informationSlab morphology in the Cascadia fore arc and its relation to episodic tremor and slip
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006053, 2010 Slab morphology in the Cascadia fore arc and its relation to episodic tremor and slip Pascal Audet,
More informationDOWNLOAD OR READ : SHALLOW SUBDUCTION ZONES SEISMICITY MECHANICS AND SEISMIC POTENTIAL PDF EBOOK EPUB MOBI
DOWNLOAD OR READ : SHALLOW SUBDUCTION ZONES SEISMICITY MECHANICS AND SEISMIC POTENTIAL PDF EBOOK EPUB MOBI Page 1 Page 2 shallow subduction zones seismicity mechanics and seismic potential shallow subduction
More informationSlab pull, slab weakening, and their relation to deep intra-slab seismicity
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L14305, doi:10.1029/2005gl022922, 2005 Slab pull, slab weakening, and their relation to deep intra-slab seismicity Susan L. Bilek Earth and Environmental Science
More informationThe Structure of the Earth and Plate Tectonics
The Structure of the Earth and Plate Tectonics Agree or Disagree? 1. The Earth if made up of 4 different layers. 2. The crust (where we live) can be made of either less dense continental crust or the more
More informationG 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society
Geosystems G 3 AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society Article Volume 10, Number 2 14 February 2009 Q02S12, doi: ISSN: 1525-2027 Click Here for Full Article
More informationWhy Does Oceanic Crust Sink Beneath Continental Crust At Convergent Boundaries
Why Does Oceanic Crust Sink Beneath Continental Crust At Convergent Boundaries What is the process by which oceanic crust sinks beneath a deep-ocean Why does oceanic crust sink beneath continental crust
More informationLearning Objectives (LO)! Lecture 11: Plate Tectonics II! No Homework!! ** Chapter 3 **! What we ll learn today:!
Learning Objectives (LO)! Lecture 11: Plate Tectonics II! No Homework!! ** Chapter 3 **! What we ll learn today:! 1. List the three types of tectonic plate boundaries! 2. Describe the processes occurring
More information