HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via ISI Web of Science (3) Citing Articles via Google Scholar Google Scholar Articles by Bockgård, N. Articles by Niemi, A. Search for Related Content PubMed Articles by Bockgård, N. Articles by Niemi, A. GeoRef GeoRef Citation Agricola Articles by Bockgård, N. Articles by Niemi, A. Related Collections Effective Parameters Fractured Rock Variably Saturated Fluid Flow

SPECIAL SECTION: RESEARCH ADVANCES IN VADOSE ZONE HYDROLOGY THROUGH SIMULATIONS WITH THE TOUGH CODES

Role of Rock Heterogeneity on Lateral Diversion of Water Flow at the Soil–Rock Interface

Dep. of Earth Sciences, Uppsala Univ., Villav. 16, SE-752 36 Uppsala, Sweden
^* Corresponding author (niclas.bockgard{at}hyd.uu.se)

Received 22 August 2003.

Groundwater recharge is an important, yet often a highly uncertain upper boundary condition in groundwater flow models for fractured rock. The net infiltration at ground surface can often be estimated from water budgets, but the redistribution of water between a soil layer and the underlying bedrock is much more difficult to quantify. To address this question, we studied the role of lateral flow diversion at the soil–rock interface and how it influences the percentage of net infiltration that becomes inflow into the bedrock. Numerical experiments were performed with different net infiltration conditions, soil–rock hydraulic conductivity contrasts, hydraulic gradients, and most notably, rock heterogeneity. The rock was first represented as a deterministic homogeneous continuum, then as a stochastic heterogeneous continuum and analyzed with Monte Carlo simulations. The heterogeneous simulations predicted a large variation in the inflow into the rock between different stochastic realizations (from 50 to 100% of the net infiltration). Furthermore, at high hydraulic gradients in the bedrock, the mean flux into the rock from the heterogeneous simulations was smaller than the corresponding results from the homogeneous medium. The heterogeneous simulations showed formation of local unsaturated zones below the groundwater table. These occurred in cases with large hydraulic gradients at locations with large hydraulic conductivity contrasts, with low-conductive regions overlying high-conductive ones. This phenomenon was not captured with the homogenous models, which thereby may overestimate flux into fractured rock in cases where hydraulic gradients are large.