Configuration of subducting Philippine Sea plate and crustal structure in the central Japan region

A seismic experiment with six explosive sources and 391 seismic stations was conducted in August 2001 in the central Japan region. The crustal velocity structure for the central part of Japan and configuration of the subducting Philippine Sea plate were revealed. A large lateral variation of the thickness of the sedimentary layer was observed, and the P-wave velocity values below the sedimentary layer obtained were 5.3–5.8 km/s. P-wave velocity values for the lower part of upper crust and lower crust were estimated to be 6.0–6.4 and 6.6–6.8 km/s, respectively. The reflected wave from the upper boundary of the subducting Philippine Sea plate was observed on the record sections of several shots. The configuration of the subducting Philippine Sea slab was revealed for depths of 20–35 km. The dip angle of the Philippine Sea plate was estimated to be 26° for a depth range of about 20–26 km. Below this depth, the upper boundary of the subducting Philippine Sea plate is distorted over a depth range of 26–33 km. A large variation of the reflected-wave amplitude with depth along the subducting plate was observed. At a depth of about 20–26 km, the amplitude of the reflected wave is not large, and is explained by the reflected wave at the upper boundary of the subducting oceanic crust. However, the reflected wave from reflection points deeper than 26 km showed a large amplitude that cannot be explained by several reliable velocity models. Some unique seismic structures have to be considered to explain the observed data. Such unique structures will provide important information to know the mechanism of inter-plate earthquakes.     

Numerical Simulation of Displacement and Stress Fields Associated with the 1993 Kushiro-oki, Japan, Earthquake

—We constructed a three-dimensional finite element model to simulate coseismic and postseismic displacement and stress fields associated with the 1993 Kushiro-oki earthquake, which was a very large intermediate-depth earthquake that occurred within the subducted Pacific plate at a depth of 107 km beneath the southeastern part of Hokkaido, Japan. Taking the configuration of the subducted Pacific plate into account, we constructed a realistic model with lateral heterogeneity of viscoelastic structure. We assigned a variable slip distribution to the fault plane, which was obtained from inversion analysis of near-field seismic waveforms. The result shows that elastic deformation associated with the faulting reflects the assigned inhomogeneous slip distribution on the fault plane near the fault region, while it does not reflect the distribution on the free surface of the model. The calculated postseismic deformation does not reflect the slip distribution, but shows symmetric spatial patterns concerning the dipping direction of the fault both near the fault region and on the model surface. For the next 20 years following the earthquake, the amount of the calculated deformation is a fraction of the coseismic deformation. The calculated coeseismic deformation is large just above and below the fault plane, reaching 1 m, while the postseismic deformation is dominant near the upper and lower material boundaries between the subducted plate and the surrounding asthenosphere. The spatial distribution of maximum shear stress near the fault plane corresponds to the assigned slip distribution, amounting to 32 MPa. The directions of principal stress-change axes represent reverse fault type in the SSE region of the fault, whereas normal fault type is dominant in the NNW region with the exception of some asymmetrical spatial patterns of the principal stress-change axes on the fault due to the inhomogeneous slip distribution. Time variations both in the amount and the directions of stresses are minor, suggesting that the coseismic state of the stress would remain unchanged for two decades after the event.     

Northern Hemisphere atmospheric blocking as simulated by 15 atmospheric general circulation models in the period 1979–1988

 As a part of the Atmospheric Model Intercomparison Project (AMIP), the behaviour of 15 general circulation models has been analysed in order to diagnose and compare the ability of the different models in simulating Northern Hemisphere midlatitude atmospheric blocking. In accordance with the established AMIP procedure, the 10-year model integrations were performed using prescribed, time-evolving monthly mean observed SSTs spanning the period January 1979–December 1988. Atmospheric observational data (ECMWF analyses) over the same period have been also used to verify the models results. The models involved in this comparison represent a wide spectrum of model complexity, with different horizontal and vertical resolution, numerical techniques and physical parametrizations, and exhibit large differences in blocking behaviour. Nevertheless, a few common features can be found, such as the general tendency to underestimate both blocking frequency and the average duration of blocks. The problem of the possible relationship between model blocking and model systematic errors has also been assessed, although without resorting to ad-hoc numerical experimentation it is impossible to relate with certainty particular model deficiencies in representing blocking to precise parts of the model formulation. Received: 16 July 1997/Accepted: 20 October 1997     

Stick–slip behavior of Indian gabbro as studied using a NIED large-scale biaxial friction apparatus

This paper reports stick–slip behaviors of Indian gabbro as studied using a new large-scale biaxial friction apparatus, built in the National Research Institute for Earth Science and Disaster Prevention(NIED), Tsukuba, Japan. The apparatus consists of the existing shaking table as the shear-loading device up to 3,600 k N, the main frame for holding two large rectangular prismatic specimens with a sliding area of 0.75 m2 and for applying normal stresses rnup to 1.33 MPa, and a reaction force unit holding the stationary specimen to the ground. The shaking table can produce loading rates v up to 1.0 m/s,accelerations up to 9.4 m/s2, and displacements d up to0.44 m, using four servocontrolled actuators. We report results from eight preliminary experiments conducted with room humidity on the same gabbro specimens at v = 0.1–100 mm/s and rn= 0.66–1.33 MPa, and with d of about 0.39 m. The peak and steady-state friction coefficients were about 0.8 and 0.6, respectively, consistent with the Byerlee friction. The axial force drop or shearstress drop during an abrupt slip is linearly proportional to the amount of displacement, and the slope of this relationship determines the stiffness of the apparatus as1.15 9 108N/m or 153 MPa/m for the specimens we used.This low stiffness makes fault motion very unstable and the overshooting of shear stress to a negative value was recognized in some violent stick–slip events. An abrupt slip occurred in a constant rise time of 16–18 ms despite wide variation of the stress drop, and an average velocity during an abrupt slip is linearly proportional to the stress drop.The use of a large-scale shaking table has a great potential in increasing the slip rate and total displacement in biaxial friction experiments with large specimens.     

Modelling basin-scale distribution of fish occurrence probability for assessment of flow and habitat conditions in rivers

Abstract

Flow regimes play an important role in sustaining biodiversity in river ecosystems. However, the effects of flow regimes on riverine fish have not been clearly described. Therefore, we propose a new methodology to quantitatively link habitat conditions (such as flow indices and physical habitat conditions) to the occurrence probability (OP) of fish species. We developed a basin-scale fish distribution model by integrating the concept of habitat suitability assessment with a distributed hydrological model in order to estimate the OP of fish, with particular attention to flow regime. A generalized linear model was used to evaluate the relationship between the probabilities of fish occurrence and major environmental factors in river sections. A geomorphology-based hydrological model was adopted to simulate river discharge, which was used to calculate 10 flow indices. The occurrence probabilities of 50 fish species in the Sagami River in Japan were modelled. For the prediction accuracy, field survey results that included at least five observations of both the presence and the absence of each species were required to obtain relatively reliable prediction (accuracy > 60%). Using the developed model, important habitat conditions for each species were identified, which showed the importance of low-flow events for more than 10 species, including Hypomesus nipponensis and Rhinogobius fluviatilis. The model also confirmed the positive effects of natural flow and the negative effect of river-crossing structures, such as dams and weirs, on the OP of most species. The suggested approach enables us to evaluate and project the ecological consequences of water resource management policy. The results demonstrate the applicability of the fish distribution model to provide quantitative information on the flow required to maintain fish communities.

Editor Z.W. Kundzewicz; Guest editor M. Acreman

Citation Sui, P., Iwasaki, A., Saavedra, V.O.C., and Yoshimura, C., 2013. Modelling basin-scale distribution of fish occurrence probability for assessment of flow and habitat conditions in rivers. Hydrological Sciences Journal, 59 (3–4), 618–628.     

Self-similar solutions for the dynamical condensation of a radiative gas layer

A new self-similar solution describing the dynamical condensation of a radiative gas is investigated under a plane-parallel geometry. The dynamical condensation is caused by thermal instability. The solution is applicable to generic flow with a net cooling rate per unit volume and time  ∝ρ² T ^α  , where  ρ,  T   and α are the density, temperature and a free parameter, respectively. Given α, a family of self-similar solutions with one parameter η is found in which the central density and pressure evolve as follows:  ρ( x = 0, t ) ∝ ( t _c− t )^{−η/(2−α)}  and   P ( x = 0, t ) ∝ ( t _c− t )^{(1−η)/(1−α)}  , where t _c is the epoch at which the central density becomes infinite. For  η∼ 0  the solution describes the isochoric mode, whereas for  η∼ 1  the solution describes the isobaric mode. The self-similar solutions exist in the range between the two limits; that is, for  0 < η < 1  . No self-similar solution is found for  α > 1  . We compare the obtained self-similar solutions with the results of one-dimensional hydrodynamical simulations. In a converging flow, the results of the numerical simulations agree well with the self-similar solutions in the high-density limit. Our self-similar solutions are applicable to the formation of interstellar clouds (H  i clouds and molecular clouds) by thermal instability.     

Diffuse component of the cosmic far UV radiation and interstellar dust grains

The diffuse far UV radiation ( 1350–1480 Å) observed in the sky region ofl ^II180°, 0°b ^II40° is analyzed in connection with the distributions of stars and dust grains as well as with optical properties of grains. Its intensity (starlight+scattered light) is about 6×10^–7 erg cm^–2 sec^–1 sr^–1 Å^–1 in the direction ofb ^II0° andl ^II180°. The latitude dependence of the intensity is in approximate agreement with the plane parallel slab model of the galaxy with a reasonable set of parameters. The interstellar scattering gives an albedo close to unity and forward phase function of about 0.6, which are not inconsistent with the model of interstellar grains of Wickramasinghe. The upper limit of the extragalactic UV is 2×10^–8 erg cm^–2 sec^–1 sr^–1 Å^–1 in the same region of wave-length.     

Assessment of factors influencing groundwater-level change using groundwater flow simulation,considering vertical infiltration from rice-planted and crop-rotated paddy fields in Japan

Assessing factors that influence groundwater levels such as land use and pumping strategy, is essential to adequately manage groundwater resources. A transient numerical model for groundwater flow with infiltration was developed for the Tedori River alluvial fan (140 km²), Japan. The main water input into the groundwater body in this area is irrigation water, which is significantly influenced by land use, namely paddy and upland fields. The proposed model consists of two models, a one-dimensional (1-D) unsaturated-zone water flow model (HYDRUS-1D) for estimating groundwater recharge and a 3-D groundwater flow model (MODFLOW). Numerical simulation of groundwater flow from October 1975 to November 2009 was performed to validate the model. Simulation revealed seasonal groundwater level fluctuations, affected by paddy irrigation management. However, computational accuracy was limited by the spatiotemporal data resolution of the groundwater use. Both annual groundwater levels and recharge during the irrigation periods from 1975 to 2009 showed long-term decreasing trends. With the decline in rice-planted paddy field area, groundwater recharge cumulatively decreased to 61 % of the peak in 1977. A paddy-upland crop-rotation system could decrease groundwater recharge to 73–98 % relative to no crop rotation.