Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

On the use of the groundwater replenishment periods approach to determining aquifer parameters

Chalk crops out from a wide belt around the Paris basin, France, covering an area of about 70 000 km². In this region, the chalk presents the most important unconfined aquifer because of its extent and the size of its resources ((11–12) × 10⁹ m³ year^?1). The assessment of underground outflow depends on the vertical feeding, the infiltration and the hydrometry. This paper analyses the regional structural map, interprets the groundwater reaction under rainfall, explains the water circulation in such media where the reservoir geometry plays an essential role on aquifer response, and determines different aquifer physical parameters. Attempts are made to identify the rapid transfer of groundwater at the level of faults and their important fissures. Based on a quantitative study of the seasonal and interannual piezometry fluctuations, it is noted that the mode and the piezometric chronology events are controlled by major parameters of geological and hydrogeologic contexts, aquifer hydraulic characteristics, the position of upstream and downstream basin limits, groundwater depth and replenishment time. This paper ends with determination of groundwater physical parameters of the diffusivity (T/S, T and S) values by examination of the groundwater replenishments periods. Copyright © 2005 John Wiley & Sons, Ltd.     

Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.     

Paraglacial response of steep,sediment‐mantled slopes to post‐‘Little Ice Age’ glacier recession in the central Swiss Alps

This research assesses the morphological consequences of recent (post‐‘Little Ice Age’) paraglacial reworking of valley‐side sediment mantles in the European Alps. It aims to identify the extent and conditioning factors of slope adjustment at sites in the Swiss Alps, model the temporal pattern, and assess the rates of sediment reworking involved. Gully systems have cut into steep, high‐level lateral moraines, and debris cones have accumulated downslope. Debris flow is the dominant agent of sediment transfer. Factors controlling the extent of this activity include moraine slope gradient, relief and moisture availability. Gullies appear to have reached their maximum dimensions within ca. 50 yr of deglaciation, after which gully relief is reduced by removal of inter‐gully slopes and gully infilling (within 80–140 yr). On the most recently deglaciated terrain, minimum erosion rates average ca. 95 mm yr^?1 since gully initiation, greatly exceeding ‘normal’ erosion rates in other environments. Mean annual accumulation of a single debris cone since ice retreat was calculated to be ca. 30 mm yr^?1. Implications of these findings are applied to patterns of paraglacial sediment‐mantled slope adjustment, conceptualising paraglacial landscape response in terms of a sediment release exhaustion model, and paraglacial landform succession. Copyright © 2005 John Wiley & Sons, Ltd.     

A tracer‐based assessment of hydrological pathways at different spatial scales in a mesoscale Scottish catchment

Geochemically based hydrograph separation techniques were used in a preliminary assessment to infer how runoff processes change with landscape characteristics and spatial scale (1–233 km²) within a mesoscale catchment in upland Scotland. A two‐component end‐member mixing analysis (EMMA) used Gran alkalinity as an assumed conservative tracer. Analysis indicated that, at all scales investigated, acidic overland flow and shallow subsurface storm flows from the peaty soils covering the catchment headwaters dominated storm runoff generation. The estimated groundwater contribution to annual runoff varied from 30% in the smallest (ca 1 km²) peat‐dominated headwater catchment with limited groundwater storage, to >60% in larger catchments (>30 km²) with greater coverage of more freely draining soils and more extensive aquifers in alluvium and other drift. This simple approach offers a useful, integrated conceptualization of the hydrological functioning in a mesoscale catchment, which can be tested and further refined by focused modelling and process‐based research. However, even as it stands, the simple conceptualization of system behaviour will have significant utility as a tool for communicating hydrological issues in a range of planning and management decisions. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of the duration and direction of storm movement on infiltrating planar flow with full areal coverage

Using kinematic wave equations, analytical solutions are derived for flow due to a storm moving up or down an infiltrating plane and covering it completely. The storm duration is assumed in two ways. First, the plane is covered everywhere for the same duration by the storm. Second, the plane is covered in a linearly decreasing manner from the beginning of its coverage of the plane to the other end of the plane. By comparing the flow due to this storm with that due to a stationary storm of the same duration, the influence of storm duration, direction and velocity on flow hydrograph is investigated. It is found that storm movement has a pronounced effect on runoff hydrograph. Copyright © 2002 John Wiley & Sons, Ltd.     

VELOCITIES OF PRECURSOR SOLITON GENERATION IN SINGLE LAYER FLOW

ImODUcrIONThegeneratingveodtyofprecmeorsolitonswasstudiedbyGrimshawetal.(l986),Wu(l987)andbeetal.(1989).TheydidnotgivetherelationbeweenthegeneratingvelocitiesandtopographyintenSityandthereexistSanunknownparameterh,inthetheoryofboetal.(l989),sotheprevioustheorycannotprediCtthegeneratingvelodtAs.ThisfurtherstudyonthetheoryofXuetal.(l996)ahetoconstIUCtatheoryonthegeneratingvefodtiesofsinglelayerflowfondbytopography.InSeCtion2,thernassandenergyoftheprecursorsolitongenerationaredetendnedint…     

A Miocene coarse volcaniclastic mass-flow deposit in the Shimane Peninsula,SW Japan: product of a deep submarine eruption?

 A subaqueous volcaniclastic mass-flow deposit in the Miocene Josoji Formation, Shimane Peninsula, is 15–16 m thick, and comprises mainly blocks and lapilli of rhyolite and andesite pumices and non- to poorly vesiculated rhyolite. It can be divided into four layers in ascending order. Layer 1 is an inversely to normally graded and poorly sorted lithic breccia 0.3–6 m thick. Layer 2 is an inversely to normally graded tuff breccia to lapilli tuff 6–11 m thick. This layer bifurcates laterally into minor depositional units individually composed of a massive, lithic-rich lower part and a diffusely stratified, pumice-rich upper part with inverse to normal grading of both lithic and pumice clasts. Layer 3 is 2.5–3 m thick, and consists of interbedded fines-depleted pumice-rich and pumice-poor layers a few centimeters thick. Layer 4 is a well-stratified and well-sorted coarse ash bed 1.5–2 m thick. The volcaniclastic deposit shows internal features of high-density turbidites and contains no evidence for emplacement at a high temperature. The mass-flow deposit is extremely coarse-grained, dominated by traction structures, and is interpreted as the product of a deep submarine, explosive eruption of vesicular magma or explosive collapse of lava. Received: 10 January 1996 / Accepted: 23 February 1996     

Inbank and overbank velocity conditions in an arid zone anastomosing river

Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q_2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.     

Water temperature and conductivity variability as indicators of groundwater behaviour in complex aquifer systems in the south‐east of Spain

Groundwater temperature and electrical conductivity logs were applied in this study to gain an understanding of the hydrogeological behaviour of a complex system. This system, compartmentalized by a large number of fractures, which may even cause geothermal anomalies, is that of the Campo de Dalías in the south‐east of the Iberian Peninsula. In addition to this complexity one must consider the different levels with aquifer properties and their varying lithologies, which in some cases are in direct contact with the sea. Vertical sections based upon the various logs allow lateral correlation to be made between the anomalies detected, and provide information on the existence of both ascending and descending vertical flows, and the inflow of sea water into the system. Analysis of environmental isotopes (¹⁸O and D) indicate the areas of recharge with the greatest flow velocity. No seasonal variation in the isotopes was observed, which suggests that the residence time of the groundwater in the aquifer is several years. Radiocarbon dating gave estimates of generally long residence times, although in some sectors it is possible to detect the current recharge based on its tritium content, from rainwater infiltration or seawater intrusion. Copyright © 2002 John Wiley & Sons, Ltd.