Numerical analysis of buoyancy effects during the dissolution and transport of oxygenated gasoline in groundwater

Abstract

Dissolution of oxygenated gasoline, as well as buoyancy-driven groundwater flow and transport of the multicomponent dissolved phase plumes, is simulated numerically in three dimensions. The simulations are based on a field experiment described by Mocanu (2007) in which three oxygenated gasoline sources were emplaced as nonaqueous phase liquids (NAPLs) below the water table of the shallow sand aquifer at Canadian Forces Base Borden, Ontario. The sources were composed of an ethanol-free gasoline mixture spiked with 9.8% methyl tert-butyl ether and 0.2% tert-butyl alcohol (GMT-E0), a gasoline with 10% ethanol (E10), and a source with 95% ethanol (E95). The numerical model includes dissolution of gasoline as a NAPL, density-dependent groundwater flow, advective-dispersive transport of the dissolved components, and ethanol cosolvency and degradation. Buoyancy effects in the dissolved plumes were compared under a homogeneous hydraulic conductivity field as well as with five realizations of spatially correlated random fields representing the Borden aquifer. The simulations showed that buoyancy was most significant in the E95 source plumes within the homogeneous system, having induced after 150 days a net upward displacement of the local peak concentrations for all but the least soluble component of approximately 1.5 m. The peak rise in ethanol from the GMT-E0 and E10 plumes was about 0.6 m. The results highlight the importance of shallow monitoring wells when monitoring high oxygenate fraction gasoline spills in groundwater and have implications for assessing mass fluxes and biodegradation rates. Copyright 2008 by the American Geophysical Union.