Les systèmes aquifères représentent une des principales ressources en eau potable actuellement, et sont directement connectés au cycle hydrologique dont ils constituent une partie souterraine. Dans le cadre du changement climatique est apparue la nécessité de développer un outil de modélisation numérique à même de rendre compte de la dynamique des systèmes aquifères de France. En Bretagne, les aquifères principaux sont sis sur un socle cristallin, peu profonds et très hétérogènes mais sont connectés directement aux rivières dont ils constituent le stock principal. Les données disponibles pour l’analyse de ces systèmes sont relativement faibles, nécessitant le développement d’une approche parcimonieuse mais à même de rendre compte de la complexité de la dynamique des systèmes aquifères situés dans les premiers mètres de roche altérée. L’approche développée ici est de modéliser les écoulements dans les versants qui sont agrégés dans la rivière dont les données de débits servent de point de comparaison à la calibration du modèle. Le modèle, couplé à des analyses de récession des débits en rivière, montre une bonne capacité à reproduire les processus d’écoulement des aquifères de socle, tout en donnant des indices sur le type de données nécessaires à la caractérisation de la physique de ces aquifères.