Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems

A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Multi-scale solution for building extraction from LiDAR and image data

This paper presents a multi-scale solution based on mathematical morphology for extracting the building features from remotely sensed elevation and spectral data. Elevation data are used as the primary data to delineate the structural information and are firstly represented on a morphological scale-space. The behaviors of elevation clusters across the scale-space are the cues for feature extraction. As a result, a complex structure can be extracted as a multi-part object in which each part is represented on a scale depending on its size. The building footprint is represented by the boundary of the largest part. Other object attributes include the area, height or number of stories. The spectral data is used as an additional source to remove vegetation and possibly classify the building roof material. Finally, the results can be stored in a multi-scale database introduced in this paper. The proposed solution is demonstrated using the data derived from a Light Detection And Ranging (LiDAR) surveying flight over Tokyo, Japan. The results show a reasonable match with reference data and prove the capability of the proposed approach in accommodation of diverse building shapes. Higher density LiDAR is expected to produce better accuracy in extraction, and more spectral sources are necessary for further classification of building roof material. It is also recommended that parallel processing should be implemented to reduce the computation time.     

Temporal and spatial variations in periglacial soil movements on alpine crest slopes

This paper describes up to ten years of continuous monitoring of frost heave, creep and associated parameters on high mountain crest slopes in the Japanese and Swiss Alps, aiming to evaluate spatial and interannual variations in the rates and controls of soil movement. Shallow frost creep re?ecting diurnal frost heave activity dominates the crest slopes that lack a vegetation mat and have a thin debris mantle with good drainage. Seasonal frost heave activity can induce slightly deeper movement where ?ne soil exists below the depth reached by diurnal freeze–thaw penetration, although the shallow bedrock impedes movements below 20 cm depth. As a result, downslope velocity pro?les display strong concavity with surface velocities of 2–50 cm a^?1. The frost creep rates vary spatially, depending on the soil texture, slope gradient, frequency of temperature cycling across 0 °C and moisture availability during freeze–thaw periods. Soil movements recur in every freeze–thaw period, although with some interannual variations affected by the length of seasonal snow cover and the occurrence of precipitation during freeze–thaw periods. The Swiss Alps encounter more signi?cant interannual variations than the Japanese Alps, re?ecting the large variability of the annual snow regime. Copyright © 2005 John Wiley & Sons, Ltd.     

Tectonic evolution along the northernmost border of the Philippine Sea plate since about 1 Ma

Recent structural, tephrochronologic and magnetostratigraphic studies conducted along the northernmost border of the Philippine Sea (PHS) plate enable us to reconstruct the precise tectonic evolution along the convergent boundary between the PHS plate and the Northeast Japan (NEJ) plate or the North American (NAM) plate since about 1 Ma. The authors of the present study split the tectonic evolution into five stages and present the characteristics of each stage. A plate tectonic interpretation is proposed, based upon the tectonic evolution, with special reference to the mode of convergent motion. In brief, our interpretations are as follows: the relative motion between the PHS and the NEJ plates was not recognized geologically within the area studied from about 1.0 to 0.9 Ma (Stage 1), suggesting either none or small influence from the coupling between the two plates during that period of time. Convergence between the PHS and the NEJ plates was possibly in N-S direction from 0.9 to 0.5 Ma (Stage II), and probably north-northwestward since 0.5 Ma (Stages III to V). The mode of the convergent motion was that of buoyant subduction in Stages II and III. The mode changed gradually from buoyant subduction during Stage IV to collision in Stage V (0.07 Ma to the present).     

Nankai Trough and Zenisu Ridge: a deep-sea submersible survey

Eight submersible dives between 3000 and 4200 m water depth were made off southern Japan in the eastern Nankai subduction zone. Benthic communities associated with chemosynthetic processes were discovered along the 800 m wide active tectonic zone, at the toe of the accretionary prism. A benthic community was also discovered along a zone of active compression, at the foot of Zenisu Ridge, 30 km south of Nankai Trough. Temperature measurements within the sediments below the benthic communities confirm that upward motion of interstitial water occurs there. Studies of water samples indicate advection of methane and light hydrocarbons. Specimens of the benthic community have been shown to have included in their shells carbonate resulting from methane consumption. Thus the benthic communities are related to overpressure-driven fluid advection along tectonic zones with active surface deformation. A 300 m high active scarp at the toe of the accretionary prism is related to relative motion in a 280° direction which is close to the 305° average direction of subduction in this area. The dives establish further that compressive deformation is presently occurring at the foot of Zenisu Ridge. The previous interpretation of the Zenisu Ridge as a zone of recent north-south intraplate shortening, 40 km south of the Nankai Trench, is confirmed. We conclude that tectonic evolution might well lead to future detachment of the Zenisu Ridge and overthrusting of this large piece of oceanic crust over the continental margin. Such a process might be an efficient one to emplace ophiolites over continents.     

Changes in tidal currents in the ariake sound due to reclamation

The Ariake Sound is a highly productive estuary located in the western part of Japan. The decline in fisheries and the frequent occurrence of red tides and hypoxia indicate a deterioration in the ecosystem of the sound. A change in tidal currents, which is one of the possible causes of the deterioration, was investigated by numerical experiments. Two major changes in the topography of the sound, which may have changed the tidal currents, are examined. These are the reclamation in the innermost part of the sound and the construction of a dyke in a subembayment called Isahaya Bay. The numerical experiments show that the reclamation caused the tidal currents to decrease by more than 10% over a large area in the innermost part of the sound. The influence of the dyke on the tidal currents is relatively local compared to the reclamation; the area where tidal currents significantly decrease is located mostly in the subembayment.     

Time variation of the X-ray spectrum and optical luminosity of SCO X-1

Results of rocket observations of SCO X-1 over the spectral range of 220 keV are presented. The observations have been performed partly in India and partly in Japan under the collaboration of the three groups. The present results are compared with results of similar observations carried out by the LRL (Lawrence Radiation Laboratory) group. Some of these X-ray observations were accompanied by simultaneous optical observations. Relationships between the hardness of the X-ray spectrum and the X-ray intensity and between the hardness and the optical luminosity are compiled. The relationships among the parameters (temperature, density and size) which characterize the postulated isothermal cloud model of SCO X-1 are given. They indicate that SCO X-1 is characterized by a temperature of about 10⁷–10⁸K, a density of about 10¹⁶–10¹⁷ cm^–3 and a radius of about 10⁸–10⁹ cm respectively. We further show that the temperature is inversely correlated with the size of the source; an increase in temperature corresponds to a decrease in the radius and an increase in density.     

Time variations of hard X-rays from Sco X-1

Simultaneous hard X-ray and optical observations of Sco X-1 were carried out on 1971 May 1 at Hyderabad, India, when Sco X-1 was optically bright. The X-ray intensity observed by balloon-borne counter telescopes increased in coincidence with optical enhancements, while the plasma temperature derived by fitting the X-ray spectrum in the energy range 20–40 keV to the thermal bremsstrahlung spectrum did not appreciably change over the whole period of observation.     

Preliminary Examinations of the Ore-forming Process by Fluid Mixing – a Test of MIX99

Abstract: Possible ore deposition by fluid mixing was preliminarily examined using MIX 99, a FORTRAN program developed for simulating water-rock interactions. We consider mixing of two fluids, the low f_O2 and high temperature source fluid and the high f_O2 and low temperature seawater. Oxygen fugacity of a mixed fluid formed by titration of seawater into the source fluid gradually decreases with decreasing temperature (model A). Sequential precipitation of ore-forming minerals was examined in this model. On the other hand, simultaneous precipitation of the minerals could be followed by simulation of instantaneous overall mixing of the two fluids (model B). Results of simulation of the both models revealed that a temporal sequence of mineralization observed in the Karuizawa mine, NE Japan, can be formed by model A, while model B is suitable for the mineralization of an active chimney found in the Rainbow hydrothermal area of the Mid-Atlantic Ridge.