Detection and Evaluation of Gas Hydrates in the Eastern Nankai Trough by Geochemical and Geophysical Methods

Abstract: Interstitial waters extracted from the sediment cores from the exploration wells, “BH‐1” and “MITI Nankai Trough”, drilled ～60 km off Omaezaki Peninsula in the eastern Nankai Trough, were analyzed for the chloride and sulfate concentrations to examine the depth profiles and occurrence of subsurface gas hydrates. Cored intervals from the seafloor to 310 mbsf were divided into Unit 1 (～70 mbsf, predominated by mud), Unit 2 (70–150 mbsf, mud with thin ash beds), Unit 3 (150–250+ mbsf, mud with thin ash and sand), and Unit 4 (275–310 mbsf, predominated by mud). The baseline level for Cl “concentrations was 540 mM, whereas low chloride anomalies (103 to 223 mM) were identified at around 207 mbsf (zone A), 234–240 mbsf (zone B), and 258–265 mbsf (zone C) in Unit 3. Gas hydrate saturation (Sh %) of sediment pores was calculated to be 60 % (zone A) to 80 % (zones B and C) in sands whereas only a few percent in clay and silt. The total amount of gas hydrates in hydrate‐bearing sands was estimated to be 8 to 10 m³ of solid gas hydrate per m², or 1.48 km³ CH₄ per 1 km². High saturation zones (A, B and C) were consistent with anomaly zones recognized in sonic and resistivity logs. 2D and high‐resolution seismic studies revealed two BSRs in the study area. Strong BSRs (BSR‐1) at ～263 mbsf were correlated to the boundary between gas hydrate‐bearing sands (zone C) and the shallower low velocity zone, while the lower BSRs (BSR‐2) at～289 mbsf corresponded to the top of the deeper low velocity zone of the sonic log. Tectonic uplift of the study area is thought to have caused the upward migration of BGHS. That is, BSR‐1 corresponds to the new BGHS and BSR‐2 to the old BGHS. Relic gas hydrates and free gas may survive in the interval between BSR‐1 and BSR‐2, and below BSR‐2, respectively. Direct measurements of the formation temperature for the top 170 m interval yield a geothermal gradient of ～4.3d?C/ 100 m. Extrapolation of this gradient down to the base of gas hydrate stability yields a theoretical BGHS at～230 mbsf, surprisingly ～35 m shallower than the base of gas hydrate‐bearing sands (zone C) and BSR‐1. As with the double BSRs, another tectonic uplift may explain the BGHS at unreasonably shallow depths. Alternatively, linear extrapolation of the geothermal gradient down to the hydrate‐bearing zones may not be appropriate if the gradient changes below the depths that were measured. Recognition of double BSRs (263 and 289 mbsf) and probable new BGHS (～230 mbsf) in the exploration wells implies that the BGHS has gradually migrated upward. Tectonically induced processes are thought to have enhanced dense and massive accumulation of gas hydrate deposits through effective methane recycling and condensation. To test the hypothetical models for the accumulation of gas hydrates in Nankai accretionary prism, we strongly propose to measure the equilibrium temperatures for the entire depth range down to the free gas zone below predicted BGHS and to reconstruct the water depths and uplift history of hydrate‐bearing area.     

Jovian magnetosphere-ionosphere current system characterized by diurnal variation of ionospheric conductance

We developed a new numerical model of the Jovian magnetosphere-ionosphere coupling current system in order to investigate the effects of diurnal variation of ionospheric conductance. The conductance is determined by ion chemical processes that include the generation of hydrogen and hydrocarbon ions by solar EUV radiation and auroral electrons precipitation. The model solves the torque equations for magnetospheric plasma accelerated by the radial currents flowing along the magnetospheric equator. The conductance and magnetospheric plasma then change the field-aligned currents (FACs) and the intensity of the electric field projected onto the ionosphere. Because of the positive feedback of the ionospheric conductance on the FAC, the FAC is the maximum on the dayside and minimum just before sunrise. The power transferred from the planetary rotation is mainly consumed in the upper atmosphere on the dayside, while it is used for magnetospheric plasma acceleration in other local time (LT) sectors. Further, our simulations show that the magnetospheric plasma density and mass flux affect the temporal variation in the peak FAC density. The enhancement of the solar EUV flux by a factor of 2.4 increases the FAC density by 30%. The maximum density of the FAC is determined not only by the relationship between the precipitating electron flux and ionospheric conductance, but also by the system inertia, i.e., the inertia of the magnetospheric plasma. A theoretical analysis and numerical simulations reveal that the FAC density is in proportion to the planetary angular velocity on the dayside and to the square of the planetary angular velocity on the nightside. When the radial current at the outer boundary is fixed at values above 30 MA, as assumed in previous model studies, the peak FAC density determined at latitude 73°-74° is larger than the diurnal variable component. This result suggests large effects of this assumed radial current at the outer boundary on the system.     

Laboratory Study of O(1S) Formation Process in the Photolysis of O3 and its Atmospheric Implications

A high-sensitive technique to detect O(¹S) atoms using vacuum ultraviolet laser-induced fluorescence (VUV-LIF) spectroscopy has been applied to study the O(¹S) production process from the UV photodissociation of O₃, N₂O, and H₂O₂. The quantum yields for O(¹S) formation from O₃ photolysis at 215 and 220 nm are determined to be (1.4 ± 0.4) × 10⁻⁴ and (5 ± 3) × 10⁻⁵, respectively. Based on thermochemical considerations, the O(¹S) formation from O₃ photolysis at 215 and 220 nm is attributed to a spin-forbidden process of O(¹S)+O₂(X³Σ_g ⁻). Analysis of the Doppler profile of O(¹S) produced from O₃ photolysis at 193 nm also indicates that the O(¹S) atoms are produced from the spin-forbidden process. In the photolysis of N₂O and H₂O₂ at 193 nm, no discernible signal of O(¹S) atoms has been detected. The upper limit values of the quantum yields for O(¹S) production from N₂O and H₂O₂ photolysis at 193 nm are estimated to be 8 × 10⁻⁵ and 3 × 10⁻⁵, respectively. Using the experimental results, the impact of the O(¹S) formation from O₃ photolysis on the atmospheric OH radical formation through the reaction of O(¹S)+H₂O has been estimated. The calculated results show that the contribution of the O(¹S)+H₂O reaction to the OH production rate is ∼2% of that of the O(¹D)+H₂O reaction at 30 km altitude in mid-latitude. Implications of the present laboratory experimental results for the terrestrial airglow of O(¹S) at 557.7 nm have also been discussed.     

Topography of the Moho and Conrad discontinuities in the Kyushu district,Southwest Japan

The first P-arrival time data from local earthquakes are inverted for two-dimensional variation of the depths to the Conrad and Moho discontinuities in the Kyushu district, southwest Japan. At the same time, earthquake hypocenters and station corrections are determined from the data. The depths to the discontinuities are estimated by minimizing the travel time residuals of first P-arrival phases for 608 earthquakes observed at 57 seismic stations. In the land area of Kyushu, the Conrad and Moho discontinuities are located within the depth ranges of 16–18 and 34–40 km, respectively. The Conrad discontinuity is not as largely undulated as the Moho discontinuity. The depth to the Moho is deep along the east coast of Kyushu, and the deepest Moho is closely related to markedly low velocity of P wave. We regard the deepest Moho as reflecting the Kyushu–Palau ridge subducting beneath the Kyushu district, together with the Philippine Sea slab. In western Kyushu, the shallow Moho is spreading in the north–northeast–south–southwest direction in the Okinawa trough region. Based on the presence of low-velocity anomaly in three-dimensional velocity structure and seismogenic stress field of shallow crustal earthquakes, the shallow Moho is interpreted as being due to lower crustal erosion associated with a small-scale mantle upwelling in the Okinawa trough region. The velocity discontinuity undulation basically has insignificant effect on hypocenter determination of the local earthquakes, but the Moho topography makes changes in focal depths of some upper mantle earthquakes. The depth variation of the Moho discontinuity has a good correlation with the Bouguer gravity anomaly map; i.e., the shallow Moho of western Kyushu and the deep Moho of eastern Kyushu closely correlate with the positive and negative Bouguer gravity anomalies, respectively.     

Free-surface implementation in a mesh-free finite-difference method for elastic wave propagation in the frequency domain

We present an original implementation of the free-surface boundary condition in a mesh-free finite-difference method for simulating elastic wave propagation in the frequency domain. For elastic wave modelling in the frequency domain, the treatment of free surfaces is a key issue which requires special consideration. In the present study, the free-surface boundary condition is directly implemented at node positions located on the free-surface. Flexible nature of the mesh-free method for nodal distribution enables us to introduce topography into numerical models in an efficient manner. We investigate the accuracy of the proposed implementation by comparing numerical results with an analytical solution. The results show that the proposed method can calculate surface wave propagation even for an inclined free surface with substantial accuracy. Next, we calculate surface wave propagation in a model with a topographic surface using our method, and compare the numerical result with that using the finite-element method. The comparison shows the excellent agreement with each other. Finally, we apply our method to the SEG foothill model to investigate the effectiveness of the proposed method. Since the mesh-free method has high flexibility of nodal distribution, the proposed implementation would deal with models of topographic surface with sufficient accuracy and efficiency.     

Effect of occurrence rate of acoustic emissions on their statistical behavior

The parameterm in Ishimoto-Iida's relation was investigated for acoustic emissions (AEs) occurring in rock samples under uniaxial compression. In the experiment, we found: 1) The large AEs are counted without serious error but the number of small AEs is systematically underestimated at high AE rates, 2) the frequency distribution of maximum AE amplitudes becomes nonlinear in logarithmic scale with increasing AE rate, and 3) there exists a strong negative correlation betweenm-value and AE rate. The miscount of small AEs was interpreted as due to overlap of the large and small AEs. We call the miscount masking effect. A statistical analysis based on the masking effect showed that them-value decreases more effectively as the AE rate increases, and thus the masking effect is a possible origin both for the nonlinear frequency distribution of maximum AE amplitudes and for the negative correlation ofm-value with AE rate. We emphasize that one should be careful of the masking effect to examine correctly the change, ofm-value. In order to eliminate the masking effect, AEs should be measured by a measurement system with low sensitivity. Even if the masking effect is eliminated, them-value decreases before the main fracture of a rock sample. Them-value is a key parameter to predict the main fracture.     

Summary of taxa and distribution of Sirenia in the North Pacific Ocean

Abstract North Pacific fossil sirenians comprise representatives of three subfamilies of the Dugongidae: Halitheriinae ( Metaxytherium arctodites , Middle Miocene, North America), Hy-drodamalinae ( Dusisiren spp., Early-Late Miocene, and Hydrodamalis spp., Late Miocene-Pleistocene, North America and Japan), and Dugonginae ( Dioplotherium allisoni , Early-Middle Miocene, North America). Indeterminate dugongid remains are also known from the Late Oligocene of Japan, and the discovery of additional taxa in the western Pacific, especially in Paleogene rocks, can be anticipated. The known North Pacific Neogene taxa apparently dispersed into the Pacific from the Caribbean. Metaxytherium gave rise in the Pacific to Dusisiren ; a series of chronospecies of the latter genus eventually culminated in Hydrodamalis , which was exterminated by humans circa AD 1768. Dioplotherium left no known descendants in the Pacific. The Recent Dugong probably entered the Pacific from the Indian Ocean. The presence in the North Pacific Miocene of at least three sympatric dugongid lineages, together with desmostylians, is evidence for a diversity of marine plants that was reduced by subsequent climatic cooling.     

Geological and microbial anomalies in the extinct submarine volcano, Shiribeshi Seamount, in the eastern margin of the Japan Sea

Abstract   The manned submersible Shinkai 2000 investigated yellow patches on the near-summit slope of Shiribeshi Seamount in the Japan Sea. Yellowish patches are often associated with seepage, and the possibility of seepage at Shiribeshi Seamount was tested by the following four lines of observation: (i) high subsurface temperature was measured at a ring-like patch, although no increase in subsurface temperature was observed at other patches; (ii) high gamma ray (γ-ray) intensity from the thorium series was recorded in the patch zone; (iii) the yellowish deposit was composed of calcite, quartz and amorphous iron compound, as seen at the yellowish patches in other seeps and volcanoes; (iv) lipid phosphate, a measure of microbial abundance, in sediments of the ring-like patch was determined, and the recorded microbial abundance was higher inside the patch than outside it. The four lines of observation are explained consistently by postulating that the seepage of warm fluid contained Fe and γ-ray sources. A hydrothermal origin of the yellow patches is not ruled out for the extinct but young (0.9 Ma) arc volcano.     

The Tiger Sulfide Chimney,Yonaguni Knoll IV Hydrothermal Field,Southern Okinawa Trough,Japan: The First Reported Occurrence of Pt–Cu–Fe‐Bearing Bismuthinite and Sn‐Bearing Chalcopyrite in an Active Seafloor Hydrothermal System

A sulfide chimney ore sampled from the flank of the active Tiger vent area in the Yonaguni Knoll IV hydrothermal field, south Okinawa trough, consists of anhydrite, pyrite, sphalerite, galena, chalcopyrite and bismuthinite. Electron microprobe analysis indicates that the chalcopyrite contains up to 2.4 wt% Sn, whereas bismuthinite contains up to 1.7 wt% Pt, 0.8 wt% Cu and 0.5 wt% Fe. The Sn‐rich chalcopyrite and Pt–Cu–Fe‐bearing bismuthinite are the first reported occurrence of such minerals in an active submarine hydrothermal system. The results confirm that Sn enters the chalcopyrite as a solid solution towards stannite by the coupled substitution of Sn⁴⁺Fe²⁺ for Fe³⁺Fe³⁺, whereas Pt, Cu and Fe enter the bismuthinite structure as a solid solution during rapid nucleation. The fluid inclusions homogenization temperatures in anhydrite (220–310°C) and measured end‐member temperature of the vent fluids on‐site (325°C) indicate that Sn‐bearing chalcopyrite and Pt–Cu–Fe‐bearing bismuthinite express the original composition of the minerals that precipitated as metastable phases at a temperature above 300°C. The result observed in this study implies that sulfides in ancient volcanogenic massive sulfide deposits have similar trace element distribution during nucleation but it is remobilised during diagenesis, metamorphism or supergene enrichment processes.     

Advanced Micro-XRF Method to Separate Sedimentary Rhythms and Event Layers in Sediments: Its Application to Lacustrine Sediment from Lake Suigetsu,Japan

Event-related sedimentary layers, which are deposited occasionally due to volcanic eruptions, earthquakes or heavy rains, are often contained in the rhythmical sequences of lacustrine and marine sediments. We have developed an analytical method for separating the sedimentary rhythms and the event layers identified using the scanning X-ray analytical microscope (SXAM) and obtained sequential profiles of seven elements Al, Si, K, Ca, Ti, Mn, and Fe in the lacustrine sediment from Lake Suigetsu, Japan. Two types of event layers could be detected from the elemental composition of 33 layers of sediment: three known volcanic ash layers and 30 clay layers containing 12 turbidites. The recurrence interval of the latter, which may potentially be initiated and archived by locally important earthquakes, is estimated to be an average of 640 ± 160 years by using Sompi event analysis (SEA) based on an autoregressive (AR) model. After removing those portions that represented event layers from the elemental profiles, we obtained event-removed (ER) temporal profiles based on the tephrochronology of the three volcanic ash layers. The ER temporal profiles of manganese and iron, probably representing the siderite content, showed a millennial-scale variation in the Holocene that corresponded well with ice-rafting events in the North Atlantic.