Name: exploring the effects of atmospheric forcings on evaporation: experimental integration of the atmospheric boundary layer and shallow subsurface
Uploaded: 2015-06-08T13:00:00.000+00:00
Duration: 13 min 27 s
Description: Watch this Scientific Journal Video about Exploring the Effects of Atmospheric Forcings on Evaporation: Experimental Integration of the Atmospheric Boundary Layer and Shallow Subsurface at JoVE.com

This research was funded by the U. S. Army Research Office Award W911NF-04-1-0169, the Engineering Research and Development Center (ERDC) and National Science Foundation grant EAR-1029069. In addition, this research was supported by a Summer Programs in Undergraduate Research grant from Colorado School of Mines. The authors wish to thank Ryan Tolene and Paul Schulte for their contributions.

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
<ol>
	<li>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.</li>
	<li>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 cm2 (Figure 1). The grid will be used to properly space the sensors within the soil tank.</li>
</ol>
<img alt="Figure 1" src="/files/ftp_upload/52704/52704fig1highres.jpg" width="700" /> 
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.
<ol start="3">
	<li>On one of the large glass planes, drill a total of twenty-five 1.9 cm (&#190; inch) diameter holes for the soil moisture sensors.
	<ol>
		<li>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.</li>
	</ol>
	</li>
	<li>Similarly, drill and tap a total of twenty-five 0.635 cm (&#188; 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.</li>
	<li>On the acrylic pane used as the bottom of the tank, drill and tap a single &#189; 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&#176; 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.</li>
	<li>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.</li>
	<li>To raise the tank off of the ground and make room for the 90&#176; 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.</li>
</ol>
2. Construction and Preparation of Climate Controlled Wind Tunnel
<ol>
	<li>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.</li>
	<li>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).</li>
</ol>
<img alt="Figure 2" src="/files/ftp_upload/52704/52704fig2highres.jpg" width="700" /> 
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.
<ol start="3">
	<li>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.</li>
	<li>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 (&#188; 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.
	<ol>
		<li>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.</li>
	</ol>
	</li>
	<li>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.</li>
	<li>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.</li>
	<li>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.</li>
	<li>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).</li>
</ol>
3. Installation of Sensors
<ol>
	<li>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.</li>
	<li>Prior to installation in the soil tank, calibrate the soil moisture sensors in accordance with the two point &#945;-mixing method developed by Sakaki et al.12.</li>
	<li>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.</li>
	<li>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.</li>
	<li>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.</li>
	<li>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.</li>
	<li>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.</li>
	<li>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.
	<ol>
		<li>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&#8217;s equation: 
		<img alt="Equation 1" fo:content-width="0.8in" src="/files/ftp_upload/52704/52704eq1.jpg" width="100" /> (1) 
		where V (m/sec) is wind speed, Pdynamic (Pa) is dynamic pressure, and &#961; (kg/m3) is the density of air.</li>
		<li>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 &#177; 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.</li>
	</ol>
	</li>
</ol>
<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr>
			<td>Sensor</td>
			<td>Sensor Measurements</td>
			<td>Number of Sensors Employed in Experimental Apparatus</td>
			<td>Sensor Sampling Frequency (min)</td>
		</tr>
		<tr>
			<td>EC-5</td>
			<td>Soil moisture</td>
			<td>25</td>
			<td>10</td>
		</tr>
		<tr>
			<td>ECT</td>
			<td>Soil/air temperature</td>
			<td>25</td>
			<td>10</td>
		</tr>
		<tr>
			<td>SH-1</td>
			<td>Thermal properties</td>
			<td>1</td>
			<td>10</td>
		</tr>
		<tr>
			<td>EHT</td>
			<td>Relative humidity/temperature</td>
			<td>5</td>
			<td>10</td>
		</tr>
		<tr>
			<td>Infrared camera</td>
			<td>Surface temperature/evaporation</td>
			<td>1</td>
			<td>1</td>
		</tr>
		<tr>
			<td>Digital camera</td>
			<td>Visualization of drying front</td>
			<td>1</td>
			<td>60</td>
		</tr>
		<tr>
			<td>Pitot static tube</td>
			<td>Wind velocity</td>
			<td>1</td>
			<td>10</td>
		</tr>
		<tr>
			<td>Weighting scale</td>
			<td>Cumulative evaporation/evaporation rate</td>
			<td>1</td>
			<td>10</td>
		</tr>
	</tbody>
</table>
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
<ol>
	<li>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.
	<ol>
		<li>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.</li>
	</ol>
	</li>
	<li>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.</li>
	<li>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.11</li>
	<li>Carefully wet-pack the soil tank using soil and deionized water.
	<ol>
		<li>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.</li>
		<li>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.</li>
	</ol>
	</li>
	<li>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/cm3) to determine the volume of the sand (Vs, cm3). Calculate the porosity (&#951;, m3/m3) of the soil in the tank according to: 
	<img alt="Equation 2" fo:content-width="0.8in" src="/files/ftp_upload/52704/52704eq2.jpg" width="100" /> (2) 
	where (VT, m3) is the total volume of the empty tank determined in Step 4.2.</li>
	<li>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.</li>
	<li>Place the tank on a weighting scale to monitor cumulative water loss which can in turn be used to calculate evaporation rate.</li>
	<li>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.</li>
</ol>
5. Start the Experiment and Begin Data Collection
<ol>
	<li>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).</li>
	<li>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).</li>
</ol>

<ol>
	<li>Verstraete, M. M., Schwartz, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Desertification+and+global+change.">Desertification and global change.</a> Vegetatio. 91, 3-13 (1991).</li><li>Warren, A., Adams, W. M., Goudie, A. S., Orme, A. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Desertification.">Desertification.</a> The Physical Geography of Africa. , 342-355 (1996).</li><li>Katata, G., Nagai, H., Ueda, H., Agam, N., Berliner, P. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Development+of+a+land+surface+model+including+evaporation+and+adsorption+processes+in+the+soil+for+the+land-air+exchange+in+arid+regions.">Development of a land surface model including evaporation and adsorption processes in the soil for the land-air exchange in arid regions.</a> J. Hydrometeorol. 8, 1307-1324 (2007).</li><li>Davarzani, H., Smits, K. M., Tolene, R., Illangasekare, T. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Study+of+the+effect+of+wind+speed+on+evaporation+from+soil+through+integrated+modeling+of+atmospheric+boundary+layer+and+shallow+subsurface.">Study of the effect of wind speed on evaporation from soil through integrated modeling of atmospheric boundary layer and shallow subsurface.</a> Water Resour. Res. 50, (2014).</li><li>Heitman, J. L., Horton, R., Ren, T., Nassar, I. N., Davis, D. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+test+of+coupled+soil+heat+and+water+transfer+prediction+under+transient+boundary+temperatures.">A test of coupled soil heat and water transfer prediction under transient boundary temperatures.</a> Soil Sci. Soc. Am. J. 72 (5), 1197-1207 (2008).</li><li>Gurr, C. G., Marshall, T. J., Hutton, J. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Movement+of+water+in+soil+due+to+a+temperature+gradient.">Movement of water in soil due to a temperature gradient.</a> Soil Sci. 74, 335-345 (1952).</li><li>Nassar, I. N., Horton, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Water+transport+in+unsaturated+non-isothermal+salty+soil:+experimental+results.">Water transport in unsaturated non-isothermal salty soil: experimental results.</a> Soil Sci. Soc. Am. J. 53, 1323-1363 (1989).</li><li>Prunty, L., Horton, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Steady-state+Temperature+Distribution+in+Nonisothermal+unsaturated+closed+soil+cells.">Steady-state Temperature Distribution in Nonisothermal unsaturated closed soil cells.</a> Soil Sci. Soc. Am. J. 58, 1358-1363 (1994).</li><li>Bachmann, J., Horton, R., Ren, T., van der Ploeg, R. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+the+thermal+properties+of+four+wettable+and+four+water-repellent+soils.">Comparison of the thermal properties of four wettable and four water-repellent soils.</a> Soil Sci. Soc. Am. J. 65 (6), 1675-1679 (2001).</li><li>Smits, K. M., Sakaki, T., Limsuwat, A., Illangasekare, T. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Thermal+conductivity+of+sands+under+varying+moisture+and+porosity+in+drainage-wetting+cycles.">Thermal conductivity of sands under varying moisture and porosity in drainage-wetting cycles.</a> Vadose Zone J. 9, 1-9 (2010).</li><li>Sakaki, T., Limsuwat, A., Smits, K. M., Illangasekare, T. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Empirical+two-point+&#945;-mixing+Model+for+calibrating+the+ECH2O+EC-5+soil+sensor+in+sands.">Empirical two-point &#945;-mixing Model for calibrating the ECH2O EC-5 soil sensor in sands.</a> Water Resources Research. 44, W00D08 (2008).</li><li>Shokri, N., Lehmann, P., Or, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Evaporation+from+layered+porous+media.">Evaporation from layered porous media.</a> J. Geophys. Res. 115, B06204 (2010).</li><li>Van Brakel, J., Mujumdar, A. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Mass+transfer+in+convective+drying.">Mass transfer in convective drying.</a> Advances in Drying. 1, 217-267 (1980).</li><li>Yiotis, A. G., Subos, A. G. B. o. u. d. o. u. v. i. s. A. K., Tsimpanogiannis, I. N., Yortsos, Y. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+effect+of+liquid+films+on+the+drying+of+porous+media.">The effect of liquid films on the drying of porous media.</a> AIChE J. 50, 2721-2737 (2004).</li><li>Ishihara, Y., Shimojima, E., Harada, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Water+vapor+transfer+beneath+bare+soil+where+evaporation+is+influenced+by+a+turbulent+surface+wind.">Water vapor transfer beneath bare soil where evaporation is influenced by a turbulent surface wind.</a> J. Hydrol. 131 (1-4), 63-104 (1992).</li><li>Lehmann, P., Assouline, S., Or, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Characteristic+lengths+affecting+evaporative+drying+of+porous+media.">Characteristic lengths affecting evaporative drying of porous media.</a> Phys. Rev. E. 77 (5 Pt 2), 056309 (2008).</li><li>Trautz, A. C., Smits, K. M., Schulte, P., Illangasekare, T. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Sensible+heat+balance+and+heat-pulse+method+applicability+to+in+situ+soil-water+evaporation.">Sensible heat balance and heat-pulse method applicability to in situ soil-water evaporation.</a> Vadose Zone J. 13, (2014).</li></ol>

The authors declare that they have no competing financial interests.

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 temperature and soil moisture). The following are some of the most critical components of the protocol presented in this study.
The tank dimensions and sensor placement were specifically chosen to maximize the number of sensors employed while accounting for the sensor&#8217;s respective sample volumes. The first row of sensors is 2.5 cm below the soil surface due to each sensor&#8217;s sample volume (defined as the volume of soil around the sensor, within which a change in environmental conditions affects the sensor readings). The sensors, placed in NPT fittings, are installed horizontally through the walls of the soil tank so that the sensor wires are not within the soil itself; all the sensor wires are outside of the tank, preventing water channeling. The installation of a large network of temperature and soil moisture sensors allows horizontal and vertical distributions of these variables to be determined at a fine spatial resolution.
Placing the soil tank on a weighting scale allows cumulative water loss and the associated evaporation rate to be determined using the water mass balance approach described above. These values can then be compared to evaporation rates obtained using other methods such as the combined heat-pulse and sensible heat balance method employed in Trautz et al.18
The wind tunnel portion of the apparatus is composed of three parts - an upstream, downstream and middle section. The upstream section is used to heat the air before it is drawn over the soil tank in the middle section with the help of a temperature control system. The middle section of the wind tunnel is outfitted with sensor technologies for the measurement of temperature and relative humidity. The downstream portion of the wind tunnel contains an in-line duct fan and damper for controlling wind velocity which is monitored using a pitot-static tube.
The applicability of the soil tank-wind tunnel apparatus described above was demonstrated in an experimental case study of the effects of wind speed on evaporation rate. Results show that increasing wind speed leads to an increased evaporation rate and shortened Stage I evaporation duration. Increasing wind speed beyond 3 m/sec however, shows little additional impact on Stage I evaporation. Stage II evaporation, governed primarily by properties of the porous medium, appears to be independent of or only slightly influenced by wind speed.
This experimental protocol is applicable to a variety of environmental conditions to include changes in soil conditions (i.e. different soils, packing configurations, vegetation, and urban environments), climate boundary conditions (temperature, wind speed, precipitation) or subsurface conditions (e.g. varying water table levels). The dimensions and sensor layout of the described apparatus can be modified to address the needs of different experiments. The packing procedure described above can similarly be modified to account for different packing configurations such as varying porosity conditions and soil heterogeneity.

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&#8217;s surface. Many weather and climate phenomena (e.g., hurricanes, El Ni&#38;#241;o, droughts, etc.) are principally driven by processes associated with atmospheric-land surface interactions1. As more than half of the land surface on the Earth is arid or semiarid2-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 interact5.
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 &#8220;lumped parameterizations&#8221; by performing well controlled experiments under transient conditions that are capable of testing heat and water transfer theory6.
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 sampling7-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 processes11.
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).

<table><tbody><tr><td>ECH2O EC-5 Soil Moisture Sensor (25)</td><td>Decagon Devices Inc. Decagon.com</td><td>40593</td><td>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 &amp;plusmn;3%, and collect data using the Em50 dataloggers</td></tr><tr><td>ECT Soil/Air Temperature Sensor (19)</td><td>Decagon Devices Inc. Decagon.com</td><td>40651</td><td>For specifications visit http://www.decagon.com/products/canopy-atmosphere/temperature/ect-air-temperature/. Sampling frequency on 10 min intervals, accuracy is &amp;plusmn;0.5&amp;nbsp;&amp;deg;C, Measure within a temperature of 5 and 40&amp;nbsp;&amp;deg;C, and collect data using the Em50 dataloggers</td></tr><tr><td>EHT Relative Humidity and Temperature Sensor (5)</td><td>Decagon Devices Inc. Decagon.com</td><td>N/A</td><td>Sampling Frequency on 10 min intervals, accuracy is &amp;plusmn;3% between 5% and 100% relative humidity, and collect data using Em50 data loggers. For more information visit decagon.com</td></tr><tr><td>Em50 Data Logger (10)</td><td>Decagon Devices Inc. Decagon.com</td><td>40800</td><td>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</td></tr><tr><td>Sartorius Weighing Scale (1)</td><td>Sartorius Corporation</td><td>11209-95</td><td>Sartorius Model 11209-95, Range = 65&amp;nbsp;kg, Resolution = &amp;plusmn;1&amp;nbsp;g</td></tr><tr><td>Infrared SalamandernCeramic Radiative Heater (1)</td><td>Mor Electric Heating Assoc., Inc. http://www.morelectricheating.com/</td><td>FTE 500-240</td><td>5 heaters needed, adjust to get the right ambient/free-flow temperature</td></tr><tr><td>2104 Temperature Control System (1)</td><td>Chromalox</td><td>2104</td><td>Controls the heaters
&lt;br /&gt;&amp;#160;&lt;br /&gt;&amp;#160;&lt;br /&gt;&amp;#160;&lt;br /&gt;&amp;#160;&lt;br /&gt;&amp;#160;&lt;br /&gt;&amp;#160;</td></tr><tr><td>Infrared Temperature Sensor Regulator (1)</td><td>Exergen Corporation</td><td>N/A</td><td>Monitors the heaters temperatures</td></tr><tr><td colspan="4">[header]</td></tr><tr><td>Stainless Steel Pitot-Static Tube (1)</td><td>Dwyer Instruments, Inc. http://www.dwyer-inst.com/</td><td>Series 160</td><td>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&#039;s equation.</td></tr><tr><td>1/2 inch Acrylic (1)</td><td>Colorado Plastics http://www.coloradoplastics.com/</td><td>N/A</td><td>Specific heat of 1,464&amp;nbsp;J&amp;nbsp;kg&lt;sup&gt;-1&lt;/sup&gt;&amp;nbsp;K&lt;sup&gt;-1&lt;/sup&gt;, thermal conductivity of 0.2&amp;nbsp;W&amp;nbsp;m&lt;sup&gt;-1&lt;/sup&gt;&amp;nbsp;K&lt;sup&gt;-1&lt;/sup&gt;, and a density of 1,150&amp;nbsp;kg&amp;nbsp;m&lt;sup&gt;-3&lt;/sup&gt;</td></tr><tr><td>Galvanized Steel Ducting Material (1)</td><td>Home Depot</td><td>N/A</td><td>Material used to build wind- tunnel, and both round and rectangular ducting were used in construction and connected using square-to-round reducer duct</td></tr><tr><td>Variable Speed Controller Connected to an In-Line Duct Fan (1)</td><td>Suncourt, Inc. http://www.suncourt.com/</td><td>VS200</td><td>15.3 cm in Diameter Placed in-line with round duct</td></tr><tr><td>Galvanized Steel Damper (1)</td><td>Home Depot</td><td>N/A</td><td>Used to control/reduce speeds in the wind tunnel for low velocity data</td></tr><tr><td>Accusand #30/40 (1)</td><td>Unimin Corporation http://www.unimin.com/</td><td>N/A</td><td>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&amp;nbsp;g/cm&lt;sup&gt;3&lt;/sup&gt;.</td></tr></tbody></table>

exploring the effects of atmospheric forcings on evaporation: experimental integration of the atmospheric boundary layer and shallow subsurface

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.

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 presented here are for a wind speed of 1.22 m/sec. Cumulative evaporation data is shown for all four experiments.
<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr>
			<td rowspan="3">Packing Conditions</td>
			<td rowspan="3">Dry Bulk Density 
			(g cm&#8211;3)</td>
			<td rowspan="3">Air Entry Pressure 
			(cm H2O)</td>
			<td rowspan="3">Residual Water Content 
			(m3 m&#8211;3)</td>
			<td colspan="2">Van Genuchten</td>
		</tr>
		<tr>
			<td colspan="2">Model Parameters *</td>
		</tr>
		<tr>
			<td>&#945; (cm&#8211;1)</td>
			<td>n (&#8211;)</td>
		</tr>
		<tr>
			<td>Tight</td>
			<td>1.79</td>
			<td>16.1/22.5</td>
			<td>0.028</td>
			<td>0.04</td>
			<td>20.53</td>
		</tr>
	</tbody>
</table>
Table 2: Key properties of experimental test soil used.
<table border="0" cellpadding="0" cellspacing="0">
	<tbody>
		<tr>
			<td rowspan="2">Experiment Run #</td>
			<td>Average Maximum Wind Velocity</td>
			<td>Initial Temperature 
			on Soil surface</td>
			<td>Final Temperature 
			on Soil Surface</td>
		</tr>
		<tr>
			<td>(m/sec)</td>
			<td>(&#176;C)</td>
			<td>(&#176;C)</td>
		</tr>
		<tr>
			<td>1</td>
			<td>0.55</td>
			<td>27</td>
			<td>31</td>
		</tr>
		<tr>
			<td>2</td>
			<td>1.22</td>
			<td>26</td>
			<td>33</td>
		</tr>
		<tr>
			<td>3</td>
			<td>3</td>
			<td>29</td>
			<td>37</td>
		</tr>
		<tr>
			<td>4</td>
			<td>3.65</td>
			<td>33</td>
			<td>44.5</td>
		</tr>
	</tbody>
</table>
Table 3: Experimental wind velocities applied.
Time-dependent relative humidity and temperature measured at the soil surface are presented in Figure 3. The relative humidity remains relatively constant at around 0.80 for approximately two days before steeply decreasing over the next four days, beyond which a stable relative humidity value of 0.35 is obtained. The temperature of the soil surface shows an increasing trend over a four-day period before stabilizing. These trends were observed in all four experiments and can be explained in terms of the soil drying. Relative humidity decreases in conjunction with a decrease in evaporation rate because there is less water vapor present over time. The temperature increases as the available water decreases (i.e., evaporation rate decreases) because the process of evaporation no longer cools the soil surface. During the first three days, the relative humidity of downstream air was higher than upstream air due to the presence of more water vapor resulting from upstream evaporation. This trend was reversed later, most likely due to the upstream sensor losing contact with the soil surface; the sensor cables are flexible and occasionally pull the sensor from the soil surface, changing the humidity reading. The relative humidity measured downstream is greater than that measured upstream because the process of evaporation along the first 21.5 cm of the tank increased the amount of moisture present in air.
<img alt="Figure 3" src="/files/ftp_upload/52704/52704fig3highres.jpg" width="700" /> 
Figure 3: Relative humidity and temperature measured on the soil surface (This figure has been modified from Davarzani et al.5).
The free flow air temperature in this experiment was set to a constant value of 40 &#176;C using the previously described temperature control system. The time-dependent temperature and relative humidity of the air in the free flow, at a height of 8.5 cm above the soil surface, are shown in Figure 4. The observed diurnal fluctuations in temperature are due to the variability of heater outputs in response to the temperature measured by the infrared temperature sensor that regulates the temperature control system (see Step 2.3). Diurnal fluctuations can be avoided, if desired, by setting the infrared temperature sensor to a set temperature value. The difference in atmospheric temperature along the length of the tank is the result of evaporative cooling (Figure 4).
<img alt="Figure 4" src="/files/ftp_upload/52704/52704fig4highres.jpg" width="700" /> 
Figure 4: Relative humidity and temperature measured 7.5 cm above the soil surface upstream and downstream of the tank (This figure has been modified from Davarzani et al.5).
In Figure 5a, the time dependent soil temperature is shown for the depths of 2.5 cm, 7.5 cm and 12.5 cm below the soil surface as well as the ambient temperature; see Figure 1 for sensor identification. As seen in Figure 5a, surface temperature and wind velocity are less influential on local temperatures at greater depths &#8212; showing no effect at depths below 12.5 cm. Figure 5b shows temperature as a function of time for three sensors located at a depth of 2.5 cm. There is a slight difference in temperature for sensors at this depth with the upstream sensor 5 showing a higher temperature than the downstream Sensor 1. This is because the free flow temperature is always higher upstream than downstream (Figure 4). The differences in temperature also results in an asymmetrical saturation profile in the soil tank as will be subsequently shown.
A 
<img alt="Figure 5a" src="/files/ftp_upload/52704/52704fig5ahighres.jpg" width="700" /> 
B 
<img alt="Figure 5b" src="/files/ftp_upload/52704/52704fig5bhighres.jpg" width="700" /> 
Figure 5: Evolution of measured soil temperature as a function of time (a) vertically at the middle of the tank and (b) horizontally at a depth of 2.5 cm (This figure has been modified from Davarzani et al.5).
Figure 6a shows the time-dependent saturation versus time at soil depths of 2.5, 7.5, 12.5, and 17.5 cm. For depths greater than 12.5 cm, the saturation remained at 100% for the duration of the experiment; closer to the soil surface however, saturation decreased over time. The saturation shown in Figure 6a can be related to the different stages of evaporation (i.e. Stage I and Stage II), defined by differences in evaporation rates, location of the drying front, and dominant transport mechanisms14. During Stage I evaporation, the drying front rapidly retreats away from the soil surface as gravitational and viscous forces begin to dominate capillary forces. This is observed in the first day by the measured decrease in soil saturation by the first row of soil moisture sensors corresponding to a depth of 2.5 cm. After Day 1, the rate at which the drying front continues to retreat slows as shown in the gradual shape of the saturation curves for sensors 6-10 located at a depth of 7.5 cm (Figure 6a). This marks the transition of evaporation to vapor diffusion limited Stage II evaporation. The initial part of the Stage II is often called the falling rate period15-17. Eventually, the saturation curves level out and change very little as the drying front reaches a depth of 12.5 cm (e.g., Sensor 13) by Day 3.
A 
<img alt="Figure 6a" src="/files/ftp_upload/52704/52704fig6ahighres.jpg" width="700" /> 
B 
<img alt="Figure 6b" src="/files/ftp_upload/52704/52704fig6bhighres.jpg" width="700" /> 
Figure 6: Time evolution of measured subsurface soil saturation (a) vertically at the middle of the tank and (b) horizontally at a depth of 2.5 cm (This figure has been modified from Davarzani et al.5).
Figure 6b shows saturation versus time for three sensors located at constant depth of 2.5 cm. The saturation curves are nearly identical and consistent across the entire length of the tank at this depth. The slight asymmetric distribution is due to the air temperature difference between the upstream and downstream sections of the wind tunnel. Since upstream temperatures were consistently a few degrees warmer, atmospheric demand, which drives evaporation, would be higher and therefore there would be a slightly faster rate of drying.
Figure 7 shows wind speed, mean value of 1.22 m/sec, as a function of time. The observed sinusoidal diurnal trend in the wind speed is the result of changes in atmospheric conditions such as barometric pressure and air density. The average wind speed was used in modeling efforts because the effects of diurnal fluctuations of atmospheric variables were not the focus of the present study. This does not mean however, that time-dependent data could not be used. As part of the series of evaporation experiments, four different mean wind velocities were applied; see Table 3 for a summary. The calculated Reynolds numbers for all experiments in this study were within the laminar and transitional flow regimes. However, it is well know that surface turbulence can affect the evaporation rates16 and should be addressed in future studies.
<img alt="Figure 7" src="/files/ftp_upload/52704/52704fig7highres.jpg" width="700" /> 
Figure 7: Time-dependent wind speed over the soil surface with mean value of 1.22 m/sec&#8211;1 (This figure has been modified from Davarzani et al.5).
The effect of the airflow in the free fluid region (i.e. atmosphere) on cumulative evaporation is shown in Figure 8. Cumulative evaporation is plotted for four different free flow average wind speeds (Vw) of 0.50, 1.20, 3.00 and 3.60 m/sec. Results demonstrate that wind speed has a very prominent effect on cumulative evaporation and the amount of water loss during the different evaporative stages. As shown in Figure 8, increasing the wind speed increases the total evaporation. By comparing the slopes of the curves, the greatest influence was on the initial evaporation rate, here referred to as stage 1. Stage 1 evaporation is often defined by high and relatively constant evaporation rates17 and is predominantly influenced by atmospheric demand rather than soil conditions. As wind speed is further increased from 3 to 3.6 m/sec, evaporation shows much less dependency on incremental changes in wind speed than was observed for changes at low wind speeds. Increasing wind speed leads to an increase in Stage I evaporation rate while simultaneously decreasing the transition time from Stage I to Stage II5. The influence of wind speed on evaporation is less significant for Stage II evaporation which is controlled predominantly by the porous medium. During this stage, evaporation is controlled by the rate at which water can be transmitted to the soil surface via diffusion rather than atmospheric demand.
<img alt="Figure 8" src="/files/ftp_upload/52704/52704fig8highres.jpg" width="700" /> 
Figure 8: The effect of different mean wind speeds on cumulative evaporation (This figure has been modified from Davarzani et al.5).

Watch this Scientific Journal Video about Exploring the Effects of Atmospheric Forcings on Evaporation: Experimental Integration of the Atmospheric Boundary Layer and Shallow Subsurface at JoVE.com

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

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.

Evaporation is directly influenced by the interactions between the atmosphere, land surface and soil subsurface. This work aims to experimentally study ...

exploring-effects-atmospheric-forcings-on-evaporation-experimental

Nanyang Technological University

Research

Education

JoVE Core

JoVE Journal

Civil Engineering

Environment

Unit 1

Leveling

SCIENTISTS IN ACTION

8.7K Views.  Colorado School of Mines. 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.

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

Surface Renewal: An Advanced Micrometeorological Method for Measuring and Processing Field-Scale Energy Flux Density Data

Advanced micrometeorological methods have become increasingly important in soil, crop, and environmental sciences. For many scientists without formal training in atmospheric science, these techniques are relatively inaccessible. Surface renewal and other flux measurement methods require an understanding of boundary layer meteorology and extensive training in instrumentation and multiple data management programs. To improve accessibility of these techniques, we describe the underlying theory of surface renewal measurements, demonstrate how to set up a field station for surface renewal with eddy covariance calibration, and utilize our open-source turnkey data logger program to perform flux data acquisition and processing. The new turnkey program returns to the user a simple data table with the corrected fluxes and quality control parameters, and eliminates the need for researchers to shuttle between multiple processing programs to obtain the final flux data. An example of data generated from these measurements demonstrates how crop water use is measured with this technique. The output information is useful to growers for making irrigation decisions in a variety of agricultural ecosystems. These stations are currently deployed in numerous field experiments by researchers in our group and the California Department of Water Resources in the following crops: rice, wine and raisin grape vineyards, alfalfa, almond, walnut, peach, lemon, avocado, and corn.

Surface renewal is a micrometeorological method that is being used increasingly to determine energy fluxes, but its technical complexity makes it inaccessible to a broad audience. We describe the steps needed to set up and calibrate a surface renewal field station, to acquire and process data, and to correctly interpret results.

Advanced micrometeorological methods have become  increasingly important in soil, crop, and environmental sciences.  For many scientists without formal ...

Evolution of Staircase Structures in Diffusive Convection

Diffusive convection (DC) occurs when the vertical stratified density is controlled by two opposing scalar gradients that have distinctly different molecular diffusivities, and the larger- and smaller- diffusivity scalar gradients have negative and positive contributions for the density distribution, respectively. The DC occurs in many natural processes and engineering applications, for example, oceanography, astrophysics and metallurgy. In oceans, one of the most remarkable features of DC is that the vertical temperature and salinity profiles are staircase-like structure, composed of consecutive steps with thick homogeneous convecting layers and relatively thin and high-gradient interfaces. The DC staircases have been observed in many oceans, especially in the Arctic and Antarctic Oceans, and play an important role on the ocean circulation and climatic change. In the Arctic Ocean, there exist basin-wide and persistent DC staircases in the upper and deep oceans. The DC process has an important effect on diapycnal mixing in the upper ocean and may significantly influence the surface ice-melting. Compared to the limitations of field observations, laboratory experiment shows its unique advantage to effectively examine the dynamic and thermodynamic processes in DC, because the boundary conditions and the controlled parameters can be strictly adjusted. Here, a detailed protocol is described to simulate the evolution process of DC staircase structure, including its generation, development and disappearance, in a rectangular tank filled with stratified saline water. The experimental setup, evolution process, data analysis, and discussion of results are described in detail.

Diffusive convection (DC) widely occurs in natural processes and engineering applications, characterized by a series of staircases with homogeneous convecting layers and stratified interfaces. An experimental procedure is described to simulate the evolution process of the DC staircase structure, including the generation, development and disappearance, in a rectangular tank.

Diffusive convection (DC) occurs when the vertical stratified density is controlled by two opposing scalar gradients that have distinctly different ...

Measurements of CO2 Fluxes at Non-Ideal Eddy Covariance Sites

This protocol is an example of utilizing the eddy covariance (EC) technique to investigate spatially and temporally averaged net CO2 fluxes (net ecosystem production, NEP), in non-typical ecosystems, on a currently reforested windthrow area in Poland. After a tornado event, a relatively narrow &#8220;corridor&#8221; was created within surviving forest stands, which complicates such kind of experiments. The application of other measuring techniques, such as the chamber method, is even more difficult under these circumstances, because especially at the beginning, fallen trees and in general great heterogeneity of the site provide a challenging platform to perform flux measurements and then to properly upscale obtained results. In comparison with standard EC measurements carried out in untouched forests, the case of windthrow areas requires special consideration when it comes to the site location and data analysis in order to ensure their representativeness. Therefore, here we present a protocol of real-time, continuous CO2 flux measurements at a dynamically changing, non-ideal EC site, which includes (1) site location and instrumentation setup, (2) flux computation, (3) rigorous data filtering and quality control, and (4) gap filling and net fluxes partitioning into CO2 respiration and absorption. The main advantage of the described methodology is that it provides a detailed description of the experimental setup and measurement performance from scratch, which can be applied to other spatially limited ecosystems. It can also be viewed as a list of recommendations on how to deal with unconventional site operation, providing a description for non-specialists. Obtained quality-checked, gap filled, half-hour values of net CO2, as well as absorption and respiration fluxes, can be finally aggregated into daily, monthly, seasonal or annual totals.

Measurements of CO2 Fluxes at Non-Ideal Eddy Covariance Sites

The presented protocol uses the eddy covariance method at non-typical locations, applicable to all types of short-canopy ecosystems with limited area, on a currently reforested windthrow site in Poland. Details of measuring site setup rules, flux calculations and quality control, and final result analysis, are described.

This protocol is an example of utilizing the eddy covariance (EC) technique to investigate spatially and temporally averaged net CO2 fluxes (net ecosystem ...

Protocol for Soil Hydraulic Property Measurement via Evaporation

Measurements of Soil Water Potential and Conductivity based on a Simple Evaporation Experiment using a Hydraulic Property Analyzer

The measurement of soil hydraulic properties is critical in understanding the physical components of soil health as well as integrated knowledge of soil systems under various management practices. Collecting reliable data is imperative for informing decisions that affect agriculture and the environment. The simple evaporation experiment described here uses instrumentation in a laboratory setting to analyze soil samples collected in the field. The soil water tension of the sample is measured by the instrument, and tension data is modeled by software to return soil hydraulic properties. This method can be utilized to measure soil water retention and hydraulic conductivity and give insight into differences in treatments or environmental dynamics over time. Initial establishment requires a user, but data acquisition is automated with the instrument. Soil hydraulic properties are not easily measured with traditional experiments, and this protocol offers a simple and optimal alternative. Interpretation of results and options for extending the data range are discussed.

Author Spotlight: Advancing Agricultural Land Ecosystem Research with a Hydraulic Property Analyzer to Assess Soil Health

This article features a simple evaporation experiment using a hydraulic property instrument for a soil sample. Through efficient means, measurements can be taken over a series of days to generate high-quality data.

The measurement of soil hydraulic properties is critical in understanding the physical components of soil health as well as integrated knowledge of soil ...

Measurement of Leaf Hydraulic Conductance and Stomatal Conductance and Their Responses to Irradiance and Dehydration Using the Evaporative Flux Method (EFM)

Water is a key resource, and the plant water transport system sets limits on maximum growth and drought tolerance. When plants open their stomata to achieve a high stomatal conductance (gs) to capture CO2 for photosynthesis, water is lost by transpiration1,2. Water evaporating from the airspaces is replaced from cell walls, in turn drawing water from the xylem of leaf veins, in turn drawing from xylem in the stems and roots. As water is pulled through the system, it experiences hydraulic resistance, creating tension throughout the system and a low leaf water potential (&Psi;leaf). The leaf itself is a critical bottleneck in the whole plant system, accounting for on average 30% of the plant hydraulic resistance3. Leaf hydraulic conductance (Kleaf = 1/ leaf hydraulic resistance) is the ratio of the water flow rate to the water potential gradient across the leaf, and summarizes the behavior of a complex system: water moves through the petiole and through several orders of veins, exits into the bundle sheath and passes through or around mesophyll cells before evaporating into the airspace and being transpired from the stomata. Kleaf is of strong interest as an important physiological trait to compare species, quantifying the effectiveness of the leaf structure and physiology for water transport, and a key variable to investigate for its relationship to variation in structure (e.g., in leaf venation architecture) and its impacts on photosynthetic gas exchange. Further, Kleaf responds strongly to the internal and external leaf environment3. Kleaf can increase dramatically with irradiance apparently due to changes in the expression and activation of aquaporins, the proteins involved in water transport through membranes4, and Kleaf declines strongly during drought, due to cavitation and/or collapse of xylem conduits, and/or loss of permeability in the extra-xylem tissues due to mesophyll and bundle sheath cell shrinkage or aquaporin deactivation5-10. Because Kleaf can constrain gs and photosynthetic rate across species in well watered conditions and during drought, and thus limit whole-plant performance they may possibly determine species distributions especially as droughts increase in frequency and severity11-14.
We present a simple method for simultaneous determination of Kleaf and gs on excised leaves. A transpiring leaf is connected by its petiole to tubing running to a water source on a balance. The loss of water from the balance is recorded to calculate the flow rate through the leaf. When steady state transpiration (E, mmol &bull; m-2 &bull; s-1) is reached, gs is determined by dividing by vapor pressure deficit, and Kleaf by dividing by the water potential driving force determined using a pressure chamber (Kleaf= E /- &Delta;&psi;leaf, MPa)15.
This method can be used to assess Kleaf responses to different irradiances and the vulnerability of Kleaf to dehydration14,16,17.

Measurement of Leaf Hydraulic Conductance and Stomatal Conductance and Their Responses to Irradiance and Dehydration Using the Evaporative Flux Method EFM

We describe a relatively rapid (30 min) and realistic method for simultaneously measurement of leaf hydraulic conductance (Kleaf) and stomatal conductance (gs) for transpiring excised leaves. The method can be modified to measure the light and dehydration responses of Kleaf and gs.

Water is a key resource, and the plant water transport system sets limits on maximum growth and drought tolerance. When plants open their stomata to ...

Biology

Boundary Layer Characteristics

When a fluid encounters a solid surface, a boundary layer forms due to the interaction between the fluid's motion and the stationary surface. This phenomenon is characterized by a thin region adjacent to the surface where viscous forces dominate, influencing the fluid's velocity profile. The development of the boundary layer begins at the leading edge of the surface and evolves as the fluid moves downstream.As the fluid flows over the surface, friction between the fluid and the wall slows down the particles nearest the surface. This effect creates a velocity gradient, where the velocity of the fluid increases from zero at the wall to the free-stream velocity farther away. The region where this velocity gradient exists defines the boundary layer. Over distance, the thickness of the boundary layer increases as more fluid is affected by the surface's resistance.In the region close to the leading edge, the boundary layer is typically laminar. Here, fluid particles move in smooth, parallel layers with minimal mixing. The flow remains stable under certain conditions, but as it progresses downstream, factors such as surface roughness, velocity increases, or pressure gradients can destabilize the flow. Once the flow reaches a critical Reynolds number, the laminar boundary layer transitions to a turbulent state. The turbulent boundary layer is characterized by chaotic mixing of fluid particles, increased energy dissipation, and a thicker boundary layer compared to its laminar counterpart.The laminar boundary layer's velocity profile can be described mathematically using the Blasius equation, a solution to the Navier-Stokes equations under specific assumptions. For practical applications, the Momentum Integral Boundary Layer equation provides a way to analyze momentum changes in the boundary layer without requiring a detailed velocity profile. These models are critical for understanding fluid dynamics near surfaces, particularly in engineering and aerodynamic applications.

When a fluid encounters a solid surface, a boundary layer forms due to the interaction between the fluid's motion and the stationary surface. This ...

Flow over Immersed Bodies

KEY TERMS AND CONCEPTS

Influence of Earth's Curvature and Atmospheric Refraction on Leveling

During leveling, the Earth's curvature and atmospheric refraction introduce deviations in the line of sight from a true horizontal reference.&nbsp;When the line of sight is leveled, it remains perpendicular to the plumb line only at a single point. Beyond this, it deviates due to the Earth’s curvature, represented by the correction C. For a sight distance D, the deviation can be derived using the relationship:<figure class="image image-style-align-center image_resized" style="display:table;margin-left:auto;margin-right:auto;width:9.53%;"><img style="aspect-ratio:197/150;" src="https://cloudfront-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/18021/18021_Equation_1_editor_18021.svg" alt="18021_Equation_1.svg"></figure>This relationship shows that the deviation increases quadratically with distance. Over a horizontal distance of 1 kilometer, it is approximately 0.0785 meters.In addition to curvature, atmospheric refraction bends the line of sight downward, further deviating it from the true horizontal. The bending effect is approximately one-seventh of the Earth's curvature deviation, amounting to about 0.0112 meters over 1 kilometer. This refraction arises from variations in the refractive index of air with altitude, causing light rays to curve slightly.&nbsp;The combined influence of Earth's curvature and atmospheric refraction is approximately 0.0675 meters over 1 kilometer. As with curvature alone, this combined effect increases with the square of the distance, which requires careful adjustments during leveling operations to maintain accuracy. Correcting for these deviations is critical for high-precision applications, such as determining benchmark elevations and establishing horizontal planes over long distances.

During leveling, the Earth's curvature and atmospheric refraction introduce deviations in the line of sight from a true horizontal reference.&nbsp;When ...

Variation of Atmospheric Pressure

Change in atmospheric pressure with height is particularly interesting. The decrease in atmospheric pressure with increasing altitude is due to the decreasing gravitational force per unit area as we move away from the surface of the earth.
Assuming the air temperature is constant at a given altitude and that the ideal gas law of thermodynamics describes the atmosphere to a good approximation, one can find the variation of atmospheric pressure with height.
Let p(y) be the atmospheric pressure at height y. The density &#961; at y, the temperature T in the Kelvin scale (K), and the mass m of a molecule of air are related to the absolute pressure by the ideal gas law, in the form:
<img alt="Equation1" src="/files/ftp_upload/14711/14711_Equation_1.svg" style="text-align: center; height: 38px;" />
Using density from the ideal gas law, the rate of variation of pressure with height is integrated from sea level, and the final expression is obtained as:
<img alt="Equation2" src="/files/ftp_upload/14711/14711_Equation_2.svg" style="text-align: center; height: 17px;" />
where,
<img alt="Equation3" src="/files/ftp_upload/14711/14711_Equation_3.svg" style="text-align: center; height: 38px;" />
Atmospheric pressure drops exponentially with height, since the y-axis points up from the ground, and y has positive values in the atmosphere above sea level. The pressure drops by a factor of 1/e when the height is 1/&#945;, which gives us a physical interpretation for &#945;. The constant 1/&#945; is a length scale that characterizes how pressure varies with height and is often referred to as the pressure scale height.

Change in atmospheric pressure with height is particularly interesting. The decrease in atmospheric pressure with increasing altitude is due to the ...

/files/media/science-education/science-education-thumbs/14711.jpg

Physics

Fluid Mechanics

Vapor Pressure of Fluid

The vapor pressure of a fluid is a crucial concept in fluid mechanics, influencing phenomena such as boiling and cavitation. Vapor pressure refers to the pressure exerted by a vapor at a state of thermodynamic equilibrium with its corresponding liquid phase at a specific temperature. It represents the tendency of molecules to escape from the fluid surface into the vapor phase.
When a liquid is placed in a closed container with a small air space, and the space is evacuated, vapor molecules will escape from the liquid into the space. As more molecules escape, the pressure in the space increases until an equilibrium is reached, where the rate of molecules escaping the liquid equals the rate of molecules returning to it. This equilibrium pressure is the vapor pressure of the liquid.
Vapor pressure is temperature-dependent. As the temperature increases, the molecules&#39; kinetic energy increases, resulting in more molecules overcoming intermolecular forces and escaping into the vapor phase. It follows that the vapor pressure increases with temperature. For instance, water at 100&#176;C has a vapor pressure of 101.3 kPa (1 atm), allowing it to boil at standard atmospheric pressure. However, at higher elevations, where atmospheric pressure is lower, water boils at lower temperatures because its vapor pressure equals the reduced atmospheric pressure at the higher elevation.

The vapor pressure, which influences boiling and cavitation, is the pressure exerted by a vapor in equilibrium with its liquid at a given temperature and ...

Fluid Properties

Vapor Pressure Lowering

The equilibrium vapor pressure of a liquid is the pressure exerted by its gaseous phase when vaporization and condensation are occurring at equal rates:
&#160;<img alt="Eq1" src="/files/ftp_upload/11367/11367_Equation_1.png" style="height: 50px;" />
Dissolving a nonvolatile substance in volatile liquid results in a lowering of the liquid&#8217;s vapor pressure. This phenomenon can be explained by considering the effect of added solute molecules on the liquid&#39;s vaporization and condensation processes. To vaporize, solvent molecules must be present at the surface of the solution. The presence of solute decreases the surface area available to solvent molecules and thereby reduces the rate of solvent vaporization. Since the rate of condensation is unaffected by the presence of solute, the net result is that the vaporization-condensation equilibrium is achieved with fewer solvent molecules in the vapor phase (i.e., at a lower vapor pressure).
While this interpretation is useful, it does not account for several important aspects of the colligative nature of vapor pressure lowering. A more rigorous explanation involves the property of entropy. For purposes of understanding the lowering of a liquid&#39;s vapor pressure, it is adequate to note that the more dispersed nature of matter in a solution, compared to separate solvent and solute phases, serves to effectively stabilize the solvent molecules and hinder their vaporization. A lower vapor pressure results, and a correspondingly higher boiling point.
The relationship between the vapor pressures of solution components and the concentrations of those components is described by Raoult&#8217;s law: The partial pressure exerted by any component of an ideal solution is equal to the vapor pressure of the pure component multiplied by its mole fraction in the solution.
&#160;<img alt="Eq2" src="/files/ftp_upload/11367/11367_Equation_2.png" style="height: 50px;" />
where PA is the partial pressure exerted by component A in the solution, P&#186;A is the vapor pressure of pure A, and XA is the mole fraction of A in the solution.
Recalling that the total pressure of a gaseous mixture is equal to the sum of partial pressures for all its components (Dalton&#8217;s law of partial pressures), the total vapor pressure exerted by a solution containing i components is
&#160;<img alt="Eq3" src="/files/ftp_upload/11367/11367_Equation_3.png" style="height: 50px;" />
A nonvolatile substance is one whose vapor pressure is negligible (P&#186; &#8776; 0), and so the vapor pressure above a solution containing only nonvolatile solutes is due only to the solvent:
&#160;<img alt="Eq4" src="/files/ftp_upload/11367/11367_Equation_4.png" style="height: 50px;" />
This text is adapted from <a href="https://openstax.org/books/chemistry-2e/pages/11-4-colligative-properties?query=vapor%20pressure&amp;target=%7B%22index%22%3A0%2C%22type%22%3A%22search%22%7D#fs-idp74224880">Openstax, Chemistry 2e, Section 11.4: Colligative Properties.</a>

Colligative Properties; Raoult&#039;s Law &amp; Vapor Pressure Lowering

The equilibrium vapor pressure of a liquid is the pressure exerted by its gaseous phase when vaporization and condensation are occurring at equal rates:
...