Evolution of Mount Fuji,Japan: Inference from drilling into the subaerial oldest volcano,pre‐Komitake

A buried, old volcanic body (pre‐Komitake Volcano) was discovered during drilling into the northeastern flank of Mount Fuji. The pre‐Komitake Volcano is characterized by hornblende‐bearing andesite and dacite, in contrast to the porphyritic basaltic rocks of Komitake Volcano and to the olivine‐bearing basaltic rocks of Fuji Volcano. K‐Ar age determinations and geological analysis of drilling cores suggest that the pre‐Komitake Volcano began with effusion of basaltic lava flows around 260 ka and ended with explosive eruptions of basaltic andesite and dacite magma around 160 ka. After deposition of a thin soil layer on the pre‐Komitake volcanic rocks, successive effusions of lava flows occurred at Komitake Volcano until 100 ka. Explosive eruptions of Fuji Volcano followed shortly after the activity of Komitake. The long‐term eruption rate of about 3 km³/ka or more for Fuji Volcano is much higher than that estimated for pre‐Komitake and Komitake. The chemical variation within Fuji Volcano, represented by an increase in incompatible elements at nearly constant SiO₂, differs from that within pre‐Komitake and other volcanoes in the northern Izu‐Bonin arc, where incompatible elements increase with increasing SiO₂. These changes in the volcanism in Mount Fuji may have occurred due to a change in regional tectonics around 150 ka, although this remains unproven.     

The role of the hyporheic zone for a benthic fish in an intermittent river: a refuge, not a graveyard

In lotic systems, the hyporheic zone has been suggested as a potential refuge for aquatic organisms during disturbances (hyporheic refuge hypothesis). However, the supporting evidence is unclear, especially regarding the survival of hyporheic refugees and their contribution to the recovery of post-disturbance populations. Moreover, few studies have focused on the importance of the hyporheic refuge for aquatic vertebrates such as fish. In this study, we present evidence that the hyporheic zone acts as a refuge for a small benthic fish (Cobitis shikokuensis) following surface drying in an intermittent river. We examined its survival during and recolonization after dry periods by direct hyporheic sampling and mark-and-recapture surveys. When the streambed dried, hyporheic sampling was conducted 58 times across 33 locations in the intermittent reach and 31 individuals of C. shikokuensis were captured from extracted hyporheic water. Mark-and-recapture surveys revealed that recolonizers after re-wetting included C. shikokuensis individuals that had survived dry periods in the hyporheic refuge. The condition factor of C. shikokuensis significantly declined after dry periods, suggesting that most recolonizers suffered from physiological stress, probably within the hyporheic refuge. These results clearly support the long-debated, hyporheic refuge hypothesis, and provide a striking example of the critical role of the hyporheic zone in population maintenance of lotic organisms.     

Theoretical and experimental investigation of nickel isotopic fractionation in species relevant to modern and ancient oceans

Nickel plays a central role as an enzyme co-factor in the metabolism of methanogenic Archaea. Methanogens can fractionate Ni isotopes during assimilation, opening the possibility of using the stable isotopic composition of Ni as a biomarker. However, other sources of Ni isotopic variations need to be evaluated before one can establish Ni isotopes as an unambiguous tracer of methanogenesis in the rock record. Equilibrium exchange of Ni between the different species present in the ocean is a potential source of isotopic fractionation. Through controlled laboratory experiments and theoretical calculations, we quantify equilibrium Ni isotope fractionation between different species relevant to the modern and ancient ocean: Ni(H₂O)₆²⁺, Ni(H₂O)₁₈²⁺, NiOH(H₂O)₅⁺, Ni(OH)₂(H₂O)₄, NiCl(H₂O)₅⁺, cis-NiCl₂(H₂O)₄, trans-NiCl₂(H₂O)₄, NiHS(H₂O)₅⁺, Ni(HS)₂(H₂O)₄, NiSO₄(H₂O)₄, NiHCO₃(H₂O)₄⁺, NiCO₃(H₂O)₄, and organic ligands (crown ether and oxalic acid). The magnitude of ligand-controlled Ni isotopic fractionation, approximately 1.25‰/amu (2.5‰ for the ⁶⁰Ni/⁵⁸Ni ratio), is similar to that previously measured in cultures of methanogenic Archaea.     

3-D internal architecture of methane hydrate-bearing turbidite channels in the eastern Nankai Trough, Japan

The reservoir architecture of methane hydrate (MH) bearing turbidite channels in the eastern Nankai Trough, offshore Japan is evaluated using a combination of 3-D seismic and well data. On the 3-D seismic section, the MH-bearing turbidite channels correspond to complex patterns of strong seismic reflectors, which show the 3-D internal architecture of the channel complex. A seismic-sequence stratigraphic analysis reveals that the channel complex can be roughly classified into three different stages of depositional sequence (upper, middle, and lower). Each depositional sequence results in a different depositional system that primarily controls the reservoir architecture of the turbidite channels. To construct a 3-D facies model, the stacking patterns of the turbidite channels are interpreted, and the reservoir heterogeneities of MH-bearing sediments are discussed. The identified channels at the upper sequence around the β1 well exhibit low-sinuosity channels consisting of various channel widths that range from tens to several hundreds of meters. Paleo-current flow directions of the turbidite channels are typically oriented along the north-northeast-to-south-southwest direction. High-amplitude patterns were identified above the channels along the north-to-south and north-northeast-to-south-southeast directions. These roughly coincide with the paleo-current flow of the turbidite channels. An interval velocity using high-density velocity analysis shows that velocity anomalies (>2000 m/s) are found on the northeastern side of the turbidite channels. The depositional stage of the northeastern side of the turbidite channels exhibits slightly older sediment stages than the depositional stages of the remaining channels. Hence, the velocity anomalies of the northeastern side of the channels are related to the different stages of sediment supply, and this may lead to the different reservoir architectures of the turbidite channels.     

Improving strategies with constraints regarding non-Gaussian statistics in a three-dimensional variational assimilation method

We assess validity of a Gaussian error assumption, the basic assumption in data assimilation theory, and propose two kinds of constraints regarding non-Gaussian statistics. In the mixed water region (MWR) off the east coast of Japan exhibiting complicated frontal structures, a probability density function (PDF) of subsurface temperature shows double peaks corresponding to the Kuroshio and Oyashio waters. The complicated frontal structures characterized by the temperature PDF sometimes cause large innovations, bringing about a non-Gaussianity of errors. It is also revealed that assimilated results with a standard three-dimensional variational (3DVAR) scheme have some issues in MWR, arising from the non-Gaussianity of errors. The Oyashio water sometimes becomes unrealistically cold. The double peaks seen in the observed temperature PDF are too smoothed. To improve the assimilated field in MWR, we introduce two kinds of constraints, J _c1 and J _c2, which model the observed temperature PDF. The constraint J _c1 prevents the unrealistically cold Oyashio water, and J _c2 intends to reproduce the double peaks. The assimilated fields are significantly improved by using these constraints. The constraint J _c1 effectively reduces the unrealistically cold Oyashio water. The double peaks in the observed temperature PDF are successfully reproduced by J _c2. In addition, not only subsurface temperature but also whole level temperature and salinity (T–S) fields are improved by adopting J _c1 and J _c2 to a multivariate 3DVAR scheme with vertical coupled T–S empirical orthogonal function modes.     

Ab initio calculations of the Fe(II) and Fe(III) isotopic effects in citrates,nicotianamine, and phytosiderophore,and new Fe isotopic measurements in higher plants

Iron is one of the most abundant transition metal in higher plants and variations in its isotopic compositions can be used to trace its utilization. In order to better understand the effect of plant-induced isotopic fractionation on the global Fe cycling, we have estimated by quantum chemical calculations the magnitude of the isotopic fractionation between different Fe species relevant to the transport and storage of Fe in higher plants: Fe(II)-citrate, Fe(III)-citrate, Fe(II)-nicotianamine, and Fe(III)-phytosiderophore. The ab initio calculations show firstly, that Fe(II)-nicotianamine is ～3‰ (⁵⁶Fe/⁵⁴Fe) isotopically lighter than Fe(III)-phytosiderophore; secondly, even in the absence of redox changes of Fe, change in the speciation alone can create up to ～1.5‰ isotopic fractionation. For example, Fe(III)-phytosiderophore is up to 1.5‰ heavier than Fe(III)-citrate₂ and Fe(II)-nicotianamine is up to 1‰ heavier than Fe(II)-citrate. In addition, in order to better understand the Fe isotopic fractionation between different plant components, we have analyzed the iron isotopic composition of different organs (roots, seeds, germinated seeds, leaves and stems) from six species of higher plants: the dicot lentil (Lens culinaris), and the graminaceous monocots Virginia wild rye (Elymus virginicus), Johnsongrass (Sorghum halepense), Kentucky bluegrass (Poa pratensis), river oat (Uniola latifolia), and Indian goosegrass (Eleusine indica). The calculations may explain that the roots of strategy-II plants (Fe(III)-phytosiderophore) are isotopically heavier (by about 1‰ for the δ⁵⁶Fe) than the upper parts of the plants (Fe transported as Fe(III)-citrate in the xylem or Fe(II)-nicotianamine in the phloem). In addition, we suggest that the isotopic variations observed between younger and older leaves could be explained by mixing of Fe received from the xylem and the phloem.     

Nuclear field shift in natural environments

The nuclear field shift (NFS) is an isotope shift in atomic energy levels caused by a combination of differences in nuclear size and shape and electron densities at the nucleus. The effect of NFS in isotope fractionation was theoretically established by Bigeleisen in 1996 [Bigeleisen J. (1996) J. Am. Chem. Soc. 118:3676–3680] and has been analytically measured in laboratory chemical exchange reactions. More recently, some isotopic variations of heavy elements (Hg, Tl, U) measured in natural systems as well as isotopic anomalies measured for lower-mass elements in meteorites have been attributed to the NFS effect. These isotopic variations open up new and exciting fields of investigations in Earth sciences. In this paper, we review the different natural systems in which NFS has been proposed to be the origin of isotopic variations.     

SHRIMP zircon and EPMA monazite dating of granitic rocks from the Maizuru terrane, southwest Japan: Correlation with East Asian Paleozoic terranes and geological implications

Abstract The Maizuru terrane, distributed in the Inner Zone of southwest Japan, is divided into three subzones (Northern, Central and Southern), each with distinct lithological associations. In clear contrast with the Southern zone consisting of the Yakuno ophiolite, the Northern zone is subdivided into the western and eastern bodies by a high-angle fault, recognized mainly by the presence of deformed granitic rocks and pelitic gneiss. This association suggests an affinity with a mature continental block; this is supported by the mode of occurrence, and petrological and isotopic data. Newly obtained sensitive high mass-resolution ion microprobe (SHRIMP) zircon U–Pb ages reveal the intrusion ages of 424 ± 16 and 405 ± 18 Ma (Siluro–Devonian) for the granites from the western body, and 249 ± 10 and 243 ± 19 Ma (Permo–Triassic) for the granodiorites from the eastern body. The granites in the western body also show inherited zircon ages of around 580 and 765 Ma. In addition, electron probe microanalysis (EPMA) monazite U–Th–total Pb dating gives around 475–460 Ma. The age of intrusion, inherited ages, mode of occurrence, and geological setting of the Siluro–Devonian granites of the Northern zone all show similarities with those of the Khanka Massif, southern Primoye, Russia, and the Hikami granitic rocks of the South Kitakami terrane, Northeast Japan. We propose that both the Siluro–Devonian and Permo–Triassic granitic rocks of the Northern zone are likely to have been juxtaposed through the Triassic–Late Jurassic dextral strike-slip movement, and to have originated from the Khanka Massif and the Hida terrane, respectively. This study strongly supports the importance of the strike-slip movement as a mechanism causing the structural rearrangement of the Paleozoic–Mesozoic terranes in the Japanese Islands, as well as in East Asia.     

Perturbed physics ensemble using the MIROC5 coupled atmosphere–ocean GCM without flux corrections: experimental design and results

In this study, we constructed a perturbed physics ensemble (PPE) for the MIROC5 coupled atmosphere–ocean general circulation model (CGCM) to investigate the parametric uncertainty of climate sensitivity (CS). Previous studies of PPEs have mainly used the atmosphere-slab ocean models. A few PPE studies using a CGCM applied flux corrections, because perturbations in parameters can lead to large radiation imbalances at the top of the atmosphere and climate drifts. We developed a method to prevent climate drifts in PPE experiments using the MIROC5 CGCM without flux corrections. We simultaneously swept 10 parameters in atmosphere and surface schemes. The range of CS (estimated from our 35 ensemble members) was not wide (2.2–3.2?°C). The shortwave cloud feedback related to changes in middle-level cloud albedo dominated the variations in the total feedback. We found three performance metrics for the present climate simulations of middle-level cloud albedo, precipitation, and ENSO amplitude that systematically relate to the variations in shortwave cloud feedback in this PPE.     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.