Optical Microwave Scintillometer Evaporation Measurements over a Saline Lake in a Heterogeneous Setting in the Atacama Desert

Felipe Lobos-Roco aMeteorology and Air Quality, Wageningen University, Wageningen, Netherlands

Search for other papers by Felipe Lobos-Roco in
Current site
Google Scholar
PubMed
Close
,
Oscar Hartogensis aMeteorology and Air Quality, Wageningen University, Wageningen, Netherlands

Search for other papers by Oscar Hartogensis in
Current site
Google Scholar
PubMed
Close
,
Jordi Vilà-Guerau de Arellano aMeteorology and Air Quality, Wageningen University, Wageningen, Netherlands

Search for other papers by Jordi Vilà-Guerau de Arellano in
Current site
Google Scholar
PubMed
Close
,
Francisca Aguirre bDepartamento de Ingeniería Hidráulica y Ambiental, Pontificia Universidad Católica de Chile, Santiago, Chile

Search for other papers by Francisca Aguirre in
Current site
Google Scholar
PubMed
Close
,
Alberto de la Fuente cDepartamento de Ingeniería Civil, Universidad de Chile, Santiago, Chile

Search for other papers by Alberto de la Fuente in
Current site
Google Scholar
PubMed
Close
, and
Francisco Suárez bDepartamento de Ingeniería Hidráulica y Ambiental, Pontificia Universidad Católica de Chile, Santiago, Chile
dCentro de Desarrollo Urbano Sustentable, Santiago, Chile
eCentro de Excelencia en Geotermia de los Andes, Santiago, Chile

Search for other papers by Francisco Suárez in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Estimating lake evaporation is a challenge due to both practical considerations and theoretical assumptions embedded in indirect methods. For the first time, we evaluated measurements from an optical microwave scintillometer (OMS) system over an open-water body under arid conditions. The OMS is a line-of-sight remote sensing technique that can be used to measure the sensible and latent heat fluxes over horizontal areas with pathlengths ranging from 0.5 to 10 km. We installed an OMS at a saline lake surrounded by a wet-salt crust in the Salar del Huasco, a heterogeneous desert landscape in the Atacama Desert. As a reference, we used eddy covariance systems installed over the two main surfaces in the OMS footprint. We performed a footprint analysis to reconstruct the surface contribution to the OMS measured fluxes (80% water and 20% wet salt). Furthermore, we investigated the applicability of the Monin–Obukhov similarity theory (MOST), which was needed to infer fluxes from the OMS-derived structure parameters to the fluxes. The OMS structure parameters and MOST were compromised, which we mitigated by fitting MOST coefficients to the site conditions. We argue that the MOST deviation from values found in the literature is due to the effects of the surface heterogeneity and the nonlocal processes induced by regional circulation. With the available dataset we were not able to rule out instrument issues, such as additional fluctuations to the scintillation signal due to absorption or the effect of vibration in high-wind conditions. The adjusted MOST coefficients lowered by a factor of 1.64 compared to using standard MOST coefficients. For H and LυE, we obtained zero-intercept linear regressions with correlations, R2, of 0.92 and 0.96, respectively. We conclude that advances in MOST are needed to successfully apply the OMS method in landscapes characterized by complex heterogeneity such as the Salar del Huasco.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Felipe Lobos-Roco, felipe.lobosroco@wur.nl, felipe.lobos.roco@gmail.com

Abstract

Estimating lake evaporation is a challenge due to both practical considerations and theoretical assumptions embedded in indirect methods. For the first time, we evaluated measurements from an optical microwave scintillometer (OMS) system over an open-water body under arid conditions. The OMS is a line-of-sight remote sensing technique that can be used to measure the sensible and latent heat fluxes over horizontal areas with pathlengths ranging from 0.5 to 10 km. We installed an OMS at a saline lake surrounded by a wet-salt crust in the Salar del Huasco, a heterogeneous desert landscape in the Atacama Desert. As a reference, we used eddy covariance systems installed over the two main surfaces in the OMS footprint. We performed a footprint analysis to reconstruct the surface contribution to the OMS measured fluxes (80% water and 20% wet salt). Furthermore, we investigated the applicability of the Monin–Obukhov similarity theory (MOST), which was needed to infer fluxes from the OMS-derived structure parameters to the fluxes. The OMS structure parameters and MOST were compromised, which we mitigated by fitting MOST coefficients to the site conditions. We argue that the MOST deviation from values found in the literature is due to the effects of the surface heterogeneity and the nonlocal processes induced by regional circulation. With the available dataset we were not able to rule out instrument issues, such as additional fluctuations to the scintillation signal due to absorption or the effect of vibration in high-wind conditions. The adjusted MOST coefficients lowered by a factor of 1.64 compared to using standard MOST coefficients. For H and LυE, we obtained zero-intercept linear regressions with correlations, R2, of 0.92 and 0.96, respectively. We conclude that advances in MOST are needed to successfully apply the OMS method in landscapes characterized by complex heterogeneity such as the Salar del Huasco.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Felipe Lobos-Roco, felipe.lobosroco@wur.nl, felipe.lobos.roco@gmail.com

1. Introduction

In deserts, water is confined to specific areas such as saline lakes, rivers, wetlands, and salt flats, all of which act as highly localized evaporative pathways (Rosen 1994). For example, the mountain regions of the Atacama Desert (Altiplano) are formed by consecutive endorheic basins. These basins catch precipitation, which then enters the hydrological system. From this system, the water wells up at specific sites and produces relatively small shallow lakes surrounded by extended salt flats (Uribe et al. 2015). In this water-limited environment, evaporation is highly relevant since it represents the principal cause of water loss in the basins.

In arid regions, measuring evaporation over small saline lakes is particularly complicated for three main reasons. First is the surface heterogeneity, with strong contrasting surface fluxes (Suarez et al. 2020), requiring the necessity of a footprint analysis to interpret the measured evaporation. Second, the atmospheric boundary layer dynamics and their influence on the surface fluxes are affected by nonlocal effects such as entrainment and advection of heat, moisture, and momentum (Lobos-Roco et al. 2021). Last, measuring evaporation over open waters always represents a challenge since installing instruments in water is technically challenging (Tanny et al. 2008; Nordbo et al. 2011).

To estimate evaporation, various direct and indirect methods with different theoretical assumptions and footprints were used and depend on the spatial scale of the subject under investigation (Abtew and Melesse 2012; Shuttleworth 2008). Within this context, we evaluated the applicability of an optical microwave scintillometer (OMS) to estimate open-water evaporation in arid regions. A scintillometer consists of a transmitter and a receiver. The transmitter emits a light beam toward the receiver, which measures the intensity of the fluctuations (the scintillation effect) caused by the turbulent transport of sensible heat and moisture crossing the instrument’s path. The OMS combines two scintillometers with different wavelengths: the large aperture scintillometer (LAS) operates an optical light beam, while the microwave scintillometer (MWS) operates a microwave beam (Green et al. 2001; Kohsiek and Herben 1983). The optical beam is primarily sensitive to temperature fluctuations induced by the sensible heat flux, whereas the microwave beam is sensitive to both heat and water vapor (Hill et al. 1980). Note that Han et al. (2019) estimated lake evaporation using a LAS, which by itself yields only the sensible heat flux, which tend to be small over water. The technique used depends on an energy balance approach where the quality of the evaporation estimate depends on the net radiation measurements and hard to determine water storage flux.

The OMS presents several advantages related to its large footprint (10–100 km2) and its capacity to integrate heterogeneous landscapes (Lagouarde et al. 2002; Beyrich et al. 2012; Isabelle et al. 2020; Zhang et al. 2021; Meijninger et al. 2006; Evans et al. 2012; Ward et al. 2015). The OMS is specifically convenient for measuring evaporation over open water since the transmitter and receiver can be installed on opposite shores, concentrating the central part of the instrument’s footprint over the water body. However, the OMS method also presents disadvantages in the flux calculation due to its dependency on the Monin–Obukhov similarity theory (MOST; Monin and Obukhov (1954)) and its underlying assumptions, such as horizontal homogeneity and the absence of outer scale processes. MOST connects turbulence-related quantities to turbulent fluxes through standard formulations determined empirically in field experiments and assumed to be universally applicable. This assumption is questionable when nonlocal contributions to the exchange process intervene (Kooijmans and Hartogensis 2016; Li et al. 2012).

The Altiplano of the Atacama Desert region presents interesting conditions to evaluate the strengths and weaknesses of the OMS method. We performed a 10-day field experiment called E-DATA, which was aimed at understanding the processes that control the evaporation of the saline lake in the Salar del Huasco in the Altiplano of the Atacama Desert (Suarez et al. 2020; Lobos-Roco et al. 2021). Part of the E-DATA experiment included the critical evaluation of the OMS method to determine open-water evaporation in a heterogeneous setting. To this end, we installed an OMS over the saline lake along with a small network of three eddy covariance (EC) systems, including one in the water body, covering the heterogeneous surfaces in the Salar del Huasco to serve as a reference.

Within this framework, we addressed the following research question: how accurate are the evaporation measurements performed by an OMS over an open-water body under arid and windy atmospheric conditions? To answer this question, we integrated into our analysis 1) the dynamic surface and atmospheric boundary layer conditions and 2) the evaluation of the MOST functions used to derive the surface fluxes and the impacts of both footprints and MOST function on the measurement of surface fluxes.

2. Methodology

This section begins with an overview of the measurement principles of the OMS and the EC systems, followed by descriptions of the calculation steps required to obtain surface fluxes from the OMS and EC measurements. Last, we describe the model used for calculating the OMS and EC footprints.

a. OMS versus EC method

Figure 1 provides an overview of some key aspects of the OMS and EC methods. Figure 1a summarizes the processing chain of the OMS to obtain surface fluxes from the OMS light intensity fluctuations, expressed as the variance of logarithm of the intensity fluctuations, σln(I)2. The OMS beams experience fluctuations in the intensity of the light due to the eddies that cross them. The statistic that captures the intensity of these fluctuations, σln(I)2, is related to the structure parameter of the refractive index Cn2 (step I in Fig. 1a), which is a spatial statistic that quantifies the turbulence intensity. The turbulent eddies that mainly contribute to the scintillation effect are of one typical size, i.e., the aperture diameter, D (= 0.15 m) for the LAS and the Fresnel length scale, F=λL, for the MWS, where λ (= 1.8 mm) is the wavelength and L is the scintillometer pathlength (0.5 < L < 10 km). Typical values for F are 1–3 m. These eddy sizes typically correspond to the spectral inertial subrange (Kolmogorov 1941). The LAS is mainly sensitive to temperature fluctuations, whereas the MWS is sensitive to both temperature and humidity fluctuations brought about by the transport of turbulent eddies. Combining the measurements from both scintillometers allows to solve the temperature and humidity fluctuations in terms of the structure parameters of temperature and humidity, CT2 and Cq2 (step II in Fig. 1a). These structure parameters, along with additional measurements of wind speed U and an estimate of the roughness length z0, are related to the fluxes of H and LυE following MOST formulations (step III in Fig. 1a). In section 2b and 2c, we make the processing steps more explicit.

Fig. 1.
Fig. 1.

Flowchart of the steps used to calculate the sensible (H) and latent heat (LυE) fluxes using the (a) optical microwave scintillation method and (b) eddy covariance method. The large arrows depict how the turbulent eddies are advected from their source area (footprint) to the scintillometer and EC sensor paths. Roman numerals I, II, and III represent the OMS calculations steps described in section 2b. A picture of every instrument is depicted next to the schemes.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Figure 1b highlights the main differences of the EC method in determining surface fluxes with respect to the OMS method: fluxes are calculated directly (i.e., do not rely on MOST) from the time series of high-frequency measurements of the 3D wind vector fluctuations (u′, υ′, w′), as well as fluctuations of air temperature (T′) and water vapor (q′) as the covariance of the vertical wind speed, w, and T for the sensible heat flux ( HwT¯), with for the latent flux ( LυEwq¯) and with for the momentum flux ( τwu¯). In contrast to the scintillometer that only considers a limited range of inertial-range eddy sizes around one dominant eddy scale in the flux estimate (e.g., Tatarski 1961; Wang et al. 1978), the EC method integrates all eddy sizes that contribute to the flux transport. The EC footprint, i.e., the source area of the turbulent eddies measured, is typically smaller than the scintillometer as it samples turbulent eddies at one point rather than over a path.

b. OMS fluxes

In the following section, we will explain the three processing steps of the OMS outlined in Fig. 1a. The first processing step relates the raw scintillation statistic σln(I)2 (—) to Cn2 (m−2/3) for the optical and microwave scintillometer (step I in Fig. 1a). This relationship, which assumes that the OMS is only sensitive to eddies in the inertial subrange, is expressed through wave propagation theory in a turbulent medium (Tatarski 1961) from which it follows that
Cn,opt2=CoptDopt7/3L3σln(I)2,and
Cn,mw2=CmwFmw7/3L3σln(I)2,
where Fmw corresponds to the Fresnel length scale ( F=λL), Dopt is the optical scintillometer aperture diameter (0.15 m), and Copt (0.907) and Cmw (3.425) are coefficients determined for each experimental setup and rely on the ratio between F and D (e.g., Ward et al. 2015).
The second processing step determines the structure parameters of temperature, CT2 (K2 m−2/3), and water vapor, Cq2 [(kg kg−1)2 m−2/3] (step II in Fig. 1a), where the wavelength-dependent structure parameter of the refractive index ( Cn,λ2), can be expressed as follows (Hill 1997):
Cn,λ2=AT,λ2CT2T2+Aq,λ2Cq2q2+2AT,λAq,λCTqTq,
where AT,λ and Aq,λ are dimensionless coefficients that depend on wavelength, pressure, temperature, and specific humidity (Ward et al. 2013). Two versions of Eq. (3) can be defined, i.e., one for the LAS (λ = 880 nm) and one for the MWS (λ = 1.8 mm). Equation (3), however, contains three unknown variables ( CT2, Cq2, and CTq), so a third expression is needed to solve the equations. Here, two approaches can be followed. The first one is a method introduced by Lüdi et al. (2005), which uses the covariance between the LAS and MWS intensity fluctuations. This defines a cross structure parameter of the refractive index between the two wavelengths used, Cn,opt,mw. Thus, a third expression of Eq. (3) is provided, which allows us to solve the three unknowns. The method provided by Lüdi et al. (2005) is quite sensitive to errors in the correlation measurement. Therefore, we used the second method, elaborated by Hill (1997). This method prescribes a value for rTq, the correlation coefficient between temperature and humidity, which is related to CT2 and Cq2. The advantage of the Hill method is that the calculations are more robust, i.e., less sensitive to measurement errors, but good representative values for rTq must be known a priori (Stoffer 2018). Following the Hill method, the solution for CT2 and Cq2 from Eq. (3) is given by
CT2=(Aq,mw2Cn,opt2+Aq,opt2Cn,mw2+2Aq,optAT,mwrTqCn,opt2Cn,mw2)T2γ2,
Cq2=(AT,mw2Cn,opt2+AT,opt2Cn,mw2+2AT,optAT,mwrTqCn,opt2Cn,mw2)q2γ2,
where γ=AT,mw2Aq,opt2AT,opt2Aq,mw2.
The third and last processing step is to connect CT2 and Cq2 to the surface fluxes in the framework of MOST (step III in Fig. 1a), which states that turbulence variables can be made dimensionless using MOST scaling variables as a function of the stability parameter z/Lob:
fCT2(z/Lob)=z2/3CT2θ*2,
fCq2(z/Lob)=z2/3Cq2q*2,
where θ* and q* are turbulent scaling variables for temperature and specific humidity and z/Lob is a dimensionless stability parameter where is the measurement height and Lob the Obukhov length. fCT2 and fCq2 are empirical functions that have the following base shape for unstable conditions (Wyngaard et al. 1971):
fCX2=c1(1c2zLob)2/3,
where X stands for T or q. Parameters c1 and c2 are coefficients that depend on stability and are considered universal when the standard assumptions for MOST of homogeneous and isotropic turbulence are met. However, no unanimity exists on the universality of fCT2 and fCq2 coefficients (Beyrich et al. 2012), which indicates that MOST assumptions are often violated. Kooijmans and Hartogensis (2016) provided an overview of the and coefficients and reported in the literature and they proposed their own values based on eleven data field experiments performed with the same instrumentation and with uniform data processing (see Table 2).
Sensible heat flux is determined as defined by the MOST framework as H=ρcpu*θ*, where ρ is the air density, cp the specific heat capacity at constant pressure, u* the friction velocity. Here, θ* is solved iteratively from Eq. (7) in conjunction with expressions for the friction velocity u* taken from standard flux profile relationships and the definition of the Obukhov length Lob, given by
u*=kUln(zUz0)Ψm(zLob)+Ψm(z0Lob),and
Lob=Tu*2kgθ*,
respectively, where k is the von Kármán constant (taken as 0.4), U is the wind speed measured at z, z0 is the surface roughness length, Ψm are stability correction functions (Businger et al. 1971; Dyer 1967), T is the air temperature, and g is the gravity acceleration. Once u* and Lob are solved as part of the sensible heat flux calculation, we can solve q* directly from Eq. (7) and from there on, the latent heat flux follows from its definition in the MOST framework ( LυE=ρLυu*q*).

The OMS data were processed from raw data to fluxes, according to the steps outlined in Fig. 1a. The raw 1-kHz intensity signals were filtered using a high-pass filter (31 s for the LAS and 71 s for the MWS for low crosswind conditions and 10 and 20 s for high crosswind conditions, respectively) to avoid the contribution of absorption fluctuations to the scintillation statistic, σln(I)2. Structure parameters of T and q were calculated according to the methodology proposed by Hill (1997) where we assumed, based on the EC measurements, the values of rTq were +0.7 for the day and +0.3 for the night. The rTq = +0.7 is a typical daytime value reflecting the strong correlation between T and q fluctuations as they are mostly being transported by the same turbulent eddies. The positive rTq values for the nighttime are atypical, which is caused by small, but positive H and LυE during the nighttime. Additionally, we used AT and Aq coefficients (Ward et al. 2015) evaluated with data from a Vaisala WXT520 all-weather sensor which was supplied with the Radiometer Physics GmbH (RPG) MWS. We used the internal data acquisition system of the RPG-MWS receiver to record the raw OMS (LAS and MWS) data (see also section 3a). The data processing was made using our own code in MATLAB and Python scripts, following the steps depicted in this section. Fluxes were computed using the MOST functions for unstable conditions defined by Kooijmans and Hartogensis (2016), but also using site-fitted MOST functions presented in the results section. Parameter z0 was solved from the wind flux profile relations [Eq. (9)] with prescribed u* from the EC systems [see section 3b(4)]. The OMS fluxes were calculated at a temporal resolution of 10 min.

c. EC fluxes

EC data were processed using the flux software package EddyPro v6.2.2 (Fratini and Mauder 2014) from LI-COR Biosciences Inc. (Lincoln, Nebraska, United States) to obtain the turbulent fluxes of latent and sensible heat, and the momentum at a temporal resolution of 10 min. All necessary standard data treatment and flux correction procedures were included, such as axis rotation with the planar-fit procedure (Wilczack et al. 2001), raw data screening including spike removal (Vickers and Mahrt 1997), interval linear detrending, and low-pass filtering correction (Massman 2000). Structure parameters from EC data were estimated using the structure-function equation (Hartogensis 2006): CX2=[x(r1)x(r2)]¯/r1,22/3 with x = q or T. The spatial separation r1,2 was approximated through the wind speed and time step (UΔt) using nominal values of 0.5 and 1 m. An average CX2 for the two separation distances was used in this analysis. For a more detailed explanation of the eddy covariance instrumentation and data processing, see Suarez et al. (2020) and Lobos-Roco et al. (2021).

d. OMS and EC footprints

To quantify the source area that determines the measurements over heterogeneous surfaces, a footprint model is a useful tool (Meijninger et al. 2002). In simple words, it represents the “field view” of an instrument defined by an upwind area as a function of measurement height (Schuepp et al. 1990), the current wind regime, surface fluxes and characteristics (Leclerc and Foken 2014). Footprints are especially relevant when measurements are performed below the so-called blending height (Wieringa 1976). The footprints are calculated using the advection–diffusion model proposed by Horst and Weil (1992), which was derived for a point measurement such as the EC method. To extend this to a scintillometer path footprint, the path between transmitter and receiver is discretized in N steps, where the point footprint at each location xi is convoluted with the path weighting function (Hartogensis et al. 2003) of the LAS and MWS (Meijninger et al. 2002) resulting in a representative source or footprint area (Spath):
Spath=i=1NWpath(xi)Spoint(xi,y,z,z0,zLob,U,συ),
where Spoint represents the point-source function (Neftel et al. 2008) as a function of the horizontal positions xi and y, the measurement height z, the Monin–Obukhov dimensionless stability parameter z/Lob, the wind speed U, and συ as the standard deviation of the lateral wind speed component that is used to model the Gaussian lateral spread. Finally, the footprints have been truncated to a 95% cumulative contribution.

3. Site and data description

This section describes the E-DATA experiment where our OMS and EC data were gathered. Moreover, we provide a detailed description of the atmospheric and surface conditions observed during the E-DATA experiment expressed by the boundary layer development, wind regime, surface heterogeneity, footprints, and roughness length.

a. E-DATA experiment

The data used in this study were gathered during the E-DATA experiment (Suarez et al. 2020; Lobos-Roco et al. 2021) performed at the saline lake of the Salar del Huasco (20.2°S, 68.8°W, 3790 m MSL) in the Altiplano of the Atacama Desert from 13 to 24 November 2018. The Salar del Huasco is representative of a typical natural hydrological system of the Altiplano region of Chile, characterized by a closed basin surrounded by complex topography of volcanic origin (Fig. 2a). Its hydrological functioning is given by the balance between occasional convective showers in the summer and a long dry season during the rest of the year (de la Fuente et al. 2021). Even so, these basins sustain permanent wetlands and feed the saline lake continuously with water through seepage from a complex hydrogeological system, making Salar del Huasco one of the few pathways of evaporation in the region (Blin et al. 2021).

Fig. 2.
Fig. 2.

Study area. (a) Salar del Huasco basin and saline lake. (b) Installation of the OMS and EC systems. Orange circles represent the OMS receiver (RX) and transmitter (TX), solid black line the OMS path and yellow circles the EC systems. Windrose indicates wind speed and direction of 15 Nov 2018.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

E-DATA was designed to study the evaporation in arid environments at local and regional scales over heterogeneous surfaces (Suarez et al. 2020). Part of the experiment was to collect data with an OMS installed over a ∼15-cm-deep saline lake in the middle of the salt flat (see Fig. 2b). The OMS was in operation from 13 to 24 November 2018. We experienced problems with the data collection for two reasons. First, there were technical problems with the storage of the raw data. The internal computer of the RPG-MWS receiver was not powerful enough to store the raw data at 1 kHz, which resulted in only 20%–35% of the samples being stored with no information on the nature of the gaps. The consequence of this is that we were not able to analyze the spectra of the raw scintillation signal to exclude potential influence of absorption fluctuations that are supposedly filtered out by the aforementioned high-pass filter and the potential influence tripod vibrations in high-wind conditions (Aguirre et al. 2022). If not filtered out adequately both these effects lead to an increase in σln(I)2 when the atmospheric conditions are favorable for these effects, i.e., a slow fluctuating humidity field and strong winds. As a result, the fluxes will be overestimated. Second, after 2.5 days into the experiment, the level of σln(I)2 increased by a factor of 4 and stayed at that level for the rest of the experiment. We hypothesize that misalignment of the beam caused ground reflections and subsequently double the amount of scintillations observed, which would increase the scintillation statistic by a factor of 4. We attribute the slight misalignment to the soft, unstable wet-salt and mud surface. The result was that we had one day with acceptable data left (15 November 2018) with overlapping OMS and EC data. This is a restriction but, as reported by Suarez et al. (2020) and Lobos-Roco et al. (2021), all days in the campaign showed very similar behavior in surface fluxes (see Fig. 5). Therefore, we are confident that what we learned from our one golden day is representative of a typical day in late spring in our study area.

Table 1 and Fig. 2b show the specific details of the OMS and EC systems installed. The OMS pathlength between the transmitter (Tx: 20.283°S, 68.870°W) and receiver (Rx: 20.280°S, 68.880°W) was 1070 m. We intended to place the OMS transmitter and receiver close to the shoreline ensuring an all-water footprint. In the end, we chose to install the OMS some distance from the shorelines where there was a somewhat higher (1 m) and firmer waterlogged wet salt. In all, the OMS path covered 300 m of wet salt and 770 m of water. The path crossed the saline lake from the northwest to southeast, offset northward by 45° with respect to the main wind speed (Fig. 2b). The transmitter and receiver were placed at 1.8 m above the wet salt, which was 0.5 m higher than the water surface, resulting in an OMS beam height of 2.3 m above the water. It is noteworthy to mention that the MWS beam size in the middle of the path, marked by the Fresnel length, reached 1.4 m, indicating that a minor misalignment could already give rise to ground reflections of the MWS beam. In addition, EC systems were installed over the lake (20.277°S, 68.882°W) and over a wet-salt surface (20.283°S, 68.875°W). Finally, to connect the analysis of the surface fluxes with the ABL dynamics, radiosondes were launched on 19 November 2018 at 0900, 1200, 1500, and 1800 LT from the lake surroundings (20.35°S–68.90°W; 3953 m MSL). The radiosondes were launched for measuring vertical profiles of potential temperature (θ) and specific humidity (q) within the first 2 km from the surface.

Table 1

Sensors and variables measured by the OMS and EC in the E-DATA experiment: large aperture scintillometer (LAS), microwave scintillometer (MWS), structure parameter of temperature ( CT2) and humidity ( Cq2), sensible heat flux (H), latent heat flux (LυE), friction velocity ( u*), air temperature (T), ∼5-cm temperature (TS), relative humidity (RH), air pressure (p), parallel wind component (u), lateral wind component (υ), vertical wind component (w), specific humidity (q), potential temperature (θ), incoming shortwave (Sw↓), outgoing shortwave (Sw↑), incoming longwave (Lw↓), outgoing longwave (Lw↑) radiation, and net radiation (Rn).

Table 1

b. E-DATA atmospheric and surface conditions

This subsection describes the atmospheric and surface conditions observed during the E-DATA experiment to provide a context for interpreting the scintillometer flux measurements. We will focus on two aspects. First, we will focus on the nonlocal processes associated with the regional atmospheric circulation that affect the atmospheric boundary layer development through advection and entrainment. Second, we will focus on the surface heterogeneity between water and wet salt in terms of near-surface temperature, net radiation, surface fluxes, roughness length, and how these areas are represented in the OMS footprint.

1) Atmospheric boundary layer

The nonlocal processes that define the land–atmosphere exchange in the Salar del Huasco are governed by a regional circulation pattern between the wet, cool, and low-elevation Pacific Ocean and the dry, hot, and high elevation of the Atacama Altiplano where the Salar del Huasco is located. This circulation is responsible for the two well-marked wind regimes observed over the entire Salar basin during the experiment (Fig. 3). The first wind regime corresponds to the relatively calm conditions that occur in the morning (0600–1200 LT) and is characterized by low wind speeds (<2 m s−1). During this regime, there is not a prevalent wind direction. The second wind regime corresponds to strong turbulent conditions that occur from noon onward (1200–2100 LT), which are characterized by high wind speeds (>10 m s−1). During this regime, there is a predominant westerly wind direction (250°–290°). The nature of this well-defined circulation has been analyzed in detail using observations and regional-scale modeling studies (Lobos-Roco et al. 2021; Suarez et al. 2020; Munoz et al. 2018; Rutllant et al. 2003).

Figure 4 shows the diurnal evolution of the potential temperature θ (top panels) and the specific humidity q (bottom panels) vertical profiles within the ABL. These profiles show a thermally stable ABL at 0900 LT, which rapidly develops to 1800 m at 1200 LT, becoming ∼15 K warmer with respect to the morning. The morning boundary layer development is dominated by the surface and entrainment fluxes leading to a reduction of the inversion layer jump to ∼1 K. The rapid boundary layer growth is interrupted after midday, where the ABL decays abruptly to 500 m. The decrease in the ABL height is associated with the advection of relatively cold air from the Pacific Ocean to the study area by a regionally driven wind regime introduced in Fig. 3 (see also Lobos-Roco et al. 2021; Suarez et al. 2020).

A similar pattern is observed for moisture in the bottom panels of Fig. 4. Here, a layer with a negative humidity gradient grows fast in the morning until it is interrupted in the early afternoon by the regional flow. The moisture profile shows that the air in the shallow boundary layer caused by the regional flow is relatively moister, i.e., we observe an increase of the specific humidity from <1 to 3.5 g kg−1.

2) Surface heterogeneity

To demonstrate the level of surface heterogeneity of the area covered by the OMS and EC measurements, we show time series of near-surface temperature, net radiation (Rn), and turbulent surface fluxes measured at the EC stations over the water and wet-salt surfaces (Fig. 5). In general terms, during the day, the near-surface (∼5-cm height) temperature is lower for the water surface (see Fig. 5b), so less longwave outgoing radiation will be lost to the atmosphere. In addition, a considerably higher albedo is observed over the wet salt (0.58) compared to the water surface (0.12) (Lobos-Roco et al. 2021). As a result of the differences in surface temperature and albedo during the afternoon, much more radiation energy is available at the water’s surface, which shows maximum values of ∼900 W m−2 compared to the maximum values of ∼500 W m−2 over the wet-salt surface. These remarkable differences in Rn have repercussions in the partitioning of the sensible and latent heat fluxes. Figures 5a and 5b also show the surface fluxes of latent heat (LυE) and sensible heat (H) over wet salt and water averaged along the E-DATA period. The first characteristic observed is the day-to-day low variability (shades) in the surface fluxes, evidence that the same diurnal cycle occurs every day. Even though both surfaces are water saturated, their flux characteristics are completely different. To begin with, the fluxes over the wet salt are smaller by a factor of ∼3 compared to the water surface given the much-reduced amount of net radiation available. In addition, H > LυE for the wet salt, whereas LυE > H for the water surface. This is remarkable since the wet salt crust is waterlogged. The hypersalt concentrations, however, impede the water escaping from the surface (e.g., Salhotra et al. 1985; Oroud 1999). As a result, most of the energy is dedicated to heating the air. Two aspects are equal for both surfaces. First is the sudden increase of fluxes at 1200 LT aligned with the arrival of the strong winds that break the virtually absent turbulent mixing seen in the morning (Lobos-Roco et al. 2021). Second is the positive sensible heat flux during both the day and night, meaning that the atmosphere is unstable even during nighttime conditions. Lobos-Roco et al. (2021) argue that the unstable nighttime condition is related to a katabatic flow that drains cold air from the surrounding mounting ridges into the Salar. However, it also might be related to the high heat capacity of the surface, causing heating of near-surface air further into the night.

Fig. 3.
Fig. 3.

Diurnal relationship between wind speed (U) and wind direction (WD) obtained by EC over water and wet-salt surfaces.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Fig. 4.
Fig. 4.

Atmospheric boundary layer conditions measured during E-DATA on 19 Nov 2018. (top) Diurnal cycle of the potential temperature radiosounding profiles. (bottom) Diurnal cycle of specific humidity radiosounding profiles.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Fig. 5.
Fig. 5.

Diurnal-averaged time series of near-surface (∼5 cm) temperature, net radiation, and turbulent fluxes over the (a) wet-salt and (b) water surfaces during the E-DATA period. The picture behind exemplifies the wet-salt and water surface during the experiment. Shades indicate the minimum and maximum values across the E-DATA period.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

3) Footprint analysis

The source area or footprint is highly dependent on the wind regime (Fig. 3). Figures 6a–d shows the footprints for both scintillometers (LAS and MWS) and the EC systems (water and wet salt) at two representative times of each wind regime (0900 and 1600 LT, respectively). The footprints have been truncated to a 95% cumulative contribution. Both scintillometer and EC footprints are considerably smaller in the morning compared to the afternoon due the different wind and related flux regimes. Their orientation changes according to the wind direction. Under all conditions, the EC footprints cover only a small portion of the surface above which they are installed (Table 1). The scintillometer footprints, on the other hand, always cover a mix of both water and wet salt. On average, about 2/3 of the LAS and MWS footprints are covered by water and 1/3 wet salt. The difference in the path weighting function means that the LAS footprints are more weighted toward the center of the path (Evans and De Bruin 2011) and therefore cover a slightly larger area of water compared to the MWS footprints.

Fig. 6.
Fig. 6.

Footprints of the (a),(c) LAS and EC and (b),(d) MWS and EC over the water and wet-salt surfaces for representative hours of wind regimes on 15 Nov 2018. Colors (given without scale) represent the percentage contribution of that area to the total footprint. The wind roses characterize the wind regimes for the same day (top) between 0700 and 1200 LT and (middle) between 1200 and 2000 LT. The background image is a georeferenced mosaic based on a drone flight on 18 Nov 2018. (e) Diurnal cycle of LAS and MWS percentage footprint contributions from the water (solid lines) and wet-salt surface (dotted lines).

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

This systematic difference can also be seen in Fig. 6e, which expresses the time-dependent footprint proportions covering water and wet salt for the two scintillometers on 15 November 2018. Driven by the wind regimes, we observe two distributions of footprints. During daytime when the winds are generally from the west and they become strong during the afternoon, the distribution is roughly 80% water and 20% wet salt. In the evening and night, the contributions of both surfaces to the scintillometer footprints are much more variable depending on the wind direction with contributions from water varying between 50% and 80% and wet salt between 20% and 50%. The diurnal cycle of the OMS water/wet-salt footprint proportion is only relevant when surface fluxes observed over both water and wet-salt surfaces are significantly different (Fig. 6e). The strong difference in flux regimes between the two surfaces demonstrates the importance of constructing a footprint-weighted composite of the EC fluxes when comparing them with those of the scintillometers.

4) Roughness length

For the OMS flux calculations, an accurate estimate of z0, is crucial given the fact the measurements were conducted close to the surface (∼2-m height) and the afternoon winds were strong. Under these conditions, mechanically generated turbulence is a key process in determining the fluxes [see Eqs. (7)(10)]. The wet-salt and water surfaces are very smooth, and an appropriate z0 estimate is, therefore, not trivial. We solved z0 from the wind flux profile relation [Eq. (9)] using the estimates of U, u*, and from the EC stations of water and wet salt. The results, presented in Fig. 7, show that for low wind speeds, z0 is ill defined and for high wind speeds, its value converges to a value indicated by the solid lines (Mahrt et al. 2001). On average, the water surface is smoother with z0 = 4 × 10−5 m against z0 = 2 × 10−4 m for the wet-salt surface. The shallow water becomes choppy in the afternoon, but this is not intense enough to cause z0 to be a function of wind speed (Charnock 1955). A weighted value considering the typical afternoon footprints over these two surfaces yields a value of z0 = 1.5 × 10−4 m that we used in the flux calculations.

Fig. 7.
Fig. 7.

Relationship between roughness length (z0) and wind speed over (a) water surface and (b) wet-salt surface. Wind speed data are taken from the EC systems. The solid lines express the values to which z0 converges for high wind speeds.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

4. Results and discussion

a. MOST functions

To assess the representativity of MOST for the special conditions in the Salar del Huasco, Fig. 8 shows the CT2 and Cq2 dimensionless groups [Eqs. (6) and (7)] as a function of z/Lob for EC-based structure parameters (Figs. 8a,b) and OMS structure parameters (Figs. 8c,d). The EC-based structure parameters are composed values from the wet-salt and water EC systems and reflect data from the entire 10-day E-DATA experiment, whereas the OMS structure parameters only include data from 15 November 2018. The data points in Fig. 8 were filtered by excluding values for which |H| < 5 W m−2. In addition, in Figs. 8a and 8b the CT2 and Cq2 dimensionless groups were filtered to the 25%–75% interquartile range. The structure parameters from the EC systems along with the EC-flux-derived variables θ*, q*, and Lob are weighted by the OMS water/wet-salt footprint proportions (wfp subscript). The data points are marked by the wind regime where, not surprisingly, the points associated with the strong wind conditions are on the neutral side of the −z/Lob range.

Fig. 8.
Fig. 8.

Dimensionless CT2 and Cq2 groups vs stability parameter (z/Lob). (a),(b) CT2 and Cq2 dimensionless group using parameters from EC water and wet salt weighted by the OMS footprint (wfp) during the E-DATA period. (c),(d) CT2 and Cq2 dimensionless group using structure parameters measured by the OMS and the remaining parameters from EC water and wet salt weighted by the OMS footprint (wfp) during 15 Nov 2018. MOST functions proposed by Kooijmans and Hartogensis (2016) and the proposed ones in this study are shown in solid and dashed lines, respectively. The experimental data points are marked by their wind regime (windy > 4 m s−1 between 1300 and 2040 LT and calm < 4 m s−1 between 2340 and 1300 LT).

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Figure 8a shows that standard MOST functions slightly overestimate the CT2 data points during the windy regime, but data are well represented during the calm wind conditions. Likewise, Fig. 8b shows that standard MOST functions represent both windy and calm regimes well for the Cq2 dimensionless group data points. Figure 8c shows that for CT2, the standard MOST functions of Kooijmans and Hartogensis (2016) cover the data points reasonably well under both wind regimes. However, for Cq2 (Fig. 8d), the data points follow the standard MOST functions during the calm regime, whereas for the windy regime, the data points are far above both the Kooijmans and Hartogensis (2016) function and the EC data points (shown in Fig. 8b). As our dataset did not allow us to discard the influence of instrumental issues that might explain the elevated OMS Cq2 we will, for the sake of the argument, assume that they do not play a role and discuss a physical process that plays a role here. The OMS and EC measure below the blending height. As the EC footprints are entirely above water and wet-salt surfaces (Fig. 6), no interaction between their respective internal boundary layers affects the measurements. However, in the case of OMS, 2/3 of the footprint is covered by water and 1/3 by wet salt, with an order of magnitude difference of LυE flux regime (Fig. 5). This high contrast contributes to having strong humidity fluctuations within the internal boundary layers at the transition between the two surfaces, elevating the measured Cq2 by the OMS. For CT2 (Fig. 8c), higher values are observed for the OMS as compared to the EC. However, the difference is less pronounced since the contrast between H flux surfaces is less as well. Note that both EC and OMS are under a setting governed by a regional advection and entrainment of air not originating from the Salar (Lobos-Roco et al. 2021). Such regional circulation contributes to 40% of the total diurnal moisture budget above the water surface and this phenomenon might cause additional humidity fluctuations. However, our data do not show evidence of this since we observe high differences between the Cq2 of EC and OMS, where the same conditions are supposed to affect both instruments. It is also under these conditions of no horizontal homogeneity and regional advection that MOST fails.

To mitigate these effects, we adopted a practical approach to eye-fit an alternative Cq2 MOST relation that inherently considers the enhanced Cq2. Our motivation to use an eye-fit rather than a statistical approach, like the one performed by Kooijmans and Hartogensis (2016), is that we only have a single day of measurements available, whereas other methods require significant amounts of data. For CT2, we use the coefficients determined by Kooijmans and Hartogensis (2016). The proposed line is depicted in Fig. 8d and the coefficients that define it are given in Table 2.

Table 2

Kooijmans and Hartogensis (2016) MOST function coefficients for unstable conditions [base function in Eq. (8)]. Numbers in parentheses are the site-adapted coefficients for the Salar del Huasco study case, 15 Nov 2018 (Fig. 8d).

Table 2

b. Impact of footprint and MOST on surface fluxes

In this section, we demonstrate that it is essential to consider the impacts of the footprint and MOST functions when comparing the performance of the OMS in a complex environment like the Salar del Huasco.

Figure 9 shows the diurnal cycle of surface fluxes measured by EC over water and wet-salt surfaces as well as the diurnal cycle combining both EC system fluxes, but weighting them using the OMS footprint (ECwfp) presented in Fig. 6e. ECwfp resembles most of the EC fluxes over water (80%) than wet salt (20% wet salt) during daytime. Given the high flux contrast between wet salt and water, the 20% wet salt seen in the OMS footprint significantly influences the composed fluxes. This is the case for LυE, where during daytime, it is ∼75 W m−2 over wet salt versus ∼550 W m−2 over water, the ECwfp LυE is ∼100 W m−2 lower than that of the water surface. For the sensible heat flux (Fig. 9c), the absolute values of the fluxes are lower and the contrast between the two surfaces is less extreme resulting in an ECwfp H that is only ∼10 W m−2 higher than that of the water surface. The orthogonal regression (Figs. 9b,d) between EC water, wet salt, and wfp shows a quantification of the differences in LυE and H of water and wet salt with respect to ECwfp.

Fig. 9.
Fig. 9.

Diurnal cycle of (a),(c) surface fluxes and (b),(d) orthogonal regression between EC over water and wet salt and an OMS footprint composed EC flux (ECwfp) for 15 Nov 2018. The “m” indicates the regression slope; gray data points refer to wet salt and black ones to water.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

When performing scintillometer measurements over heterogenous terrain under the blending height (Wieringa 1976), the structure parameters of the individual surfaces within the footprint will be sensed by the instrument (Meijninger et al. 2002). In our study, the OMS was installed close to the surface (2.3 m), and it is safe to assume we were below the blending height (Meijninger et al. 2002). Due to the nonlinear relationship between structure parameters and surface fluxes, there will be a difference between the EC flux weighted by the OMS footprint (ECwfp), and the EC flux calculated using structure parameters weighted by the OMS footprint ( ECCX2,wfp)(Lagouarde et al. 2002; Meijninger et al. 2006). The latter more accurately simulates the way of measuring the OMS. To quantify the difference between the two approaches, we constructed the OMS footprint-weighted, indirect EC flux estimate through CX2 ( ECCX2,wfp) following these steps: 1) obtain CT2 and Cq2 from the EC fluxes over water and wet salt applying an inverse MOST procedure using standard MOST function coefficients (Kooijmans and Hartogensis 2016); 2) construct simulated OMS-path-weighted structure parameters using the OMS-weighted footprint factors; and finally 3) calculate simulated OMS-path-weighted fluxes based from the path averaged structure parameters of step 2 using the same MOST functions of step 1.

Figure 10 shows the difference between the two approaches described above. The big contrast in Cq2 between water and wet salt, connected to their contrast in LυE (see Figs. 5a,b) and the nonlinear relationship between Cq2 and LυE, causes a difference of ∼13% between ECwfp and ECCX2,wfp. So, this is a significant factor to consider whenvalidating OMS derived fluxes with a footprint composite of EC fluxes. For the sensible heat fluxes, the difference is negligible, which is related to the much smaller fluxes with less contrast between the two surfaces (see Fig. 5).

Fig. 10.
Fig. 10.

Relationships of (a) latent and (b) sensible heat fluxes between OMS-footprint-weighted, direct EC flux estimates (ECwfp) and indirect EC flux estimates through CX2 ( ECCX2,wfp).

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Next, we will look at the difference in OMS-derived fluxes with standard MOST functions [Eq. (8)] taken from Kooijmans and Hartogensis (2016) and site-fitted MOST functions as determined in Fig. 8 and Table 2. Figure 11 shows an overestimation of LυE (factor of 1.64) if standard MOST functions are used. This is a clear indication that MOST assumptions are violated in the Salar del Huasco basin. In part, this has to do with the surface heterogeneity, noting that we measured below the blending height and measurements were thus sensitive to the non-MOST turbulence behavior at the transition between the two contrasting surfaces in the OMS footprint. In addition, regional advection of dry air has a strong influence on the local humidity (Lobos-Roco et al. 2021). Consequently, the relationship between the measured humidity fluctuations and local fluxes do not follow the standard relationships.

Fig. 11.
Fig. 11.

(a) Diurnal cycle and (b) orthogonal regression of latent heat flux from OMS data calculated with standard and site-fitted MOST functions. m indicates the regression slope.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

Figure 12 brings together the results presented so far, i.e., OMS fluxes based on site-fitted MOST functions (OMSsite-fitted-MOST), and ECwfp and ECCX2,wfp fluxes that serve as references. In general, we observe a reasonable agreement. The orthogonal regressions in Figs. 12b and 12d are based on ECCX2,wfp as this more closely resembles the OMS flux estimate given the fact that we measured below the blending height. After all the careful and guided treatment of the data, we arrive at zero-intercept linear regression with slope coefficients of 0.95 and 1.07 and correlation coefficients R2 of 0.92 and 0.96 for LυE and H, respectively.

Fig. 12.
Fig. 12.

Diurnal cycle of (a),(c) surface fluxes and (b),(d) orthogonal regression between the OMS with site-fitted MOST functions and OMS-footprint-weighted, direct EC flux estimates (ECwfp) and indirect EC flux estimates through CX2 ( ECCX2,wfp). The “m” indicates the regression slope.

Citation: Journal of Hydrometeorology 23, 6; 10.1175/JHM-D-21-0100.1

5. Conclusions

This paper studies the performance of an optical microwave scintillometer (OMS) in measuring evaporation over a saline lake in a heterogeneous setting in the Atacama Desert. Our research question, “How accurate can the measurements of an OMS over an open-water body be under arid conditions?”, was inspired by the possibility of the OMS to be installed at the borders of an open-water body. We conclude that the OMS is capable to measure open-water evaporation without the need to install instrumentation in the water. However, care must be taken to not include (unwillingly) surrounding areas with contrasting flux regimes in the OMS footprint. We argue that with measurements below the blending height an enhanced level of temperature and humidity fluctuations is observed. This is due to the passage of the scintillometer beam though the internal boundary layer that defines the transition between the different surfaces. As a result, Monin–Obukhov similarity theory (MOST) relations, needed to infer the sensible and latent heat fluxes from the measured structure parameters of temperature and humidity, can potentially fail. The OMS dependence on the MOST, which relies on the assumption of locally driven, undisturbed turbulence is a weakness that requires further study. We should also note that there also more trivial, instrumental issues that can cause MOST to fail, in the case of the OMS most notably are absorption fluctuations and tower or tripod vibrations.

For our study focusing on a saline lake in the Salar del Huasco, the OMS deployment proved to be challenging due to the unstable ground near the water edge. We therefore extended the path such that the instrument transmitter and receiver were installed on stable ground with a waterlogged salt surface. The resulting scintillometer setup was defined by an installation height of 2.3 m, which can be considered below the blending height and a pathlength of 1070 m covering a mixed footprint of open water and a waterlogged salt surface with a daytime footprint coverage of ∼80% and ∼20%, respectively. As a reference, we installed an EC system over the open water and one over the wet-salt surface. Even though both surfaces in the footprint were water saturated, they were strongly heterogeneous in their flux regimes due to their different evaporation properties. In short, the amount of available energy was much less over the wet-salt surface compared to the water surface (maximum Rn values of 900 W m−2 over water versus 500 W m−2 over the wet salt) due to the larger albedo and lower surface temperature. The salt also impeded evaporation leading to H > LυE in contrast to the water surface where LυE > H. The large contrast in fluxes was especially visible in LυE values (maximum values of 75 W m−2 over wet salt versus 550 W m−2 over water). For H values, the contrast was also notable but much less pronounced (maximum values of 200 W−2 over wet salt versus 100 W m−2 over water). In addition to the local heterogeneity, the fluxes were governed by a regional wind-driven turbulence regime with a low-wind, low-flux regime during the night and morning and strong-wind, high-flux regime during the afternoon.

The mixed footprint setup including a transition between the heterogeneous surfaces marked by disturbed turbulence made that our evaluation of the OMS fluxes against those measured by EC focused on the following three points. In our evaluation of the OMS fluxes against those measured by EC we addressed the heterogenous setting of the OMS measurements by critically evaluating the MOST functions, estimating an OMS footprint composing EC flux and discussing the regional-scale-driven turbulence regime. Concerning these three issues our main findings are as follows.

First, we evaluated the MOST functions used to infer the OMS measured structure parameters to fluxes. We did this for both the scintillometer as well as the EC data-derived structure parameters for temperature CT2 and humidity Cq2. Here, we found that for CT2, our data correspond well to the standard MOST functions as determined by Kooijmans and Hartogensis (2016), which was based on many experiments. For Cq2, however, we observe much higher OMS values than those based on the EC data. We attribute this to the aforementioned internal boundary layer disturbances at the transition between the surfaces. The OMS CT2 values are also somewhat higher than those from the EC data, but the difference is limited as the sensible heat flux contrast between the two surfaces is much less than that for Cq2. Again, instrument issues could be an explanation as well. We were, however, not able to fully discard this possibility as the recorded raw data stream missed 70% of the samples making a spectral analysis impracticable. To compensate for the elevated OMS Cq2 values, we adopted a practical approach by using site-adapted MOST coefficients. We should note that our approach to adjust Cq2 MOST functions to account for disturbances at the transition of the two heterogeneous surfaces assumes that the instrument was otherwise performing well. We applied a strict filtering the available data. In the end, we only used one day of data as the other days showed unphysically high scintillation intensities, likely due to a minor misalignment that might have caused surface reflections of the beam. Also, we strictly high-pass filter the raw 1-kHz data to exclude potential absorption fluctuations. Here, we adopted filter constants depending on the wind regimes.

Second, following the approach by Meijninger et al. (2002) for scintillometer measurements performed below the blending height, we validated the OMS against a composed flux consisting of OMS-footprint-weighted EC fluxes over the two surfaces. In addition, we considered the nonlinear relationship between structure parameters and surface fluxes in the footprint analysis that lead to an increase in EC composed LυE of 15%. All things combined, we reach a zero-intercept linear regression with an R2 of 0.92 and a slope of 0.95 when comparing the evaporation of the OMS against the OMS footprint composed EC estimates.

Third, we also discussed the influence of the regionally driven wind regime that controls the boundary layer development and turbulent mixing and is responsible for the advection and entrainment of dry air to the study area. These aspects influence both the EC and OMS measurements. Judging by the differences in structure parameters derived from the OMS and EC and how these influence the dimensionless MOST relations, it seems that the local MOST functions, in either the water or the wet salt, are not really affected by these nonlocal disturbances.

Finally, the weakness of the OMS method is its dependence on the Monin–Obukhov similarity theory, which relies on the assumption of locally driven fluxes. This assumption is often violated in (semi)arid regions where strong spatial contrasts in the temperature and humidity lead to temperature and humidity fluctuations that are driven by nonlocal processes such as regional circulations and internal boundary layers. Our approach was to introduce site-fitted MOST functions. We do not see this as a general solution but merely as a way to quantify the impact of the MOST violation. This is a challenge that we must overcome for the OMS to be used as a standard method for evaporation over open-water bodies.

Acknowledgments.

This research received financial support from the Chilean National Agency of Science and Development through the project ANID/FONDECYT/1210221. Support for Felipe Lobos-Roco was provided by the Wageningen University Ph.D. Sandwich Project 5160957644. F. Suárez acknowledges the support from the Centro de Desarrollo Urbano Sustentable (CEDEUS - ANID/FONDAP/15110020) and from the Centro de Excelencia en Geotermia de los Andes (CEGA - ANID/FONDAP/15200001).

Data availability statement.

Data available at http://dx.doi.org/10.17632/fkdzr7b5s6.1.

REFERENCES

  • Abtew, W., and A. Melesse, 2012: Evaporation and Evapotranspiration: Measurements and Estimations. Springer, 219 pp.

  • Aguirre, F., O. Hartogensis, F. Meza, and F. Suárez, 2022: Refinements and analysis of the optical-microwave scintillometry method applied to measurements over a vineyard in Chile. Water, 14, 474, https://doi.org/10.3390/w14030474.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Beyrich, F., and Coauthors, 2012: Towards a validation of scintillometer measurements: The LITFASS-2009 experiment. Bound.-Layer Meteor., 144, 83112, https://doi.org/10.1007/s10546-012-9715-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blin, N., M. Hausner, S. Leray, C. Lowry, and F. Suárez, 2021: Potential impacts of climate change on an aquifer in the arid Altiplano, northern Chile: The case of the protected wetlands of the Salar del Huasco basin. J. Hydrol. Reg. Stud., 39, 100996, https://doi.org/10.1016/j.ejrh.2022.100996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Businger, J. A., J. Wyngaard, Y. Izumi, and E. F. Bradley, 1971: Flux-profile relationships in the atmospheric surface layer. J. Atmos. Sci., 28, 181189, https://doi.org/10.1175/1520-0469(1971)028<0181:FPRITA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Charnock, H., 1955: Wind stress on a water surface. Quart. J. Roy. Meteor. Soc., 81, 639640, https://doi.org/10.1002/qj.49708135027.

  • de la Fuente, A., C. Meruane, and F. Suarez, 2021: Long-term spatiotemporal variability in high Andean wetlands in northern Chile. Sci. Total Environ., 756, 143830, https://doi.org/10.1016/j.scitotenv.2020.143830.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dyer, A., 1967: The turbulent transport of heat and water vapour in an unstable atmosphere. Quart. J. Roy. Meteor. Soc., 93, 501508, https://doi.org/10.1002/qj.49709339809.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Evans, J., and H. De Bruin, 2011: The effective height of a two-wavelength scintillometer system. Bound.-Layer Meteor., 141, 165177, https://doi.org/10.1007/s10546-011-9634-0.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Evans, J., D. McNeil, J. Finch, T. Murray, R. Harding, H. Ward, and A. Verhoef, 2012: Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain. Agric. For. Meteor., 166, 221233, https://doi.org/10.1016/j.agrformet.2012.07.010.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fratini, G., and M. Mauder, 2014: Towards a consistent eddy-covariance processing: An intercomparison of EddyPro and TK3. Atmos. Meas. Tech., 7, 22732281, https://doi.org/10.5194/amt-7-2273-2014.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Green, A., M. Astill, K. McAneney, and J. Nieveen, 2001: Path-averaged surface fluxes determined from infrared and microwave scintillometers. Agric. For. Meteor., 109, 233247, https://doi.org/10.1016/S0168-1923(01)00262-3.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Han, P.-F., X.-S. Wang, and J.-Z. Wang, 2019: Using large-aperture scintillometer to estimate lake-water evaporation and heat fluxes in the Badain Jaran Desert, China. Water, 11, 2575, https://doi.org/10.3390/w11122575.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hartogensis, O. K., 2006: Exploring scintillometry in the stable atmospheric surface layer. Ph.D. thesis, Wageningen University, 240 pp.

    • Search Google Scholar
    • Export Citation
  • Hartogensis, O. K., C. Watts, J. Rodriguez, and H. De Bruin, 2003: Derivation of an effective height for scintillometers: La Poza experiment in Northwest Mexico. J. Hydrometeor., 4, 915928, https://doi.org/10.1175/1525-7541(2003)004<0915:DOAEHF>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hill, R., 1997: Algorithms for obtaining atmospheric surface-layer fluxes from scintillation measurements. J. Atmos. Oceanic Technol., 14, 456467, https://doi.org/10.1175/1520-0426(1997)014<0456:AFOASL>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hill, R., S. Clifford, and R. S. Lawrence, 1980: Refractive-index and absorption fluctuations in the infrared caused by temperature, humidity, and pressure fluctuations. J. Opt. Soc. Amer., 70, 11921205, https://doi.org/10.1364/JOSA.70.001192.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Horst, T., and J. Weil, 1992: Footprint estimation for scalar flux measurements in the atmospheric surface layer. Bound.-Layer Meteor., 59, 279296, https://doi.org/10.1007/BF00119817.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Isabelle, P.-E., D. F. Nadeau, A. O. Perelet, E. R. Pardyjak, A. N. Rousseau, and F. Anctil, 2020: Application and evaluation of a two-wavelength scintillometry system for operation in a complex shallow boreal-forested valley. Bound.-Layer Meteor., 174, 341370, https://doi.org/10.1007/s10546-019-00488-7.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kohsiek, W., and M. Herben, 1983: Evaporation derived from optical and radio-wave scintillation. Appl. Opt., 22, 25662570, https://doi.org/10.1364/AO.22.002566.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kolmogorov, A. N., 1941: The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. C. R. Acad. Sci. URSS, 30, 301305.

    • Search Google Scholar
    • Export Citation
  • Kooijmans, L. M., and O. K. Hartogensis, 2016: Surface-layer similarity functions for dissipation rate and structure parameters of temperature and humidity based on eleven field experiments. Bound.-Layer Meteor., 160, 501527, https://doi.org/10.1007/s10546-016-0152-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lagouarde, J.-P., J.-M. Bonnefond, Y. Kerr, K. McAneney, and M. Irvine, 2002: Integrated sensible heat flux measurements of a two-surface composite landscape using scintillometry. Bound.-Layer Meteor., 105, 535, https://doi.org/10.1023/A:1019631428921.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leclerc, M. Y., and T. Foken, 2014: Footprints in Micrometeorology and Ecology. Springer, 258 pp.

  • Li, D., E. Bou-Zeid, and H. A. De Bruin, 2012: Monin–Obukhov similarity functions for the structure parameters of temperature and humidity. Bound.-Layer Meteor., 145, 4567, https://doi.org/10.1007/s10546-011-9660-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lobos-Roco, F., O. Hartogensis, J. Vilà-Guerau de Arellano, A. de la Fuente, R. Muñoz, J. Rutllant, and F. Suárez, 2021: Local evaporation controlled by regional atmospheric circulation in the Altiplano of the Atacama Desert. Atmos. Chem. Phys., 21, 91259150, https://doi.org/10.5194/acp-21-9125-2021.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lüdi, A., F. Beyrich, and C. Matzler, 2005: Determination of the turbulent temperature–humidity correlation from scintillometric measurements. Bound.-Layer Meteor., 117, 525550, https://doi.org/10.1007/s10546-005-1751-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, J. Sun, N. O. Jensen, H. Jørgensen, E. Pardyjak, and H. Fernando, 2001: Determination of the surface drag coefficient. Bound.-Layer Meteor., 99, 249276, https://doi.org/10.1023/A:1018915228170.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Massman, W., 2000: A simple method for estimating frequency response corrections for eddy covariance systems. Agric. For. Meteor., 104, 185198, https://doi.org/10.1016/S0168-1923(00)00164-7.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meijninger, W., O. Hartogensis, W. Kohsiek, J. Hoedjes, R. Zuurbier, and H. De Bruin, 2002: Determination of area-averaged sensible heat fluxes with a large aperture scintillometer over a heterogeneous surface–Flevoland field experiment. Bound.-Layer Meteor., 105, 3762, https://doi.org/10.1023/A:1019647732027.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meijninger, W. M. L., F. Beyrich, A. Ludi, W. Kohsiek, and H. D. Bruin, 2006: Scintillometer-based turbulent fluxes of sensible and latent heat over a heterogeneous land surface–A contribution to LITFASS-2003. Bound.-Layer Meteor., 121, 89110, https://doi.org/10.1007/s10546-005-9022-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Monin, A. S., and A. M. Obukhov, 1954: Basic laws of turbulent mixing in the surface layer of the atmosphere. Tr. Geofiz. Inst., Akad. Nauk SSSR, 24, 163187.

    • Search Google Scholar
    • Export Citation
  • Munoz, R. C., M. J. Falvey, M. Arancibia, V. I. Astudillo, J. Elgueta, M. Ibarra, C. Santana, and C. Vasquez, 2018: Wind energy exploration over the Atacama Desert: A numerical model–guided observational program. Bull. Amer. Meteor. Soc., 99, 20792092, https://doi.org/10.1175/BAMS-D-17-0019.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Neftel, A., C. Spirig, and C. Ammann, 2008: Application and test of a simple tool for operational footprint evaluations. Environ. Pollut., 152, 644652, https://doi.org/10.1016/j.envpol.2007.06.062.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nordbo, A., S. Launiainen, I. Mammarella, M. Lepparanta, J. Huotari, A. Ojala, and T. Vesala, 2011: Long-term energy flux measurements and energy balance over a small boreal lake using eddy covariance technique. J. Geophys. Res., 116, D02119, https://doi.org/10.1029/2010JD014542.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Oroud, I., 1999: Temperature and evaporation dynamics of saline solutions. J. Hydrol., 226, 110, https://doi.org/10.1016/S0022-1694(99)00138-9.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rosen, M. R., 1994: The importance of groundwater in playas: A review of playa classifications and the sedimentology and hydrology of playas. Paleoclimate and Basin Evolution of Playa Systems, GSA Special Papers, Vol. 289, Geological Society of America, https://doi.org/10.1130/SPE289-p1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rutllant, J. A., H. Fuenzalida, and P. Aceituno, 2003: Climate dynamics along the arid northern coast of Chile: The 1997–1998 Dinaica del Clima de la region de Antofagasta (DICLIMA) experiment. J. Geophys. Res., 108, 4538, https://doi.org/10.1029/2002JD003357.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Salhotra, A. M., E. E. Adams, and D. R. Harleman, 1985: Effect of salinity and ionic composition on evaporation: Analysis of Dead Sea evaporation pans. Water Resour. Res., 21, 13361344, https://doi.org/10.1029/WR021i009p01336.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Schuepp, P., M. Leclerc, J. MacPherson, and R. Desjardins, 1990: Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation. Bound.-Layer Meteor., 50, 355373, https://doi.org/10.1007/BF00120530.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shuttleworth, W. J., 2008: Evapotranspiration measurement methods. Southwest Hydrol., 7, 2223.

  • Stoffer, R., 2018: Revisiting raw data processing of combined optical-microwave scintillometers. M.S. thesis, Wageningen University, 59 pp.

    • Search Google Scholar
    • Export Citation
  • Suarez, F., F. Lobos-Roco, A. De La Fuente, J. Vila-Guerau de Arellano, A. Prieto, C. Meruane, and O. Hartogensis, 2020: E-data: A comprehensive field campaign to investigate evaporation enhanced by advection in the hyper-arid Altiplano. Water, 12, 745, https://doi.org/10.3390/w12030745.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tanny, J., S. Cohen, S. Assouline, F. Lange, A. Grava, D. Berger, B. Teltch, and M. Parlange, 2008: Evaporation from a small water reservoir: Direct measurements and estimates. J. Hydrol., 351, 218229, https://doi.org/10.1016/j.jhydrol.2007.12.012.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tatarski, V. I., 1961: Wave Propagation in a Turbulent Medium. McGraw-Hill, 285 pp.

  • Uribe, J., J. F. Munoz, J. Gironas, R. Oyarzun, E. Aguirre, and R. Aravena, 2015: Assessing groundwater recharge in an Andean closed basin using isotopic characterization and a rainfall runoff model: Salar del Huasco basin, Chile. Hydrogeol. J., 23, 15351551, https://doi.org/10.1007/s10040-015-1300-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Vickers, D., and L. Mahrt, 1997: Quality control and flux sampling problems for tower and aircraft data. J. Atmos. Oceanic Technol., 14, 512526, https://doi.org/10.1175/1520-0426(1997)014<0512:QCAFSP>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, T.-i., G. Ochs, and S. Clifford, 1978: A saturation-resistant optical scintillometer to measure CN2. J. Opt. Soc. Amer., 68, 334338, https://doi.org/10.1364/JOSA.68.000334.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ward, H., J. Evans, O. Hartogensis, A. Moene, H. De Bruin, and C. Grimmond, 2013: A critical revision of the estimation of the latent heat flux from two-wavelength scintillometry. Quart. J. Roy. Meteor. Soc., 139, 19121922, https://doi.org/10.1002/qj.2076.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ward, H., J. G. Evans, C. S. B. Grimmond, and J. Bradford, 2015: Infrared and millimetre-wave scintillometry in the suburban environment–Part 1: Structure parameters. Atmos. Meas. Tech., 8, 13851405, https://doi.org/10.5194/amt-8-1385-2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wieringa, J., 1976: An objective exposure correction method for average wind speeds measured at a sheltered location. Quart. J. Roy. Meteor. Soc., 102, 241253, https://doi.org/10.1002/qj.49710243119.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wilczack, J., S. Oncley, and S. Stage, 2001: Sonic anemometers tilt corrections algorithms. Bound.-Layer Meteor., 99, 127150, https://doi.org/10.1023/A:1018966204465.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wyngaard, J., Y. Izumi, and S. A. Collins, 1971: Behavior of the refractive-index-structure parameter near the ground. J. Opt. Soc. Amer., 61, 16461650, https://doi.org/10.1364/JOSA.61.001646.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhang, G., J. Zhang, and P. Meng, 2021: Estimation of kilometer-scale heat fluxes over a hilly area in northern China using an optical-microwave scintillometer. Agric. Water Manage., 244, 106582, https://doi.org/10.1016/j.agwat.2020.106582.

    • Crossref
    • Search Google Scholar
    • Export Citation
Save
  • Abtew, W., and A. Melesse, 2012: Evaporation and Evapotranspiration: Measurements and Estimations. Springer, 219 pp.

  • Aguirre, F., O. Hartogensis, F. Meza, and F. Suárez, 2022: Refinements and analysis of the optical-microwave scintillometry method applied to measurements over a vineyard in Chile. Water, 14, 474, https://doi.org/10.3390/w14030474.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Beyrich, F., and Coauthors, 2012: Towards a validation of scintillometer measurements: The LITFASS-2009 experiment. Bound.-Layer Meteor., 144, 83112, https://doi.org/10.1007/s10546-012-9715-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blin, N., M. Hausner, S. Leray, C. Lowry, and F. Suárez, 2021: Potential impacts of climate change on an aquifer in the arid Altiplano, northern Chile: The case of the protected wetlands of the Salar del Huasco basin. J. Hydrol. Reg. Stud., 39, 100996, https://doi.org/10.1016/j.ejrh.2022.100996.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Businger, J. A., J. Wyngaard, Y. Izumi, and E. F. Bradley, 1971: Flux-profile relationships in the atmospheric surface layer. J. Atmos. Sci., 28, 181189, https://doi.org/10.1175/1520-0469(1971)028<0181:FPRITA>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Charnock, H., 1955: Wind stress on a water surface. Quart. J. Roy. Meteor. Soc., 81, 639640, https://doi.org/10.1002/qj.49708135027.

  • de la Fuente, A., C. Meruane, and F. Suarez, 2021: Long-term spatiotemporal variability in high Andean wetlands in northern Chile. Sci. Total Environ., 756, 143830, https://doi.org/10.1016/j.scitotenv.2020.143830.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Dyer, A., 1967: The turbulent transport of heat and water vapour in an unstable atmosphere. Quart. J. Roy. Meteor. Soc., 93, 501508, https://doi.org/10.1002/qj.49709339809.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Evans, J., and H. De Bruin, 2011: The effective height of a two-wavelength scintillometer system. Bound.-Layer Meteor., 141, 165177, https://doi.org/10.1007/s10546-011-9634-0.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Evans, J., D. McNeil, J. Finch, T. Murray, R. Harding, H. Ward, and A. Verhoef, 2012: Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain. Agric. For. Meteor., 166, 221233, https://doi.org/10.1016/j.agrformet.2012.07.010.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fratini, G., and M. Mauder, 2014: Towards a consistent eddy-covariance processing: An intercomparison of EddyPro and TK3. Atmos. Meas. Tech., 7, 22732281, https://doi.org/10.5194/amt-7-2273-2014.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Green, A., M. Astill, K. McAneney, and J. Nieveen, 2001: Path-averaged surface fluxes determined from infrared and microwave scintillometers. Agric. For. Meteor., 109, 233247, https://doi.org/10.1016/S0168-1923(01)00262-3.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Han, P.-F., X.-S. Wang, and J.-Z. Wang, 2019: Using large-aperture scintillometer to estimate lake-water evaporation and heat fluxes in the Badain Jaran Desert, China. Water, 11, 2575, https://doi.org/10.3390/w11122575.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hartogensis, O. K., 2006: Exploring scintillometry in the stable atmospheric surface layer. Ph.D. thesis, Wageningen University, 240 pp.

    • Search Google Scholar
    • Export Citation
  • Hartogensis, O. K., C. Watts, J. Rodriguez, and H. De Bruin, 2003: Derivation of an effective height for scintillometers: La Poza experiment in Northwest Mexico. J. Hydrometeor., 4, 915928, https://doi.org/10.1175/1525-7541(2003)004<0915:DOAEHF>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hill, R., 1997: Algorithms for obtaining atmospheric surface-layer fluxes from scintillation measurements. J. Atmos. Oceanic Technol., 14, 456467, https://doi.org/10.1175/1520-0426(1997)014<0456:AFOASL>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hill, R., S. Clifford, and R. S. Lawrence, 1980: Refractive-index and absorption fluctuations in the infrared caused by temperature, humidity, and pressure fluctuations. J. Opt. Soc. Amer., 70, 11921205, https://doi.org/10.1364/JOSA.70.001192.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Horst, T., and J. Weil, 1992: Footprint estimation for scalar flux measurements in the atmospheric surface layer. Bound.-Layer Meteor., 59, 279296, https://doi.org/10.1007/BF00119817.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Isabelle, P.-E., D. F. Nadeau, A. O. Perelet, E. R. Pardyjak, A. N. Rousseau, and F. Anctil, 2020: Application and evaluation of a two-wavelength scintillometry system for operation in a complex shallow boreal-forested valley. Bound.-Layer Meteor., 174, 341370, https://doi.org/10.1007/s10546-019-00488-7.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kohsiek, W., and M. Herben, 1983: Evaporation derived from optical and radio-wave scintillation. Appl. Opt., 22, 25662570, https://doi.org/10.1364/AO.22.002566.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kolmogorov, A. N., 1941: The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. C. R. Acad. Sci. URSS, 30, 301305.

    • Search Google Scholar
    • Export Citation
  • Kooijmans, L. M., and O. K. Hartogensis, 2016: Surface-layer similarity functions for dissipation rate and structure parameters of temperature and humidity based on eleven field experiments. Bound.-Layer Meteor., 160, 501527, https://doi.org/10.1007/s10546-016-0152-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lagouarde, J.-P., J.-M. Bonnefond, Y. Kerr, K. McAneney, and M. Irvine, 2002: Integrated sensible heat flux measurements of a two-surface composite landscape using scintillometry. Bound.-Layer Meteor., 105, 535, https://doi.org/10.1023/A:1019631428921.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leclerc, M. Y., and T. Foken, 2014: Footprints in Micrometeorology and Ecology. Springer, 258 pp.

  • Li, D., E. Bou-Zeid, and H. A. De Bruin, 2012: Monin–Obukhov similarity functions for the structure parameters of temperature and humidity. Bound.-Layer Meteor., 145, 4567, https://doi.org/10.1007/s10546-011-9660-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lobos-Roco, F., O. Hartogensis, J. Vilà-Guerau de Arellano, A. de la Fuente, R. Muñoz, J. Rutllant, and F. Suárez, 2021: Local evaporation controlled by regional atmospheric circulation in the Altiplano of the Atacama Desert. Atmos. Chem. Phys., 21, 91259150, https://doi.org/10.5194/acp-21-9125-2021.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lüdi, A., F. Beyrich, and C. Matzler, 2005: Determination of the turbulent temperature–humidity correlation from scintillometric measurements. Bound.-Layer Meteor., 117, 525550, https://doi.org/10.1007/s10546-005-1751-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, J. Sun, N. O. Jensen, H. Jørgensen, E. Pardyjak, and H. Fernando, 2001: Determination of the surface drag coefficient. Bound.-Layer Meteor., 99, 249276, https://doi.org/10.1023/A:1018915228170.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Massman, W., 2000: A simple method for estimating frequency response corrections for eddy covariance systems. Agric. For. Meteor., 104, 185198, https://doi.org/10.1016/S0168-1923(00)00164-7.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meijninger, W., O. Hartogensis, W. Kohsiek, J. Hoedjes, R. Zuurbier, and H. De Bruin, 2002: Determination of area-averaged sensible heat fluxes with a large aperture scintillometer over a heterogeneous surface–Flevoland field experiment. Bound.-Layer Meteor., 105, 3762, https://doi.org/10.1023/A:1019647732027.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meijninger, W. M. L., F. Beyrich, A. Ludi, W. Kohsiek, and H. D. Bruin, 2006: Scintillometer-based turbulent fluxes of sensible and latent heat over a heterogeneous land surface–A contribution to LITFASS-2003. Bound.-Layer Meteor., 121, 89110, https://doi.org/10.1007/s10546-005-9022-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Monin, A. S., and A. M. Obukhov, 1954: Basic laws of turbulent mixing in the surface layer of the atmosphere. Tr. Geofiz. Inst., Akad. Nauk SSSR, 24, 163187.

    • Search Google Scholar
    • Export Citation
  • Munoz, R. C., M. J. Falvey, M. Arancibia, V. I. Astudillo, J. Elgueta, M. Ibarra, C. Santana, and C. Vasquez, 2018: Wind energy exploration over the Atacama Desert: A numerical model–guided observational program. Bull. Amer. Meteor. Soc., 99, 20792092, https://doi.org/10.1175/BAMS-D-17-0019.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Neftel, A., C. Spirig, and C. Ammann, 2008: Application and test of a simple tool for operational footprint evaluations. Environ. Pollut., 152, 644652, https://doi.org/10.1016/j.envpol.2007.06.062.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nordbo, A., S. Launiainen, I. Mammarella, M. Lepparanta, J. Huotari, A. Ojala, and T. Vesala, 2011: Long-term energy flux measurements and energy balance over a small boreal lake using eddy covariance technique. J. Geophys. Res., 116, D02119, https://doi.org/10.1029/2010JD014542.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Oroud, I., 1999: Temperature and evaporation dynamics of saline solutions. J. Hydrol., 226, 110, https://doi.org/10.1016/S0022-1694(99)00138-9.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rosen, M. R., 1994: The importance of groundwater in playas: A review of playa classifications and the sedimentology and hydrology of playas. Paleoclimate and Basin Evolution of Playa Systems, GSA Special Papers, Vol. 289, Geological Society of America, https://doi.org/10.1130/SPE289-p1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Rutllant, J. A., H. Fuenzalida, and P. Aceituno, 2003: Climate dynamics along the arid northern coast of Chile: The 1997–1998 Dinaica del Clima de la region de Antofagasta (DICLIMA) experiment. J. Geophys. Res., 108, 4538, https://doi.org/10.1029/2002JD003357.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Salhotra, A. M., E. E. Adams, and D. R. Harleman, 1985: Effect of salinity and ionic composition on evaporation: Analysis of Dead Sea evaporation pans. Water Resour. Res., 21, 13361344, https://doi.org/10.1029/WR021i009p01336.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Schuepp, P., M. Leclerc, J. MacPherson, and R. Desjardins, 1990: Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation. Bound.-Layer Meteor., 50, 355373, https://doi.org/10.1007/BF00120530.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shuttleworth, W. J., 2008: Evapotranspiration measurement methods. Southwest Hydrol., 7, 2223.

  • Stoffer, R., 2018: Revisiting raw data processing of combined optical-microwave scintillometers. M.S. thesis, Wageningen University, 59 pp.

    • Search Google Scholar
    • Export Citation
  • Suarez, F., F. Lobos-Roco, A. De La Fuente, J. Vila-Guerau de Arellano, A. Prieto, C. Meruane, and O. Hartogensis, 2020: E-data: A comprehensive field campaign to investigate evaporation enhanced by advection in the hyper-arid Altiplano. Water, 12, 745, https://doi.org/10.3390/w12030745.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tanny, J., S. Cohen, S. Assouline, F. Lange, A. Grava, D. Berger, B. Teltch, and M. Parlange, 2008: Evaporation from a small water reservoir: Direct measurements and estimates. J. Hydrol., 351, 218229, https://doi.org/10.1016/j.jhydrol.2007.12.012.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tatarski, V. I., 1961: Wave Propagation in a Turbulent Medium. McGraw-Hill, 285 pp.

  • Uribe, J., J. F. Munoz, J. Gironas, R. Oyarzun, E. Aguirre, and R. Aravena, 2015: Assessing groundwater recharge in an Andean closed basin using isotopic characterization and a rainfall runoff model: Salar del Huasco basin, Chile. Hydrogeol. J., 23, 15351551, https://doi.org/10.1007/s10040-015-1300-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Vickers, D., and L. Mahrt, 1997: Quality control and flux sampling problems for tower and aircraft data. J. Atmos. Oceanic Technol., 14, 512526, https://doi.org/10.1175/1520-0426(1997)014<0512:QCAFSP>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, T.-i., G. Ochs, and S. Clifford, 1978: A saturation-resistant optical scintillometer to measure CN2. J. Opt. Soc. Amer., 68, 334338, https://doi.org/10.1364/JOSA.68.000334.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ward, H., J. Evans, O. Hartogensis, A. Moene, H. De Bruin, and C. Grimmond, 2013: A critical revision of the estimation of the latent heat flux from two-wavelength scintillometry. Quart. J. Roy. Meteor. Soc., 139, 19121922, https://doi.org/10.1002/qj.2076.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ward, H., J. G. Evans, C. S. B. Grimmond, and J. Bradford, 2015: Infrared and millimetre-wave scintillometry in the suburban environment–Part 1: Structure parameters. Atmos. Meas. Tech., 8, 13851405, https://doi.org/10.5194/amt-8-1385-2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wieringa, J., 1976: An objective exposure correction method for average wind speeds measured at a sheltered location. Quart. J. Roy. Meteor. Soc., 102, 241253, https://doi.org/10.1002/qj.49710243119.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wilczack, J., S. Oncley, and S. Stage, 2001: Sonic anemometers tilt corrections algorithms. Bound.-Layer Meteor., 99, 127150, https://doi.org/10.1023/A:1018966204465.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wyngaard, J., Y. Izumi, and S. A. Collins, 1971: Behavior of the refractive-index-structure parameter near the ground. J. Opt. Soc. Amer., 61, 16461650, https://doi.org/10.1364/JOSA.61.001646.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhang, G., J. Zhang, and P. Meng, 2021: Estimation of kilometer-scale heat fluxes over a hilly area in northern China using an optical-microwave scintillometer. Agric. Water Manage., 244, 106582, https://doi.org/10.1016/j.agwat.2020.106582.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Flowchart of the steps used to calculate the sensible (H) and latent heat (LυE) fluxes using the (a) optical microwave scintillation method and (b) eddy covariance method. The large arrows depict how the turbulent eddies are advected from their source area (footprint) to the scintillometer and EC sensor paths. Roman numerals I, II, and III represent the OMS calculations steps described in section 2b. A picture of every instrument is depicted next to the schemes.

  • Fig. 2.

    Study area. (a) Salar del Huasco basin and saline lake. (b) Installation of the OMS and EC systems. Orange circles represent the OMS receiver (RX) and transmitter (TX), solid black line the OMS path and yellow circles the EC systems. Windrose indicates wind speed and direction of 15 Nov 2018.

  • Fig. 3.

    Diurnal relationship between wind speed (U) and wind direction (WD) obtained by EC over water and wet-salt surfaces.

  • Fig. 4.

    Atmospheric boundary layer conditions measured during E-DATA on 19 Nov 2018. (top) Diurnal cycle of the potential temperature radiosounding profiles. (bottom) Diurnal cycle of specific humidity radiosounding profiles.

  • Fig. 5.

    Diurnal-averaged time series of near-surface (∼5 cm) temperature, net radiation, and turbulent fluxes over the (a) wet-salt and (b) water surfaces during the E-DATA period. The picture behind exemplifies the wet-salt and water surface during the experiment. Shades indicate the minimum and maximum values across the E-DATA period.

  • Fig. 6.

    Footprints of the (a),(c) LAS and EC and (b),(d) MWS and EC over the water and wet-salt surfaces for representative hours of wind regimes on 15 Nov 2018. Colors (given without scale) represent the percentage contribution of that area to the total footprint. The wind roses characterize the wind regimes for the same day (top) between 0700 and 1200 LT and (middle) between 1200 and 2000 LT. The background image is a georeferenced mosaic based on a drone flight on 18 Nov 2018. (e) Diurnal cycle of LAS and MWS percentage footprint contributions from the water (solid lines) and wet-salt surface (dotted lines).

  • Fig. 7.

    Relationship between roughness length (z0) and wind speed over (a) water surface and (b) wet-salt surface. Wind speed data are taken from the EC systems. The solid lines express the values to which z0 converges for high wind speeds.

  • Fig. 8.

    Dimensionless CT2 and Cq2 groups vs stability parameter (z/Lob). (a),(b) CT2 and Cq2 dimensionless group using parameters from EC water and wet salt weighted by the OMS footprint (wfp) during the E-DATA period. (c),(d) CT2 and Cq2 dimensionless group using structure parameters measured by the OMS and the remaining parameters from EC water and wet salt weighted by the OMS footprint (wfp) during 15 Nov 2018. MOST functions proposed by Kooijmans and Hartogensis (2016) and the proposed ones in this study are shown in solid and dashed lines, respectively. The experimental data points are marked by their wind regime (windy > 4 m s−1 between 1300 and 2040 LT and calm < 4 m s−1 between 2340 and 1300 LT).

  • Fig. 9.

    Diurnal cycle of (a),(c) surface fluxes and (b),(d) orthogonal regression between EC over water and wet salt and an OMS footprint composed EC flux (ECwfp) for 15 Nov 2018. The “m” indicates the regression slope; gray data points refer to wet salt and black ones to water.

  • Fig. 10.

    Relationships of (a) latent and (b) sensible heat fluxes between OMS-footprint-weighted, direct EC flux estimates (ECwfp) and indirect EC flux estimates through CX2 ( ECCX2,wfp).

  • Fig. 11.

    (a) Diurnal cycle and (b) orthogonal regression of latent heat flux from OMS data calculated with standard and site-fitted MOST functions. m indicates the regression slope.

  • Fig. 12.

    Diurnal cycle of (a),(c) surface fluxes and (b),(d) orthogonal regression between the OMS with site-fitted MOST functions and OMS-footprint-weighted, direct EC flux estimates (ECwfp) and indirect EC flux estimates through CX2 ( ECCX2,wfp). The “m” indicates the regression slope.

All Time Past Year Past 30 Days
Abstract Views 372 0 0
Full Text Views 677 362 43
PDF Downloads 517 250 42