A regional-scale conceptual and numerical groundwater flow model in fluvio-glacial sediments for the Milan Metropolitan area (Northern Italy)

Study region: The Milan metropolitan area lies on one of the most important aquifer in Italy, heavily exploited for public and industrial water supply. The area, covering 3135 km2 in the Po Plain (Northern Italy) with a continental climate, is bounded by the Po, the Adda and the Ticino rivers and by the prealpine foothills. Regional hydrology is characterised by a network of natural and man-made elements, and lowland springs. The sedimentary sequence, from bottom to top, is formed by meandering river plain deposits, the distal fringe of the glacial outwash plains and proximal braid-plain deposits. Study focus: This study proposes a general approach for aquifer geometry reconstruction and hydrodynamic parametrization of hydrofacies in fluvio-glacial deposits, and their implementation into a 3D regional groundwater flow model. This approach is based on sedimentologically-defined lithofacies/hydrofacies and their correlation in space to obtain nearly homogeneous subunits starting from available data (i.e. 8628 borehole logs, grain size distributions, well tests) and sedimentological knowledge. New hydrological insights for the region: The calibrated 3D FEM groundwater model allows quantifying the main components of the hydrogeological budget at the regional scale, and the fluxes among the different hydro-stratigraphic units. A sensitivity analysis of groundwater levels to the main recharge components suggests importance of anthropogenic disturbances with respect to natural recharge, and that land-use change may impact water resources more than climate change.