Groundwater isoscapes in a montane headwater catchment show dominance of well‐mixed storage

We conducted an integrated groundwater–surface water monitoring programme in a 3.2‐km² experimental catchment in the Scottish Highlands by sampling all springs, seepages, and wells in six, spatially extensive synoptic surveys over a 2‐year period. The catchment has been glaciated, with steep hillslopes and a flat valley bottom. There is around 70% glacial drift cover in lower areas. The solid geology, which outcrops at higher elevations, is granite and metamorphic schist. The springs and seepages generally occur at the contact between the solid geology and drift or at breaks of slopes in the valley bottom. Samples were analysed for stable isotopes, Gran alkalinity and electrical conductivity. Despite the surveys encompassing markedly different antecedent conditions, the isotopic composition of groundwater at each location exhibited limited temporal variability, resulting in a remarkable persistence of spatial patterns indicating well‐mixed shallow, groundwater stores. Moreover, line‐conditioned excess values derived from the isotope data indicated no evidence of fractionation affecting the groundwater, which suggests that most recharge occurs in winter. The alkalinity and electrical conductivity of groundwater reflected geological differences in the catchment, being highest where more weatherable calcareous rocks outcrop at higher altitudes in the catchment. Springs draining these areas also had the most variable isotope composition, which indicated that they have shorter residence times than the drift covered part of the catchment. The study showed that even in geologically heterogeneous upland catchments, groundwater can be characterized by a consistent isotopic composition, reflecting rapid mixing in the recharge zone. Our work, thus, emphasizes the critical role of groundwater in upland catchments and provides tracer data that can help constrain quantitative groundwater models.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

Toward a new conceptual model for groundwater flow in merokarst systems: Insights from multiple geophysical approaches

Merokarst aquifers — relatively thin (<1–2 m) karstified carbonate units interbedded between mudstone, shale, or sandstone — constitute a significant proportion of carbonate terrain and underlie a large portion of the west- and south-central USA, yet few advances have been made in our understanding of porosity development and flow-path generation in these complex systems in decades. Toward this end, we used a multi-geophysical approach at the well-studied Konza Prairie Biological Station (KPBS), a part of the larger Flint Hills (25,734 km²), underlain by thin limestone units (1–2 m thick) interbedded with mudstone/shale units (2–4 m thick), to elucidate hydrologic connectivity and potential controls on known groundwater flow directions. We combined electrical resistivity tomography (ERT), surface and borehole nuclear magnetic resonance (NMR), and ground penetrating radar (GPR) measurements across a low order catchment where over 25 boreholes and groundwater wells sampling perched aquifers could be used to constrain interpretation of lithology, potential flow paths, and permeability. Data revealed that groundwater export may be an unappreciated component of lateral-flow-dominated models used to represent merokarst in that: (a) potentiometric surfaces from two limestone units showed groundwater flows toward a hydrologic depression, opposite the direction of stream flow, in the upstream portion of the catchment, (b) long term measures of groundwater levels revealed a greater variance and overall water storage in this same upstream area compared to wells near the outlet, and (c) ERT and NMR results indicate pronounced lateral heterogeneity within a given unit, suggestive of a greater degree of vertical hydrological connectivity than usually considered for horizontally-layered merokarst. Our data suggest vertical connectivity can shunt water to depth in these “sandwiched” merokarst aquifers, yielding atypical groundwater flow directions and unrealized deep export of weathering solutes and carbon.     

Application of electrical and electromagnetic methods to study sedimentary covers in high mountain areas

The results of geophysical studies conducted with selected electrical and electromagnetic methods in the Kondratowa Valley in the Tatra Mountains (the Carpathian Mountains, Poland) are presented in the article. The surveys were performed with the following methods: electrical resistivity tomography (ERT), georadar (GPR) and conductivity meter (CM). The objective of the noninvasive geophysical measurements was to determine the thickness of the Quaternary postglacial sediments that fill the bottom of the valley and to designate the accumulation of boulders deposited on Quaternary sediments. The results of ERT surveys conducted along the axis of the valley allowed to determine the changeability of the thickness of the postglacial sediments and allowed to designate a few areas of occurrence of boulders. The ERT, GPR and CM surveys conducted across the valley allowed to designate with high accuracy the thickness of the accumulation of boulders sliding down the valley bottom from the couloirs surrounding the valley.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial organization of groundwater dynamics and streamflow response from different hydropedological units in a montane catchment

Groundwater dynamics play an important role in runoff generation and hydrologic connectivity between hillslopes and streams. We monitored a network of 14 shallow groundwater (GW) wells in a 3.2 km² experimental catchment in the Scottish Highlands. Wells were placed in three contrasting landscape units with different hydropedological characteristics and different topographic positions relative to the stream network, encompassing a catena sequence from freely draining podzols on steeper hillslopes to increasingly thick peats (histosols) in the valley bottom riparian zone. GW dynamics were characterized by statistical analyses of water table fluctuations, estimation of variabilities in lag times and hysteresis response in relation to streamflow. The three landscape units had distinct storage–discharge relationships and threshold responses with a certain GW level above which lateral flow dominates. Steeper hillslopes with freely draining podzols were characterized by GW fluctuations of around 150 cm in the underlying drift. GW usually showed peak response up to several hours after stream flow. During persistent wet periods the water table remained in the soil profile for short spells and connected shallow flow paths in the near surface horizons to the lower hillslopes. In the peaty gleys in the lower foot slopes, GW was characterized by a water table generally within 20 cm of the soil surface, though at some locations this could fall to 50 cm in extreme dry periods. GW responses were usually a few hours prior to the stream responses. In riparian peats, the water table was also usually less than 20 cm deep and responded several hours before the stream. These riparian peat soils remain at, or very near saturation with near‐continuous GW–surface water connectivity. In contrast, the steeper slopes remain disconnected for prolonged periods and need large recharge events to overcome storage thresholds. GW responses vary seasonally, and landscape controls on the spatial organization of GW dynamics are strongest at low flows and in small events. During wettest periods, limited storage and extensive saturation weaken such controls. This study demonstrated that montane catchments can have highly dynamic GW stores, which are important in generating both storm flows and baseflows. Copyright © 2016 John Wiley & Sons, Ltd.     

Using high resolution tracer data to constrain water storage,flux and age estimates in a spatially distributed rainfall‐runoff model

Models simulating stream flow and conservative tracers can provide a representation of flow paths, storage distributions and mixing processes that is advantageous for many predictive purposes. Compared with models that only simulate stream flow, tracer data can be used to investigate the internal consistency of model behaviour and to gain insight into model performance. Here, we examine the strengths and weaknesses of a data‐driven, spatially distributed tracer‐aided rainfall‐runoff model. The model structure allowed us to assess the influence of landscape characteristics on the routing and mixing of water and tracers. The model was applied to a site in the Scottish Highlands with a unique tracer data set; ~4 years of daily isotope ratios in stream water and precipitation were available, as well as 2 years of weekly soil and ground water isotopes. The model structure was based on an empirically based, lumped tracer‐aided model previously developed for the catchment. The best model runs were selected from Monte Carlo simulations based on dual calibration criteria using objective functions for both stream isotopes and discharge at the outlet. Model performance for these criteria was reasonable (Nash–Sutcliffe efficiencies for discharge and isotope ratios were ~0.4–0.6). The model could generally reproduce the variable isotope signals in the soils of the steeper hill slopes where storage was low, and damped isotope responses in valley bottom cells with high storage. The model also allowed us to estimate the age distributions of internal stores, water fluxes and stream flow. Average stream water age was ~1.6 years, integrating older groundwater in the valley bottom and dynamic younger soil waters. By tracking water ages and simulating isotopes, the model captured the changes in connectivity driven by distributed storage dynamics. This has substantially improved the representation of spatio‐temporal process dynamics and gives a more robust framework for projecting environmental change impacts. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

The topographic signature of a major typhoon

In August 2009, the typhoon Morakot, characterized by a cumulative rainfall up to 2884 mm in about three days, triggered thousands of landslides in Taiwan. The availability of LiDAR surveys before (2005) and after (2010) this event offers a unique opportunity to investigate the topographic signatures of a major typhoon. The analysis considers the comparison of slope–area relationships derived by LiDAR digital terrain models (DTMs). This approach has been successfully used to distinguish hillslope from channelized processes, as a basis to develop landscape evolution models and theories, and understand the linkages between landscape morphology and tectonics, climate, and geology. We considered six catchments affected by a different degree of erosion: three affected by shallow and deep‐seated landslides, and three not affected by erosion. For each of these catchments, 2 m DTMs were derived from LiDAR data. The scaling regimes of local slope versus drainage area suggested that for the catchments affected by landslides: (i) the hillslope‐to‐valley transitions morphology, for a given value of drainage area, is shifted towards higher value of slopes, thus indicating a likely migration of the channelized processes and erosion toward the catchment boundary (the catchment head becomes steeper because of erosion); (ii) the topographic gradient along valley profiles tends to decrease progressively (the valley profile becomes gentler because of sediment deposition after the typhoon). The catchments without any landslides present a statistically indistinguishable slope–area scaling regime. These results are interesting since for the first time, using multi‐temporal high‐resolution topography derived by LiDAR, we demonstrated that a single climate event is able to cause significant major geomorphic changes on the landscape, detectable using slope–area scaling analysis. This provides new insights about landscape evolution under major climate forcing. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrological response to thinning a small jarrah (Eucalyptus marginata) forest catchment

The jarrah (Eucalyptus marginata) forest on a small catchment was thinned early in 1983 to study the effect on catchment hydrology. The thinning reduced canopy cover, basal area and stocking about two-thirds. Rainfall during the pre and post-treatment periods was 21% and 10%, respectively, below the long-term average. Streamflow increased from 0.5% of rainfall (4.3 mm) before thinning to 7.6% of rainfall (90 mm) 9 years after thinning. Streamflow duration increased, with the largest increases in streamflow in the wet winter months of June-October. The deep groundwater level at a midslope location increased by 8 m and at a valley location by 4 m in the 8 years after thinning. There was no indication of a new equilibrium being reached for either streamflow or groundwater, thus further increases in streamflow and groundwater level are likely.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

A multi‐method field experiment to determine local groundwater flow in a glacier forefield

We implemented multiple independent field techniques to determine the direction and velocity of groundwater flow at a specific stream reach in a glacier forefield. Time‐lapse experiments were conducted using two electrical resistivity tomography (ERT) lines installed in a cross pattern. A circular array of groundwater tubes was also installed to monitor groundwater flow via discrete salt injections. Both inter‐borehole and ERT results confirmed this stream section as a losing reach and enabled quantification of the flow direction. Both techniques yielded advection velocities varying between 5.7 and 21.8 m/day. Estimates of groundwater flow direction and velocity indicated that groundwater infiltrates from the stream nearby and not from the adjacent lateral moraine. Groundwater age estimated from radon concentration measurements supported this hypothesis. Despite uncertainties inherent to each of the methods deployed, the combination of multiple field techniques allowed drawing consistent conclusions about local groundwater flow. We thus regard our multi‐method approach as a reliable way to characterize the two‐dimensional groundwater flow at sites where more invasive groundwater investigation techniques are difficult to carry out and local heterogeneities can make single measurements unreliable. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the impact of extensive irrigation on the groundwater resources

Drastic groundwater resource depletion due to excessive extraction for irrigation is a major concern in many parts of India. In this study, an attempt was made to simulate the groundwater scenario of the catchment using ArcSWAT. Due to the restriction on the maximum initial storage, the deep aquifer component in ArcSWAT was found to be insufficient to represent the excessive groundwater depletion scenario. Hence, a separate water balance model was used for simulating the deep aquifer water table. This approach is demonstrated through a case study for the Malaprabha catchment in India. Multi‐site rainfall data was used to represent the spatial variation in the catchment climatology. Model parameters were calibrated using observed monthly stream flow data. Groundwater table simulation was validated using the qualitative information available from the field. The stream flow was found to be well simulated in the model. The simulated groundwater table fluctuation is also matching reasonably well with the field observations. From the model simulations, deep aquifer water table fluctuation was found very severe in the semi‐arid lower parts of the catchment, with some areas showing around 60 m depletion over a period of eight years. Copyright © 2012 John Wiley & Sons, Ltd.     

Electrical resistivity characterization and defect detection on a geosynthetic clay liner (GCL) on an experimental site

In this paper we analyze the onsite characterization of a geosynthetic clay liner (GCL) that serves to ensure the impermeability of a landfill cap by DC electrical methods. The imaging of the GCL geoelectrical properties is a challenging problem because it is a very thin (between 4 and 7 mm thick) and resistive layer (from 100,000 to 2,000,000 Ω·m) depending on meteorological conditions and aging. We compare results obtained using electrical resistivity tomography (ERT) using two different kinds of arrays (dipole–dipole DD and Wenner–Schlumberger) on an experimental site with engineered defects. To confirm these results and to find the real onsite GCL resistivity we have performed sampling of the posterior distribution of this parameter using vertical electrical sounding (VES) inversions. Different VES methods were extracted from ERT with DD array and converted into a Schlumberger array.As a main conclusion the dipole–dipole array provides a better resistivity resolution of the defects than the Wenner–Schlumberger array. On ERT images, the defect detection seems to be impossible if the GCL has very high resistivity, as it happened when it was put in place. Taking into account the equivalence rules, the inversions are in both cases (ERT and VES) compatible. The GCL resistivity estimated from PSO (particle swarm optimization) varies from 3.0 10⁵ to 1.10⁶ Ω·m depending on saturation conditions during the twenty first months of its placing. Then, the resistivity dropped to 4.10⁴–9.10⁴ Ω·m, indicating a probable chemical damage of the GCL due to aging. Finally the fact that the VES inversions are solved via PSO sampling allows for the detection of a very thin and resistive layer and opens the possibility of performing micro VES surveys along the landfill to detect possible GCL defects.     

Tracing solute infiltration using a combined method of dye tracer test and electrical resistivity tomography in an undisturbed forest soil profile

An accurate prediction of solute infiltration in a soil profile is important in the area of environmental science, groundwater and civil engineering. We examined the infiltration pattern and monitored the infiltration process using a combined method of dye tracer test and electrical resistivity tomography (ERT) in an undisturbed field soil (1 m × 1 m). A homogeneous matrix flow was observed in the surface soil (A horizon), but a preferential flow along macropores and residual rock structure was the dominant infiltration pattern in the subsurface soil. Saturated interflow along the slopping boundaries of A and C1 horizons and of an upper sandy layer and a lower thin clay layer in the C horizon was also observed. The result of ERT showed that matrix flow started first in A horizon and then the infiltration was followed by the preferential flows along the sloping interfaces and macropores. The ERT did not show as much detail as the dye‐stained image for the preferential flow. However, the area with the higher staining density where preferential flow was dominant showed a relatively lower electrical resistivity. The result of this study indicates that ERT can be applied for the monitoring of solute transportation in the vadose zone. Copyright © 2010 John Wiley & Sons, Ltd.     

Geophysical and Hydrologic Studies of Lake Seepage Variability

Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.     

Hydrologic functioning of the deep critical zone and contributions to streamflow in a high‐elevation catchment: Testing of multiple conceptual models

High‐elevation mountain catchments are often subject to large climatic and topographic gradients. Therefore, high‐density hydrogeochemical observations are needed to understand water sources to streamflow and the temporal and spatial behaviour of flow paths. These sources and flow paths vary seasonally, which dictates short‐term storage and the flux of water in the critical zone (CZ) and affect long‐term CZ evolution. This study utilizes multiyear observations of chemical compositions and water residence times from the Santa Catalina Mountains Critical Zone Observatory, Tucson, Arizona to develop and evaluate competing conceptual models of seasonal streamflow generation. These models were tested using endmember mixing analysis, baseflow recession analysis, and tritium model “ages” of various catchment water sources. A conceptual model involving four endmembers (precipitation, soil water, shallow, and deep groundwater) provided the best match to observations. On average, precipitation contributes 39–69% (55 ± 16%), soil water contributes 25–56% (41 ± 16%), shallow groundwater contributes 1–5% (3 ± 2%), and deep groundwater contributes ~0–3% (1 ± 1%) towards annual streamflow. The mixing space comprised two principal planes formed by (a) precipitation‐soil water‐deep groundwater (dry and summer monsoon season samples) and (b) precipitation‐soil water‐shallow groundwater (winter season samples). Groundwater contribution was most important during the wet winter season. During periods of high dynamic groundwater storage and increased hydrologic connectivity (i.e., spring snowmelt), stream water was more geochemically heterogeneous, that is, geochemical heterogeneity of stream water is storage‐dependent. Endmember mixing analysis and ³H model age results indicate that only 1.4 ± 0.3% of the long‐term annual precipitation becomes deep CZ groundwater flux that influences long‐term deep CZ development through both intercatchment and intracatchment deep groundwater flows.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Scale variation of post‐glacial sediment yield in Chilliwack Valley,British Columbia

The Holocene volumetric sediment budget is estimated for coarse textured sediments (sand and gravel) in a large, formerly glaciated valley in southwest British Columbia. Erosion is estimated by compiling volumetric loss estimated in digital elevation models (DEMs) of gullied topography and by applying a non‐linear diffusion model on planar, undissected hillslopes. Estimates of steepland yield are based on estimates of post‐glacial deposition volumes in fans, cones and deltas at the outlets of low‐order tributary catchments. Erosion of post‐glacial fans and tributary valley fills is estimated by reconstructing formerly continuous surfaces. Results are classed by catchment order and compared across scales of contributing area, revealing declining specific sediment yield (in m³ km^?2 a^?1) with catchment area for the smaller tributaries (<10 km²) and increasing specific sediment yield for larger tributaries and Chilliwack Valley itself. Approximately 60% of mobilized sediment is redeposited in first‐ to third‐order catchments, with lesser proportions stored at the outlets of higher order catchments. A simple network routing model emphasizes the significant sediment flux contributions from colluvium, drift blankets and gullies in steeper terrain. As this material is deposited at junctions within the lower drainage network, an increasing proportion of material is derived from remnant valley fills and para‐glacial fans in the major valleys. Yield from lower‐order, steepland catchments tends to remain in storage, indefinitely sequestered on footslopes. These observations have implications for modelling the post‐glacial sediment balance amongst catchments of varying size. After 10⁴ years, the system remains in disequilibrium. The critical linkage lies between low‐order, hillslope catchments (相似文献   

Arctic rockwall retreat rates estimated using laboratory‐calibrated ERT measurements of talus cones in Longyeardalen,Svalbard

Holocene rockwall retreat rates quantify integral values of rock slope erosion and talus cone evolution. Here we investigate Holocene rockwall retreat of exposed arctic sandstone cliffs in Longyeardalen, central Svalbard and apply laboratory‐calibrated electrical resistivity tomography (ERT) to determine talus sediment thickness. Temperature–resistivity functions of two sandstone samples are measured in the laboratory and compared with borehole temperatures from the talus slope. The resistivity of the higher and lower‐porosity sandstone at relevant borehole permafrost temperatures defines a threshold range that accounts for the lithological variability of the dominant bedrock and debris material. This helps to estimate the depth of the transition from higher resistivities of ice‐rich debris to lower resistivities of frozen bedrock in the six ERT transects. The depth of the debris–bedrock transition in ERT profiles is confirmed by a pronounced apparent resistivity gradient in the raw data plotted versus depth of investigation. High‐resolution LiDAR‐scanning and ERT subsurface information were collated in a GIS to interpolate the bedrock surface and to calculate the sediment volume of the talus cones. The resulting volumes were referenced to source areas to calculate rockwall retreat rates. The rock mass strength was estimated for the source areas. The integral rockwall retreat rates range from 0.33 to 1.96 mm yr^–1, and are among the highest rockwall retreat rates measured in arctic environments, presumably modulated by harsh environmental forcing on a porous sandstone rock cliff with a comparatively low rock mass strength. Here, we show the potential of laboratory‐calibrated ERT to provide accurate estimates of rockwall retreat rates even in ice‐rich permafrost talus slopes. Copyright © 2012 John Wiley & Sons, Ltd.