This dataset includes the results of the application of the SHYFEM numerical model to the Po River-Delta-Sea (RDS) system (Italy). The hydrodynamics and the interaction between riverine and marine waters in the system were studied through the application of the unstructured finite elements 3D numerical model to a domain comprising the Po river branches, 7 coastal lagoons and the shelf sea. The vertical structure of the model is z-levels, with 1 m thickness up to 10 m depth, 2 m thickness up to 26 and finally 4 m thickness to 34 m depth. The model was forced with real data of water levels, meteorological parameters and freshwater inputs. The simulations included the baroclinic gradients in the hydrodynamics calculations induced by mixing of riverine and sea waters, and by heat fluxes with the atmosphere.Two simulations were run, one for the model calibration (March - September 2009) and a two year long simulation (2010-2011) for the model validation and analysis. The SHYFEM parameter files of the two simulations are included in the data archive.The dataset includes the time-averaged hydrodynamic and T/S fields of the two simulations and of a relevant river flood in May 2009. For the 2010-2011 simulation we also provide the Water Renewal Time distribution and its range of variability.The model values calculated on the triangular finite elements grid (scalar values on the nodes and velocities on the elements) were converted into a regular grid netcdf format with horizontal spacing of 100 m. The reference coordinate system is EPSG 3003.

Supplement to: Maicu, Francesco; De Pascalis, Francesca; Ferrarin, Christian; Umgiesser, Georg (2018): Hydrodynamics of the Po River‐Delta‐Sea system. Journal of Geophysical Research: Oceans, 123