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.     

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.     

Depositional and post‐depositional history of warm stage deposits at Knocknacran,Co. Monaghan,Ireland: implications for preservation of Irish last interglacial deposits

Organic‐rich deposits, uncovered during overburden removal from mantled gypsum karst at Knocknacran opencast gypsum mine, Co. Monaghan, are the best candidate to date for a last interglacial record in Ireland. The two till and organic‐rich deposits (preserved at different quarry elevations) were emplaced on to a Tertiary dolerite surface during high‐energy flood events and subsequently folded and faulted by movement towards sinkholes in underlying gypsum. Uranium–thorium disequilibrium dating suggests that the organic‐rich deposits in the upper section were hydrologically isolated at ca. 41 ka and those in the lower section at ca. 86 ka. Interpretation of the pollen content, although tentative because of the depositional and post‐depositional history of the material, suggests that the organic material originated in a warm stage possibly warmer than the post‐Eemian interstadials. The unusual setting of preservation may indicate that in situ, last interglacial deposits have generally been removed by erosion in Ireland. Copyright © 2004 John Wiley & Sons, Ltd.     

Himalayan magmatism and porphyry copper–molybdenum mineralization in the Yulong ore belt, East Tibet

Summary ?The NW–SE-trending Yulong porphyry Cu–Mo ore belt, situated in the Sanjiang0 area of eastern Tibet, is approximately 400 km long and 35 to 70 km wide. Complex tectonic and magmatic processes during the Himalayan epoch have given rise to favorable conditions for porphyry-type Cu–Mo mineralization. Porphyry masses of the Himalayan epoch in the Yulong ore belt are distributed in groups along regional NW–SE striking tectonic lineaments. They were emplaced mainly into Triassic and Lower Permian sedimentary-volcanic rocks. K–Ar und U–Pb isotopic datings give an intrusion age range of 57–26 Ma. The porphyries are mainly of biotite monzogranitic and biotite syenogranitic compositions. Geological and geochemical data indicate that the various porphyritic intrusions in the belt had a common or similar magma source, are metaluminous to peraluminous, Nb–Y–Ba-depleted, I-type granitoids, and belong to the high-K calc-alkaline series. Within the Yulong subvolcanic belt a number of porphyry stocks bear typical porphyry type Cu–Mo alteration and mineralization. The most prominent porphyry Co–Mo deposits include Yulong, Malasongduo, Duoxiasongduo, Mangzong and Zhanaga, of which Yulong is one of the largest porphyry Cu (Mo) deposits in China with approximately 8 × 10⁶ tons of contained Cu metal. Hydrothermal alteration at Yulong developed around a biotite–monzogranitic porphyry stock that was emplaced within Upper Triassic limestone, siltstone and mudstone. The earliest alteration was due to the effects of contact metamorphism of the country rocks and alkali metasomatism (potassic alteration) within and around the porphyry body. The alteration of this stage was accompanied by a small amount of disseminated and veinlet Cu–Mo sulfide mineralization. Later alteration–mineralization zones form more or less concentric shells around the potassic zone, around which are distributed a phyllic or quartz–sericite–pyrite zone, a silicification and argillic zone, and a propylitic zone. Fluid inclusion data indicate that three types of fluids were involved in the alteration–mineralization processes: (1) early high temperature (660–420 °C) and high salinity (30–51 wt% NaCl equiv) fluids responsible for the potassic alteration and the earliest disseminated and/or veinlet Cu–Mo sulfide mineralization; (2) intermediate unmixed fluids corresponding to phyllic alteration and most Cu–Mo sulfide mineralization, with salinities of 30–50 wt% NaCl equiv and homogenization temperatures of 460–280 °C; and (3) late low to moderate temperature (300–160 °C) and low salinity (6–13 wt% NaCl equiv) fluids responsible for argillic and propylitic alteration. Hydrogen and oxygen isotopic studies show that the early hydrothermal fluids are of magmatic origin and were succeeded by increasing amounts of meteoric-derived convective waters. Sulfur isotopes also indicate a magmatic source for the sulfur in the early sulfide mineralization, with the increasing addition of sedimentary sulfur outward from the porphyry stock. Received August 29, 2001; revised version accepted May 1, 2002 Published online: November 29, 2002     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.