Significance of groundwater flux on contaminant concentration and mass discharge in the nonaqueous phase liquid (NAPL) contaminated zone

•We develop models of how groundwater flux variations impact NAPL dynamics.•Any forms of temporal variations of groundwater flux can be incorporated.•Flux affects contaminant mass discharge reduction (CMDR) and mass reduction (MR).•Models based on variable groundwater flux conditions compare better with field data.•The study can guide field campaign efforts in sites contaminated by NAPLs.

Groundwater flowing through residual nonaqueous phase liquid (NAPL) source zone will cause NAPL dissolution and generate large contaminant plume. The use of contaminant mass discharge (CMD) measurements in addition to NAPL aqueous phase concentration to characterize site conditions and assess remediation performance is becoming popular. In this study, we developed new and generic numerical models to investigate the significance of groundwater flux temporal variations on the NAPL source dynamics. The developed models can accommodate any temporal variations of groundwater flux in the source zone. We examined the various features of groundwater flux using a few selected functional forms of linear increase/decrease, gradual smooth increase/decrease, and periodic fluctuations with a general trend. Groundwater flux temporal variations have more pronounced effects on the contaminant mass discharge dynamics than the aqueous concentration. If the groundwater flux initially increases, then the reduction in contaminant mass discharge (CMDR) vs. NAPL mass reduction (MR) relationship is mainly downward concave. If the groundwater flux initially decreases, then CMDR vs. MR relationship is mainly upward convex. If the groundwater flux variations are periodic, the CMDR vs. MR relationship tends to also have periodic variations ranging from upward convex to downward concave. Eventually, however, the CMDR vs. MR relationship approaches 1:1 when majority of the NAPL mass becomes depleted.