Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan

The fine-scale seismic structure of the central Mexico, southern Peru, and southwest Japan subduction zones is studied using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. The morphology of the transition from flat to normal subduction is explored in central Mexico and southern Peru, while in southwest Japan the spatial coincidence of a thin ultra-slow velocity layer (USL) atop the flat slab with locations of slow slip events (SSEs) is explored. This USL is also observed in central Mexico and southern Peru, where its lateral extent is used as one constraint on the nature of the flat-to-normal transitions.

In western central Mexico, I find an edge to this USL which is coincident with the western boundary of the projected Orozco Fracture Zone (OFZ) region. Forward modeling of the 2D structure of the subducted Cocos plate using a finite-difference algorithm provides constraints on the velocity and geometry of the slab’s seismic structure in this region and confirms the location of the USL edge. I propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, by a process analogous to that which occurred when the Rivera plate separated from the proto-Cocos plate 10 Ma.

In eastern central Mexico, observations of a sharp transition in slab dip near the abrupt end of the Trans Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of the USL is found to be coincident with these features and with the western boundary of a zone of decreased seismicity, indicating a change in structure which I interpret as evidence of a possible tear. Analysis of intraslab seismicity patterns and focal mechanism orientations and faulting types provides further support for a possible tear in the South Cocos slab. This potential tear, together with the tear along the projection of the OFZ to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments.

In southern Peru, observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest a smooth contortion of the slab. Concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry. The lateral extent of the USL atop the horizontal Nazca slab is found to be coincident with the margin of the projected linear continuation of the subducting Nazca Ridge, implying a causal relationship, but not a slab tear. Waveform modeling of the 2D structure in southern Peru provides constraints on the velocity and geometry of the slab’s seismic structure and confirms the absence of any tears in the slab.

In southwest Japan, I estimate the location of a possible USL along the Philippine Sea slab surface and find this region of low velocity to be coincident with locations of SSEs that have occurred in this region. I interpret the source of the possible USL in this region as fluids dehydrated from the subducting plate, forming a high pore-fluid pressure layer, which would be expected to decrease the coupling on the plate interface and promote SSEs.

Seismic structure above and below the core-mantle boundary

Seismic structure above and below the core-mantle boundary (CMB) has been studied through use of travel time and waveform analyses of several different seismic wave groups. Anomalous systematic trends in observables document mantle heterogeneity on both large and small scales. Analog and digital data has been utilized, and in many cases the analog data has been optically scanned and digitized prior to analysis.

Differential travel times of S - SKS are shown to be an excellent diagnostic of anomalous lower mantle shear velocity (V s) structure. Wavepath geometries beneath the central Pacific exhibit large S- SKS travel time residuals (up to 10 sec), and are consistent with a large scale 0(1000 km) slower than average V_s region (≥3%). S - SKS times for paths traversing this region exhibit smaller scale patterns and trends 0(100 km) indicating V_s perturbations on many scale lengths. These times are compared to predictions of three tomographically derived aspherical models: MDLSH of Tanimoto [1990], model SH12_WM13 of Suet al. [1992], and model SH.10c.17 of Masters et al. [1992]. Qualitative agreement between the tomographic model predictions and observations is encouraging, varying from fair to good. However, inconsistencies are present and suggest anomalies in the lower mantle of scale length smaller than the present 2000+ km scale resolution of tomographic models. 2-D wave propagation experiments show the importance of inhomogeneous raypaths when considering lateral heterogeneities in the lowermost mantle.

A dataset of waveforms and differential travel times of S, ScS, and the arrival from the D" layer, Scd, provides evidence for a laterally varying V_s velocity discontinuity at the base of the mantle. Two different localized D" regions beneath the central Pacific have been investigated. Predictions from a model having a V_s discontinuity 180 km above the CMB agree well with observations for an eastern mid-Pacific CMB region. This thickness differs from V_s discontinuity thicknesses found in other regions, such as a localized region beneath the western Pacific, which average near 280 km. The "sharpness" of the V_s jump at the top of D", i.e., the depth range over which the V_s increase occurs, is not resolved by our data, and our data can in fact may be modeled equally well by a lower mantle with the increase in V_s at the top of D" occurring over a 100 krn depth range. It is difficult at present to correlate D" thicknesses from this study to overall lower mantle heterogeneity, due to uncertainties in the 3-D models, as well as poor coverage in maps of D" discontinuity thicknesses.

P-wave velocity structure (V_p) at the base of the mantle is explored using the seismic phases SKS and SPdKS. SPdKS is formed when SKS waves at distances around 107° are incident upon the CMB with a slowness that allows for coupling with diffracted P-waves at the base of the mantle. The P-wave diffraction occurs at both the SKS entrance and exit locations of the outer core. SP_dKS arrives slightly later in time than SKS, having a wave path through the mantle and core very close to SKS. The difference time between SKS and SP_dKS strongly depends on V_p at the base of the mantle near SK Score entrance and exit points. Observations from deep focus Fiji-Tonga events recorded by North American stations, and South American events recorded by European and Eurasian stations exhibit anomalously large SP_dKS - SKS difference times. SKS and the later arriving SP_dKS phase are separated by several seconds more than predictions made by 1-D reference models, such as the global average PREM [Dziewonski and Anderson, 1981] model. Models having a pronounced low-velocity zone (5%) in V_p in the bottom 50-100 km of the mantle predict the size of the observed SP_dK S-SKS anomalies. Raypath perturbations from lower mantle V_s structure may also be contributing to the observed anomalies.

Outer core structure is investigated using the family of SmKS (m=2,3,4) seismic waves. SmKS are waves that travel as S-waves in the mantle, P-waves in the core, and reflect (m-1) times on the underside of the CMB, and are well-suited for constraining outermost core V_p structure. This is due to closeness of the mantle paths and also the shallow depth range these waves travel in the outermost core. S3KS - S2KS and S4KS - S3KS differential travel times were measured using the cross-correlation method and compared to those from reflectivity synthetics created from core models of past studies. High quality recordings from a deep focus Java Sea event which sample the outer core beneath the northern Pacific, the Arctic, and northwestern North America (spanning 1/8th of the core's surface area), have SmKS wavepaths that traverse regions where lower mantle heterogeneity is pre- dieted small, and are well-modeled by the PREM core model, with possibly a small V_p decrease (1.5%) in the outermost 50 km of the core. Such a reduction implies chemical stratification in this 50 km zone, though this model feature is not uniquely resolved. Data having wave paths through areas of known D" heterogeneity (±2% and greater), such as the source-side of SmKS lower mantle paths from Fiji-Tonga to Eurasia and Africa, exhibit systematic SmKS differential time anomalies of up to several seconds. 2-D wave propagation experiments demonstrate how large scale lower mantle velocity perturbations can explain long wavelength behavior of such anomalous SmKS times. When improperly accounted for, lower mantle heterogeneity maps directly into core structure. Raypaths departing from homogeneity play an important role in producing SmKS anomalies. The existence of outermost core heterogeneity is difficult to resolve at present due to uncertainties in global lower mantle structure. Resolving a one-dimensional chemically stratified outermost core also remains difficult due to the same uncertainties. Restricting study to higher multiples of SmKS (m=2,3,4) can help reduce the affect of mantle heterogeneity due to the closeness of the mantle legs of the wavepaths. SmKS waves are ideal in providing additional information on the details of lower mantle heterogeneity.

^(18)O/^(16)O and D/H studies on the interactions between heated meteoric ground waters and igneous intrusions: Western San Juan Mountains, Colorado and the Isle of Skye, Scotland

Detailed oxygen, hydrogen and carbon isotope studies have been carried out on igneous and metamorphic rocks of the Stony Mountain complex, Colorado, and the Isle of Skye, Scotland, in order to better understand the problems of hydrothermal meteoric water-rock interaction.

The Tertiary Stony Mountain stock (~1.3 km in diameter), is composed of an outer diorite, a main mass of biotite gabbro, and an inner diorite. The entire complex and most of the surrounding country rocks have experienced various degrees of ¹⁸O depletion (up to 10 per mil) due to interaction with heated meteoric waters. The inner diorite apparently formed from a low-¹⁸O magma with δ¹⁸O ≃ +2.5, but most of the isotopic effects are a result of exchange between H₂O and solidified igneous rocks. The low-¹⁸O inner diorite magma was probably produced by massive assimilation and/or melting of hydrothermally altered country rocks. The δ¹⁸O values of the rocks generally increase with increasing grain size, except that quartz typically has δ¹⁸O = +6 to +8, and is more resistant to hydrothermal exchange than any other mineral studied. Based on atom % oxygen, the outer diorites, gabbros, and volcanic rocks exhibit integrated water/rock ratios of 0.3 ± 0.2, 0.15 ± 0.1, and 0.2 ± 0.1, respectively. Locally, water/rock ratios attain values greater than 1.0. Hydrogen isotopic analyses of sericites, chlorites, biotites, and amphiboles range from -117 to -150. δD in biotites varies inversely with Fe/Fe+Mg, as predicted by Suzuoki and Epstein (1974), and positively with elevation, over a range of 600 m. The calculated δD of the mid-to-late-Tertiary meteoric waters is about -100. Carbonate δ¹³C values average -5.5 (PDB), within the generally accepted range for deep-seated carbon.

Almost all the rocks within 4 km of the central Tertiary intrusive complex of Skye are depleted in ¹⁸O. Whole-rock δ¹⁸O values of basalts (-7. 1 to +8.4), Mesozoic shales (-0.6 to + 12.4), and Precambrian sandstones (-6.2 to + 10.8) systematically decrease inward towards the center of the complex. The Cuillin gabbro may have formed from a ¹⁸O-depleted magma (depleted by about 2 per mil); δ¹⁸O of plagioclase (-7.1 to + 2.5) and pyroxene (-0.5 to + 3.2) decrease outward toward the margins of the pluton. The Red Hills epigranite plutons have δ¹⁸O quartz (-2.7 to + 7.6) and feldspar (-6.7 to + 6.0) that suggest about 3/4 of the exchange took place at subsolidus temperatures; profound disequilibrium quartz-feldspar fractionations (up to 12) are characteristic. The early epigranites were intruded as low-¹⁸O melts (depletions of up to 3 per mil) with δ¹⁸O of the primary, igneous quartz decreasing progressively with time. The Southern Porphyritic Epigranite was apparently intruded as a low-¹⁸O magma with δ¹⁸O ≃ -2.6. A good correlation exists between grain size and δ¹⁸O for the unique, high-¹⁸O Beinn an Dubhaich granite which intrudes limestone having a δ¹⁸O range of +0.5 to +20.8, and δ¹³C of -4.9 to -1.0. The δD values of sericites (-104 to -107), and amphiboles, chlorites, and biotites (-105 to -128) from the igneous rocks , indicate that Eocene surface waters at Skye had δD ≃ -90. The average water/rock ratio for the Skye hydrothermal system is approximately one; at least 2000 km³ of heated meteoric waters were cycled through these rocks.

Thus these detailed isotopic studies of two widely separated areas indicate that (1) ¹⁸O-depleted magmas are commonly produced in volcanic terranes invaded by epizonal intrusions; (2) most of the ¹⁸O-depletion in such areas are a result of subsolidus exchange (particularly of feldspars); however correlation of δ¹⁸O with grain size is generally preserved only for systems that have undergone relatively minor meteoric hydrothermal exchange; (3) feldspar and calcite are the minerals mos t susceptible to oxygen isotopic exchange, whereas quartz is very resistant to oxygen isotope exchange; biotite, magnetite, and pyroxene have intermediate susceptibilities; and (4) basaltic country rocks are much more permeable to the hydrothermal convective system than shale, sandstone, or the crystalline basement complex.

Studies toward the total synthesis of ritterazine B : the investigation of alkynylation reactions for use in the synthesis of the western fragment

The ritterazine and cephalostatin natural products have biological activities and structures that are interesting to synthetic organic chemists. These products have been found to exhibit significant cytotoxicity against P388 murine leukemia cells, and therefore have the potential to be used as anticancer drugs. The ritterazines and cephalostatins are steroidal dimers joined by a central pyrazine ring. Given that the steroid halves are unsymmetrical and highly oxygenated, there are several challenges in synthesizing these compounds in an organic laboratory.

Ritterazine B is the most potent derivative in the ritterazine family. Its biological activity is comparable to drugs that are being used to treat cancer today. For this reason, and the fact that there are no reported syntheses of ritterazine B to date, our lab set out to synthesize this natural product.

Herein, efforts toward the synthesis of the western fragment of ritterazine B are described. Two different routes are explored to access a common intermediate. An alkyne conjugate addition reaction was initially investigated due to the success of this key reaction in the synthesis of the eastern fragment. However, it has been found that a propargylation reaction has greater reactivity and yields, and has the potential to reduce the step count of the synthesis of the western fragment of ritterazine B.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.