We propose a thermo-fluid-dynamics model to study some recent uplift episodes occurred in the period 2008–2013 at Campi Flegrei caldera (Italy). Accounting for eight overpressure sources (from 5 to 40 MPa) in the hydrothermal system, our model solves for heat and momentum balances to obtain fluid velocities responsible for the observed ground displacement. For a validation of the model we use a dataset from seven continuous GNSS stations of the Neapolitan Volcanoes Continuous GPS network (NeVoCGPS), belonging to a geodetic monitoring system covering the Neapolitan volcanic area, and operated by the Istituto Nazionale di Geofisica e Vulcanologia. We compare the observed and modelled vertical displacements to assess “threshold” values for the...