Modeling the subglacial hydrology of the James Lobe of the Laurentide Ice Sheet

To investigate the drainage conditions that might be expected to develop beneath soft-bedded ice sheets, we modeled the subglacial hydrology of the James Lobe of the Laurentide Ice Sheet from Hudson Bay to the Missouri River. Simulations suggest the James Lobe had little effect on regional groundwater flow because the poorly conductive Upper-Cretaceous shale that occupies the upper layer of the bedrock would have functioned as a regional aquitard. This implies that general northward groundwater flow out of the Williston Basin has likely persisted throughout the Quaternary. Moreover, the simulations indicate that the regional aquifer system could not have drained even the minimum amount of basal meltwater that might have been produced from at the glacier bed. Therefore, excess drainage must have occurred by some sort of channelized drainage network at the ice–till interface. Using a regional groundwater model to determine the hydraulic conductivity for an equivalent porous medium in a 1-m thick zone between the ice and underlying sediment, and assuming conduit dimensions from previous theoretical work, we use a theoretical karst aquifer analog as a heuristic approach to estimate the spacing of subglacial conduits that would have been required at the ice–till interface to evacuate the minimum water flux. Results suggest that for conduits assumed to be on the order of a tenth of a meter deep and up to a meter wide, inter-conduit spacing must be on the order of tens–hundreds of meters apart to maintain basal water pressures below the ice overburden pressure while evacuating the hypothesized minimum meltwater flux.