The dynamic behavior of groundwater flow and salt transport is affected by tide and pumping in coastal multilayered aquifers. In this paper, two groups of experiments were conducted considering different constant head inland boundaries. The fluctuation of the groundwater level and the process of seawater intrusion in the multilayered aquifers were observed. A two-dimensional SEAWAT model is developed to simulate the seawater intrusion to coastal aquifers under the influences of tidal fluctuation and groundwater exploitation. The hydrogeological parameters in the model are calibrated by the records of the groundwater level and salinity measurements. The results showed that the simulated groundwater level and salt concentration match the observation well. The groundwater level has the characteristics of periodic fluctuation with tide. The lag time of the groundwater level fluctuation in each monitoring point increases slightly with the increasing distance from the saltwater chamber. For the low tide, the inland freshwater recharge has main effect on groundwater level fluctuation. The rising tide has a negative effect on the drawdown of the groundwater level induced by pumping. For the high tide, the tide plays a major role on groundwater level fluctuation, compared with the inland freshwater recharge. Compared with the condition of high head of inland recharge, larger saltwater intrusion lengths and area have been observed and simulated in the aquifer, which means that faster inland motion of the seawater wedge would occur when the inland recharge is small in the coastal aquifers. It revealed that inland recharge plays a major role in the seawater intrusion for the same pumping rate of groundwater in different seasons. The analysis provides insights into how the tide fluctuation, groundwater pumping, and inland recharge effect on the area and rates of seawater intrusion.

The safety of groundwater resources in coastal areas is related to the sustainable development of the national economy and society. Seawater intrusion is a serious problem threatening the groundwater environment in coastal areas. Climate change, tidal effects and groundwater exploitation may destroy the balance between salt water and fresh water in coastal aquifers, leading to seawater intrusion. The threat of seawater intrusion has attracted close attention, especially in the coastal areas of northern China, and the accuracy and efficiency of seawater intrusion monitoring need to be improved. The aim of this study was to fill the blanks in seawater intrusion research in the coastal aquifer of the Daqing River Catchment, northeastern Liaodong Bay, China, and determine the extent and evolutionary characteristics of seawater intrusion in this area. In this study, historical chloride concentration data were used to trace the evolution of the salinization, and electrical resistivity tomography (ERT) was used to supplement the data in areas with limited hydrochemical data and to detect the saltwater/freshwater interface, especially in the area near the Xihai Sluice. The results show that seawater intrusion in the Daqing River Catchment is mainly caused by overexploitation of groundwater. Since 2012, strict controls have been placed on the groundwater exploitation rate, and the chloride concentration of 250 mg/L has receded year by year, with saltwater being significantly reduced by 2018. The Daqing River plays an important role in the saltwater distribution. The Xihai Sluice, located in the lower reaches of the Daqing River, intercepts and controls the seawater intrusion in a certain range by raising the level of fresh groundwater to intercept and control saltwater intrusion within a certain range. The research results also confirmed that a combination of geophysical and geochemical methods is of great value in studying seawater intrusion, especially in areas with limited available hydrochemical data. A monitoring network with ERT instruments and wells should be established to collect regular measurements of the electrical resistivity distribution, as well as the groundwater level and chemical composition.

In the context of global warming, rising sea levels are intensifying seawater intrusion in coastal areas. Due to the complex hydrodynamic conditions and increasing groundwater over-extraction in these regions, understanding the patterns of seawater intrusion is crucial for effective prevention and control. This study employed a sandbox model to investigate both vertical and horizontal seawater intrusion into a coastal unconfined aquifer with an impermeable dam under varying conditions of sea level rise, coastal slope, and groundwater pumping rate. Additionally, a two-dimensional SEAWAT model was developed to simulate seawater intrusion under these experimental conditions. The results indicate that sea level rise significantly increases the extent and intensity of seawater intrusion. When sea level rises by 3.5 cm, 4.5 cm, and 5.5 cm, the areas of the saline wedge reached 362 cm², 852 cm², and 1240 cm², respectively, with both horizontal and vertical intrusion ranges expanding considerably. When groundwater extraction is superimposed, vertical seawater intrusion is notably intensified. At an extraction rate of 225 cm³/min, the vertical intrusion areas corresponding to sea level rises of 3.5 cm, 4.5 cm, and 5.5 cm were 495 cm², 1035 cm², and 1748 cm², respectively, showing significant expansion, and this expansion becomes more pronounced as sea levels rise. In contrast, slope variations had a significant impact only on vertical seawater intrusion. As the slope decreased from tanα = 1/5 to tanα = 1/9, the upper saline wedge area expanded from 525 cm² to 846 cm², considerably increasing the vertical intrusion range. Finally, the combined effects of groundwater extraction and sea level rise exacerbate seawater intrusion more severely than either factor alone, presenting greater challenges for coastal water resource management.

The coupled effect of seawater intrusion and inland freshwater recharge plays an important role in contamination transport in coastal heterogeneous aquifer. In this study, the effects of seawater intrusion and inland recharge on contamination transport were investigated by conducting laboratory experiments and numerical simulations. The laboratory tests were conducted in a sand tank considering two scenarios, namely the conditions of landward and seaward hydraulic gradients. The SEAWAT software was applied for validating the contaminant transport in coastal heterogeneous aquifer. The results indicated that the simulated seawater wedge and contours of the saltwater contaminant matched the observed ones well. The length of the seawater wedge in the scenario of seaward hydraulic gradient was smaller than that in the scenario of landward hydraulic gradient, which reflected that the large quantity of inland recharge have a negative effect on the invasion process of seawater. The plume moved mainly downward in the heterogeneous unconfined aquifer for both scenarios. The pollution plume became concave at the interface between each two layers, which was because the velocity of contaminant plume migration increased gradually from the upper layer to lower layer. The migration direction of the front of the plume was consistent with the direction of hydraulic gradient, which indicated that it was influenced by the water flowing. The maximum area of plume in the scenario of seaward hydraulic gradient was slightly smaller than that in the scenario of landward hydraulic gradient. The maximum area and vertical depth of the pollutant plume were sensitive to the hydraulic conductivity, dispersivity and contamination concentration. This study was of great significance to the controlling of pollution and utilization of freshwater resources in coastal areas.

In this work, the influence of the Hurst exponent and Peclet number (Pe) on the temporal mixing behavior of a conservative solute in the self-affine fractures with variable-aperture fracture and constant-aperture distributions were investigated. The mixing was quantified by the scalar dissipation rate (SDR) in fractures. The investigation shows that the variable-aperture distribution leads to local fluctuation of the temporal evolution of the SDR, whereas the temporal evolution of the SDR in the constant-aperture fractures is smoothly decreasing as a power-law function of time. The Peclet number plays a dominant role in the temporal evolution of mixing in both variable-aperture and constant-aperture fractures. In the constant-aperture fracture, the influence of Hurst exponent on the temporal evolution of the SDR becomes negligible when the Peclet number is relatively small. The longitudinal SDR can be related to the global SDR in the constant-aperture fracture when the Peclet number is relatively small. As the Peclet number increases the longitudinal SDR overpredicts the global SDR. In the variable-aperture fractures, predicting the global SDR from the longitudinal SDR is inappropriate due to the non-monotonic increase of the longitudinal concentration second moment, which results in a physically meaningless SDR.