Numerical Groundwater Flow Model Research Articles

Impacts of groundwater dynamics around a macro-tidal river on agricultural soil salinity

Estuaries are vulnerable to oceanic and atmospheric climate change. Much of the research investigating climate change impacts on estuaries is focused on saltwater intrusion within surface water due to drought and rising sea levels, with implications for ecosystems and humans. Groundwater and soil near estuaries may also be influenced, as estuary salinity and hydraulic head changes can impact soils and aquifers not previously at risk of salinization. This study was conducted to address knowledge gaps related to present and future groundwater salinity distribution in a groundwater system connected to a macro-tidal estuary. The studied estuary experiences a tidal bore due to its hydraulic connection to the Bay of Fundy in Nova Scotia, Canada. A parcel of agricultural land adjacent to the estuary was selected to assess the groundwater response to episodic fluctuations in estuary water levels and salinity. Groundwater monitoring and electromagnetic surveys were conducted to map soil and groundwater salinity patterns. A numerical model of groundwater flow and solute transport informed by field data was used to investigate how varying estuary salinity due to droughts and sea-level rise could impact groundwater salinity. Results showed that, in contrast to salt wedges observed along marine coasts, the saline groundwater existed as a plume immediately around the estuary. Model simulations showed that short-term droughts had an insignificant impact on the adjacent groundwater salinity. However, permanent increases in salinity caused by sea-level rise increased the plume volume by 86 %, or an additional ∼11 m horizontally and ∼ 4.5 m vertically. Our results suggest that increased river salinity in this setting would not result in widespread salinization of porewater and agricultural soils, but more extensive salinization may be experienced in permeable aquifers or along more saline estuarine zones. Findings may inform land management decisions in regions exposed to increased salinity in the future.

Optimal pumping policy from well field connected with a stream

Conjunctive use management of water resources is of paramount importance in the backdrop of the growing demand for freshwater resources. With this viewpoint, this study is aimed to present a methodological approach utilizing the combined use of simulation and optimization models for evolving strategies or policies for the management of water resources in a well-field connected to a stream. The proposed methodology incorporates the coupling of the groundwater simulation model with the optimization model for decision-making through the response matrix method. In this study, a numerical groundwater flow model using MODFLOW was used to simulate the response of stream-aquifer interactions. A linear optimization model was formulated to account for the best optimal groundwater withdrawal taking into consideration the allowable drawdown, permissible stream flow depletion, and satisfaction of the water demands. Such management options permit the investigation of alternative strategies for groundwater abstraction in a well-field hydraulically connected to a stream for sustainable development of water resources. The model was applied in a hypothetical well field study area hydraulically connected to a stream and satisfactory results were obtained. The methodology was also very important to understand the tradeoff between groundwater withdrawal and the depletion of flow to the nearby stream. This approach was applied to the Aynalem well field stream-aquifer system which is located in Northern Ethiopia. Results show that groundwater extraction from the Aynalem well field triggers stream flow depletion and it is possible to decrease the stream flow depletion by 3.5–4.5% by adjusting current groundwater withdrawal schedules and configuration of existing wells.

Slope stabilization through groundwater management with limited hydrogeological data: a case study from Majes, southern Peru

Water table rise near a cliff may trigger a landslide due to the associated increase pore pressure and decrease in frictional resistance. One main cause of water table rise is intense irrigation for agriculture in arid and semi-arid regions. One such case is in Majes, southern Peru, where a landslide has evolved near an intensively irrigated agricultural area. Mitigation strategies for landslides exist, such as physical strengthening of the cliff, but can be expensive. We describe a groundwater management approach to reduce the pore pressure in the vicinity of the cliff to either slow the propagation of an existing landslide or prevent the initiation of a new landslide. A 3D numerical groundwater flow model was built for the Majes area which employs the limited data existing on the local hydrogeology. Simulations were run to understand the connection between the hydraulic properties and the water table level change due to irrigation and pumping. Results show that through a series of pumping wells near the cliff edge, the pore pressure can be decreased effectively. Moreover, decreasing the water table via pumping can be accomplished in 25–35% of the time it took to elevate the water table level by irrigation. In addition, the pumping can capture water that could be reused for irrigation. Thus, based on our analysis, we conclude that wells could provide a groundwater management approach that keeps the pore pressure at low levels to mitigate landslide processes and simultaneously supplies water for irrigation existing and future irrigation-heavy agriculture in semi-arid environments.

Numerical groundwater flow modeling under future climate change in the Central Rift Valley Lakes Basin; Ethiopia

Study area: Katar and Meki subbasins, Rift Valley Lakes Basin, Ethiopia.Study focus: This research was carried out to characterized the recharge mechanism and quantify the steady-state groundwater balance and its sensitivity to future climate change. A groundwater simulation model was constructed and calibrated using a hydro-geo spatial dataset. Three regional climate models were used to assess the potential impact of changes in future precipitation on the recharge rate and groundwater balance components.New Hydrogeological Insight: Groundwater potential assessment depends on accurate estimation of the recharge rate. Precipitation contributed 11.95% and 11.96% to groundwater recharge in the Katar and Meki subbasins, respectively. The steady-state numerical groundwater model was calibrated and the model performed in the ranges of R2: 0.95–0.99; RMSE: 16.17–25.18; and MAE: 12.69–24.55, demonstrating 'excellent' model performance. In particular, the model exhibited high sensitivity to changes in the recharge rate and horizontal hydraulic conductivity. Future change in precipitation caused a reduction in groundwater potential in the range of 6.24–40.32% by the 2040 s and 2070 s, respectively, in the Katar subbasin. Likewise, the Meki Subbasin will experience a reduction in groundwater potential in the range of 0.29–37.17% by the 2040 s and 2070 s, respectively. These results emphasize how crucial it is for future water resource development initiatives to take into account climate variability for sustainable groundwater development.

Evaluating groundwater resources trends through multiple conceptual models and GRACE satellite data.

In this research, three numerical groundwater flow models, developed and calibrated from three equally plausible conceptual models over the Nasia Basin, have been used to assess groundwater resources variations over a transient period. The use of multiple numerical models reduces the effect of uncertainties in conceptual model formulation. All the three calibrated numerical models indicate an increasing trend of groundwater recharge and storage over the period of the groundwater level monitoring. This suggests that the prevailing erratic climatic conditions in the area are conducive for increasing groundwater recharge and storage in the terrain. The high-intensity, short duration rainfall patterns, attending climate change in the basin, enhance high levels of infiltration and percolation, leading to steadily increasing groundwater recharge. Groundwater recharge estimates from each of the models over the transient period appear to reflect the pattern of seasonal variations in rainfall in the region. Data from the models indicates a significant role of baseflow in sustaining perennial streamflow in the area. This presents a significant development in terms of groundwater-based adaptation projects, especially in agriculture. The trend of groundwater recharge in the Nasia Basin is in sync with regional groundwater storage variations estimated from the Gravity Recovery and Climate Experiment (GRACE) satellite data collected and processed over the Volta Basin. At the Volta Basin level, groundwater storage variations indicate a strong positive trend of increasing groundwater recharge from 2002 (beginning of the GRACE mission) to 2022 (end point of the data used for this research). Analysis of the GRACE data suggests that there is a cumulative increase in groundwater storage by 30cm, representing approximately 120 km3 of groundwater over the period in the basin. This translates into approximately 15mm/year of groundwater storage increase. Thus, at both the regional and local levels, groundwater appears to be responding positively to the impacts of erratic rainfall patterns observed in the area recently. The high-intensity, short duration rainfall patterns appear to favor significant groundwater recharge, resulting in a strong positive groundwater storage signal. The high positive groundwater storage signal suggests increasing groundwater resources potential in the area, indicating promising opportunities for groundwater-based climate change adaptation interventions.

ERT data assimilation to characterize aquifer hydraulic conductivity heterogeneity through a heat‐tracing experiment

AbstractGeothermal energy systems, such as heat pumps relying on aquifers, use renewable sources of energy that are accessible in urban areas. It is necessary to characterize the subsurface hydraulic properties prior to the installation of such systems. In this context, a heat‐tracing experiment is a typical field test that can help with the characterization of the subsurface. During a heat‐tracing experiment, monitoring with downhole temperature sensors, water‐level pressure transducers and electrical resistivity tomography (ERT) can be used to help characterize the hydrogeological properties. Previous monitoring tools have shortcomings, such as low‐resolution data and over‐smoothing; thus, they fail to reproduce the heterogeneity of hydrogeological properties. Ensemble Kalman filter (EnKF) is a promising tool that can overcome the over‐smoothing problem to replicate the hydrogeological property heterogeneity. In this work, we proposed a new procedure to assimilate time‐lapse cross‐borehole ERT data into a numerical model of groundwater flow and heat transfer, where the groundwater is extracted and heated water is reinjected into an unconfined sandy‐gravel aquifer. The finite element model (FEFLOW 7.3) of groundwater flow and heat transfer is integrated with petrophysical relationship and electrical forward modelling (ResIPy) to estimate cross‐borehole ERT measurements. Then, the estimated apparent resistivity is assimilated to update the hydraulic conductivity model using EnKF. The results of the application of the proposed approach to an experimental site located in Quebec City (Canada) demonstrate that the heterogeneity of K is correctly reproduced as the updated K model is reasonably consistent with the lithological log. In addition, the proposed approach was able to replicate the cross‐borehole ERT field and temperature measurements. The comparison between prior and posterior distribution of K with slug test results shows that the EnKF made the final (assimilated) distribution of K move towards K values inferred with slug tests.

Numerical simulation of the interaction between mine water drainage and recharge: A case study of Wutongzhuang coal mine in Heibei Province, China

Water inflow poses a critical safety concern in coal mining. The pressing challenge in current mine water inflow research is the need to enhance the prediction of dynamic changes in mine water inflow across different mining stages. In this study, the hydrogeological conditions and aquifer structural characteristics of the Wutongzhuang Coal Mine in Hebei Province, China were analyzed, and the results were combined with pumping test and drilling data. A numerical model of groundwater flow in the study area was established using Visual MODFLOW and the model was calibrated and verified by comparing its results with the measured values to ensure that it realistically reflected the state of groundwater movement in the study area. The results of numerical simulation indicated that: (1) The water level of the Ordovician aquifer varies significantly by season. A large exposed area of Ordovician limestone in the western mountainous region receives a large amount of rainfall infiltration recharge in the summer rainy period that raises the water level of the Ordovician aquifer in the study area. (2) The groundwater system in the study area is in an equilibrium state, with the water inflow approximately equal to the outflow. (3) The recharge of mine drainage water has a slight positive effect on the Ordovician water level and the floor water inflow, with the impact on the water level of the aquifer decreasing as the distance from the recharge well and the import of Ordovician limestone water increase. In cases in which there is mine water recharge, the water level of the Ordovician aquifer is approximately 0.33 to 1.37 m higher than in cases without recharge, resulting in an increase in water inflow of approximately 10.6 m3/d. (4) In cases in which there is recharge of all drainage water, the maximum water inflow into the fifth mining area is 3,491.3 m3/d, which is 12.5 m3/d higher than in the case without recharge. These results provide important technical references for the formulation of mining schemes and water prevention measures in mining areas.

Numerical modelling of groundwater flow in Nambiyar River Basin

Numerical modelling of groundwater flow in Nambiyar River Basin

Three-dimensional inter-basin groundwater flow toward a groundwater-fed stream: Identification, partition, and quantification

Inter-basin groundwater flow (IGF), defined as groundwater flows across the topographic divides, occurs in nature as a common phenomenon, and seriously affect water balance and solute transport in a basin. Although the existence of IGF has been checked in drainage basins with indirect methods, such as using information of hydrogeochemistry, few studies are conducted to directly identify and quantify the IGF in the three-dimensional (3D) space using numerical modeling of groundwater flow. In this study, we identify, partition and quantify IGF based on a regional-scale 3D groundwater flow model covering several traditional surface water basins. The investigation is focused on the Hailiutu River, a groundwater-fed stream in the Ordos Plateau, China. It is the first time that the source zones and 3D paths of IGF as well as its contribution to the streamflow are fully captured for a groundwater-fed stream. Results indicate that the IGF contributes 36%∼59% to the baseflow and 32%∼52% to the streamflow of the Hailiutu River, through different groups of circulation paths with the maximum depth up to ∼ 1000 m. Uncertainty analyses of the model show that the contribution of IGF to the streamflow is positively correlated to the infiltration recharge and hydraulic conductivities of the aquifer media. Several indirect methods are also checked for comparison, but they could only be used to assess the potential contribution of IGF, not the IGF paths. This study will facilitate in deeply understanding the evolution mechanism of the inter-basin groundwater circulation and the interaction between surface water and groundwater.

Quantifying aquifer interaction using numerical groundwater flow model evaluated by environmental water tracer data: Application to the data-scarce area of the Bandung groundwater basin, West Java, Indonesia

Study Region:Bandung groundwater basin, Indonesia. Study focus:Groundwater abstraction of various magnitudes, pumped out from numerous depths in a multitude of layers of aquifers, stimulates different changes in hydraulic head distribution, including ones under vertical cross-sections. This generates groundwater flow in the vertical direction, where groundwater flows within its storage from the shallow to the underlying confined aquifers. In the Bandung groundwater basin, previous studies have identified such processes, but quantitative evaluations have never been conducted, with data scarcity mainly standing as one of the major challenges. In this study, we utilize the collated (1) environmental water tracer data, including major ion elements (Na+/K+, Ca2＋, Mg2+, Cl−, SO42−,HCO3−), stable isotope data (2H and δ18O), and groundwater age determination (14C), in conjunction with (2) groundwater flow modeling to quantify the aquifer interaction, driven mainly by the multi-layer groundwater abstraction in the Bandung groundwater basin, and demonstrate their correspondence. In addition, we also use the model to quantify the impact of multi-layer groundwater abstraction on the spatial distribution of the groundwater level changes. New hydrological insights for the region:In response to the limited calibration data availability, we expand the typical model calibration that makes use of the groundwater level observations, with in-situ measurement and a novel qualitative approach using the collated environmental water tracers (EWT) data for the model evaluation. The analysis in the study area using EWT data and quantitative methods of numerical groundwater flow modeling is found to collaborate with each other. Both methods show agreement in their assessment of (1) the groundwater recharge spatial distribution, (2) the regional groundwater flow direction, (3) the groundwater age estimates, and (4) the identification of aquifer interaction. On average, the downwelling to the deeper aquifer is quantified at 0.110 m/year, which stands out as a significant component compared to other groundwater fluxes in the system. We also determine the unconfined aquifer storage volume decrease, calculated from the change in the groundwater table, resulting in an average declining rate of 51 Mm3/year. This number shows that the upper aquifer storage is dwindling at a rate disproportionate to its groundwater abstraction, hugely influenced by losses to the deeper aquifer. The outflow to the deeper aquifer contributes to 60.3% of the total groundwater storage lost, despite representing only 32.3% of the total groundwater abstraction. This study shows the possibility of quantification of aquifer interaction and groundwater level change dynamics driven by multi-layer groundwater abstraction in a multi-layer hydrogeological setting, even in a data-scarce environment. Applying such methods can assist in deriving basin-scale groundwater policies and management strategies under the changing anthropogenic and climatic factors, thereby ensuring sustainable groundwater management.

Numerical simulation of different pollutant control measures around an old landfill contaminated site: A field scale study

The remediation of contaminated soils is a great challenge for global environmental sciences and engineering. The landfill was a kind of infrastructure to deal with waste from different sources while it would also cause the threat to groundwater. Cut-off walls and pumping wells were usually applied in the landfill to prevent the spread of pollutants to wider areas. However, the combination of using both of methods was rarely analyzed, especially using field data for calibrating and fitting groundwater flow and pollutant transport. 7 monitoring wells were arranged in the study area to survey the subsurface seepage. The pollution monitoring was carried out for a period of 50 days, covering 31 types of inorganic and organic pollutants. The concentration of 2,4,6-trichlorophenol (TCP) was 556.7 times greater than the standard concentration. A coupled numerical model of groundwater flow and pollutant transport was developed to assess the effectiveness of various control methods. Three options were tested, including the implementation of a single cut-off wall as well as a combination of a cut-off wall and a pumping well, for preventing the discharge of pollutants from landfills. The combination of a cut-off wall and a pumping well is the best strategy for removal of TCP. The combination approaches lead to a reduction of pollution plumes by a factor of 11 compared to the case without pollution control measures. The research findings may provide a basis and reference for the application of cutoff walls and pumping well in landfill sites or contaminated groundwater.

Rock core VOC profiles diagnostic of aquitard occurrence and integrity in a multi-layered sedimentary rock aquifer flow system

Aquitard assessment poses a challenge for practitioners and researchers managing dense non-aqueous phase liquid (DNAPL) contaminants in sedimentary rock aquifer systems. Within these systems, DNAPL migration and bulk hydraulic properties are controlled by fracture networks which are difficult to characterize at relevant scales with conventional methods. A previous study at a DNAPL contaminated site used contrasts in high-resolution vertical gradients to delineate four aquitards and two surfaces with poor vertical connectivity within a thick (>200 m) sedimentary rock aquifer system primarily composed of layered sandstones. In this study, we corroborate that method by using the existing DNAPL contamination as a tracer to identify and characterize the aquitards in the sequence and assess aquitard integrity. Closely spaced samples of rock from four continuous cores were collected and analyzed for volatile organic compounds (VOCs). The cores were collected along a transect perpendicular to groundwater flow and 75 m downgradient of the DNAPL source zone. Dissolved phase contaminant concentrations were collected using high-resolution multilevel systems placed along the plume centerline, downgradient from the source. The rock core and groundwater VOC concentration profiles suggest a primary aquitard limits vertical migration of both DNAPL and the associated dissolved phase contaminants, despite relatively large, downward gradients across its lower boundary. A calibrated, three-dimensional numerical groundwater flow model indicates an anisotropy factor of 1,000 for this aquitard. The Kv values for all four aquitards are only slightly lower than the values for the aquifers. However, the aquitard Kv values are 1–3 orders of magnitude less than their Kh values, and the aquitard Kh values are often similar to or even greater than the Kh values for the aquifers. Therefore, large anisotropy, rather than exceptionally low Kv values, distinguish the aquitards from the aquifers in this system. Fracture network characterization of nearby outcrops showed the anisotropy of the aquitards can be attributed to laterally extensive, bedding parallel partings coupled with short, frequently terminating, and poorly connected vertical fractures. In addition, the rock VOC profiles combined with fracture network data show that fracture termination surfaces can be powerful inhibitors of downward DNAPL migration providing additional evidence for the importance of these features in sedimentary rock flow systems. Ultimately, this study reveals the important influence of anisotropic aquitards on DNAPL migration, accumulation, and persistence, as well as associated plume migration. Furthermore, it demonstrates that subtle contrasts critical to aquitard integrity may be easily missed with long screened wells and other conventional, low resolution site characterization methods.

Numerical modeling of groundwater flow and nitrate transport in karst and wastewater irrigated agricultural and forest landscapes, PA, USA

AbstractUnderstanding the influence of treated wastewater reuse for irrigation on nitrate groundwater contamination in karstic landscape is critical for developing long‐term management plans and assessing the sustainability of these practices. The Pennsylvania State University has been operating a secondary treated wastewater effluent irrigation site—the “Living Filter”—since the 1960s. The objective of this study is to develop a regional groundwater flow and nitrate transport model to assess the impact of treated wastewater and disposal of large volumes of wastewater for irrigation of agricultural and forested lands on nitrate groundwater contamination in a karst landscape in the Spring Creek Watershed in central Pennsylvania. The modular three‐dimensional finite‐difference ground‐water flow model (MODFLOW) was used to construct the groundwater flow model of the local aquifer and simulate three‐dimensional regional changes in the water table over time. The simulation of nitrate fate and transport in the aquifer was conducted by developing a groundwater mass transport model using MODFLOW and a modular three‐dimensional transport multi‐species model (MT3DMS). The calibrated groundwater model yielded a root mean square error of 5.41 ft for the area of irrigation and 24.67 ft for the entire modeled area. Results of the calibrated model show the extension of a nitrate plume over most of the study area with concentrations below 10 mg·L−1 under the current 2015 setting and the long‐term 2035 scenario of increased irrigation rates. The nitrate migration was dominated by the presence of two high conductivity groundwater troughs adjacent to the study area resulting in a near steady‐state total mass since the 1990s.

A Novel Simulation-Optimization Model Built by FloPy: Pollutant Traceability in a Chemical Park in China

Heavy metal pollution of groundwater will not only destroy the ecological environment but also negatively affect the functioning of the human liver. Tracing the source of groundwater pollution is an important way to protect groundwater resources. FloPy is promoting the use of big data in the groundwater field, especially in groundwater resource planning and management and contaminant traceability. This paper takes Mn as an example and codes a simulation-optimization model for solving the groundwater pollutant traceability problem using FloPy. The Bayesian optimization and strengthen elitist genetic algorithm (SEGA) algorithms are then used to optimize the hydraulic conductivity and pollutant sources in the study area. The results show that the model runs in 411 s, which is an acceptable amount of time spent, the slope of the fitted curve between the model-calculated water level and the actual observed water level is 0.914, and the contaminant traceability results can successfully locate the contaminant sources in real engineering problems. The numerical groundwater flow model and solute transport model can be quickly built, modified, and run by writing code, and can be easily and efficiently coupled with various optimization algorithms with FloPy.

Numerical modeling of site-scale groundwater flow with stochastic parameterized hydraulic conductivity fields for geological disposal of high-level radioactive waste in China

Deep geological disposal is internationally accepted as a feasible and safe approach for high-level radioactive waste (HLW) disposal. To evaluate hydrogeological conditions in support of safety assessments for HLW disposal at the Xinchang preselected site in China, stochastic inverse models of site-scale groundwater flow were developed in this study to generate highly parameterized hydraulic conductivity fields for the fractured medium by adopting pilot point and null-space Monte Carlo methods. The hydraulic influence zones of the faults were carefully considered to have homogenous parameter values across the zones. A total of 120 multiple random realizations were generated to characterize the 3D spatial variations in hydraulic conductivity by conditioning to 1218 fixed pilot points. The simulated hydraulic heads of the realizations were in good agreement with observational data. The hydraulic conductivity fields close to the boreholes could be characterized in detail and had relatively low uncertainties. Furthermore, the hydraulic conductivity (K) values around underground research laboratory (URL) at the HLW disposal zone depth were less than 10−9 m/s. The K values obtained for most of the influence zones of these faults at this depth were consistent with hydraulic test data and ranged from 10−10 to 10−9 m/s, thereby indicating that the surrounding rock has low permeability and is suitable for HLW disposal. These findings offer new insights into hydrogeological conditions by providing key information on the permeability characteristics that need further focus on HLW disposal.

Permeability Tests and Numerical Simulation of Argillaceous Dolomite in the Jurong Pumped-Storage Power Station, China

Due to its poor hydro-physical properties and other characteristics, argillaceous dolomite is susceptible to seepage failure under high water pressure, affecting the seepage stability of a rock mass. To ensure the safety of the project, when the argillaceous dolomite is present, it is necessary to study the conditions pertaining to its seepage failure. Taking the argillaceous dolomite of Jurong Pumped Storage Power Station as the research object, the spatial distribution, occurrence, scale, degree of weathering, and mechanical and hydrogeological characteristics of the argillaceous dolomite were studied. Through on-site water pressure tests and laboratory variable head tests, the permeability characteristics of argillaceous dolomite were analyzed, and the hydraulic conductivity of the argillaceous dolomite in the upper reservoir and underground powerhouse areas was quantified. The argillaceous dolomite specimens were collected, and seepage failure tests were conducted to determine the critical water pressure for its seepage failure. Based on the results of the laboratory tests, a numerical model of groundwater flow was established. By changing the water level of the upper reservoir and the measures of the anti-seepage and drainage, the seepage stability of the argillaceous dolomite was discussed. The actual water pressure of argillaceous dolomite in the underground powerhouse area was identified during the operation of the Jurong pumped-storage power station. The calculations show that when fully enclosed anti-seepage and drainage measures are taken for the underground powerhouse, the maximum head of water is 98 m, which is lower than the critical water pressure of seepage failure for the argillaceous dolomite. Therefore, no seepage failure will occur. The results provide a scientific basis for the anti-seepage and drainage design of the underground powerhouse area.

Assessing the Connectivity of a Regional Fractured Aquifer Based on a Hydraulic Conductivity Field Reversed by Multi-Well Pumping Tests and Numerical Groundwater Flow Modeling

Assessing the Connectivity of a Regional Fractured Aquifer Based on a Hydraulic Conductivity Field Reversed by Multi-Well Pumping Tests and Numerical Groundwater Flow Modeling

Effects of a Gravel Pit Lake on Groundwater Hydrodynamic

In Europe, 1132 Mt of sand and gravel were mined in 2019, which causes major changes to the hydrogeological cycle. Such changes may lead to significantly raised or lowered groundwater levels. Therefore, the aggregate sector has to ensure that impacts on existing environmental and water infrastructures are kept to a minimum in the post-mining phase. Such risk assessments are often made by empirical methods, which are based on assumptions that do not meet real aquifer conditions. To investigate this effect, predictions by empirical and numerical methods about hydraulic head changes caused by a pit lake were compared. Wrobel’s equation, which is based on Sichardt’s equation, was used as the empirical method, while a numerical groundwater flow model has been solved by means of the finite-element method in FEFLOW. The empirical method provides significantly smaller ranges of increased/decreased groundwater levels caused by the gravel pit lake as the numerical method. The underestimation of the empirical results was related to the finding that field measurements during pumping tests show a larger extent of groundwater drawdown than calculations with the Sichardt’s equation. Simplifications of the 2D model approach have been evaluated against hydraulic head changes derived from a 3D groundwater model. Our results clearly show that the faster and cheaper empirical method—Wrobel’s equation, which is often preferred over the more expensive and time-consuming numerical method, underestimates the drawdown area. This is especially critical when the assignment of mining permits is based on such computations. Therefore, we recommend using numerical models in the pre-mining phase to accurately compute the extent of a gravel/sand excavation’s impacts on hydraulic head and hence more effective protection of groundwater and other related environmental systems.

Numerical Modelling of Groundwater Flow and Contamination Transport (TDS) in Kaniqurzala Area- Central Basin – Erbil City- Kurdistan - Iraq

Numerical Modelling of Groundwater Flow and Contamination Transport (TDS) in Kaniqurzala Area- Central Basin – Erbil City- Kurdistan - Iraq

Integrated Three‐Dimensional Geological and Numerical Groundwater Model Development

AbstractConceptual site model (CSM) development is a foundational step for numerical groundwater model construction. Tools for reviewing and combining disparate data to develop a CSM have been furthered through software that efficiently integrates multiple data streams, especially when a large amount of information is available at complex sites. Three‐dimensional (3D) interfaces for data visualization and geological modeling are becoming commonly used to compile data, support development of a robust CSM, and create effective visuals that enable project stakeholders to understand site complexity. This article describes a 3D interface workflow applied to a geological and numerical groundwater modeling case study where numerical modeling (i.e., MODFLOW family of codes) was used to assess remedy effectiveness. Specifically, the numerical model was used to overview an in situ solidification (ISS) application, where soil mixing will be conducted with cement, for possible adverse effects to an existing groundwater treatment system. 3D geological modeling allowed for the existing numerical groundwater flow model to be reconstructed to include previously overlooked data on buried infrastructure, which led to a detailed depiction of high‐permeability channels and a more accurate representation of site conditions. Modeling showed that the ISS remedy would not adversely impact the existing groundwater treatment system. Integration of existing data within a 3D interface led to hydrogeological insights that would have been difficult to obtain from two‐dimensional visualization of information, improved numerical modeling efficiency, and informed site remedy decisions. Additionally, the 3D interface created effective visuals that increased internal and external project team conceptual understanding.