Global land water underpins livelihoods, socioeconomic development, and ecosystems. It remains unclear how water availability has changed in recent decades. Using an ensemble of observations, we quantified global land water availability over the past two decades. We show that the Southern Hemisphere has dominated the declining trend in global water availability from 2001 to 2020. The significant decrease occurs mainly in South America, southwestern Africa, and northwestern Australia. In the Northern Hemisphere, the complex regional increasing and decreasing trends cancel each other, resulting in a negligible hemispheric trend. The variability and trend in water availability in the Southern Hemisphere are largely driven by precipitation associated with climate modes, particularly the El Niño-Southern Oscillation. This study highlights their dominant role in controlling global water availability.

Surface downward radiation (SDR), including shortwave downward radiation (SWDR) and longwave downward radiation (LWDR), is of great importance to energy and climate studies. Considering the lack of reliable SDR data with a high spatiotemporal resolution in the East Asia–Pacific (EAP) region, we derived SWDR and LWDR at 10-min and 0.05° resolutions for this region from 2016 to 2020 based on the next-generation geostationary satellite Himawari-8 (H-8). The SDR product is unique in terms of its all-sky features, high accuracy, and high-resolution levels. The cloud effect is fully considered in the SDR product, and the influence of high aerosol loadings and topography on the SWDR are considered. Compared to benchmark products of the radiation, such as Clouds and the Earth’s Radiant Energy System (CERES) and the European Centre for Medium-Range Weather Forecasts (ECMWF) next-generation reanalysis (ERA5), and the Global Land Surface Satellite (GLASS), not only is the resolution of the new SDR product notably much higher, but the product accuracy is also higher than that of those products. In particular, hourly and daily root-mean-square errors of the new SWDR are 104.9 and 31.5 W m⁻², respectively, which are much smaller than those of CERES (at 121.6 and 38.6 W m⁻², respectively), ERA5 (at 176.6 and 39.5 W m⁻², respectively), and GLASS (daily of 36.5 W m⁻²). Meanwhile, RMSEs of hourly and daily values of the new LWDR are 19.6 and 14.4 W m⁻², respectively, which are comparable to that of CERES and ERA5, and even better over high-altitude regions.

Land surface temperature (LST) plays a critical role in land surface processes. However, as one of the effective means for obtaining global LST observations, remote sensing observations are inherently affected by cloud cover, resulting in varying degrees of missing data in satellite-derived LST products. Here, we propose a solution. First, the data interpolating empirical orthogonal functions (DINEOF) method is used to reconstruct invalid LSTs in cloud-contaminated areas into ideal, clear-sky LSTs. Then, a cumulative distribution function (CDF) matching-based method is developed to correct the ideal, clear-sky LSTs to the real LSTs. Experimental results prove that this method can effectively reconstruct missing LST data and guarantee acceptable accuracy in most regions of the world, with RMSEs of 1–2 K and R values of 0.820–0.996 under ideal, clear-sky conditions and RMSEs of 4–7 K and R values of 0.811–0.933 under all weather conditions. Finally, a spatiotemporally continuous MODIS LST dataset at 0.05° latitude/longitude grids is produced based on the above method.

The development and continuity of consistent long-term data records from similar overlapping satellite observations is critical for global monitoring and environmental change assessments. We developed an empirical approach for inter-calibration of satellite microwave brightness temperature (T_b) records over land from the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) and Microwave Scanning Radiometer 2 (AMSR2) using overlapping T_b observations from the Microwave Radiation Imager (MWRI). Double Differencing (DD) calculations revealed significant AMSR2 and MWRI biases relative to AMSR-E. Pixel-wise linear relationships were established from overlapping T_b records and used for calibrating MWRI and AMSR2 records to the AMSR-E baseline. The integrated multi-sensor T_b record was largely consistent over the major global vegetation and climate zones; sensor biases were generally well calibrated, though residual T_b differences inherent to different sensor configurations were still present. Daily surface air temperature estimates from the calibrated AMSR2 T_b inputs also showed favorable accuracy against independent measurements from 142 global weather stations (R² ≥ 0.75, RMSE ≤ 3.64 °C), but with slightly lower accuracy than the AMSR-E baseline (R² ≥ 0.78, RMSE ≤ 3.46 °C). The proposed method is promising for generating consistent, uninterrupted global land parameter records spanning the AMSR-E and continuing AMSR2 missions.

Snow water equivalent (SWE) is a key parameter in the Earth’s energy budget and water cycle. It has been demonstrated that SWE can be retrieved using active microwave remote sensing from space. This necessitates the development of forward models that are capable of simulating the interactions of microwaves and the snow medium. Several proposed models have described snow as a collection of sphere- or ellipsoid-shaped ice particles embedded in air, while the microstructure of snow is, in reality, more complex. Natural snow usually forms a sintered structure following mechanical and thermal metamorphism processes. In this research, the bi-continuous vector radiative transfer (bi-continuous-VRT) model, which firstly constructs snow microstructure more similar to real snow and then simulates the snow backscattering signal, is used as the forward model for SWE estimation. Based on this forward model, a parameterization scheme of snow volume backscattering is proposed. A relationship between snow optical thickness and single scattering albedo at X and Ku bands is established by analyzing the database generated from the bi-continuous-VRT model. A cost function with constraints is used to solve effective albedo and optical thickness, while the absorption part of optical thickness is obtained from these two parameters. SWE is estimated after a correction for physical temperature. The estimated SWE is correlated with the measured SWE with an acceptable accuracy. Validation against two-year measurements, using the SnowScat instrument from the Nordic Snow Radar Experiment (NoSREx), shows that the estimated SWE using the presented algorithm has a root mean square error (RMSE) of 16.59 mm for the winter of 2009–2010 and 19.70 mm for the winter of 2010–2011.