A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

Dera Ismail Khan (DIK) is situated in the Lower Indus Basin of Pakistan. The land use has been changed in the canal command area due to irrigation activities near the Indus River. To check the current status and predict the groundwater levels in the area, the unconfined aquifer has been simulated in Visual MODFLOW for a period of 35 years, i.e., from 1985 to 2020. The 2900-km² area has been modeled with a grid of 500 by 500 m and the depth set to 100 m. The aquifer in the study area has been divided vertically and laterally into three and ten zones, respectively, for the characterization. Water wells and streams were used as the sinks and hydrologic boundaries, respectively. The model was successfully calibrated in steady and the non-steady state. The simulation revealed that the whole simulation can be divided into two phases, i.e., before and after the construction of the Chashma Right Bank Canal (CRBC), whereas the results were summarized in the form of water table depth maps and groundwater budget calculations. To determine the groundwater sustainability, a conjunctive use scenario has been employed to simulate the aquifer dynamics till 2020. The simulation revealed incremental drawdowns till the end.     

Erratum to: A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

The Gour Oumelalen area exposes Paleoproterozoic (1.9 Ga) marbles and calc-silicate granulites. Some marbles show a specific mineralogy characterized by the presence of a highly aluminous clinopyroxene with Al₂O₃ content exceeding16 wt%. This clinopyroxene shows a marked zoning with a hedenbergitic core rimmed by fassaite. Phase relations are expressed by spectacular reaction textures in calc-silicate granulites as Opx + Cpx + Pl + H₂O  == > Grt + Qtz ± Am and Cpx + Ilm + Pl ± Mt = > Grt + Qtz + Spn. In olivine-bearing marbles, clinopyroxene and dolomite occur around olivine and calcite. According to thermodynamic modeling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (NCFMASHTO) system, the first stage of metamorphism is located at 800 °C and 6–7 kbar, which is followed by an increase in pressure at 9–10 kbar and 800 °C and an isobaric cooling at 680 °C. The deduced anticlockwise P–T path is consistent with a granulitic metamorphism occurring in an active continental margin context.