Mathematical modelling of groundwater flow at Sellafield, UK

Sellafield in West Cumbria was a potential site for the location of the UK's first underground repository for radioactive, intermediate level waste (ILW). The repository was to lie around 650 m beneath the ground surface within rocks of the Borrowdale volcanic group (BVG), a thick suite of SW dipping, fractured, folded and metamorphosed Ordovician meta-andesites and ignimbrites. These are overlain by an onlapping sequence of Carboniferous and Permo-Triassic sediments. In situ borehole measurements showed that upward trending fluid pressure gradients exist in the area of the potential repository site, and that there are three distinct fluid types in the subsurface; fresh, saline and brine (at depth, to the west of the site). Simulations of fluid flow in the Sellafield region were undertaken with a 2D, steady-state, coupled fluid and heat flow simulation code (OILGEN). In both simplified and geologically complex models, topographically driven flow dominated the regional hydrogeology. Fluids trended persistently upwards through the potential repository site. The dense brine to the west of the site promoted upward deflection of topographically driven groundwaters. The inclusion in hydrogeological models of faults and variably saline sub-surface fluids was essential to the accurate reproduction of regional hydraulic head variations. Sensitivity analyses of geological variables showed that the rate of groundwater flow through the potential repository site was dependent upon the hydraulic conductivity of the BVG, and was unaffected by the hydraulic conductivity of other hydrostratigraphic units. Calibration of the model was achieved by matching simulated subsurface pressures to those measured in situ. Simulations performed with BVG hydraulic conductivity 100 times the base case median value provided the “best-fit” comparison between the calculated equivalent freshwater head and that measured in situ, regardless of the hydraulic conductivity of other hydrostratigraphic units. Transient mass transport simulations utilising the hydraulic conductivities of this “best fit” simulation showed that fluids passing through the potential repository site could reach the surface in 15 000 years. Simple safety case implications drawn from the results of the study showed that the measured BVG hydraulic conductivity must be less than 0.03 m year⁻¹ to be simply declared safe. Recent BVG hydraulic conductivity measurements showed that the maximum BVG hydraulic conductivity is around 1000 times this safety limit.     

A model for Holocene retreat of the West Antarctic Ice Sheet

Marine ice sheets are grounded on land which was below sea level before it became depressed under the ice-sheet load. They are inherently unstable and, because of bedrock topography after depression, the collapse of a marine ice sheet may be very rapid. In this paper equations are derived that can be used to make a quantitative estimate of the maximum size of a marine ice sheet and of when and how rapidly retreat would take place under prescribed conditions. Ice-sheet growth is favored by falling sea level and uplift of the seabed. In most cases the buttressing effect of a partially grounded ice shelf is a prerequisite for maximum growth out to the edge of the continental shelf. Collapse is triggered most easily by eustatic rise in sea level, but it is possible that the ice sheet may self-destruct by depressing the edge of the continental shelf so that sea depth is increased at the equilibrium grounding line.Application of the equations to a hypothetical “Ross Ice Sheet” that 18,000 yr ago may have covered the present-day Ross Ice Shelf indicates that, if the ice sheet existed, it probably extended to a line of sills parallel to the edge of the Ross Sea continental shelf. By allowing world sea level to rise from its late-Wisconsin minimum it was possible to calculate retreat rates for individual ice streams that drained the “Ross Ice Sheet.” For all the models tested, retreat began soon after sea level began to rise (～15,000 yr B.P.). The first 100 km of retreat took between 1500 and 2500 yr but then retreat rates rapidly accelerated to between 0.5 and 25 km yr^?1, depending on whether an ice shelf was present or not, with corresponding ice velocities across the grounding line of 4 to 70 km yr^?1. All models indicate that most of the present-day Ross Ice Shelf was free of grounded ice by about 7000 yr B.P. As the ice streams retreated floating ice shelves may have formed between promontories of slowly collapsing stagnant ice left behind by the rapidly retreating ice streams. If ice shelves did not form during retreat then the analysis indicates that most of the West Antarctic Ice Sheet would have collapsed by 9000 yr B.P. Thus, the present-day Ross Ice Shelf (and probably the Ronne Ice Shelf) serves to stabilize the West Antarctic Ice Sheet, which would collapse very rapidly if the ice shelves were removed. This provides support for the suggestion that the 6-m sea-level high during the Sangamon Interglacial was caused by collapse of the West Antarctic Ice Sheet after climatic warming had sufficiently weakened the ice shelves. Since the West Antarctic Ice Sheet still exists it seems likely that ice shelves did form during Holocene retreat. Their effect was to slow and, finally, to halt retreat. The models that best fit available data require a rather low shear stress between the ice shelf and its sides, and this implies that rapid shear in this region encouraged the formation of a band of ice with a preferred crystal fabric, as appears to be happening today in the floating portions of fast bounded glaciers.Rebound of the seabed after the ice sheet had retreated to an equilibrium position would allow the ice sheet to advance once more. This may be taking place today since analysis of data from the Ross Ice Shelf indicates that the southeast corner is probably growing thicker with time, and if this persists then large areas of ice shelf must become grounded. This would restrict drainage from West Antarctic ice streams which would tend to thicken and advance their grounding lines into the ice shelf.     

An integrated hydrogeological study to support sustainable development and management of groundwater resources: a case study from the Precambrian Crystalline Province,India

The rapid expansion of agriculture, industries and urbanization has triggered unplanned groundwater development leading to severe stress on groundwater resources in crystalline rocks of India. With depleting resources from shallow aquifers, end users have developed resources from deeper aquifers, which have proved to be counterproductive economically and ecologically. An integrated hydrogeological study has been undertaken in the semi-arid Madharam watershed (95 km²) in Telangana State, which is underlain by granites. The results reveal two aquifer systems: a weathered zone (maximum 30 m depth) and a fractured zone (30–85 m depth). The weathered zone is unsaturated to its maximum extent, forcing users to tap groundwater from deeper aquifers. Higher orders of transmissivity, specific yield and infiltration rates are observed in the recharge zone, while moderate orders are observed in an intermediate zone, and lower orders in the discharge zone. This is due to the large weathering-zone thickness and a higher sand content in the recharge zone than in the discharge zone, where the weathered residuum contains more clay. The NO₃ ^? concentration is high in shallow irrigation wells, and F^? is high in deeper wells. Positive correlation is observed between F^? and depth in the recharge zone and its proximity. Nearly 50 % of groundwater samples are unfit for human consumption and the majority of irrigation-well samples are classed as medium to high risk for plant growth. Both supply-side and demand-side measures are recommended for sustainable development and management of this groundwater resource. The findings can be up-scaled to other similar environments.     

Modeling of groundwater flow at depth in crystalline rock beneath a moving ice-sheet margin, exemplified by the Fennoscandian Shield, Sweden

On-going geological disposal programs for spent nuclear fuel have generated strong demands for investigation and characterization of deep-lying groundwater systems. Because of the long time scales for which radiological safety needs to be demonstrated in safety assessment applications, an analysis of the hydrogeological performance of the geosphere system during glacial climate conditions is needed. Groundwater flow at depth in crystalline rock during the passage of an ice-sheet margin is discussed based on performed groundwater-flow modeling of two bedrock sites, Forsmark and Laxemar, in the Fennoscandian Shield, Sweden. The modeled ice sheet mimics the Weichselian ice sheet during its last major advance and retreat over northern Europe. The paper elaborates and analyzes different choices of top boundary conditions at the ice sheet–subsurface interface (e.g. ice-sheet thickness and ice-margin velocity) and in the proglacial area (presence or lack of permafrost) and relates these choices to available groundwater-flow-model hydraulic output and prevailing conceptual hydrogeochemical models of the salinity evolution at the two sites. It is concluded that the choice of boundary conditions has a strong impact on results and that the studied sites behave differently for identical boundary conditions due to differences in their structural-hydraulic properties.     

Groundwater balance in the Khor Arbaat basin, Red Sea State, eastern Sudan

The Khor Arbaat basin is the main source of potable water supply for the more than 750,000 inhabitants of Port Sudan, eastern Sudan. The variation in hydraulic conductivity and storage capacity is due to the heterogeneity of the sediments, which range from clay and silt to gravely sand and boulders. The water table rises during the summer and winter rainy seasons; it reaches its lowest level in the dry season. The storage capacity of the Khor Arbaat aquifer is estimated to be 21.75?×?10⁶ m³. The annual recharge through the infiltration of flood water is about 1.93?×?10⁶ m³. The groundwater recharge, calculated as underground inflow at the ‘upper gate’, is 1.33?×?10⁵ m³/year. The total annual groundwater recharge is 2.06?×?10⁶ m³. The annual discharge through underground outflow at the ‘lower gate’ (through which groundwater flows onto the coastal plain) is 3.29?×?10⁵ m³/year. Groundwater discharge due to pumping from Khor Arbaat basin is 4.38?×?10⁶ m³/year on average. The total annual groundwater discharge is about 4.7?×?10⁶ m³. A deficit of 2.6?×?10⁶ m³/year is calculated. Although the total annual discharge is twice the estimated annual recharge, additional groundwater flow from the fractured basement probably balances the annual groundwater budget since no decline is observed in the piezometric levels.     

Consideration of the use of direct recharge in analytical models of flood-wave response: a case study of the Merguellil alluvial aquifer,central Tunisia

In a stream–aquifer system, adequate knowledge of the aquifer hydraulic parameters is paramount for the analysis of groundwater/surface-water interaction. In such systems, analytical flood-wave response models are commonly used to assess these parameters. However, when an important means of direct recharge (precipitation) occurs simultaneously with the stream flood waves, the classical flood-wave response model (i.e. when only stream flood waves are considered as an input signal) leads to overestimated and variable values of the aquifer hydraulic parameters. In this study, an analytical flood-wave response model was implemented to determine hydraulic parameters of the Merguellil alluvial aquifer (central Tunisia) taking into account the recharge rate as a second input signal. An array of three piezometers and a river-stage gauge equipped with an automated Diver data logger were installed in the floodplain of the Merguellil Wadi. The results showed a good agreement between the observed and modelled hydrographs when the recharge rate was considered, which allows the model to provide the same value of diffusivity (50–70 m² h^?1) in the three piezometers and for all flood events being tested. However, when the recharge rate was not taken into account, the model gave an overestimated diffusivity with a very poor hydrograph match. Overall, the methodology described in this paper can be applicable to any area in the world because of its simplicity and its better estimation of aquifer hydraulic parameters.     

Provenance of buried esker groundwater: the case of Vars-Winchester esker aquifer,Eastern Ontario,Canada

An innovative mode of groundwater recharge to a buried esker aquifer is considered. The current conceptual model affords a natural safeguard to underlying aquifers from the overlying muds. A hypothesis of groundwater recharge to a buried esker aquifer via preferential pathways across its overlying muds is tested here by heuristic numerical one-dimensional and two-dimensional modeling simulations. The hypothesis has been tested against two other conventionally accepted scenarios involving: (1) distal esker outcrop areas and (2) remote shallow-bedrock recharge areas. The main evidence comes from documented recharge pressure pulses in the overlying mud aquitard and in the underlying esker hydraulic-head time series for the Vars-Winchester esker aquifer in Eastern Ontario, Canada. These perturbations to the potentiometric surface are believed to be the aquifer response to recharge events. The migration rate of these pressure pulses is directly related to the hydraulic diffusivity of the formation. The measured response time and response amplitude between singular radar precipitation events and well hydrographs constituted the heuristic model calibration targets. The main evidence also includes mud-layering deformation (water escape features) which was observed in seismic surveys of the over-esker muds. These disturbed stratigraphic elements provide a realistic mechanism for migrating water to transit through the muds. The effective hydraulic conductivities of these preferential pathways in the muds were estimated to be between 2?×?10^?6 and 7?×?10^?6 m/s. The implications of these findings relate to the alleged natural safeguard of these overlying muds.     

Characterization of the shallow groundwater system in an alpine watershed: Handcart Gulch,Colorado, USA

Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10^?5 m/s. Discharge was estimated at 1.28×10^?3 km³. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10^?5–2.0×10^?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10^?9–2.0×10^?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.     

What can be learned from sequential multi-well pumping tests in fracture-karst media? A case study in Zhangji, China

A fracture-karst aquifer is a karst aquifer with a fractured rock matrix, and its parameters are difficult to determine. Two sequential pumping tests in a fracture-karst aquifer system at the Zhangji well field of China are considered, one carried out before (in 2000) and one after (in 2005) the operation of a pumping station in the well field (2003–2005). The sequential tests serve multiple purposes. First, they provide a cross check of the parameters obtained. Second, they can be used to assess the effect of long-term groundwater exploitation of the aquifer. A three-dimensional finite-element transient flow model has been developed to simulate groundwater flow at the site. Generally good agreement has been found between the simulated and observed hydraulic heads for both tests. The hydraulic parameters obtained from the 2005 test are generally consistent with their counterparts from the 2000 test. However, a small but steady increase of hydraulic conductivities from 2000 to 2005 at the site has been observed. A 10-year prediction of groundwater resources has been made and indicates that the well field can accommodate the proposed 8.0?×?10⁴ m³/day exploitation rate under relative drought conditions without causing a steady decline of groundwater levels.     

Groundwater flow dynamics of weathered hard-rock aquifers under climate-change conditions: an illustrative example of numerical modeling through the equivalent porous media approach in the north-western Pyrenees (France)

A numerical groundwater model of the weathered crystalline aquifer of Ursuya (a major water source for the north-western Pyrenees region, south-western France) has been computed based on monitoring of hydrological, hydrodynamic and meteorological parameters over 3 years. The equivalent porous media model was used to simulate groundwater flow in the different layers of the weathered profile: from surface to depth, the weathered layer (5?·?10^?8?≤?K?≤?5?·? 10^?7 m s^?1), the transition layer (7?·?10^?8?≤?K?≤?1?·? 10^?5 m s^?1, the highest values being along major discontinuities), two fissured layers (3.5?·?10^?8?≤?K?≤?5?·?? 10^?4 m s^?1, depending on weathering profile conditions and on the existence of active fractures), and the hard-rock basement simulated with a negligible hydraulic conductivity (K = 1 10 ^?9 ). Hydrodynamic properties of these five calculation layers demonstrate both the impact of the weathering degree and of the discontinuities on the groundwater flow. The great agreement between simulated and observed hydraulic conditions allowed for validation of the methodology and its proposed use for application on analogous aquifers. With the aim of long-term management of this strategic aquifer, the model was then used to evaluate the impact of climate change on the groundwater resource. The simulations performed according to the most pessimistic climatic scenario until 2050 show a low sensitivity of the aquifer. The decreasing trend of the natural discharge is estimated at about ?360 m³ y^?1 for recharge decreasing at about ?5.6 mm y^?1 (0.8 % of annual recharge).     

Residence times of shallow groundwater in West Africa: implications for hydrogeology and resilience to future changes in climate

Although shallow groundwater (<50 mbgl) sustains the vast majority of improved drinking-water supplies in rural Africa, there is little information on how resilient this resource may be to future changes in climate. This study presents results of a groundwater survey using stable isotopes, CFCs, SF₆, and ³H across different climatic zones (annual rainfall 400–2,000 mm/year) in West Africa. The purpose was to quantify the residence times of shallow groundwaters in sedimentary and basement aquifers, and investigate the relationship between groundwater resources and climate. Stable-isotope results indicate that most shallow groundwaters are recharged rapidly following rainfall, showing little evidence of evaporation prior to recharge. Chloride mass-balance results indicate that within the arid areas (<400 mm annual rainfall) there is recharge of up to 20 mm/year. Age tracers show that most groundwaters have mean residence times (MRTs) of 32–65 years, with comparable MRTs in the different climate zones. Similar MRTs measured in both the sedimentary and basement aquifers suggest similar hydraulic diffusivity and significant groundwater storage within the shallow basement. This suggests there is considerable resilience to short-term inter-annual variation in rainfall and recharge, and rural groundwater resources are likely to sustain diffuse, low volume abstraction.     

Drainage beneath ice sheets: groundwater–channel coupling,and the origin of esker systems from former ice sheets

The nature of the drainage system beneath ice sheets is crucial to their dynamic behaviour but remains problematic. An experimentally based theory of coupling between groundwater and major channel systems is applied to the esker systems in the area occupied the last ice sheet in Europe, which we regard as a fossil imprint of major longitudinal drainage channels. We conclude that the large-scale distribution and spacing of major eskers is consistent with the theory of groundwater control, in which esker spacing is partly controlled by the transmissivity of the bed.It is concluded that esker patterns reflect the large-scale organisation of the subglacial drainage pattern in which channel development is coupled to groundwater flow and to the ice sheet's dynamic regime. The theory is then used to deduce: basal meltwater recharge rates and their spatial variability from esker spacing in an area in which the ice sheet was actively streaming during its final retreat; patterns of palaeo-groundwater flow and head distribution; and the seasonally varying magnitude of discharge from stream tunnels at the retreating ice sheet margin. Major channel/esker systems appear to have been stable at least over several hundred of years during the retreat of the ice sheet, although major dynamic events are demonstrably associated with major shifts in the hydraulic regime.Modelling suggests: that glaciation can stimulate deep groundwater circulation cells that are spatially linked to channel locations, with groundwater flow predominantly transverse to ice flow; that the circulation pattern has the potential to create large-scale anomalies in groundwater chemistry; and that the spacing of channels will change through the glacial cycle, influencing water pressures in stream tunnels, subglacial hydraulic gradients and effective pressure. If the latter is reduced sufficiently, it could trigger enhanced bed deformation, thus coupling drainage to ice sheet movement. It suggests the possibility of distinctive phases of sediment deformation and drumlin mobilisation during a glacial cycle.     

Relative importance and chemical effects of diffuse and focused recharge in an eogenetic karst aquifer: an example from the unconfined upper Floridan aquifer, USA

Karst aquifer studies often focus on allogenic water inputs and large conduit flow. However, diffuse recharge can be significant, particularly in unconfined eogenetic karst aquifers that retain high matrix permeability. This study examines an unconfined region of the upper Floridan aquifer (USA) that hosts a sinking stream, its resurgence, and a large conduit system. Daily diffuse recharge was approximated using a water-budget method and ranged from 17% of precipitation during a low precipitation year to >53% during the highest precipitation year, illustrating the highly variable nature of diffuse recharge in this region. The total allogenic input via the sinking stream over the 5 years of the study was significantly larger than the volume of diffuse recharge. However, only about 2% of the allogenic recharge flows from the conduit into the surrounding aquifer. That flow is restricted to storm events when hydraulic heads in the conduits exceed those in the surrounding aquifer. The estimated volume of dissolution is similar for allogenic recharge and diffuse recharge to the unconfined region surrounding the conduits, but dissolution from the diffuse recharge is distributed over a larger area than dissolution from allogenic recharge. These results exemplify how recharge type impacts flow and water–rock interactions in eogenetic karst aquifers.     

The formation of natural cryogenic brines

The source of salts in the Ca-chloridic, hypersaline brines (up to 190 g Cl L⁻¹) occurring in crystalline basement rocks in the Canadian, Fennoscandian and Bohemian Shields and their evolution have been investigated and reported. The Cl-Br-Na relationship indicates that these waters have been concentrated from seawater, by freezing during glacial times. The Na/Cl ratio (0.25 to 0.35) in the more saline fluids is compatible with cooling down to −30°C, where the most saline waters have been concentrated by a factor of 25 to 30 relative to the parent seawater.The brines formed from seawater within cryogenic troughs, along the subarctic continental margins, around ice sheets. The depressions within which the brines formed are the cryogenic analogues of the classic, evaporitic lagoon. One million years suffice to saturate with brine a 2000km-radius by 1km-depth rock volume at an H₂O removal rate of only 2.8 mm/yr. Density-induced brine migration on a continental scale takes place via fissures below the ice.Our calculations, that were performed on a hypothetical ice sheet with dimensions compatible with the Laurentide ice sheet, demonstrate that during 1m.y., a 60m thick cryogenic sediment section could have formed. However, the precipitated minerals (mirabilite and hydrohalite) are repeatedly dispersed by the advance and retreat of the ice sheet, dissolved by melt water-seawater mixtures, and eroded during postglacial uplift, leaving almost no trace in the geological record.The cryogenic brines formed intermittently during and between glacial periods. The repeating advance and retreat of the ice sheets exerted a major control on the direction and intensity of brine flow. The cryogenic concentration of seawater and the migration of brine towards the center of the glaciostatic depression occurred mainly during the build up of the ice sheet, while reversal of the water flow from the center of the cryogenic basin outwards happened upon deglaciation. The flow of the waters in the subsurface was, inevitably, accompanied by significant dilution with melt water from the ice sheets.Using a “granitic” U concentration of 4 ppm and a (Ca-Mg mass balance based) rock/water ratio anywhere between 3.4 and 6.8 kg L⁻¹, a few hundred thousand years of brine-rock interaction are sufficient for the growth of ¹²⁹I in the most saline Canadian Shield brine to its present concentration (3.4×10⁸ atoms ¹²⁹I L⁻¹). Hence, both the formation of the saline fluids and their emplacement in their present sites occurred most likely within the Pleistocene.The young age calculated for cryogenic brines in crystalline shields and the dynamic water flow therein should raise concern about the planning and construction of high-grade nuclear waste repositories in such rocks, which are already under way.     

利用RQD估算岩体不同深度的平均渗透系数和平均变形模量

渗透系数是进行裂隙岩体渗流模拟的必备参数,变形模量是工程岩体数值模拟的必备参数。大量研究中忽视了渗透系数和变形模量随深度变化这一重要规律,从而影响模拟结果的可靠性。为此,探讨了利用极易获取的RQD（岩石质量指标）资料估算不同深度的渗透系数和变形模量的可行性。通过以某花岗岩体为例,研究发现,RQD均值随深度增大,渗透系数均值随深度减小,其相关性很好。因此,利用RQD估算不同深度的平均渗透系数是可行的。根据RQD随深度的变化,利用经验公式估算了不同深度的变形模量均值和变化范围。估算得到的变形模量与实测结果基本一致。     

Ribbed moraine stratigraphy and formation in southern Finnish Lapland

Characteristics of ribbed moraines, the dominating moraine type in southern Finnish Lapland, have been studied in detail. The ridges are composed of several till units, of which the bottommost units consist of mature basal tills and the surficial parts are enriched with local, short‐transport rock fragments and boulders in till and at the surface of ridges. As a result of this re‐examination a two‐step model of the formation process of ribbed moraines is presented. In the first stage, while cold‐based conditions prevailed, both the bottommost part of the ice sheet and the frozen, substrate fractured under compressive ice flow. Following glacial transport of fractured blocks and formation of the transverse ridge morphology, erosion between the ridges continued owing to freeze–thaw process under variable pressure conditions. In the areas with a low pre‐existing till sheet, the process caused quarrying of the bedrock surface and subsequent deposition of rock fragments and boulders under high pressure on the next ridge. The most suitable conditions for ribbed moraine formation existed during Late Weichselian deglaciation, after the Younger Dryas when the climate warmed very quickly, leading to an imbalance between a warm glacier surface and a cold base. Copyright © 2006 John Wiley & Sons, Ltd.     

The role of meltwater in glacial processes

Water plays a dominant role in many glacial processes and the erosional, depositional and climatic significance of meltwaters and associated fluvioglacial processes cannot be overemphasized. At its maximum extent c. 20,000 years ago, the volume of the Laurentide ice sheet was 33 × 10⁶ km³ (about the same as the volume of all ice present today on planet Earth). The bulk of this was released as water in little more than 10,000 years. Pulses of meltwater flowing to the Atlantic Ocean from large ice dammed lakes altered thermohaline circulation of the world's oceans and global climate. One such discharge event via Hudson Bay at 8200 years BP released 160,000 km³ of water in 12 months. Global sea levels recovered from glacial maximum low stands reached at about 20,000 years ago at an average rate of 15 m per thousand years but estimates of shorter term rates suggest as much as 20 m sea level rise in 1000 years and for short periods, rates as high as 4 m per hundred years. Meltwaters played a key role in lubricating ice sheet motion (and thus areal abrasion) across the inner portions of the ice sheet where it slid over rigid crystalline bedrock of the Canadian Shield. The recharge of meltwater into the ice sheets bed was instrumental in generating poorly sorted diamict sediments (till) by sliding-induced shearing and deformation of overpressured sediment and soft rock. The transformation of overpressured till into hyperconcentrated slurries in subglacial channels may have generated a highly effective erosional tool for selective overdeepening and sculpting of bedrock substrates. Some workers credit catastrophic subglacial ‘megafloods’ with the formation of drumlins and flutes on till surfaces. Subglacial melt river systems were instrumental in reworking large volumes of glaciclastic sediment to marine basins; it has been estimated that less than 6% of the total volume of glaciclastic sediment produced during the Pleistocene remains on land. Fluvioglacial and glaciolacustrine sediments and landforms dominate large tracts of the ‘glacial’ landscape in North America. The recharge of subglacial meltwater into underlying bedrock and sediment aquifers created transient reversals in the long-term equilibrium flow directions of basinal fluids. With regard to pre-Pleistocene glacial record, meltwaters moved enormous volumes of terrestrial ‘glaciclastic’ sediment to marine basins and thus played a key role in preserving a record of glaciation, a record otherwise almost entirely lost on land.     

Hydrogeology of a montane headwater groundwater system downgradient of a coal-mine waste rock dump: Elk Valley,British Columbia,Canada

Steelmaking-coal waste rock placed in mountain catchments in the Elk Valley, British Columbia, Canada, drain constituents of interest (CIs) to surface water downgradient of the waste rock dumps. The role of groundwater in transporting CIs in the headwaters of mountain catchments is not well understood. This study characterizes the physical hydrogeology of a portion of a 10-km² headwater catchment (West Line Creek) downgradient of a 2.7-km² waste rock dump placed over a natural headwater valley-bottom groundwater system. The study site was instrumented with 13 monitoring wells. Drill core samples were collected to determine subsurface lithology and geotechnical properties. The groundwater system was characterized using field testing and water-level monitoring. The valley-bottom sediments were composed of unconsolidated glacial and meltwater successions (<64 m thick) deposited as a series of cut and fill structures overlying shale bedrock. An unconfined basal alluvial aquifer located above fractured bedrock was identified as the primary conduit for groundwater flow toward Line Creek (650 m from the toe of the dump). Discharge through the basal alluvial aquifer was estimated using the geometric mean hydraulic conductivity (±1 standard deviation). These calculations suggest groundwater discharge could account for approximately 15% (ranging from 2 to 60%) of the total water discharged from the watershed. The residence time from the base of the waste rock dump to Line Creek was estimated at <3 years. The groundwater system was defined as a snowmelt (i.e., nival) regime dominated by direct recharge (percolation of precipitation) across the catchment.     

The hydraulic conductivity of the fractures intersecting Cambrian sandstone rock masses, central Jordan

This paper outlines the hydraulic characteristics of fractured rock masses and their implication in engineering works. The hydraulic behavior of subsurface fracture systems has been evaluated by means of hydraulic testing using packer tests and by fracture analysis. A comparison of the borehole results with those of surface fracture mapping provides a reasonable correlation between the two methods of measuring fractured rock hydraulic conductivity. The mean hydraulic conductivity value obtained from the boreholes is 36.5 LU (9.26᎒^-5 m/s), while the mean value of hydraulic conductivity obtained from field mapping of fracture data is in the order of 1᎒^-5 m/s. Based on the hydraulic conductivity values the sandstone rock mass can be considered medium to highly conductive; nevertheless, it seems to be almost impervious at greater depth. The empirical relationships which have been derived between hydraulic conductivity and both rock quality designation (RQD) and rock mass rating (RMR) indices indicated that the mean value of hydraulic conductivity of the rock mass could be estimated to be in the order of 10^-5 m/s, which is confirmed by the packer tests.