Conceptual model of the geometry and physics of water flow in a fractured basalt vadose zone

A conceptual model of the geometry and physics of water flow in a fractured basalt vadose zone was developed based on the results of lithological studies and a series of ponded infiltration tests conducted at the Box Canyon site near the Idaho National Engineering and Environmental Laboratory. The infiltration tests included one 2‐week test in 1996, three 2‐day tests in 1997, and one 4‐day test in 1997. For the various tests, initial infiltration rates ranged from 4.1 cm/d (4.75 ×10−7 m/s) to 17.7 cm/d (2.05×10−7 m/s) and then decreased with time, presumably because of mechanical or microbiological clogging of fractures and esicular basalt in the near‐surface zone, as well as the effect of entrapped air. The subsurface moisture redistribution was monitored with tensiometers, neutron logging, time domain reflectrometry, and ground‐penetrating radar. A conservative tracer, potassium bromide, was added to the pond water at a concentration of 3 g/L to monitor water flow with electrical resistivity probes and water sampling. Analysis of the data shows evidence of preferential flow rather than the propagation of a uniform wetting front. We propose a conceptual model describing the saturation‐desaturation behavior of the basalt, in which rapid preferential flow occurs through the largest vertical fractures, followed by a gradual wetting of other fractures and the basalt matrix. Fractures that are saturated early in the tests may become desaturated thereafter, which we attribute to the redistribution of water between fractures and matrix. Lateral movement of water takes place within horizontal fracture and rubble zones, enabling development of perched water bodies.