19Nitrate transport in heterogeneous bedrock aquifers is influenced by mechanisms that operate at 20 different spatial and temporal scales. To understand these mechanisms in a fractured sandstone 21 aquifer with high porosity, a groundwater-flow and nitrate transport model -reproducing multiple 22 2010). In these studies, parameters are varied and the amount of nitrate leached below the root 59 zone is monitored and modelled to assess the impact in term of groundwater nitrate contamination. 60The transferability of small-scale results at watershed scale is however rather limited. Distributed 61 physically-based hydrological models are thus commonly used for a more comprehensive 62 4 understanding of nitrate contamination at the watershed scale (e.g., Vaché et al., 2002; Hattermann 63 et al., 2005; Bracmort et al., 2006;Rode et al., 2009). The representation of aquifers in hydrological 64 models is generally over-simplified, with flow and transport processes modelled through simple 65 transfer functions, which for complex aquifer systems may not be adequate to simulate the spatial 66 and temporal distribution of nitrate contamination in groundwater and its transport to downstream 67 surface-water receivers. For this key task, understanding aquifer dynamics, especially groundwater 68 residence times, is crucial to support agricultural and water management plans, given that the 69 contamination intensity and time needed for water quality to change as a result of BMPs both 70 depend to a great extent on hydrogeological characteristics. Therefore, many researchers have 71 taken into account the dynamics of aquifers using numerical groundwater models (e.g., Lunn et al., 72 1996;Vinten and Dunn, 2001;Molénat et al., 2002;Wriedt and Rode, 2006; Dimitriou et al., 2010; 73 Paradis et al., 2016). Numerical modeling of groundwater flow and nitrate transport in aquifers 74 can provide quantitative insight into water quality sustainability as long as model parameters (e.g., 75hydraulic property values and structures) and inputs (e.g., groundwater recharge, nitrate flux) are 76 reasonably well established (Michael and Voss 2009). Estimation of model parameters and inputs 77 can be difficult, notably for large aquifers with limited data. To narrow the range of model 78 uncertainties, model calibration with multiple data sets is generally proposed, with each data set 79 providing different and complementary constraints on the aquifer model (e.g., Kim et al., 1999). 80
81In this paper, the control exerted by the dynamics of a sandstone aquifer on groundwater 82 contamination by nitrate, in a watershed with intensive agriculture, is examined using a 83 groundwater flow and nitrate transport model. The heterogeneous and dual-porosity nature of the 84 aquifer revealed by field and laboratory data were expected to influence water quality, but the 85 5 magnitude and trends of the contamination still remained unknown. To address this shortcoming, 86 the nitrate transport mechanisms within the aquifer were inferred from the grou...