Modeling groundwater flow in bedrock can be particularly challenging due to heterogeneities associated with fracture zones. However, fracture zones can be difficult to map, particularly in forested areas where tree cover obscures land surface features. This study presents the evidence of fracture zones in a small, snowmelt-dominated mountain headwater catchment and explores the significance of these fracture zones on groundwater flow in the catchment. A newly acquired bare earth image acquired using LiDAR identifies a previously undetected linear erosion zone that passes near a deep bedrock well at low elevation in the catchment. Borehole geophysical logs indicate more intense fracturing in this well compared to two wells at higher elevation. The well also exhibited a linear flow response during a pumping test, which is interpreted to reflect the influence of a nearby vertical fracture zone. The major ion chemistry and stable isotope composition reveal a slightly different chemical composition and a more depleted isotopic signature for this well compared to other groundwaters and surface waters sampled throughout the catchment. With this evidence of fracturing at the well scale, an integrated land surface – subsurface hydrologic model is used to explore four different model structures at the catchment scale. The model is refined in steps, beginning with a single homogeneous bedrock layer, and progressively adding 1) a network of large-scale fracture zones within the bedrock, 2) a weathered bedrock zone, and 3) an updated LiDAR-derived digital elevation model, to gain insight into how increasing subsurface geological complexity and land surface topography influence model fit to observed data and the various water balance components. Ultimately, all of the models are considered plausible, with similar overall fit to observed data (snow, streamflow, pressure heads in piezometers, and groundwater levels) and water balance results. However, the models with fracture zones and a weathered zone had better fits for the low elevation well. These models contributed slightly more baseflow (~14% of streamflow) compared to models without a weathered zone (~1%). Thus, in the watershed scale model, including a weathered bedrock zone appears to more strongly influence the hydrology than only including fracture zones.
Member of collection