Effects of Three-dimensional Inhomogeneous Viscoelastic Structures on Postseismic Surface Deformations Associated with the Great 1946 Nankaido Earthquake

—We investigated the effects of various viscoelastic structures on postseismic surface displacement and principal strain fields associated with the great 1946 Nankaido earthquake, which occurred on the plate boundary between the subducting Philippine Sea plate and the continental Eurasian plate. For this purpose, we constructed two kinds of three-dimensional structural models using the finite element method one is the Layered Model, in which a semi-infinite Maxwell viscoelastic material is underlying an elastic layer, and the other is the more realistic Plate Model, in which the three-dimensional configuration of the subducted Philippine Sea plate is taken into account. We also considered two cases with different thicknesses of the elastic layer (50 and 33km) for the respective models. The difference between the two models in postseismic surface deformations is significant for the case with the thinner elastic layer. In this case the horizontal surface displacement and principal strain for the Layered Model is two to three times larger than those for the Plate Model. Downward surface deformation tends to be dominant for the Layered Model, while the change in the pattern for the Plate Model is less marked. The spatial extent of uplift and subsidence for the Plate Model is broader than that for the Layered Model. Postseismic vertical displacements in Shikoku were found to be strongly dependent on the viscoelastic structures. From these results, we suggest that the estimates of the viscosity of the uppermost mantle, interplate coupling, and the area and the amount of after-slip following the 1946 Nankaido earthquake, which have been estimated based on simple layered viscoelastic models, should be re-evaluated using realistic three-dimensionally heterogeneous viscoelastic structures.     

Bathymetric and geochemical investigation of Kawah Ijen Crater Lake, East Java, Indonesia

A bathymetric survey of Kawah Ijen crater lake was conducted by acoustic sounding in 1996 to compare the lake morphology with those measured in 1922, 1925 and 1938, and to calculate the present lake volume. Even though the lake experienced several hydrothermal eruptions, the maximum depth became shallower (182 m) than before (200 m), resulting in a reduced lake volume (3.0×10⁷ m³).Fifty-two major and minor constituents including rare earth elements and polythionates (PT) of the lake waters at various depths were determined by ICP-AES, ICP-MS and HPLC, respectively. These ions except for several volatile elements are taken up by lake fringe through congruent dissolution of pyroclastics of Kawah Ijen volcano. Most ions are homogeneously distributed throughout the lake, although PT showed a considerable vertical variation. Rare earth elements (REE) in the Kawah Ijen water as well as those from other hyper-acidic crater lakes show distribution patterns likely due to the three rock dissolution (preferential, congruent and residual) types, and their logarithmic concentrations linearly depend upon the pH values of the lake waters.Using the PT degradation kinetics data, production rates of PT, injection rates of SO₂ and H₂S into the lake were estimated to be 114, 86 and 30 tons/day, respectively. Also travel time of the spring water at the Banyupahit Riverhead from Kawah Ijen was estimated to be 600–1000 days through the consideration of decreasing rates of PT. Molten sulfur stocks containing Sn, Cu, Bi sulfides and Pb-barite exposed on the inner crater slope were presumed to be extinct molten sulfur pools at the former lake bottom. This was strongly supported by the barite precipitation temperature estimated through the consideration of the temperature dependence of Pb-chlorocomplex formation.     

Comparison of runoff characteristics of two adjacent basins in a tropical rainforest using a modified hydrologic cycle model with outflow

We propose a new runoff model including an outflow process that was applied to two adjacent basins (CL, TL) located in Lambir Hills National Park in north‐central Sarawak, Malaysia. Rainfall, runoff, topography, and soil layer thickness were observed. About 19% of annual runoff was observed in the CL basin (21.97 ha), whereas about 46% was observed in the TL basin (23.25 ha). It was inferred that the CL basin has an outflow because of low base flow, small runoff peak, and excessive water loss. By incorporating the outflow process into the HYdrological CYcle MODEL, good agreement between the data generated by the model and that observed was shown, with the exception of the data from the rainless period. Then, the fitting parameters for each basin were exchanged, except for the outflow parameter, and the characteristics of each basin were compared by calculating virtual runoff. As a result, the low base flow of the CL basin was estimated by the movement of the rainwater that escaped from the basin as deep percolation or lateral flow (11% of rainfall). The potential of the CL basin for mitigating flood and drought appeared to be higher than that of the TL basin. This is consistent with the topographic characteristics of the CL basin, which has a gentler slope than the TL basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Viscosity of peridotite liquid up to 13 GPa: Implications for magma ocean viscosities

The viscosity of synthetic peridotite liquid has been investigated at high pressures using in-situ falling sphere viscometry by combining a multi-anvil technique with synchrotron radiation. We used a newly designed capsule containing a small recessed reservoir outside of the hot spot of the heater, in which a viscosity marker sphere is embedded in a forsterite + enstatite mixture having a higher solidus temperature than the peridotite. This experimental setup prevents spheres from falling before a stable temperature above the liquidus is established and thus avoids difficulties in evaluating viscosities from velocities of spheres falling through a partially molten sample.

Experiments have been performed between 2.8 and 13 GPa at temperatures ranging from 2043 to 2523 K. Measured viscosities range from 0.019 (± 0.004) to 0.13 (± 0.02) Pa s. At constant temperature, viscosity increases with increasing pressure up to  8.5 GPa but then decreases between  8.5 and 13 GPa. The change in the pressure dependence of viscosity is likely associated with structural changes of the liquid that occur upon compression. By combining our results with recently published 0.1 MPa peridotite liquid viscosities [D.B. Dingwell, C. Courtial, D. Giordano, A. Nichols, Viscosity of peridotite liquid, Earth Planet. Sci. Lett. 226 (2004) 127–138.], the experimental data can be described by a non-Arrhenian, empirical Vogel-Fulcher-Tamman equation, which has been modified by adding a term to account for the observed pressure dependence of viscosity. This equation reproduces measured viscosities to within 0.08 log₁₀-units on average. We use this model to calculate viscosities of a peridotitic magma ocean along a liquid adiabat to a depth of  400 km and discuss possible effects on viscosity at greater pressures and temperatures than experimentally investigated.     

Effects of inertial and kinematic interaction on seismic behavior of pile with embedded foundation

Effects of inertial and kinematic forces on pile stresses are studied based on large shaking table tests on pile-structure models with a foundation embedded in dry and liquefiable sand deposits. The test results show that, if the natural period of the superstructure, T_b, is less than that of the ground, T_g, the ground displacement tends to be in phase with the inertial force from the superstructure, increasing the shear force transmitted to the pile. In contrast, if T_b is greater than T_g, the ground displacement tends to be out of phase with the inertial force, restraining the pile stress from increasing. With the effects of earth pressures on the embedded foundation and pile incorporated in, pseudo-static analysis is conducted to estimate maximum moment distribution in pile. It is assumed that the maximum moment is equal to the sum of the two stresses caused by the inertial and kinematic effects if T_b<T_g or the square root of the sum of the squares of the two if T_b>T_g. The estimated pile stresses are in good agreement with the observed ones regardless of the occurrence of soil liquefaction.     

Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

A Late Cretaceous mammalian dentary from the Ashizawa Formation (Futaba Group), Fukushima,northeastern Japan

A mammalian dentary discovered in the Coniacian Ashizawa Formation (Fukushima, northeastern Japan) is described. The specimen is a fragment of the horizontal ramus of a left edentulous dentary with five alveoli, the distal four of which are plugged with broken roots. Based on the morphologies of the dentary and the roots, it is considered to be of a therian mammal. This constitutes the first discovery of a Mesozoic mammal in northeastern Japan and highlights the potential for future mammal discoveries in the Cretaceous System in northeastern Japan, which will be significant for disclosure of the mammalian faunal evolution in East Asia during the Late Cretaceous.     

Chemical characteristics of chromian spinel in plutonic rocks: Implications for deep magma processes and discrimination of tectonic setting

We summarize chemical characteristics of chromian spinels from ultramafic to mafic plutonic rocks (lherzolites, harzburgites, dunites, wehrlites, troctolites, olivine gabbros) with regard to three tectonic settings (mid‐ocean ridge, arc, oceanic hotspot). The chemical range of spinels is distinguishable between the three settings in terms of Cr# (= Cr/(Cr + Al) atomic ratio) and Ti content. The relationships are almost parallel with those of chromian spinels in volcanic rocks, but the Ti content is slightly lower in plutonics than in volcanics at a given tectonic environment. The Cr# of spinels in plutonic rocks is highly diverse; its ranges overlap between the three settings, but extend to higher values (up to 0.8) in arc and oceanic hotspot environments. The Ti content of spinels in plutonics increases, for a given lithology, from the arc to oceanic hotspot settings by mid‐ocean ridge on average. This chemical diversity is consistent with that of erupted magmas from the three settings. If we systematically know the chemistry of chromian spinels from a series of plutonic rocks, we can estimate their tectonic environments of formation. The spinel chemistry is especially useful in dunitic rocks, in which chromian spinel is the only discriminating mineral. Applying this, discordant dunites cutting mantle peridotites were possibly precipitated from arc‐related magmas in the Oman ophiolite, and from an intraplate tholeiite in the Lizard ophiolite, Cornwall.     

Physical parameters identification for full-scale ten-story reinforced concrete building with degrading tri-linear model by modal iterative error correction method

To assess the post-earthquake seismic safety of buildings, it is crucial to predict seismic response, and it is necessary to set the appropriate physical parameters of the response analysis model. Numerous methods have been proposed to identify physical parameters. However, most of them are limited to linear systems, and previous researches on nonlinear systems have difficulties in practical applications. In this paper, a nonlinear response analysis model is identified for a full-scale ten-story reinforced concrete building with the degrading tri-linear stiffness model by the modal iterative error correction (MIEC) method, and the accuracy of this technique is discussed by comparing with the shaking table test.     

Three-dimensional EM computer simulation on sprite initiation above a horizontal lightning discharge

The purpose of this paper is to study the three-dimensional (3D) effects of the source current and the electromagnetic (EM) pulse on the distribution of upper atmospheric electric field and ionization created. A new lightning model has been employed; i.e., we have included a horizontal channel in addition to the conventional vertical channel. The effects of the horizontal lightning channel are summarized as follows: (1) the effect of a vertical channel plays a fundamental role in the whole view of sprite initiation; (2) the position of a sprite is shifted from the position of its parent vertical channel in response to the length of a horizontal channel; and (3) we observe very fine structures with local maxima and minima in the “reduced” electric field (the electric field divided by the neutral gas density). These theoretical inferences are discussed and compared with the observational facts so far reported (such as lateral shift of sprites, morphological difference of sprites, etc.).