Yong Ge Zang

This study employed a coupled water‐air two‐phase flow and salt water transport model to analyze the behaviors of generated airflow in unsaturated zones and the fluctuations of salinity at the salt–fresh water interface in a two‐layered unconfined aquifer with a sloping beach surface subjected to tidal oscillations. The simulation results show that as the new dynamic steady state including effects of tidal fluctuations is reached through multiple tidal cycles, the dispersion zone in the lower salt water wedge is broadened because fresh water/salt water therein flows continuously landward or seaward during tidal cycles. The upper salt–fresh water interface exhibits more vulnerable to the tidal fluctuations, and the variation of salinity therein is periodic, which is irrelevant to the hydraulic head but is influenced by the direction and velocity of surrounding water‐flow. With the tidal level fluctuating, airflow is mainly concentrated in the lower permeable layer due to the restraint of the upper semi‐permeable layer, and the time‐lag between the pore‐air pressure and the tidal level increases with distance from the coastline. The effect of airflow in unsaturated zones can be transmitted downward, causing both the magnitude of salinity and its amplitude in the upper salt–fresh water interface to be smaller for the case with airflow than without airflow due to the resistance of airflow to water‐flow. Sensitivity analysis reveal that distributions of airflow in unsaturated zones are affected by the permeability of the upper/lower layer and the van Genuchten parameter of the lower layer, not by the van Genuchten parameter of the upper layer, whereas the salinity fluctuations in the salt–fresh water interface are affected only by soil parameters of the lower layer.

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.