Solute transport characteristics of columnar volumetric contraction networks with mega column structure and aperture variability

Numerical simulations explore for the first time the role of mega columns and aperture variability on particle transport through mature volumetric contraction networks as informed by a unique synthesis of network propagation and maturity. Columnar fracture patterns are generated by updating a series of Voronoi centers to the midpoint of a generated polygon over many iterations, creating 250 network realizations. A DFN simulator solves for fluid flow and tracks conservative particle trajectories within each network. Dominant fracture attributes impacting fluid flow and solute transport in volumetric contraction networks are fracture orientation, density, and aperture/transmissivity. Ensemble plume snapshots generated by networks with equal fracture transmissivity define a baseline-level of dispersion that is solely attributed to network structure and connectivity. Longitudinal and transverse dispersion increase and the center of plume mass becomes delayed relative to the baseline case when fracture transmissivity is varied according to a lognormal distribution. The incorporation of highly-transmissive, large-aperture mega column fractures leads to plume snapshots with a more pronounced leading edge and an order of magnitude faster average breakthrough times. The breakthrough curves contain three peaks reflecting contrasting transport pathways in which particles are: (i) initially placed in mega column fractures and remain in these features until exiting the model domain, (ii) initially placed into small column fractures, incur additional time to migrate and enter a mega column fracture, and remain within those mega columns, and (iii) initially placed in small column fractures and remain in these fractures. Incorporating variability in fracture transmissivity for both small column and mega column fractures disrupts the binary distinction between small column and mega column fracture velocities and leads to dispersed breakthroughs over long time scales with a single peak. These results demonstrate that preferential flow paths emerge in volumetric contraction networks due to contracts in fracture transmissivity, not fracture connectivity as observed in tectonic networks.