A boundary element method for identifying orthotropic material parameters

A new approach is proposed for identifying elastic constants for orthotropic material bodies by using the boundary element method. The material parameters which must be determined are two Young's moduli, a shear modulus and two Poisson's ratios. The method is applied to geotechnical engineering problems.     

Crystal chemistry of wadsleyite II and water in the Earth’s interior

Wadsleyite II is a variably hydrous magnesium-iron silicate phase similar to spinelloid IV and a potential host for H in the Transition Zone of the Earths mantle. Two separate samples of wadsleyite II synthesized at 17.5 GPa and 1400°C and at 18 GPa and 1350°C have been characterized by electron microprobe, single-crystal X-ray diffraction, visible, IR, Raman, and Mössbauer spectroscopies, and transmission electron microscopy including electron energy-loss spectroscopy. The two samples have the following chemical formulae: Mg_1.71Fe_0.18Al_0.01H_0.33Si_0.96O₄ and Mg_1.60Fe_0.22Al_0.01 H_0.44Si_0.97O₄. Mössbauer spectroscopy and electron energy loss spectroscopy (EELS) indicate that about half of the iron present is ferric. Refinement of the structures shows them to be essentially the same as spinelloid IV. Calculated X-ray powder diffraction patterns show only subtle differences between wadsleyite and wadsleyite II. The hydration mechanism appears to be protonation of the non-silicate oxygen (O2) and possibly the oxygens surrounding the partially vacant tetrahedral site Si2, charge-balanced by cation vacancies in Si2, M5 and M6. The unit cell volume of this phase and its synthesis conditions indicate that it may be an intermediate phase occurring between the fields of wadsleyite and ringwoodite, if sufficient trivalent cations are available. The unit cell parameters have been refined at pressures up to 10.6 GPa by single-crystal X-ray diffraction in the diamond anvil cell. The refined bulk modulus for the sample containing 2.8 wt% H₂O is 145.6 ± 2.8 GPa with a K of 6.1 ± 0.7. Similar to wadsleyite and ringwoodite, hydration has a large effect on the bulk modulus. The presence of this phase in the mantle could serve to obscure the seismic expression of the phase boundary between wadsleyite and ringwoodite near 525 km. The large apparent effect of hydration on bulk modulus is consistent with hydration having a larger effect on seismic velocities than temperature in the Transition Zone.     

Field survey report on tsunami disasters caused by the 1993 Southwest Hokkaido earthquake

Detailed field work at Okushiri Island and along the southwest coast of Hokkaido has revealed quantitatively (1) the advancing direction of tsunami on land, (2) the true tsunami height (i.e., height of tsunami, excluding its splashes, as measured from the ground) and (3) the flow velocity of tsunami on land, in heavily damaged areas. When a Japanese wooden house is swept away by tsunami, bolts that tie the house to its concrete foundation resist until the last moment and become bent towards the direction of the house being carried away. The orientations of more than 850 of those bent bolts and iron pipes (all that can be measured, mostly at Okushiri Island) and fell-down direction of about 400 trees clearly display how tsunami behaved on land and caused serious damage at various places. The true tsunami height was estimated by using several indicators, such as broken tree twigs and a window pane. The flow velocity of tsunami on land was determined by estimating the hydrodynamic force exerted on a bent handrail and a bent-down guardrail by the tsunami throughin situ strength tests.Contrary to the wide-spread recognition after the tsunami hazard, our results clearly indicate that only a few residential areas (i.e., Monai, eastern Hamatsumae, and a small portion at northern Aonae, all on Okushiri Island) were hit by a huge tsunami, with true heights reaching 10 m. Southern Aonae was completely swept away by tsunami that came directly from the focal region immediately to the west. The true tsunami height over the western sea wall of southern Aonae was estimated as 3 to 4 m. Northern Aonae also suffered severe damage due to tsunami that invaded from the corner zone of the sand dune (8 m high) and tide embankment at the northern end of the Aonae Harbor. This corner apparently acted as a tsunami amplifier, and tide embankment or breakwater can be quite dangerous when tsunami advances towards the corner it makes with the coast. The nearly complete devastation of Inaho at the northern end of Okushiri Island underscored the danger of tsunami whose propagation direction is parallel to the coast, since such tsunami waves tend to be amplified and tide embankment or breakwater is constructed low towards the coast at many harbors or fishing ports. Tsunami waves mostly of 2 to 4 m in true height swept away Hamatsumae on the southeast site of Okushiri Island where there were no coastal structures. Coastal structures were effective in reducing tsunami hazard at many sites. The maximum flow velocity at northern Aonae was estimated as 10 to 18 m/s (Tsutsumi et al., 1994), and such a high on-land velocity of tsunami near shore is probably due to the rapid shallowing of the deep sea near the epicentral region towards Okushiri Island. If the advancing direction, true height, and flow velocity of tsunami can be predicted by future analyses of tsunami generation and progagation, the analyses will be a powerful tool for future assessment of tsunami disasters, including the identification of blind spots in the tsunami hazard reduction.     

Simulating bubble number density of rhyolitic pumices from Plinian eruptions: constraints from fast decompression experiments

Decompression experiments of a crystal-free rhyolitic liquid with ≈ 6.6 wt. % H₂O were carried out at a pressure range from 250 MPa to 30–75 MPa in order to characterize effects of magma ascent rate and temperature on bubble nucleation kinetics, especially on the bubble number density (BND, the number of bubbles produced per unit volume of liquid). A first series of experiments at 800°C and fast decompression rates (10–90 MPa/s) produced huge BNDs (≈ 2 × 10¹⁴ m⁻³ at 10 MPa/s ; ≈ 2 × 10¹⁵ m⁻³ at 90 MPa/s), comparable to those in natural silicic pumices from Plinian eruptions (10¹⁵–10¹⁶ m⁻³). A second series of experiments at 700°C and 1 MPa/s produced BNDs (≈ 9×10¹² m⁻³) close to those observed at 800°C and 1 MPa/s (≈ 6 × 10¹² m⁻³), showing that temperature has an insignificant effect on BNDs at a given decompression rate. Our study strengthens the theory that the BNDs are good markers of the decompression rate of magmas in volcanic conduits, irrespective of temperature. Huge number densities of small bubbles in natural silicic pumices from Plinian eruptions imply that a major nucleation event occurs just below the fragmentation level, at which the decompression rate of ascending magmas is a maximum (≥ 1 MPa/s).     

Thermal Energy Balance Analysis of the Tokyo Metropolitan Area Using a Mesoscale Meteorological Model Incorporating an Urban Canopy Model

The summer climate around the Tokyo metropolitan area has been analysed on an urban scale, and the regional characteristics of the thermal energy balance of a bayside business district in the centre of Tokyo (Otemachi) have been compared with an inland residential district (Nerima), using a mesoscale meteorological model incorporating an urban canopy model. From the results of the analysis, the mechanism of diurnal change in air temperature and absolute humidity in these areas is quantitatively demonstrated, with a focus on the thermal energy balance. Moreover, effective countermeasures against urban heat-islands are considered from the viewpoint of each region’s thermal energy balance characteristics. In addition to thermal energy outflux by turbulent diffusion, advection by sea-breezes from Tokyo Bay discharges sensible heat in Otemachi. This mitigates temperature increases during the day. On the other hand, because all sea-breezes must first cross the centre of Tokyo, it has less of a cooling effect in Nerima. As a result, the air temperature during the day in Nerima is higher than that in Otemachi.     

Multiscale characterization and modeling of granular materials through a computational mechanics avatar: a case study with experiment

Through a first-generation computational mechanics avatar that directly links advanced X-ray computed tomographic experimental techniques to a discrete computational model, we present a case study where we made the first attempt to characterize and model the grain-scale response inside the shear band of a real specimen of Caicos ooids subjected to triaxial compression. The avatar has enabled, for the first time, the transition from faithful representation of grain morphologies in X-ray tomograms of granular media to a morphologically accurate discrete computational model. Grain-scale information is extracted and upscaled into a continuum finite element model through a hierarchical multiscale scheme, and the onset and evolution of a persistent shear band is modeled, showing excellent quantitative agreement with experiment in terms of both grain-scale and continuum responses in the post-bifurcation regime. More importantly, consistency in results across characterization, discrete analysis and continuum response from multiscale calculations are found, achieving the first and long sought-after quantitative breakthrough in grain-scale analysis of real granular materials.     

Onshore tsunami deposits caused by the 1993 Southwest Hokkaido and 1983 Japan Sea earthquakes

Onshore tsunami deposits resulting from the 1993 Southwest Hokkaido and 1983 Japan Sea earthquakes were described to evaluate the feasibility of tsunami deposits for inferring paleoseismic events along submarine faults. Tsunami deposits were divided into three types, based on their composition and aerial distribution: (A) deposits consisting only of floating materials, (B) locally distributed siliclastic deposits, and (C) widespread siliclastic deposits. The most widely distributed tsunami deposits consist of the first two types. Type C deposits are mostly limited to areas where the higher tsunami runup was observed. The scale of tsunami represented by vertical tsunami runup is an important factor controlling the volume of tsunami deposits. The thickest deposits, about 10 cm, occur behind coastal dunes. To produce thick siliclastic tsunami deposits, a suitable source area, such as sand bar or dune, must be available in addition to sufficient vertical tsunami runup. Estimation of the amounts of erosion and deposition indicates that tsunami deposits were derived from both onshore and shoreface regions. The composition and grain size of the tsunami deposits strongly reflect the nature of the sedimentary materials of their source area. Sedimentary structures of the tsunami deposits suggest both low and high flow régimes. Consequently, it seems very difficult to identify tsunami deposits based only on grain size distribution or sedimentary structure of a single site in ancient successions.     

The stability of slopes and underground openings against flexural toppling and their stabilisation

Summary The stability of slopes and underground openings during and after excavation is always of great concern in the field of rock engineering. Depending upon the geologic conditions and material properties of rock and discontinuities, and the geometry of excavations and topography, various kinds of instabilities are likely to be encountered. One of these is the flexural toppling failure which has become to be known recently. A stability analysis method for slopes and underground openings under various loading conditions against the flexural toppling failure is proposed. In addition, for the stabilisation of structures, a method is suggested to take into account the reinforcement effect of fully grouted rockbolts. The applicability and validity of the proposed method is checked through model tests carried out in laboratory under well controlled conditions.