Numerical Investigation into the Evolution of Groundwater Flow and Solute Transport in the Eastern Qaidam Basin since the Last Glacial Period

Abstract

A complete understanding of groundwater circulation as well as the transport and distribution of solutes in arid-semiarid basin regions is a prerequisite for the safe use of groundwater resources. The distributions of the groundwater flow systema and solutes are affected by the basin morphology, lithology, and climate variations; therefore, they can change over geologic time. In this study, we performed a case study of the Qaidam Basin in the northeastern Tibetan Plateau, in which we utilized reactive solute transport simulations to build a numerical model in TOUGHREACT for a typical section of the eastern Qaidam Basin since the last glacial period. The results show that the groundwater in the eastern Qaidam Basin developed into a three-level groundwater flow system and that the seepage velocity of the local water flow system is significantly higher than that of the intermediate and regional water flow systems. Although groundwater in the discharge region has been continuously concentrated and enriched since the last glacial period, the distributions of the groundwater flow system and solutes have been greatly affected by climate variations. During warm periods, the centres of groundwater discharge and solute concentration shifted to areas with more groundwater recharge; in contrast, both centres shifted to the central basin during drought periods. The groundwater in the basin mainly contains Na+ and Cl− ions, which vary significantly from the recharge region to the discharge region. Evaporation of groundwater results in increases in the concentrations of most of the components except HCO3−. The groundwater in the discharge region is currently in the stage of carbonate precipitation and is far from gypsum and halite precipitation.