The Surface Soil Moisture (SSM) products derived from microwave remote sensing have a coarse spatial resolution, therefore downscaling is required to obtain accurate SSM at high spatial resolution. An effective way to handle the stratified heterogeneity is to model for various stratifications, however the number of samples is often limited under each stratification, influencing the downscaling accuracy. In this study, a machine learning-based geostatistical model, which combines various ancillary information at fine spatial scale, is developed for spatial downscaling. The proposed support vector area-to-area regression kriging (SVATARK) model incorporates support vector regression and area-to-area kriging by considering the nonlinear relationships among variables for various stratifications. SVATARK also considers the change of support problem in the downscaling interpolation process as well as for solving the small sample size in trend prediction. The SVATARK method is evaluated in the Naqu region on the Tibetan Plateau, China to downscale the European Space Agency's (ESA) 25-km-resolution SSM product. The 1-km-resolution SSM predictions have been produced every 8 days over a six-year period (2010-2015). Compared with other five downscaling methods, the downscaled predictions from the SVATARK method performs the best with in-situ observations, resulting in a 24.4 percent reduction in root mean square error with 0.08 m3.m-3 and a 8.2 percent increase in correlation coefficient with 0.72, on average. Additionally, anomalously low SSM values, an indicator of drought, had a record low anomaly in mid-July for 2015, as noted by previous studies, indicating that SVATARK could be utilized for drought monitoring.

The distribution of the population is an essential aspect of addressing social, economic, and environmental problems. Gridded population data can provide more detailed information than census data, and multisource data from remote sensing and geographic information systems have been widely used for population estimation studies. However, due to spatial heterogeneity, the population has different distribution characteristics and variation patterns at different scales, while the relationships between multiple variables also vary with scale. This article presents a stepwise downscaling approach in that the random forest regression kriging technique is used to downscale census data to multi-resolution gridded population datasets. Using Nanjing, China, as the experimental case, population distribution maps were generated at 100 m, 500 m, and 1 km spatial resolution, and compared with the other three downscaling methods and three population products. The results demonstrated the produced gridded population maps by the proposed approach have higher accuracy and more accurate details of population distribution with the smallest mean absolute error (MAE) and root mean squared error (RMSE) values of 1.590 and 2.189 ten thousand people (over 40% reduction). The artificial land and road data are the two most important indicators of population distribution for the regional random forest modeling in Nanjing. Our proposed method can be a valuable tool for population mapping and has the potential to monitor sustainable development goals.

Accurate spatial distribution of gridded gross domestic product (GDP) data is crucial for revealing regional disparities within administrative units, thus facilitating a deeper understanding of regional economic dynamics, industrial distribution, and urbanization trends. The existing GDP spatial models often rely on prediction residuals for model evaluation or utilize residual distribution to improve the final accuracy, frequently overlooking the modifiable areal unit problem within residual distribution. This paper introduces a hybrid downscaling model that combines random forest and area-to-area kriging to map gridded GDP. Employing Thailand as a case study, GDP distribution maps were generated at a 1 km spatial resolution for the year 2015 and compared with five alternative downscaling methods and an existing GDP product. The results demonstrate that the proposed approach yields higher accuracy and greater precision in detailing GDP distribution, as evidenced by the smallest mean absolute error and root mean squared error values, which stand at USD 256.458 and 699.348 ten million, respectively. Among the four different sets of auxiliary variables considered, one consistently exhibited a higher prediction accuracy. This particular set of auxiliary variables integrated classification-based variables, illustrating the advantages of incorporating such integrated variables into modeling while accounting for classification characteristics.