Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

Consolidation of double‐layered ground with vertical drains

Vertical drains are usually installed in subsoil consisting of several layers. Due to the complex nature of the problem, over the past decades, the consolidation properties of multi‐layered ground with vertical drains have been analysed mainly by numerical methods. An analytical solution for consolidation of double‐layered ground with vertical drains under quasi‐equal strain condition is presented in this paper. The main steps for the computation procedure are listed. The convergence of the series solution is discussed. The comparisons between the results obtained by the present analytical method and the existing numerical solutions are described by figures. The orthogonal relation for the system of double‐layered ground with vertical drains is proven. Finally, some consolidation properties of double‐layered ground with vertical drains are analysed. Copyright © 2001 John Wiley & Sons, Ltd.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Relict mountain slope deposits of northern Portugal: facies,sedimentogenesis and environmental implications

This paper deals with the formation processes and the palaeoenvironmental significance of relict slope deposits located on the uppermost part of the north Portugal mountains. For this purpose, seven key sites representative of the different lithofacies have been selected and analysed in detail. The data show that three main dynamic processes are responsible for the emplacement of regional fossil slope deposits: runoff, debris flows and dry grain flows. The ubiquity of these processes and the lack of frost‐related features or landforms do not support the existence of severe Pleistocene climates in this part of the lberian Peninsula as postulated by previous work. Pedological data gathered at one of the study sites show that a subalpine environment was probably present at 700–800 m altitude between 29 and 14 kyr. Using data from the Pyrenees Mountains, a 6.5 to 12°C depression in mean annual temperature has been tentatively postulated for this Pleniglacial period. Copyright © 2003 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Multistep matched-field inversion for broad-band data from ASIAEX2001

This paper discusses the results of geoacoustic inversion carried out using explosive charge data from the Asian Seas International Acoustic Experiment (ASIAEX) East China Sea (ECS) Experiment. A multifrequency incoherent matched-field inversion processor and a genetic algorithm (GA) are used for the inversion. A multistep matched field inversion approach is presented, which makes use of the varying sensitivities of wave fields at various frequencies to reduce the inversion problem into a sequence of smaller inversions with fewer unknowns to estimate at each stage. Different parameters are estimated using data at different frequencies according to their sensitivities. Inversion results for different areas in the ECS region are summarized and compared with core data.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Abundance and isotopic composition of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite

Abstract— Abundances and isotopic compositions of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite were measured. The abundance ratios of spallogenic components in metal reveal a ³He deficiency which is due to the diffusive loss of parent isotopes, that is, tritium (Tilles, 1963; Schultz, 1967). The diffusive loss likely has been induced by thermal heating by the Sun during cosmic‐ray exposure (～160 Ma; Lavielle et al, 1999). Thermal process such as impact‐induced partial loss may have affected the isotopic composition of spallogenic Ne. The ¹²⁹Xe/¹³¹Xe ratio of cosmogenic components in the metal indicates an enhanced production of epi‐thermal neutrons. The abundance ratios of spallogenic components in the graphite reveal that it contained small amounts of metal and silicates. The isotopic composition of heavy noble gases in graphite itself was obtained from graphite treated with HF/HCl. The isotopic composition of the etched graphite shows that it contains two types of primordial Xe (i.e., Q‐Xe and El Taco Xe). The isotopic heterogeneity preserved in the Bohumilitz graphite indicates that the Bohumilitz graphite did not experience any high‐temperature event and, consequently, must have been emplaced into the metal at subsolidus temperatures. This situation is incompatible with an igneous model as well as the impact melting models for the IAB‐IIICD iron meteorites as proposed by Choi et al. (1995) and Wasson et al (1980).