Highly alkaline,high-temperature hydrothermal fluids in the early Archean ocean

Based on the petrology of hydrothermally altered Archean basaltic greenstones, thermodynamic calculations of phase equilibria were conducted to estimate the composition of a high-temperature (∼350 °C) hydrothermal fluid in an Archean subseafloor basalt-hosted hydrothermal system. The results indicate that the hydrothermal fluid was highly alkaline attributed to the presence of calcite in the alteration minerals under a high-CO₂ condition, and predict a generation of SiO₂-rich, Fe-poor hydrothermal fluids in the Archean subseafloor hydrothermal system. The chemically reactive mixing zones between alkaline hydrothermal fluids and slightly acidic-neutral seawater are characterized by inverse pH and chemical polarity to modern hydrothermal systems, leading to extensive precipitation of silica and iron oxyhydroxides on/under the seafloor. Such processes can be responsible for the abiotic formation of voluminous chert and subseafloor silica dike, the mechanism of silicification, and the pH-controlled generation of banded iron formation that has been arising mainly from the redox chemistry in the Archean ocean. Such high-temperature alkaline fluids could have had a significant role not only in the early ocean geochemical processes but also in the early evolution of life.     

Influence of interaction between small asperities on various types of slow earthquakes in a 3-D simulation for a subduction plate boundary

Recently, the occurrence of slow earthquakes such as low-frequency earthquakes and very low-frequency earthquakes have been recognized at depths of about 30 km in southwest Japan and Cascadia. These slow earthquakes occur sometimes in isolation and sometimes break into chain-reaction, producing tremor that migrates at a speed of about 5–15 km/day and suggesting a strong interaction among nearby small asperities. In this study, we formulate a 3-D subduction plate boundary model with two types of small asperities chained along the trench at the depth of 30 km. Our simulation succeeds in representing various types of slow earthquakes including low-frequency earthquakes and rapid slip velocity in the same asperity, and indicates that interaction between asperities may cause the very low-frequency earthquakes. Our simulation also shows chain reaction along trench with propagation speed that can be made consistent with observations by adjusting model parameters, which suggests that the interactions also explain the observed migration of slow earthquakes.     

High precision slip distribution of the 1944 Tonankai earthquake inferred from tsunami waveforms: Possible slip on a splay fault

Inversion of tsunami waveforms is a well-established technique for estimating the slip distributions of subduction zone earthquakes, with some of the most detailed results having been obtained for earthquakes in the Nankai Trough, SW Japan. The present study, although it uses a method and tsunami waveform data set almost identical to previous study, aims to improve on previous work by using a more precise specification of initial conditions for the calculation of tsunami Green's functions. Specifically, we incorporated four improvements in the present study: (1) we used a realistic plate model based only on seismic survey results, and assumed it to be the fault plane of the 1944 Tonankai earthquake; (2) the smallest subfaults consistent with the long wavelength approximation were used in the tsunami inversion analysis; (3) we included the effect of horizontal displacement of the ocean bottom on tsunami generation; and (4) we performed a checkerboard resolution test. As obtained in previous studies, a zone of high slip (> 2.0 m) was resolved off the Shima Peninsula. However, the more precise calculation of tsunami Green's functions has revealed additional detail that was not evident in previous studies, which we demonstrate is resolvable and correlates with the position of known faults in the accretionary prism. While there was little or no slip near the trench axis in the eastern part of the rupture zone, there was up to 1.5 m of slip resolved within 30 km of the trough axis in the western part, along the coast of the Kii Peninsula. This troughward slip zone coincides with the position of a large splay fault mapped in multichannel reflection surveys. Furthermore, it is also clear that the upper edge of the Enshu fault off Shima and Atsumi peninsulas is consistent with the up-dip limit of slip in the eastern part of our model. We tested the possibility that slip occurred on the former splay fault instead of on the plate interface during the 1944 Tonankai earthquake, and find that slip on this splay fault is also consistent with the data, although we cannot distinguish whether slip was dominant on the splay fault or on the plate interface. We further suggest that the position of the Enshu fault may be determined by the subduction of topographic highs, and that such faults may have an important influence on the up-dip rupture limit of the 1944 Tonankai and, potentially, other subduction zone earthquakes.     

Software Dedicated for the Curation of Geochemical Data Sets in Analytical Laboratories

Software designed for analytical laboratories to guarantee traceability and accessibility of rocks with their geochemical properties has been developed. The software documents the sample origin, current sample location and the location of any sample subsets (e.g., thin sections, solutions, etc.), and archives all associated geochemical data sets. The software can be installed on a personal computer so is available for use in any laboratory and allows curation before and after publication. The software will be of use in integrating and sharing geological reference materials within and among institutes. In this article, the system design and implementation are detailed. All source codes for the software are available at http://dream.misasa.okayama-u.ac.jp/ .     

Boron cycling by subducted lithosphere; insights from diamondiferous tourmaline from the Kokchetav ultrahigh-pressure metamorphic belt

Subduction of lithosphere, involving surficial materials, into the deep mantle is fundamental to the chemical evolution of the Earth. However, the chemical evolution of the lithosphere during subduction to depth remains equivocal. In order to identify materials subjected to geological processes near the surface and at depths in subduction zones, we examined B and Li isotopes behavior in a unique diamondiferous, K-rich tourmaline (K-tourmaline) from the Kokchetav ultrahigh-pressure metamorphic belt. The K-tourmaline, which includes microdiamonds in its core, is enriched in ¹¹B relative to ¹⁰B (δ¹¹B = −1.2 to +7.7) and ⁷Li relative to ⁶Li (δ⁷Li = −1.1 to +3.1). It is suggested that the K-tourmaline crystallized at high-pressure in the diamond stability field from a silicate melt generated at high-pressure and temperature conditions of the Kokchetav peak metamorphism. The heavy isotope signature of this K-tourmaline differs from that of ordinary Na-tourmalines in crustal rocks, enriched in the light B isotope (δ¹¹B = −16.6 to −2.3), which experienced isotope fractionation through metamorphic dehydration reactions. A possible source of the heavy B-isotope signature is serpentine in the subducted lithospheric mantle. Serpentinization of the lithospheric mantle, with enrichment of heavy B-isotope, can be produced by normal faulting at trench-outer rise or trench slope regions, followed by penetration of seawater into the lithospheric mantle. Serpentine breakdown in the lithospheric mantle subducted in subarc regions likely provided fluids with the heavy B-isotope signature, which was acquired during the serpentinization prior to subduction. The fluids could ascend and cause partial melting of the overlying crustal layer, and the resultant silicate melt could inherit the heavy B-isotope signature. The subducting lithospheric mantle is a key repository for modeling the flux of fluids and associated elements acquired at a near the surface into the deep mantle.     

Deformable backstop as seaward end of coseismic slip in the Nankai Trough seismogenic zone

Wide-angle seismic surveys performed in the last decade have clarified the 3-D crustal structure along the Nankai Trough. The geometry and velocity structure of the southwestern Japan subduction zone are now well constrained. Comparing these observations with the rupture distribution of historic great thrust earthquakes, it appears that the coseismic rupture occurred along plate boundaries deeper than the wedge/backstop boundary (the boundary between the Neogene-Quaternary accretionary wedge and the crust forming the backstop). From the view of spatial relationship, both rupture distributions of the last two large events and the crust forming the backstop are considerably retreated from the trough axis in the west and east off the Kii Peninsula. In both areas, seamount or ridge subduction is apparent in seismic results, geomorphological data and geomagnetic data. The landward indentation of the deformable backstop, which corresponds to the crustal block of old accreted sediments, may be formed by seamount subduction according to published results of sandbox modeling. In particular, the subducted seamount may be a structural factor affecting the recession of the crustal block forming the backstop.     

Temporal Variations of Seismic Wave Velocity Associated with 1998 M6.1 Shizukuishi Earthquake

—?We attempt to detect temporal variations of seismic wave velocity before and after 1998 M6.1 Shizukuishi, northeastern Japan, earthquake by using waveform data from explosions and earthquake doublets spanning the period immediately before and after the earthquake. Direct P waves of the second explosion are delayed by ～20 ms at observation stations with epicentral distances less than 15 km. This tendency does not change if the analysis frequency band is changed. Our result suggests that the P-wave velocity decreased by at least 1% in the extremely shallow region of the hanging wall of the M6.1 thrust event after its occurrence. On the other hand, there was the frequency dependence of the coda wave delays for both artificial sources and for natural events. At 5–10 Hz, immediate sharp increases by more than 20 ms in time delays and lower coherency were observed at several stations. We estimated the region in which P-wave velocity might have decreased after the M6.1 earthquake. Maximum depth of the region is 13 km. The region includes the aftershock area of the M6.1 earthquake, but is eccentric to the earthquake in the west. Considering the frequency band analyzed (5–10 Hz), the scale of the spatial inhomogeneity which led to the coda wave delay is several hundreds meters. We investigated candidates for the cause of the direct P-wave and coda wave delay. Observed direct P-wave delay can be partly explained by the stress changes caused by coseismic fault slip. However, the coda wave delay cannot be explained by the stress changes that are limited to the superficial area. Crustal heterogeneity should have changed at coseismic time in the deeper area where aftershocks of the M6.1 earthquake occurred.     

Role of precursory less-viscous mixed magma in the eruption of phenocryst-rich magma: evidence from the Hokkaido-Komagatake 1929 eruption

During the 1929 activity of Hokkaido-Komagatake volcano, the Plinian eruption of a phenocryst-rich andesite was preceded by a small eruption of more mafic magma formed by magma mixing. A similar eruption sequence has been reported for some other eruptions (Pallister et al. 1996; Venezky and Rutherford 1997), suggesting that eruption of a mixed magma is a precursor of phenocryst-rich magmas. For the purpose of understanding the tapping processes of the phenocryst-rich magma chamber, we investigated the temporal variation in the erupted magma and estimated the viscosity and density of the end-member and mixed magmas with constraints drawn from petrography. For the precursory mixed magma we estimate 33dž vol.% phenocrysts, andesitic-dacitic melt composition, 3 wt.% H₂O content, and temperature of 1040°C. In comparison, for the climactic, silicic end-member magma we estimate 48Dž vol.% phenocryst, high-silica rhyolitic melt, 3 wt.% H₂O, and temperature of 950°C, respectively. The mafic end-member magma, which was not erupted, is thought to be an almost aphyric basaltic-andesitic magma, based on mass balance calculation of the phenocryst content. The proportion of the mafic end-member magma component in the mixed magma was calculated to be 20-40 wt.%. On the basis of these data, we estimate magma viscosities of 10^3.9, 10^6.9, and 10^2.0 Pa s for the mixed, silicic end-member, and mafic end-member magmas, respectively. The calculated density differences among these magmas are inconsequential when possible errors are considered. We calculate the minimum excess pressure required for dike propagation to be 31 MPa for the silicic end-member magma and 8 MPa for the mixed magma, using the estimated viscosity and dike propagation model of Rubin (1995). If we assume that excess pressure is limited by the wall rock strength of the magma chamber, excess pressure retainable in the magma chamber is less than ca. 20 MPa. This suggests that the mixed magma was able to ascend to the surface without freezing, whereas the viscous silicic end-member magma could not. The formation and precursory eruption of the mixed magma are, therefore, effective and necessary initiation processes for the phenocryst-rich, viscous magma eruption.     

Helium and carbon isotopes in fluorites: implications for mantle carbon contribution in an ancient subduction zone

The concentrations of helium and carbon in fluorite associated with Cretaceous to Neogene (90–13 Ma) granitic magmatism in the Japanese arc have been measured. Concentrations of Li, U, Th and Gd were measured to correct for secondary generated ³He. The CO₂/³He of fluorites are almost uniform (1.5×10¹⁰–4×10¹⁰) and in fair agreement with the range of present island arc volcanic gases. The calculated mantle C contribution in the Mesozoic subduction zone appear to have been identical to the present one (7–19%) indicating that the C flux from the mantle in supra-subduction zone environments has remained fairly constant during the past 70 million years.     

Semidiurnal tides from the extended Canadian Middle Atmosphere Model (CMAM) and comparisons with TIMED Doppler interferometer (TIDI) and meteor radar observations

The extended Canadian Middle Atmosphere Model (extended CMAM) is a general circulation model, which extends from the surface to about 210 km. Spatial complex spectral analysis is applied to horizontal winds simulated by the extended CMAM to obtain semidiurnal tidal amplitudes and phases (from e5 to w5) in the mesosphere and lower thermosphere (MLT) region. The dominant w2 migrating component and the presence of eight nonmigrating tides (w3, w4, w5, e1, e2, e3, e4 and e5) in the mid-latitudes are identified. Components w1 and s0, which tend to maximize at high latitudes, will be discussed separately in a later paper. The migrating semidiurnal tide (w2) has amplitudes reaching over 20 m s⁻¹ for both zonal and meridional winds in the mid-latitude region. Its form compares well to the published results. The amplitudes of nonmigrating semidiurnal tides are non-negligible compared with the migrating semidiurnal tides. The amplitudes for w3 and e2 exceed 12 and 8 m s⁻¹, respectively.Comparisons are made with four nonmigrating semidiurnal components (w3, w4, e1 and e2) derived from the TIMED Doppler interferometer (TIDI) wind measurements between 85 and 105 km altitude and between 45°S and 45°N latitude. Overall, the basic CMAM and TIDI latitudinal structures of the amplitudes agree well and the agreement between the annual mean amplitudes varies with component. Relative to the TIDI results, the CMAM seasonal variations of w4 are in good agreement, of e2 are in reasonable agreement, of w3 are in partial agreement and of e1 are in poor agreement.The 11 semidiurnal components from the model are superimposed to generate the total semidiurnal winds at Jakarta (6°S, 106°E) and Kototabang (0°, 100°E) and are compared with measurements from two equatorial meteor radar stations at these sites. The relative contributions of components to the reconstructed amplitude vary from month to month. The CMAM reconstructions are generally larger than the radar results by a factor varying between one and two. The phases in the radar data are typically stationary with respect to height, whereas they generally decrease with height in the CMAM reconstruction.