Mapping soil moisture across the UK: assimilating cosmic-ray neutron sensors, remotely sensed indices, rainfall radar and catchment water balance data in a Bayesian hierarchical model

Abstract. Soil moisture is important in many hydrological and ecological processes. However, data sets which are currently available have issues with accuracy and resolution. To translate remotely sensed data to an absolute measure of soil moisture requires mapped estimates of soil hydrological properties and estimates of vegetation properties, and this introduces considerable uncertainty. We present an alternative methodology for producing daily maps of soil moisture over the UK at 2 km resolution (“SMUK”). The method is based on a simple linear statistical model, calibrated with 5 years of daily data from cosmic-ray neutron sensors at ∼ 40 sites across the country. The model is driven by precipitation, humidity, a remotely sensed “soil water index” satellite product and soil porosity. The spatial variation in the parameter describing the soil water retention (and thereby the response to precipitation) was estimated using daily water balance data from ∼ 1200 catchments with good coverage across the country. The model parameters were estimated by Bayesian calibration using a Markov chain–Monte Carlo method, so as to characterise the posterior uncertainty in the parameters and predictions. The approach reduces uncertainty by integrating multiple data sources, all of which have weaknesses but together act as a better constraint on the true soil moisture. The model explains around 70 % of the variance in the daily observations with a root-mean-square error of 0.05 m3 m−3, better than results from more complex process-based models. Given the high resolution of the inputs in time and space, the model can predict the very detailed variation in soil moisture which arises from the sporadic nature of precipitation events, including the small-scale and short-term variations associated with orographic and convective rainfall. Predictions over the period 2016 to 2023 demonstrated realistic patterns following the passage of weather fronts and prolonged droughts. The model has negligible computation time, and inputs and predictions are updated daily, lagging approximately 1 week behind real time.