The Kiryu Experimental Watershed: 50-years of rainfall-runoff data for a forest catchment in central Japan

The Kiryu Experimental Catchment (KEW) is a small (5.99 ha) forest catchment located in Shiga Prefecture, central Japan (34°58′ N, 136°00′ E; www.bluemoon.kais.kyoto-u.ac.jp/kiryu/contents.html ). Around this area, forest devastation occurred from ca. 1250 to ca. 150 years ago because of overuse of forest and timbers. Then, hillside forestation was carried out for more than 100 years to prevent soil erosion and support the timber industry, and consequently, most of this area is now covered with plantation forests mainly by Chamaecyparis obtusa Sieb. et Zucc. (Japanese cypress) planted around 1960's. This plantation forest is not actively managed. The KEW is one of the leading experimental forests with long-term monitoring data in Japan. Research in the KEW began in 1967 to elucidate the hydrological and biogeochemical processes in the forested catchment in relation to climate, geology, soil, and vegetation growth. Since then, the long-term hydrological data of precipitation, runoff and sediment transport are continuously monitoring. In this study, we provide the data and preliminarily discuss the rainfall–runoff patterns and sediment transport through 50 years in the KEW. The annual precipitation and the maximum daily rainfall have been greater than the average over the last decade. In response to the rainfall patterns, the ratio of annual direct runoff to precipitation was also larger in the last decade. The sediment transport in this decade was consequently larger than the preceding decades. Our data presented here suggest that a close relationship exists between the climate condition, rainfall–runoff response, sediment dynamics, as well as a slowly progressing change of forest condition.     

Molecular dynamics study of the crystal structure and phase relation of the GeO2 polymorphs

A two-body interatomic potential model for GeO₂ polymorphs has been determined to simulate the structure change of them by semi-empirical procedure, total lattice energy minimization of GeO₂ polymorphs. Based on this potential, two polymorphs of GeO₂; α-quartz-type and rutile-type, have been reproduced using the molecular dynamics (MD) simulation techniques. Crystal structures, bulk moduli, volume thermal expansion coefficients and enthalpies of these polymorphs of GeO₂ were simulated. In spite of the simple form of the potential, these simulated structural values, bulk moduli and thermal expansivities are in excellent agreement with the reliable experimental data in respect to both polymorphs. Using this potential, MD simulation was further used to study the structural changes of GeO₂ under high pressure. We have investigated the pressure-induced amorphization. As reported in previous experimental studies, quartz-type GeO₂ undergoes pressure-induced crystalline-to-amorphous transformation at room temperature, the same as other quartz compounds; SiO₂, AlPO₄. Under hydrostatic compression, in this study, α-quartz-type GeO₂ transformed to a denser amorphous state at 7.4 GPa with change of the packing of oxygen ions and increase of germanium coordination. At higher pressure still, rutile-type GeO₂ transformed to a new phase of CaCl₂-type structure as a post-rutile candidate. Received: 29 July 1996 / Revised, accepted: 30 April 1997     

Eccentric instability of a two-dimensional fluid disc with a central massive object

In this paper we investigate, by linear modal analysis, the one-armed dynamical instability of a two-dimensional fluid disc that has a massive object at its centre. The model of the disc is chosen to avoid the artificial instabilities that originate from the unrealistic disc configurations that have been adopted in previous studies. We find a one-armed instability for which the central massive object is displaced from the centre, which is generally called the 'eccentric instability'. However, to excite the eccentric instability, the mass of the central object should be appreciably smaller than that of the disc, and this mass ratio is far smaller than what was originally proposed. The instability shown in this paper is likely to be excited in a stellar system with a central massive object, e.g. a galactic nucleus harbouring a massive black hole, and further studies are desirable via techniques such as numerical simulations.     

Computational model of the structural and elastic properties of the ilmenite and perovskite phases of MgSiO3

The structural and elastic properties of the ilmenite and perovskite phases of MgSiO₃ are investigated with a computational model based on energy minimization. The potential energies of these two crystals are approximated by the sum of Coulomb, van der Waals, and repulsion terms between atoms. Required energy parameters are derived by fitting the parameters to the observed crystal structures of these two phases as well as to the measured elastic constants of the ilmenite phase. The resulting potential model is applied to predicting the elastic constants of the perovskite phase. The calculated bulk modulus of the perovskite phase compares favorably with the data obtained from volume-compression experiments as well as the values estimated from empirical elasticity systematics of perovskite type compounds. The predicted shear modulus of the perovskite phase is also in reasonable agreement with the values proposed from similar empirical elasticity systematics. Subsequently, the model is used to simulate the high pressure behaviors of the crystal structures and elastic constants of these two phases.     

Comparison between thermal–viscous coupling and frictional sliding

We investigated the similarity between thermal–viscous coupling (TVC) and frictional sliding, proposed by Kameyama and Kaneda [Pure Appl. Geophys. 159 (2002) 2011]. We consider a one-dimensional layer composed of viscous material, which is sandwiched and sheared by two thick elastic layers. The rate of viscous deformation depends on the temperature T_c in the viscous layer as well as shear stress τ. The temperature T_c changes owing to heating by viscous dissipation and conductive cooling. We carried out velocity-stepping tests for the steady-state deformation both numerically and analytically, and compared the temporal evolution of small perturbations with that of the spring-block model with rate- and state-dependent friction (RSF). We found that, as is the case of frictional slip stability, the manner of temporal evolution is classified into four regimes depending on whether it is stable or not and whether it is monotonous or oscillatory with time. By further interpreting TVC in terms of general RSF theory by Ruina [J. Geophys. Res. 88 (1983) 10359], we obtained the relations between the parameters appearing in the phenomenological RSF law and the nondimensional parameters which characterize the nature of TVC. A further improvement of this approach might be important for estimating the actual values of frictional constitutive parameters at the deeper portion of seismogenic faults of interplate or inland earthquakes where a ductile deformation is expected to be significant.     

Stable hydrogen isotope measurement of archaeal ether-bound hydrocarbons

Compound-specific hydrogen isotope analysis of ether-bound isoprenoid hydrocarbons from archaeal membranes has been developed using chemical degradation and gas chromatography/pyrolysis/isotope ratio mass spectrometry. The ether-bound hydrocarbons are quantitatively converted to saturated hydrocarbons by cleavage of ether bonds with HI followed by H₂ reduction in the presence of PtO₂. The δD value of ether-bound hydrocarbon moieties are corrected by way of isotopic mass balance calculation for the hydrogen incorporated during the hydrogenation. The method was successfully applied to determination of the δD values of biphytane moieties in glycerol dialkyl glycerol tetraethers (GDGTs) from a Sulfolobus sp. culture and a marine sediment.     

The influence of the brightness of the asteroidal dust bands on the gegenschein

The position and shape of the Gegenschein’s maximum brightness provide information on the structure of the interplanetary dust cloud. We show that the asteroidal dust bands, extended near the anti-solar point, play an important role in determining both the position of the maximum brightness and the shape of the Gegenschein. After removing the asteroidal dust bands from an observation of the Gegenschein on November 2, 1997, it was found that the maximum brightness point shifted −0.4° in ecliptic latitude, i.e., to the south of the ecliptic plane, at an ecliptic longitude of 180°, in contrast to a latitude value of +0.1° when the dust bands were included. Furthermore, the part of the Gegenschein to the south of the ecliptic plane was brighter than the northern part at the time of observation. Referring to the cloud model of T. Kelsall et al. (1998, Astrophy. J. 508, 44-73), it can be estimated that the ascending node of the symmetry plane of the dust cloud is 57°^−3°_+7° when its inclination is 2.03° ? 0.50°.     

High-temperature Raman spectroscopy and quasi-harmonic lattice dynamic simulation of diopside

We investigated the lattice vibrational properties and lattice dynamical behaviour of diopside by combining laser micro-Raman spectroscopic measurements with quasi-harmonic lattice dynamic simulation using a transferable interatomic potential. We obtained polarized Raman spectra from a Fe-poor natural diopside and the temperature dependencies of the Raman modes to 1125?K from high-temperature Raman spectra of a Fe-poor and a Fe-rich natural diopside. The various modes display different temperature dependencies: from ?0.021?cm^?1/K to ?0.004?cm^?1/K. The temperature shift of low frequency modes is generally higher. A comparison of experimentally determined frequencies and symmetries of vibrational modes of the optical type (Raman and infrared) obtained in this and earlier studies with those calculated by us suggests that a consistent characterization of the vibrational properties was achieved. The good agreement between the experimental and simulated data on the temperature-dpendencies of the Raman modes (within 5%), crystal structure (2%), bulk modulus (5%), volume thermal expansivity (6%), and constant volume heat capacity (0.2%) testifies to the applicability of the transferable interatomic potential and the lattice dynamic model to predicting the vibrational, physical, and thermodynamic properties. The simulated properties from the lattice dynamic calculations are very similar to those obtained by molecular dynamic calculations with the same potential model.