H21F-1450
A Numerical Model of Deuterium and Oxygen-18 Diffusion in the Confined Lower Wilcox Aquifer of the Lower Mississippi Valley (USA)

Tuesday, 15 December 2015
Poster Hall (Moscone South)
Benjamin James Currens1, Audrey H Sawyer2, Alan E Fryar1, Thomas M Parris3 and Junfeng Zhu1, (1)University of Kentucky, Lexington, KY, United States, (2)Ohio State University Main Campus, Columbus, OH, United States, (3)Kentucky Geologic Survey, Lexington, KY, United States
Abstract:
Deuterium and oxygen-18 are routinely used with noble gases and radioisotopes (e.g., 2H, 14C, 36Cl) to infer climate during groundwater recharge. However, diffusion of 2H and 18O between a confined aquifer and bounding aquitards could alter total isotope concentrations and the inferred temperature during recharge if groundwater flow is sufficiently slow. Hendry and Schwartz (WRR 24(10), 1988) explained anomalous 2H and 18O enrichment in the Milk River aquifer of Alberta by analytically modeling isotope diffusion between the lower bounding aquitard and the aquifer. Haile (PhD dissertation, U. Kentucky, 2011) inferred the same mechanism to explain 2H and 18O enrichment along a flowpath in the confined Lower Wilcox aquifer of the northern Gulf Coastal Plain in Missouri and Arkansas. Based on the geologic and hydraulic properties of the Lower Wilcox aquifer, a numerical model has been constructed to determine how diffusion may influence 2H and 18O concentrations in regional aquifers with residence times on the order of 104 to 105 years. The model combines solutions for a 1D forward-in-time, finite-difference groundwater flow equation with an explicit-implicit Crank-Nicholson algorithm for advection and diffusion to solve for flow velocity and isotope concentration. Initial results are consistent with the analytical solution of Hendry and Schwartz (1988), indicating diffusion as a means of isotopic enrichment along regional groundwater flowpaths.