Kalman filter finite element method applied to dynamic ground motion

The purpose of this paper is to investigate the estimation of dynamic elastic behavior of the ground using the Kalman filter finite element method. In the present paper, as the state equation, the balance of stress equation, the strain–displacement equation and the stress–strain equation are used. For temporal discretization, the Newmark ¼ method is employed, and for the spatial discretization the Galerkin method is applied. The Kalman filter finite element method is a combination of the Kalman filter and the finite element method. The present method is adaptable to estimations not only in time but also in space, as we have confirmed by its application to the Futatsuishi quarry site. The input data are the measured velocity, acceleration, etc., which may include mechanical noise. It has been shown in numerical studies that the estimated velocity, acceleration, etc., at any other spatial and temporal point can be obtained by removing the noise included in the observation. Copyright © 2008 John Wiley & Sons, Ltd.     

Hortonian overland flow from Japanese forest plantations—an aberration,the real thing,or something in between?

There is a growing opinion that poorly managed plantation forests in Japan are contributing to increased storm runoff and erosion. Here we present evidence to the contrary from runoff plots at two scales (hillslope and 0·5 × 2 m plots) for several forest conditions in the Mie and Nariki catchments. Runoff coefficients from small plots in untended hinoki forests were variable but typically higher than from better managed or deciduous forests during small storms at Nariki; at Mie, runoff during small events was highly variable from all small plots but runoff coefficients were similar for hinoki plots with and without understory vegetation, while the deciduous plot had lower runoff coefficients. Storm runoff was less at the hillslope scale than the plot scale in Mie; these results were more evident at sites with better ground cover. During the largest storms at both sites, differences in runoff due to forest condition were not evident regardless of scale. Dynamic soil moisture tension measurements at Nariki indicated that during a large storm, flow in the upper organic‐rich and root‐permeated soil horizons was 3·2 times higher than measured overland runoff from a small hinoki plot with poor ground cover and 8·3 times higher than runoff from a deciduous forest plot. On the basis of field observations during storms, at least a portion of the monitored ‘Hortonian overland flow’ was actually occurring in this near‐surface ‘biomat’. Therefore our field measurements in both small and large plots potentially included biomat flow in addition to short‐lived Hortonian runoff. Because overland flow decreased with increasing scale, rill erosion did not occur on hillslopes. Additionally, runoff coefficients were not significantly different among cover conditions during large storms; thus, the ‘degraded’ forest conditions appear not to greatly enhance peak flows or erosion potential at larger scales, especially when biomat flow is significant. Copyright © 2007 John Wiley & Sons, Ltd.     

Electrokinetic phenomena associated with a water injection experiment at the Nojima fault on Awaji Island, Japan

Abstract Self-potential variations were measured to estimate the magnitude of electrokinetic and hydrological parameters (zeta potential and permeability) of the Nojima Fault zone in Awaji, Japan. The study observed self-potential variations that seemed to be associated with water flow from the injection well to the fracture zone, which were induced by turning the injection on and off. Amplitudes of the variations were a few to 0.03 V across 320–450 m dipoles. These variations can be explained well with an electrokinetic model. The quantity k/ζ (permeability/zeta potential) is in the range 1.6 × 10⁻¹³− 5.4 × 10⁻¹³ m²/V. Permeability of the Nojima fault zone can be estimated as approximately 10⁻¹⁶–10⁻¹⁵ m² on the assumption that the zeta potential is in the range –0.01 to –0.001 V.     

Synthesis of γ-Mg2SiO4

Magnesium orthosilicate with spinel structure (γ-Mg₂SiO₄) was synthesized at about 250 kbar and 1000°C. Unit cell dimension was established to be 8.076 ± 0.001Å. X-ray powder diffraction pattern revealed a significant difference between γ-Mg₂SiO₄ and other γ-M₂SiO₄ spinels (M = Fe, Co, and Ni) in the intensities of (111) and (331) reflections, both of which are virtually absent in the Mg₂SiO₄ spinel. This feature could be thoroughly understood by the calculation of the intensities for several silicate spinels.     

Young,olivine xenocryst-bearing alkali-basalt from the oceanward slope of the Japan Trench

Young (6 Ma) alkali-basalts were collected from the toe of the oceanward slope of the northern Japan Trench. Two types of olivine are present in these lavas, xenocrysts with reaction rims and magmatic. The forsterite (Fo) (no. 91–92) values and NiO contents (0.3–0.5 wt%) of the xenocrysts are similar in composition to those of the depleted mantle peridotite. The groundmass olivines have relatively lower Fo values (no. 81–88) and NiO contents (0.1–0.5 wt%). Reaction rims and the vicinity of the silicate inclusion in xenocrysts show the intermediate compositions between the xenocryst and magmatic olivines. Chromian spinel inclusions in the xenocrysts also show the depleted composition in the range of abyssal peridotite. CO₂ fluid inclusions in the xenocryst records pressures before entrainment into the host magma up to 0.4 GPa, which corresponds to a depth of up to 14 km of lithospheric mantle. These data indicate that the xenocrysts originate from MORB-depleted mantle.     

Geographical Distribution of 3He/4He Ratios in the Chugoku District, Southwestern Japan

We have collected 34 hot spring and mineral spring gases and waters in the Chugoku and Kansai districts, Southwestern Japan and measured the ³He/⁴He and ⁴He/²⁰Ne ratios by using a noble gas mass spectrometer. Observed ³He/⁴He and ⁴He/²⁰Ne ratios range from 0.054 R_atm to 5.04 R_atm (where R_atm is the atmospheric ³He/⁴He ratio of 1.39 × 10⁻⁶) and from 0.25 to 36.8, respectively. They are well explained by a mixing of three components, mantle-derived, radiogenic, and atmospheric helium dissolved in water. The ³He/⁴He ratios corrected for air contamination are low in the frontal arc and high in the volcanic arc regions, which are consistent with data of subduction zones in the literature. The geographical contrast may provide a constraint on the position of the volcanic front in the Chugoku district where it was not well defined by previous works. Taking into account the magma aging effect, we cannot explain the high ³He/⁴He ratios of the volcanic arc region by the slab melting of the subducting Philippine Sea plate. The other source with pristine mantle material may be required. More precisely, the highest and average ³He/⁴He ratios of 5.88 R_atm and 3.8±1.6 R_atm, respectively, in the narrow regions near the volcanic front of the Chugoku district are lower than those in Kyushu and Kinki Spot in Southwestern Japan, but close to those in NE Japan. This suggests that the magma source of the former may be related to the subduction of the Pacific plate, in addition to a slight component of melting of the Philippine Sea slab.     

Identification of critical ground motions for seismic performance assessment of structures

A method is established to identify critical earthquake ground motions that are to be used in physical testing or subsequent advanced computational studies to enable seismic performance to be assessed. The ground motion identification procedure consists of: choosing a suitable suite of ground motions and an appropriate intensity measure; selecting a computational tool and modelling the structure accordingly; performing Incremental Dynamic Analysis on a non‐linear model of the structure; interpreting these results into 50th (median) and 90th percentile performance bounds; and identifying the critical ground motions that are close to these defining probabilistic curves at ground motion intensities corresponding to the design basis earthquake and the maximum considered earthquake. An illustrative example of the procedure is given for a reinforced concrete highway bridge pier designed to New Zealand specifications. Pseudodynamic tests and finite element based time history analyses are performed on the pier using three earthquake ground motions identified as: (i) a Design Basis Earthquake (10% probability in 50 years) with 90 percent confidence of non‐exceedance; (ii) a Maximum Considered Event (2% probability in 50 years) representing a median response; and (iii) a Maximum Considered Event representing 90 percent confidence of non‐exceedance. Copyright © 2006 John Wiley & Sons, Ltd.     

Kalman filtering in groundwater flow modelling: problems and prospects

The popularity of applying filtering theory in the environmental and hydrological sciences passed its first climax in the 1970s. Like so many other new mathematical methods it was simply the fashion at the time. The study of groundwater systems was not immune to this fashion, but neither was it by any means a prominent area of application. The spatial-temporal characteristics of groundwater flow are customarily described by analytical or, more frequently, numerical, physics-based models. Consequently, the state-space representations associated with filtering must be of a high order, with an immediately apparent computational over-burden. And therein lies part of the reason for the but modest interest there has been in applying Kalman filtering to groundwater systems, as reviewed critically in this paper. Filtering theory may be used to address a variety of problems, such as: state estimation and reconstruction, parameter estimation (including the study of uncertainty and its propagation), combined state-parameter estimation, input estimation, estimation of the variance-covariance properties of stochastic disturbances, the design of observation networks, and the analysis of parameter identifiability. A large proportion of previous studies has dealt with the problem of parameter estimation in one form or another. This may well not remain the focus of attention in the future. Instead, filtering theory may find wider application in the context of data assimilation, that is, in reconstructing fields of flow and the migration of sub-surface contaminant plumes from relatively sparse observations. Received: October 27, 1997     

Fine structure of laboratory wind-wave surfaces studied using an optical method

The fine structure of laboratory wind-wave surfaces was investigated using an optical method. Several characteristic structures, closely associated with wind speed and the stage of development of the waves, are described. A parallel light, incident on the wind-wave surface, was backscattered by specular facets perpendicular to the incident angle. These specular facets, which represented the fine structures, were photographed and quantitatively analyzed by image processing techniques.For wind speeds ranging from 3.6 to 13.6 m s^–1, different types of small-scale undulations appeared: a three-dimensional rhombic structure; a train of capillary waves on the forward face with wavelengths gradually decreasing with distance from the crest, and a streaky structure on the backward face in the direction of the wind; and random and wide-spreading irregularities covering the whole wind-wave surface. The spatial scales and temporal variation of the undulations are described.