Interpretable modeling of spatiotemporal variations of groundwater total hardness across China