This work deals with the PNRR project “campionamento delle matrici ambientali”, whitin we particularly focused on soil and surface water. For these matrices we studied sampling method and procedure for pH analysis. For soil matrix, the sampling uncertainty was estimated using the ANOVA method.