Adaptation of ecosystems’ root zones to climate change critically affects drought resilience and vegetation productivity. However, a global quantitative assessment of this mechanism is missing. In this study, we analyzed high-quality observation-based data to find that the global average root zone water storage capacity (S_R) increased by 11%, from 182 to 202 mm in 1982–2020. The total increase of S_R equals to 1652 billion m³ over the past four decades. S_R increased in 9 out of 12 land cover types, while three relatively dry types experienced decreasing trends, potentially suggesting the crossing of ecosystems’ tipping points. Our results underscore the importance of accounting for root zone dynamics under climate change to assess drought impacts.

Assimilating remote sensing observations into land surface models has become common practice to improve the accuracy of terrestrial water storage (TWS) estimates such as soil moisture and groundwater, for understanding the land surface interaction with the climate system, as well as assessing regional and global water resources. Such remote sensing observations include soil moisture information from the L-band Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions, and TWS information from the Gravity Recovery And Climate Experiment (GRACE). This study evaluates the benefit of assimilating them into the Community Atmosphere and Biosphere Land Exchange (CABLE) land surface model. The evaluation is conducted in the Goulburn River catchment, South-East Australia, where various in situ soil moisture and groundwater level data are available for validating data assimilation (DA) approaches. It is found that the performance of DA mainly depends on the type of observations that are assimilated. The SMOS/SMAP-only assimilation (SM DA) improves the top soil moisture but degrades the groundwater storage estimates, whereas the GRACE-only assimilation (GRACE DA) improves only the groundwater component. Assimilating both observations (multivariate DA) results in increased accuracy of both soil moisture and groundwater storage estimates. These findings demonstrate the added value of multivariate DA for simultaneously improving different model states, thus leading to a more robust DA system.

The spatial heterogeneity and temporal variation of soil moisture and surface heat fluxes are key to many geophysical and environmental studies. It has been demonstrated that they can be mapped by assimilating soil thermal and wetness information into surface energy balance models. The aim of this work is to determine whether enhancing the spatial resolution or temporal sampling frequency of soil moisture data could improve soil moisture or surface heat flux estimates. Two experiments are conducted in an area mainly covered by grassland, and land surface temperature (LST) observations from the Geostationary Operational Environmental Satellite (GOES) mission are assimilated together with either an enhanced L-band passive soil moisture product (9 km, 2-3 days) from the Soil Moisture Active Passive (SMAP) mission or a merged product (36 km, quasi-daily) from the SMAP and the Soil Moisture Ocean Salinity (SMOS) mission. The results suggest that the availability of soil moisture observations is increased by 41% after merging data from the SMAP and the SMOS missions. A comparison with results from a previous study that assimilated a coarser SMAP soil moisture product (36 km, 2-3 days) suggests that enhancing the temporal sampling frequency of soil moisture observations leads to improved soil moisture estimates at both the surface and root zone, and the largest improvement is seen in the bias metric (0.008 and 0.007m 3 m -3 on average at the surface and root zone, respectively). Enhancing the spatial resolution, however, does not significantly improve soil moisture estimates, particularly at the surface. Surface heat flux estimates from assimilating soil moisture data of different spatial or temporal resolutions are very similar.