Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Isotope ratios of basalts from ODP Leg 115 holes (Table 1)

Strontium, neodymium, and lead isotope ratios are reported for 13 Leg 115 basalts as well as 3 basalts from Texaco drill hole SM-1 on the Mascarene Plateau. The 87Sr/86Sr ratios and eNd range from 0.70330 to 0.70439 and 5.5 to 7.4, respectively, although 87Sr/86Sr ratios higher than 0.70383 are found only in SM-1 basalts. The high 87Sr/86Sr values are thought to reflect seawater Sr in secondary phases, although all samples were strongly leached in HC1 before analysis. 206Pb/204Pb ratios range from 18.53 to 18.80, and sho high 207Pb/204Pb and 208Pb/204Pb ratios relative to 206Pb/204Pb ratios, typical of Indian Ocean mid-ocean ridge (MORB) and oceanic-island basalts (OIB). Isotopic compositions of Leg 115 basalts generally fall between fields for MORB and Reunion Island basalts, consistent with the conclusion drawn from geochronological studies that Deccan flood basalt volcanism, the Chagos-Laccadive Ridge, and the Mascarene Plateau are all products of the Reunion mantle plume. Isotopic compositions of magmas produced by this plume have varied systematically with time in the direction of less "depleted," less MORB-like isotopic signatures. This compositional change has been accompanied by a decrease in eruption rate. We interpret Deccan volcanism as the voluminous beginning of the plume. Reduced entrainment of asthenosphere following melting of the plume head resulted in less MORB-like isotope ratios in magmas and a decrease in eruptive activity with time.

40Ar-39Ar plateau and isochron ages from the Ninetyeast and Broken Ridges (Table 1)

Concordant plateau and isochron ages are obtained from 40Ar-39Ar incremental heating experiments on volcanic rocks recovered by drilling at three Leg 121 sites along the Ninetyeast Ridge and two dredge locations on the southern scarp of the Broken Ridge, eastern Indian Ocean. The new data confirm a northerly increase in the age of volcanism along the Ninetyeast Ridge, from 38 to 82 Ma; this lineament links current hotspot volcanism near the Kerguelen islands with the Rajmahal flood basalt eruptions at M0 time (117 ± 1 Ma). The Broken Ridge was formed over the same hotspot at 88-89 Ma, but later experienced rift-related volcanism in Paleocene time (63 Ma). The geometry and distribution of ages along these prominent volcanic ridges and the Mascarene-Chagos-Laccadive-Maldive ridge system in the western Indian Ocean are most compatible with plate motions over fixed hotspots near Kerguelen and Reunion islands, respectively.

Primary mineralogy of ODP Leg 115 basalts

Basement rocks were recovered at four sites on Leg 115 along the Reunion hotspot track in the western Indian Ocean. Plate tectonic reconstructions indicate that the drilled structures formed in three different volcanic environments. Sites 706 and 713 from the eastern side of the Saya de Malha Bank and the northern end of the Chagos Bank, respectively, are on a large volcanic platform analogous to Iceland on the Mid-Atlantic Ridge. Lavas at Site 707 on the northwestern side of the Saya de Malha Bank erupted during the early stages of rifting of the Seychelles from India. Basalts from Site 715 were erupted onto an isolated oceanic island that was distant from ocean ridges and continents much as Reunion Island is today. Many of the rocks were examined in thin section and found to be primarily augite-plagioclase basalts with minor olivine and rare opaque oxides. Site 715 is unusual in that it contains a variety of basalts including olivine-rich and aphyric Fe-Ti basalts. At each of the four sites the rocks were grouped into chemical types (units) on the basis of ship- board bulk-rock analyses and at least one thin section from each chemical unit was analyzed by electron microprobe. The plagioclase and augite chemistry reflects the bulk-rock chemistry and, in general, these minerals were in equilibrium with their host magmas at the time the basalts were quenched. Olivine was rarely preserved, but where it is still present it also appears to have crystallized in equilibrium with the host magma. At three of the drill sites plagioclase phenocrysts or megacrysts that crystallized from a primitive magma are also present. The one site (715) that does not contain these primitive plagioclase phenocrysts is also the site that appears to have been influenced the least by ocean- ridge or Deccan-type magmas. Site 715, furthermore, has a mineralogy that is dominated by olivine as compared with the plagioclase-rich lavas of the other sites.

(Table 1) 40Ar/39Ar plateau and isochron ages for basalts from the Réunion hotspot track

Concordant plateau and isochron ages were calculated from 40Ar/39Ar incremental heating experiments on volcanic rocks recovered by drilling at four Leg 115 sites and two industry wells along the volcanic lineament connecting Reunion Island to the Deccan flood basalts, western Indian Ocean. The new ages provide unequivocal evidence that volcanic activity migrated southward along this sequence of linear ridges. The geometry and age distribution of volcanism are most compatible with origin above a stationary hotspot centered beneath Reunion. The hotspot became active with rapid eruption of the Deccan flood basalts, western India, and subsequent volcanic products record the northward motion of the Indian and African plates over the hotspot through Tertiary time. The radiometric ages are in general accord with basal biostratigraphic age estimates, although some adjustments in current magnetobiostratigraphic time scales may be required.

Eocene-Oligocene diatoms in the western Indian Ocean

The occurrence of diatom species in the Eocene-Oligocene sections of Ocean Drilling Program (ODP) Leg 115 sites and Deep Sea Drilling Project (DSDP) Sites 219 and 236 in the low-latitude Indian Ocean are investigated. Diatoms are generally rare and poorly preserved in the Paleogene sequences we studied. The best-preserved assemblages are found close to ash layers in early Oligocene sediments. The low-latitude diatom zonation established for the Atlantic region by Fenner in 1984 is fully applicable to the Paleogene sequences of the western Indian Ocean. Correlation of the diatom zones to the calcareous nannofossil stratigraphy of the sites places the Coscinodiscus excavatus Zone of Fenner within calcareous nannofossil Subzone CP16b. For the Mascarene Plateau and the Chagos Ridge, the times when the sites studied, together with the areas upslope from them, subsided to below the euphotic zone are deduced from changes in the relative abundance between the group of benthic, shallow-water species and Grammatophora spp. vs. the group of fully planktonic diatom species. The Eocene section of Site 707, on the Mascarene Plateau, is characterized by the occurrence of benthic diatoms (approximately 10% of the diatom assemblage). These allochthonous diatoms must have originated from shallow-water environments around volcanic islands that existed upslope from ODP Site 707 in Eocene times. In Oligocene and younger sediments of Sites 707 and 706, occurrences of benthic diatoms are rare and sporadic and interpreted as reworked from older sediments. This indicates that the area upslope from these two Mascarene Plateau sites had subsided below the euphotic zone by the early Oligocene. Only Grammatophora spp., for which a neritic but not benthic habitat is assumed, continues to be abundant throughout the Oligocene sequences. The area of the Madingley Rise sites (Sites 709-710) and nearby shallower areas subsided below the euphotic zone already in middle Eocene times, as benthic diatoms are almost absent from these Eocene sections. Only sites located on abyssal plains, and which intermittently received turbidite sediments (e.g., Sites 708 and 711), contain occasionally single, benthic diatoms of Oligocene age. The occurrence of the freshwater diatom Aulacosira granulata in a few samples of late early Oligocene and late Oligocene age at Sites 707, 709, and 714 is interpreted as windblown. Their presence indicates at least seasonally arid conditions for these periods in the source areas of eastern Africa and India. Three new species and two new combinations are defined: Chaetoceros asymmetricus Fenner sp. nov.; Hemiaulus gracilis Fenner, sp. nov.; Kozloviella meniscosa Fenner, sp. nov.; Cestodiscus demergitus (Fenner) Fenner comb, nov.; and Rocella princeps (Jouse) Fenner comb. nov.

Annotated record of the detailed examination of Mn deposits from ANTIPODE Expedition stations

Deep-sea sediment cores from Scripps Institution of Oceanography's ANTIPODE Expedition were described to identify visually distinct units based on color, texture, or other feature, sedimentary structures, lithology and abundance of component grains, and paleontology. Sixty-eight cores were examined, of which 34 are large diameter piston cores. Photographs and graphic lithology legs are included as PLATES 1-48. ANTIPODE Expedition recovered cores from: the Monterey-Ascension Fan, the Northeast Pacific, the Aleutians, the Northwest Pacific, the Philippine Sea, Indonesia, the Tonga Ridge, the Seychelles, Chagos Archipelago, the Mid-Indian Ridge, the Bay of Bengal, near Sumatra, and near the Cocos Island in the Indian Ocean. The purpose of this report is to present sufficient basic data on ANTIPODE cores for invesiigators to choose samples for their own research.