Exploring the Effects of Atmospheric Forcings on Evaporation: Experimental Integration of the Atmospheric Boundary Layer and Shallow Subsurface

Podsumowanie

A protocol for the design and construction of a soil tank interfaced to a small climate controlled wind tunnel to study the effects of atmospheric forcings on evaporation is presented. Both the soil tank and wind tunnel are instrumented with sensor technologies for the continuous in situ measurement of environmental conditions.

Streszczenie

Evaporation is directly influenced by the interactions between the atmosphere, land surface and soil subsurface. This work aims to experimentally study evaporation under various surface boundary conditions to improve our current understanding and characterization of this multiphase phenomenon as well as to validate numerical heat and mass transfer theories that couple Navier-Stokes flow in the atmosphere and Darcian flow in the porous media. Experimental data were collected using a unique soil tank apparatus interfaced with a small climate controlled wind tunnel. The experimental apparatus was instrumented with a suite of state of the art sensor technologies for the continuous and autonomous collection of soil moisture, soil thermal properties, soil and air temperature, relative humidity, and wind speed. This experimental apparatus can be used to generate data under well controlled boundary conditions, allowing for better control and gathering of accurate data at scales of interest not feasible in the field. Induced airflow at several distinct wind speeds over the soil surface resulted in unique behavior of heat and mass transfer during the different evaporative stages.

Wprowadzenie

Understanding the interaction between the land and atmosphere is paramount to our understanding of many current world problems such as leaking of geologically-sequestered carbon dioxide in soil, climate change, water and food supply, the accurate detection of landmines, and the remediation of ground water and soil. In addition, the primary exchanges of heat and water that drive global and regional meteorological conditions occur at the Earth’s surface. Many weather and climate phenomena (e.g., hurricanes, El Niño, droughts, etc.) are principally driven by processes associated with atmospheric-land surface interactions¹. As more than half of the land surface on the Earth is arid or semiarid^2-4, accurately describing the water cycle in these regions on the basis of heat and water exchanges between the atmospheric air and the soil surface is critical to improving our understanding of the aforementioned issues, particularly in regions vulnerable to extended drought and desertification. However, despite decades of research, there still remain many knowledge gaps in the current understanding of how the shallow subsurface and atmosphere interact⁵.

Transport processes involving liquid water, water vapor, and heat in soil are dynamic and strongly coupled with respect to interactions with the soil and enforced boundary conditions (i.e., temperature, relative humidity, thermal radiation). Numerical heat and mass transfer models commonly oversimplify or overlook a number of these complexities due in part to a lack of testing and refinement of existing theories resulting from a paucity of high temporal and spatial resolution data. Datasets developed for model validation are oftentimes lacking critical atmospheric or subsurface information to properly test the theories, resulting in numerical models that do not properly account for important processes or depend on the use of poorly understood parameters that are adjusted or fitted in the model. This approach is widely used due to its simplicity and ease of use and has in some applications shown much merit. However, this approach can be improved upon by better understanding the physics behind these “lumped parameterizations” by performing well controlled experiments under transient conditions that are capable of testing heat and water transfer theory⁶.

Careful experimentation in the laboratory allows precision datasets to be generated that can subsequently be used to validate numerical models. Data available from field sites are often incomplete and costly to obtain, and the degree of control needed to obtain a fundamental understanding of processes and to generate data for model validation could be considered inadequate in some cases. Laboratory experimentation of natural phenomena such as soil evaporation allows atmospheric conditions (i.e., temperature, relative humidity, wind speed) and soil conditions (i.e., soil type, porosity, packing configuration) to be carefully controlled. Many laboratory techniques used to study soil evaporation and soil thermal and hydraulic properties use destructive sampling^7-10. Destructive sampling methods require that a soil sample be unpacked to obtain point data, preventing the measurement of transient behavior and disrupting soil physical properties; this approach introduces error and uncertainty to the data. Nondestructive measurements, like the method presented here, allow for more accurate determination and study of the interdependency of soil properties and processes¹¹.

The goal of this work is to develop a soil tank apparatus and associated protocol for the generation of high spatial and temporal resolution data pertaining to the effects of changes in atmospheric and subsurface conditions on bare-soil evaporation. For this work, a small wind tunnel capable of maintaining a constant wind speed and temperature is interfaced with a soil tank apparatus. The wind tunnel and soil tank are instrumented with a suite of state of the art sensor technologies for autonomous and continuous data collection. Wind speed is measured using a stainless steel pitot-static tube attached to a pressure transducer. Temperature and relative humidity are monitored in the atmosphere using two types of sensors. Relative humidity and temperature are also monitored at the soil surface. Sensors in the subsurface measure soil moisture and temperature. Weight measurements of the tank apparatus are used to determine evaporation through a water mass balance. To demonstrate the applicability of this experimental apparatus and protocol, we present an example of bare-soil evaporation under varying wind speed conditions. The soil tank, packed homogeneously with a well characterized sand, was initially fully saturated and allowed to evaporate freely under carefully controlled atmospheric conditions (i.e. temperature, wind speed).

Protokół

Note: Laboratory testing is performed using a two-dimensional bench scale tank interfaced with a climate controlled wind tunnel apparatus. Both the bench scale tank and wind tunnel are instrumented with various sensor technologies. The following protocol will first discuss the construction and preparation of the soil tank, followed by a discussion of the wind tunnel and the instrumentation of both. The tank dimensions, wind tunnel dimensions, number of sensors, and sensor technology type presented can be modified to suit the needs of a specific experimental set-up. The protocol presented below was used to experimentally study the effects of wind speed on bare-soil evaporation.

1. Construction and Preparation of Porous Media Soil Tank

1. Cut a large piece of 1.2 cm thick acrylic glass into five individual panes. Assemble these panes into an open-topped soil tank with internal length, width and height of 25, 9.1 and 55 cm, respectively. Acrylic glass allows processes in the subsurface to be visually observed.
2. Draw a 5 x 5 grid that is 25 cm by 25 cm on each of the two large glass panes (length 25 cm and height 55 cm) as shown in Figure 1. Ensure that each square within the grid has an area of 25 cm² (Figure 1). The grid will be used to properly space the sensors within the soil tank.

Figure 1: Schematic front and side views of the soil tank used for the experimental set-up (dimensions are in centimeters). (a) The front view of the soil tank displaying the grid system consisting of twenty-five 5 cm x 5 cm squares. (b) The side view of the soil tank, showing the installed temperature, relative humidity and soil moisture sensor network as a function of depth. Note that the schematics are not drawn to scale.

3. On one of the large glass planes, drill a total of twenty-five 1.9 cm (¾ inch) diameter holes for the soil moisture sensors.
   1. Drill each hole in the center of every square in the grid established in Step 1.2 so that the centers of the holes of two abutting squares are 5 cm apart; the first set of holes is 2.5 cm below the top of the tank. Use appropriately sized taps to cut threads into each of the newly created holes. The 5 cm spacing between sensors ensures that each sensor is outside of the sampling volume of the next closest sensor.
4. Similarly, drill and tap a total of twenty-five 0.635 cm (¼ inch) diameter holes in the center of each grid box created during Step 1.2. Make sure that the center of each hole is spaced 5 cm apart with the first row of holes located 2.5 cm below the top of the soil tank. The 5 cm spacing between sensors ensures that each sensor is outside of the sampling volume of the next closest sensor.
5. On the acrylic pane used as the bottom of the tank, drill and tap a single ½ inch diameter hole in the middle of the pane. Glue a mesh screen (finer than the test soils to be used) over the hole on the internal side of the glass. On the external side of the bottom plane, install a 90° elbow that is attached to flexible tubing with an adjustable valve. This valve and tubing is used to drain water from the tank at the termination of an experiment or as a way to install constant head devices for maintaining constant water table depths.
6. Use marine grade glue or similar water resistant polymer adhesive to attach and seal the tank together as shown in Figure 1. Allow the adhesive to cure for one day.
7. To raise the tank off of the ground and make room for the 90° elbow (Figure 1), attach two additional pieces of 1.2 cm thick acrylic glass with length 12 cm and height 5 cm to the bottom of the tank.

2. Construction and Preparation of Climate Controlled Wind Tunnel

1. Construct the 215 cm long upstream portion of the wind tunnel out of rectangular galvanized steel ducting material that has a width of 8.5 cm and a height of 26 cm. Surround the outside of the duct with polystyrene insulation.
2. Drill a small hole in the side of the duct work near the downstream exit of the upstream portion of the wind tunnel for the insertion of a relative humidity-temperature sensor (Figure 2).

Figure 2: Complete experimental set-up, including tank, ductwork, sensors grid (dimensions are in centimeters). Complete experimental set up of the combined wind tunnel and soil-tank apparatus. The wind tunnel is elevated and sits flush with the surface of the soil tank. The soil tank is instrumented with a network of sensors used to measure a variety of subsurface and atmospheric variables. The grid circles represent the locations for inserting these sensors. A heating control system and an in-line duct fan are used to control temperature and wind speed, respectively. The pitot-static tube is used to measure wind speed. The entire apparatus sits on a weighting scale to obtain a mass balance during experimentation. Note that the schematic is not drawn to scale.

3. Install five ceramic infrared heating elements positioned in parallel within a reflector along the length of the upstream portion of the wind tunnel. Connect the infrared heating elements to a temperature control system regulated by an infrared temperature sensor.
4. Construct the mid-section of the wind tunnel out of two 1.2 cm thick acrylic panels with a length and height of 25 cm and 26 cm respectively. Drill two 0.635 cm (¼ inch) diameter holes in one of the mid-section panels to insert temperature and/or relative humidity-temperature sensors at the locations shown in Figure 2.
   1. Secure the acrylic panels to the top of the soil tank side walls (i.e., panels with dimensions 25 cm x 55 cm) using a strong adhesive tape, ensuring that the wind tunnel and soil tank panels sit flush with one another.
5. Construct the first 50 cm of the downstream portion of the wind tunnel out of the same size rectangular ducting material described in Step 2.1. On the terminating side, reduce the rectangular ducting material to a 15.3 cm diameter round duct with length of 170 cm. Install a galvanized steel damper, used to adjust wind speeds, at the far downstream end of the round duct for aid in wind speed control.
6. As in Step 2.2, drill one 0.635 cm diameter hole in the side of the downstream rectangular duct near the entrance for the insertion of a relative humidity-temperature sensor. Drill a second 0.635 cm diameter hole from the top of the rectangular duct along the centerline of the wind tunnel.
7. Install an in-line duct fan in the middle of the round duct (i.e., 85 cm downstream from reduction described in Step 2.4) oriented to expel air from the downstream portion of the wind tunnel. Interface the fan with a variable speed controller for more precise control of rotational frequency and as a result wind speed.
8. Use weldment material and adjustable shelving units to elevate and secure the wind tunnel apparatus. Ensure that the bottom of the upstream and downstream ductwork are flush with the top of the soil tank (Figure 2).

3. Installation of Sensors

1. Prior to installation within the soil tank, secure each soil moisture and temperature sensor within a threaded NPT housing (1.9 cm and 0.635 cm housings, respectively) and seal with flashing sealant to prevent moisture intrusion. Do not use silicone-based sealant products as they can interfere with the electronics within some sensors. Cure the sensors for approximately one week.
2. Prior to installation in the soil tank, calibrate the soil moisture sensors in accordance with the two point α-mixing method developed by Sakaki et al.¹².
3. Wrap the threads of each NPT housing with plumbers tape prior to installation in the tank to help provide a better seal between the NPT threading and acrylic glass.
4. Install a total of 25 soil moisture and temperature sensors each horizontally through the walls of the tank at the locations discussed in Step 1.2. Twist the sensor cables in sync with the NPT fitting/housing so as not damage the internal wiring within the cables. Do not over-torque the NPTs so as to prevent the glass from cracking. Connect the soil moisture sensors and temperature sensors to their designated data loggers.
5. Install 3 relative humidity-temperature sensors on the soil surface at distances of 2.5, 12.5 and 21.5 cm from the leading edge of the tank. Place the sensors in good contact with the soil surface so that the relative humidity readings reflect the conditions at the soil surface rather than the surrounding air. Connect the sensors to data loggers.
6. To obtain the requisite air temperature and relative humidity measurements in the atmosphere, install relative humidity-temperature sensors in the free flow section of the wind tunnel, using the holes drilled through the upstream and downstream sections of the wind tunnel as well as the panels.
7. Install a pitot-static tube directly downstream of the soil tank through the 0.635 cm hole drilled in the top of the downstream wind tunnel section. Hold the pitot-static tube at a height of 13 cm from the floor of the section. Connect the pitot-static tube to a differential pressure transducer.
8. Calibrate the differential pressure transducer. The pitot-static tube measures dynamic pressure which is defined as the difference of the stagnation and static pressures. The pressure differential is interpreted by the pressure transducer as a voltage differential.
   1. Measure the voltage under no flow conditions (voltage should be approximately equal to 0) and for flow of a known dynamic pressure; this allows a linear relationship to be established between dynamic pressure and voltage. Determine wind speed by applying Bernoulli’s equation:
      (1)
      where V (m/sec) is wind speed, P_dynamic (Pa) is dynamic pressure, and ρ (kg/m³) is the density of air.
   2. Compare the velocity calculated using Equation (1) with another measurement device. Here, compare the pitot-static tube differential pressure transducer with Laser Doppler Velocimetry (LDV) measurements which has an accuracy of ± 0.01 m/sec.
      Note: A summary of sensors employed and their associated sampling frequencies can be found in Table 1. For sensor specifications and other information, refer to the enclosed materials/equipment list.

Sensor	Sensor Measurements	Number of Sensors Employed in Experimental Apparatus	Sensor Sampling Frequency (min)
EC-5	Soil moisture	25	10
ECT	Soil/air temperature	25	10
SH-1	Thermal properties	1	10
EHT	Relative humidity/temperature	5	10
Infrared camera	Surface temperature/evaporation	1	1
Digital camera	Visualization of drying front	1	60
Pitot static tube	Wind velocity	1	10
Weighting scale	Cumulative evaporation/evaporation rate	1	10

Table 1: Summary of sensors used in experimental portion of present study.

4. Pack the Soil Tank and Prepare for the Start of the Experiment

1. Prior to packing the tank with soil, test its integrity by performing a leak test. Fill the tank with water and wait for 4-6 hr to ensure that no leaks in the structure or sensors have developed.
   1. If leaks develop, drain the tank, let it dry overnight and fix the leaks using the same marine adhesive used during original construction. If no leaks develop, drain the soil tank and prepare for the steps below.
2. Determine the total volume of the tank with sensors in place. Carefully fill the tank with water using a graduated cylinder, making sure to record the amount of water added. Convert the recorded total volume to cubic centimeters for use in Step 4.5.
3. Obtain dry soil to pack the soil tank. Characterize the hydraulic and thermal properties of the selected soil separately in accordance with the methods discussed in Smits et al.¹¹
4. Carefully wet-pack the soil tank using soil and deionized water.
   1. To wet-pack the soil tank, first pour approximately 5 cm of water into the tank. Slowly add dry soil to the water in the tank, using a scoop, in 2.5 cm depth increments. Record the weight of the sand added during each lift so the porosity of the soil packing can be calculated.
   2. Upon completion of each layer, repeatedly tap the tank walls using a rubber mallet, 100-200 times, to obtain a uniform bulk density throughout. While tapping, avoid contact with the sensors and sensor wires. The use of vibratory devices should be avoided so as not to damage the network of sensitive sensors.
5. Upon completion of packing the tank, sum together the weights of each soil layer (see Step 4.4) to get the total mass of soil. Divide the total mass by the bulk density of soil (e.g. the bulk density of quartz sand is 2.65 g/cm³) to determine the volume of the sand (V_s, cm³). Calculate the porosity (η, m³/m³) of the soil in the tank according to:
   (2)
   where (V_T, m³) is the total volume of the empty tank determined in Step 4.2.
6. Once the tank is fully packed, place a plastic cover such as saran wrap over the tank until the experiment is ready to start to prevent the onset of evaporation.
7. Place the tank on a weighting scale to monitor cumulative water loss which can in turn be used to calculate evaporation rate.
8. Calculate the hourly evaporation rate by dividing the hourly weight loss by the product of the density of water and the cross sectional area of the evaporating surface.

5. Start the Experiment and Begin Data Collection

1. Once the set-up is complete, determine the desired atmospheric conditions (i.e. temperature, wind speed). Ensure that the data loggers and other data acquisition systems are turned on and set to the correct sampling intervals (e.g., every 10 min).
2. Start the fan and temperature control system. Allow the climate conditions to equilibrate before removing the plastic cover on the surface of the soil tank. Run the experiment for the desired length of time (e.g., 15 days).

Wyniki

The objective of the experiment presented here was to study the effect of wind speed on evaporation from bare soil. Key properties of the test soil used in the present study are summarized in Table 2. A series of experiments were performed in which different boundary conditions at the soil surface (i.e., wind speed and temperature) were applied (Table 3). Although four experiments at different wind speeds and temperatures were performed, the majority of experimental results pres...

Dyskusje

The purpose of this protocol was to develop an experimental apparatus and associated procedures for the generation of high spatial and temporal resolution data required for studying land-atmospheric interactions with respect to heat and mass transfer processes. The experimental apparatus described consisted of a soil tank and a small wind tunnel, both of which were outfitted with an array of sensors for the measurement of pertinent soil and atmospheric variables (e.g., wind speed, relative humidity, soil and air...

Materiały

Name	Company	Catalog Number	Comments
ECH2O EC-5 Soil Moisture Sensor (25)	Decagon Devices Inc. Decagon.com	40593	For specifics visit: http://www.decagon.com/products/soils/volumetric-water-content-sensors/ec-5-soil-moisture-small-area-of-influence/. Sampling frequency on 10 minute intervals, accuracy is ±3%, and collect data using the Em50 dataloggers
ECT Soil/Air Temperature Sensor (19)	Decagon Devices Inc. Decagon.com	40651	For specifications visit http://www.decagon.com/products/canopy-atmosphere/temperature/ect-air-temperature/. Sampling frequency on 10 min intervals, accuracy is ±0.5 °C, Measure within a temperature of 5 and 40 °C, and collect data using the Em50 dataloggers
EHT Relative Humidity and Temperature Sensor (5)	Decagon Devices Inc. Decagon.com	N/A	Sampling Frequency on 10 min intervals, accuracy is ±3% between 5% and 100% relative humidity, and collect data using Em50 data loggers. For more information visit decagon.com
Em50 Data Logger (10)	Decagon Devices Inc. Decagon.com	40800	For specifics visit http://www.decagon.com/products/data-management/data-loggers/em50-digital-analog-data-logger/. ECH2O decagon devices, pulls data from the ECT, EC-5, and EHT sensors, and each data logger has 5 sensor connections and a com port that connects from the logger to USB to computer
Sartorius Weighing Scale (1)	Sartorius Corporation	11209-95	Sartorius Model 11209-95, Range = 65 kg, Resolution = ±1 g
Infrared SalamandernCeramic Radiative Heater (1)	Mor Electric Heating Assoc., Inc. http://www.morelectricheating.com/	FTE 500-240	5 heaters needed, adjust to get the right ambient/free-flow temperature
2104 Temperature Control System (1)	Chromalox	2104	Controls the heaters
Infrared Temperature Sensor Regulator (1)	Exergen Corporation	N/A	Monitors the heaters temperatures
[header]
Stainless Steel Pitot-Static Tube (1)	Dwyer Instruments, Inc. http://www.dwyer-inst.com/	Series 160	For specifics visit http://www.dwyer-inst.com/Product/%20TestEquipment/PitotTubes/Series160. Sensor sampling frequency is every 10 minutes, must be connected to differential pressure transducer and anemometer, and convert the pressure data collected into win velocities using Bernoulli's equation.
1/2 inch Acrylic (1)	Colorado Plastics http://www.coloradoplastics.com/	N/A	Specific heat of 1,464 J kg-1 K-1, thermal conductivity of 0.2 W m-1 K-1, and a density of 1,150 kg m-3
Galvanized Steel Ducting Material (1)	Home Depot	N/A	Material used to build wind- tunnel, and both round and rectangular ducting were used in construction and connected using square-to-round reducer duct
Variable Speed Controller Connected to an In-Line Duct Fan (1)	Suncourt, Inc. http://www.suncourt.com/	VS200	15.3 cm in Diameter Placed in-line with round duct
Galvanized Steel Damper (1)	Home Depot	N/A	Used to control/reduce speeds in the wind tunnel for low velocity data
Accusand #30/40 (1)	Unimin Corporation http://www.unimin.com/	N/A	This sand is silica sand and is 99.8% quartz, its grain shape is classified as rounded, the uniformity coefficient is approximately 1.2, and the grain density is 2.66 g/cm3.