Downscaling WGHM-Based Groundwater Storage Using Random Forest Method: A Regional Study over Qazvin Plain, Iran

Abstract

Climate change, urbanization, and a growing population have led to a rapid increase in groundwater (GW) use. As a result, monitoring groundwater changes is essential for water managers and decision-makers. Due to the lack of reliable and insufficient in situ information, remote sensing and hydrological models may be counted as alternative sources to assess GW storage changes on regional and global scales. However, often, these hydrological models have a low spatial resolution for water-related applications on a small scale. Therefore, the main purpose of this study is to downscale the GW storage anomaly (GWSA) of the WaterGAP Global Hydrology Model (WGHM) from a coarse (0.5 degrees) to a finer spatial resolution (0.1 degrees) using fine spatial resolution auxiliary datasets (0.1 degrees), such as evaporation (E), surface (SRO), subsurface runoff (SSRO), snow depth (SD), and volumetric soil water (SWVL), from the ERA5-Land model, as well as the global precipitation (Pre) measurement (GPM-IMERG) product. The Qazvin Plain in central Iran was selected as the case study region, as it faces a severe decline in GW resources. Different statistical regression models were tested for the GWSA downscaling to find the most suitable method. Moreover, since different water budget components (such as precipitation or storage) are known to have temporal lead or lag relative to each other, the approach also incorporates a time shift factor. The most suitable regression model with the highest skill score during the training-validation was selected and applied to predict the final 0.1-degree GWSA. The downscaled results showed high agreement with the in situ groundwater levels over the Qazvin Plain on both interannual and monthly time scales, with a correlation coefficient of 0.989 and 0.62, respectively. Moreover, the downscaled product represents clear proof that the developed downscaling technique is able to learn from high-resolution auxiliary data to capture GWSA features at a higher spatial resolution. The major benefit of the proposed method lies in the utilization of only the auxiliary data that are available with global coverage and are free of charge, while not requiring in situ GW records for training or prediction. Therefore, the proposed downscaling technique can potentially be applied at a global scale and to aquifers in other geographical regions.