<p>Knowledge of the variances and covariances of gridded terrestrial water storage anomalies (TWS) as observed with GRACE and GRACE-FO is crucial for many applications thereof. For example, data assimilation into different models, trend estimations, or combinations with other data set require reliable estimations of the variances and covariances. Today, the Level-2 Stokes coefficients are provided with formal variance-covariance matrices which can yield variance-covariance matrices of the gridded data after a labourious variance propagation through all post-processing steps, including filtering and spherical harmonic synthesis. Unfortunately, this is beyond the capabilities of many, if not most, users.</p><p><br>This is why, we developed a spatial covariance model for gridded TWS data. The covariance model results in non-homogeneous, non-stationary, and anisotropic covariances. This model also accommodates a wave-like behaviour in latitudinal-directed correlations caused by residual striping errors. The model is applied to both VDK3 filtered GFZ RL06 and ITSG-Grace2018 TWS data.&#160;</p><p><br>With thus derived covariances it is possible to estimate the uncertainties of mean TWS time series for any arbitrary region such as river basins. On the other hand, such time series uncertainties can also be derived from the afore mentioned formal covariance matrices. Here, only the formal covariance matrices of ITSG-Grace2018 are used which are also filtered with the VDK3 filter. All together, we are able to compare globally the time series uncertainties of both the modelled and formal approach. Further, the modelled uncertainties are compared to empirical standard deviations in arid regions in the Arabian, Sahara, and Gobi desert where residual hydrological signal can be neglected. Both in the temporal and spatial domain they show a very satisfying agreement proving the usefulness of the covariance model for the users.&#160;</p>
Abstract We evaluate trends in terrestrial water storage over 1950–2100 in CMIP6 climate models against a new global reanalysis from assimilating GRACE and GRACE-FO satellite observations into a hydrological model. To account for different timescales in our analysis, we select regions in which the influence of interannual variability is relatively small and observed trends are assumed to be representative of the development over longer periods. Our results reveal distinct biases in drying and wetting trends in CMIP6 models for several world regions. Specifically, we see high model consensus for drying in the Amazon, which disagrees with the observed wetting. Other regions show a high consensus of models and observations suggesting qualitatively correctly simulated trends, e.g., for the Mediterranean and parts of Central Africa. A high model agreement might therefore falsely indicate a robust trend in water storage if it is not assessed in light of the observed developments. This underlines the potential use of maintaining an adequate observational capacity of water storage for climate change assessments.