Assessment of contamination dispersion in a porous aquifer environment using explicit and implicit methods in the MODFLOW–MODPATH model

Abstract

This study presents a meticulously edited article that assesses the uncertainty associated with contamination dispersion in a porous aquifer environment. The article employs a conceptual model incorporating finite difference flow modeling and particle tracking methods grounded in particle movement theory within groundwater flow. The groundwater flow model utilizes an automatic calibration technique, PES, and independent verification, to minimize statistical discrepancies in the optimized parameters. The statistical summary of the model output reveals the average pollution concentration in the study area, specifically the landfill aquifer in northern Iran. The results indicate that the initial contamination value at the boundaries is zero, whereas 91,802 units are discharged into the sea. Moreover, approximately 3,317,290,003 pollution units have been extracted from pumping wells along the contamination plume over 10 years. Furthermore, the drainage network releases 1,778,680,001 pollution units, resulting in a net outflow of 5,082,740,007 units from the same drainage network. Additionally, 9,912,180,003 units were discharged from the leaky boundaries. It should be noted that these statistics can be reversed depending on the specified period. The difference between the input and output pollution values in the study area indicates a net accumulation of 8192 units. To analyze the sensitivity of contamination dispersion in the porous aquifer environment and the variations in coastal drainage patterns obtained from the finite difference simulation, the explicit and implicit methods, as well as the random walk particle tracking method, were investigated in a portion of the basin area using the South Side Implicit method. One of the main objectives of this sensitivity analysis is to evaluate the results of the uncalibrated qualitative model and examine the particle dispersion in the adjacent area of the contaminated landfill region to assess the error of the proposed model by developing a limited test. The calculated error or computational difference between the explicit and implicit methods and the applet model is considered negligible. For example, the calculated error between the explicit and backward methods and the applet model ranges from 0.086 to 0.372 units for a concentration interval of 0.512 units and 0.680 units for a time interval 0.537. Therefore, the computational differences are not significant. Furthermore, based on theoretical sensitivity analysis of contamination dispersion, the laboratory model demonstrates that the flow direction in this modeling is from west to east, as expected.