Soil moisture is key for quantifying soil-atmosphere interactions and the ESA-CCI (European Space Agency-Climate Change Initiative) provides historical (>30 years) satellite soil moisture global grids with spatial resolution of ~27km. This dataset is incomplete (contains gaps) due to conditions such as dense vegetation or extremely dry surfaces. Here we provide a framework to increase the spatial resolution and fill gaps (reporting associated uncertainty) of the ESA-CCI (v4.5) soil moisture dataset. The outcome is a new dataset of gap-free global mean annual soil moisture and uncertainty
for 28 years (1991-2018) across 15km grids. We compare the performance of machine learning odels using only terrain parameters (e.g., slope, wetness index) against predictions using terrain parameters, bioclimatic information, and soil type classes. We use independent field information from the International Soil Moisture Network (ISMN, n=13376) and in-situ precipitation records (n=171) only for model evaluation purposes. Using only terrain parameters to predict soil moisture results in a parsimonious approach comparable with a more complex model that includes additional bioclimatic and soil information. The correlation between observed and predicted soil moisture values varies from r=0.69 to r=0.87 with root mean squared errors (RMSE) around 0.03 and 0.04 m3/m3. Our soil moisture predictions improve: (a) the correlation with the ISMN (when compared with the original ESA-CCI product) from r=0.30 (RMSE=0.09 m3/m3 ) to r=0.66 (RMSE=0.05 m3/m3 ); and (b) the
correlation with local precipitation records across boreal (from r=<0.3 up r=0.49) or tropical areas (from r=<0.3 to r=0.46) which are currently poorly represented in the ISMN. Temporal trends show a decline of global annual soil moisture using: a) data from the ISMN (-1.5 [-1.8, -1.24]%, b) associated locations from the original ESA-CCI dataset (- 0.87[-1.54, -0.17]%), c) associated locations from predictions based on terrain parameters (-0.85[-1.01, -0.49]%), and d)associated locations from predictions including bioclimatic and soil type classes (-0.68[-0.91, -0.45]%). Our parsimonious downscaled soil moisture predictions are independent of climate variables and vegetation indexes, to avoid potential spurious correlations in future research, and they complement information about soil moisture dynamics worldwide.