Quantifying the spatiotemporal variations of basin-scale surface water (SW)–groundwater (GW) interactions is vital for the conjunctive management of water resources in the basin. In this study, an integrated hydrological model (SWAT-MODFLOW) is used to simulate the SW–GW system in the Huashan Basin. The numerical model was calibrated and validated using the streamflow observations of the watershed outlet and the groundwater levels of the long-term monitoring wells from 2016 to 2020 in the study area. The model results show that the SWAT–MODFLOW can achieve a better fit for the streamflow discharge, compared with the results in the single SWAT model, with R² (coefficient of correlation) and NSE (Nash-Sutcliffe efficiency coefficient) of 0.85 and 0.83, respectively. The water table fitting results indicate that R² and RMSE can reach 0.95 and 0.88, respectively. The water budgets analysis demonstrates that the average rate (0.5281 m³/s) of GW abstraction to SW is larger than the rate (0.1289 m³/s) of SW recharge to GW. Moreover, the exchange rate of SW and GW gradually reaches a peak value from June to August, and the lowest value is shown in April, for each hydrological year. Based on the IPPC6 CanESM5 dataset supplied by the Canadian Climate Centre, the regional precipitation scenario subject to climate change was predicted by the ASD (Auto Statistical Downscaling Model) a statistical downscaling method, under the climate scenarios of SSP2_4.5 and SSP5_8.5. The SW–GW interaction pattern was modeled under the future scenarios in the study area. The current (2016–2020) average annual rate of the SW–GW interaction is considered as the base value. Subject to the SSP2_4.5 scenario, the average exchange rate of the SW recharge to GW is 0.1583 m³/s, which is an increase of 22.8%. The average exchange rate of the GW discharge to SW is 0.5189 m³/s which is a reduction of 0.017%. Subject to the SSP5_8.5 scenario, the average exchange rate of SW recharge to GW is 0.1469 m³/s, which is an increase of 14.7%. The average exchange rate of the GW discharge to SW is 0.5953 m³/s, which is an increases of 12.7%. The results can assist in water resource management in the basin, by identifying potential locations of nutrient transport from the aquifer to the river, as well as changes in spatial variability under future climatic conditions.

The combined hydrochemical analysis, factor analysis, and isotopic signals of water and nitrate were applied to explore the hydrochemical origin and identify the sources and transformation of nitrate in river water and groundwater in the Huashan watershed. Additionally, a Bayesian isotope mixing model (SIAR) was employed for quantitative assessment of the nitrate sources. The results indicated that both river water and groundwater were dominated by HCO₃-Ca and HCO₃-Ca·Mg types; both originated from precipitation and were influenced by evaporation. The main constituent ions in the river water and groundwater primarily originated from carbonate and silicate dissolution, with the presence of cation exchange in the groundwater. The water chemistry of river water was greatly influenced by physicochemical factors, while that of groundwater was mainly controlled by water–rock interaction. NO₃⁻ in river water was mainly influenced by soil nitrogen (SN) and manure and septic wastes (MSWs), while NO₃⁻ in groundwater was jointly affected by ammonium fertilizers (AF), SN, and MSWs. With the exception of denitrification observed in the groundwater at the watershed outlet, denitrification was absent in both groundwater in the piedmont area and in river water. The SIAR model results demonstrated that the contribution rates of atmospheric precipitation (AP), AF, SN, and MSWs to river water were 12%, 21%, 25%, and 42%, respectively, while to groundwater, they were 16%, 27%, 10%, and 47%, respectively. Overall, MSWs were the main sources of nitrate in the river water and groundwater. It is necessary to prevent the leakage of MSWs when managing water resources.

The measurement of river discharge is essential for sustainable water resource management. The velocity–area approach is the most common method for calculating river discharge. Although several velocity measurement methods exist, they often have varying degrees of technical issues attributed to their operational complexity, time effectiveness, accuracy, and environmental impact. To address these issues, we propose a three-dimensional (3D) hydrodynamic model coupled with data assimilation (DA) for velocity measurement with improved accuracy and efficiency. We then apply this model to the Lanxi River reach in Zhejiang Province, China. The experimental results confirm that the obtained assimilated velocities using our proposed algorithm are much closer to the observed velocities than the simulated velocities. Our results show that when using the proposed method, the RMSE is decreased by 78%, and the SKILL and DASS values are 0.96 and 0.92, respectively. These confirm that the DA scheme of the flow velocity measurement is effective and capable of significantly improving the accuracy of the velocity with lower computational complexity.

To better understand processes involved in the evaporative drying of the soil, simulations on the dynamics of the evaporation zone, condensation zone, and dry surface layer (DSL) were conducted during a 10-day drying event under diurnal atmospheric conditions. Simulated water contents and soil temperatures matched well with the measured data in the lysimeter. Surface evaporation predominantly occurred during the early period each day, while subsurface evaporation dominated during the remaining part of the day. The evaporation zone presented a distinctly diurnal pattern, moving toward the deeper soil layer during the daytime and back toward the soil surface during the nighttime. The DSL and condensation zone, located immediately above and below the evaporation zone, respectively, also presented diurnal patterns following those of the evaporation zone. As soil drying progressed, both the position of the evaporation zone within the profile and the DSL width exhibited an overall increasing trend, reaching about 4.9 mm by the end of the study period. The occurrence of condensation zones was limited to the daytime when there was a downward surface temperature gradient present. Diurnal patterns observed in both evaporation zones and DSL could potentially be determined by quantifying changes in the near-surface profile’s soil water content, relative humidity, pressure head, and vapor density.