the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A Multi-chain Surrogate-assisted Hybrid Optimization Framework for Joint Identification of Groundwater Contaminant Sources and Hydrogeological Parameters
Abstract. Rapid and accurate identification of groundwater contaminant information and hydrogeological parameters is crucial for effective groundwater remediation and risk management. Within a simulation-optimization framework, this task is inherently posed as a mixed-variable optimization problem involving discrete parameters (e.g., source locations) and continuous ones (e.g., hydraulic heads, conductivities, and release fluxes). However, several challenges arise in this context. First, conventional optimization algorithms often exhibit slow convergence and unstable performance. Second, they typically require thousands of simulations to adequately explore the complex parameter space, resulting in prohibitive computational costs. To address these issues, this study develops a surrogate-assisted hybrid algorithm that integrates the Cooperative Search Algorithm (CSA) and Tabu Search (TS) within a synergistic multi-chain optimization framework, termed SA-CSA-TS. In each iteration, individual chains first perform independent CSA-based optimization to promote broad global exploration, after which they collaboratively refine source locations through a neighbourhood search guided by a shared tabu list. In addition, surrogate models equipped with a reconstruction strategy partially replace groundwater simulations, thereby substantially reducing the computational burden. Case studies reveal that the Radial Basis Function (RBF) outperforms other mainstream surrogate models in both accuracy and stability. Furthermore, comparative experiments confirm that the proposed SA-CSA-TS framework not only achieves higher solution accuracy but also significantly reduces computational demand, demonstrating strong potential for efficient groundwater contamination diagnosis.
- Preprint
(2915 KB) - Metadata XML
-
Supplement
(202 KB) - BibTeX
- EndNote
Status: open (until 10 Feb 2026)
- CC1: 'Comment on egusphere-2025-6140', Giacomo Medici, 19 Dec 2025 reply
-
CC2: 'Comment on egusphere-2025-6140', Nima Zafarmomen, 25 Dec 2025
reply
The study introduces a significant advancement in the simulation-optimization (S−O) framework for groundwater contamination diagnosis. The novelty of the SA-CSA-TS framework lies in its synergistic multi-chain architecture. Unlike traditional single-population algorithms, the integration of the Cooperative Search Algorithm (CSA) for global exploration and Tabu Search (TS) for local refinement—guided by a shared tabu list—effectively addresses the "equifinality" and multimodality inherent in mixed-variable groundwater problems.
Furthermore, the systematic evaluation of surrogate models and the implementation of a dynamic reconstruction strategy (achieving an 85–88% reduction in computational demand) provides a highly practical blueprint for real-world remediation efforts where time and computing resources are limited.
Minor Comments:
1. In Section 3.4, the authors describe the update rules for the tabu list. It would be beneficial to briefly clarify the "tenure" or size of the tabu list. Does the list have a maximum capacity, or does it grow indefinitely throughout the FEmaxiterations?
2. While the authors mention using default settings in UQPyL, the performance of Kriging and Gaussian Processes is often highly sensitive to the choice of kernel/correlation functions. A brief sentence justifying the choice of the Cubic RBF kernel over others (like Thin Plate Spline) would add more depth to the surrogate comparison section.
3. In the discussion of the "parameter-compensation effect" (Section 7.2), the authors correctly identify that multiple locations can yield similar concentrations. It might be helpful to suggest how monitoring well placement (optimal experimental design) could potentially reduce this equifinality in future iterations of the framework.
4. In Figure 16 and Figure 18 (Radar Charts), the overlap of the GA and CSA lines can be difficult to distinguish. Consider using slightly different line textures (e.g., dashed vs. dotted) to improve accessibility for the reader.
5. To broaden the impact of the study, the authors should consider how this framework interacts with broader hydrological cycles and diverse data sources. I strongly recommend the authors consider and potentially reference studies such as: "Assimilation of sentinel‐based leaf area index for modeling surface‐ground water interactions in irrigation districts" This would help contextualize how satellite-derived data and surface-water interactions might provide additional constraints to the groundwater simulation models, potentially refining the identification of hydrogeological parameters.
Citation: https://doi.org/10.5194/egusphere-2025-6140-CC2
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 140 | 82 | 20 | 242 | 20 | 33 | 36 |
- HTML: 140
- PDF: 82
- XML: 20
- Total: 242
- Supplement: 20
- BibTeX: 33
- EndNote: 36
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General comments
Good research that needs some improvement. See my specific comments that should improve the manuscript.
Specific comments
Lines 34-35. “Groundwater contamination has become an increasingly critical issue, posing significant risks to environmental safety and public health”. Insert recent literature on groundwater contamination with an evident worldwide angle:
- Agbotui, P. Y., Firouzbehi, F., Medici, G. 2025. Review of effective porosity in sandstone aquifers: insights for representation of contaminant transport. Sustainability, 17(14), 6469.
- Sauvé, S.,Desrosiers, M. 2014. A review of what is an emerging contaminant. Chemistry Central Journal, 8(1), 15.
Line 91. You need to disclose the general aim of the research.
Line 91. You need to describe the specific objectives of your research by using numbers (e.g., i, ii, and iii).
Line 92-onwards. You need to add more information on the boundary conditions.
Line 92-onwards. Add more detail on the nature of the geological material modelled.
Line 109. Overall, 9 equations in the manuscript are too many, not all of the are necessary. Equation 2 is very well known.
Line 175. Equations on kriging (very well-known method) not necessary.
Line 515. Assign a number to this equation.
Figures and tables
Figure 1-5. Room to make the figures larger.
Figure 8. You need to discuss boundary conditions in more detail in the main body.
Figure 8. Increase the graphic resolution of the figure.