Abstract. The Newmarket Till forms a regionally extensive aquitard separating two major aquifer systems in the Greater Toronto area, Canada. The till is incised, and sometimes eroded entirely, by a network of sand- and gravel-filled channels forming productive aquifers and, locally, high-conductivity windows between aquifer systems. Leakage through the till may also be substantial in places. This study investigates the spatial variability of aquitard leakance in order to assess the relative importance of recharge processes to the lower aquifers. With a large database derived from water-well records and containing both hard and soft information, the Sequential Indicator Simulation method is used to generate maps of aquitard thickness and window probability. These can be used for targeting channel aquifers and for identifying potential areas of recharge to the lower aquifers. Conductivities are modeled from sparse data assuming that their correlation range is much smaller than the grid spacing. Block-scale leakances are obtained by upscaling nodal values based on simulated conductivity and thickness fields. Under the "aquifer-flow'' assumption, upscaling is performed by arithmetic spatial averaging. Histograms and maps of upscaled leakances show that heterogeneities associated with aquitard windows have the largest effect on regional groundwater flow patterns.
Rsum. La moraine glaciaire de Newmarket constitue un impermable d'extension rgionale sparant deux systmes aquifres dans la rgion du Grand Toronto (Canada). La moraine est entaille, et parfois entirement rode, par un rseau de chenaux combls de sables et de graviers formant des aquifres productifs et, localement, des fentres, zones forte conductivit hydraulique reliant les systmes aquifres. Une drainance au travers de la moraine peut galement tre significative par endroits. Cette tude s'intresse la variabilit spatiale de la drainance au travers de l'impermable, dans le but d'valuer l'importance relative des processus d'alimentation des aquifres infrieurs. partir d'une vaste base de donnes constitue par les mesures faites dans les puits et contenant des informations la fois certaines et incertaines, la mthode de simulation par indicateur squentiel est utilise pour crer des cartes d'paisseur de l'impermable et de probabilit d'existence des fentres. Ces cartes peuvent tre utilises pour mettre en vidence les aquifres de chenaux et pour identifier les zones potentielles de recharge des aquifres infrieurs. Les conductivits hydrauliques sont modlises partir de donnes clairsemes en supposant que leur gamme de corrlation est beaucoup plus faible que le pas de la grille. La drainance l'chelle des blocs est obtenue par accroissement du niveau d'chelle des valeurs nodales bas sur les champs de conductivit simule et des paisseurs. partir de l'hypothse d'coulement dans l'aquifre, l'accroissement d'chelle est ralis en faisant une moyenne arithmtique spatiale. Des histogrammes et des cartes de drainance avec accroissement d'chelle montrent que les htrognits associes aux fentres dans l'impermable constituent l'effet le plus important dans l'organisation des coulements souterrains rgionaux.
Resumen. El Till de Newmarket forma un acuitardo de extensin regional que separa dos sistemas acuferos principales en la zona de Greater Toronto (Canad). El till est horadado, y a veces completamente erosionado, por una red de canales rellenos de arena y grava. Estos constituyen acuferos productivos y, localmente, conexiones de alta permeabilidad entre sistemas acuferos. El goteo a travs del till puede ser fundamental en ciertos lugares. Este estudio investiga la variabilidad espacial del goteo desde el acuitardo con el fin de establecer la importancia relativa de los procesos de recarga hacia los acuferos inferiores. Se ha utilizado el mtodo de la Simulacin Indicadora Secuencial, soportado por una gran base de datos que contiene registros de sondeos e incluye informacin blanda y dura, para generar mapas de espesor del acuitardo y mapas de probabilidad de discontinuidades laterales en el acuitardo. Estos resultados se pueden emplear para identificar canales y las zonas de recarga potencial a los acuferos inferiores. Se han modelado las conductividades a partir de datos dispersos, suponiendo que la distancia de correlacin es mucho menor que el espaciado de la malla. Se han calculado valores de goteo a escala de bloque por medio del escalado de los valores nodales, basndose en los campos simulados de conductividad y espesor. El escalado se calcula mediante un promedio aritmtico espacial, bajo la hiptesis de "flujo en el acufero". Los histogramas y mapas de goteos escalados muestran que las heterogeneidades asociadas a las discontinuidades en el acuitardo producen el efecto ms notable en la distribucin regional del flujo de aguas subterrneas.