A density-dependent numerical model was set up to quantify the actual and future (2050) salinization of a coastal aquifer in the Po Delta (Italy). SEAWAT 4.0 was used for this purpose, while PEST was used for automatic inverse parameter calibration. The use of high-resolution multi-level sampling allowed obtaining insights in the salinity stratification within the aquifer and on the vertical fluxes induced by the reclamation drainage network. The calibrated model was employed to build up a scenario investigating the effects of the projected relative sea level rise on groundwater salinity by 2050. This scenario allowed to identify the zones that resulted more affected by the relative sea level rise and to quantify the increase in salinizatio...