The aim was to study density-driven groundwater flow and analyse groundwater mixing because of seasonal changes in groundwater temperature. Here, density-driven convection in groundwater was studied by numerical simulations in a subarctic climate, i.e. where the water temperature was <4 °C. The effects of soil permeability and groundwater temperature (i.e. viscosity and density) were determined. The influence of impermeable obstacles in otherwise homogeneous ground was also studied. An initial disturbance in the form of a horizontal groundwater flow was necessary to start the convection. Transient solutions describe the development of convective cells in the groundwater and it took 22 days before fully developed convection patterns were formed. The thermal convection reached a maximum depth of 1.0 m in soil of low permeability (2.71 · 10⁻⁹ m²). At groundwater temperature close to its density maximum (4 °C), the physical size (in m) of the convection cells was reduced. Small stones or frost lenses in the ground slightly affect the convective flow, while larger obstacles change the size and shape of the convection cells. Performed simulations show that “seasonal groundwater turnover” occurs. This knowledge may be useful in the prevention of nutrient leakage to underlying groundwater from soils, especially in agricultural areas where no natural vertical groundwater flow is evident. An application in northern Sweden is discussed.

Temperature, discharge, and stable isotope ratios of five karst springs in a mountainous area of Zigui County, Hubei Province, Central China, were analyzed. The purpose was to illustrate the heat exchanges linked to circulation depth in the exposed karst water systems through the development of a method for estimating heat input and heat flux during a rainstorm. Meteorological water in the study area conformed to a local meteoric water line (δD = 8.37 δ¹⁸O + 12.99) with a mean δ¹⁸O elevation gradient of −4.0‰ km⁻¹, which was used to estimate mean circulation depths of 209–686 m. The mean spring temperatures defined a vertical gradient of −5.4 °C km⁻¹, which resembled that of the stable atmosphere of the Earth, indicating that the thermal response patterns are mainly controlled by surface air temperature. Thermal convection after rainfall events dominated heat exchange between baseflow and recharge water, leading to a warmer and colder recharge during summer and winter, respectively, whereas thermal conduction dominated the heat exchange only between groundwater, surrounding geology, and the interface air under a condition of no rainfall, resulting in only small temperature variations of the baseflow. Successful application of the method for estimating heat exchange showed that the characteristics of shallow circulation, strong karstification, and well-developed epikarst readily allowed disruption of the thermal balance of the Yuquandong system, resulting in a poor heat regulation capacity, a larger variation of heat input, a lower mean heat flux, and lower baseflow temperatures compared to those of the Dayuquan system.