[1] Diffusive mass transfer in rock fractures is strongly affected by fluid flow in addition to material properties. The flow-dependence of matrix diffusion is quantified by a random variable (''transport resistance'') denoted as b [T/L] and computed from the flow field by following advection trajectories. The numerical methodology for simulating fluid flow is mesh-free, using Fup basis functions. A generic statistical model is used for the transmissivity field, featuring three correlation structures: (i) highly connected non-multiGaussian; (ii) poorly connected (or disconnected) non-multi-Gaussian ; and (iii) multiGaussian. The moments of b are shown to be linear with distance, irrespective of the structure, after approximately 10 integral scales of ln T. Percentiles of b are found to be linear with the mean b when considering all three structures. Taking advantage of this property, a potentially useful relationship is presented between b percentiles and the fracture mean water residence time that integrates all structures with high variability; it can be used in discrete fracture network simulations where T statistical data on individual fractures are not available.Citation: Cvetkovic, V., and H. Gotovac (2013), Flow-dependence of matrix diffusion in highly heterogeneous rock fractures, Water Resour. Res., 49,[7587][7588][7589][7590][7591][7592][7593][7594][7595][7596][7597]