Abstract A numerical model accounting for variable density flow and transport was built up to quantify the actual and future (2050) salinization trends of a coastal aquifer in the Po delta (Northern Italy). SEAWAT 4.0 was employed to model the interaction between the surface drainage system and the underlying aquifer. PEST was employed for inverse parameter calibration using hydraulic heads and groundwater salinities as constraints. The calibrated model was used to predict the behavior of the coastal aquifer using a multiple scenario approach: increase in evapotranspiration induced by temperature increase; increase in the frequency of extreme high rainfall events; extreme drought conditions; and irrigation canals dewatering due to salinizat...