The role of dew and radiation fog inputs in the local water cycling of a temperate grassland during dry spells in central Europe
- 1Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
- 2Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland
Correspondence: Yafei Li (email@example.com) and Werner Eugster (firstname.lastname@example.org)
During dry spells, non-rainfall water (hereafter NRW) mostly formed from dew and fog potentially plays an increasingly important role in temperate grassland ecosystems with ongoing global warming. Dew and radiation fog occur in combination during clear and calm nights, and both use ambient water vapor as a source. Research on the combined mechanisms involved in NRW inputs to ecosystems is rare, and distillation of water vapor from the soil as a NRW input pathway for dew formation has hardly been studied. Furthermore, eddy covariance (EC) measurements are associated with large uncertainties on clear, calm nights when dew and radiation fog occur. The aim of this paper is thus to use stable isotopes as tracers to investigate the different NRW input pathways into a temperate Swiss grassland at Chamau during dry spells in summer 2018. Stable isotopes provide additional information on the pathways from water vapor to liquid water (dew and fog) that cannot be measured otherwise. We measured the isotopic composition (δ18O, δ2H, and ) of ambient water vapor, NRW droplets on leaf surfaces, and soil moisture and combined them with EC and meteorological observations during one dew-only and two combined dew and radiation fog events. The ambient water vapor d was found to be strongly linked with local surface relative humidity (), highlighting the dominant role of local moisture as a source for ambient water vapor in the synoptic context of the studied dry spells. Detailed observations of the temporal evolution of the ambient water vapor and foliage NRW isotopic signals suggest two different NRW input pathways: (1) the downward pathway through the condensation of ambient water vapor and (2) the upward pathway through the distillation of water vapor from soil onto foliage. We employed a simple two-end-member mixing model using δ18O and δ2H to quantify the NRW inputs from these two different sources. With this approach, we found that distillation contributed 9 %–42 % to the total foliage NRW, which compares well with estimates derived from a near-surface vertical temperature gradient method proposed by Monteith in 1957. The dew and radiation fog potentially produced 0.17–0.54 mm d−1 NRW gain on foliage, thereby constituting a non-negligible water flux to the canopy, as compared to the evapotranspiration of 2.7 mm d−1. Our results thus underline the importance of NRW inputs to temperate grasslands during dry spells and reveal the complexity of the local water cycle in such conditions, including different pathways of dew and radiation fog water inputs.
The role of dew and fog inputs in the hydrological cycle is well understood in desert areas, where rainfall totals are small (Malek et al., 1999; Jacobs et al., 2002; Kidron et al., 2002; Agam and Berliner, 2006; del Prado and Sancho, 2007; Pan et al., 2010; Ucles et al., 2013; McHugh et al., 2015). Such water inputs are, however, mostly neglected in regions where average rainfall is abundant, and thus the expected water gains from dew of up to 0.7–0.8 L m−2 d−1 during nights with perfect clear-sky conditions (Beysens, 2018), or fog providing on the order of 8.5 L m−2 d−1 in tropical montane cloud forests (Bruijnzeel et al., 2006) appear to be small and negligible in comparison to average precipitation rates. However, during dry spells, especially during the warm season when daily evapotranspiration rates are high, it can be expected that, although small, non-rainfall water (hereafter NRW) inputs from various sources (see below) may become essential for the vegetation to alleviate stress (Tuller and Chilton, 1973). This may even be the case in temperate climates, where average annual precipitation typically balances or exceeds actual annual evapotranspiration. Grasslands tend to be the first to suffer from prolonged dry spells and droughts (Wolf et al., 2013). Here we investigate the small-scale processes of how fog and dew water influence the water cycling over a grassland at a central European temperate climate site during representative warm-season nights.
Rainfall measurements with conventional rain gauges collect liquid and solid precipitation (Glickman and Zenk, 2000) and thus the vast amount of above-ground water entering the vegetation canopy in wet climates, but in temperate and even more pronounced in dry climates, some important components of the hydrological cycle are missed, e.g., NRW inputs. NRW inputs include a number of components: (1) dew formation (Monteith, 1957); (2) fog deposition (Dawson, 1998); (3) water vapor adsorption (Agam and Berliner, 2006); (4) rime ice deposition (Hindman et al., 1983); (5) hoar frost (Monteith and Unsworth, 2013); and (6) guttation (Long, 1955). During extended periods without rainfall, it is well known that mainly dew and fog (out of the long list of NRW components) are essential water sources for plants in (1) arid and semi-arid regions (Malek et al., 1999; Jacobs et al., 2002; Kidron et al., 2002; Agam and Berliner, 2006; del Prado and Sancho, 2007; Kidron and Temina, 2013; Ucles et al., 2013; He and Richards, 2015; McHugh et al., 2015; Tomaszkiewicz et al., 2017); (2) Mediterranean coastal regions (Beysens et al., 2007); (3) temperate ecosystems (Jacobs et al., 2006); and (4) tropical climates (Clus et al., 2008). On clear, calm nights when dew and radiation fog occur, the atmospheric boundary layer becomes stably stratified, leading to a shallow stable nocturnal boundary layer (hereafter NBL) with a depth on the order of no more than 50–100 m (Garratt, 1992). Dew and radiation fog occur at the bottom of the NBL (Stull, 1988; Garratt, 1992; Oke, 2002; Monteith and Unsworth, 2013). Both dew and radiation fog are formed due to the cooling of the Earth's surface after sunset by long-wave radiation losses on clear nights (Oke, 2002). This radiative cooling is a process due to which a body loses heat by long-wave thermal radiation, whereby its surface cools below the dew point of the adjacent air. Under such conditions, dew can form on plant surfaces while fog forms on activated aerosol particles in the near-surface atmosphere.
NRW inputs contribute to the water budget across many ecosystems including croplands (Atzema et al., 1990; Wen et al., 2012; He and Richards, 2015; Meng and Wen, 2016; Tomaszkiewicz et al., 2017), grasslands (Jacobs et al., 2006; Wen et al., 2012; He and Richards, 2015), and forests (Fritschen and Doraiswamy, 1973; Dawson, 1998; Lai and Ehleringer, 2011; Hiatt et al., 2012; Berkelhammer et al., 2013). As compared to forests, grasslands present favorable conditions for dew and radiation fog formation: (1) a cooler surface due to a higher albedo and thus lower net solar radiation input (Moore, 1976), as well as higher evapotranspiration (Kelliher et al., 1993; Williams et al., 2012); (2) weaker aerosol particle deposition due to shorter roughness length of grasslands (Gallagher et al., 2002), and thus more aerosol particles remaining in the near-surface atmosphere, which consequently results in better conditions for radiation fog formation over grasslands. From the perspective of ecological functioning, small amounts of NRW inputs have a more important influence on grasslands than forests because of a reduced capability to increase crop water use efficiency (WUE), defined as gross carbon uptake per unit water lost, when water availability is low (Wolf et al., 2013), but also due to lower soil water availability and shallower rooting depth in grasslands. At the beginning of drought stress in ecosystems, forests increase their WUE by closing their stomata, which increases stomatal resistance and thus reduces evapotranspiration, while grasslands maintain their evapotranspiration as long as the soil water is available to supply the evaporative demand (e.g., Wolf et al., 2013). Therefore, grasslands are more prone to suffer from soil water scarcity. In addition, as opposed to the deep-rooted systems for forest plants, grassland plants take up water from the topsoil, where scarcity of soil water occurs more frequently in the absence of precipitation; therefore grasslands tend to anticipate lower soil water availability compared to forests.
Ambient water vapor is the main vapor source for both dew and radiation fog; therefore, dew and radiation fog usually occur in combination. Because of the variability of temperature and humidity conditions, a single NRW night may transit from dew only to intermittent dew and radiation fog in combination. Before the atmospheric humidity reaches saturation at the standard measurement height at 2 m a.g.l., dew can only form if the surface temperature drops below air temperature. When the ambient water vapor reaches saturation or even super-saturation, dew and radiation fog can form intermittently. Kaseke et al. (2017) used hydrogen and oxygen stable isotope regression to separate the different types of dew and fog, but they focused on dew and fog events separately. Research that focusses on relevant phase change processes during dew and radiation fog in combination is thus rare.
The moisture movement in the soil–plant–atmosphere continuum has been well understood by the eddy covariance (hereafter EC) technique, but the reliability of the method suffers during nighttime with weak turbulence (Berkelhammer et al., 2013) when dew and radiation fog occur. In principle, downward water vapor flux measured by EC should provide a quantitative estimate of dew formation on the vegetation surface (termed “phantom dew” by Gay et al. 1996). The results by Jacobs et al. (2006), however, showed that dew formation quantified by EC was less than one-third of the dew amount measured by a lysimeter (estimated from Fig. 1 in their paper). Moreover, katabatic cold-air drainage flows in non-flat topography lead to advective fluxes that are not directly captured by EC measurements (e.g., Eugster and Siegrist 2000, and Sun et al., 2006), which typically leads to a gap in the local energy budget , with Rn the net all-wave radiation, H the sensible heat flux, LE the latent heat flux, and ΔQ the energy flux to close the budget (see also Wilson et al., 2002, and Franssen et al., 2010), which makes estimates of dew formation during calm nights highly uncertain and unreliable if ΔQ≠0, and they are thus not further addressed in this paper.
Monteith (1957) identified two input pathways for dew formation: (1) the downward pathway through the condensation of ambient water vapor on the plants and/or on soil surface and (2) the upward pathway through distillation of water vapor from soil onto plant surfaces. Soil vapor diffusion from the soil to the atmosphere is driven by the temperature gradient between the soil and the atmosphere and between different depths of the soil (Monteith, 1957; Oke, 1970). The temperature gradient generally reaches a maximum at the soil–atmosphere interface, where the soil surface is roughly 2–5 ∘C warmer than the adjacent air at 1 cm a.g.l. for short grass cover according to Monteith (1957). The diffusing soil vapor can therefore condense onto cooler foliage. Since Monteith (1957) quantified the downward and upward components of dew formation by absorbing NRW on foliage with filter paper, research has rarely focused on distinguishing these two pathways of dew formation. Furthermore, Monteith (1957) distinguished the two pathways by collecting NRW on separate nights when only one or the other of the two pathways was assumed to occur. In Monteith (1957), distillation of water vapor from soil as one component of NRW was quantified on very calm nights with a 2 m wind speed (hereafter u2m) of less than 0.5 m s−1, whereas the maximum NRW condensing from ambient water was assumed to occur on slightly windy nights with u2m in the range of 2–3 m s−1. However, for clear, calm nights with u2m between 0.5 and 2 m s−1, condensation of ambient water vapor and soil-diffusing vapor can occur in combination, with NRW on the foliage being a mix from these two pathways. Research focusing on distinguishing and quantifying the ratio of these two NRW components, i.e., NRW from ambient water vapor and distillation, is scarce.
When the condensation of ambient water vapor and distillation occur simultaneously, the partitioning of NRW into these two components becomes difficult because there is no direct measurement possible to quantify distillation amounts. Hydrometric approaches, e.g., using lysimeters, can easily quantify the condensation amount of ambient water vapor but cannot quantify the distillation amount, if the water vapor condensing on the above-ground parts (e.g., leaf surfaces) of the lysimeter stems from the below-ground part (soil) of the same lysimeter without a net change in lysimeter weight. Monteith (1957) provided the equations to calculate the distillation rate through measuring the soil surface temperature and air temperature at 1 cm a.g.l. (see Sect. 3.2.5), which he compared with filter paper measurements and interpreted that the “agreement was reasonable”, with a mean ratio of observed vs. calculated distillation of 0.76 (i.e., calculated distillation was ≈32 % higher than the observed distillation; see Monteith, 1957). The disagreement according to Monteith (1957) is not only related to the unknown collection efficiency of the filter paper he used but may have arisen from errors or uncertainties in the following three assumptions: (1) the assumption of purely molecular and thus nonturbulent transfer, (2) the assumption of linear (not curvilinear or exponential) temperature gradient, and (3) the assumption of saturation at the soil surface that air in direct contact with the soil may be undersaturated if the 1 cm temperature is lower than the soil surface temperature. To overcome the above-mentioned challenges of quantifying distillation with traditional methods such as EC, filter paper or the vertical temperature gradient method by Monteith (1957), a useful approach to quantify the ratio between condensation and distillation is the use of stable isotopes: NRW inputs from ambient water vapor and from distillation carry different isotopic signatures due to their different sources, i.e., the atmosphere and the soil moisture respectively. Therefore, a two-end-member mixing model using stable isotopes in water (Keeling, 1958; Dawson, 1998; Phillips et al., 2005) can be employed to quantify the individual contributions of these two sources (see details in Sect. 3.2.4).
Our aim was thus to (1) investigate the isotopic fractionations during dew-only and dew–fog combined events and (2) estimate the contribution of NRW from atmospheric vapor and from soil-diffusing vapor. We carried out three 24 h observation campaigns during summer 2018 using stable isotopes combined with EC and meteorological measurements to characterize the meteorological conditions, to analyze the isotope fractionation of dew and radiation fog formation, to quantify the NRW contribution from ambient water vapor and soil-diffusing vapor, and to explore the potential role of dew and radiation fog during dry spells in temperate grasslands.
2.1 Hydrogen and oxygen isotopes
Hydrogen and oxygen stable isotopes are a useful research tool to investigate the dynamics of the continental water cycle (Aemisegger et al., 2014; Huang and Wen, 2014; Delattre et al., 2015; Parkes et al., 2017) and can therefore be used to trace dew formation and radiation fog deposition into ecosystems (Spiegel et al., 2012; Wen et al., 2012; Delattre et al., 2015; He and Richards, 2015; Parkes et al., 2017). The isotopic composition of a water sample is expressed in terms of the abundance of hydrogen (2H and 1H) or oxygen (18O and 16O) isotopes by using the delta notation (hereafter δ) as ‰, where Rsample and Rstandard are the molar ratios of either 2H 1H or 18O 16O for the sample and standard, respectively. The standard is the Vienna Standard Mean Ocean Water (V-SMOW) controlled and distributed by the International Atomic Energy Agency (IAEA, 2009). With this definition, δ18O and δ2H are expressed as per mil (‰) discriminations from the standard. Water molecules with different isotopes are termed “isotopologues”. Three isotopologues, i.e., 1HO, 1HO, and 1H2H16O, are the most abundant in the water cycle. During phase changes such as evaporation and condensation, heavier isotopologues (i.e., 1HO and 1H2H16O) become enriched in the liquid phase and depleted in the gaseous phase, which thus causes an increase of δ2H and δ18O in the liquid phase and a decrease of δ2H and δ18O in the gaseous phase. During the evaporation and condensation processes, equilibrium fractionation always occurs at the interface between two phases and results in a ratio of 1:8 between the variability of δ18O and δ2H. When the ambient air is unsaturated, a deviation from the 1:8 ratio becomes measurable due to non-equilibrium fractionation (Dansgaard, 1964) driven by faster molecular diffusivity of the lighter isotopologue (i.e., 1HO) than its heavier counterparts (i.e., 1HO and 1H2H16O). The second-order parameter deuterium excess (hereafter d), defined as according to Dansgaard (1964), is a useful measure of non-equilibrium fractionation and provides information complementary to δ2H and δ18O. The d is often used as a tracer for the water vapor source of a given water pool in the water cycle (Gat, 1996; Yakir and Sternberg, 2000; Yepez et al., 2003; Welp et al., 2012; Aemisegger et al., 2014; Galewsky et al., 2016). For example, at the local scale, as compared to the higher d vapor of entrained free tropospheric air, local evapotranspiration is a vapor source with lower d because soil water vapor at the evaporation front has a lower d value (Delattre et al., 2015; Parkes et al., 2017). The diurnal cycle of deuterium excess in a well-mixed convective boundary layer has been studied previously (e.g., Lai and Ehleringer 2011), whereas relevant processes affecting d in the NBL are much less well known, in particular over grasslands.
2.2 Excluding the confusion of guttation
Long (1955) pointed out that guttation droplets distributed on the edges of plant leaves and can easily be mistaken by observers for dew droplets. Dew is however distinct from guttation, which is the exudation of drops of liquid from the hydathodes of the leaves of grasses driven by root pressure (Long, 1955; Stocking, 1956; Hughes and Brimblecombe, 1994). Both dew and guttation occur under high relative humidity. A soil water content near field capacity is favorable for guttation, whilst dew can also occur at very low soil water contents. In our study, we exclusively focused on the role of NRW during warm-season dry spells when soil water content in the main rooting zone was rather low, closer to the wilting point than to the field capacity, and hence guttation can be neglected here. Furthermore, guttation could easily be distinguished from dew by analyzing the stable isotopes of the respective water component: guttation stems from plant-internal water, whilst dew is plant-external water condensed from ambient water vapor or distilled from vapor related to the soil water isotopic signals. Consequently, the isotopic composition of guttation droplets should vary by species in parallel with the plant-internal water because no isotopic fractionation is expected during the guttation process. In all our samples, however, the isotopic composition of dew water was not related to the plant species from the surfaces of which the water was collected, which allowed us to exclude guttation as a relevant process during dry-spell periods.
3.1 Eddy covariance and meteorological measurements
The Chamau site (47∘12′36.8′′ N, 8∘24′37.6′′ E) is an intensively managed temperate grassland (4–6 cuts per year) at 393 m a.s.l., located in a valley bottom in Switzerland. The EC and meteorological measurement station (Fig. A1 in Appendix A) have been operational since 2005. The EC measurement setup consisted of a 3-D sonic anemometer (Gill R3, Gill Instruments Ltd., Lymington, UK) and an open-path infrared gas analyzer (IRGA, Li-7500, Li-Cor, Lincoln, NE, USA). The center of the sonic anemometer axis was at 2.4 m a.g.l. (see Fuchs et al., 2018, for more details). The EC measurements at 20 Hz were processed to 30 min averages using EddyPro Version 7.0.6 (LI-COR, 2017) and following established community guidelines (Aubinet et al., 2012; see also Appendix B) for horizontal wind speed (hereafter u2m, in m s−1), atmospheric specific humidity (hereafter qa2m, in g kg−1), dew point temperature (hereafter Td, in ∘C), turbulent latent heat flux (hereafter LE, in W m−2), turbulent sensible heat flux (hereafter H, in W m−2), and net radiation flux (hereafter Rn, in W m−2); negative fluxes denote a downward flux, whilst positive values stand for upward fluxes. Evapotranspiration (ET, in mm h−1) was derived from LE (see Appendix B). Ground heat flux (hereafter G, in W m−2) was measured at 0.02 m depth with two heat flux plates (HFP01 heat flux sensor, Hukseflux, Delft, the Netherlands).
The meteorological instruments were installed at 2.0 m a.g.l. (see Fuchs et al., 2018, for more details). Measurements were taken every 10 s and then aggregated to 30 min averages for air temperature (hereafter Ta2m, in ∘C), relative humidity (hereafter RH, in %) (a shaded, sheltered, and ventilated HydroClip S3, Rotronic AG, Basserdorf, Switzerland), as well as long-wave outgoing and ingoing radiation (hereafter LWout and LWin, in W m−2; obtained from a ventilated four-way CNR1 radiometer, Kipp & Zonen B.V., Delft, Netherlands, that also provided all-wave net radiation, Rn). The horizontal visibility (in km) was measured every 10 s with a fog sensor (MiniOFS, Optical Sensors Inc., Göteborg, Sweden) and a present weather detector (PWD10, Vaisala Oyj, Helsinki, Finland). The meteorological measurements were processed to 30 min averages for Ta2m, RH, and LWout and to 1 min averages for visibility. The vegetation surface temperature (T0, in ∘C) was determined following Stefan–Boltzmann's law as (Moene and van Dam, 2014)
where an emissivity (hereafter ε) of 0.98 was used to calculate temperatures for wet leaf surfaces (hereafter index w; T0=T0w), and a value of 0.96 was used for dry leaf surfaces (hereafter index d; T0=T0d) according to López et al. (2012); σ is the Stefan–Boltzmann constant at W m−2 K−1. The LWsurface was derived as suggested by Moene and van Dam (2014) as the difference between measured upwelling long-wave radiation LWout corrected for the first-order reflection of downwelling long-wave radiation LWin; i.e., LW.
The saturation specific humidity (q0, in g kg−1) and the relative humidity (h0) with respect to surface temperature T0 for wet and dry vegetation surfaces was calculated following Tetens' formula (Buck, 1981; Campbell and Norman, 1998) (see the equations in Appendix B).
Flux measurements were also used to assess the local surface energy budget as
where ΔQ is the energy budget closure term remaining when all other components (Rn, H, LE, G; in W m−2) are measured. A deviation of ΔQ from 0 W m−2 is typically a result of inaccuracies in determining the components of the energy budget, differences in footprint areas covered by the three different types of measurements (Rn: radiation flux; H and LE: turbulent fluxes; G: molecular flux), or advection of sensible and latent heat. Here we make the assumption that inaccuracies of the individual measurements do not change substantially over each field campaign, and variations of footprint areas mostly relate to H and LE, with smaller footprints during daytime and larger ones at night, whereas advective influences should be best detectable on the hourly timescale during the day–night transition around sunrise and sunset.
3.2 Experiment setup during the three 24 h observation campaigns
Three 24 h observation campaigns were carried out during expected dew/fog events on 25–26 July (event 1), 20–21 August (event 2), and 9–10 September (event 3) 2018. The time series were all recorded in CET (UTC+1). The precipitation at the Chamau site was 870 mm in 2018, which was 297 mm (about 25 %) less than the multiyear average over 2006–2017. The year-to-date precipitation before the three events was 393, 474, and 536 mm, respectively, which was 311 mm (−44 %), 359 mm (−43 %), and 367 mm (−41 %) less than the corresponding 2006–2017 averages (Fig. 1a). From April to September 2018, the average temperature was 17.3 ∘C, which was 1.8 ∘C higher than the corresponding 2006–2017 average (Fig. 1b). The corresponding consecutive rainless periods were 23–27 July, 18–21 August, and 8–12 September 2018 respectively. The daily average ET during the rainless periods was 2.7 mm (Fig. 1c).
Because of the extreme summer drought in 2018, no harvesting of the grassland was carried out during the three campaigns, but harvests were possible 46 d before event 1 on 9 June 2018 and 1 d after event 3 on 10 September 2018. The leaf area index was 1.5–2.5 m2 m−2 as measured 7 d before events 1 and 2 with a LAI-2000 plant canopy analyzer (LI-COR Biosciences, Lincoln, NE, USA). The mean vegetation height (zc) was roughly 0.2–0.3 m during the three campaigns. The wilting point, field capacity, and saturation water content (all in volumetric soil water content) were 12 %–14 %, 27 %–30 %, and 47 %–49 %, respectively, according to the soil texture reported by Roth (2006) and the equations by Saxton et al. (1986) (see details in Appendix C). The volumetric soil water content (SWC) was measured at 10, 20, 30, and 50 cm respectively (ML2x sensors, Delta-T Devices Ltd., Cambridge, UK). The plant roots were mainly distributed in the top 0–15 cm of the soil (Prechsl et al., 2015), and SWC in this layer was 17 %–20 % during the three events (Fig. 1d). The rainfall after event 1 was not sufficient to refill the deficient soil water storage, which explains why the observed SWC remained low until event 3.
3.2.1 Isotopic composition of non-rainfall water on foliage, leaf water, and soil water
To analyze the isotopic composition of NRW on foliage (hereafter fNRW), leaf water, and soil water, the sampling was carried out on a grassland area of 100×130 m2 around the EC and meteorological installations (Fig. A1 in Appendix A). NRW droplets on foliage (fNRW) were absorbed in triplicates with cotton balls from the leaf surfaces of randomly selected plants for Lolium sp. with long and narrow leaves; taller vegetation Taraxacum sp. with long and wide leaves; and shorter vegetation Trifolium spp. with short and wide leaves, as well as both shorter and taller vegetation. The fNRW samples were taken at the end of the nights of events 1 and 3 (one sampling per event) but every 2 h during the night of event 2 (i.e., four samplings in event 2). Simultaneously, leaf samples were taken in triplicates from the randomly selected plants for the three species after softly drying the leaf surfaces with tissue paper. To prevent the disturbance of destructive sampling on the effect of dew and fog formation, the NRW droplets and leaf samples were taken from different plants of the same species in the sampling area. The soil cores were taken with a soil auger and were then cut into slabs to separate the soil depths of 0–5, 5–10, 10–15, 15–20, and 20–40 cm. Soil samples in event 1 were taken without replicate within 2 h before sunset and at the end of the night; soil samples in event 2 was taken without replicate within 2 h before sunset, as well as every 2 h (i.e., four samplings in event 2) during the night; soil samples in event 3 were taken in triplicates within 2 h before sunset and at the end of the night.
After collection, NRW droplets on foliage (fNRW) and leaf and soil samples were immediately transferred into gas-tight 12 mL exetainers (Labco Exetainer® vials, High Wycombe, UK) and stored in a portable cooling box filled with ice blocks. Before extracting the water in a cryogenic vacuum extraction system (Prechsl et al., 2015), the samples were stored at −19∘ C. The isotopic composition of extracted water samples for fNRW (hereafter δ18OfNRW and δ2HfNRW), leaf water (δ18Oleaf, and δ2Hleaf), and soil water (hereafter δ18Os, and δ2Hs) was measured using an isotope ratio mass spectrometer (IRMS, DELTAplusXP, Finnigan MAT, Bremen, Germany). The measured uncertainties of δ18O and δ2H using IRMS are ±0.1 ‰ and better than ±1.0 ‰, respectively (Werner and Brand, 2001; Gehre et al., 2004).
3.2.2 Isotopic composition of ambient water vapor and non-rainfall water condensed from this vapor
The isotopic composition and the volumetric mixing ratio of ambient water vapor were measured at 0.5–1 Hz using a cavity ring-down laser absorption spectrometer (L2130-i, Picarro Inc., Santa Clara, CA, USA). The L2130-i was placed in a house 200 m away from the EC and meteorological measurements (Fig. A1 in Appendix A). Ambient air was pulled into the instrument through a PTFE intake hose, with an outer diameter of in., and a PTFE-filter inlet (FS-15-100 and TF50, Solberg International Ltd., Itasca, IL, USA) fixed at 6 m a.g.l. The intake hose was thermally isolated and heated using a resistive heating wire (Raychem 5BTV2-CT, Von Rotz, Kerns, Switzerland) that was wrapped around the entire length of the intake tube to prevent condensation and minimize the response time of the inlet system. An external membrane pump (N022, KNF Neuberger GmbH, Munzingen, Freiburg, Germany) with a flow rate of 9 L min−1 was used to maintain turbulent flow (Reynolds number Re >2900) in the tube to minimize memory effects within the inlet system. The isotopic composition of ambient water vapor (hereafter δa) and the volumetric ambient water vapor mixing ratio (hereafter wa) were measured using a flow split with a flow rate of 300 mL min−1 through the L2130-i cavity. The instrument's response time in this setup was found to be on the order of 10 s in Aemisegger et al. (2012).
To correct for instrument drifts and to normalize the data to the international VSMOW–SLAP (Vienna Standard Mean Ocean Water–Standard Light Antarctic Precipitation) scale, the raw data were calibrated using a standard delivery module (SDM; A0101, Picarro Inc., Santa Clara, CA, USA) by performing two-point calibrations every 12 h (Aemisegger et al., 2012) using two liquid standards (standard 1: δ18O = −11.43 ‰, δ2H = −81.84 ‰, d= 9.64 ‰; standard 2: δ18O = −40.66 ‰, δ2H = −325.67 ‰, ‰ measured by an IRMS). The δ18O and δ2H of the standards thus bracket the range of the measured δ18Oa and δ2Ha. Laser spectrometric measurements are known to be affected by a bias dependent on the water vapor mixing ratio due to spectroscopic effects (absorption peak fitting and baseline effects). In our study, all measurements were performed at wa>12 mmol mol−1; therefore no mixing-ratio-dependent isotope bias correction was necessary (see more details in Aemisegger et al., 2012). The L2130-i was calibrated using a dew point generator (LI-610, Li-Cor Inc., Lincoln, NE, USA) following the procedure by Thurnherr et al. (2020). Calibrated δ18Oa and δ2Ha were then averaged over 30 min intervals. The second-order parameter d of ambient water vapor (hereafter da) was calculated with the calibrated δ18Oa and δ2Ha. The overall random uncertainties of δ18O and δ2H measurements were 0.2 ‰ and 0.8 ‰ respectively (for more details about the uncertainty quantification, see Aemisegger et al., 2012).
To analyze the correlation between da and surface humidity, the surface relative humidity (RH0 in %) computed from water vapor mixing ratio wa and surface saturation specific humidity (q0w; see Eq. B1 in Appendix B) was calculated as
where Mv=0.018015 kg mol−1 is mole weight for water vapor, Md=0.028965 kg mol−1 is mole weight for dry air, a is a unit conversion factor (10−3 mol mmol−1 ×103g kg−1), and qa_L2130i is specific humidity (in g kg−1) computed from wa.
Ambient water vapor is one source of NRW on foliage (fNRW) which experiences fractionation during the condensation process. With the assumption of equilibrium fractionation, the isotopic composition of equilibrium liquid (hereafter aNRW, and its isotopic composition δ18OaNRW and δ2HaNRW) formed from ambient water vapor δ18Oa and δ2Ha was calculated using the temperature-dependent equilibrium fractionation factors following Horita and Wesolowski (1994) as
where and were equilibrium fractionation factors calculated as (Horita and Wesolowski, 1994)
An approach to calculate the NRW isotope composition from ambient vapor, which considers both equilibrium and non-equilibrium fractionation in the laminar sublayer of the leaf boundary layer, has been proposed by Wen et al. (2012). The isotope composition of the NRW formed from ambient vapor under such conditions (hereafter naNRW, and its isotopic composition δ18OnaNRW and δ2HnaNRW) was calculated as follows:
where δnaNRW is either δ18OnaNRW or δ2HnaNRW, ϵk is the non-equilibrium fractionation factor in permil, calculated from ‰, given , following Merlivat (1978), and m=0.67 for laminar flow following Dongmann et al. (1974); ϵeq is equilibrium fractionation factor in permil calculated from ( ‰ in Eqs. (6) and (7).
3.2.3 Determination of the atmospheric layer heights and assessment of eddy covariance setup height
The isotopic fractionation during phase change at the Earth surface is linked to the micrometeorological layers near the surface (Fig. 2). The inclusion of a zero-plane displacement (or fluid dynamic height origin, zd) (Fig. 2) in wind profiles allows us to separate the downward flux from ambient water vapor and the upward flux from soil-diffusing vapor. The average wind speed is zero at zd+z0, where z0 is aerodynamic roughness length (z0). The roughness length z0 at the Chamau site was 0.03 m on average. It was computed by solving the logarithmic wind profile equation for z0 using measured horizontal wind speed u2m and friction velocity u*,
during neutral atmospheric stratification (e.g., Panofsky, 1984; see data in Appendix D), with z2m the measurement height (2 m) and κ the von Kármán constant (0.40). The zero-plane displacement zd can be approximated as two-thirds of vegetation height (Stull, 1988; Oke, 2002), i.e., 0.13–0.20 m. With respect to –0.23 m, we consider three pathways of NRW inputs onto the foliage of grasslands for dew and radiation fog: (1) the downward component of dew formation condensing from ambient water vapor, (2) the upward component of dew formation via distillation of water vapor from soil, and (3) radiation fog deposition.
The top of NBL is difficult to quantify because in many cases the NBL does not have a strong demarcation at its top. Therefore, many definitions of the NBL are based on relative comparisons of the stable boundary layer state aloft to near-surface state (Stull, 1988). We determined the top of the NBL as the lowest height where the vertical stratification of the atmosphere becomes isothermal, i.e., ∂T/ 0 (Stull, 1988; Garratt, 1992), where T is air temperature extracted from the hourly COSMO-1 model (Consortium for Small-scale Modeling) with a resolution of 1.1 km (meridional) ×1.1 km (zonal) over Switzerland (MeteoSwiss, 2021; Westerhuis et al., 2020) and 80 vertical levels. During the three events in this study, the NBL top was at 114, 55, and 193 m a.g.l., respectively (Fig. 3). Therefore, the EC measurement setup at 2.4 m a.g.l. is expected to have captured roughly 98 % of the expected flux (Eugster and Merbold, 2015). The roughness sublayer (1–3 times the vegetation height according to Oke, 2002) was at 0.2–0.9 m at the Chamau site; therefore the EC instruments were installed well above the roughness sublayer. Here we simply use NBL as background information on atmospheric stability but did not use it for NBL budgets (Denmead et al., 1996), as was done by Stieger et al. (2015) at this exact same site, and thus the uncertainty in the exact value extracted for the NBL top from the COSMO-1 model output has no influence on our dew estimates.
3.2.4 Partitioning of non-rainfall water inputs using a two-end-member mixing model
We partitioned the contribution of NRW input pathways into the two main processes: (1) the downward component of dew formation and fog droplet deposition (aNRW) and (2) the distillation of soil-diffusing vapor on plant leaves. With unsaturated conditions, NRW on foliage (fNRW) was a mix of aNRW and distillation, while with saturated conditions, fNRW originated from dew or from fog (aNRW), which could lead to a mixture of water from both sources over the course of a night when dew and fog occur intermittently. “Unsaturated conditions” in this context refers to the standard meteorological measurements at 2 m a.g.l. level. Dew forming in unsaturated conditions is a mixture of aNRW and distillation but lacks contribution from fog deposition. Thus, the isotopic signature of NRW resulting from the isotopic composition of distillation (hereafter δ18Odistillation and δ2Hdistillation) and the proportion of distillation (hereafter fdistillation) in fNRW can be expressed as
where faNRW is the proportion of aNRW in fNRW. The four parameters δ18Odistillation, δ2Hdistillation, fdistillation, and faNRW are unknown. Therefore, solving for four unknowns with only three equations (Eqs. 10–12) requires two time points of measurements (here we used 23:00 and 01:00 CET in event 2) to obtain empirical estimates for the four unknowns. By doing so, we implicitly assumed that δ18Odistillation and δ2Hdistillation were constant over time (i.e., did not change within this 2 h interval during event 2), and only fdistillation and faNRW were allowed to change between these two sampling times. For δfNRW, the median value for each sampling was taken, and for δaNRW, the period between two measurements was computed from 30 min data. Consequently, the three equations (Eqs. 10–12) can be expanded to six equations via the inclusion of two sampling times (τ and τ+1) as
which can be solved for the six unknowns δ18Odistillation, δ2Hdistillation, fdistillation_τ, fdistillation_τ+1, faNRW_τ, and faNRW_τ+1 using the “limSolve::Solve” function in R (Venables and Ripley, 2002).
3.2.5 Partitioning of non-rainfall water inputs using Monteith (1957) approach (M57)
To assess the results from our mixing model by Eqs. (13)–(18), the partitioning of NRW components was also performed using the Monteith (1957) approach (hereafter M57), i.e., partitioning the NRW components from the amount of aNRW and distillation. The amount of NRW from soil diffusing vapor was calculated as follows based on the near-ground vertical temperature gradient:
where Kv is the diffusion coefficient given m2 s−1 (Monteith, 1957); Ts1cm in ∘C is the soil temperature measured at 1 cm in depth; Ta1cm in ∘C is the air temperature at 1 cm a.g.l., which was computed from the simulated wet vegetation surface temperature T0w and measured soil temperature Ts1cm:
and the saturated absolute humidity χs1cm and χa1cm at soil temperature at 1 cm depth (Ts1cm) and air temperature at 1 cm (Ta1cm) were calculated following Parish and Putnam (1977) as
where T is substituted by either Ts1cm or Ta1cm to calculate χs1cm and χa1cm with relative humidity h=100 %, respectively.
The condensation rate of ambient water vapor was calculated as (Pasquill, 1949; Monteith, 1957)
where ∂χ/∂z is the gradient of absolute humidity from Ta at z=2 m and from T0w at ; thus )]; Φ is the stability parameter proportional to the Richardson number Ri with numerous semi-empirical forms (Garratt, 1992); here we followed Monteith (1957) given ), in which σ is a proportionality factor associated with thermal stratification assumed to be on the order of 10 (Pasquill, 1949; Monteith, 1957); Ri is calculated as (Wyngaard, 2010)
where L in meters is the Monin–Obukhov length calculated following Monin and Obukhov (1954); other semi-empirical forms of Φ and its effect on NRW amount estimates are given in Appendix E.
3.3 Statistics and imaging
We report means ± SD (standard deviation), unless specified differently. For the isotopic composition of NRW on foliage (δ18OfNRW, δ2HfNRW, and dfNRW) and leaf water (δ18Oleaf, δ2Hleaf, and dleaf), we report the interquartile range (25 % and 75 % quantile) together with the median to account for the unknown empirical distribution of destructive sampling of individual plants. The statistical significance of among-species differences was assessed with Tukey's honest significant difference test using the “agricolae:: HSD.test” function in R. All analyses were performed with R version 4.0.4 (R Core Team, 2020). Orthogonal regression was performed using the “mcr::mcreg” function in R (total least-squares fit; Gat, 1981) for all linear regression analyses.
4.1 Environmental conditions during dew and radiation fog events
Dew and radiation fog generally form during clear-sky nights with low wind speed and weak turbulence. During the three field campaigns presented in this study, wind speed (u2m) and latent heat flux (LE) showed an abrupt weakening from around 17:00 CET onwards (Fig. 4a, b). With nightfall, u2m remained below 0.7 m s−1 (Fig. 4a), and LE was very low (−26 to 14 W m−2; Fig. 4b), indicating a vanishing of turbulent fluxes. These are favorable conditions for dew and radiation fog formation.
The three events with dew or radiation fog were characterized by high relative humidity (RH) with respect to air temperature measured at 2 m a.g.l. From around 17:00 CET, RH increased rapidly and reached 100 % around 03:00 CET during event 2 and around 20:30 CET during event 3 (Fig. 4c). These saturated conditions led to the formation of fog characterized by a horizontal visibility <1 km (Fig. 4d). Fog appeared around 05:00 CET during event 2, lasting for less than an hour until sunrise, whilst the onset of fog was much earlier during event 3 (around 23:00 CET), lasting for a longer period until dissipation around sunrise. The visibility was always >1 km in event 1, indicating that fog was absent during event 1. Therefore, event 1 can be considered to be a dew-only event, whilst events 2 and 3 were characterized by a combination of dew and the partial influence of radiation fog.
Dew or radiation fog occurred when the surface cooled below dew point. Both grassland surfaces and ambient air started to cool from around 17:00 CET onwards, due to substantial net long-wave radiation loss (−36 W m−2 at sunset; Fig. 5a). The vegetation surfaces of the grassland cooled more rapidly than the near-surface atmosphere; thus with nightfall, the vegetation surface temperature T0 derived from radiation measurement remained cooler than air temperature Ta2m at 2 m a.g.l., although both gradually decreased (Fig. 5b). The first sign of condensation occurred when the leaf surfaces cooled below dew point temperature (i.e., T0<Td; Fig. 5b). The level of computed dry surface temperature T0d became lower than dew point Td at around 00:30 CET in event 1, 21:30 CET in event 2, and 19:00 CET in event 3 (Fig. 5b), determining the time when the first signs of condensation can be expected. During event 3, the surface already cooled below the dew point rapidly after sunset (i.e., T0<Td in Fig. 5b), indicating that condensation already started with nightfall.
Dew and radiation fog were characterized by a decrease in specific humidity (Fig. 5c). But before the formation of dew and fog set in, the specific humidity of the air (qa2m) steeply increased by 2.0–3.5 g kg−1 from around 17:00 CET until sunset (Fig. 5c), suggesting the mixing of moisture from local evapotranspiration into a shallow inversion layer. With nightfall, qa2m reached a nighttime maximum of 9.6–12.5 g kg−1 (Fig. 5c). In particular, in events 1 and 2, before starting to decrease, qa2m fluctuated for a short period from sunset until the first sign of condensation (Fig. 5c). When condensation started (T0<Td, Fig. 5b), qa2m gradually decreased (Fig. 5c). With the saturation specific humidity at surface temperature (q0) falling to values below qa2m (Fig. 5c), computed theoretical surface relative humidity h0 exceeded 100 % (Fig. 4c). The decrease of qa2m was much faster in event 3 (0.4 g kg−1 h−1; Fig. 5c) than that in events 1 and 2 (0.2 and 0.3 g kg−1 h−1; Fig. 5c), indicating stronger condensation of ambient water vapor.
According to the variability of environmental conditions, water vapor and thermal dynamics of dew and radiation fog events can be separated into four periods from 17:00 CET until sunrise: (1) pre-condensation period (hereafter P1 period) with the gradual weakening of turbulence, and warmer surface above the dew point (T0>Td; Fig. 5b); and (2) condensation period (hereafter P2 period) with a cooler surface below dew point (T0<Td; Fig. 5b). The P1 period was further separated into the P1a period starting around 17:00 CET until sunset, with the weakening of turbulence and the increase of specific humidity qa2m, and the P1b period from sunset until the first sign of condensation, with short-term fluctuations of specific humidity (qa2m). The P2b period was further split into the P2a period with dew only in the conditions of RH <100 % and the P2b period with combined dew and radiation fog in the conditions of RH = 100 %.
4.2 Isotopic dynamics of ambient water vapor during dew and fog events
The four periods of water vapor and thermal dynamics defined in Sect. 4.1 are reflected in the temporal evolution of the volumetric water vapor mixing ratio (wa) and the isotopic composition of ambient water vapor (δ18Oa, δ2Ha, and da; Fig. 6). From 17:00 CET until sunset (P1a period), when the turbulence was weakening and the surface was cooling, wa, δ18Oa, and δ2Ha showed a steep increase by 0.3–0.4 mmol mol−1, 2.0 ‰–3.2 ‰, and 7.4 ‰–12.5 ‰, respectively (Fig. 6a, b, c), whilst da showed a steep decrease by 11.6 ‰–16.9 ‰ (Fig. 6d). The decrease in da and increase in δ18Oa and δ2Ha were due to the effect of local evapotranspiration under the conditions of reduced entrainment from the free troposphere. The vapor sourced from local evapotranspiration features a lower da and higher δ18Oa and δ2Ha than the free troposphere. With nightfall, wa, δ18Oa, and δ2Ha reached a plateau, with 15.5 to 17.8 mmol mol−1 in wa (Fig. 6a), −15.5 ‰ to −14.3 ‰ in δ18Oa (Fig. 6b), and −128.0 ‰ to −113.2 ‰ in δ2Ha (Fig. 6c).
The start of condensation then caused a decrease of δ18Oa and δ2Ha because heavier water isotopologues have a lower partial vapor pressure at saturation than their lighter counterpart (psat[1H2H16O] < psat[1HO] <psat[1HO]) and thus preferentially prevail in the phase with stronger bonds (liquid > vapor; Bigeleisen, 1961). During the condensation period with RH <100 % (P2a period), wa steeply decreased by 0.8–5.5 mmol mol−1 (Fig. 6a), δ2Ha decreased by 3.3 ‰–16.7 ‰ (Fig. 6c), and da reached its minimum at −11.8 ‰ to −4.7 ‰ (Fig. 6d). During the condensation period with RH =100 % (P2b period), the decreasing rate of δ2Ha in event 3 (1.6 ‰ δ2Ha h−1) was almost double compared to that in events 1 and 2 (0.8 and 1.0 δ2Ha h−1 respectively, Fig. 6c), suggesting stronger condensation in event 3. Note that the changes of δ18Oa and da (Fig. 6b, d) depended on the humidity dynamics and the occurrence of dew and fog (Fig. 4c, d). During the dew-only periods (P2a period) in events 1 and 2 (Fig. 4d), δ2Ha decreased by 3.3 ‰–5.7 ‰ (Fig. 6c), and da slightly decreased by 3.4 ‰–3.7 ‰ (Fig. 6d), while δ18Oa showed fluctuations around the maximum reached 4 h and 2 h after nightfall of events 1 and 2 respectively (−15.5 ‰ to −14.3 ‰; Fig. 6b). The slight fluctuation of δ18Oa and decrease of da during the P2a period were a result of concurrent evaporation, which leads to an additional non-equilibrium fractionation with variations of δ18Oa:δ2Ha deviating from 1:8. Furthermore, as condensation was stronger than evaporation (i.e., net condensation), this caused a decrease of wa and δ2Ha (Fig. 6a, c). Because δ18Oa is more sensitive to evaporation than δ2Ha due to the higher partial vapor pressure of 1HO than 1H2H16O, evaporation accompanying condensation is the likely reason for the fluctuations of δ18Oa (Fig. 6b) but had only a minor effect on the variability of δ2Ha (Fig. 6c). During P2b periods in events 2 and 3 with dew and fog in combination, both δ18Oa and δ2Ha gradually decreased (by 0.3 ‰–1.5 ‰, and 2.1 ‰–12.8 ‰ respectively) with a ratio of around 1:8 (Fig. 6b, c); hence da was relatively constant during the nighttime minimum (−6.0 ‰ to −4.7 ‰; Fig. 6d), with only small fluctuations. In this saturated condition, evaporation was negligible, and condensation was the dominant process. This is confirmed by the constant values of da during P2b (Fig. 6d), showing that this period was dominated by equilibrium fractionation.
4.3 Isotopic composition of different non-rainfall water components
As one of the components of NRW on foliage (fNRW), the isotopic composition of NRW equilibrium liquid from ambient water vapor (aNRW) was comparable with the isotopic composition of fNRW. The isotopic composition of aNRW was −5.0 ‰ to −4.3 ‰ for δ18OaNRW, −47.4 ‰ to −38.6 ‰ for δ2HaNRW, and −12.1 ‰ to −2.4 ‰ for daNRW (Fig. 7a, b, c). For comparison, NRW on foliage (fNRW) was −6.1 ‰ to −1.5 ‰ for δ18OfNRW, −64.3 ‰ to −35.6 ‰ for δ2HfNRW, and −33.8 ‰ to 8.0 ‰ for dfNRW (Fig. 7a, b, c). The isotopic composition of fNRW varied over time with gradually decreasing δ18OfNRW (Fig. 7a) but gradually increasing δ2HfNRW (Fig. 7b) and dfNRW (Fig. 7c). The relationship between the isotopic composition of fNRW and aNRW was related to humidity conditions. With unsaturated conditions when dew formation occurred, δ18OaNRW ( ‰; Fig. 7a) was lower than δ18OfNRW (−3.8 ‰; Fig. 7a), while δ2HaNRW ( ‰; Fig. 7b) was higher than δ2HfNRW (−47.7 ‰; Fig. 7b), and daNRW ( ‰; Fig. 7c) was higher than dfNRW (−20.5 ‰; Fig. 7c). With saturated conditions at 03:00 and 05:00 CET of event 2, the isotopic composition of aNRW ( ‰ in δ18OaNRW, ‰ for δ2HaNRW, and ‰ for daNRW; Fig. 7) was identical to the isotopic composition of fNRW (−4.7 ‰ for δ18OfNRW, −43.0 ‰ for δ2HfNRW, and −5.4 ‰ for dfNRW; Fig. 7). In particular, with saturated condition at 05:00 CET in event 3 when radiation fog occurred, δ18OfNRW and δ2HfNRW were lowered by 0.7 ‰ and 1.4 ‰ with respect to δ18OaNRW and δ2HaNRW respectively (Fig. 7a, b), and daNRW was 5.5 ‰ higher than dfNRW (Fig. 7c).
The isotopic composition of the distillation component, i.e., NRW from soil-diffusing vapor, was computed with a two-end-member mixing model using the values from 23:00 to 01:00 CET in event 2. In unsaturated conditions, with respect to aNRW, δ18OfNRW and δ2HfNRW deviated to the higher and lower sides of δ18OaNRW and δ2HaNRW, respectively (Fig. 7a, b). This is in contrast to the effect that evaporation would have had; both δ18OfNRW and δ2HfNRW would be higher than δ18OaNRW and δ2HaNRW. Therefore, we assumed that the observed deviations of δ18OfNRW and δ2HfNRW with respect to δ18OaNRW and δ2HaNRW were caused by the mixed source of NRW on foliage, i.e., the mixing of NRW from ambient water vapor and soil-diffusing vapor (i.e., distillation). Based on the measurements from 23:00 to 01:00 CET in event 2, the averages of δ18Odistillation, δ2Hdistillation, and ddistillation during this 2 h period were computed as −1.0 ‰, −71.8 ‰, and −63.4 ‰ respectively (Fig. 7) via the mixing model. As a comparison, from 1 h before sunset till sunrise in event 2, the isotopic composition of soil water in 0–40 cm varied in the range of −10.4 ‰ to 5.5 ‰ for δ18Os, −78.8 ‰ to −8.5 ‰ for δ2Hs, and −52.4 ‰ to 4.1 ‰ for ds (Fig. 8). The computed distillation δ18Odistillation and δ2Hdistillation fell in the range of the soil water δ18Os and δ2Hs (Fig. 8a, b), whilst ddistillation was lower than the soil water ds (Fig. 8c) probably derived from the uncertainty of δ18Odistillation and δ2Hdistillation estimates (see in Sect. 5.3).
The relationships of δ2HfNRW–δ18OfNRW and δ2HaNRW–δ18OaNRW with respect to the local meteoric water line (LMWL: δ2H = 7.68×δ18O + 6.97; Prechsl et al., 2014) suggested that the local vapor is the primary source for dew and radiation fog during all three events (Fig. 9). Both δ2HfNRW–δ18OfNRW and δ2HaNRW–δ18OaNRW fell to the right-hand sides of the LMWL, suggesting lower d from NRW inputs as compared to local precipitation. When we only considered the condensation of ambient water vapor under equilibrium fractionation, δ2HfNRW and δ18OfNRW pairs fell on the δ2HaNRW–δ18OaNRW regression line (for the sampling at 03:00 and 05:00 CET in event 2; Fig. 9). However, with the mix of the component condensing from soil-diffusing vapor (distillation) in the conditions of RH <100 %, the δ2HfNRW– δ18OfNRW pairs fell to the right-hand sides of the δ2HaNRW–δ18OaNRW regression line (for the sampling at 03:00 CET in event 1 and the samplings at 23:00 and 01:00 in event 2; Fig. 9), and dfNRW was lower than daNRW ( ‰ for dfNRW, and ‰ for daNRW; Fig. 7c). This suggested that the soil-diffusing vapor was a lower d vapor source as compared to the ambient water vapor daNRW (Fig. 9), which corresponded to the fact that the soil water ds ( ‰, Fig. 8c) was lower than daNRW, whereas with the mix of the component from radiation fog deposition in the conditions of RH = 100 %, δ2HfNRW–δ18OfNRW pairs fell to the left-hand sides of the δ2HaNRW–δ18OaNRW regression line (for the sampling at 05:00 CET in event 3); hence the corresponding dfNRW was higher than daNRW (Fig. 9).
The condensation of ambient water vapor for dew formation can be approximated as an equilibrium fractionation process (e.g., Wen et al., 2012, and Delattre et al., 2015); the condensation of ambient water vapor to form radiation fog can cause lower δ18O and δ2H of NRW on foliage compared to NRW equilibrium liquid obtained from ambient water vapor. When considering non-equilibrium fractionation, the isotopic composition of NRW from ambient water vapor (δ18OnaNRW and δ2HnaNRW; Fig. 7a, b) was much lower than the isotopic composition of NRW on foliage (δ18OfNRW and δ2HfNRW; Fig. 7a, b), and the lowering of δ18OnaNRW and δ2HnaNRW was more severe with the increase of the computed relative humidity (h0; Fig. 4c) at surface temperature. The lowering of δ18OnaNRW and δ2HnaNRW with respect to δ18OfNRW and δ2HfNRW was most likely due to the overestimate of the non-equilibrium fractionation factor when computed h0 exceeded 100 % (going up to 132 %; Fig. 4c). Non-equilibrium fractionation is usually considered to be negligible above −10 ∘C in the process of vapor condensing to liquid in clouds (Jouzel et al., 1987). However, non-equilibrium fractionation driven by molecular diffusion might have played an important role in a laminar fog boundary layer (FBL) (Castillo and Rosner, 1989; Epstein et al., 1992), which led to lower δ18OfNRW and δ2HfNRW than δ18OaNRW and δ2HaNRW at 05:00 CET in event 3 (Fig. 7a, b) when radiation fog occurred (Fig. 4d). Heavier isotopologues move more slowly than their lighter counterpart in air (molecular diffusivity: D[1HO] < D[1H2H16O] < D[1HO]; Merlivat, 1978); hence the rate at which heavy isotopologues (1HO and 1H2H16O) in ambient air pass through the laminar FBL to be condensed at the liquid–vapor interface is smaller than the rate of condensation of their lighter counterpart (1HO). Therefore, δ18OfNRW and δ2HfNRW can become lower than δ18OaNRW and δ2HaNRW. Fog lasted from 23:00 CET until sunrise of event 3 and appeared around 5:00 CET within half an hour before sunrise in event 2 (Fig. 4d). However, we only observed a lower δ18OfNRW and δ2HfNRW than δ18OaNRW and δ2HaNRW in event 3 (Fig. 7a, b), suggesting that the lowering of δ18OfNRW and δ2HfNRW might also be related to the duration of radiation fog.
4.4 Contribution of distillation in the total non-rainfall water on foliage
The contribution of distillation in the total NRW on foliage (fNRW) was computed via the mixing model using the values from 23:00 to 01:00 CET in event 2 and the M57 approach, respectively (Fig. 10). Based on this model we estimated a contribution of 28 % and 9 % of foliage NRW (fNRW) from distillation at 23:00 CET and 01:00 CET of event 2, respectively (Fig. 10b; Table 1); hence 72 % and 91 %, respectively, was dew condensed from ambient water vapor. A linear extrapolation of fdistillation to the beginning of dew formation at 21:30 CET of event 2 increased the contribution of distillation to 42 % (Fig. 10b; Table 1), and thus the contribution of aNRW was 58 %. Similarly, when using the values of δ18Odistillation and δ2Hdistillation computed from event 2 for estimating fdistillation during event 1, the contribution of distillation was around 18 %–31 %, and thus the contribution of aNRW was around 69 %–82 % for our sampling at 03:00 CET of event 1 (vertical whiskers in Fig. 10b; Table 1). For comparison, the contribution of distillation was also calculated using the M57 approach (Eqs. 19 and 22; Fig. 10a). The dew and radiation fog potentially produced 0.17–0.54 mm d−1 NRW gain on foliage, which, compared to evapotranspiration water loss of on average 2.7 mm d−1, constitutes a non-negligible water flux into the canopy. The computed dew water gain from aNRW (0.12–0.50 mm) was generally larger than the internal redistribution via distillation (0.04–0.05 mm) (Fig. 10a). As the nights progressed, the contribution of distillation (fdistillation) to NRW on foliage (fNRW) decreased from 76 % at 00:30 CET to 27 % before dawn in event 1 and from 45 % at 21:30 CET to 19 % before dawn in event 2. Overall lower fdistillation was observed in events 2 and 3 as compared to that of event 1. No clear trend was observed for fdistillation in event 3, with slight variations around 6 %–8 % (Fig. 10b). The fdistillation estimate from the mixing model during events 1 and 2 agrees well with the M57 approach (compare black and grey data in Fig. 10b).
5.1 Diurnal patterns of isotopic composition in ambient water vapor
The diurnal patterns of d for ambient water vapor were mainly affected by the entrainment from the free troposphere and local evapotranspiration (Lee et al., 2006; Lai and Ehleringer, 2011; Welp et al., 2012; Huang and Wen, 2014; Delattre et al., 2015; Parkes et al., 2017). Moreover, the effect of local evapotranspiration might be enhanced by density-driven katabatic drainage flow down the slopes of the local topography (Drobinski et al., 2003; Whiteman et al., 2010; Nadeau et al., 2013; Duine et al., 2016) and by the regional thermodynamic conditions with weak large-scale influence during clear and calm nights (Eugster and Siegrist, 2000; Goulden et al., 2006; Eugster and Merbold, 2015). Entrainment from the free troposphere played a dominant role in midday atmospheric water vapor dynamics, whilst local evapotranspiration was the main driver of atmospheric water vapor dynamics in the late afternoon when entrainment from free troposphere was already reduced. Entrainment from the free troposphere is a vapor source with lower δ18O and δ2H and higher d, whilst local evapotranspiration is a vapor source with higher δ18O and δ2H and lower d (Parkes et al., 2017). Consequently, as compared to the nighttime periods, we observed a higher da (Fig. 6d) and a decrease in δ18Oa and δ2Ha (Fig. 6b, c) during 13:00–17:00 CET. Although evapotranspiration is stronger at midday as compared to late afternoon, evapotranspiration is not the main factor controlling δ18Oa and δ2Ha variabilities at midday. On the contrary, during the periods of turbulence weakening and surface cooling from around 17:00 CET to sunset with the reduced entrainment from the free troposphere (weakened u2m and reduced LE in Fig. 4a, b), local evapotranspiration became the main driver of isotopic dynamics of ambient water vapor for the three events in our study. This combination of weakening entrainment and evapotranspiration into a shallower mixed layer hence caused a steep decrease in da (Fig. 6d) and increases in δ18Oa and δ2Ha (Fig. 6b, c) during the P1a period, which is in accordance with previous studies by Huang and Wen (2014) and Parkes et al. (2017). The soil moisture at 0–5 cm in a short period before sunset (e.g., within 1 h before sunset) showed extremely varied isotopic composition from −8.5 ‰ to 5.9 ‰ for δ18Os_0–5cm, from −72.8 ‰ to −8.5 ‰ for δ2Hs_0–5cm, and from −52.4 ‰ to −4.1 ‰ for ds_0–5cm (Fig. 8), which is in accordance with the report by Welp et al. (2012) that soil evaporation showed very large variability of isotopic signals. The chamber experiment by Parkes et al. (2017) showed that the soil water vapor at the evaporation front had much higher δ18O and δ2H and much lower d as compared to the soil water at 0–5 cm. This much higher δ18O and δ2H and lower d vapor source at the soil evaporation front might have caused an enhanced variability observed in δ18Oa, δ2Ha, and da in ambient water vapor in the P1a period. As the Chamau site studied here is located in a valley bottom, the relative energy budget closure differed slightly from zero (Fig. 11e) in period P1a, suggesting that the effect of local evapotranspiration on the isotopic dynamics of ambient water vapor might have been slightly accompanied by cold-air drainage towards the valley bottom. Non-equilibrium fractionation is intrinsically dominant in the processes of evaporation with unsaturated ambient air (RH <100 % at 2 m a.g.l.), which induced a slight decrease of da during the condensation period P2a with RH <100 % (Fig. 6d). During the dew and radiation fog period P2b with RH =100 % at 2 m a.g.l., the condensation of ambient water vapor could essentially be described by an equilibrium fractionation process, with da remaining constant at a low nighttime minimum level (Fig. 6d), which is in accordance with the results by Huang and Wen (2014) and Delattre et al. (2015).
Isotopic signals in ambient water vapor provide information on the strength of continental moisture recycling (Aemisegger et al., 2014). In particular, the da has been shown to be a useful tracer for moisture source conditions and to be strongly anticorrelated with the surface relative humidity RH0 (computed from wa using Eq. 3) at the moisture source location (Craig and Gordon, 1965; Pfahl and Wernli, 2008; Welp et al., 2012; Aemisegger et al., 2014). The physical foundation for this strong link is the sensitivity of da to the non-equilibrium fractionation effect. The lower the surface relative humidity (RH0), the stronger the non-equilibrium fractionation and the higher da becomes. Spiegel et al. (2012) found an exceptionally high d in fog droplets after the passage of a cold front in central Europe with important moisture advected from the subpolar North Atlantic with anomalously low RH0. In Aemisegger et al. (2014), synoptic events were classified into events with a remote or local moisture source based on backward trajectories and a detailed correlation analysis between da and surface relative humidity. They found that events dominated by local sources show a strong anticorrelation between da and local surface relative humidity. In our study, da shows a strong anticorrelation with RH0 (; Fig. 12), suggesting that dew and radiation fog was formed from local moisture as a vapor source. The slope of the da–RH0 relation found here (−0.26 ‰ %−1) is similar to the relations found at another Swiss grassland site in dry summer periods (−0.17 ‰ %−1 by Aemisegger et al., 2014). From this analysis, we conclude that during the studied events, the isotopic signals were dominated by local moisture and that large-scale advection with the weak synoptic-scale flow in the context of central European anticyclones likely had a negligible influence.
5.2 Processes affecting non-rainfall water on foliage
Besides the main contribution of NRW from ambient water vapor to dew formation and radiation fog deposition, NRW on foliage (fNRW) can also be affected by three additional processes: (1) re-evaporation of NRW on foliage (He and Richards, 2015); (2) distillation (Monteith, 1957); and (3) guttation (Hughes and Brimblecombe, 1994; Xu et al., 2019). The role of distillation was quantified in Sect. 4.4, and in the following we argue why the other two additional processes at most had a minor influence on dew formation during all three events investigated here. Re-evaporation should have caused both δ18OfNRW and δ2HfNRW being higher than δ18OaNRW and δ2HaNRW (e.g., He and Richards, 2015), which was not the case in our study: we observed higher δ18OfNRW but lower δ2HfNRW as compared to δ18OaNRW and δ2HaNRW (Fig. 7a, b). Re-evaporation of NRW droplets on foliage might have occurred but was not the dominant process that could have led to the observed isotopic differences between fNRW and aNRW. Guttation, the exudation at leaf edges, is a process without a phase change of liquid water and thus does not involve isotopic fractionation. Hence, δ18O and δ2H of guttation water should be identical to δ18O and δ2H of leaf water. In our study, we found significant among-species differences in δ18O and δ2H of leaf water (Table 2), most likely resulting from species-specific leaf water evaporation and root water uptake, which contrasts with the insignificance of among-species differences in δ18O and δ2H of fNRW. This suggests that plant water only has a minor effect on δ18O and δ2H of fNRW. Furthermore, when the soil water content is much lower than field capacity, as was the case during all three events studied here, guttation hardly occurs (Long, 1955). During all three events SWC was very low (17 %–20 %) and thus close to the wilting point (12 %–14 %) and much lower than field capacity (27 %–30 %) in the main rooting zone in 0–15 cm soil depth. From these considerations, we conclude that re-evaporation and guttation are of no concern at our site during dry spells, and only distillation constitutes an important component for NRW on foliage besides the dominant NRW from ambient water vapor during the events in our study.
5.3 Uncertainty assessment of partitioning non-rainfall water components
The uncertainty of partitioning non-rainfall water components arises from the difficulties of measuring or calculating the distillation amount, although the NRW amount from ambient water vapor can be easily and accurately measured by a hydrometric approach, e.g., using a lysimeter (Jacobs et al., 2006). Distillation is an internal recycling of water from soil to plant surfaces (Monteith, 1957), which cannot be captured by a lysimeter because the latter device measures the water budget of plant and soil monoliths (Agam and Berliner, 2006) and thus does not distinguish between water in the soil and water on plant leaves. The EC method is widely used to investigate the water flux dynamics in ecosystem, but its suitability for quantitative NRW estimates can be questioned when an open-path IRGA is used to measure LE in clear and calm nights with dew and radiation fog occurrences. As soon as fog occurs or dew drips to the optical windows of the IRGA, LE measurements become unrealistic and cannot be analyzed quantitatively. The use of a closed-path IRGA that does not suffer from this problem may be a solution but could not be tested at the Chamau site in this study. But even when LE measurements appear to be of high quality, the EC-derived NRW estimates from ambient water vapor may not be very accurate, as shown by Jacobs et al. (2006), who found that the EC approach obtained less than one-third of the NRW amount as compared to the NRW amount obtained by a lysimeter. Monteith (1957) gave the equations of calculating the distillation amount, but a reanalysis of the data he published revealed that only the order of magnitude of distillation (reported as 1–2 mg cm−2 h−1, which corresponds with 0.01–0.02 mm h−1) agreed reasonably with observations, and large uncertainties remained, most likely as a result of untestable assumptions that have to be made about molecular transfer, linear temperature gradient, and saturated vapor at the soil surface for the M57 method as shown in Eq. (19) to be valid (Monteith, 1957).
To overcome this problem, we used a two-end-member mixing model as an alternative approach to partition NRW components. We compared the results of partitioning NRW components by our mixing model with estimates computed using the M57 approach. In general, with our mixing model we obtain lower values of fdistillation in the second half of condensation periods (i.e., 03:00 CET of event 1, and 01:00 CET of event 2; Fig. 10b) in comparison with the results obtained from the M57 approach. In particular, in the second half condensation periods of events 2 and 3, when the ambient air reached saturation, the mixing model was not applicable for partitioning NRW into aNRW and distillation fractions because δ18O and δ2H were too similar between aNRW and fNRW. Under these conditions the NRW amount approach yielded 6 % offdistillation at minimum. This could be explained by the isotopic exchanges in the soil–plant–atmosphere continuum that attenuated the δ18O and δ2H differences among different water sources driven by molecular diffusivity (Farquhar and Lloyd, 1993; Dawson et al., 2002). During the nights of events 1 and 2, we observed a very large variability of soil water δ18O and δ2H (Fig. 8), suggesting that distillation is not a vapor source with constant δ18O and δ2H. Furthermore, the spatial variability of shallow soil water content and its isotopic composition might enlarge the variability of δ18O and δ2H for distillation. This might explain why fdistillation obtained from the mixing model differed more substantially from the M57 approach at 01:00 CET of event 2 than that at 23:00 of event 2. One reason might be the shortcoming that we had to assume constant values of δ18O and δ2H of distillation for estimating fdistillation during 23:00 and 01:00 CET of event 2. Therefore, we recommend more intensive sampling of NRW on foliage in future studies, e.g., every 30 min for δ18OfNRW and δ2HfNRW in Eqs. (13)–(18). This should help to improve the accuracy of the fdistillation estimate by the mixing model. Another reason might be that the NRW droplets we took from the different leaves represent the cumulated NRW, while the temporal variability of NRW droplets on foliage might enlarge the uncertainty of NRW partitioning.
The M57 approach was not accurate enough for computing fdistillation during the events studied here because air temperature at 1 cm a.g.l. (Ta1cm) had to be computed from soil and surface temperatures, as it was not directly measured. Consequently, the large uncertainty in our Ta1cm estimates translates to increased uncertainty in distillation estimates computed via Eq. (19). Thus, we have to assume that our calculated distillation rates using the M57 approach are even more than 32 % of the levels that Monteith (1957) measured via his filter paper sampling. In the M57 approach as shown in Eq. (22), the stability term ) was used. However, the stability term is sometimes written as for Ri<−0.1 and for as, e.g., in Monteith and Unsworth (2013), which would cause higher condensation rate using Eq. (22) (see Fig. E1 in Appendix E), hence a lower relative contribution of distillation in the total NRW than given the term ). In future research, we recommend combining isotopic composition measurements with lysimetric measurements to partition NRW from ambient water vapor and distillation. This would provide a useful benchmark to better evaluate the isotope-based estimates of NRW inputs. The NRW amount from ambient water vapor can be measured directly by a lysimeter as the net water gain of the soil and plant monoliths (Kaseke et al., 2012; Ucles et al., 2013), while distillation is an indirect estimate based on stable water isotope data of the transfer of moisture from one part of the surface (soil surface) to another (foliage) within grassland ecosystems.
5.4 Potential effects of non-rainfall water on local water cycling
From the perspective of ecological relevance, distillation might be more important than previously thought, although it has no large-scale hydrological significance as it is simply a moisture transfer from the soil to the atmosphere (Monteith, 1957) and was thought to be detrimental by enlarging soil water loss via providing a shortcut for the water transfer (Chaney, 1981). The ecological relevance of distillation can be expected if the transfer of moisture is from one hydrological pool that is inaccessible to plants to another that is actually accessible to plants. For example, distillation could transfer soil-diffusing vapor from layers deeper than the effective rooting zone of grassland to droplets forming or depositing on leaf surfaces or surface soil where it can be accessed by the fine roots. Wang et al. (2017) observed that 0.0092 mm of water was transferred from deeper soil layers to the surface by vapor diffusion in a grassland plot, although it was debated whether the water went onto foliage or was absorbed by the topsoil. The process of vapor diffusion from deeper soil layers to the surface strongly depends on soil properties and thus might differ from site to site. For event 2 in our study, the nighttime variability of soil water δ18O and δ2H was not only observed in the top 0–15 cm, but also in the deeper soil layers at 15–40 cm (Fig. 8), suggesting that isotopic exchange occurs between deeper soil and topsoil layers. Furthermore, when soil water content was low and close to the wilting point, soil vapor diffusion is expected to become more important for distillation than capillary rise. Theoretically, the soil vapor diffusion rate increases with the decrease of soil water content on the conditions of volumetric soil water content >10 % (Philip and De Vries, 1957; Barnes and Turner, 1998), which makes soil vapor diffusion more important under such conditions. This corresponded to our results of the isotopic composition of soil water that δ18Os and δ2Hs variability (Fig. 8) was stronger in events 1 and 2 (SD of δ18Os and δ2Hs was 4.1 ‰ and 20.1 ‰ respectively) with a slightly lower SWC than that in event 3 (1 %–3 % lower of SWC, Fig. 1d; SD of δ18Os and δ2Hs was 2.2 ‰ and 12.4 ‰ respectively) with a slightly higher SWC. Therefore, future research focusing on the continuous measurements of the isotopic composition (δ18O and δ2H) of soil vapor is expected to give more quantitative insights on vapor transfer in soils during nights with dew and radiation fog. According to the results from our mixing model, distillation contributed up to 42 % of the total NRW inputs (Fig. 10b) and was the important pathway of NRW inputs during very calm nights (u2m<0.7 m s−1; see also Monteith, 1957, for u2m<0.5 m s−1) besides the condensation of ambient water vapor. According to our mixing model, distillation contributed up to 42 % of the total NRW on foliage (Fig. 10b); according to the M57 approach, the distillation amount was 0.04–0.05 mm per night (Fig. 10a) as compared to the NRW from ambient water vapor (0.12–0.50 mm per night) for the three events in our study.
We investigated the small-scale processes of how fog and dew formation influence the water cycling over a central European temperate grassland site during representative warm-season nights. Our results revealed different input pathways for dew and radiation fog during three dry intensive observation periods in summer 2018. Dew and radiation fog occurred on clear, calm nights with very low wind speed (u2m<0.7 m s−1) and weak turbulence, with LE ranging from −26 to 14 W m−2. Three primary pathways of NRW gains during the night were investigated in detail: (1) the condensation of atmospheric water vapor to plants, which constitutes a net water gain and might be important during dry spells or droughts; (2) internal recycling by distillation of water vapor from soil onto foliage, thereby redistributing the water within the ecosystem with no net water gain; and (3) radiation fog droplet deposition, which also leads to a net water gain. Condensation of ambient water vapor during dew and radiation fog was found to be predominantly an equilibrium fractionation process, which was deduced from the rather constant da during NRW nights. A decrease of 0.8 ‰–1.6 ‰ δ2Ha h−1 was observed in ambient water vapor, induced by condensation under equilibrium conditions during dew and radiation fog. Local evapotranspiration at high relative humidity from 17:00 CET until sunset caused the lowering of da to values in the range of 2.4 ‰ to 4.8 ‰ as compared to the higher daytime da (12.2 ‰ to 18.0 ‰). A further decrease to da values in the range of −11.8 ‰ to −4.7 ‰ was observed during the occurrence of dew and radiation fog in the night. Dew only formed with unsaturated conditions with a mixed NRW condensing from ambient water vapor and soil-diffusing vapor (distillation). The comparison between the foliage NRW δ18OfNRW and δ2HfNRW and the equilibrium NRW δ18OaNRW and δ2HaNRW of ambient water vapor allowed us to trace the source of NRW input pathways during dew formation. Distillation contributed 9 %–42 % to the total foliage NRW computed from a two-end-member mixing model. The dew and radiation fog produced 0.17–0.54 mm d−1 NRW gain on foliage, computed from a near-surface vertical temperature gradient method proposed by Monteith in 1957, which constitutes a non-negligible water flux into the canopy compared to evapotranspiration of on average 2.7 mm d−1. The strong anticorrelation between da and local RH0 suggested only a minor influence of large-scale air advection and highlighted the dominant role of local moisture as a source for ambient water vapor. Our results thus underline the importance of NRW inputs to temperate grasslands during dry spells and reveal the complexity of the local water cycle in such conditions, including different pathways of dew and radiation fog water inputs. In future studies, more intensive and continuous isotope measurements of foliage NRW, ambient water vapor, and soil vapor should be combined with direct lysimeter measurements to partition the NRW components from ambient water vapor and soil-diffusing vapor.
The eddy covariance fluxes were calculated using the software EddyPro (version 7.0.6, LI-COR, 2017) and following established community guidelines (Aubinet et al., 2012). Eddy covariance raw data were despiked and screened following Vickers and Mahrt (1997). Wind data were rotated (2D rotation; Wilczak et al., 2001), and time lags between the turbulent wind and H2O data were compensated for using covariance maximization. For spectral corrections, fluxes were corrected for high-pass and low-pass filtering effects (Moncrieff et al., 2005; Fratini et al., 2012) and instrument separation (Horst and Lenschow, 2009). Processed H2O fluxes were rejected from further analyses (1) if they were found outside a physically plausible range (between −20 and 50 mmol H2O m−2 s−1) and (2) if they failed the tests for stationarity and well-developed turbulence (e.g., Foken et al., 2005).
The saturation specific humidity (q0, g kg−1) for wet (q0w) and dry (q0d) vegetation surfaces was calculated as (Garratt, 1992)
where p (in hPa) is air pressure, and es0 (in hPa) is saturation vapor pressure at T0 calculated as (Garratt, 1992)
The dew point temperature (Td, ∘C) was calculated as (Garratt, 1992)
where esa (in hPa) is saturation vapor pressure at Ta calculated as (Garratt, 1992)
The evapotranspiration rate (in mm h−1) was calculated from the turbulent latent heat flux (LE in W m−2) as (Stull, 1988)
where (Stull, 1988), kg m−3 is water density, and b is a unit conversion factor (3.6×106 mm m−1 s h−1). Negative values of ET indicate dew formation below the eddy covariance flux instrumentation.
The soil water potential (ψs in kPa) was calculated from soil water content (θs) as (Saxton et al., 1986)
where (% sand) and (% clay) are the percentage of sand and clay, respectively; and wilting point and field capacity were calculated given kPa and kPa, respectively (Rai et al., 2017).
The saturated water content was calculated as (Saxton et al., 1986)
On relatively windy nights when the leaf surfaces remained dry (Monteith, 1957; Garratt, 1992), the roughness length (z0) at the Chamau site was computed from wind speed (u2m) and friction velocity (u*) following Monin–Obukhov similarity theory (Monteith, 1957; Garratt, 1992) as shown in Eq. (9). We selected two relatively windy nights, i.e., neutral atmospheric stratification (Panofsky, 1984), from 22 to 23 June 2018, and from 1 to 2 July 2018, to calculate roughness length z0. During these two nights, no precipitation occurred, and the latent heat flux (LE) was purely upward (i.e., no condensation); therefore leaf surfaces remained dry. The average of roughness length z0 was thus 0.03 m (Table D1). No harvest occurred between these two nights and the three events (see Sect. 3.2), and the grassland height was 0.2–0.3 m; therefore the grassland growth causes a minor change of z0.
The stability parameter Φ is proportional to the Richardson number Ri with numerous semi-empirical forms (Garratt, 1992). In our results, we followed the term ) as suggested in the original publication by Monteith (1957) that we refer to. However, in Monteith and Unsworth (2013), Φ is given as
In Fig. E1, we assess the effect that a replacement of the original approach by Monteith (1957) (Fig. E1a,b) has if it is replaced by the variant found in Monteith and Unsworth (2013) (Fig. E1c,d). The change does not affect our estimates for distillation (red curves in Fig. E1a,c) but increases the cumulative NRW amount gained over each of the three events studied (blue curves in Fig. E1a,c) by 34 %, 75 %, and 43 %, respectively.
Data were deposited at the ETH Zurich research collection at https://doi.org/10.3929/ethz-b-000465064 (Li, 2020).
YL, AR, WE, and FA designed the project. YL and AR performed the field experiment. YL carried out the laboratory work. YL performed statistical analyses. FA wrote the code of data calibration for the ambient water vapor isotopes. FA and WE commented on the results of the data analysis. YL wrote and revised the manuscript, with contributions and feedback by FA, WE, AR, and NB.
The authors declare that they have no conflict of interest.
The Picarro L2130-i analyzer was provided by the Atmospheric Dynamics group (Heini Wernli) at ETH Zurich. We thank Iris Feigenwinter, ETH Zurich, for the quality check of EC and meteorological data as well as measurements of LAI; Kathrin Fuchs, KIT Garmisch-Partenkirchen (Germany), for her introduction to the Chamau (Chamau) field site; Lukas Hörtnagl, ETH Zurich, for the method introduction of EC data; and Stephanie Westerhuis for the coding of extracting COSMO-1 model data. The Canton of Zug staff are thanked for maintaining the agricultural management of the Chamau field site. The authors also acknowledge MeteoSwiss for the access to the COSMO-1 model data products. We also thank the editor and two anonymous reviewers, who provided very insightful and detailed comments that considerably improved this final version.
This research has been supported by the Swiss National Science Foundation (grant no. 175733).
This paper was edited by Ryan Teuling and reviewed by two anonymous referees.
Aemisegger, F., Sturm, P., Graf, P., Sodemann, H., Pfahl, S., Knohl, A., and Wernli, H.: Measuring variations of δ18O and δ2H in atmospheric water vapour using two commercial laser-based spectrometers: an instrument characterisation study, Atmos. Meas. Tech., 5, 1491–1511, https://doi.org/10.5194/amt-5-1491-2012, 2012.
Aemisegger, F., Pfahl, S., Sodemann, H., Lehner, I., Seneviratne, S. I., and Wernli, H.: Deuterium excess as a proxy for continental moisture recycling and plant transpiration, Atmos. Chem. Phys., 14, 4029–4054, https://doi.org/10.5194/acp-14-4029-2014, 2014.
Agam, N. and Berliner, P. R.: Dew formation and water vapor adsorption in semi-arid environments – A review, J. Arid Environ., 65, 572–590, https://doi.org/10.1016/j.jaridenv.2005.09.004, 2006.
Atzema, A. J., Jacobs, A. F. G., and Wartena, L.: Moisture distribution within a maize crop due to dew, Neth. J. Agr. Sci., 38, 117–129, https://doi.org/10.18174/njas.v38i2.16599, 1990.
Aubinet, M., Vesala, T., and Papale, D.: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, Springer Science & Business Media, Dordrecht, The Netherlands, 438 pp., https://doi.org/10.1007/978-94-007-2351-1, 2012.
Barnes, C. J. and Turner, J. V.: Isotopic exchange in soil water, in: Isotope tracers in catchment hydrology, edited by: Kendall, C. and McDonnell, J. J., Elsevier, Amsterdam, The Netherlands, 137–163, https://doi.org/10.1016/B978-0-444-81546-0.50012-4, 1998.
Berkelhammer, M., Hu, J., Bailey, A., Noone, D. C., Still, C. J., Barnard, H., Gochis, D., Hsiao, G. S., Rahn, T., and Turnipseed, A.: The nocturnal water cycle in an open-canopy forest, J. Geophys. Res.-Atmos., 118, 10225–10242, https://doi.org/10.1002/jgrd.50701, 2013.
Beysens, D.: Dew water, River Publishers, Denmark, 2018.
Beysens, D., Clus, O., Mileta, M., Milimouk, I., Muselli, M., and Nikolayev, V. S.: Collecting dew as a water source on small islands: the dew equipment for water project in Bisevo (Croatia), Energy, 32, 1032–1037, https://doi.org/10.1016/j.energy.2006.09.021, 2007.
Bigeleisen, J.: Statistical mechanics of isotope effects on the thermodynamic properties of condensed systems, J. Chem. Phys., 34, 1485–1493, https://doi.org/10.1063/1.1701033, 1961.
Bruijnzeel, L., Eugster, W., and Burkard, R.: Fog as a hydrologic input, in: Encyclopedia of hydrological sciences, John Wiley & Sons, Ltd, 1–24, https://doi.org/10.1002/0470848944.hsa041, 2006.
Buck, A. L.: New equations for computing vapor pressure and enhancement factor, J. Appl. Meteorol., 20, 1527–1532, https://doi.org/10.1175/1520-0450(1981)020<1527:NEFCVP>2.0.CO;2, 1981.
Campbell, G. S. and Norman, J. M.: Water vapor and other gases, in: An Introduction to Environmental Biophysics, edited by: Campbell, G. S. and Norman, J. M., Springer New York, New York, USA, 37–51, https://doi.org/10.1007/978-1-4612-1626-1_3, 1998.
Castillo, J. L. and Rosner, D. E.: Theory of surface deposition from a unary dilute vapor-containing stream allowing for condensation within the laminar boundary-layer, Chem. Eng. Sci., 44, 925–937, https://doi.org/10.1016/0301-9322(89)90088-8, 1989.
Chaney, W. R.: Sources of water, in: Woody plant communities, edited by: Kozlowski, T. T., Academic Press, 1–47, New York, USA, https://doi.org/10.1016/B978-0-12-424156-5.50007-X, 1981.
Clus, O., Ortega, P., Muselli, M., Milimouk, I., and Beysens, D.: Study of dew water collection in humid tropical islands, J. Hydrol., 361, 159–171, https://doi.org/10.1016/j.jhydrol.2008.07.038, 2008.
Craig, H. and Gordon, L. I.: Deuterium and Oxygen 18 variations in the ocean and the marine atmosphere, Consiglio nazionale delle richerche, Laboratorio de geologia nucleare Pisa, Spoleto, Italy, 9–130, 1965.
Dansgaard, W.: Stable isotopes in precipitation, Tellus, 16, 436–468, https://doi.org/10.1111/j.2153-3490.1964.tb00181.x, 1964.
Dawson, T. E.: Fog in the California redwood forest: ecosystem inputs and use by plants, Oecologia, 117, 476–485, https://doi.org/10.1007/s004420050683, 1998.
Dawson, T. E., Mambelli, S., Plamboeck, A. H., Templer, P. H., and Tu, K. P.: Stable isotopes in plant ecology, Annu. Rev. Ecol. Syst., 33, 507–559, https://doi.org/10.1146/annurev.ecolsys.33.020602.095451, 2002.
del Prado, R. and Sancho, L. G.: Dew as a key factor for the distribution pattern of the lichen species Teloschistes lacunosus in the Tabernas Desert (Spain), Flora, 202, 417–428, https://doi.org/10.1016/j.flora.2006.07.007, 2007.
Delattre, H., Vallet-Coulomb, C., and Sonzogni, C.: Deuterium excess in the atmospheric water vapour of a Mediterranean coastal wetland: regional vs. local signatures, Atmos. Chem. Phys., 15, 10167–10181, https://doi.org/10.5194/acp-15-10167-2015, 2015.
Denmead, O. T., Raupach, M. R., Dunin, F. X., Cleugh, H. A., and Leuning, R.: Boundary layer budgets for regional estimates of scalar fluxes, Global Change Biol., 2, 255–264, https://doi.org/10.1111/j.1365-2486.1996.tb00077.x, 1996.
Dongmann, G., Nurnberg, H. W., Forstel, H., and Wagener, K.: On the enrichment of HO in the leaves of transpiring plants, Radiat. Environ. Bioph., 11, 41–52, https://doi.org/10.1007/BF01323099, 1974.
Drobinski, P., Haeberli, C., Richard, E., Lothon, M., Dabas, A. M., Flamant, P. H., Furger, M., and Steinacker, R.: Scale interaction processes during the MAP IOP 12 south fohn event in the Rhine Valley, Q. J. Roy. Meteor. Soc., 129, 729–753, https://doi.org/10.1256/qj.02.35, 2003.
Duine, G. J., Hedde, T., Roubin, P., and Durand, P.: A simple method based on routine observations to nowcast down-valley flows in shallow, narrow valleys, J. Appl. Meteorol. Clim., 55, 1497–1511, https://doi.org/10.1175/JAMC-D-15-0274.1, 2016.
Epstein, M., Hauser, G. M., Fauske, H. K., Grolmes, M. A., Henry, R. E., and Leung, J. C.: Fog formation and deposition within laminar and turbulent natural-convection boundary-layers along cold vertical plates, Chem. Eng. Commun., 118, 163–187, https://doi.org/10.1080/00986449208936092, 1992.
Eugster, W. and Merbold, L.: Eddy covariance for quantifying trace gas fluxes from soils, SOIL, 1, 187–205, https://doi.org/10.5194/soil-1-187-2015, 2015.
Eugster, W. and Siegrist, F.: The influence of nocturnal CO2 advection on CO2 flux measurements, Basic Appl. Ecol., 1, 177–188, https://doi.org/10.1078/1439-1791-00028, 2000.
Farquhar, G. D. and Lloyd, J.: Carbon and oxygen isotope effects in the exchange of carbon dioxide between terrestrial plants and the atmosphere, in: Stable Isotopes and Plant Carbon-water Relations, edited by: Ehleringer, J. R., Hall, A. E., and Farquhar, G. D., Academic Press, San Diego, USA, 47–70, https://doi.org/10.1016/B978-0-08-091801-3.50011-8, 1993.
Foken, T., Göockede, M., Mauder, M., Mahrt, L., Amiro, B., and Munger, W.: Post-field data quality control, in: Handbook of micrometeorology: a guide for surface flux measurement and analysis, edited by: Lee, X., Massman, W., and Law, B., Springer, Dordrecht, The Netherlands, 181–208, https://doi.org/10.1007/1-4020-2265-4_9, 2005.
Franssen, H. J. H., Stockli, R., Lehner, I., Rotenberg, E., and Seneviratne, S. I.: Energy balance closure of eddy-covariance data: A multisite analysis for European FLUXNET stations, Agr. Forest Meteorol., 150, 1553–1567, https://doi.org/10.1016/j.agrformet.2010.08.005, 2010.
Fratini, G., Ibrom, A., Arriga, N., Burba, G., and Papale, D.: Relative humidity effects on water vapour fluxes measured with closed-path eddy-covariance systems with short sampling lines, Agr. Forest Meteorol., 165, 53–63, https://doi.org/10.1016/j.agrformet.2012.05.018, 2012.
Fritschen, L. J. and Doraiswamy, P.: Dew: An dddition of hydrologic balance of Douglas Fir, Water Resour. Res., 9, 891–894, https://doi.org/10.1029/WR009i004p00891, 1973.
Fuchs, K., Hörtnagl, L., Buchmann, N., Eugster, W., Snow, V., and Merbold, L.: Management matters: testing a mitigation strategy for nitrous oxide emissions using legumes on intensively managed grassland, Biogeosciences, 15, 5519–5543, https://doi.org/10.5194/bg-15-5519-2018, 2018.
Galewsky, J., Steen-Larsen, H. C., Field, R. D., Worden, J., Risi, C., and Schneider, M.: Stable isotopes in atmospheric water vapor and applications to the hydrologic cycle, Rev. Geophys., 54, 809–865, https://doi.org/10.1002/2015RG000512, 2016.
Gallagher, M. W., Nemitz, E., Dorsey, J. R., Fowler, D., Sutton, M. A., Flynn, M., and Duyzer, J.: Measurements and parameterizations of small aerosol deposition velocities to grassland, arable crops, and forest: Influence of surface roughness length on deposition, J. Geophys. Res.-Atmos., 107, AAC 8-1–AAC 8-10, https://doi.org/10.1029/2001JD000817, 2002.
Garratt, J. R.: The atmospheric boundary layer, Earth-Sci. Rev., 89–134, Cambridge, England., 336 pp., 1992.
Gat, J. R.: Stable isotope hydrology: Deuterium and oxygen-18 in the water cycle, International Atomic Energy Agency (IAEA), available at: https://inis.iaea.org/collection/NCLCollectionStore/_Public/13/677/13677657.pdf?r=1 (last access: 4 May 2021), 1981.
Gat, J. R.: Oxygen and hydrogen isotopes in the hydrological cycle, Annu. Rev. Earth Pl. Sc., 24, 225–262, https://doi.org/10.1146/annurev.earth.24.1.225, 1996.
Gay, L. W., Vogt, R., Bernhofer, C., and Blanford, J. H.: Flux agreement above a Scots pine plantation, Theor. Appl. Climatol., 53, 33–48, https://doi.org/10.1007/Bf00866409, 1996.
Gehre, M., Geilmann, H., Richter, J., Werner, R. A., and Brand, W. A.: Continuous flow 2HH and 18OO analysis of water samples with dual inlet precision, Rapid Commun. Mass Sp., 18, 2650–2660, https://doi.org/10.1002/rcm.1672, 2004.
Glickman, T. S. and Zenk, W.: Glossary of meteorology, American Meteorological Society (AMS), Boston, USA, 855 pp., 2000.
Goulden, M. L., Miller, S. D., and da Rocha, H. R.: Nocturnal cold air drainage and pooling in a tropical forest, J. Geophys. Res.-Atmos., 111, D08S04, https://doi.org/10.1029/2005JD006037, 2006.
He, S. and Richards, K.: The role of dew in the monsoon season assessed via stable isotopes in an alpine meadow in Northern Tibet, Atmos. Res., 151, 101–109, https://doi.org/10.1016/j.atmosres.2014.02.014, 2015.
Hiatt, C., Fernandez, D., and Potter, C.: Measurements of fog water deposition on the California central coast, Atmospheric and Climate Sciences, 2, 525–531, https://doi.org/10.4236/acs.2012.24047, 2012.
Hindman, E. E., Borys, R. D., and Demott, P. J.: Hydrometeorological significance of rime ice deposits in the colorado rockies, Water Resour. Bull., 19, 619–624, https://doi.org/10.1111/j.1752-1688.1983.tb02779.x, 1983.
Horita, J. and Wesolowski, D. J.: Liquid-vapor fractionation of oxygen and hydrogen isotopes of water from the freezing to the critical temperature, Geochim. Cosmochim. Ac., 58, 3425–3437, https://doi.org/10.1016/0016-7037(94)90096-5, 1994.
Horst, T. W. and Lenschow, D. H.: Attenuation of scalar fluxes measured with spatially-displaced sensors, Bound.-Lay. Meteorol., 130, 275–300, https://doi.org/10.1007/s10546-008-9348-0, 2009.
Huang, L. and Wen, X.: Temporal variations of atmospheric water vapor δD and δ18O above an arid artificial oasis cropland in the Heihe River Basin, J. Geophys. Res.-Atmos., 119, 11456–11476, https://doi.org/10.1002/2014JD021891, 2014.
Hughes, R. N. and Brimblecombe, P.: Dew and guttation – formation and environmental significance, Agr. Forest Meteorol., 67, 173–190, https://doi.org/10.1016/0168-1923(94)90002-7, 1994.
IAEA: Reference Sheet for VSMOW2 and SLAP2 International Measurement Standards, International Atomic Energy Agency (IAEA), Vienna, Austria, 8 pp., 2009.
Jacobs, A. F. G., Heusinkveld, B. G., and Berkowicz, S. M.: A simple model for potential dewfall in an arid region, Atmos. Res., 64, 285–295, https://doi.org/10.1016/S0169-8095(02)00099-6, 2002.
Jacobs, A. F. G., Heusinkveld, B. G., Kruit, R. J. W., and Berkowicz, S. M.: Contribution of dew to the water budget of a grassland area in the Netherlands, Water Resour. Res., 42, W03415, https://doi.org/10.1029/2005WR004055, 2006.
Jouzel, J., Russell, G. L., Suozzo, R. J., Koster, R. D., White, J. W. C., and Broecker, W. S.: Simulations of the HDO and HO atmospheric cycles using the NASA GISS general circulation model: The seasonal cycle for present-day conditions, J. Geophys. Res.-Atmos., 92, 14739–14760, https://doi.org/10.1029/JD092iD12p14739, 1987.
Kaseke, K. F., Mills, A. J., Brown, R., Esler, K. J., Henschel, J. R., and Seely, M. K.: A method for direct assessment of the “Non rainfall” atmospheric water cycle: input and evaporation from the soil, Pure Appl. Geophys., 169, 847–857, https://doi.org/10.1007/s00024-011-0328-9, 2012.
Kaseke, K. F., Wang, L. X., and Seely, M. K.: Nonrainfall water origins and formation mechanisms, Science Advances, 3, e1603131, https://doi.org/10.1126/sciadv.1603131, 2017.
Keeling, C. D.: The concentration and isotopic abundances of atmospheric carbon dioxide in rural areas, Geochim. Cosmochim. Ac., 13, 322–334, https://doi.org/10.1016/0016-7037(58)90033-4, 1958.
Kelliher, F. M., Leuning, R., and Schulze, E. D.: Evaporation and canopy characteristics of coniferous forests and grasslands, Oecologia, 95, 153–163, https://doi.org/10.1007/BF00323485, 1993.
Kidron, G. J. and Temina, M.: The Effect of Dew and Fog on Lithic Lichens Along an Altitudinal Gradient in the Negev Desert, Geomicrobiol. J., 30, 281–290, https://doi.org/10.1080/01490451.2012.672542, 2013.
Kidron, G. J., Herrnstadt, I., and Barzilay, E.: The role of dew as a moisture source for sand microbiotic crusts in the Negev Desert, Israel, J. Arid Environ., 52, 517–533, https://doi.org/10.1006/jare.2002.1014, 2002.
Lai, C. T. and Ehleringer, J. R.: Deuterium excess reveals diurnal sources of water vapor in forest air, Oecologia, 165, 213–223, https://doi.org/10.1007/s00442-010-1721-2, 2011.
Lee, H., Smith, R., and Williams, J.: Water vapour 18OO isotope ratio in surface air in New England, USA, Tellus B, 58, 293–304, https://doi.org/10.1111/j.1600-0889.2006.00191.x, 2006.
LI-COR: Eddy covariance processing software, Version 7.0.6 [Software], LI-COR, Inc, available at: https://www.licor.com/env/products/eddy_covariance/software.html (last access: 4 May 2021), 2017.
Li, Y.: Data from: The role of dew and radiation fog inputs in the local water cycling of a temperate grassland during dry spells in Central Europe, ETH Zurich, https://doi.org/10.3929/ethz-b-000465064, 2020.
Long, I. F.: Dew and guttation, Weather, 10, 128–128, https://doi.org/10.1002/j.1477-8696.1955.tb00170.x, 1955.
López, A., Molina-Aiz, F. D., Valera, D. L., and Peña, A.: Determining the emissivity of the leaves of nine horticultural crops by means of infrared thermography, Sci. Hortic.-Amsterdam, 137, 49–58, https://doi.org/10.1016/j.scienta.2012.01.022, 2012.
Malek, E., McCurdy, G., and Giles, B.: Dew contribution to the annual water balances in semi-arid desert valleys, J. Arid Environ., 42, 71–80, https://doi.org/10.1006/jare.1999.0506, 1999.
McHugh, T. A., Morrissey, E. M., Reed, S. C., Hungate, B. A., and Schwartz, E.: Water from air: an overlooked source of moisture in arid and semiarid regions, Sci. Rep.-UK, 5, 13767, https://doi.org/10.1038/srep13767, 2015.
Meng, Y. and Wen, X. F.: Characteristics of dew events in an arid artificial oasis cropland and a sub-humid cropland in China, J. Arid Land, 8, 399–408, https://doi.org/10.1007/s40333-016-0006-y, 2016.
Merlivat, L.: Molecular diffusivities of HO, HD16O, and HO in gases, J. Chem. Phys., 69, 2864–2871, https://doi.org/10.1063/1.436884, 1978.
MeteoSwiss Operational Applications within COSMO, available at: http://www.cosmo-model.org/content/tasks/operational/meteoSwiss/default.htm, last access: 4 May 2021.
Moene, A. F. and van Dam, J. C.: Transport in the Atmosphere-Vegetation-Soil Continuum, Cambridge University Press, Cambridge, UK, https://doi.org/10.1017/CBO9781139043137, 2014.
Moncrieff, J., Clement, R., Finnigan, J., and Meyers, T.: Averaging, detrending, and filtering of eddy covariance time series, in: Handbook of Micrometeorology: A Guide for Surface Flux Measurement and Analysis, edited by: Lee, X., Massman, W., and Law, B., Springer, Dordrecht, The Netherlands, 7–31, https://doi.org/10.1007/1-4020-2265-4_2, 2005.
Monin, A. S. and Obukhov, A. M.: Osnovnye zakonomernosti turbulentnogo peremesivanija v prizemnom sloe atmosfery (in Russian), Trudy Geofizicheskogo Instituta Akademiya Nauk SSSR, 24, 163–187, 1954.
Monteith, J. L.: Dew, Q. J. Roy. Meteor. Soc., 83, 322–341, https://doi.org/10.1002/qj.49708335706, 1957.
Monteith, J. L. and Unsworth, M. H. (Eds.): Principles of environmental physics, Principles of environmental physics, edn. 4, Academic Press, Boston, Oxford, UK, 422 pp., https://doi.org/10.1016/B978-0-12-386910-4.00018-4, 2013.
Moore, C. J.: A comparative study of radiation balance above forest and grassland, Q. J. Roy. Meteor. Soc., 102, 889–899, https://doi.org/10.1002/qj.49710243416, 1976.
Nadeau, D. F., Pardyjak, E. R., Higgins, C. W., Huwald, H., and Parlange, M. B.: Flow during the evening transition over steep Alpine slopes, Q. J. Roy. Meteor. Soc., 139, 607–624, https://doi.org/10.1002/qj.1985, 2013.
Oke, T. R.: Temperature Profile near Ground on Calm Clear Nights, Q. J. Roy. Meteor. Soc., 96, 14–23, https://doi.org/10.1002/qj.49709640703, 1970.
Oke, T. R.: Boundary layer climates, edn. 2, Routledge, London, UK, https://doi.org/10.4324/9780203407219, 2002.
Pan, Y. X., Wang, X. P., and Zhang, Y. F.: Dew formation characteristics in a revegetation-stabilized desert ecosystem in Shapotou area, Northern China, J. Hydrol., 387, 265–272, https://doi.org/10.1016/j.jhydrol.2010.04.016, 2010.
Panofsky, H. A.: Atmospheric turbulence, John Wiley & Sons, New York, USA, 1984.
Parish, O. O. and Putnam, T. W.: Equations for the determination of humidity from dewpoint and psychrometric data, NASA, Washington D.C., USA, available at: https://ntrs.nasa.gov/citations/19770009916 (last access: 4 May 2021), 1977.
Parkes, S. D., McCabe, M. F., Griffiths, A. D., Wang, L., Chambers, S., Ershadi, A., Williams, A. G., Strauss, J., and Element, A.: Response of water vapour D-excess to land–atmosphere interactions in a semi-arid environment, Hydrol. Earth Syst. Sci., 21, 533–548, https://doi.org/10.5194/hess-21-533-2017, 2017.
Pasquill, F.: Eddy diffusion of water vapour and heat near the ground, Proc. R. Soc. Lon. Ser.-A, 198, 116–140, https://doi.org/10.1098/rspa.1949.0090, 1949.
Pfahl, S. and Wernli, H.: Air parcel trajectory analysis of stable isotopes in water vapor in the eastern Mediterranean, J. Geophys. Res.-Atmos., 113, D20104, https://doi.org/10.1029/2008JD009839, 2008.
Philip, J. R. and De Vries, D. A.: Moisture movement in porous materials under temperature gradients, Eos-Transactions American Geophysical Union, 38, 222–232, https://doi.org/10.1029/TR038i002p00222, 1957.
Phillips, D. L., Newsome, S. D., and Gregg, J. W.: Combining sources in stable isotope mixing models: alternative methods, Oecologia, 144, 520–527, https://doi.org/10.1007/s00442-004-1816-8, 2005.
Prechsl, U. E., Gilgen, A. K., Kahmen, A., and Buchmann, N.: Reliability and quality of water isotope data collected with a lowbudget rain collector, Rapid Commun. Mass Sp., 28, 879–885, https://doi.org/10.1002/rcm.6852, 2014.
Prechsl, U. E., Burri, S., Gilgen, A. K., Kahmen, A., and Buchmann, N.: No shift to a deeper water uptake depth in response to summer drought of two lowland and sub-alpine C3-grasslands in Switzerland, Oecologia, 177, 97–111, https://doi.org/10.1007/s00442-014-3092-6, 2015.
R Core Team: A language and environment for statistical computing, available at: https://www.R-project.org/ (last access: 4 May 2021), 2020.
Rai, R. K., Singh, V. P., and Upadhyay, A.: Soil analysis, in: Planning and evaluation of irrigation projects, edited by: Rai, R. K., Singh, V. P., and Upadhyay, A., Academic Press, Texas, USA, 505–523, https://doi.org/10.1016/B978-0-12-811748-4.00017-0, 2017.
Roth, K.: Bodenkartierung und GIS-basierte Kohlenstoffinventur von Graslandböden: Untersuchungen an den ETH-Forschungsstationen Chamau und Früebüel (ZG, Schweiz), University of Zurich, Zurich, Switzerland, 2006 (in German).
Saxton, K. E., Rawls, W. J., Romberger, J. S., and Papendick, R. I.: Estimating generalized soil-water characteristics from texture, Soil Sci. Soc. Am. J., 50, 1031–1036, https://doi.org/10.2136/sssaj1986.03615995005000040039x, 1986.
Schreel, J. D. M. and Steppe, K.: Foliar Water Uptake in Trees: Negligible or Necessary?, Trends Plant Sci., 25, 590–603, https://doi.org/10.1016/j.tplants.2020.01.003, 2020.
Spiegel, J. K., Aemisegger, F., Scholl, M., Wienhold, F. G., Collett Jr., J. L., Lee, T., van Pinxteren, D., Mertes, S., Tilgner, A., Herrmann, H., Werner, R. A., Buchmann, N., and Eugster, W.: Temporal evolution of stable water isotopologues in cloud droplets in a hill cap cloud in central Europe (HCCT-2010), Atmos. Chem. Phys., 12, 11679–11694, https://doi.org/10.5194/acp-12-11679-2012, 2012.
Stieger, J., Bamberger, I., Buchmann, N., and Eugster, W.: Validation of farm-scale methane emissions using nocturnal boundary layer budgets, Atmos. Chem. Phys., 15, 14055–14069, https://doi.org/10.5194/acp-15-14055-2015, 2015.
Stocking, C. R.: Guttation and bleeding, in: Pflanze und Wasser/Water relations of plants, edited by: Adriani, M. J., Aslyng, H. C., Burström, H., Geiger, R., Gessner, F., Härtel, O., Huber, B., Hülsbruch, M., Kalle, K., Kern, H., Killian, C., Kisser, J. G., Kramer, P. J., Lemée, G., Levitt, J., Meyer, B. S., Mothes, K., Pisek, A., Ruttner, F., Stålfelt, M. G., Stiles, W., Stocker, O., Stocking, C. R., Straka, H., Thornthwaite, W. C., Troll, C., Ullrich, H., and Veihmeyer, F. J., Springer Berlin Heidelberg, Berlin, Heidelberg, Germany, 489–502, https://doi.org/10.1007/978-3-642-94678-3_25, 1956.
Stull, R. B. (Ed.): Stable boundary layer, in: An Introduction to Boundary Layer Meteorology, Springer, Dordrecht, Dordrecht, The Netherlands, 499–543, https://doi.org/10.1007/978-94-009-3027-8_12, 1988.
Sun, H. Z., Clark, T. L., Stull, R. B., and Black, T. A.: Two-dimensional simulation of airflow and carbon dioxide transport over a forested mountain – Part I: Interactions between thermally-forced circulations, Agr. Forest Meteorol., 140, 338–351, https://doi.org/10.1016/j.agrformet.2006.03.023, 2006.
Thurnherr, I., Kozachek, A., Graf, P., Weng, Y., Bolshiyanov, D., Landwehr, S., Pfahl, S., Schmale, J., Sodemann, H., Steen-Larsen, H. C., Toffoli, A., Wernli, H., and Aemisegger, F.: Meridional and vertical variations of the water vapour isotopic composition in the marine boundary layer over the Atlantic and Southern Ocean, Atmos. Chem. Phys., 20, 5811–5835, https://doi.org/10.5194/acp-20-5811-2020, 2020.
Tomaszkiewicz, M., Abou Najm, M., Zurayk, R., and El-Fadel, M.: Dew as an adaptation measure to meet water demand in agriculture and reforestation, Agr. Forest Meteorol., 232, 411–421, https://doi.org/10.1016/j.agrformet.2016.09.009, 2017.
Tuller, S. E. and Chilton, R.: The role of dew in the seasonal moisture balance of a summer-dry climate, Agr. Meteorol., 11, 135–142, https://doi.org/10.1016/0002-1571(73)90057-5, 1973.
Ucles, O., Villagarcia, L., Canton, Y., and Domingo, F.: Microlysimeter station for long term non-rainfall water input and evaporation studies, Agr. Forest Meteorol., 182, 13–20, https://doi.org/10.1016/j.agrformet.2013.07.017, 2013.
Venables, W. N. and Ripley, B. D.: Modern applied statistics with S, edn. 4, Springer, New York, USA, 498 pp., https://doi.org/10.1007/978-0-387-21706-2, 2002.
Vickers, D. and Mahrt, L.: Quality control and flux sampling problems for tower and aircraft data, J. Atmos. Ocean. Tech., 14, 512–526, https://doi.org/10.1175/1520-0426(1997)014<0512:QCAFSP>2.0.CO;2, 1997.
Wang, C., Cen, Y., Liu, M., and Bowler, P.: Formation and influencing factors of dew in sparse elm woods and grassland in a semi-arid area, Acta Ecologica Sinica, 37, 125–132, https://doi.org/10.1016/j.chnaes.2017.06.004, 2017.
Welp, L. R., Lee, X., Griffis, T. J., Wen, X.-F., Xiao, W., Li, S., Sun, X., Hu, Z., Val Martin, M., and Huang, J.: A meta-analysis of water vapor deuterium-excess in the midlatitude atmospheric surface layer, Global Biogeochem. Cy., 26, GB3021, https://doi.org/10.1029/2011GB004246, 2012.
Wen, X. F., Lee, X., Sun, X. M., Wang, J. L., Hu, Z. M., Li, S. G., and Yu, G. R.: Dew water isotopic ratios and their relationships to ecosystem water pools and fluxes in a cropland and a grassland in China, Oecologia, 168, 549–561, https://doi.org/10.1007/s00442-011-2091-0, 2012.
Werner, R. A. and Brand, W. A.: Referencing strategies and techniques in stable isotope ratio analysis, Rapid Commun. Mass Sp., 15, 501–519, https://doi.org/10.1002/rcm.258, 2001.
Westerhuis, S., Fuhrer, O., Cermak, J., and Eugster, W.: Identifying the key challenges for fog and low stratus forecasting in complex terrain, Q. J. Roy. Meteor. Soc., 146, 3347–3367, https://doi.org/10.1002/qj.3849, 2020.
Whiteman, C. D., Hoch, S. W., Lehner, M., and Haiden, T.: Nocturnal cold-air intrusions into a closed basin: observational evidence and conceptual model, J. Appl. Meteorol. Clim., 49, 1894–1905, https://doi.org/10.1175/2010JAMC2470.1, 2010.
Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic anemometer tilt correction algorithms, Bound.-Lay. Meteorol., 99, 127–150, https://doi.org/10.1023/A:1018966204465, 2001.
Williams, C. A., Reichstein, M., Buchmann, N., Baldocchi, D., Beer, C., Schwalm, C., Wohlfahrt, G., Hasler, N., Bernhofer, C., Foken, T., Papale, D., Schymanski, S., and Schaefer, K.: Climate and vegetation controls on the surface water balance: Synthesis of evapotranspiration measured across a global network of flux towers, Water Resour. Res., 48, W06523, https://doi.org/10.1029/2011WR011586, 2012.
Wilson, K., Goldstein, A., Falge, E., Aubinet, M., Baldocchi, D., Berbigier, P., Bernhofer, C., Ceulemans, R., Dolman, H., Field, C., Grelle, A., Ibrom, A., Law, B. E., Kowalski, A., Meyers, T., Moncrieff, J., Monson, R., Oechel, W., Tenhunen, J., Valentini, R., and Verma, S.: Energy balance closure at FLUXNET sites, Agr. Forest Meteorol., 113, 223–243, https://doi.org/10.1016/S0168-1923(02)00109-0, 2002.
Wolf, S., Eugster, W., Ammann, C., Häni, M., Zielis, S., Hiller, R., Stieger, J., Imer, D., Merbold, L., and Buchmann, N.: Contrasting response of grassland versus forest carbon and water fluxes to spring drought in Switzerland, Environ. Res. Lett., 8, 035007, https://doi.org/10.1088/1748-9326/8/3/035007, 2013.
Wyngaard, J. C.: Turbulence in the Atmosphere, Cambridge University Press, Cambridge, https://doi.org/10.1017/CBO9780511840524, 2010.
Xu, Y. Y., Yi, Y., Yang, X., and Dou, Y. B.: Using stable hydrogen and oxygen isotopes to distinguish the sources of plant leaf surface moisture in an urban environment, Water-Sui, 11, 2287, https://doi.org/10.3390/w11112287, 2019.
Yakir, D. and Sternberg, L. d. S. L.: The use of stable isotopes to study ecosystem gas exchange, Oecologia, 123, 297–311, https://doi.org/10.1007/s004420051016, 2000.
Yepez, E. A., Williams, D. G., Scott, R. L., and Lin, G.: Partitioning overstory and understory evapotranspiration in a semiarid savanna woodland from the isotopic composition of water vapor, Agr. Forest Meteorol., 119, 53–68, https://doi.org/10.1016/S0168-1923(03)00116-3, 2003.