Due to climate and population evolution, increased pressure is put on the groundwater resource, which calls for better understanding and models. In this paper, we describe a novel workflow to model the geological heterogeneity of coastal aquifers and apply it to the Roussillon plain (southern France). The main strength of the workflow is its capability to model aquifer heterogeneity when only sparse data are available while honoring the local geological trends and quantifying uncertainty.