Abstract. Potential impacts of global climate changes on the groundwater are largely unknown, especially for densely populated areas where groundwater is heavily exploited for public and industrial supply. Hence, to better plan and manage the groundwater resources, medium-long term numerical modelling of groundwater flow, which takes into account climate change, population growth, and industrial and agricultural activities, is needed. The objective of this paper is to tackle three main issues: (1) the development of a robust hydro-stratigraphic model of a multi-aquifer system from a well logs database by means of a novel multi-dimensional approach which includes a hierarchical classification of the lithologies, the interpretation of cross-sections, and the interpolation of aquifer boundary surfaces; (2) the parametrization and calibration of both a steady state and a transient groundwater flow model, starting from empirical relationships (for unconfined aquifer) and step-drawdown and well tests (for semi-confined and confined aquifers) to define equivalent homogenous sub-units; and (3) the simulation of steady state and transient scenarios based on projections about global climate change and variation in abstraction and recharge rates. These issues are illustrated for the Milan metropolitan area (Northern Italy) and the conterminous Po Plain portion. The results of the model allow to analyse the major components of the regional groundwater system (i.e. public supply wells withdrawals, discharge to gaining rivers and springs, recharge from irrigation networks and vegetated areas, flow transfer between aquifers). The groundwater level rising observed in the last decades caused serious problems in the urban areas and a progressive increase in the base-flow towards the gaining rivers. Simulations including effects of future climate scenarios (2017–2030) indicate a further increase in groundwater level in the next decades at a lower rate (ca. 0.3 m/year) with respect to that of the 1970–2016 period (ca. 1 m/year), due to the combined action of decreasing withdrawals and recharge.