Effects of sedimentary heterogeneity on groundwater flow in a Quaternary pro-glacial delta environment: joining facies analysis and numerical modelling

We present the sedimentological survey of Quaternary sediments developed in glacio-fluvial to glacio-lacustrine delta environments and exposed at an a bandoned quarry in Northern Italy. This outcrop is considered an analogue o f sedimentary structures that could form a real aquifer. We recognise the f ollowing sedimentary units: (I) colluvial sediments with human artefacts (l ess than 2 m thick); (II) gravelly sands with oversized pebbles (2.5-5.5 m thick), bounded at the base by a roughly flat composite disconformity surfa ce; (m) and (TV) gravelly-sand foresets which form two bodies with differen t dips that are separated by an inclined, almost planar, erosional surface. From past observations in the surrounding area and a geoelectrical survey, we estimate that the minimum total thickness of the delta system is 20 m. We have analysed the sedimentary characteristics of two well exposed areas, which are located within units II, III and IV and covering areas of 115 an d 95 m(2), using detailed sedimentological surveys along six vertical secti ons. With the aid of a topographic survey we have reconstructed the distrib ution of textural units based on sedimentary facies and including informati on about grain-size distribution, sorting and packing. We have constructed a sedimentary model for the two areas using orthogonal cells of 10 x 10 cm, attributing the dominant textural unit to each cell. The distribution of s edimentary textural units is used to estimate the distribution of hydraulic conductivity for the two areas of this aquifer analogue; we have determine d the hydraulic conductivity by assigning different values to the different textural units according to laboratory measurements of porosity and grain- size distribution, Kozeny-Carman's equation and literature data. We have mo delled groundwater how in the two areas approximating the balance equation with finite differences, using a regular grid with 10 cm spacing, and we ha ve performed some numerical experiments. First, we have evaluated the equiv alent conductivity tensor for the two areas. Secondly we have considered a coarse grid, with 1.5 m spacing, and at this scale we have analysed the val idation of the discrete form of Darcy's law usually introduced in finite di fference modelling. The results of the numerical experiments show the aniso tropic behaviour of the medium, especially for unit II composed of sandy be ds, some of which lead to a preferential path for groundwater flow in the h orizontal direction. The deeper sedimentary units m and TV are characterise d by a lower anisotropy ratio for hydraulic conductivity, because the layer ed structures of these units have different dips.