The effect of terrace geology on ground‐water movement and on the interaction of ground water and surface water on a mountainside near Mirror Lake,New Hampshire,USA

The west watershed of Mirror Lake in the White Mountains of New Hampshire contains several terraces that are at different altitudes and have different geologic compositions. The lowest terrace (FSE) has 5 m of sand overlying 9 m of till. The two next successively higher terraces (FS2 and FS1) consist entirely of sand and have maximum thicknesses of about 7 m. A fourth, and highest, terrace (FS3) lies in the north‐west watershed directly adjacent to the west watershed. This highest terrace has 2 m of sand overlying 8 m of till. All terraces overlie fractured crystalline bedrock. Numerical models of hypothetical settings simulating ground‐water flow in a mountainside indicated that the presence of a terrace can cause local ground‐water flow cells to develop, and that the flow patterns differ based on the geologic composition of the terrace. For example, more ground water moves from the bedrock to the glacial deposits beneath terraces consisting completely of sand than beneath terraces that have sand underlain by till. Field data from Mirror Lake watersheds corroborate the numerical experiments. The geology of the terraces also affects how the stream draining the west watershed interacts with ground water. The stream turns part way down the mountainside and passes between the two sand terraces, essentially transecting the movement of ground water down the valley side. Transects of water‐table wells were installed across the stream's riparian zone above, between, and below the sand terraces. Head data from these wells indicated that the stream gains ground water on both sides above and below the sand terraces. However, where it flows between the sand terraces the stream gains ground water on its uphill side and loses water on its downhill side. Biogeochemical processes in the riparian zone of the flow‐through reach have resulted in anoxic ground water beneath the lower sand terrace. Results of this study indicate that it is useful to understand patterns of ground‐water flow in order to fully understand the flow and chemical characteristics of both ground water and surface water in mountainous terrain. Copyright © 2007 John Wiley & Sons, Ltd.     

Mathematical and numerical filtration–advection–dispersion model of miscible grout propagation in saturated porous media

The development of a predictive model of behaviour of porous media during injection of miscible grout, taking into account convection, dilution and filtration of grout solution with interstitial water, as well as consolidation aspects, is presented. Model assumptions are reviewed and discussed first. During the establishment of the model, we insist on surface terms and their physical relevance in expressing adsorption effects. Constitutive laws such as Fick's law for diffusive mass transport, hydrodynamic dispersion tensor dealing with miscibility, are modified by taking into account filtration effects. A new surface term appears in mass balance equations as a consequence of filtration. According to the filtration laws used, an initial filtration rate is estimated on the basis of a one‐dimensional experimental campaign. The field equations are discretized by using Galerkin finite element and θ‐scheme standard method. For transport equation, Streamline Upwind Petrov Galerkin method is employed to prevent numerical oscillations. Lastly, confrontation of numerical results with laboratory experiments constitutes a first step to validate the model on a realistic basis. Copyright © 2001 John Wiley & Sons, Ltd.     

Moraine development at a small High‐Arctic valley glacier: Rieperbreen,Svalbard

Ice‐cored lateral and frontal moraine complexes, formed at the margin of the small, land‐based Rieperbreen glacier, central Svalbard, have been investigated through field observations and interpretations of aerial photographs (1936, 1961 and 1990). The main focus has been on the stratigraphical and dynamic development of these moraines as well as the disintegration processes. The glacier has been wasting down since the ‘Little Ice Age’ (LIA) maximum, and between 1936 and 1990 the glacier surface was lowered by 50–60 m and the front retreated by approximately 900 m. As the glacier wasted, three moraine ridges developed at the front, mainly as melting out of sediments from debris‐rich foliation and debris‐bands formed when the glacier was polythermal, probably during the LIA maximum. The disintegration of the moraines is dominated by wastage of buried ice, sediment gravity‐flows, meltwater activity and some frost weathering. A transverse glacier profile with a northward sloping surface has developed owing to the higher insolation along the south‐facing ice margin. This asymmetric geometry also strongly affects the supraglacial drainage pattern. Lateral moraines have formed along both sides of the glacier, although the insolation aspect of the glacier has resulted in the development of a moraine 60 m high along its northern margin. Copyright © 2001 John Wiley & Sons, Ltd.     

Late Glacial fluvial response of the Niers‐Rhine (western Germany) to climate and vegetation change

The Niers valley was part of the Rhine system that came into existence during the maximum Saalian glaciation and was abandoned at the end of the Weichselian. The aim of the study was to explain the Late Pleniglacial and Late Glacial fluvial dynamics and to explore the external forcing factors: climate change, tectonics and sea level. The sedimentary units have been investigated by large‐scale coring transects and detailed cross‐sections over abandoned channels. The temporal fluvial development has been reconstructed by means of geomorphological relationships, pollen analysis and ¹⁴C dating. The Niers‐Rhine experienced a channel pattern change from braided, via a transformational phase, to meandering in the early Late Glacial. This change in fluvial style is explained by climate amelioration at the Late Pleniglacial to Late Glacial transition (at ca. 12.5 k ¹⁴C yr BP) and climate‐related hydrological, lithological and vegetation changes. A delayed fluvial response of ca. 400 ¹⁴C yr (transitional phase) was established. The channel transformations are not related to tectonic effects and sea‐level changes. Successive river systems have similar gradients of ca. 35–40 cm km^?1. A meandering river system dominated the Allerød and Younger Dryas periods. The threshold towards braiding was not crossed during the Younger Dryas, but increased aeolian activity has been observed on the Younger Dryas point bars. The final abandonment of the Niers‐Rhine was dated shortly after the Younger Dryas to Holocene transition. Traces of Laacher See pumice have been found in the Niers valley, indicating that the Niers‐Rhine was still in use during the Younger Dryas. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling the effects of land‐use change on water and salt delivery from a catchment affected by dryland salinity in south‐east Australia

A comprehensive framework for the assessment of water and salt balance for large catchments affected by dryland salinity is applied to the Boorowa River catchment (1550 km²), located in south‐eastern Australia. The framework comprised two models, each focusing on a different aspect and operating on a different scale. A quasi‐physical semi‐distributed model CATSALT was used to estimate runoff and salt fluxes from different source areas within the catchment. The effects of land use, climate, topography, soils and geology are included. A groundwater model FLOWTUBE was used to estimate the long‐term effects of land‐use change on groundwater discharge. Unlike conventional salinity studies that focus on groundwater alone, this study makes use of a new approach to explore surface and groundwater interactions with salt stores and the stream. Land‐use change scenarios based on increased perennial pasture and tree‐cover content of the vegetation, aimed at high leakage and saline discharge areas, are investigated. Likely downstream impacts of the reduction in flow and salt export are estimated. The water balance model was able to simulate both the daily observed stream flow and salt load at the catchment outlet for high and low flow conditions satisfactorily. Mean leakage rate of about 23·2 mm year^?1 under current land use for the Boorowa catchment was estimated. The corresponding mean runoff and salt export from the catchment were 89 382 ML year^?1 and 38 938 t year^?1, respectively. Investigation of various land‐use change scenarios indicates that changing annual pastures and cropping areas to perennial pastures is not likely to result in substantial improvement of water quality in the Boorowa River. A land‐use change of about 20% tree‐cover, specifically targeting high recharge and the saline discharge areas, would be needed to decrease stream salinity by 150 µS cm^?1 from its current level. Stream salinity reductions of about 20 µS cm^?1 in the main Lachlan River downstream of the confluence of the Boorowa River is predicted. The FLOWTUBE modelling within the Boorowa River catchment indicated that discharge areas under increased recharge conditions could re‐equilibrate in around 20 years for the catchment, and around 15 years for individual hillslopes. Copyright © 2004 John Wiley & Sons, Ltd.     

Extension of the Griffith's fracture criteria to saturated clays

Inglis [1] has solved the problem of distribution of stress in an elastic plate around an elliptical hole. His works clarify the role of cracks in the failure of an elastic material. However, his solution cannot be applied to saturated clay because he considers only total stresses, while, in saturated clay, the criterion of rupture should be expressed in terms of effective and not total stresses. The solution of Atkinson and Craster [2] using Biot's poroelasticity theory, shows that there is no high pore pressure in the vicinity of the crack tips for saturated clay. The major difference between this approach and the Biot's theory of is that, in saturated clay, strain is a function of the variation of the effective stress [3], while, in poroelastic media, strain is only a function of the variation of the total stress [4, Equation 2.2]. Also in their solution there is continuity between the pore fluid and the inner fluid in the crack. Their solution is valid for poroelastic media involving a movement of the pore fluid. In our solution there is no movement of the pore fluid (Undrained condition). In this paper we have solved the same problem as Inglis [1], but for the particular case of saturated clay obeying elastic law. By solving this problem we obtained the expressions for pore pressure, effective stress, total stress and displacements. The results show that not only the total stress but also the pore pressure and the effective stress are also high in the vicinity of the crack tips. A new failure criterion, based on Griffith's strain energy principle [5] and maximum tensile stress [6], valid for saturated clay is developed in this paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Displacement‐controlled method and its applications to material non‐linearity

For the analysis of non‐linear problems, the displacement‐controlled method (DCM) has a more extensive application scope and more powerful abilities than the load‐controlled method (LCM). However, difficulties of the DCM's procedure not amenable to most finite element implementations of the conventional LCM have restricted its applications in geomechanics. By means of Sherman–Morrison's theorem, the solution of DCM is improved. The improved procedure is characterized by high efficiency, good numerical stability and a programme structure similar to LCM. Two aspects of applications of DCM are illustrated. The first application is to compute the response of a structure under a given load level like the conventional finite element analysis. The second application is to trace the equilibrium path of a structure under a given load distribution type. A simple but effective algorithm is presented for automatically adjusting the step length in tracing the equilibrium path. Examples illustrate that the proposed procedures are suited for modelling complicated non‐linear problems in geomechanics. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of property variation and modelling on convection in a fluid overlying a porous layer

An accurate computational analysis is presented for the onset of thermal convection in a two‐layer system which is comprised of a saturated layer of porous material described by Darcy's law, over which lies a layer of the same saturating fluid. The two‐layer system is heated from below and the upper (fluid) surface is allowed to be fixed or stress free. The onset of convection may have a bi‐modal nature in which convection may be dominated by the porous medium or by the fluid depending on the depths of the relative layers, but this is strongly controlled by material parameters. The effect of variation of relevant fluid and porous material properties is investigated in detail, as is the effect of the interface boundary condition between the fluid and the porous medium. Copyright © 2001 John Wiley & Sons, Ltd.     

Three‐dimensional finite element analysis of the interaction between tunneling and pile foundations

This paper concerns analysis of the impact of construction of urban tunnels on adjacent pile foundations. It is carried out using an elastoplastic three‐dimensional finite element modelling. Numerical simulations are performed in two stages, which concern, respectively, the application of the pile axial loading and the construction of the tunnel in presence of the pile foundations. Analysis is carried out for both single piles and groups of piles. Results of numerical simulations show that tunneling induces significant internal forces in adjacent piles. The distribution of internal forces depends mainly on the position of the pile tip regarding the tunnel horizontal axis and the distance of the pile axis from the centre of the tunnel. Analysis of the interaction between tunneling and a group of piles reveals a positive group effect with a high reduction of the internal forces in rear piles. Copyright © 2002 John Wiley & Sons, Ltd.