Application of analytical diffusion models to outcrop observations: Implications for mass transport by fluid flow through fractures

A pavement outcrop with excellent exposure of spatial relationships among joints, veins, small offset normal faults, and associated alteration halos (redox fronts) provided an opportunity to compare predictions of analytical models for reaction front propagation in a fracture-matrix system with a real-field situation. The results have important implications for fluid flow and pollutant transport through a fractured medium. The alteration halos observed suggest that all joints of different sets and most small faults are conductive to meteoric water at shallow depth. On the other hand, veins are local barriers to mass transport by diffusion. By using petrologic and petrophysical data, analytical modeling, and the width of the alteration halos, it was possible to estimate when the fracture network was open to fluid flow. The inferred time for fluid flow and diffusion through the fracture network is sensitive to the porosity n of the rock matrix used in the analytical solutions: 2200 ± 500 years with n = 0.08, 4600 ± 900 years with n = 0.05, and 16,000 ± 4000 with n = 0.02. The second and third age determinations are consistent with the landscape evolution of the area since the end of the last Wűrmian ice age and with the timing required to fill the fractures observed in outcrop. We suggest that analytical modeling is an important tool for the determination of transport and reaction time scales in fractured formations where it is constrained by a robust petrophysical and chemical properties data set.