Mixed effects models were used to estimate within- and between-borehole variances.
The estimated variance components helped determining the optimal number of samples.
The vertical distribution of contaminants was assessed using mixed effects models.
Mixed effects models showed to have a superior fit compared to classical models.