In this research, univariate and bivariate statistical methods were applied to rainfall, river and piezometric level datasets belonging to 24 years long time series (1986-2009). These methods, that often are used to understand the effects of precipitation on rivers and karstic springs discharge, have been used to assess, piezometric level response to rainfall and river level fluctuations in a porous aquifer. A rain gauge, a river level gauge and three wells, located in Central Italy along the lower Pescara river valley in correspondence of its important alluvial aquifer, provided the data. The statistical analysis has been used within a known hydrogeological framework, which has been refined by mean of a photo-interpretation and a GPS survey. Water-groundwater relationships were identified following the autocorrelation and cross-correlation analyses; the spectral analysis and mono-fractal features of time series were assessed, in order to provide information on multy-year variability, data distributions, their fractal dimension and the distribution return time within the historical time series. The statistical-mathematical results were interpreted through field work that identified distinct groundwater flowpaths within the aquifer and enabled the implementation of a conceptual model, improving the knowledge on water resources management tools.