Phase and group velocities and Q of mantle Love and Rayleigh waves of the first two modes and their azimuthal dependences for the 1963 Kurile Islands earthquake

Phase and group velocities and Q of mantle Love and Rayleigh waves from the 1963 Kurile Islands earthquake (M_w = 8.5) were determined over 37 great circle paths by a time variable filtering technique, in a period range 100–500 s for the fundamental modes and 100–275 s for the first higher modes. The preliminary reference Earth model (PREM) explains reasonably well the average dispersion results for the fundamental Love and Rayleigh waves. There exists a small, but significant inconsistency between the observation and the model for the first higher Love and Rayleigh waves. The Q structure of PREM is inconsistent with the observation for the fundamental Love waves, but explains other observations reasonably well. The dispersion of each mode shows a clear azimuthal dependence from which the four azimuthal windows were established. The phase and group velocity measurements for each window were, in general, shown to be mutually consistent. The azimuthal variations are largest for the first higher Rayleigh waves, indicating strong lateral heterogeneity in the structure of the low velocity zone. The first of the four windows is characterized by the largest fraction of Precambrian shields and the second window by the largest fraction of normal oceans. A comparison of these two windows may give some insight into deep lateral heterogeneity between continents and oceans. The observed phase and group velocities of the first window are systematically higher than those of the second window for the fundamental Love and Rayleigh waves at periods up to 400 s, and for the first higher Love and Rayleigh waves up to 175 s. Their differences are greatest for the first higher Rayleigh waves and least for the fundamental Rayleigh waves. Although the fundamental Rayleigh waves show the least velocity differences, their persistence up to a period of longer than 300 s is in striking contrast with some of the pure path phase velocities derived earlier for continents and oceans. A set of models for continents and oceans. PEM-C and PEM-O are not consistent with our observation. The third azimuthal window is characterized by trench-marginal seas and the fourth window by mountainous areas, typically the Asian high plateaus from northern China to the Middle East through Tibet. A comparison of these two windows gives some information about deep structural differences between subduction zones and continental collision zones, both belonging to plate convergence zones. For the fundamental and the first higher Love waves, the phase and group velocities for the third window are markedly low, whereas those for the fourth window are somewhat comparable to those for the second window. Slow Rayleigh waves are evident for two windows, with the fourth window apparently being the slowest for the fundamental Rayleigh above 200 s and for the first higher Rayleigh. For the fundamental Rayleigh waves, the third window is very slow below 200 s, but becomes progressively fast as the period increases and tends to be the fastest window around 400 s, suggesting a deep seated high velocity anomaly beneath trench-marginal seas. The dispersion characteristics of the fourth window indicate a thick high velocity lid with an extensive low velocity zone beneath it. The shield-like lithosphere, coupled with an extensive low velocity zone, may be a characteristic feature of continental collision zones. The particle motion of the fundamental Love waves was found not to be purely transverse to a great-circle connecting the epicenter to a station. The departure from the purely transverse motion is systematic among different periods, different G arrivals (G₂, G₃,…) and different stations, which may be interpreted as being due to lateral refraction.     

汉诺坝玄武岩中地幔岩捕掳体REE和Sr,Nd同位素地球化学

本文报道汉诺坝玄武岩中地幔岩捕掳体的ＲＥＥ丰度和Ｓｒ、Ｎｄ同位素组成。不同岩石类型的ＲＥＥ配分模式和同位素组成反映地幔部分熔融程度和交代作用过程。二辉橄榄岩亏损轻稀土,是原始地幔经不同程度部分熔融的残留体。方辉橄榄岩具Ｕ型ＲＥＥ配分模式,是强烈亏损的地幔岩被熔体非化学平衡交代的结果。二辉岩脉状体富轻、中稀土,它同与脉状体接触的二辉橄榄岩可达化学平衡或近于化学平衡,而二辉岩脉状体的形成与玄武岩岩浆无成因关系。据对二辉岩脉状体和不含脉状体橄榄岩的Ｓｍ－Ｎｄ同位素定年,这种脉状体形成于３００Ｍａ左右。     

A methodology for estimating business interruption loss caused by flood disasters: insights from business surveys after Tokai Heavy Rain in Japan

This study proposes a probabilistic methodology for estimating the business interruption loss of industrial sectors as an extension of current methodology. The functional forms and parameters are selected and calibrated based on survey data obtained from businesses located in the inundated area at the time of the 2000 Tokai Heavy Rain in Japan. The Tokai Heavy Rain was a rare event that hit a densely populated and industrialized area. In the estimation of business interruption losses, functional fragility curves and accelerated failure time models are selected to estimate the extent of damage to production capacity and production recovery time. Significant explanatory variables, such as inundation depth, distinct vulnerability, and the resilience characteristics of each sector, as well as the accuracy of fit of the model, are analyzed in the study. The function obtained and the estimated parameters can be utilized as benchmarks in estimating the probabilistic distribution of business interruption losses, especially in the case of urban flood disasters.

     

Excess bottom222Rn profiles and their implications in the northwestern Pacific Ocean

Vertical profiles of excess bottom²²²Rn and potential temperature were measured at 23 stations in the northwestern Pacific Ocean. The Rn profiles were classified into the following three types: quasi-exponential (type E), benthic boundary layer (type B), and horizontally disturbed (type H). The ratio among types E, B and H was approximately 2 : 3 : 1.An apparent vertical eddy diffusivity (K) was calculated by applying a one-dimensional diffusion model to the Rn profiles of types E and B. Type E had K values ranging from 15 to 180 cm² s^?1 (average: 70 cm² s^?1). As to type B, K values for the benthic boundary layer (4.5–260 cm² s^?1, average: 120 cm² s^?1) were always more than an order of magnitude larger than those for the upper layer (0.2–35 cm² s^?1, average: 7 cm² s^?1), indicating more active vertical mixing in the benthic boundary layer than in the upper layer.Rn profiles were measured in regions where the bottom topography is known. It was verified that the occurrence of type H related closely with local bottom topographic features accompanied by lateral transient supply of Rn-rich or Rn-poor water.A couple of Rn profiles at the same location, measured at time intervals of several years, were compared with each other for three locations. The general characteristics of Rn profiles were shown to remain unaltered with time, while the fine structure of Rn profiles may have short-term variations caused by local bottom topography and fluctuations of bottom current as indicated in type H.     

Water wedge advancing along the interface between two homogeneous layers

The behavior of the water wedge advancing along the interface between two homogeneous layers is investigated experimentally. The horizontally intruding wedge is strongly influenced by the detailed structure of the interface, even when the interface is sharp enough for its position to be visually determined. Within our experimental range, the intrusion velocity is very slow and depends on the thickness of the interface along which the water wedge intrudes. The wedge length is well described by the relation which is derived using an analogy to the linearly stratified case. The results could be useful for understanding the generation mechanism of the oceanic microstructure.     

GIS-based weights-of-evidence modelling of rainfall-induced landslides in small catchments for landslide susceptibility mapping

Landslide susceptibility mapping is a vital tool for disaster management and planning development activities in mountainous terrains of tropical and subtropical environments. In this paper, the weights-of-evidence modelling was applied, within a geographical information system (GIS), to derive landslide susceptibility map of two small catchments of Shikoku, Japan. The objective of this paper is to evaluate the importance of weights-of-evidence modelling in the generation of landslide susceptibility maps in relatively small catchments having an area less than 4 sq km. For the study area in Moriyuki and Monnyu catchments, northeast Shikoku Island in west Japan, a data set was generated at scale 1:5,000. Relevant thematic maps representing various factors (e.g. slope, aspect, relief, flow accumulation, soil depth, soil type, land use and distance to road) that are related to landslide activity were generated using field data and GIS techniques. Both catchments have homogeneous geology and only consist of Cretaceous granitic rock. Thus, bedrock geology was not considered in data layering during GIS analysis. Success rates were also estimated to evaluate the accuracy of landslide susceptibility maps and the weights-of-evidence modelling was found useful in landslide susceptibility mapping of small catchments.     

Temperature coefficients of elastic constants of single crystal MgO between 80 and 1,300 K

Elastic constants of single crystal MgO have been measured by the rectangular parallelepiped resonance (RPR) method at temperatures between 80 and 1,300 K. Elastic constants C _ij (Mbar=10³ kbar) and their temperature coefficients (kbar/K) are: $$\begin{gathered} {\text{ }}C_{{\text{11}}} {\text{ }}C_{{\text{12}}} {\text{ }}C_{{\text{44}}} {\text{ }}K_s {\text{ }}C_s \hfill \\ C_{ij} {\text{ 300 K 2}}{\text{.966 0}}{\text{.959 1}}{\text{.562 1}}{\text{.628 1}}{\text{.004}} \hfill \\ \partial C_{ij} {\text{/}}\partial T{\text{100 K }} - {\text{0}}{\text{.259 0}}{\text{.013 }} - {\text{0}}{\text{.072 }} - {\text{0}}{\text{.078 }} - {\text{0}}{\text{.136}} \hfill \\ {\text{ 300K }} - {\text{0}}{\text{.596 0}}{\text{.068 }} - {\text{0}}{\text{.122 }} - {\text{0}}{\text{.153 }} - {\text{0}}{\text{.332}} \hfill \\ {\text{ 800 K }} - {\text{0}}{\text{.619 0}}{\text{.009 }} - {\text{0}}{\text{.152 }} - {\text{0}}{\text{.200 }} - {\text{0}}{\text{.314}} \hfill \\ {\text{ 1,300 K }} - {\text{0}}{\text{.598 0}}{\text{.036 }} - {\text{0}}{\text{.130 }} - {\text{0}}{\text{.223 }} - {\text{0}}{\text{.218}} \hfill \\ \end{gathered} $$ By combining the present results with the previous data on the thermal expansivity and specific heat, the thermodynamic properties of magnesium oxide are presented and discussed. The elastic parameters of MgO at very high temperatures in the earth's lower mantle are also clarified.     

Emergent radiation from the atmosphere bounded by the two half Lambert surfaces composed of different albedo,II

For the evaluation of the effect of the non-uniform surface albedo on the emergent radiation from the atmosphere, the emergent radiation from the atmosphere bounded by the two half Lambert surfaces composed of different albedo is computed. This paper is the improved version of the previous paper (Takashima and Masuda, 1991). The atmosphere is assumed to be homogeneous, which is composed of aerosol, molecules, and absorbent gases. Their optical thicknesses are (1) 0.25, 0.23, and 0.02, and (2) 0.75, 0.23, and 0.02, respectively. The model aerosol is of the oceanic and water soluble types.In the computational procedure, the emergent radiation is calculated approximately by the contributions due to the multiple scattering in the atmosphere, and due to the diffusely or directly transmitted radiation through the atmosphere which is reflected by the surfaces once (4 interactive radiative modes between atmosphere and surface). Furthermore, to perform the hemispherical integration processing the radiative interaction, the transmission function based on the single scattering in the atmosphere is introduced and then the transmission function is averaged over the hemisphere with weighting function. The numerical simulation exhibits the extraordinary effect near the two half surface boundary of different albedoes. The effect decreases exponentially with the distance from the boundary. The effect depends on the atmospheric aerosol type, optical thickness, and surface albedo. The present version enables us to quantitatively discuss the radiative transfer trend near the boundary of two half surfaces. The upward radiance would simply be evaluated using the present scattering approximation method if the surface albedo is less than 0.3. The present method is thought of as a first step extending the one-dimensional radiative transfer model to two-dimensional using the doubling-adding method.     

Simulation analysis of ground liquefaction induced by earthquake

A dynamic effective stress analysis with the finite element method has long been recommended to predict the liquefaction phenomena of sandy soil by authors and Zienkiewicz et al. as well as the similar approaches by the others. Our approach of the analysis is summarized in the first.

Until recently, however, these approaches has not commonly been used as the means of design, although its capability of prediction is appreciated by geotechnicians. This method has been neglected because of the lack of verification studies of soil models of sand and mathematical formulation for boundary problems of liquefaction phenomena. Therefore the verification of the numerical method to evaluate liquefaction potential are urgent requirement of the recent engineering practice. To respond this requirement, extensive numerical studies on the liquefation simulations are performed by DIANA program for shaking table tests which have been conducted by the authors. The test models are soil-structure type models with combination of homogenous ground and partially improved ground by compaction.

The good performance of our approach is proved by the results of numerical simulation showing good agreement with experimental data in terms of response acceleration, excess pore pressure, and deformation profile. It is also demonstrated that the numerical results can provide substantial information to understand the mechanisms of soil ground behavior which is not easily obtained by experiments.

The procedure to identify soil constants for the reflecting surface model is also reported in details.