Hydrological Characteristics of Columnar Basalt Aquifers: Measuring and Modeling Skaftafellsheiði, Iceland

Basalt with columnar jointing can act as a good groundwater conductor. In areas with limited water resources in sedimentary rock, such as the Deccan Traps in India and the Columbia River basalt formations in Washington State (USA), large quantities of groundwater are abstracted from such basalt formations for drinking water supply and irrigation. The hydraulic properties of basaltic formations are difficult to quantify. To obtain a better understanding of their hydraulic properties, intensive field campaigns in Iceland were combined with a conceptual groundwater model in MODFLOW. The field experiments enabled us to derive the upper boundary conditions, like precipitation surplus, and obtain reliable ranges for the kh (0.01–0.3 m d−1) and kv (0.01–10 m d−1) of the basalt formations. The main objective was to test the concept of representative elementary volumes (REVs) for such basaltic regions. Precipitation excess for the Vestragil and Eystragil catchments was calculated by taking into account the orographic effect of precipitation. It was found that at higher elevations (600 m + msl) the precipitation was twice the amount compared to the base camp rain gauge at 100 m + msl. Calculated evapotranspiration (1–2 mm d−1) is in line with the literature. In the MODFLOW model, best results were obtained when the top layer (organic soil, peat, and regolith) was considered to be most conductive (up to 10 m d−1), with a gradual reduction in hydraulic conductivity with depth in the basaltic aquifers. This study shows that, when larger elementary volumes are used, a good model representation of basaltic regions can be created.