Given the inherent heterogeneity of subsurface formations, how can we accurately model and predict the spatial and temporal variability of groundwater flow and contaminant transport in fractured aquifers, where traditional Darcy-based flow models fail to account for anisotropy, preferential flow paths, and micro-scale interactions between fracture networks and porous media?