2014 GSA Annual Meeting in Vancouver, British Columbia (19–22 October 2014)

Paper No. 337-1
Presentation Time: 1:00 PM


KHAN, Mahfuzur R.1, KONESHLOO, Mohammad2, KNAPPETT, Peter S.K.3, MOZUMDER, M. Rajib Hassan4, MAILLOUX, Brian5, AHMED, Kazi Matin6, BOSTICK, Benjamin C.7, HARVEY, Charles F.8, VAN GEEN, Alexander4 and MICHAEL, Holly A.1, (1)Department of Geological Sciences, University of Delaware, 255 Academy Street, Newark, DE 19716, (2)Department of Geological Sciences, University of Delaware, Newark, DE 19716, (3)Geology and Geophysics, Texas A&M University, College Station, TX 77843, (4)Lamont-Doherty Earth Observatory, Columbia University, 61 Route 9w, Palisades, NY 10964, (5)Environmental Sciences, Barnard College, 3009 Broadway, New York, NY 10027, (6)Department of Geology, University of Dhaka, Dhaka, Dhaka 1000, Bangladesh, (7)Lamont-Doherty Earth Observatory, Columbia University, 61 Route 9W, PO Box 1000, Palisades, NY 19064, (8)Dept of Civil and Environmental Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA 02139, mahfuz@udel.edu

Groundwater pumping beneath Dhaka, a mega-city in Bangladesh, has created an extensive drawdown cone. Outside the city, shallow groundwater is often contaminated with dangerous concentrations of arsenic (As) and deep groundwater is considered a reliable source of As-safe drinking water. Unfortunately, vertical head gradients induced by Dhaka pumping may threaten the deep aquifer with downward migration of shallow groundwater elevated in As and/or dissolved organic carbon capable of causing in-situ As release. This study explores the role of geologic heterogeneity in this downward flow. Five realizations of heterogeneous hydraulic conductivity were simulated with sequential indicator simulation conditioned to 400 borehole logs within the study area. A regional MODFLOW model with high resolution in the Dhaka area was developed for flow simulation through both heterogeneous and equivalent homogeneous systems. Advective particle tracking was used to determine the transport time and flowpaths of shallow groundwater. Results indicate that deep Dhaka pumping has induced vertical groundwater flow from the rivers and shallow aquifers in the surrounding regions. Although the overall shape and extent of the Dhaka pumping cone of depression is the same in both homogeneous and heterogeneous systems, there is substantial local variability in the vertical hydraulic gradient and travel times in heterogeneous systems. Both the minimum and maximum travel times in heterogeneous systems differ by an order of magnitude compared to that of equivalent homogeneous systems. Moreover, in the heterogeneous systems, although high gradients exist across clays, vertical flow is largely impeded; travel times to deep aquifer zones are thousands of years. In contrast, areas with thick sandy sediments showing near-zero vertical gradients are vulnerable to vertical flow of shallow groundwater in the deeper part of the aquifer system within a few decades. The results have implications for managing the deep groundwater resource in the Bengal Delta and in similar areas with contaminated shallow groundwater.