An efficient and robust soil moisture (SM) sampling scheme that can capture the spatial variability of SM is required for the accurate calibration and validation of satellite-based SM retrievals. Often, this process requires numerous sampling points, consuming a significant amount of time. Therefore, it is crucial to develop efficient sampling methods for the improvement of satellite-based SM estimations. The objectives of this study were to define an efficient sampling strategy that could be beneficial for the validation of satellite SM estimations; investigate the role of RS covariates in developing such a strategy; and evaluate the performance of the new sampling scheme over various spatial and temporal domains. In this study, we used the conditioned Latin hypercube sampling (cLHS) algorithm to define an efficient sampling strategy. To this end, remote sensing (RS) raster and digital elevation models (DEM) were used to identify numerous environmental covariates to locate sampling points for characterizing spatial variability of SM at the agricultural field scale. A random forest-based technique, the Boruta algorithm, was also applied to select the most important covariates for utilization into the cLHS algorithm. We used the statistical moments (mean and standard deviation, SD) of the field to select the efficient sample size that can best represent SM status in the field. To evaluate the new sampling scheme, a second data set obtained during a different month for the same agricultural field was used. However, because of the potential for high spatial and temporal correlations between training and test covariates when obtained for the same region, we also used different test datasets in New Zealand to evaluate the sampling scheme. Results showed that the RS covariates obtained from SAR and optical imagery were among the most significant covariates for capturing the spatial variability of SM even if they were not acquired on the day of collection. Also, the new sampling scheme could capture the SM spatial pattern of the field for both test datasets with RMSE less than 4% volumetric SM, which is within the range of the expected performance for most satellite SM products. The evaluation of the new sampling scheme on the New Zealand datasets confirmed the functionality of the proposed sampling scheme for a different temporal and spatial domain.

• Most conceptual bucket models have an upper limit on simulated soil moisture deficit. • Problems arise when the bucket “empties” because ET drops to unrealistic (low) levels. • Alternatives include bottomless buckets or deficit-based soil moisture accounting. • Here, we switch to a deficit-based scheme while keeping everything else constant. • Tested over historic drought, model performance and realism are enhanced. Rainfall-runoff models based on conceptual “buckets” are frequently used in climate change impact studies to provide runoff projections. When these buckets approach empty, the simulated evapotranspiration approaches zero, which places an implicit limit on the soil moisture deficit that can accrue within the model. Such models may cease to properly track the moisture deficit accumulating in reality as dry conditions continue, leading to overestimation of subsequent runoff and possible long-term bias under drying climate. Here, we suggest that model realism may be improved through alternatives which remove the upper limit on simulated soil moisture deficit, such as “bottomless” buckets or deficit-based soil moisture accounting. While some existing models incorporate such measures, no study until now has systematically assessed their impact on model realism under drying climate. Here, we alter a common bucket model by changing the soil moisture storage to a deficit accounting system in such a way as to remove the upper limit on simulated soil moisture deficit. Tested on 38 Australian catchments, the altered model is better able to track the decline in soil moisture at the end of seasonal dry periods, which leads to superior performance over varied historic climate, including the 13-year “Millennium” drought. However, groundwater and GRACE data reveal long-term trends that are not matched in simulations, indicating that further changes may be required. Nonetheless, the results suggest that a broader adoption of bottomless buckets and/or deficit accounting within conceptual rainfall runoff models may improve the realism of runoff projections under drying climate.

• Most conceptual bucket models have an upper limit on simulated soil moisture deficit. • Problems arise when the bucket “empties” because ET drops to unrealistic (low) levels. • Alternatives include bottomless buckets or deficit-based soil moisture accounting. • Here, we switch to a deficit-based scheme while keeping everything else constant. • Tested over historic drought, model performance and realism are enhanced. Rainfall-runoff models based on conceptual “buckets” are frequently used in climate change impact studies to provide runoff projections. When these buckets approach empty, the simulated evapotranspiration approaches zero, which places an implicit limit on the soil moisture deficit that can accrue within the model. Such models may cease to properly track the moisture deficit accumulating in reality as dry conditions continue, leading to overestimation of subsequent runoff and possible long-term bias under drying climate. Here, we suggest that model realism may be improved through alternatives which remove the upper limit on simulated soil moisture deficit, such as “bottomless” buckets or deficit-based soil moisture accounting. While some existing models incorporate such measures, no study until now has systematically assessed their impact on model realism under drying climate. Here, we alter a common bucket model by changing the soil moisture storage to a deficit accounting system in such a way as to remove the upper limit on simulated soil moisture deficit. Tested on 38 Australian catchments, the altered model is better able to track the decline in soil moisture at the end of seasonal dry periods, which leads to superior performance over varied historic climate, including the 13-year “Millennium” drought. However, groundwater and GRACE data reveal long-term trends that are not matched in simulations, indicating that further changes may be required. Nonetheless, the results suggest that a broader adoption of bottomless buckets and/or deficit accounting within conceptual rainfall runoff models may improve the realism of runoff projections under drying climate.

Evidence suggests that catchment state variables such as groundwater can exhibit multiyear trends. This means that their state may reflect not only recent climatic conditions but also climatic conditions in past years or even decades. Here we demonstrate that five commonly used conceptual “bucket” rainfall‐runoff models are unable to replicate multiyear trends exhibited by natural systems during the “Millennium Drought” in south‐east Australia. This causes an inability to extrapolate to different climatic conditions, leading to poor performance in split sample tests. Simulations are examined from five models applied in 38 catchments, then compared with groundwater data from 19 bores and Gravity Recovery and Climate Experiment data for two geographic regions. Whereas the groundwater and Gravity Recovery and Climate Experiment data decrease from high to low values gradually over the duration of the 13‐year drought, the model storages go from high to low values in a typical seasonal cycle. This is particularly the case in the drier, flatter catchments. Once the drought begins, there is little room for decline in the simulated storage, because the model “buckets” are already “emptying” on a seasonal basis. Since the effects of sustained dry conditions cannot accumulate within these models, we argue that they should not be used for runoff projections in a drying climate. Further research is required to (a) improve conceptual rainfall‐runoff models, (b) better understand circumstances in which multiyear trends in state variables occur, and (c) investigate links between these multiyear trends and changes in rainfall‐runoff relationships in the context of a changing climate.