Modeling the interaction between fluid flow and particle movement with discontinuous deformation analysis (DDA) method

Coupled hydro‐mechanical processes in granular media represent the interaction between solid particles movement and fluid flow during external and/or internal loading. The processes attract attention in geotechnical engineering since they cause many disasters such as landslide, slope collapse, boiling/quicksand and soil liquefaction, etc. This study presents a numerical method combining discontinuous deformation analysis (DDA) for mechanical calculation and finite element method for fluid flow simulation to model the interaction between solid particles' movement and fluid flow from microscopic point of view. The term of hydraulic pressure was newly formulated and introduced to the original DDA. The simplified examples are given to verify the new method, and the computational results correlate well with the theoretical calculations, further development is also considered. Copyright © 2010 John Wiley & Sons, Ltd.     

Interannual variability of North Pacific eastern subtropical mode water formation in the 1990s derived from a 4-dimensional variational ocean data assimilation experiment

The interannual variability of eastern subtropical mode water (ESTMW) formation in the North Pacific is examined using a new ocean dataset constructed by a 4-dimensional variational data assimilation experiment covering the decade of the 1990s. The volume of newly formed ESTMW varies due to interannual variability in the following three physical processes taking place in the surface layer: (1) convergence in the transport of surface saline water induced by Ekman flow in the vicinity of the formation region, (2) thermal stratification in the preconditioning phase in association with the insolation anomaly induced largely by low-level cloud coverage, and (3) wintertime surface cooling in the eastern subtropics. We find that, in addition to the surface forcing, the properties of both the ESTMW and the upper mixed-layer water are broadly controlled by the volume of the new ESTMW component, and that the variations in the upper mixed-layer water affect the properties of ESTMW formed in the following winter. Due to the combined effect of these processes, the ESTMW subducts down to subsurface layers with a wide range of σ_θ values lying between 24.8 and 25.4 and with significant interannual variation in water mass formation.     

2011年日本东北近海地震震源上方的位移

2011年3月11日,一次大的板间地震(MW9.0)发生在日本东北宫城县附近的板块交界处。根据余震分布推测的震源区为约500km长、200km宽的近海区(Japan Meteo-rological Agency,2011)。     

Geochemistry and petrogenesis of carbonatites from South Nam Xe,Lai Chau area,northwest Vietnam

This paper presents a study of the petrography, mineral chemistry, geochemistry, and Sr–Nd–Pb–C–O isotope systematics of carbonatite dykes and associated rocks from the northeastern part of the Song Da intracontinental rift in South Nam Xe (northwest Vietnam) aimed at constraining the origin of the carbonatite magmas. The carbonatites are characterized by SiO₂ < 12.18 wt.% and by wide ranges in FeO, MgO and CaO content that define them as calciocarbonatite and ferrocarbonatite. On U–Th–Pb isochron diagrams, whole rocks and mineral separates from the ferrocarbonatites form linear arrays corresponding to ages of 30.2–31.6 Ma (Rupelian, Oligocene). The South Nam Xe carbonatites are extremely enriched in Sr, Ba, and light rare earth elements (LREE), and depleted in high field strength elements (HFSE) (e.g. Ti, Nb, Ta, Zr and Hf). The age–corrected Sr–Nd–Pb isotope ratios and C isotope data are relatively uniform (⁸⁷Sr/⁸⁶Sr_(t) = 0.708193–0.708349; ¹⁴³Nd/¹⁴⁴Nd_(t) = 0.512250–0.512267; ε_Nd(t) = ?6.46 to ?6.80; ²⁰⁶Pb/²⁰⁴Pb_(t) = 18.26–18.79; ²⁰⁷Pb/²⁰⁴Pb_(t) = 15.62–15.64; ²⁰⁸Pb/²⁰⁴Pb_(t) = 38.80–39.38; δ¹³C_V-PDB = –2.7?‰ to ?4.1?‰). These isotopic compositions indicate source contamination that occurred before the production of the carbonatite magmas, and did not change noticeably during or after emplacement. The variation in oxygen isotopes is consistent with the change in mineral compositions and trace element abundances: the lower δ¹⁸O values (9.1–11.0?‰) coupled with Sr-rich, Mn-poor calcite, and igneous textures such as triple junctions among calcite grain boundaries, define a magmatic origin. However, the elevated δ¹⁸O values of the ferrocarbonatites (12.0–13.3?‰) coupled with a volatile-bearing mineral assemblages (including REE fluorcarbonates, sulfates, sulfides and fluorite) may be due to interaction with meteoric water during low-temperature alteration. High δ¹³C values and Sr–Pb ratios, and low Rb/Sr (0.00014–0.00301), Sm/Nd (0.089–0.141) and ¹⁴³Nd/¹⁴⁴Nd ratios, coupled with very high Sr-Nd concentrations, suggest the involvement of an enriched mantle component, which probably resulted from metasomatism due to the migration of subducted material. Because of the lack of tectonic data and the limited number of samples studied, this conclusion is still ambiguous and requires further study.     

A development of the discontinuous deformation analysis for rock fall analysis

Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental study on countermeasure against liquefaction-induced flotation of manhole using recycled materials packed in sandbags

The 2004 Niigata-ken Chuetsu earthquake in Japan caused serious damage to sewage facilities such as uplift of manholes and settlement of pavement above backfill soil for pipes.This paper deals with shake table tests in a 1 g gravity field on application of recycled materials for ground improvement to mitigation of liquefaction-induced flotation of manhole during earthquakes.The recycled materials used in tests were tire chips made of waste tires and crushed gravels made of waste reinforced concrete,and they were packed in sandbags.From the test results,it was confirmed that the recycled materials packed in sandbags could be treated as one of the countermeasures to restrain the flotation of manholes and settlement of ground surrounded by sandbags.     

Damage of 1975 Haicheng <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>7.3 earthquake in the Korean Peninsula and its effect on Kyushu,Japan

The 1975 Haicheng M _S7.3 earthquake affected whole Korean Peninsula and even Kyushu, Japan. In this paper isoseismal map of this earthquake was drawn for these areas by sorting out and analyzing data systematically collected, furthermore, the characteristics of the earthquake damage were discussed. In the northern Korean Peninsula, the intensity is V in great majority of this area except local area with intensity VI, and the damage is basically in accordance with typical characteristics of corresponding intensity grade. In the southern Korean Peninsula, the intensity is IV in large part of the area, however, in some big cities such as Seoul etc., not only high buildings shook strongly, but also some people appeared perceptions such as feeling dizzy and so on, and electrical service was interrupted on one or two districts because of transformer tripping. These phenomena could be caused by the effect of long-period wave generated from a large earthquake on far-field.     

Robust and exploratory analysis of active mesoscale tectonic zones in Japan utilizing the nationwide GPS array

A monitoring GPS array recently developed in Japan can yield nationwide maps of active inland tectonic zones (ATZs) on a mesoscale, approximately 70 to several hundred kilometers in lateral extent. But it has been difficult to characterize ATZs in Japan, as they are in fact operational on multiple scales and our efforts are often hindered by various irregularities in the data. The key to overcoming these problems would be to gain an insight into the available data before any precise kinematic modeling is performed with indefinite assumptions. In this study, horizontal velocity fields, deduced from the nationwide GPS array, were treated with a set of techniques in robust smoothing and exploratory data analysis that brought out exceptionally powerful mesoscale ATZs, and made them easier to characterize. The resolved ATZs were then retrospectively monitored to study their regional and temporal variations, using a set of approx. 840 observation stations, about 30 km apart, for a 4-year series of fixed observation time-intervals, 810 days each. The smoothing operation involved three steps: (1) imputation of the velocity fields for the purpose of anti-aliasing, (2) robust smoothing of the velocity fields with the median operative, and (3) visualization of deformation-rate distributions in several coordinate independent parameters, and post-filtering. The geometrical resolvability of mesoscale ATZs was confirmed by calibrating the smoothing scheme against synthetic tectonic boundary models before it was applied to the case study in Japan. ATZs in Japan, which are essentially visible as systematic deviations in the velocity fields on the International Terrestrial Reference Frame (ITRF) and as strain rate anomalies, were highlighted sharply along some known tectonic zones, chains of active volcanoes, and areas above low seismic velocity anomalies in the crust and upper mantle, all of which generally paralleled the offshore trench axes. The geometrical agreements among the mapped ATZs and the physical anomalies in the crust are presumably due to their common structural weakness on the mesoscale. In the four main islands of Japan, all but 30–40% of the strain rate anomalies persisted during the entire 6 years of the case study period, while the rest sporadically appeared or disappeared in a period from several months to a few years. The transient shifts in the deformation rates were remarkably synchronous with some nearby major tectonic episodes: large earthquakes and slow events. Differential plate coupling strengths along the subduction zones can also be inferred from the persistent pattern of rotational strain rate anomalies forming clockwise and counterclockwise pairs along the Pacific. Our empirical observations suggest that the first-order features of interseismic crustal deformations in Japan can be characterized as collateral processes behaving in response to fluctuations of the tectonic stresses on multiple scales, likely influenced by changes of plate coupling strengths on the contiguous subduction faults.