Teleseismic array analysis of upper mantle compressional velocity structure

Large quantities of teleseismic short-period seismograms recorded at SCARLET provide travel time, apparent velocity and waveform data for study of upper mantle compressional velocity structure. Relative array analysis of arrival times from distant (30° < Δ < 95°) earthquakes at all azimuths constrains lateral velocity variations beneath southern California. We compare dT/dΔ back azimuth and averaged arrival time estimates from the entire network for 154 events to the same parameters derived from small subsets of SCARLET. Patterns of mislocation vectors for over 100 overlapping subarrays delimit the spatial extent of an east-west striking, high-velocity anomaly beneath the Transverse Ranges. Thin lens analysis of the averaged arrival time differences, called 'net delay' data, requires the mean depth of the corresponding lens to be more than 100 km. Our results are consistent with the PKP-delay times of Hadley and Kanamori (1977), who first proposed the high-velocity feature, but we place the anomalous material at substantially greater depths than their 40-100 km estimate.

Detailed analysis of travel time, ray parameter and waveform data from 29 events occurring in the distance range 9° to 40° reveals the upper mantle structure beneath an oceanic ridge to depths of over 900 km. More than 1400 digital seismograms from earthquakes in Mexico and Central America yield 1753 travel times and 58 dT/dΔ measurements as well as high-quality, stable waveforms for investigation of the deep structure of the Gulf of California. The result of a travel time inversion with the tau method (Bessonova et al., 1976) is adjusted to fit the p(Δ) data, then further refined by incorporation of relative amplitude information through synthetic seismogram modeling. The application of a modified wave field continuation method (Clayton and McMechan, 1981) to the data with the final model confirms that GCA is consistent with the entire data set and also provides an estimate of the data resolution in velocity-depth space. We discover that the upper mantle under this spreading center has anomalously slow velocities to depths of 350 km, and place new constraints on the shape of the 660 km discontinuity.

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for a comparative investigation of the upper mantle beneath the Cascade Ranges-Juan de Fuca region, an ocean-continent transit ion. These data consist of 853 seismograms (6° < Δ < 42°) which produce 1068 travel times and 40 ray parameter estimates. We use the spreading center model initially in synthetic seismogram modeling, and perturb GCA until the Cascade Ranges data are matched. Wave field continuation of both data sets with a common reference model confirms that real differences exist between the two suites of seismograms, implying lateral variation in the upper mantle. The ocean-continent transition model, CJF, features velocities from 200 and 350 km that are intermediate between GCA and T7 (Burdick and Helmberger, 1978), a model for the inland western United States. Models of continental shield regions (e.g., King and Calcagnile, 1976) have higher velocities in this depth range, but all four model types are similar below 400 km. This variation in rate of velocity increase with tectonic regime suggests an inverse relationship between velocity gradient and lithospheric age above 400 km depth.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Regional variations in upper mantle structure beneath North America

Several types of seismological data, including surface wave group and phase velocities, travel times from large explosions, and teleseismic travel time anomalies, have indicated that there are significant regional variations in the upper few hundred kilometers of the mantle beneath continental areas. Body wave travel times and amplitudes from large chemical and nuclear explosions are used in this study to delineate the details of these variations beneath North America.

As a preliminary step in this study, theoretical P wave travel times, apparent velocities, and amplitudes have been calculated for a number of proposed upper mantle models, those of Gutenberg, Jeffreys, Lehman, and Lukk and Nersesov. These quantities have been calculated for both P and S waves for model CIT11GB, which is derived from surface wave dispersion data. First arrival times for all the models except that of Lukk and Nersesov are in close agreement, but the travel time curves for later arrivals are both qualitatively and quantitatively very different. For model CIT11GB, there are two large, overlapping regions of triplication of the travel time curve, produced by regions of rapid velocity increase near depths of 400 and 600 km. Throughout the distance range from 10 to 40 degrees, the later arrivals produced by these discontinuities have larger amplitudes than the first arrivals. The amplitudes of body waves, in fact, are extremely sensitive to small variations in the velocity structure, and provide a powerful tool for studying structural details.

Most of eastern North America, including the Canadian Shield has a Pn velocity of about 8.1 km/sec, with a nearly abrupt increase in compressional velocity by ~ 0.3 km/sec near at a depth varying regionally between 60 and 90 km. Variations in the structure of this part of the mantle are significant even within the Canadian Shield. The low-velocity zone is a minor feature in eastern North America and is subject to pronounced regional variations. It is 30 to 50 km thick, and occurs somewhere in the depth range from 80 to 160 km. The velocity decrease is less than 0.2 km/sec.

Consideration of the absolute amplitudes indicates that the attenuation due to anelasticity is negligible for 2 hz waves in the upper 200 km along the southeastern and southwestern margins of the Canadian Shield. For compressional waves the average Q for this region is > 3000. The amplitudes also indicate that the velocity gradient is at least 2 x 10^-3 both above and below the low-velocity zone, implying that the temperature gradient is < 4.8°C/km if the regions are chemically homogeneous.

In western North America, the low-velocity zone is a pronounced feature, extending to the base of the crust and having minimum velocities of 7.7 to 7.8 km/sec. Beneath the Colorado Plateau and Southern Rocky Mountains provinces, there is a rapid velocity increase of about 0.3 km/sec, similar to that observed in eastern North America, but near a depth of 100 km.

Complicated travel time curves observed on profiles with stations in both eastern and western North America can be explained in detail by a model taking into account the lateral variations in the structure of the low-velocity zone. These variations involve primarily the velocity within the zone and the depth to the top of the zone; the depth to the bottom is, for both regions, between 140 and 160 km.

The depth to the transition zone near 400 km also varies regionally, by about 30-40 km. These differences imply variations of 250 °C in the temperature or 6 % in the iron content of the mantle, if the phase transformation of olivine to the spinel structure is assumed responsible. The structural variations at this depth are not correlated with those at shallower depths, and follow no obvious simple pattern.

The computer programs used in this study are described in the Appendices. The program TTINV (Appendix IV) fits spherically symmetric earth models to observed travel time data. The method, described in Appendix III, resembles conventional least-square fitting, using partial derivatives of the travel time with respect to the model parameters to perturb an initial model. The usual ill-conditioned nature of least-squares techniques is avoided by a technique which minimizes both the travel time residuals and the model perturbations.

Spherically symmetric earth models, however, have been found inadequate to explain most of the observed travel times in this study. TVT4, a computer program that performs ray theory calculations for a laterally inhomogeneous earth model, is described in Appendix II. Appendix I gives a derivation of seismic ray theory for an arbitrarily inhomogeneous earth model.

Part I. Approximate Hartree-Fock wavefunctions one-electron properties and electronic structure of the water molecule. Part II. Perturbation-variational calculation of the nuclear spin-spin isotope coupling constant in HD

Part I

Several approximate Hartree-Fock SCF wavefunctions for the ground electronic state of the water molecule have been obtained using an increasing number of multicenter s, p, and d Slater-type atomic orbitals as basis sets. The predicted charge distribution has been extensively tested at each stage by calculating the electric dipole moment, molecular quadrupole moment, diamagnetic shielding, Hellmann-Feynman forces, and electric field gradients at both the hydrogen and the oxygen nuclei. It was found that a carefully optimized minimal basis set suffices to describe the electronic charge distribution adequately except in the vicinity of the oxygen nucleus. Our calculations indicate, for example, that the correct prediction of the field gradient at this nucleus requires a more flexible linear combination of p-orbitals centered on this nucleus than that in the minimal basis set. Theoretical values for the molecular octopole moment components are also reported.

Part II

The perturbation-variational theory of R. M. Pitzer for nuclear spin-spin coupling constants is applied to the HD molecule. The zero-order molecular orbital is described in terms of a single 1s Slater-type basis function centered on each nucleus. The first-order molecular orbital is expressed in terms of these two functions plus one singular basis function each of the types e^-r/r and e^-r ln r centered on one of the nuclei. The new kinds of molecular integrals were evaluated to high accuracy using numerical and analytical means. The value of the HD spin-spin coupling constant calculated with this near-minimal set of basis functions is J_HD = +96.6 cps. This represents an improvement over the previous calculated value of +120 cps obtained without using the logarithmic basis function but is still considerably off in magnitude compared with the experimental measurement of J_HD = +43 0 ± 0.5 cps.

Vertical temperature structure of 250 m. water layer around Sri Lanka at the tail-end of the north-east monsoon

Data collected during cruises of the Hoyo Maru in Jan-Mar 1975 are analysed. Data is tabulated to show vertical temp profile, surface water temp and temp gradient. Each of these features is discussed. Thermoclines are shown to be established off the coast of India, their depth varying according to time of year. Upwelling off the Cochin coast is discussed. This occurs during Oct-Nov. Surface temp is considerably influenced by the north-east monsoon. The 'clockwise current' (or 'transparent current) characterised by high salinity, high transparency, and rich nutrient conchs, and which prevails in Jan-May brings oceanic water into the Bay of Bengal and sweeps along the Ceylon coast.

Molecular modeling on some human interleukins sharing gamma c of interleukin-2 receptor, and structure-function relationships

Five models for human interleukin-7 (HIL-7), HIL-9, HIL-13, HIL-15 and HIL-17 have been generated by SYBYL software package. The primary models were optimized using molecular dynamics and molecular mechanics methods. The final models were optimized using a steepest descent algorithm and a subsequent conjugate gradient method. The complexes with these interleukins and the common gamma chain of interleukin-2 receptor (IL-2R) were constructed and subjected to energy minimization. We found residues, such as Gln127 and Tyr103, of the common gamma chain of IL-2R are very important. Other residues, e.g. Lys70, Asn128 and Glu162, are also significant. Four hydrophobic grooves and two hydrophilic sites converge at the active site triad of the gamma chain. The binding sites of these interleukins interaction with the common gamma chain exist in the first helical and/or the fourth helical domains. Therefore, we conclude that these interleukins binds to the common gamma chain of IL-2R by the first and the fourth helix domain. Especially at the binding sites of some residues (lysine, arginine, asparagine, glutamic acid and aspartic acid), with a discontinuous region of the common gamma chain of IL-2R, termed the interleukins binding sites (103-210). The study of these sites can be important for the development of new drugs. (C) 2000 Elsevier Science B.V. All rights reserved.

MnSb/porous silicon hybrid structure prepared by physical vapor deposition

MnSb/porous silicon hybrid structure was prepared by physical vapor deposition technique. The structure and surface morphology of the MnSb films were analyzed by X-ray diffraction and scanning electron microscope, respectively. The magnetic hysteresis loops were obtained by an alternative gradient magnetometer. Based on the measurements, only MnSb phase was found and the surface morphology was rough and island-like. MnSb thin films show ferromagnetism at room temperature. (C) 2003 Elsevier B.V. All rights reserved.

Investigation on a pressure-gradient fiber laser hydrophone

In this paper, a pressure-gradient fiber laser hydrophone is demonstrated. Two brass diaphragms are installed at the end of a metal cylinder as sensing elements. A distributed feedback fiber laser, fixed at the center of the two diaphragms, is elongated or shortened due to the acoustic wave. There are two orifices at the middle of the cylinder. So this structure can work as a pressure-gradient microphone in the acoustic field. Furthermore, the hydrostatic pressure is self-compensated and an ultra-thin dimension is achieved. Theoretical analysis is given based on the electro-acoustic theory. Field trials are carried out to test the performance of the hydrophone. A sensitivity of 100 nm MPa-1 has been achieved. Due to the small dimensions, no directivity is found in the test.

band Structure calculations on orthorhombic bulk and nanocrystalline SrY2O4

The electronic structure of SrY2O4 is calculated by using a density functional method, and the exchange and correlation have been treated by using a the generalized gradient approximation (GGA) within the scheme due to Perdew, Burke, and Ernzerhof (PBE). SrY2O4 is predicted to be a direct-gap material because the top of the valence band and the bottom of the conduction band are along the same direction at G. The bond length and the bond covalency are also calculated by using a chemical bond method.

Magnetic structure and orbital ordering in tetragonal and monoclinic KCrF3 from first-principles calculations

KCrF3 has been systematically investigated by using the full-potential linearized augmented plane wave plus local orbital method within the generalized gradient approximation and the local spin density approximation plus the on-site Coulomb repulsion approach. The total energies for ferromagnetic and three different antiferromagnetic configurations are calculated in the high-temperature tetragonal and low-temperature monoclinic phases, respectively.

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

The variation of upper ocean structure and paleoproductivity in the Kuroshio source region during the last 200 kyr

A core from the source region of the Kuroshio warm current (east of the Luzon Island) was analyzed using several proxies in order to study the variability of the Western Pacific Warm Pool (WPWP) during the last two glacial-interglacial cycles. Primary productivity (PP) variations were deduced from variations in the coccolith flora. Primary productivity was higher during glacial periods (the end of Marine Isotope Stage [MIS] 3, some periods in MIS 2 and 6), and decreased during interglacial periods (MIS 7, MIS Se and probably MIS 5c-5d), with the lowest PP in MIS 5e. variations in the delta C-13 difference in benthic and bulk carbonate, thus in the vertical gradient of delta C-13 in dissolved inorganic carbon (Delta delta C-13(c). (wuellerstorfi-N. dutertrei) and Delta delta C-13(c.) (wuellerstorfi-coccolith)) Coincided With the PP Changes, showing that export productivity was low during interglacial periods (MIS 7, MIS 5e and Holocene) and high during glacial periods (MIS 6, probably MIS 5c-5d, late MIS 4 and late MIS 3). Comparison of foraminiferal carbonate dissolution indicators and PP changes reveals that nannofossil assemblage in core Ph05-5 is not sensitive to carbonate dissolution intensity. The depth of the thermocline (DOT) was estimated from planktonic forminiferal assemblages, and was relatively greater during interglacial periods (MIS 7, MIS 5e, probably MIS 5c and Holocene) than during glacials (middle MIS 6, probably MIS 5b and 5d, some periods in MIS 4, MIS 3 and MIS 2). Good coherence between the paleoproductivity records and the DOT suggests that the DOT changes could be the primary control factor in changes of paleoproductivity, and the glacial high productivity in the Kuroshio source region could be associated with a global increase of nutrient concentration in the intermediate waters that upwelled into the photic zone. The low CO2 values derived for intervals of high productivity and a relatively shallow DOT suggest that the changes in biological productivity and DOT in the equatorial Pacific could have modified atmospheric CO2 concentrations. High Sea Surface Temperatures (SSTs) during the warm MIS 5e in combination with intensified monsoonal rain fall could have resulted in a more intense stratification of the upper waters, resulting in low nutrient supply to the surface waters and a resulting decrease in productivity. (C) 2010 Elsevier B.V. All rights reserved.

Reconstruction of vertical thermal structure from several subsurface temperatures in the China Seas and adjacent waters

Empirical Orthogonal Function (EOF) analysis is used in this study to generate main eigenvector fields of historical temperature for the China Seas (here referring to Chinese marine territories) and adjacent waters from 1930 to 2002 (510 143 profiles). A good temperature profile is reconstructed based on several subsurface in situ temperature observations and the thermocline was estimated using the model. The results show that: 1) For the study area, the former four principal components can explain 95% of the overall variance, and the vertical distribution of temperature is most stable using the in situ temperature observations near the surface. 2) The model verifications based on the observed CTD data from the East China Sea (ECS), South China Sea (SCS) and the areas around Taiwan Island show that the reconstructed profiles have high correlation with the observed ones with the confidence level > 95%, especially to describe the characteristics of the thermocline well. The average errors between the reconstructed and observed profiles in these three areas are 0.69A degrees C, 0.52A degrees C and 1.18A degrees C respectively. It also shows the model RMS error is less than or close to the climatological error. The statistical model can be used to well estimate the temperature profile vertical structure. 3) Comparing the thermocline characteristics between the reconstructed and observed profiles, the results in the ECS show that the average absolute errors are 1.5m, 1.4 m and 0.17A degrees C/m, and the average relative errors are 24.7%, 8.9% and 22.6% for the upper, lower thermocline boundaries and the gradient, respectively. Although the relative errors are obvious, the absolute error is small. In the SCS, the average absolute errors are 4.1 m, 27.7 m and 0.007A degrees C/m, and the average relative errors are 16.1%, 16.8% and 9.5% for the upper, lower thermocline boundaries and the gradient, respectively. The average relative errors are all < 20%. Although the average absolute error of the lower thermocline boundary is considerable, but contrast to the spatial scale of average depth of the lower thermocline boundary (165 m), the average relative error is small (16.8%). Therefore the model can be used to well estimate the thermocline.

Phytoplankton biomass size structure and its regulation in the Southern Yellow Sea (China): Seasonal variability

Phytoplankton size structure plays a significant role in controlling the carbon flux of marine pelagic ecosystems. The mesoscale distribution and seasonal variation of total and size-fractionated phytoplankton biomass in surface waters. as measured by chlorophyll a (Chl a), was studied in the Southern Yellow Sea using data from four cruises during 2006-2007. The distribution of Chl a showed a high degree of spatial and temporal variation in the study area. Chl a concentrations were relatively high in the summer and autumn, with a mean of 142 and 1.27 mg m(-3), respectively. Conversely, in the winter and spring. the average Chl a levels were only 098 and 0.99 mg m(-3) Total Chl a showed a clear decreasing gradient from coastal areas to the open sea in the summer, autumn and winter cruises. Patches of high Chl a were observed in the central part of the Southern Yellow Sea in the spring due to the onset of the phytoplankton bloom. The eutrophic coastal waters contributed at least 68% of the total phytoplankton biomass in the surface layer. Picophytoplankton showed a consistent and absolute dominance in the central region of the Southern Yellow Sea (>40%) in all of the cruises, while the proportion of microphytoplankton was the highest in coastal waters The relative proportions of pico- and nanophytoplankton decreased with total biomass, whereas the proportion of the micro-fraction increased with total biomass. Relationships between phytoplankton biomass and environmental factors were also analysed. The results showed that the onset of the spring bloom was highly dependent on water column stability. Phytoplankton growth was limited by nutrient availability in the summer due to the strong thermocline. The combined effects of P-limitation and vertical mixing in the autumn restrained the further increase of phytoplankton biomass in the Surface layer. The low phytoplankton biomass in winter was caused by vertical dispersion due to intense mixing. Compared with the availability of nutrients. temperature did not seem to cause direct effects on phytoplankton biomass and its size structure. Although interactions of many different environmental factors affected phytoplankton distributions. hydrodynamic conditions seemed to be the dominant factor. Phytoplankton size structure was determined mainly by the size-differential capacity in acquiring resource. Short time scale events, such as the spring bloom and the extension of Yangtze River plume, can have substantial influences, both on the total Chl a concentration and on the size structure of the phytoplankton. (C) 2009 Elsevier Ltd. All rights reserved.