An efficient method to compute the dynamic response of a fluid-filled crack

The vibration of a fluid-filled crack is considered to be one of the most plausible source mechanisms for the long-period events and volcanic tremors occurring around volcanoes. As a tool for the quantitative interpretation of source process of such volcanic seismic signals, we propose a method to numerically simulate the dynamic response of a fluid-filled crack. In this method, we formulate the motions of the fluid inside and the elastic solid outside of the crack, using boundary integrals in the frequency domain and solve the dynamic interactions between the fluid and the elastic solid using the point collocation method. The present method is more efficient compared with the time-domain finite difference method, which has been used in simulations of a fluid-filled crack and enables us to study the dynamics of a fluid-filled crack over a wide range of physical parameters. The method also allows us direct calculation of the attenuation quality factor of the crack resonance, which is an indispensable parameter for estimating the properties of the fluid inside the crack. The method is also designed to be flexible to many applications, which may be encountered in volcano seismology, and thus, extensions of the method to more complicated problems are promising.     

Gas hydrate saturation at Site C0002, IODP Expeditions 314 and 315, in the Kumano Basin,Nankai trough

The degree of gas hydrate saturation at Integrated Ocean Drilling Program (IODP) Site C0002 in the Kumano Basin, Nankai Trough, was estimated from logging‐while‐drilling logs and core samples obtained during IODP Expeditions 314 and 315. Sediment porosity data necessary for the calculation of saturation were obtained from both core samples and density logs. Two forms of the Archie equation (‘quick‐look’ and ‘standard’) were used to calculate gas hydrate saturation from two types of electrical resistivity log data (ring resistivity and bit resistivity), and a three‐phase Biot‐type equation was used to calculate gas hydrate saturation from P‐wave velocity log data. The gas hydrate saturation baseline calculated from both resistivity logs ranges from 0% to 35%, and that calculated from the P‐wave velocity log ranges from 0% to 30%. High levels of gas hydrate saturation (>60%) are present as spikes in the ring resistivity log and correspond to the presence of gas hydrate concentrations within sandy layers. At several depths, saturation values obtained from P‐wave velocity data are lower than those obtained from bit resistivity data; this discrepancy is related to the presence of free gas at these depths. Previous research has suggested that gas from deep levels in the Kumano Basin has migrated up‐dip towards the southern and seaward edge of the basin near Site C0002. The high saturation values and presence of free gas at site C0002 suggest that a large gas flux is flowing to the southern and seaward edge of the basin from a deeper and/or more landward part of the Kumano Basin, with the southern edge of the Kumano Basin (the location of site C0002) being the main area of fluid accumulation.     

Preferred orientation of the trace fossil Entradichnus meniscus in eolian dune strata (Djadokhta Formation) at Tugrikiin Shiree,southern Mongolia and its paleoecological implications

Eolian sand dune deposits of the Upper Cretaceous Djadokhta Formation at Tugrikiin Shiree, southern Mongolia, yield not only dinosaur skeletal remains but also numerous trace fossils produced by invertebrates. This paper describes the trace fossil Entradichnus meniscus, a long unlined and unbranched trail that is filled with meniscate laminae and occurs characteristically in positive epirelief. The trail is straight to gently meandering, parallel to the foreset laminae of the eolian dunes, and exhibits a significant preferred orientation parallel to the depositional dip of the cross‐stratification laminae. In addition, almost all the crescentic internal laminae of the trail show concave down‐dips. These features indicate that the trails were produced beneath the slipface of eolian dunes by the downward burrowing of the trace‐makers. This occurrence mode of E. meniscus of the Mongolian Cretaceous is very similar to that described from the Jurassic eolian dune deposits in North America. Hence, the downward burrowing of the E. meniscus animal might be a common feature in arid eolian dune deposits at least during the Jurassic and Cretaceous, and possibly reflecting a behavioral response to the morphology of large sand dunes under an arid climate.     

Spatial variations of Sr and Nd isotope ratios of Cretaceous-Paleogene granitoid rocks,Southwest Japan Arc

Cretaceous-Paleogene granitoid rocks and contemporaneous volcanic rocks are widely distributed in the Inner Zone of Southwest Japan. This intense intermediate to felsic magmatism is considered to have taken place on the eastern margin of the Eurasian Continent, before the Southwest Japan Arc drifted away from the continent in the middle Miocene, resulting in the opening of the Japan Sea. The granitoid rocks show regional variations in terms of their radiometric age, petrography, Sr, Nd and O isotope ratios. Based on Sr and Nd isotope ratios, granitoid rocks can be divided into three zones (South, Transitional and North) between the Median Tectonic Line and the Japan Sea. Granitoid rocks and associated gabbros of the North Zone have low initial Sr isotope ratios (0.7048 to 0.7068) and high initial Nd values (+3 to-2.2), whereas granitoid rocks and gabbros from the South Zone have high initial Sr isotope ratios (0.7070 to 0.7088) and low initial Nd values (-3.0to-8.0). Most granitoid rocks from the Transitional Zone have Sr and Nd isotope ratios that lie between those of the North and South Zones, although there is some overlap. Contamination of magmas by upper crust cannot explain this geographical variation in Sr and Nd isotopes. Instead, the regional variation is attributed to compositionally different, magma sources (probably upper mantle and lower crust), beneath the North and South Zones. This is supported by the Sr and Nd isotopic ratios of upper mantle and lower crustal xenoliths included in Cenozoic volcanic rocks in the North and South Zones. These ratios are similar to those of the granitoid rocks in the respective zones. It is suggested that a micro-continent or island arc consisting of continental crust was underthrust beneath the South Zone before or during the Cretaceous, resulting in compositionally distinct sources for granitoid rocks of the North and South Zones. The large variation observed in Sr and Nd isotope ratios for Transitional Zone granitoid rocks is explained by variable proportions of the two different crustal and upper mantle components.     

Isotopic composition of gas hydrates in subsurface sediments from offshore Sakhalin Island, Sea of Okhotsk

Hydrate-bearing sediment cores were retrieved from recently discovered seepage sites located offshore Sakhalin Island in the Sea of Okhotsk. We obtained samples of natural gas hydrates and dissolved gas in pore water using a headspace gas method for determining their molecular and isotopic compositions. Molecular composition ratios C₁/_C2+ from all the seepage sites were in the range of 1,500–50,000, while δ¹³C and δD values of methane ranged from ?66.0 to ?63.2‰ VPDB and ?204.6 to ?196.7‰ VSMOW, respectively. These results indicate that the methane was produced by microbial reduction of CO₂. δ¹³C values of ethane and propane (i.e., ?40.8 to ?27.4‰ VPDB and ?41.3 to ?30.6‰ VPDB, respectively) showed that small amounts of thermogenic gas were mixed with microbial methane. We also analyzed the isotopic difference between hydrate-bound and dissolved gases, and discovered that the magnitude by which the δD hydrate gas was smaller than that of dissolved gas was in the range 4.3–16.6‰, while there were no differences in δ¹³C values. Based on isotopic fractionation of guest gas during the formation of gas hydrate, we conclude that the current gas in the pore water is the source of the gas hydrate at the VNIIOkeangeologia and Giselle Flare sites, but not the source of the gas hydrate at the Hieroglyph and KOPRI sites.     

Enhanced predictions of wave impact pressure by improved incompressible SPH methods

A new criterion is proposed for a more efficient assessment of free-surface particles in a particle-based simulation. Enhanced wave impact simulations are carried out by improved Incompressible SPH (ISPH) methods. The first improvement is the same as that in the Corrected ISPH (CISPH; [Khayyer A, Gotoh, H, Shao SD. Corrected incompressible SPH method for accurate water-surface tracking in breaking waves, Coast Eng 2008; 55 (3): 236–250]) method and is proposed for the improvement of momentum conservation. The second improvement is achieved by deriving and employing a higher order source term based on a more accurate differentiation to obtain a less fluctuating and more accurate pressure field. The enhanced performance of improved ISPH methods is demonstrated through the simulation of several fluid impact simulations in comparison with the experimental data and simulation results by other numerical methods.     

Upwelling of oxygen-depleted water (Sumishio) in Omura Bay,Japan

Wind-induced circulation and the distribution of hypoxia corresponding to the upwelling of oxygen-depleted water (called “Aoshio” in Japan and “Sumishio” locally in Omura Bay) in Omura Bay, Japan, was examined field observations and by three-dimensional modeling. During the calm weather in summer, well-mixed strait water, in rich oxygen at the mouth of the bay intruded into the middle layer of the bay, southward and northward along the west and east coast, respectively, forming basin-scale cyclonic circulation. A stagnant water mass was formed below the center of this cyclonic circulation, and it became hypoxic water. As a result of the prevailing strong southeast (northeast) wind, the bottom hypoxia moved in a southeasterly (northeasterly) direction. This induces the upwelling of hypoxic water, accompanied by mass mortality of marine organisms.