The problem of identifying and reproducing the hydrological behaviour of groundwater systems can often be set in terms of ordinary differential equations relating the inputs and outputs of their physical components under simplifying assumptions. Conceptual linear and nonlinear models described as ordinary differential equations are widely used in hydrology and can be found in several studies. Groundwater systems can be described conceptually as an interlinked reservoir model structured as a series of nonlinear tanks, so that the groundwater table can be schematized as the water level in one of the interconnected tanks. In this work, we propose a methodology for inferring the dynamics of a groundwater system response to rainfall, based on recorded time series data. The use of evolutionary techniques to infer differential equations from data in order to obtain their intrinsic phenomenological dynamics has been investigated recently by a few authors and is referred to as evolutionary modelling. A strategy named Evolutionary Polynomial Regression (EPR) has been applied to a real hydrogeological system, the shallow unconfined aquifer of Brindisi, southern Italy, for which 528 recorded monthly data over a 44-year period are available. The EPR returns a set of non-dominated models, as ordinary differential equations, reproducing the system dynamics. The choice of the representative model can be made both on the basis of its performance against a test data set and based on its incorporation of terms that actually entail physical meaning with respect to the conceptualization of the system. Key words groundwater; conceptual model; ordinary differential equations; evolutionary modelling; shallow aquifer InfĂ©rer la dynamique du systĂšme hydrogĂ©ologique Ă partir de sĂ©ries des donnĂ©es hydrologiques RĂ©sumĂ© Le problĂšme de l'identification et de la reproduction du comportement hydrologique des systĂšmes hydrogĂ©ologiques peut souvent ĂȘtre posĂ© en termes d'Ă©quations diffĂ©rentielles ordinaires relatives aux entrĂ©es et aux sorties de leurs composantes physiques, avec des hypothĂšses simplificatrices. Des modĂšles conceptuels linĂ©aires et non-linĂ©aires dĂ©crits sous forme d'Ă©quations diffĂ©rentielles ordinaires sont largement utilisĂ©s en hydrologie et peuvent ĂȘtre trouvĂ©s dans plusieurs Ă©tudes. Les systĂ©mes hydrogĂ©ologiques peuvent ĂȘtre dĂ©crits sur le plan conceptuel par un modĂšle de rĂ©servoirs interdĂ©pendants structurĂ© comme une sĂ©rie de rĂ©servoirs nonlinĂ©aires, de sorte que le niveau de la nappe peut ĂȘtre schĂ©matisĂ© comme Ă©tant le niveau d'eau dans lĂ©un des rĂ©servoirs interconnectĂ©s. Dans ce travail, nous proposons une mĂ©thodologie permettant d'infĂ©rer la dynamique de la rĂ©ponse d'un systĂ©me hydrogĂ©ologique aux prĂ©cipitations, sur la base de donnĂ©es temporelles enregistrĂ©es. L'utilisation de techniques Ă©volutives pour dĂ©duire les Ă©quations diffĂ©rentielles Ă partir de donnĂ©es afin d'obtenir leur dynamique phĂ©nomĂ©nologique intrinsĂ©que a Ă©tĂ© Ă©tudiĂ©e rĂ©cemment par quelques auteurs et est appelĂ©e modĂ©lisation Ă©volutive. Une stratĂ©gie...