Geochemical composition of dykes along the Cameroon Line (CL): Petrogenesis and similarities with the Central Atlantic Magmatic Province

We report whole-rock geochemistry and Sr–Nd–Pb isotopic compositions of mafic dykes intruded in the Precambrian granito-gneissic basement complex, exposed at Nyos, Batibo, Dschang and Foumban on the Cameroon Line. The dykes are alkaline (Batibo), transitional (Foumban), and subalkaline (Nyos, Batibo and Dschang) with SiO₂ of 45–54 wt% and MgO of 2–9 wt%, similar to dykes reported in other areas of the Cameroon Line (CL) and the Central Atlantic Magmatic Province (CAMP). The abundances of rare earth elements (REE) and the Primitive Mantle normalised patterns for the Nyos, Batibo and Dschang dykes are similar to those of MORB, indicating that the dykes formed at shallower depths by a higher degree of partial melting relative to the Foumban dykes and the alkaline lavas of the CL. The transitional basaltic dykes with steeper REE patterns have their sources at deeper levels in the lithospheric mantle, possibly the garnet-spinel transition zone and were generated by a lower degree partial melting of the lithospheric and plume components. The Nyos and Batibo subalkaline dykes show similar isotopic compositions with a spectrum extending from depleted (DMM-like) to enriched (EM1-like) mantle, indicating the similarity in their source components. The Dschang dykes show distinct isotopic characteristics with relatively unradiogenic Nd-Pb isotope compositions compared to the Batibo and Nyos dykes. The Foumban transitional dykes with characteristic wide ranges in Sr-Nd-Pb isotopic compositions reveal varying contributions from enriched mantle components (EM1 and EM2) in addition to its plume signature similar to those of CL lavas. The Nyos and Batibo dykes alongside other dykes on the CL have low TiO₂ abundances (<2 wt%), negative PM-normalised Nb-anomalies, and moderately to strongly enriched REE patterns, and isotopic composition that overlaps with those of CAMP, suggesting a similar lithospheric origin.     

Finite-element simulations of long-period ground motions: Japanese subduction-zone earthquakes and the 1906 San Francisco earthquake

Large earthquakes at shallow depths commonly excite long-period ground motions in distant sedimentary basins, thereby inflicting damage upon large-scale structures. For example, the 2003 Tokachi-oki earthquake, Japan, damaged oil tanks in the Yufutsu Basin, located 250 km from the epicenter. Similar long-range effects were also observed during the 2004 earthquake off Kii Peninsula, Japan. In this study, we present the results of simulations of these earthquakes undertaken using the finite element method (FEM) with a voxel mesh. In addition, to examine whether the 1906 San Francisco earthquake excited long-period ground motions in the Los Angeles-area basins, we performed long-period ground motion simulations of most of the California region. The FEM simulations confirmed the importance of path effects for the development of long-period ground motions.     

Temporal changes in the chemistry of lake water within Yugama Crater,Kusatsu-Shirane Volcano,Japan: Implications for the evolution of the magmatic hydrothermal system

The chemical composition and temperature of lake water within Yugama Crater, Kusatsu-Shirane Volcano, Japan, have fluctuated markedly over the period 1966 to 2005. From 1966 to 1981, water temperature was relatively low, and Cl⁻ and SO₄²⁻ concentrations gradually decreased. These trends are explained by the growth of a sealing zone that formed between the region of hot rock surrounding the degassing magma and cool crust saturated with groundwater. The sealing zone acted to restrict the emission of magmatic volatiles.     

Volcanic tremor during eruptions: Temporal characteristics,scaling and constraints on conduit size and processes

We investigated characteristics of eruption tremor observed for 24 eruptions at 18 volcanoes based on published reports. In particular, we computed reduced displacements (D_R) to normalize the data and examined tremor time histories. We observed: (a) maximum D_R is approximately proportional to the square root of the cross sectional area of the vent, however, with lower than expected slope; (b) about one half of the cases show approximately exponential increases in D_R at the beginnings of eruptions, on a scale of minutes to hours; (c) one half of the cases show a sustained maximum level of tremor; (d) more than 90% of the cases show approximately exponential decay at the ends of eruptions, also on a scale of minutes to hours; and (e) exponential increases, if they occur, are commonly associated with the first large stage of eruptions. We estimate the radii of the vents using several methods and reconcile the topographic estimates, which are systematically too large, with those obtained from D_R itself and theoretical considerations. We compare scaling of tremor D_R with that for explosions and find that explosions have large absolute pressures and scale with vent radius squared, whereas tremor consists of pressure fluctuations that have lower amplitudes than the absolute pressure of explosions, and the scaling is different. We explore several methods to determine the appropriate scaling. This characteristic helps us to distinguish the type of eruptions: explosive (Vulcanian or Strombolian) eruptions versus sustained or continuous ash (e.g. Plinian) eruptions. Average eruption discharge, estimated from the total volume of tephra and the total duration of eruption tremor, is well correlated with peak discharge calculated from cross sectional area of the vent and velocity of volcanic ejecta. These results suggest similar scaling between different eruption types and the overall usefulness of monitoring tremor for evaluating volcanic activity.     

Misfit functionals in Laplace‐Fourier domain waveform inversion,with application to wide‐angle ocean bottom seismograph data

In seismic waveform inversion, non‐linearity and non‐uniqueness require appropriate strategies. We formulate four types of L2 normed misfit functionals for Laplace‐Fourier domain waveform inversion: i) subtraction of complex‐valued observed data from complex‐valued predicted data (the ‘conventional phase‐amplitude’ residual), ii) a ‘conventional phase‐only’ residual in which amplitude variations are normalized, iii) a ‘logarithmic phase‐amplitude’ residual and finally iv) a ‘logarithmic phase‐only’ residual in which the only imaginary part of the logarithmic residual is used. We evaluate these misfit functionals by using a wide‐angle field Ocean Bottom Seismograph (OBS) data set with a maximum offset of 55 km. The conventional phase‐amplitude approach is restricted in illumination and delineates only shallow velocity structures. In contrast, the other three misfit functionals retrieve detailed velocity structures with clear lithological boundaries down to the deeper part of the model. We also test the performance of additional phase‐amplitude inversions starting from the logarithmic phase‐only inversion result. The resulting velocity updates are prominent only in the high‐wavenumber components, sharpening the lithological boundaries. We argue that the discrepancies in the behaviours of the misfit functionals are primarily caused by the sensitivities of the model gradient to strong amplitude variations in the data. As the observed data amplitudes are dominated by the near‐offset traces, the conventional phase‐amplitude inversion primarily updates the shallow structures as a result. In contrast, the other three misfit functionals eliminate the strong dependence on amplitude variation naturally and enhance the depth of illumination. We further suggest that the phase‐only inversions are sufficient to obtain robust and reliable velocity structures and the amplitude information is of secondary importance in constraining subsurface velocity models.     

Seismic Wavefields in the Deep Seafloor Area from a Submarine Landslide Source

We use the finite difference method to simulate seismic wavefields at broadband land and seafloor stations for a given terrestrial landslide source, where the seafloor stations are located at water depths of 1,900–4,300 m. Our simulation results for the landslide source explain observations well at the seafloor stations for a frequency range of 0.05–0.1 Hz. Assuming the epicenter to be located in the vicinity of a large submarine slump, we also model wavefields at the stations for a submarine landslide source. We detect propagation of the Airy phase with an apparent velocity of 0.7 km/s in association with the seawater layer and an accretionary prism for the vertical component of waveforms at the seafloor stations. This later phase is not detected when the structural model does not consider seawater. For the model incorporating the seawater, the amplitude of the vertical component at seafloor stations can be up to four times that for the model that excludes seawater; we attribute this to the effects of the seawater layer on the wavefields. We also find that the amplification of the waveform depends not only on the presence of the seawater layer but also on the thickness of the accretionary prism, indicating low amplitudes at the land stations and at seafloor stations located near the trough but high amplitudes at other stations, particularly those located above the thick prism off the trough. Ignoring these characteristic structures in the oceanic area and simply calculating the wavefields using the same structural model used for land areas would result in erroneous estimates of the size of the submarine landslide and the mechanisms underlying its generation. Our results highlight the importance of adopting a structural model that incorporates the 3D accretionary prism and seawater layer into the simulation in order to precisely evaluate seismic wavefields in seafloor areas.     

Geochemical characteristics of Carboniferous greenstones in the Inner Zone of Southwest Japan

Abstract Greenstones, representing remnants of paleo-oceanic crust, occur in Permian and Jurassic accretionary complexes of the Inner Zone in the Southwestern Japan arc. The formation age of most of the greenstones is early Carboniferous, based on fossil ages for overlying limestones and Sm-Nd isotope ages of the greenstones themselves. The geochemistry of such greenstones is similar to those of present-day oceanic islands. Greenstones of the Permian accretionary complex (Akiyoshi belt) are alkalic and tholeiitic in composition. Some alkali basalts show peculiar features from an EM-1 mantle source, such as the Gough Island and Tristan da Chunha basalts in the South Atlantic. Greenstones of the Jurassic accretionary complex (Tamba belt) are also alkali and tholeiitic basalts with both basalt types in the northern part of the Tamba belt coming from strongly depleted characters similar to a mid-ocean ridge basalt source mantle. The variable geochemistry of the oceanic basalts is explained by hypothesis on existence of a Carboniferous mantle plume below the spreading ridge which divides the Farallon and Izanagi plates. The Akiyoshi belt seamounts and/or oceanic islands of the Farallon plate and Tamba belt seamounts and/or oceanic islands of the Izanagi plate formed simultaneously by the upwelling of the thermal plume. Some part of the Akiyoshi belt basalts originated locally from an EM-1 mantle source, while basalts from the northern parts of the Tamba belt have a normal-type mid-ocean ridge basalt (N-MORB) source component. Existence of an N-MORB signature is consistent with the presence of a spreading center in a Carboniferous 'Pacific Ocean' that caused separation of the Farallon and Izanagi plates. Disparity in accretion ages of the basaltic rocks in the Permian and Jurassic may have been caused by differences in the relative motion of the two plates.     

Dynamic rupture and stress change in a normal faulting earthquake in the subducting Cocos plate

A large nearly vertical, normal faulting earthquake ( M _w = 7.1) took place in 1997 in the Cocos plate, just beneath the ruptured fault zone of the great 1985 Michoacan thrust event ( M _w = 8.1). Dynamic rupture and resultant stress change during the 1997 earthquake have been investigated on the basis of near-source strong-motion records together with a 3-D dynamic model.
Dynamically consistent waveform inversion reveals a highly heterogeneous distribution of stress drop, including patch-like asperities and negative stress-drop zones. Zones of high stress drop are mainly confined to the deeper, southeastern section of the vertical fault, where the maximum dynamic stress drop reaches 280 bars (28 MPa). The dynamically generated source time function varies with location on the fault, and yields a short slip duration, which is caused by a short scalelength of stress-drop heterogeneities. The synthetic seismograms calculated from the dynamic model are generally consistent with the strong-motion velocity records in the frequency range lower than 0.5 Hz.
The pattern of stress-drop distribution appears, in some sense, to be consistent with that of coseismic changes in shear stress resulting from the 1985 thrust event. This consistency suggests that the stress transfer from the 1985 event to the subducting plate could be one of the possible mechanisms that increased the chance of the occurrence of the 1997 earthquake.     

Estimation of Oxygen Consumption Rate Using T-DO Diagram in the Benthic Layer of Ohmura Bay, Kyushu, Japan

We have observed the temporal variation of oxygen deficient water with short time scale (less than a few days) in the central area of Ohmura Bay, Kyushu, Japan, in summer, 1995 and 1996. The vertical profiles of temperature were similar to those of dissolved oxygen. We noticed a linear relation between temperature and dissolved oxygen in the bottom layer, and applied the T-DO relation to estimate the net oxygen consumption rate, rather than direct evaluation of the advection and diffusion. Oxygen consumption rate just above the bottom was estimated to be about 0.21 g O₂ m^–3day^–1 in July 1995, and about 0.28 g O₂ m^–3day^–1 in August 1996. The net oxygen consumption rate estimated for the bottom layer below the second thermocline was about 0.61 g O₂ m^–3day^–1 with variability from 0.55 to 0.66 g O₂ m^–3day^–1 during July 25 to 29, 1995. This is was about 0.64 g O₂ m^–3day^–1 with variability from 0.18 to 1.4 g O₂ m^–3day^–1 during August 22 to 30, 1996. The net oxygen consumption rates are about half of those measured with a closed system in the Seto Inland Sea.