Backward probabilistic model of groundwater contamination in non-uniform and transient flow

Backward location and travel time probabilities, which provide information about the former location of contamination in an aquifer, can be used to identify unknown contamination sources. Backward location probability describes the possible upgradient positions of contamination at a known time in the past, and backward travel time probability describes the time required for contamination to travel from a known upgradient location to an observation point. These probabilities are related to adjoint states of resident concentration, and their governing equation is the adjoint of a forward contaminant transport model. Using adjoint theory to obtain the appropriate governing equation, we extend the backward probability model for conservative solutes to more general non-uniform and transient flow fields. In particular, we address three important extensions, spatially-varying porosity, transient flow and temporally-varying porosity, and internal distributed sources and sinks of solute and water. For the first time we learn that forward and backward location and travel time probabilities are not necessarily equivalent to adjoint states, but are related to them. The extensions are illustrated using a vertically-integrated groundwater model, creating transient flow by a step change in pumping and using areal recharge as an internal distributed source. Both the movement and spread of probabilities are affected. With internal sources of water, there are two interpretations of backward probability, depending on whether or not the source of water is also a source of solute. The results demonstrate how the backward probability model can be applied to other, perhaps more important, non-uniform and transient flow conditions, with time- and space-varying water storage, such as time-varying pumping or unsaturated (or saturated–unsaturated) flow and transport with spatially- and temporally-varying moisture content.