Pollution loads pose a major threat to the health of the marine environment and the long-term viability of the coastal economy. The present study developed a coupling model to simulate the chemical oxygen demand (COD) transport in upper rivers (1D) and subsequent diffusion in the coastal zone (2D) in Liaodong Bay, based on the HydroInfo system. Three main seagoing rivers, including Daliao, Liao, and Daling Rivers, were selected and investigated for hydrodynamic and hydrochemical analyses. The mathematical model was evaluated by monitoring data from state-controlled cross-sections scattered along the three rivers, and the observation data showed good agreement with simulated values, confirming the model’s accuracy in terms of spatial and temporal distribution. The transport and propagation process of COD in inlet rivers, such as Daliao, Liao, and Daling, including the sea area of Liaodong Bay, were simulated and analyzed. Simulated results revealed that the pollution range of COD in Liaodong Bay was 258–391 km² in different seasons. The pollutant leakage scenarios for the three rivers entering the sea were simulated utilizing the developed mathematical model. The study simulated and predicted that, in the event of a sudden water pollution accident (e.g., sneak discharge and leakage at various sections of sea-entering rivers, such as Daliao, Liao, and Daling), pollutants might take 2–11 days to reach the sea-entering mouth, and the sea area would take 8–32 days to reach the maximum pollution range. Our numerical modeling may be used to analyze and make decisions on pollution control in Liaodong Bay and major sea-entry rivers, and be useful to water environment management in sea-entry rivers and Liaodong Bay, and water pollution emergency responses.

The recharge efficiency during artificial groundwater recharge (AGR) is reduced primarily by clogging that is triggered by suspended particles. However, there are loopholes in the current standards of recharge-water quality for clogging control during AGR, and the threshold values of turbidity to prevent clogging have not been reasonably determined. In this study, secondary effluents from wastewater treatment plants (WWTPs) were injected into saturated sand columns to simulate the process of AGR. Batch experiments under different turbidity conditions were conducted, and the numerical modeling of particle transport and deposition was performed to assess the clogging processes. Theories of single-collector contact and interfacial interaction energy were applied to elucidate possible microcosmic mechanisms. The results showed that the diluted secondary effluent (SE) with turbidities of 0.540 ± 0.050, 1.09 ± 0.050, and 1.84 ± 0.060 NTU caused considerable clogging in the porous media, which decreased the relative hydraulic conductivities (K/K₀) by 13.2%, 17.6%, and 83.6%, respectively. The filtered SE with a turbidity of 0.160 NTU did not cause clogging, and K/K₀ was reduced by only 1.70%. The clogging was attributed to the deposition of suspended particles in the sand matrix because they have a high collision efficiency (0.007–1.98) and attachment efficiency (0.029–0.589 k_BT). Finally, this paper recommends that the turbidity of the recharge water should not exceed 0.500 NTU during AGR practices.