Stable hydrogen and oxygen isotopes reveal aperiodic non-river evaporative solute enrichment in the solute cycling of rivers in arid watersheds

• Published in: The Science of the total environment Vol. 856; p. 159113
• Main Authors: Letshele, Kesego P., Atekwana, Eliot A., Molwalefhe, Loago, Ramatlapeng, Goabaone J., Masamba, Wellington R.L.
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 15-01-2023
• Subjects: Botswana; d-Excess; Endorheic river basin; Hydrologic perturbations; Okavango Delta; Total dissolved ions; Hydrologic perturbations; Endorheic river basin; d-Excess; Botswana; Okavango Delta; Total dissolved ions
• ISSN: 0048-9697; 1879-1026
• DOI: 10.1016/j.scitotenv.2022.159113

We investigated the spatial and temporal variations of the stable isotope composition of hydrogen (δD) and oxygen (δ18O) and the total dissolved ions (TDI) concentrations in the Okavango River in the middle Kalahari Desert. We aimed to elucidate the role of evaporation in controlling river solute enrichment from samples collected at a one- to two-month frequency from nine stations along a ∼460 km river transect for one year. We found that the δD and δ18O composition and the TDI concentrations increased downriver. Seasonal increases in the δD and δ18O composition and TDI concentrations during the hot, rainy season were subdued or decreased during the cool, dry season from pulse flooding. The δD and δ18O values of the samples plot along the Okavango Delta Evaporation Line consistent with evaporation. The effect of evaporation during river transit was related to the mean δD (δD = 0.07*River distance (km) − 37.9; R2 = 0.98) and mean d-excess (d-excess = −0.04*River distance (km) + 9.9; R2 = 0.94). The effect of evaporation on the river solute behavior is characterized by the mean d-excess and TDI concentrations (d-excess = −0.29*TDI (mg/L) + 15.0; R2 = 0.97). Some samples from this study and those compiled from published studies plot at greater than one sigma standard deviation above and below the mean TDI concentration vs. mean d-excess regression model line. We use these marked deviations from the mean TDI concentration vs. the mean d-excess regression model to suggest that additional solutes from river-floodplain-wetland-island interaction driven by pulse flooding are delivered into the river. While our findings support an evaporation-dominated solute enrichment during river transit at the seasonal scale, we conclude that intermittent hydrology (pulse flooding) plays an important role in the sub-seasonal spatiotemporal behavior of solutes in rivers in arid watersheds, which must be considered in solute cycling models.
[Display omitted]
•The δD, δ18O and total dissolved ions (TDI) were investigated along the Okavango River in arid Botswana.•The δ18O vs. δD of river samples plot along the Okavango Delta Evaporation Line due to evaporation.•Downriver d-excess was depleted by −0.04 ‰ km−1 (R2 = 0.94) and correlates to −0.29 mean TDI (R2 = 0.97).•A mean TDI-d-excess model supports evaporation and reveals aperiodic solute input into the river.•The aperiodic solute input driven by annual pulse flooding is spatially and temporally variable.