Abstract
Abstract. The stable isotope compositions of soil water (δ2H and δ18O) carry important information about the prevailing soil hydrological conditions and for constraining ecosystem water budgets. However, they are highly dynamic, especially during and after precipitation events. In this study, we present an application of a method based on gas-permeable tubing and isotope-specific infrared laser absorption spectroscopy for in situ determination of soil water δ2H and δ18O. We conducted a laboratory experiment where a sand column was initially saturated, exposed to evaporation for a period of 290 days, and finally rewatered. Soil water vapor δ2H and δ18O were measured daily at each of eight available depths. Soil liquid water δ2H and δ18O were inferred from those of the vapor considering thermodynamic equilibrium between liquid and vapor phases in the soil. The experimental setup allowed for following the evolution of soil water δ2H and δ18O profiles with a daily temporal resolution. As the soil dried, we could also show for the first time the increasing influence of the isotopically depleted ambient water vapor on the isotopically enriched liquid water close to the soil surface (i.e., atmospheric invasion). Rewatering at the end of the experiment led to instantaneous resetting of the stable isotope profiles, which could be closely followed with the new method. From simple soil δ2H and δ18O gradients calculations, we showed that the gathered data allowed one to determinate the depth of the evaporation front (EF) and how it receded into the soil over time. It was inferred that after 290 days under the prevailing experimental conditions, the EF had moved down to an approximate depth of −0.06 m. Finally, data were used to calculate the slopes of the evaporation lines and test the formulation for kinetic isotope effects. A very good agreement was found between measured and simulated values (Nash and Sutcliffe efficiency, NSE = 0.92) during the first half of the experiment, i.e., until the EF reached a depth of −0.04 m. From this point, calculated kinetic effects associated with the transport of isotopologues in the soil surface air layer above the EF provided slopes lower than observed. Finally, values of kinetic isotope effects that provided the best model-to-data fit (NSE > 0.9) were obtained from inverse modeling, highlighting uncertainties associated with the determinations of isotope kinetic fractionation and soil relative humidity at the EF.
Highlights
Stable isotopologues of water, namely, 1H2H16O and 1H128O, are powerful tools used in a wide range of research disciplines at different and complementary temporal and spatial scales
From simple soil δ2H and δ18O gradients calculations, we showed that the gathered data allowed one to determinate the depth of the evaporation front (EF) and how it receded into the soil over time
For an unsaturated soil, assuming in a first approximation that isotope movement occurs in the vapor phase above the soil evaporation front (EF) and strictly in the liquid phase below it, the maximal soil water δ2H and δ18O values are no longer observed at the surface but at the depth of the EF
Summary
1H2H16O and 1H128O, are powerful tools used in a wide range of research disciplines at different and complementary temporal and spatial scales. Braud et al (2005), Haverd and Cuntz (2010), Rothfuss et al (2012), Singleton et al (2004) and Sutanto et al (2012) implemented the description of the transport of 1H2H16O and 1H128O in physically based soil–vegetation– atmosphere transfer (SVAT) models (HYDRUS 1D, SiSPATIsotope, soil–litter–iso, TOUGHREACT). In these models, movement of soil 1H2H16O and 1H128O occurs in both phases below and above the EF, and heat and water transports are properly coupled
Sign-in/Register to view highlights & summary Sign-in/Register to access full text options 
Free) 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
