Determination of Evaporative Fluxes Using a Bench-Scale Atmosphere Simulator

: An accurate determination of evaporative ﬂuxes is critical for efﬁcient water management in semi-arid climates such as in the Canadian Prairies. The main achievements of this research are the design and operation of a bench-scale atmosphere simulator, performance evaluation using selected weather scenarios pertaining to regional atmospheric conditions, validation using established empirical correlations, and estimation of evaporation rates and the amount for a typical local water body. Results indicate that the measured data achieved the target values for the various parameters and the data were found to be stable during the 3-h test duration. The vapour ﬂux was found to have large variation during summer (0.120 g · s − 1 · m − 2 during the day and 0.047 g · s − 1 · m − 2 at night), low variation during spring (0.116 g · s − 1 · m − 2 during the day and 0.062 g · s − 1 · m − 2 at night), and negligible change during fall (0.100 g · s − 1 · m − 2 during the day and 0.076 g · s − 1 · m − 2 at night). The measured vapour ﬂux was generally within one standard deviation of the equality line when compared with that predicted by both the mass-transfer equations and the combination equations. The average evaporation ranged from 4 mm · d − 1 to 8 mm · d − 1 during the day and decreased to 1 mm · d − 1 to 3 mm · d − 1 at night. The 24-h evaporation was found to be 8 ± 1 mm · d − 1 from late April through late October. Likewise, the cumulative annual evaporation was found to be 1781 mm, of which 82% occurs during the day and 18% at night.


Introduction
Potential evaporation from open water bodies is a complex phenomenon because of the interactions between meteorological and physiological factors [1]. The net upward movement of vapour from an exposed water surface to the atmosphere depends on the interaction between atmospheric and surface parameters. The constitutive relationships are based on the laws of conservation for momentum, mass, and energy [2]. Several methods to determine evaporative fluxes from water surfaces are available in the literature [3][4][5]. The resulting empirical relationships for predicting potential evaporation require accurate field and/or laboratory measurements [6]. Field data is affected by spatial and temporal variations in atmospheric parameters and physiographic features and, as such, the resulting correlations are site-specific and assume constant climatic conditions [7]. For example, the variation in pan evaporation and estimated evaporation is reported to be governed by solar irradiation in an open lake in Australia [8]. Conversely, laboratory determination can create a simplified environment by isolating selected influencing parameters provided these are adequately replicated [9]. Most climate chambers only partially capture the various atmospheric components because of the difficulty arising from the interactions. Therefore, there is a need to understand the phenomenon of evaporation from a fundamental perspective.
The semi-arid Canadian Prairies have the highest water demand-to-availability ratio in Canada [10]. Low and spatiotemporally variable precipitation [11], a reliance on seasonally variable glacial runoff [12], interprovincial water use agreements [13], and competing Water 2021, 13, 84 2 of 14 municipal and industrial requirements [14] pose an acute risk to water security in the area. More than 70 storage facilities have been constructed to manage the regional water supply [15]. The large surface areas of these shallow water bodies facilitate evaporation [16]. In particular, the surface area of Lake Diefenbaker varies from 102-430 km 2 , such that a 10 mm elevation drop results in a water loss of 1-4.3 × 10 6 m 3 [17]. The annual evaporation in this man-made reservoir is estimated to be up to 1000 mm [18] thereby accounting for more than all other withdrawals combined [19]. Furthermore, the upcoming irrigation project proposes to extract 850 × 10 6 m 3 of water from the lake which is almost twice the estimated annual evaporation [20]. To ensure uninterrupted availability of water for all purposes, there is an exigent need to understand evaporation from Lake Diefenbaker.
The main objective of this paper is to determine evaporative fluxes using a de novo bench-scale atmosphere simulator (BAS). First, the design and operation of the climate chamber and the supporting modules are presented. Second, the performance of BAS was evaluated using selected weather scenarios pertaining to regional atmospheric conditions. Third, the measured and analysed data were validated using established empirical correlations. Fourth, potential evaporation rates and the amount for a typical local water body were estimated.

Design and Operation of the BAS
The BAS was designed to capture the main atmospheric parameters affecting evaporation particularly those operating in the region: (i) solar irradiance because the surface receives 220 W/m 2 for up to 16 h/d in the spring, 240 W/m 2 in the summer for 17 h/d, and 105 W/m 2 for 13 h/d in fall in Saskatchewan, and the area receives the highest annual radiant energy (154 W/m 2 in Regina) in Canada [21]; (ii) wind effects because of turbulent transport due to a flat area devoid of physical buffers such that wind speeds generally vary between 16 and 21 km/h (measured at 2 m above the surface) throughout the year [11], and; (iii) air temperature and humidity because of the warm and dry climate (Dfb according to Köppen-Geiger classification) and the significance of fetch effects [22]. Figure 1 gives a schematic of the BAS. The system layout comprises the climate chamber ( Figure 2) and supporting modules ( Figure 3) that were configured to govern the atmospheric parameters of air velocity, humidity, temperature, and solar irradiance. The BAS worked on the modular operation of the individual components. To create steadystate meteorological conditions during testing, the various modules were adjusted and the desired atmospheric parameters were achieved within a conditioning period of 0.5 to 1.5 h. The semi-arid Canadian Prairies have the highest water demand-to-availability ratio in Canada [10]. Low and spatiotemporally variable precipitation [11], a reliance on seasonally variable glacial runoff [12], interprovincial water use agreements [13], and competing municipal and industrial requirements [14] pose an acute risk to water security in the area. More than 70 storage facilities have been constructed to manage the regional water supply [15]. The large surface areas of these shallow water bodies facilitate evaporation [16]. In particular, the surface area of Lake Diefenbaker varies from 102-430 km 2 , such that a 10 mm elevation drop results in a water loss of 1-4.3 × 10 6 m 3 [17]. The annual evaporation in this man-made reservoir is estimated to be up to 1000 mm [18] thereby accounting for more than all other withdrawals combined [19]. Furthermore, the upcoming irrigation project proposes to extract 850 × 10 6 m 3 of water from the lake which is almost twice the estimated annual evaporation [20]. To ensure uninterrupted availability of water for all purposes, there is an exigent need to understand evaporation from Lake Diefenbaker.
The main objective of this paper is to determine evaporative fluxes using a de novo bench-scale atmosphere simulator (BAS). First, the design and operation of the climate chamber and the supporting modules are presented. Second, the performance of BAS was evaluated using selected weather scenarios pertaining to regional atmospheric conditions. Third, the measured and analysed data were validated using established empirical correlations. Fourth, potential evaporation rates and the amount for a typical local water body were estimated.

Design and Operation of the BAS
The BAS was designed to capture the main atmospheric parameters affecting evaporation particularly those operating in the region: (i) solar irradiance because the surface receives 220 W/m 2 for up to 16 h/d in the spring, 240 W/m 2 in the summer for 17 h/d, and 105 W/m 2 for 13 h/d in fall in Saskatchewan, and the area receives the highest annual radiant energy (154 W/m 2 in Regina) in Canada [21]; (ii) wind effects because of turbulent transport due to a flat area devoid of physical buffers such that wind speeds generally vary between 16 and 21 km/h (measured at 2 m above the surface) throughout the year [11], and; (iii) air temperature and humidity because of the warm and dry climate (Dfb according to Köppen-Geiger classification) and the significance of fetch effects [22]. Figure 1 gives a schematic of the BAS. The system layout comprises the climate chamber ( Figure 2) and supporting modules ( Figure 3) that were configured to govern the atmospheric parameters of air velocity, humidity, temperature, and solar irradiance. The BAS worked on the modular operation of the individual components. To create steadystate meteorological conditions during testing, the various modules were adjusted and the desired atmospheric parameters were achieved within a conditioning period of 0.5 to 1.5 h.   The climate chamber was designed to measure a number of parameters in the atmospheric and surface layers. The sensor dimensions, technical details, measuring principles, accuracy and precision are not given in this paper. The various sensors were arranged around the sample surface to preclude interferences and data were usually recorded at 10-s intervals. The air velocity was measured 65 mm behind the sample centre and 30 mm above the surface by a propeller-type anemometer (60 mm diameter). Air turbulence above the surface was reduced and assumed negligible by the incorporation of a flow straightener [23] and draft shield [24]. Air pressure in the chamber was recorded by a wallmounted barometer. Air temperature and relative humidity, recorded by thermometers and hygrometers, were measured at four points: low (25 mm above surface); high (75 mm above surface); upwind (50 mm before sample centre); and downwind (50 mm after sample centre). Values of absolute humidity are also reported to represent a measure of water in air irrespective of temperature. Sample mass was recorded by an analytical balance. Surface temperature was measured approximately 95 mm above the surface using a pyrometer. Incoming and outgoing shortwave radiant energy (285 nm to 3 μm), recorded before and after each experiment by a pyranometer and albedometer, were measured at surface level and about 125 mm above. The corresponding longwave radiant energy (8 to 14 μm) were estimated [25]. Atmospheric sensor data were displayed on digital readouts in realtime, thereby allowing parameter re-adjustments when required.
The air supply ( Figure 1) design was based on an open-air and non-return wind tunnel [26] with air sucked by a fan at a controlled flow rate [27]. The simple open-air and nonreturn design completely isolated the BAS from interfering with ambient humidity and temperature.
Likewise, the fan-driven suction precluded possible pumping effects [28]. An ultrasonic humidifier (Figure 3a) was used to increase vapour content in the air [29]. In contrast, the cooler/dehumidifier (Figure 3b) used a glycol chiller to decrease vapour content and remove thermal energy from the air [9]. Finally, the heater ( Figure 3c) used an electric element to add thermal energy to the air [30].   The solar irradiation module ( Figure 3d) used a 275 W quartz-tungsten-halogen (QTH) bulb (~225 mm above the surface) to mimic the irradiant shortwave energy and spectral pattern of the sun [31]. The externally mounted module ensured clear separation of shortwave and longwave bands thereby closely replicating actual conditions. The thermal heat generated by the QTH bulb was removed by a vacuum-fan suction system [32]. The power controller was adjusted to ensure that the sample received the desired shortwave irradiance, as measured by the pyranometer. The climate chamber was designed to measure a number of parameters in the atmospheric and surface layers. The sensor dimensions, technical details, measuring principles, accuracy and precision are not given in this paper. The various sensors were arranged around the sample surface to preclude interferences and data were usually recorded at 10-s intervals. The air velocity was measured 65 mm behind the sample centre and 30 mm above the surface by a propeller-type anemometer (60 mm diameter). Air turbulence above the surface was reduced and assumed negligible by the incorporation of a flow straightener [23] and draft shield [24]. Air pressure in the chamber was recorded by a wallmounted barometer. Air temperature and relative humidity, recorded by thermometers and hygrometers, were measured at four points: low (25 mm above surface); high (75 mm above surface); upwind (50 mm before sample centre); and downwind (50 mm after sample centre). Values of absolute humidity are also reported to represent a measure of water in air irrespective of temperature. Sample mass was recorded by an analytical balance. Surface temperature was measured approximately 95 mm above the surface using a pyrometer. Incoming and outgoing shortwave radiant energy (285 nm to 3 µm), recorded before and after each experiment by a pyranometer and albedometer, were measured at surface level and about 125 mm above. The corresponding longwave radiant energy (8 to 14 µm) were estimated [25]. Atmospheric sensor data were displayed on digital readouts in real-time, thereby allowing parameter re-adjustments when required.

Research Methodology
The air supply ( Figure 1) design was based on an open-air and non-return wind tunnel [26] with air sucked by a fan at a controlled flow rate [27]. The simple open-air and non-return design completely isolated the BAS from interfering with ambient humidity and temperature.
Likewise, the fan-driven suction precluded possible pumping effects [28]. An ultrasonic humidifier (Figure 3a) was used to increase vapour content in the air [29]. In contrast, the cooler/dehumidifier (Figure 3b) used a glycol chiller to decrease vapour content and remove thermal energy from the air [9]. Finally, the heater ( Figure 3c) used an electric element to add thermal energy to the air [30].
The solar irradiation module ( Figure 3d) used a 275 W quartz-tungsten-halogen (QTH) bulb (~225 mm above the surface) to mimic the irradiant shortwave energy and spectral pattern of the sun [31]. The externally mounted module ensured clear separation of shortwave and longwave bands thereby closely replicating actual conditions. The thermal heat generated by the QTH bulb was removed by a vacuum-fan suction system [32]. The power controller was adjusted to ensure that the sample received the desired shortwave irradiance, as measured by the pyranometer.

Research Methodology
Vapour flux tests were conducted based on atmospheric parameters in the study area. For each test, the sample container was cleaned with distilled water, allowed to air dry, and placed on the analytical scale balance. The balance was tared with the empty container and approximately 15 mL of room temperature distilled water was added through a syringe. Data logging sensors were turned on and the climate chamber lid was closed, and modules relevant to the active weather scenario were engaged. After the desired atmospheric conditions in the chamber were attained, the lid was re-opened and~15 mg of water added. The evaporation tests were conducted for 3.0 h, similar to [33], to obtain sufficient data while maintaining constant meteorological conditions in the BAS. This generated 1100 measurements at 10-s intervals. The average vapour flux over the course of each experiment was determined using the change in sample mass over time, and the surface area of the sample. Table 1 provides the average atmospheric parameters for the selected weather scenarios (spring, summer, and fall) showing both daytime and nighttime values. Winter was not included because freezing conditions are prevalent during this season. The atmospheric parameters (hourly land-based measurements) were obtained from the Canadian Weather Energy and Engineering Datasets (CWEEDS). Based on data from 1998 to 2014, the variations of atmospheric parameters were plotted in the form of mean values and standard deviation over the year. The selected value of each parameter was horizontally extended and a distinct change in standard deviation was identified as the cut-off date for a given season. Figure 4 presents typical annual daytime conditions in the study area with the summer weather scenario highlighted. The atmospheric parameters represent the weighted average values between the stations surrounding Lake Diefenbaker: Regina (172 km), Swift Current (117 km), Wynyard (196 km), and Kindersley (170 km). The measured values of atmospheric parameters at a station were divided by their respective distances from Lake Diefenbaker and averaged using the sum of the inverse distances following the method [34]. The presented data were used to simulate evaporative conditions at the water surface. The air velocity (Figure 4a) data were downscaled using the method described by [35].
Likewise, the air humidity (Figure 4b), presented as vapour density, was calculated from partial vapour pressure using dewpoint temperature and air temperature, as given in Figure 4c [36]. Solar irradiance (Figure 4d) represents the sum of direct and diffuse radiant energy received at the ground surface. Finally, the water surface temperature (Figure 4e) of a 1.0 km 2 section on Lake Diefenbaker was obtained using the Oak Ridge National Laboratory (ORNL) MODIS Land Product Subset Tool-MYD21A2 that contains 17-years (2002 to 2019) of 8-day water-based measurements [37].     Table 2 provides a summary of the average atmospheric and water surface parameter measurements in the BAS for the investigated weather scenarios. The standard error over the test duration for each of the parameters was found to be negligible, thereby demonstrating that the BAS is able to maintain steady-state atmospheric conditions. The target air velocity (Table 1) was achieved in the BAS. Depending on the weather scenario, the air pressure measured from about 93 kPa to 96 kPa. The variation in relative humidity and air temperature among the four sensors is attributed to their fetch location with respect to the sample. Generally, the upwind and high sensors measured low relative humidity and high air temperature compared to their downwind and low counterparts. This is attributed to wind direction (which laterally moves the water vapour released from the sample surface [38]) and elevation (which retards the vertical movement of water vapour released from the sample [39]). Based on the average of the lower sensors, the target air temperature (Table 1) was achieved in the BAS. With the target incoming shortwave flux ( Table 1) achieved in the BAS, the outgoing shortwave flux was found to be lower than the suggested albedo coefficient for open water bodies [5]. This is primarily attributed to the stationary and perpendicular flux source as opposed to the moving and angular direction of the sun [40]. The rate of mass change due to evaporation over 3-h was obtained using ten-degree polynomial regression curves that best-fitted the measured data. The resulting values were found to range between 7 × 10 −3 g·s −1 to 18 × 10 −3 g·s −1 . Finally, the surface temperature followed the expected seasonal and day-night trends that may be partially attributed to the specific heat capacity of water. Given the short test duration and small sample size, the measured data do not capture the long-term heat storage effects of deep water bodies [41].

Results and Discussion
Note. The ± indicates the standard error (SE) of measurement within the 3-h test period (n = 1100) of the weather scenario a Applies to the fitted high-order polynomial curve used to determine rate of mass change (∆M). Figure 5 presents the typical atmospheric and water surface conditions measured over the test duration in the BAS for a summer day. The air velocity (Figure 5a) matched the target value of 1.3 m·s −1 with minor fluctuations. The downward trend in air humidity (Figure 5b) from the target value of 9 g·m −3 is attributed to possible interference with the unregulated laboratory humidity. Such interferences were not observed with air temperature (Figure 5c), which was controlled in the laboratory. The target air temperature of 19 • C was achieved, as measured by the lower sensors. The variation in measured relative humidity and air temperature is due to the interrelationship of the two parameters [42]. The solar irradiance (Figure 5d) measurements before and after testing showed no deviation from the set calibration voltage and achieved the target value. Finally, the water surface temperature (Figure 5e) was noted to be stable during testing with the sporadic spikes attributable to the relatively lower precision (1 • C) of the pyrometer. Table 3 gives a summary of the analysed data of average atmospheric and water surface parameters for the investigated weather scenarios. The corresponding transient data (not given in this paper) were found to be steady. From a momentum perspective, the aerodynamic resistance is a pertinent parameter affecting evaporation in the BAS. This parameter captures the atmospheric turbulent diffusion of mass and energy and is inversely related to wind velocity and cover height [5]. The values of aerodynamic resistance were found to range from 41 s·m −1 to 47 s·m −1 (typical for water surfaces) such that the lower values in spring and fall were associated with high wind velocity and the reverse was true for summer.  Table 3 gives a summary of the analysed data of average atmospheric and water surface parameters for the investigated weather scenarios. The corresponding transient data (not given in this paper) were found to be steady. From a momentum perspective, the aerodynamic resistance is a pertinent parameter affecting evaporation in the BAS. This parameter captures the atmospheric turbulent diffusion of mass and energy and is inversely related to wind velocity and cover height [5]. The values of aerodynamic resistance were found to range from 41 s•m −1 to 47 s•m −1 (typical for water surfaces) such that the lower values in spring and fall were associated with high wind velocity and the reverse was true for summer. The target air humidity (Table 1) was achieved in the BAS as represented by vapour density. The vapour pressure deficit characterises the ease of evaporative mass transfer based on the residual water vapour capacity in the atmosphere that, in turn, depends on the temperature source and its measurement location [43]. The vapour pressure deficits were found to be generally high during the day as compared to those at night. In the BAS atmosphere, the data indicated large variation in vapour pressure deficit during summer (1047 Pa during the day and 389 Pa at night), low variation during spring (588 Pa during the day and 483 Pa at night), and negligible change fall (549 Pa during the day and 559 Pa at night). Likewise, the surface water data indicated large variations in vapour pressure deficit during spring (1037 Pa during the day and 353 Pa at night) and summer (898 Pa during the day and 190 Pa at night), and low variation during fall (559 Pa during the day and 439 Pa at night). The separate determination of air and surface deficits allows for their subsequent use in the various empirical equations, as described later in this paper.
Assuming an infinitely thin surface with no heat storage, the available energy (Q) was calculated according to the following equation [44]: In the above equation, net radiant heat flux (R) was determined by measuring shortwave energy ( Table 2) and estimating longwave energy (Table 3). Likewise, sensible heat flux (H) was determined using the Bowen Ratio and evaporative heat flux (λE). The ground heat flux (G) was determined by rearranging Equation (1) and solving using the measured and analysed variables. The available energy was found to be highest during spring (328 J·s −1 ·m −2 during the day and 67 J·s −1 ·m −2 at night), followed by summer (252 J·s −1 ·m −2 during the day and 30 J·s −1 ·m −2 at night), and then by fall (241 J·s −1 ·m −2 during the day and 102 J·s −1 ·m −2 at night).
The vapour flux was obtained from the measured rate of change in mass and the surface area. The data followed seasonal patterns similar to the atmospheric vapour pressure deficit, namely: a large variation during summer (0.120 g·s −1 ·m −2 during the day and 0.047 g·s −1 ·m −2 at night), a low variation during spring (0.116 g·s −1 ·m −2 during the day and 0.062 g·s −1 ·m −2 at night), and a negligible change during fall (0.100 g·s −1 ·m −2 during the day and 0.076 g·s −1 ·m −2 at night). Figure 6 presents the effect of measured and analysed parameters on vapour flux. For measured parameters, the vapour flux showed opposite patterns for the night when compared to the day. This highlights the significance of capturing solar irradiance in the BAS, albeit more data is required to validate the observed patterns. Despite data scarcity, vapour flux generally increased with an increase in both vapour pressure deficit and energy, such that the night values were lower than the day values. The scatter in aerodynamic resistance is attributed to that in measured air velocity. The vapour flux was obtained from the measured rate of change in mass and the surface area. The data followed seasonal patterns similar to the atmospheric vapour pressure deficit, namely: a large variation during summer (0.120 g•s −1 •m −2 during the day and 0.047 g•s −1 •m −2 at night), a low variation during spring (0.116 g•s −1 •m −2 during the day and 0.062 g•s −1 •m −2 at night), and a negligible change during fall (0.100 g•s −1 •m −2 during the day and 0.076 g•s −1 •m −2 at night). Figure 6 presents the effect of measured and analysed parameters on vapour flux. For measured parameters, the vapour flux showed opposite patterns for the night when compared to the day. This highlights the significance of capturing solar irradiance in the BAS, albeit more data is required to validate the observed patterns. Despite data scarcity, vapour flux generally increased with an increase in both vapour pressure deficit and energy, such that the night values were lower than the day values. The scatter in aerodynamic resistance is attributed to that in measured air velocity.   Table 4 gives the established empirical equations for potential evaporation relevant to the collected BAS data. The mass-transfer models used equations requiring variables of momentum and mass. These are based on Dalton's law of partial pressures and the eddy motion transfer of vapour from water surfaces into the atmosphere [45]. In contrast, the combination models included energy variables. As such, these models are designed to accommodate the effects of solar radiation and ground heat flux [6].
The measured vapour flux was generally within one standard deviation of the equality line when compared with that predicted by both the mass-transfer equations ( Figure  7a) and the combination equations (Figure 7b). The underestimated predictions are at-   Figure 7 gives a comparison of the measured values with estimates based on empirical relationships and Table 4 gives the established empirical equations for potential evaporation relevant to the collected BAS data. The mass-transfer models used equations requiring variables of momentum and mass. These are based on Dalton's law of partial pressures and the eddy motion transfer of vapour from water surfaces into the atmosphere [45]. In contrast, the combination models included energy variables. As such, these models are designed to accommodate the effects of solar radiation and ground heat flux [6].  Tables 2 and 3. RMSE is root mean square error, R 2 is correlation coefficient, and SI is scatter index.
tributed to a downscaled air velocity in the BAS and/or the selected constants in the empirical equations. Among the former set, the equation by [46] was found to be closest to the measured data owing to a better prediction of the wind function. Likewise, the combination equation by [47] best matched the measured values due to the inclusion of atmospheric vapour pressure deficit. Finally, the vapour flux measured up to 0.12 g•s −1 •m −2 and this number does not pertain to any limitations of BAS.  The measured vapour flux was generally within one standard deviation of the equality line when compared with that predicted by both the mass-transfer equations (Figure 7a) and the combination equations (Figure 7b). The underestimated predictions are attributed to a downscaled air velocity in the BAS and/or the selected constants in the empirical equations. Among the former set, the equation by [46] was found to be closest to the measured data owing to a better prediction of the wind function. Likewise, the combination equation by [47] best matched the measured values due to the inclusion of atmospheric vapour pressure deficit. Finally, the vapour flux measured up to 0.12 g·s −1 ·m −2 and this number does not pertain to any limitations of BAS. Figure 8 gives predicted evaporation for the study area and Table 5 (Table 3) and duration (Table 1) along with water density. During the day time, the BAS data predicted 1459 mm/year of evaporation that was 9% higher than that estimated by the mass-based Meyer equation but identical to the combinationbased Monteith equation. During the night time, the BAS data predicted 322 mm/year of evaporation, which was over predicted by 39% and 27%, respectively. This might be attributed to less adequate chamber insulation for relatively cooler scenarios. Nonetheless, the average daily evaporation ranged from 4 mm·day −1 to 8 mm·day −1 during the day and decreased to 1 mm·day −1 to 3 mm·day −1 at night. The combined 24-h evaporation was found to be 8 ± 1 mm·day −1 from late April through late October. This is similar to the estimated data reported by [52] for Lake Diefenbaker. Furthermore, the cumulative annual evaporation was found to be 1781 mm, of which 82% occurs during the day and 18% during the night. Overall, the BAS over predicted the cumulative annual evaporation from the above-mentioned estimates by 15% and 4%, respectively.

Summary and Conclusions
Knowledge of evaporation is critical for efficient water management in semi-arid climates such as in the Canadian Prairies. This research focused on determining evaporative fluxes using a de novo bench-scale atmosphere simulator that specifically captured the shortwave energy and spectral pattern of the sun. The system is capable of imposing atmospheric conditions on a sample of distilled water. While there are inherent limitations in laboratory investigations, this undertaking allowed a clear understanding of complicated interactions of mass and energy at water surfaces. This means that the findings correlate well with empirical relationships and, as such, are useful in water resources management. Based on regional atmospheric conditions, seasonal weather scenarios representing spring, summer, and fall, were simulated under controlled air velocity, humidity, temperature, and solar irradiance. The main conclusions of this study are given below:

•
All of the measured data achieved the target values for the various parameters and the data were found to be stable during the 3-h test duration. The slight downward trend in air humidity with respect to the target value is attributed to interferences with the unregulated laboratory humidity.  The BAS effectively captured horizontal advection (fetch effects) through an adequate sensor placement. The analysed data must be understood in the context of land-surface transitions at shorelines, heat advection from inflowing water, and variable heat storage across the reservoir [53]. The seasonal values are partly affected by the long-term heat storage in the investigated water body [41]. In the absence of surface-atmosphere monitoring data, laboratory simulations can only be correlated with relevant empirical equations. Although the BAS cannot measure three-dimensional air velocity (eddy covariance) or be used for long-term studies, it can simulate various weather parameters that are required for the estimation of evaporation. This will help in informed decision making related to water storage and usage in the region, especially when coupled with field observations.

Summary and Conclusions
Knowledge of evaporation is critical for efficient water management in semi-arid climates such as in the Canadian Prairies. This research focused on determining evaporative fluxes using a de novo bench-scale atmosphere simulator that specifically captured the shortwave energy and spectral pattern of the sun. The system is capable of imposing atmospheric conditions on a sample of distilled water. While there are inherent limitations in laboratory investigations, this undertaking allowed a clear understanding of complicated interactions of mass and energy at water surfaces. This means that the findings correlate well with empirical relationships and, as such, are useful in water resources management. Based on regional atmospheric conditions, seasonal weather scenarios representing spring, summer, and fall, were simulated under controlled air velocity, humidity, temperature, and solar irradiance. The main conclusions of this study are given below:

•
All of the measured data achieved the target values for the various parameters and the data were found to be stable during the 3-h test duration. The slight downward trend in air humidity with respect to the target value is attributed to interferences with the unregulated laboratory humidity.

•
The vapour flux was found to have large variation during summer (0.120 g·s −1 ·m −2 during the day and 0.047 g·s −1 ·m −2 at night), low variation during spring (0.116 g·s −1 ·m −2 during the day and 0.062 g·s −1 ·m −2 at night), and negligible change during fall (0.100 g·s −1 ·m −2 during the day and 0.076 g·s −1 ·m −2 at night).

•
The measured vapour flux was generally within one standard deviation of the equality line when compared with that predicted by both the mass-transfer equations and the combination equations. The underestimated predictions are attributed to a downscaled air velocity in the simulator and/or the selected constants in the empirical equations.

•
The average evaporation ranged from 4 mm·d −1 to 8 mm·d −1 during the day and decreased to 1 mm·d −1 to 3 mm·d −1 at night. The 24-h evaporation was found to be 8 ± 1 mm·d −1 from late April through late October. Likewise, the cumulative annual evaporation was found to be 1781 mm, of which 82% occurs during the day and 18% at night. Informed Consent Statement: Not applicable.

Data Availability Statement:
The authors can provide access to sensor data and test data upon request.