Crust and subduction zone structure of Southwestern Mexico
|
|
- Christiana French
- 5 years ago
- Views:
Transcription
1 PUBLICATIONS Journal of Geophysical Research: Solid Earth RESEARCH ARTICLE Key Points: Geometry of Rivera and Cocos plate from crustal study High V p /V s ratio beneath trenchward edge of magmatism Low shear velocity layer on top subduction plate Crust and subduction zone structure of Southwestern Mexico Sandy Kurniawan Suhardja 1, Stephen P. Grand 1, David Wilson 2, Marco Guzman-Speziale 3, Juan Martin Gomez-Gonzalez 3, Tonatiuh Dominguez-Reyes 4, and James Ni 5 1 Jackson School of Geosciences, University of Texas at Austin, Austin, Texas, USA, 2 U.S. Geological Survey Albuquerque Seismological Laboratory, Albuquerque, New Mexico, USA, 3 Centro de Geociencias, Universidad Nacional Autónoma de México Juriquilla Campus, Queretero, Mexico, 4 Observatorio Vulcanologico, Universidad de Colima, Colima, Mexico, 5 New Mexico State University, Las Cruces, New Mexico, USA Correspondence to: S. P. Grand, steveg@jsg.utexas.edu Citation: Suhardja, S. K., S. P. Grand, D. Wilson, M. Guzman-Speziale, J. M. Gomez- Gonzalez, T. Dominguez-Reyes, and J. Ni (2015), Crust and subduction zone structure of Southwestern Mexico, J. Geophys. Res. Solid Earth, 120, , doi:. Received 3 SEP 2014 Accepted 4 JAN 2015 Accepted article online 9 JAN 2015 Published online 11 FEB 2015 Abstract Southwestern Mexico is a region of complex active tectonics with subduction of the young Rivera and Cocos plates to the south and widespread magmatism and rifting in the continental interior. Here we use receiver function analysis on data recorded by a 50 station temporary deployment of seismometers known as the MARS (MApping the Rivera Subduction zone) array to investigate crustal structure as well as the nature of the subduction interface near the coast. The array was deployed in the Mexican states of Jalisco, Colima, and Michoacan. Crustal thickness varies from 20 km near the coast to 42 km in the continental interior. The Rivera plate has steeper dip than the Cocos plate and is also deeper along the coast than previous estimates have shown. Inland, there is not a correlation between the thickness of the crust and topography indicating that the high topography in northern Jalisco and Michoacan is likely supported by buoyant mantle. High crustal V p /V s ratios (greater than 1.82) are found beneath the trenchward edge of magmatism including below the Central Jalisco Volcanic Lineament and the Michoacan-Guanajuato Volcanic Field implying a new arc is forming closer to the trench than the Trans Mexican Volcanic Belt. Elsewhere in the region, crustal V p /V s ratios are normal. The subducting Rivera and Cocos plates are marked by a dipping shear wave low-velocity layer. We estimate the thickness of the low-velocity layer to be 3 to 4 km with an unusually high V p /V s ratio of 2.0 to 2.1 and a drop in S velocity of 25%. We postulate that the low-velocity zone is the upper oceanic crust with high pore pressures. The low-velocity zone ends from 45 to 50 km depth and likely marks the basalt to eclogite transition. 1. Introduction The tectonic evolution of western North America has been strongly influenced by successive fragmentation events of the ancient Farallon plate, as various segments of the East Pacific rise approached the paleotrench offshore. When a spreading center encounters a subduction zone, a major change in stress and plate boundaries occurs. The change in stress affects the tectonics of the overriding plate which may include unusual volcanism, fragmentation of the overriding plate, and microplate capture [Stock and Lee, 1994]. However, the details of how this occurs and the controlling dynamics are still poorly understood. Currently, the East Pacific Rise is close to encountering the Southwestern Mexico subduction zone near the Colima Rift (Figure 1). In this area the Rivera plate detached from the Cocos plate since 10 Myr [Demets and Traylen, 2000] and is subducting beneath Jalisco, Mexico. The exact boundary between the Cocos and Rivera plates is uncertain although it is likely associated with the El Gordo Graben offshore from the Colima Rift [Bandy et al., 1995]. The oceanic plates subducting have ages between 10 and 13 Myr except for a small section near the El Gordo Graben that may be as young as 4 Myr [Lawver et al., 2013]. Associated with the subduction of the Rivera plate is the overriding Jalisco Block. The Jalisco Block is bounded by the Colima Rift on the east and the Tepic-Zacoalco Rift to the north. Luhr et al. [1985] proposed that the Jalisco Block was produced by an eastward jump of the East Pacific Rise and that the Colima Rift may evolve to an ocean spreading center in the future with the Jalisco Block captured by the Pacific Plate. This model predicts right lateral strike-slip motion along the Tepic-Zacoalco Rift. Rosas-Elguera et al. [1996], however, find evidence of extension along the Tepic-Zacoalco Rift and postulate that the relative motion of the Jalisco Block with respect to the North America plate is due to applied stresses at the plate boundaries as well as the differing motions of the Rivera and Cocos plates, respectively. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1020
2 Figure 1. Map of Jalisco Block (JB) and adjacent regions, showing plate boundaries, volcanoes (open circles with number: 1, Volcan Colima; 2, Navado De Colima; 3, Volcan Cantaro; 4, Sierra La Primavera; 5, Volcan Tequila; and 6, Volcan Ceboruco); the Central Jalisco Volcanic Lineament (CJVL) [Bandy et al., 2001](redtriangles:A:Ayutla,LV:LosVolcanoes,TA:Talpa de Allende, M: Mascota, and SS: San Sebastian, from southeast to northwest). Black circles are MARS stations. Red circles are Colima Volcano Deep Seismic Experiment (CODEX) station. Brown circle are Mexican stations. SCR, Southern Colima Rift; NCG, Northern Colima Graben; TMVB, Trans-Mexican Volcanic Belt; MGVF, Michoacán-Guanajuato Volcanic Field; MAT, Middle American Trench; RI, Rivera plate; CO, Cocos Plate; and PA, Pacific Plate. Widespread magmatism is also present in Southwestern Mexico. Within the Tepic-Zacoalco Rift are several stratovolcanoes that make up the Trans-Mexican Volcanic Belt (Figure 1). An exception to this trend is Colima Volcano that lies closer to the trench. Magmatism has been migrating trenchward during the past 10 Ma [Gómez-Tuena et al., 2007; Ferrari et al., 2012] such that the current front of the magmatic arc is 60 to 80 km south of the Tepic-Zacoalco Rift. In the Jalisco Block, the arc front is marked by a line of volcanic fields known as the Central Jalisco Volcanic Lineament (CJVL). Bandy et al. [2001] have dated volcanic rocks along the CJVL and found a northwest trend of younger rocks ranging from 3 Myr near Ayutla to Holocene ages in the Mascota fields (Figure 1). East of the Jalisco Block is the Michoacan-Guanajuato Volcanic Field (MGVF), a recent volcanic field with more than 1000 Quaternary eruptive centers [Gómez-Tuena et al., 2007]. In spite of the unique tectonic and magmatic activity within Southwestern Mexico, there has been little detailed work done on the crustal and lithospheric structure of the region. A detailed knowledge of crustal structure provides constraints for the dynamic link between the tectonics of the overriding plate and the evolution of the downgoing Rivera and Cocos slabs. Regional surface wave studies [Gaite et al., 2012] provide estimates of crustal thickness although such studies commonly have strong trade-offs between Moho depth and velocity. Gravity studies have been the main source of information on crustal thickness in Southwestern Mexico [Urrutia-Fucugauchi and Flores-Ruiz, 1996] but also rely on estimates of densities of different rock units to determine crustal thickness. In this study we apply receiver function (RF) analysis to a large data set recorded by two temporary seismic arrays deployed in Southwestern Mexico as well as two stations of the Mexican National Network. We use the H-κ stacking method developed by Zhu and Kanamori [2000] to estimate both crustal thickness (H)aswellastheaverageV p /V s ratio (κ) of the crust beneath each station. This approach worked well in the interior of the region where the Moho is relatively flat but did not produce good results near the coast where clear dipping structures are evident. We also use Common Conversion Point (CCP) analysis to improve images near the coast where the dipping oceanic plates produce clear signals in the RF s. The aim of this study is to determine the crustal structure for both continental and oceanic crust beneath Southwestern Mexico. 2. Data Receiver functions were computed from data acquired from the MApping the Rivera Subduction Zone (MARS) array. The MARS array consisted of 50 broadband sensors that were deployed throughout Southwestern Mexico (Figure 1) for 18 months beginning in January Stations were spaced from 35 to 50 km apart with a combination of Strekeisen STS-2 and Guralp 3T and ESP sensors. Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) data acquisition systems recorded the data at 40 sps (samples per second). The MARS experiment involved two U.S. institutions, the University of Texas at Austin and New Mexico State University, in collaboration with two Mexican institutions, Centro de Geociencias, Universidad Nacional Autónoma de México, and the Volcanic Observatory at the Universidad de Colima. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1021
3 Figure 2. Distribution of earthquakes used for the teleseismic receiver function study. Red triangle shows station locations. Black circles show earthquake epicenters. We also use data collected by the Colima Volcano Deep Seismic Experiment (CODEX). CODEX consisted of a deployment of 20 short-period instruments around Colima Volcano (Figure 1) and overlapped the MARS experiment for 5 months in time [Gardine et al., 2007]. Finally, we also used two broadband seismic stations from the Mexican National Seismic network located within Southwestern Mexico. 3. Method Ideally, receiver functions are time series that show pulses corresponding to converted P to S waves recorded by a particular station. Since the conversion of P to S occurs at sharp boundaries in shear velocity in the subsurface, receiver functions can map out subsurface layering beneath seismic stations. A simple approach to isolate the converted S waves on a seismogram is to deconvolve the vertical from the radial component [Langston, 1977; Ammon, 1991]. Data from earthquakes with magnitudes greater than 5.7 and epicentral distances between 30 and 100 were collected from the seismic arrays discussed above. The locations of the hypocenters and the station locations are shown in Figure 2. Most of the good quality data come from Tonga, South America, and the Aleutian islands. Seismic traces were cut from 20 s before the incoming P wave to 100 s after the P wave to ensure all converted phases to a depth of 100 km are included. The data were then rotated to radial and tangential components. A signal-to-noise ratio check was done by comparing the power in the seismic traces 20 s before and after the predicted arrival time of the P wave. Only seismograms with signal-to-noise ratios higher than 2 for both the P and SV components were used to ensure high-quality data before deconvolution. Overall, the signal-to-noise ratio criteria removed 30 40% of all collected traces. For example, at station MA35, from 120 events collected, 44 were removed due to high noise levels. The deconvolution was done in the frequency domain by using the water-level stabilization method and a low-pass Gaussian filter to remove high-frequency noise [Langston, 1977]. The receiver function H(ω) is calculated from the following: RðωÞZ ðωþ HðωÞ ¼ max ZðωÞZ GðωÞ (1) f ðωþ; c*maxðzðwþz* ðwþg GðωÞ ¼ exp ω2 4α 2 (2) where ω is angular frequency, Z(ω) is the Fourier transform of the P component waveform, R(ω)isthetransform of the SV component, and Z*(ω) is the complex conjugate of Z(ω). G(ω) is a Gaussian filter that has zero phase distortion and a lack of sidelobes [Langston, 1977]. It is added to eliminate high-frequency noise in the RF. The frequency content of the Gaussian filter is controlled by the parameter α. Finally, the denominator of equation (1) is either Z(ω)Z * (ω) or a constant c times the maximum of Z(ω)Z * (ω) depending on which value is larger. This is necessary to account for holes in the spectra of Z(ω) and therefore division by a small number which is inherently unstable. The constant c is called the water level and serves to stabilize the deconvolution. The values of α and c were chosen by trial and error where we tried to make the receiver function as sharp as possible but also tried to minimize noise. All receiver functions were computed using a water-level parameter c of and a Gaussian smoothing parameter α of 3.5 resulting in receiver functions with a dominant period near 3 s. A final visual check was also performed. Good receiver functions are identified by having a sharp SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1022
4 Figure 3. Receiver functions for station MA18 with time on the vertical axis and epicentral distance on the horizontal axis. On the right is a linear stack of the receiver functions in the depth domain using a 1-D velocity model to adjust for timing variations due to different incidence angles. P wave signal with little energy arriving earlier. Low-quality RF s tend to have anomalously high-amplitude signals at later times or very wide sidelobes. We eliminated these data before further data processing. An example of receiver functions for station MA18, plotted as a function of epicentral distance, is shown in Figure 3. Most of the receiver functions have similar signals with a peak at 0 s (the P wave) followed by negative side lobes. Almost all RFs from MA18 show a strong positive pulse at about 5 s. After 5 s, there is variability in the RFs, but we can still observe some continuity of arrivals versus distance. Station MA18 is located far from the coast and is expected to have a relatively flat Moho. At this station, we collected 35 high-quality RFs and the clear arrival at about 5 s is interpreted as the P to S wave conversion from the Moho. We also performed a stack of all the MA18 receiver functions in the depth domain. Starting with a 1-D velocity model and horizontal slowness information, each RF can be interpolated into the depth domain to correct for move out and then stacked. On the right side of Figure 3, a stacked RF in the depth domain shows high amplitude at 0 km followed by a strong positive amplitude at 39 km depth interpreted as an arrival from the Moho. Strong signals are also seen at 150 km (16.3 s in the time domain) and 180 km (22 s in the time domain) which are likely PpPs and PsPs/PpSs multiples from the Moho but which could be due to conversions from deeper mantle discontinuities Receiver Function Imaging Ideally, receiver functions show pulses as a function of time that are due to converted P to S waves from interfaces in the subsurface beneath a seismic station. To convert these pulses to depths to interfaces beneath the station requires knowledge of the P and S velocities in the rocks above the interface. Zhu and Kanamori [2000] introduced a method (the H-κ method) that can minimize the ambiguity due to the trade-off between depth and velocity. The time difference between the converted wave and direct P wave depends on the angle of incidence of the incoming P wave, the subsurface P and S velocities and the depth of the interface. The timing of multiple reflected waves (PpPs and PpSs + PsPs) from the conversion depth also depends on incidence angle but differ in their dependence on depth of converter and velocity. Assuming a flat layer over a half-space and a uniform velocity model, the estimated arrival time for the Ps wave as well as the multiples can be calculated. The essence of the H-κ method is that for each station, all receiver functions are summed at the times corresponding to the Ps arrival time as well as the arrival times for the multiples for various choices of thickness H of the layer as well as V p /V s ratio [Zhu and Kanamori, 2000]. For correct values of H and V p /V s ratio, the summation should be a maximum. The approach is essentially a grid search but assumes the structure beneath a station is a single layer (the crust) overlying a half-space. The H-κ method stacks all receiver function amplitudes at one station including the Ps converted wave and its two multiples. The direct Ps converted wave has the largest amplitude and is the least sensitive to dips in the crust-mantle boundary; thus, the Ps converted wave is often given heavier weight in the stacking process. The stacking can be written as follows: SH; ð κþ ¼ X w 1 rt ð 1 Þþw 2 rt ð 2 Þ w 3 rt ð 3 Þ (3) where H is crustal thickness, κ is V p /V s ratio, r i (t) is the ith receiver function at times t 1, t 2, and t 3 which are the predicted times for T Ps, T PpPs, and T PsSs + PsPs. We used a weighting of 0.5, 0.3, and 0.2 for w 1, w 2, and w 3, respectively. Zhu and Kanamori [2000] used 0.7, 0.2, and 0.1 for the weighting in their study of Southern California. We tried several weighting schemes and found a slightly better result with the weights given above. However, overall, we did not see a dramatic change in results using different weighting schemes. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1023
5 Figure 4. (a) Contour of stack amplitudes for station MA18 as a function of crustal thickness H and V p /V s ratio. The grid search calculates a stack amplitude of receiver functions for all ranges of crustal thickness (20 45 km) and V p /V s ratio (1.65 2). The final result is taken by choosing the highest amplitude from the contour, and uncertainty is calculated by measuring the flatness of the contour peak. (b) A plot of all receiver functions from all events for this station The green line is the predicted Ps times assuming a crustal thickness of 39.2 km and a V p /V s ratio of Similarly, the red line is predicted time for PpPs, and the blue line is for PpSs. Crustal thickness and V p /V s ratio of the crust beneath Southwestern Mexico were analyzed using the H-κ stacking method described above. The bounds for the grid search were set to be 20 to 55 km for Moho depth and 1.65 to 2.00 for V p /V s ratio. Figure 4 illustrates the results for station MA18 which had 35 high-quality receiver functions. We contoured the value of the stacked receiver functions as a function of Moho depth and V p /V s ratio with the highest value being the likely true parameters. The contour plot shows a clear maximum with a realistic number for crustal thickness (39.2 km) and crustal V p /V s ratio (1.82). The predicted Moho Ps arrival times agree with the receiver function signals showing a strong positive converted wave at 5 s. Predicted times for the multiples are also plotted. The contour plots using the H-κ method show strong peaks and thus good resolution for most inland stations far from the coast but show poorer results for stations closer to the coast. This is not surprising as it is likely there are strongly dipping structures near the coast where the subducting plate is near the surface. For the coastal stations the peaks in the H-κ contour plots are not sharp and just taking the maximum value from the plots results in unrealistic values for Moho depth and V p /V s ratio. The H-κ method can be extended to dipping layers [e.g., Rossi et al., 2006], but given the complicated nature of the crustal structure near the coast (see below), we did not attempt more sophisticated inversions. For the coastal stations where the H-κ method failed, we used a velocity model with a V p /V s ratio of 1.78, the global crustal average according to Chevrot and van der Hilst [2000], to stack receiver functions in the depth domain. The peak in the stacks was then measured to give an estimate of crustal thickness. Near the coast the first peak in the RFs is usually a negative pulse indicating a drop in velocity with depth. This has been observed in other subduction zones [Bostock, 2013] and has been interpreted as overpressured oceanic basalt crust. Thus, we picked these negative peaks as the Moho. The RF s with initial negative pulses are discussed further below. The contour of the H-κ stack amplitude allows one to calculate uncertainties by measuring the flatness at the maximum point as discussed in Zhu and Kanamori [2000]. For coastal stations the H-κ method did not give satisfactory results. Therefore, to estimate uncertainties for all data we employed a bootstrapping technique [Efron and Tibshirani, 1991] to estimate the 2σ uncertainty for each stacked receiver function. We sampled 100 random populations of receiver functions for each station, for this purpose, and then processed them as discussed in the text. For stations where the H-κ worked, we processed the bootstrapped data with that technique. For stations where the H-κ method did not work, we processed data using a V p /V s ratio of SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1024
6 Figure 5. An interpolation of crustal thickness measurements using the H-κ method and picked Ps arrivals. The measured crustal thickness beneath individual stations is also given with units of kilometers. Crustal thickness varies between ~18 km on the coast and 43 km in land. CO is Colima Volcano, and MGVF is an active magmatic region known as the Michoacán- Guanajuato Volcanic Field. The red triangles show magmatic centers along the Central Jalisco Volcanic Lineament and 1.85 to give uncertainty in Moho depth due to uncertainties in V p /V s ratio. We also performed a bootstrapping analysis with a constant V p /V s ratio and added the two uncertainties together. The results of the RF analysis are shown in Figures 5 and 6 and listed in Table 1 along with the 2σ errors. In Figure 5 point measurements of crustal thickness are interpolated and plotted in map view. Figure 6 shows the interpolated V p /V s ratios in the region where we were able to make reliable measurements. Crustal thickness varies from 20 km along the coast to a maximum of 42 km inland. There is a difference in crustal structure between the Jalisco Block to the northwest and Michoacan to the southeast. Beneath the Jalisco Block, which is underthrust by the Rivera plate, the crust is km thick along the coast and rapidly Figure 6. An interpolated map of V p /V s ratio of the crust as well as the individual measurements. The Colima and Tepic-Zacoalco Rifts are shown as dashed lines. Young volcanoes along the Jalisco Volcanic Lineament are shown as red triangles. CO is Colima Volcano, and MGVF is the Michoacán-Guanajuato Volcanic Field. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1025
7 Table 1. Crustal Thickness and V p /V s Ratio With 2σ Uncertainties a A B Station Name Longitude Latitude Thickness (km) V p /V s Ratio Station Name Longitude Latitude Thickness (km) MA ± ± 0.04 MA ± 1.9 MA ± ± 0.04 MA ± 1.4 MA ± ± 0.02 MA ± 2.2 MA ± ± 0.07 MA ± 3.5 MA ± ± 0.01 MA ± 3.7 MA ± ± 0.04 MA ± 3.0 MA ± ± 0.02 MA ± 2.4 MA ± ± 0.04 MA ± 0.8 MA ± ± 0.09 MA ± 3.6 MA ± ± 0.07 MA ± 5.0 MA ± ± 0.02 MA ± 1.1 MA ± ± 0.16 MA ± 1.1 MA ± ± 0.03 MA ± 2.9 MA ± ± 0.13 MA ± 3.0 MA ± ± 0.10 MA ± 2.1 MA ± ± 0.05 MA ± 1.7 MA ± ± 0.14 MA ± 2.0 MA ± ± 0.06 MA ± 1.2 MA ± ± 0.15 MA ± 2.5 MA ± ± 0.08 MA ± 4.0 MA ± ± 0.03 MA ± 1.1 MA ± ± 0.02 MA ± 1.5 MA ± ± 0.05 MA ± 2.1 MA ± ± 0.19 CJIG ± 0.9 MA ± ± 0.08 MMIG ± 1.0 a For stations where the H-κ method worked, we give both crustal thickness and Vp /V s ratio (A). For the other stations we have no independent control on V p /V s ratio and just give the crustal thickness with associated uncertainties (B). reaches 40 km thick within 50 km of the coast. To the southeast of the Colima Rift, however, the crustal thickness along the coast is km thick and gradually thickens inland reaching 40 km more than 100 km from the coast. In Figure 7 we plot the crustal thickness estimates at each station on a map of topography. Note that there is not a good correlation between crustal thickness and topography. The thickest crust is Figure 7. Map of the topography of Southwestern Mexico with crustal thickness measurements at individual stations in kilometers. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1026
8 roughly in the middle of our study region, parallel to the strike of the subduction zone, but has rather subdued topography. Further inland, the topography increases to near 2 km, yet the crust actually thins a little to km thick. This indicates the high topography is likely supported by buoyant mantle as has been observed in other back-arc regions [Hyndman et al., 2005]. This observation also supports a model of slab rollback and perhaps tearing as proposed by Ferrari [2004] and later by Yang et al. [2009]. In this model hot asthenospheric mantle fills the space previously occupied by slab to relatively shallow depths providing a source of buoyancy in the mantle beneath the high topography. If this interpretation is correct, then the high topography would be relatively recent, i.e., following the time of the rollback estimated to be during the last 5 Myr or so [Ferrari, 2004]. This model is also consistent with observations of current and recent tectonic activity in the northern Jalisco Block including high rock uplift rates determined from stream incision profiles [Castillo et al., 2014], seismic activity in the northern Jalisco Block [Pacheco et al., 1999; Nunez-Cornu et al., 2002], and Pliocene active deformation seen from structural mapping [Ferrari and Rosas-Elguera, 2000]. Another interesting aspect of the crustal thickness map is the lack of significant thinning of crust in the Colima Rift (Figure 5). We do not have the station density to look at the Rift crustal structure in detail, but stations close to the Rift such as MA26 and MA24 do not show thinner crust than adjacent stations further from the Rift. We conclude that there has been little extension in the Rift. This agrees with the conclusion of Rosas-Elguera et al. [1996] that the southern Colima Rift is actually a broad slowly extending zone that has been active only since the late Pliocene. We do see a trend to the north of thinning crust from km thick to km thick. However, given that our deployment did not reach the surface manifestation of the Tepic-Zoacala Rift, it is difficult to place constraints on the overall extension across that Rift. We found the average V p /V s ratio of the crust to be 1.81 (Figure 6) which is slightly higher than the global average of 1.78 [Christensen, 1996; Chevrot and Van der Hilst, 2000]. There is a large range in values, however, from 1.72 to Two regions, in particular, have abnormally high crustal V p /V s ratios. The first region is located in the central part of the Jalisco Block, close to the volcanoes shown as red triangles in Figure 6. Four stations located at the southwestern edge of the region have V p /V s ratios of 1.85 or greater. It is interesting that to the north of these stations the V p /V s ratio is closer to normal. The second region, located to the northeast of Colima Volcano, shows a band of high V p /V s ratios ranging from 1.85 to The eastern part of this band is within the Michoacan-Guanajuato Volcanic Field (MGVF), a wide region of volcanic activity that began about 2.8 Ma and continues today [Gómez-Tuena et al., 2007]. The band of high V p /V s crust extends farther west than the MGVF, toward the northern Colima Rift. Bandy et al. [1995] located the Rivera-Cocos boundary at depth beneath this region based on the occurrence of thermal springs and crustal seismicity in the region although Yang et al. [2009] place the boundary farther west. Average crustal V p /V s can be used to interpret the petrology and physical state of the crust. Christensen [1996] showed from laboratory experiments that V p /V s ratio does not vary much with changes in temperature and pressure for pressures greater than MPa. The main factor that controls the V p /V s ratio in the crust is the presence of melting or fluids and the mineralogy. The relative abundance of quartz and plagioclase feldspar has a dominant effect on V p /V s [Christensen, 1996]: for felsic quartz-rich rocks such as granite, V p /V s is 1.71; intermediate rocks have a V p /V s ratio near 1.78 and mafic plagioclase-rich rocks such as gabbro have a V p /V s ratio near The average composition for continental crust is close to andesite or diorite [Anderson, 1989], and laboratory measurements by Carmichael [1982] confirmed that V p /V s for diorite at crustal pressures ranges from 1.75 to We found the region just to the east of Colima volcano and the northern part of the Jalisco Block have normal continental crust V p /V s ratios (Figure 6). The high V p /V s ratio regions could indicate a very mafic crust or that the crust has high pore pressure fluids or partial melt. The surface geology of the region is largely granitic or siliceous volcanic rock with no indication of a particularly mafic crustal composition [Valencia et al., 2013]. Thus, we conclude that the high V p /V s ratios we observe along two bands in our study area are due to partial melt or high fluid content within the crust. The high V p /V s ratios in the Jalisco Block are just seaward of a trend of volcanoes (red triangles in Figure 6) known as the Central Jalisco Volcanic Lineament (CJVL). It has been known for some time that there has been an overall trenchward migration of magmatism in Jalisco for the past 10 Ma [Ferrari et al., 2001]. The CJVL forms the front of this migrating magmatism. Ferrari [2004] and Yang et al. [2009] proposed that this migration was due to slab rollback of the Rivera plate allowing asthenosphere to warm the lithosphere progressively trenchward. Arc magmas and fluids could also be rising from the deepening slab. The highest SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1027
9 Figure 8. The groups of stations used for the Rivera and Cocos Line RF plots. R-Rʹ is the line used for CCP imaging of the Rivera plate, and C-Cʹ is used for CCP imaging of the Cocos plate. V p /V s ratios we find in the Jalisco Block are just trenchward of the CJVL volcanoes and are not collocated with any volcanoes. This indicates that the crust is being broadly heated, with possible partial melting, or has had the addition of fluids trenchward of the recent volcanism. Thermal modeling of Rivera plate subduction shows dehydration of the slab just trenchward of the volcanic front [Ferrari et al., 2012] supporting the latter model. We further note that magnetotelluric results show high crustal conductivity in this region [Corbo-Camargo et al., 2013] consistent with high crustal fluid content. The second region with anomalously high crustal V p /V s is located in the northeast of our study area partly within the Michoacan-Guanajuato Volcanic Field. This region has late Pliocene-Quaternary mafic and intermediate volcanism less than 3 Myr [Gómez-Tuena et al., 2007]. It is a 40,000 km 2 area with more than 1000 Quaternary eruptive centers, and thus, a crust with distributed partial melt is not surprising [Gómez-Tuena et al., 2007]. We find high crustal V p /V s ratios to the west of the MGVF, as well. There is no volcanism in this region, so high fluid content is likely the cause which is also consistent with the location of the beginning of slab dehydration in the thermal model derived by Ferrari et al. [2012] Cross Sections In order to study the downdip variation of the Rivera and Cocos subduction zone, we plot RF s along two linear profiles relative to an estimate of depth to the subducting Rivera and Cocos plates. We divide the data into two groups (Figure 8) corresponding to the Jalisco Block with the Rivera plate subducting beneath it and stations to the east below which is the subducting Cocos plate. We make this division because of the difference in crustal structure discussed above as well as the difference in slab dip reported by Pardo and Suarez [1995] between the two regions. We use estimates of slab interface depth from Pardo and Suarez [1995]. For each of our stations, we plotted the stacked receiver function as a function of depth to slab (Figure 9). The stacking was done assuming the V p /V s ratio found in our H-κ analysis. For stations that had a null result from the H-κ analysis, we assumed a V p /V s ratio of An estimated depth from H-κ analysis is also placed in the figure. Overall, the stacked RF s show good agreement with the H-κ analysis results although there are slight differences that may be due to a more complex crustal velocity structure than the simple homogenous layer assumed in H-κ analysis. Both panels in Figure 9 show a northerly dipping structure in the southwest of the array that is visible to near 60 km depth whereas in the northeast there is a relatively simple Moho near 40 km depth with no visible SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1028
10 Figure 9. (a) Plot of stacked RF s with respect to depth to the slab for the Rivera Line. (b) Plot of stacked RF s with respect to depth of the slab for the Cocos Line. The vertical axis is depth, and horizontal axis is the depth to the slab taken from Pardo and Suarez [1995]. The dashed red line is our interpretation of the dipping oceanic crust. The red stars are the estimates of continental crust depth using the H-κ method. dipping structures. We plot the RFs to 80 km depth because at deeper depths crustal multiples are visible and make interpretation problematic. We interpret the dipping structures as the subducting Rivera and subducting Cocos plates (Figures 9a and 9b), respectively. Note that for the subducting Cocos plate, the top of the descending plate agrees with the results of Pardo and Suarez [1995] but that we find the Rivera plate to be about 10 km systematically deeper than Pardo and Suarez [1995] claim. This discrepancy is not surprising in that seismicity is sparse in the Rivera subduction zone and instrumentation was lacking in the past. The techniques used above, H-κ analysis and stacking of receiver functions, work well for flat-lying structures but are problematic for dipping structures. Our analyses to this point have assumed all the arrivals from a given boundary at a single station are produced at the same depth. To account for the different locations of conversion points corresponding to different back azimuths to a station, we have also processed our receiver functions using the Common Conversion Point (CCP) stacking method. Individual receiver functions were back projected along the raypath corresponding to the receiver function. The back-projected amplitudes are stacked in lateral and vertical bins resulting in a 3-D image of convertors in the subsurface. This approach has been used in numerous studies [e.g., Dueker and Sheehan, 1997; Schulte-Pelkum et al., 2005] but requires a high density of stations with numerous receiver functions such that each bin is sampled by numerous traces from different directions. The MARS array is relatively sparse, so we collapsed our images onto 2-D lines so as to better compare the results with Figure 10. We divided the data into the same two groups as before (Figure 8) and back projected the receiver functions onto lines RRʹ and CCʹ, respectively. This was done by projecting each convergence point to its correct location and then projecting that location onto the respective lines along a perpendicular to the lines. Figure 10 shows the results of our CCP stacking along the two lines, RRʹ and CCʹ. We used a bin size radius of 12 km for the CCʹ line and 14 km for the RRʹ line. The results can be compared to Figure 9 although it should be noted that we use distance from the trench on the horizontal axis here, but depth to slab was used in Figure 9 as the horizontal axis. The red color in Figure 10 corresponds to a positive pulse in the receiver SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1029
11 functions and thus to a jump in shear velocity with increasing depth. The blue color corresponds to a decrease in shear velocity with depth. The CCP results are similar to the stacking results shown in Figure 9 in that a clear dipping structure is seen in the southwest and a sharp relatively flat continental Moho is seen in the northeast for both profiles. The dipping slab structure is seen to about 50 km depth and then disappears Slab Structure The cross sections shown in Figures 9 and 10 show several interesting features. First, the dipping structure in the southwest is marked by a negative pulse over a positive pulse (blue over red in Figure 10 and dashed lines in Figure 9) indicating that there is a shear velocity drop with depth associated with the subducting slabs underlain by a sharp jump in velocity. Figure 11 shows the Figure 10. (a) The final CCP image result above the Rivera plate line (line raypaths and a hypothetical velocity R-Rʹ in Figure 8). The vertical axis is depth in kilometers, and the horizontal model for this situation, as well as an axis is the distance calculated along the line. Note the starting point is close to the coast. Blue colors show negative receiver function pulses, and red example of a receiver function response. colors show positive receiver function pulses. (b) The final CCP image for Here we assume the dipping crustal the Cocos plate (line C-Cʹ). layer has a low-velocity layer underlain by a jump in velocity to normal crustal velocity and then finally another jump in velocity representing the oceanic Moho. The Ps converted wave going from low velocity to high velocity has a negative polarity. The converted waves from the deeper layers go from higher velocity at depth to lower velocity above and thus have positive polarities. There are also multiple reflections and conversions within the layers. All these waves arrive close in time resulting in a two-sided pulse in the receiver function similar to what we observe in the coastal RF data. One might interpret the slow shear velocity layer to be subducting sediment: however, the gravity study of Manea et al. [2003] suggests that the thickness of the sedimentary column over the Rivera plate probably does not exceed ~20 m and that it gradually increases eastward along the trench. Moreover, the Deep Sea Drilling Project at site 487, located ~11 km offshore of Guerrero State (latitude N and longitude W), found a sedimentary column of ~100 m of Quaternary hemipelagic sediments which overlay ~70 m of late Miocene to Pliocene pelagic. A layer just a few 100 m thick would be too thin to resolve at the periods of the RFs. A slow velocity layer associated with subducting oceanic lithosphere has been observed in several other subduction zones [Bostock, 2013]. Bostock et al. [2002] used a scattered wave inversion technique on seismic data recorded over the Juan de Fuca plate in central Oregon. Their images show a layer of slow velocity that is associated with the subducting plate to a depth near 45 km. Bostock et al. [2002] interpret the slow anomaly as the basaltic oceanic crust subducting beneath continental fore arc. In Figure 9b it can be seen that the RF peak corresponding to the top of the slow velocity occurs at the depth predicted for the subducting plate by Pardo and Suarez [1995] and is even deeper in Figure 9a. Thus, we interpret the drop in velocity found in the receiver functions to mark the top of the subducting oceanic crust in accord with Bostock et al. [2002]. The estimate of thickness of the low-velocity layer depends on the V p /V s ratio of P and S waves within the slow layer. Using normal values of V p /V s ratio, the thickness of the layer in Bostock et al. [2002] corresponds to a layer roughly 8 km thick, and thus, they associated the layer with the entire oceanic basalt crust. However, Audet et al. [2009], using multiple reflections within the slow layer for data collected above Cascadia, found SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1030
12 Figure 11. (a) Synthetic seismogram response for the four-layer velocity model shown in Figure 11b. The red line in Figure 11a shows impulsive spikes corresponding to different converted arrivals, while the blue line shows the corresponding receiver function at the longer periods used in this study. (c) (h) The raypaths for various converted and reflected waves with labels corresponding to the arrivals in Figure 11a. that the V p /V s ratio is abnormally high within the low-velocity layer, 2.4 to 2.8, and thus, the low-velocity layer is thinner than previously thought and is on the order of 3 to 5 km thick. Recently, Kim et al. [2010] and Song et al. [2009] examined data from a dense linear profile across central Mexico. They also found a thin 2 4 km thick layer of anomalously slow S velocity with anomalously high V p /V s ratio dipping beneath the continent. Audet et al. [2009] as well as Kim et al. [2010] and Song et al. [2009] identify the slow layer as the top half of the subducting oceanic crust. To better constrain the structure producing the RFs, we calculated synthetic receiver functions for various models using a reflectivity code [Levin and Park, 1997]. The three parameters we investigated are the drop in shear velocity at the basalt interface, the V p /V s ratio within the low-velocity layer, and the thickness of the layer. The velocity model we used has four layers, a continental crust, a low-velocity layer identified with the subducting oceanic crust, a lower oceanic crust, and the mantle. We constrain the continental crust, lower SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1031
13 Figure 12. The left figures show synthetic receiver function match real stack RF data at different thickness model and the right figures show the optimized velocity model to match the data. oceanic crust, and mantle to have a V p /V s ratio of 1.78 and vary the thickness, shear velocity, and V p /V s ratio of the low-velocity layer to match our data. The velocity of the lower oceanic crust and mantle are taken as normal, 6.8 km/s V p and 3.8 km/s V s for deep ocean crust and 8.0 km/s V p and 4.5 km/s V s for the mantle. Our receiver functions along the coast vary considerably with different amounts of noise such that no single model will fit all the data. For the purposes of waveform modeling, we stacked all coastal RFs with negative first motions near the Moho arrival time with back azimuths near 300. The stacked RF is then an average of converted waves from the subducting Rivera and Cocos plates from one direction. Figure 12 shows the results of our simulations with varying thicknesses of the low-velocity layer compared to the stacked coastal RF. We varied the thickness of the low-velocity zone from 1 to 6 km thick. For each thickness, we adjusted the S velocity drop in the low-velocity zone as well as the V p /V s ratio to best match the stacked data. The S velocity drop controls the amplitude of the first arriving negative pulse, and the V p /V s ratio controls the timing of the following positive pulses. Our optimal model is a 3 km thick layer with a drop in shear velocity of 25% and a V p /V s ratio within the low-velocity layer of , well above the range for normal rocks at this depth. This is similar to the results of Audet et al. [2009] although we clearly do not have tight constraints on the slow velocity layer s properties due to variable receiver functions in our profile as well as the fact that multiple variables affect the results. SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1032
14 Basalt is usually considered a high-velocity component of crust, particularly within continental crust; therefore, it is surprising that it is showing up as a slow velocity anomaly in our study as well as in the work of Bostock et al. [2002]. Audet et al. [2009] studied this issue and concluded that the basalt would have the observed low velocities if it had pervasive water present that was overpressured. They propose that the top of the oceanic crust formed of pillow basalts and sheeted dikes of gabbro contains significant water as well as hydrated minerals. As the plate subducts some of the minerals dehydrate creating high pore fluid pressures that cause the dramatic drop in shear velocity observed. This would imply that the boundary between the subducting oceanic crust and the overriding continental material is impermeable as well as the deeper crust. Kim et al. [2010] and Song et al. [2009] proposed a similar interpretation for their images in central Mexico although Kim et al. [2010] also suggest hydrous minerals such as talc may be an important contributor to the velocity drop. We observe the disappearance of the slow dipping converted phase at about 45 km (Figure 10). We interpret this to be the depth where basalt begins to transform to eclogite. When this occurs, there is a relatively large volume change that may make the subducting crust and its boundary more permeable and thus allow it to lose water to the overlying mantle wedge. If the slow velocity layer is due to high pore pressure fluids, then its disappearance may show the depth at which the water can escape into the mantle wedge. Eclogite also has high seismic velocity, and thus, its contrast with the surrounding mantle will be less than the contrast of the slow basalt with mantle causing the disappearance of the positive converted wave at the oceanic Moho. The slow velocity layer found by Bostock et al. [2002] also disappears at a similar depth. The similarity of subduction zone structure between Cascadia and Southwestern Mexico is striking but both regions have young ocean (5 Myr to 15 Myr) subducting beneath continent. One difference between the regions is the thickness of sedimentary cover in the subduction zone. Offshore Southwestern Mexico estimates of the sedimentary thickness range from 20 to 200 m, as mentioned above. However, the subduction zone offshore Cascadia has a thick sedimentary prism and far more flux of sediment into the subduction zone [Rea and Ruff, 1996]. The similarity of the two regions shows that sediment flux is not a significant factor in the permeability of the oceanic crust interface nor the temperature structure of the slab as the eclogite phase transition seems to occur at similar depths as well. It also supports the model for the slow velocity layer being basalt crust and not a layer of subducting sediment. A second unusual feature of the CCP images of Southwestern Mexico (Figure 10) is the lack of any strong Moho signal in the middle of the images, both for the Rivera plate system and the Cocos plate system. A clear continental Moho is visible to the east and clear dipping convertors are seen to the west, but in the middle there are no clear signals corresponding to the Moho. Above the Rivera slab, we measure a 20 km width of weak Moho signal, and above the Cocos slab, we measure about a 50 km width without a clear Moho. Again, this is similar to what is observed beneath Oregon. Bostock et al. [2002] explained this observation by postulating that the mantle wedge at this location is serpentinized due to the water release from the subducting slab. Serpentine has very low shear velocity and in fact has slower shear velocity than typical lower crustal mineral assemblages [Christensen, 1996]. Bostock et al. [2002] show that a peridotite with 50% serpentization will have shear velocity similar to lower crustal rocks. With no shear velocity contrast between crust and mantle, no converted P to S wave will be created at the Moho. The implication is that this mantle does not participate in the mantle flow associated with the subducting slab and is hydrated by released water from the subducting crust. Our data support this interpretation. It is interesting that the width of the zone with weak Moho is quite different between the Rivera and Cocos subduction zones. The Rivera plate has a steeper subduction angle than the Cocos plate thus creating a thicker mantle wedge as a function of distance inland. A thicker mantle wedge can be entrained in the downward flow easier than a thin wedge. These observations may provide constraints on viscosity of the mantle wedge. 4. Conclusions We have analyzed the lithospheric structure beneath Southwestern Mexico using the receiver function technique. Along the coast, receiver functions show the dipping Rivera and Cocos plates subducting beneath Mexico. We find the Rivera plate along the coast is about 10 km deeper than previously estimated by Pardo and Suarez [1995]. The receiver functions also show that the subducting slabs have a 3 4 km thick extremely slow shear velocity layer that is similar to what has been found in several other subduction zones. The SUHARDJA ET AL American Geophysical Union. All Rights Reserved. 1033
Plate 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 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 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 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 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 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 informationEarthquakes. Earthquakes are caused by a sudden release of energy
Earthquakes Earthquakes are caused by a sudden release of energy The amount of energy released determines the magnitude of the earthquake Seismic waves carry the energy away from its origin Fig. 18.1 Origin
More informationLate 20 th Century Tests of the Continental Drift Hypothesis
Late 20 th Century Tests of the Continental Drift Hypothesis 5 Characteristics of the Ocean Trenches Unless otherwise noted the artwork and photographs in this slide show are original and by Burt Carter.
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 informationSeismic structure beneath the Rivera subduction zone from finite-frequency seismic tomography
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005830, 2009 Seismic structure beneath the Rivera subduction zone from finite-frequency seismic tomography Ting Yang, 1 Stephen P. Grand, 2
More informationUSU 1360 TECTONICS / PROCESSES
USU 1360 TECTONICS / PROCESSES Observe the world map and each enlargement Pacific Northwest Tibet South America Japan 03.00.a1 South Atlantic Arabian Peninsula Observe features near the Pacific Northwest
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 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 informationTopic 12: Dynamic Earth Pracatice
Name: Topic 12: Dynamic Earth Pracatice 1. Earth s outer core is best inferred to be A) liquid, with an average density of approximately 4 g/cm 3 B) liquid, with an average density of approximately 11
More informationDynamic Earth A B1. Which type of plate boundary is located at the Jordan Fault? (1) divergent (3) convergent (2) subduction (4) transform
Dynamic Earth A B1 1. The edges of most lithospheric plates are characterized by (1) reversed magnetic orientation (2) unusually rapid radioactive decay (3) frequent volcanic activity (4) low P-wave and
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 informationTheory of Plate Tectonics
Plate Tectonics Theory of Plate Tectonics Lithosphere is made of sections called plates that move around called continental drift Plates move because they float on the asthenosphere Most plates consist
More informationNC Earth Science Essential Standards
NC Earth Science Essential Standards EEn. 2.1 Explain how processes and forces affect the Lithosphere. EEn. 2.1.1 Explain how the rock cycle, plate tectonics, volcanoes, and earthquakes impact the Lithosphere.
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 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 information1. In the diagram below, letters A and B represent locations near the edge of a continent.
1. In the diagram below, letters A and B represent locations near the edge of a continent. A geologist who compares nonsedimentary rock samples from locations A and B would probably find that the samples
More informationCopyright McGraw-Hill Education. All rights reserved. No reproduction or distribution without the prior written consent of McGraw-Hill Education
Copyright McGraw-Hill Education. All rights reserved. No reproduction or distribution without the prior written consent of McGraw-Hill Education Tibetan Plateau and Himalaya -southern Asia 11.00.a VE 10X
More informationBeneath our Feet: The 4 Layers of the Earty by Kelly Hashway
Beneath our Feet: The 4 Layers of the Earty by Kelly Hashway The Earth is more than a giant ball made up of dirt, rocks, and minerals. The Earth may look like a giant ball from when looking at it from
More informationStructure of the Earth and the Origin of Magmas
Page 1 of 12 EENS 2120 Petrology Tulane University Prof. Stephen A. Nelson Structure of the Earth and the Origin of Magmas This document last updated on 23-Jan-2015 Magmas do not form everywhere beneath
More informationCaptain s Tryouts 2017
Captain s Tryouts 2017 Dynamic Planet Test Written by: Araneesh Pratap (Chattahoochee High School) Name: Date: Answer all questions on the answer sheet. Point values are given next to each question or
More informationPlate Tectonics Tutoiral. Questions. Teacher: Mrs. Zimmerman. Plate Tectonics and Mountains Practice Test
Teacher: Mrs. Zimmerman Print Close Plate Tectonics and Mountains Practice Test Plate Tectonics Tutoiral URL: http://www.hartrao.ac.za/geodesy/tectonics.html Questions 1. Fossils of organisms that lived
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 informationContinental Drift and Plate Tectonics
Continental Drift and Plate Tectonics Continental Drift Wegener s continental drift hypothesis stated that the continents had once been joined to form a single supercontinent. Wegener proposed that the
More informationQuestions and Topics
Plate Tectonics and Continental Drift Questions and Topics 1. What are the theories of Plate Tectonics and Continental Drift? 2. What is the evidence that Continents move? 3. What are the forces that
More informationGLY 155 Introduction to Physical Geology, W. Altermann. Grotzinger Jordan. Understanding Earth. Sixth Edition
Grotzinger Jordan Understanding Earth Sixth Edition Chapter 4: IGNEOUS ROCKS Solids from Melts 2011 by W. H. Freeman and Company Chapter 4: Igneous Rocks: Solids from Melts 1 About Igneous Rocks Igneous
More informationSubduction zones 3 arc magmatism
5. 3 Subduction zones 3 arc magmatism Where can we observe magmatic/volcanic activities along subduction zones? Characteristics of arc magmatism (vs. mid-ocean ridge/intraplate magmatism) Structure of
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 informationPlate Tectonics Practice Test
Plate Tectonics Practice Test 1. What is the main idea Alfred Wegner proposed in the Theory of Continental Drift that he published in 1915? a. The continents float on a liquid layer that allows them to
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 informationChapter 4 Rocks & Igneous Rocks
Chapter 4 Rocks & Igneous Rocks Rock Definition A naturally occurring consolidated mixture of one or more minerals e.g, marble, granite, sandstone, limestone Rock Definition Must naturally occur in nature,
More informationDR
DR2003071 0 0 270 0 30 0 90 0 60 0 120 0 150 0 90 0 180 0 180 0 A) RadialReceiverFunctions B ackazimuth (in degrees relative to north) -135-90 -45 0 45 90 135 180-5.0-2.5 Tangential R eceiver Functions
More informationImaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges
2017 SCEC Proposal Report #17133 Imaging sharp lateral velocity gradients using scattered waves on dense arrays: faults and basin edges Principal Investigator Zhongwen Zhan Seismological Laboratory, California
More informationOCN 201: Seafloor Spreading and Plate Tectonics I
OCN 201: Seafloor Spreading and Plate Tectonics I Revival of Continental Drift Theory Kiyoo Wadati (1935) speculated that earthquakes and volcanoes may be associated with continental drift. Hugo Benioff
More informationChapter 7 Plate Tectonics
Chapter 7 Plate Tectonics Earthquakes Earthquake = vibration of the Earth produced by the rapid release of energy. Seismic Waves Focus = the place within the Earth where the rock breaks, producing an earthquake.
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 informationDynamic Subsidence and Uplift of the Colorado Plateau. Supplementary Material
GSA DATA REPOSITORY 2010177 Liu and Gurnis Dynamic Subsidence and Uplift of the Colorado Plateau Supplementary Material Lijun Liu and Michael Gurnis Seismological Laboratory California Institute of Technology
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 informationPlate Tectonics Lab II: Background Information
Plate Tectonics Lab II: Background Information This lab is based on a UW ESS101 Lab. Note: Hand in only the Answer Sheet at the back of this guide to your Instructor Introduction One of the more fundamental
More information4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution
4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution PI: Lijun Liu Department of Geology, University of Illinois at Urbana-Champaign Corresponding author: Lijun Liu, ljliu@illinois.edu
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
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 informationGeoffrey Abers. Water and volcanism. Imaging. Lessons from Central America. Oceanic. Terrestrial
Imaging Circum-Pacific Subduction Zones with Earthquakes: Fluid Pathways and the Origins of Volcanic Arcs Geoffrey Abers Oceanic Terrestrial Water and volcanism Imaging Lessons from Central America Wet
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 informationThe Dynamic Crust 2) 4) Which diagram represents the most probable result of these forces? 1)
1. The diagrams below show cross sections of exposed bedrock. Which cross section shows the least evidence of crustal movement? 1) 3) 4. The diagram below represents a section of the Earth's bedrock. The
More informationPlate Tectonics 3. Where Does All the Extra Crust Go?
Plate Tectonics 3 Where Does All the Extra Crust Go? Unless otherwise noted the artwork and photographs in this slide show are original and by Burt Carter. Permission is granted to use them for non-commercial,
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 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 informationThe continental lithosphere
Simplicity to complexity: The continental lithosphere Reading: Fowler p350-377 Sampling techniques Seismic refraction Bulk crustal properties, thickness velocity profiles Seismic reflection To image specific
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 informationMoho (Mohorovicic discontinuity) - boundary between crust and mantle
Earth Layers Dynamic Crust Unit Notes Continental crust is thicker than oceanic crust Continental Crust Thicker Less Dense Made of Granite Oceanic Crust Thinner More Dense Made of Basalt Moho (Mohorovicic
More informationCHAPTER 2 THE WAY THE EARTH WORKS: EXAMINING PLATE TECTONICS
CHAPTER 2 THE WAY THE EARTH WORKS: EXAMINING PLATE TECTONICS Coverage of plate tectonics is required early in introductory Geology as context for the mineralogy, petrology, structure, internal processes,
More informationIgneous Rocks. Igneous Rocks. Genetic Classification of
Igneous Rocks Fig. 5.1 Genetic Classification of Igneous Rocks Intrusive: crystallized from slowly cooling magma intruded within the Earth s crust; e.g. granite, gabbro 1 Fig. 5.2 Genetic Classification
More informationReceiver Function (RF) Estimation Using Short Period Seismological Data
Receiver Function (RF) Estimation Using Short Period Seismological Data Shantanu Pandey Department of Earth Sciences, Kurukshetra University, Kurukshetra-136119 Abstract The advent of the receiver function
More informationUnit Topics. Topic 1: Earth s Interior Topic 2: Continental Drift Topic 3: Crustal Activity Topic 4: Crustal Boundaries Topic 5: Earthquakes
The Dynamic Earth Unit Topics Topic 1: Earth s Interior Topic 2: Continental Drift Topic 3: Crustal Activity Topic 4: Crustal Boundaries Topic 5: Earthquakes Topic 1: Earth s Interior Essential Question:
More informationEarth s Interior and Geophysical Properties. Chapter 13
Earth s Interior and Geophysical Properties Chapter 13 Introduction Can we just go there? Deep interior of the Earth must be studied indirectly Direct access only to crustal rocks and upper mantle fragments
More information12.2 Plate Tectonics
12.2 Plate Tectonics LAYERS OF THE EARTH Earth is over 1200 km thick and has four distinct layers. These layers are the crust, mantle (upper and lower), outer core, and inner core. Crust outer solid rock
More informationDivergent Boundaries: Origin and Evolution of the Ocean Floor Earth - Chapter 13 Stan Hatfield Southwestern Illinois College
Divergent Boundaries: Origin and Evolution of the Ocean Floor Earth - Chapter 13 Stan Hatfield Southwestern Illinois College Mapping the Ocean Floor Depth was originally measured by lowering weighted lines
More informationChapter. Graphics by Tasa Graphic Arts. Inc.
Earth Chapter Plate Science 9 Tectonics Graphics by Tasa Graphic Arts. Inc. 1 I. Earth s surface is made up of lithospheric plates. A. Lithospheric plates are composed of the crust and part of the upper
More informationStandard 2, Objective 1: Evaluate the source of Earth s internal heat and the evidence of Earth s internal structure.
Standard 2: Students will understand Earth s internal structure and the dynamic nature of the tectonic plates that form its surface. Standard 2, Objective 1: Evaluate the source of Earth s internal heat
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationFigures S1 S4 show the measurements taken from the synthetic vespagrams where a)
Figures S1 S4 show the measurements taken from the synthetic vespagrams where a) is the differential travel time versus the Dʺ discontinuity height, b) is the differential travel time versus δv S, c) is
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 informationA. rift valley B. deep trench C. volcanic cone D. uplifted mountain. A. a volcanic island arc B. a large watershed area
DC Practice Test Dynamic Crust PRCTICE TEST Zimm's Earth Science 2013 1. Which feature would most likely be forming at converging continental tectonic plate boundaries?. rift valley. deep trench C.
More informationRocky Mountain evolution: Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005726, 2009 Rocky Mountain evolution: Tying Continental Dynamics of the Rocky Mountains and Deep Probe seismic experiments with receiver functions
More informationNAME HOMEWORK ASSIGNMENT #4 MATERIAL COVERS CHAPTERS 19, 20, 21, & 2
NAME HOMEWORK ASSIGNMENT #4 MATERIAL COVERS CHAPTERS 19, 20, 21, & 2 Assignment is due the beginning of the class period on December 14, 2004. Mark answers on a scantron sheet, which will be provided.
More informationOrigin and Evolution of the Ocean Floor
Origin and Evolution of the Ocean Floor Outline Mapping the Ocean Floor Continental Margins Origin of Oceanic Lithosphere Structure of Ocean Crust Mapping the ocean floor Depth originally measured by lowering
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 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 informationEarth. Temp. increases with depth, the thermal gradient is 25 o C/km. Pressure and density also increase with depth.
Plate Tectonics Earth Earth overall average density = 5.5 g/cm 3. Temp. increases with depth, the thermal gradient is 25 o C/km. Pressure and density also increase with depth. Spheroid: with a longer major
More informationLecture 24: Convergent boundaries November 22, 2006
Lecture 24: Convergent boundaries November 22, 2006 Convergent boundaries are characterized by consumption of oceaninc crust in subduction zones, with attendant arc volcanism, metamorphism and uplift.
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 informationEarth Dynamics. Landforms at Plate Boundaries
Earth Dynamics CHAPTER 8 LESSON 2 Landforms at Boundaries What do you think? Read the two statements below and decide whether you agree or disagree with them. Place an A in the Before column if you agree
More informationEarthquake patterns in the Flinders Ranges - Temporary network , preliminary results
Earthquake patterns in the Flinders Ranges - Temporary network 2003-2006, preliminary results Objectives David Love 1, Phil Cummins 2, Natalie Balfour 3 1 Primary Industries and Resources South Australia
More informationD) outer core B) 1300 C A) rigid mantle A) 2000 C B) density, temperature, and pressure increase D) stiffer mantle C) outer core
1. In which area of Earth's interior is the pressure most likely to be 2.5 million atmospheres? A) asthenosphere B) stiffer mantle C) inner core D) outer core Base your answers to questions 2 and 3 on
More informationCrustal Boundaries. As they move across the asthenosphere and form plate boundaries they interact in various ways. Convergent Transform Divergent
Name: Date: Period: Plate Tectonics The Physical Setting: Earth Science CLASS NOTES Tectonic plates are constantly moving and interacting As they move across the asthenosphere and form plate boundaries
More informationThe importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru
Chapter 7 The importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru Abstract Flat subduction near Peru occurs only where the thickened crust of
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 informationEarth Science Lesson 3 Dynamic Earth Processes (Grades 8-12 ) Instruction 3-3 Plate Boundaries
As we told you in a recent Instruction, much of the Earth's lithosphere (its Crust and the top layer of its Outer Mantle) is made up of moving masses of rock called tectonic plates (tectonic means moving).
More informationSeismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan
Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan Thesis by Sara L. Dougherty In Partial Fulfillment of the Requirements for the Degree
More informationPrentice Hall EARTH SCIENCE
Prentice Hall EARTH SCIENCE Tarbuck Lutgens Chapter 9 Plate Tectonics 9.1 Continental Drift An Idea Before Its Time Wegener s continental drift hypothesis stated that the continents had once been joined
More informationPrentice Hall EARTH SCIENCE
Prentice Hall EARTH SCIENCE Tarbuck Lutgens Chapter 9 Plate Tectonics 9.1 Continental Drift An Idea Before Its Time Wegener s continental drift hypothesis stated that the continents had once been joined
More informationPractice Questions: Plate Tectonics
Practice Questions: Plate Tectonics 1. Base your answer to the following question on The block diagram below shows the boundary between two tectonic plates. Which type of plate boundary is shown? A) divergent
More informationGEOLOGY MEDIA SUITE Chapter 12
UNDERSTANDING EARTH, SIXTH EDITION GROTZINGER JORDAN GEOLOGY MEDIA SUITE Chapter 12 Volcanoes 2010 W.H. Freeman and Company Plate tectonics explains the global pattern of volcanism. Key Figure 12.20 (page
More informationPhysical Geology, 15/e
Lecture Outlines Physical Geology, 15/e Plummer, Carlson & Hammersley Plate Tectonics: The Unifying Theory Physical Geology 15/e, Chapter 19 Plate Tectonics Plate Tectonics Earth s surface is composed
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 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 informationMost mafic magmas come from the upper mantle and lower crust. This handout will address five questions:
Geology 101 Origin of Magma From our discussions of the structure of the interior of the Earth, it is clear that the upper parts of the Earth (crust and mantle) are mostly solid because s-waves penetrate
More informationIntroduction to Subduction Zones
PAGEOPH, Vol. 128, Nos. 3/4 (1988) 0033~4553/88/040449-0551.50 + 0.20/0 9 1988 Birkh/iuser Verlag, Basel Introduction to Subduction Zones LARRY J. RUFF j and HIROO KANAMORI 2 Subduction zones present many
More information5. Convergent boundaries produce a relatively low number of earthquakes compared to other boundaries. a. True
1. Earth s crust is thinner than its mantle. ANSWER: True 2. The concept of isostacy states that high-density rock will stand higher than low-density rock, which explains the formation of subduction zones.
More informationAnswers: Internal Processes and Structures (Isostasy)
Answers: Internal Processes and Structures (Isostasy) 1. Analyse the adjustment of the crust to changes in loads associated with volcanism, mountain building, erosion, and glaciation by using the concept
More informationEarth s Continents and Seafloors. GEOL100 Physical Geology Ray Rector - Instructor
Earth s Continents and Seafloors GEOL100 Physical Geology Ray Rector - Instructor OCEAN BASINS and CONTINENTAL PLATFORMS Key Concepts I. Earth s rocky surface covered by of two types of crust Dense, thin,
More informationDirected Reading. Section: How Mountains Form MOUNTAIN RANGES AND SYSTEMS. Skills Worksheet
Skills Worksheet Directed Reading Section: How Mountains Form 1. How high is Mount Everest? a. about 1980 km above sea level b. more than 8 km below sea level c. more than 8 km above sea level d. more
More informationPlate tectonics, rock cycle
Dikes, Antarctica Rock Cycle Plate tectonics, rock cycle The Rock Cycle A rock is a naturally formed, consolidated material usually composed of grains of one or more minerals The rock cycle shows how one
More informationNorthern Tanzanian Earthquakes: Fault orientations, and depth distribution
Northern Tanzanian Earthquakes: Fault orientations, and depth distribution Stewart Rouse (NC A&T Physics) Penn State University SROP Mentors: Dr. Andy Nyblade & Dr. Rick Brazier July 27, 2005 1.0 Introduction
More informationThe Lithosphere and the Tectonic System. The Structure of the Earth. Temperature 3000º ºC. Mantle
The Lithosphere and the Tectonic System Objectives: Understand the structure of the planet Earth Review the geologic timescale as a point of reference for the history of the Earth Examine the major relief
More information