Next Article in Journal
The Impact of Land Cover on Selected Water Quality Parameters in Polish Lowland Streams during the Non-Vegetative Period
Next Article in Special Issue
The Influence of pH on Subsurface Denitrification Stimulated with Emulsified Vegetable Oil
Previous Article in Journal
Characterization of the Propagation of Meteorological Drought Using the Copula Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Evaporation and Seasonal Effects on the Water Content in the Unsaturated Zone: A Multi-Year Laboratory Study

by
Aderonke O. Adegbule
1,2 and
Tohren C. G. Kibbey
1,*
1
School of Civil Engineering and Environmental Science, University of Oklahoma, Norman, OK 73019, USA
2
Office of Chemical Safety and Pollution Prevention, Environmental Protection Agency, Washington, DC 20460, USA
*
Author to whom correspondence should be addressed.
Water 2022, 14(20), 3294; https://doi.org/10.3390/w14203294
Submission received: 30 September 2022 / Revised: 14 October 2022 / Accepted: 16 October 2022 / Published: 19 October 2022

Abstract

:
The unsaturated zone of the subsurface plays a critical role in the environmental fate and transport of a wide range of environmental contaminants, and is also important for plant growth and agriculture. Quantitative prediction of processes in the unsaturated zone requires knowledge of how water content varies with elevation above the water table, a relationship known as the capillary pressure (Pc)–saturation (S) relationship. While the PcS relationship is conventionally thought of as being primarily a property of the porous medium and fluids, previous work found evidence suggesting that it actually results from a dynamic equilibrium between capillary forces and phenomena driven by evaporation. The focus of this work was on gaining a further understanding of the role of evaporation on the PcS relationship. The work made use of a tall instrumented laboratory column connected to an external reservoir for maintaining water table height. Following an initial imbibition experiment, the column was saturated and allowed to drain, and then water content was monitored in the column as a function of height over 1207 days (3.31 years). While initial imbibition and drainage were rapid, on the order of hours, redistribution and evaporation effects became apparent over longer time scales (hundreds of hours). A drying front moving downward in the column was apparent from changes in the slopes of saturation vs. time curves as it passed individual sensors; unlike previous experiments with a more fully-vented column, the evaporation front appeared to stall, balanced by the capillary-driven upward flow of water. Over the full multi-year duration of the experiment, seasonal trends in water saturation were apparent, with significant, reversible variations observed that closely followed atmospheric conditions. Specifically, saturation above the water table appeared to increase during the spring and summer months and decrease during the fall and winter months, despite the constant water table location, consistent with a changing driving force for evaporation. This result may suggest the likelihood of seasonal effects in the long-term transport and fate of contaminants in the unsaturated zone.

1. Introduction

The unsaturated zone is a region of the subsurface above the water table. Both air and water are present in pores in the unsaturated zone, and their relative quantities change dynamically over time as a result of external inputs (e.g., infiltration from rainfall) or outputs (evaporation and evapotranspiration), as well as due to spatial redistribution. At equilibrium, the relationship between the water content and elevation is described by the capillary pressure–saturation (PcS) relationship. In the absence of gas pressure gradients in a porous medium, capillary pressure (Pc) is equal to the negative water pressure; if vertical flow does not occur, Pc in height units is the same as the elevation above the water table. Saturation (S), a measure of water content, is the ratio of the volume of water in a porous medium to the volume of pores; a saturation of 1.0 indicates a fully saturated medium, while a saturation of 0.0 is a fully dry porous medium. The PcS relationship is a function of the porous medium, with finer-grained porous media exhibiting greater S at a given Pc (or elevation). Furthermore, the PcS relationship is hysteretic, meaning that the observed S at any given point depends not only on the Pc, but also on the wetting/drying history. In general, when water drains from a previously saturated porous medium, the S is higher at a given Pc (or elevation) than would be the case when water imbibes into a previously dry porous medium.
The capillary fringe is a region of the unsaturated zone just above the water table where the porous medium is near full saturation. Finer-grained porous media have thicker capillary fringes, while coarser-grained porous media have thinner capillary fringes. The ability to predict and control the water content in the unsaturated zone and the capillary fringe is of critical importance for many practical applications, including agriculture and plant growth [1,2], as well as the movement of solutes and colloids in the environment [3,4,5,6,7,8,9]. Understanding and simulating water movement during these applications requires knowledge of the PcS relationship for the porous medium of interest. While methods of measuring or estimating PcS parameters for a particular porous medium are available (e.g., [10,11,12]), little is known about how environmental effects, including evaporation, influence the PcS relationship. While the dynamics of evaporation and its effects on water content near the ground surface are of great interest and have been widely studied (e.g., [13,14,15,16,17]), less information is available about its effect on the PcS relationship itself.
Previous work by the authors explored the hypothesis that the equilibrium PcS relationship is the result of a dynamic equilibrium between the various factors that influence water content in unsaturated media, including capillary forces, evaporative water loss, and the upward flow of water to replenish evaporated water [18]. In that study, the PcS relationship was studied in both a coarse sand and a thin nylon membrane. The sand was studied in two parallel systems packed between two parallel surfaces, in one case with one surface open to the atmosphere, in the other with both surfaces closed to the atmosphere (in both cases, the top of the porous media was open to the atmosphere). Drainage curves were imaged over 117 days, and it was observed that evaporation caused the saturation to continually drop over time in the open system, with the observed PcS curves gradually moving towards (although not reaching) the imbibition PcS relationship, and the thickness of the capillary fringe steadily decreasing. In contrast, the PcS relationships in the closed system remained relatively constant over time over much of the system, only exhibiting saturation decreases near the open top of the column, where a drying front induced by evaporation continued to move downward over the 117 days of the experiment. The membrane system studied in that work explored an extreme case where evaporation dominated the observed PcS behavior due to the large external area available for evaporative mass transfer. In that system, decreasing the driving force for evaporation by covering the system (thus increasing relative humidity in the air in contact with the membrane) led to rapid increases in the water content in the membrane.
The purpose of this work was to further explore the effects of evaporation on the moisture content in the unsaturated zone and the capillary fringe, with emphasis on understanding how evaporation might impact water content over a much longer time period in a tall laboratory column. The experiment conducted for the work extended over 1207 days (3.31 years), more than ten times the duration of the previous sand experiment, and was conducted in a conventional laboratory sand column with only minimal venting.

2. Materials and Methods

Figure 1 shows the experimental configuration of the sand column used for this work. The sand column was constructed from clear polyvinyl chloride plastic with an inner diameter of 2.5 cm, and a final packed height of approximately 87 cm. The column was outfitted with 12 custom capacitive water content sensors and four vents based on 1/8″ NPT quick connects for ¼″ tubing. Sand at the bottom of the column was supported by a #80 mesh stainless steel screen, and connected to an external reservoir using the same fittings described in earlier work [4]. Sensors were glued into milled slots at 5 cm spacings starting 30.5 cm from the bottom of the column using RTV silicone. Vents were positioned 2.5 cm below sensors 1, 4, 7 and 10 (as numbered from the bottom). The column was dry-packed to allow spontaneous imbibition measurements to be conducted. Packing was conducted a few cm at a time, and the material was compressed by tapping with a long rod after each sand addition. Sensors were glued in sequentially from the bottom as the level of added sand reached each sensor location; this allowed the sand to be packed without interference from sensors present higher in the column.
Sensors for the work were fabricated as custom printed circuits. Each sensor was double-sided, with 36 interdigitated traces on each side, and approximate dimensions of 2 cm × 2 cm for the active region. Water content was determined by measuring capacitance between the interdigitated traces. Capacitance was measured using an Arduino Uno, by way of a multiplexor to allow sensors to be read sequentially. All sensors were read at 10 s intervals for all but the later portion of the imbibition experiment, during which sensors were read at 325 s intervals; observed noise in the sensor data during imbibition led to reversion to 10 s intervals for the remainder of the experiment to permit greater amounts of averaging smoothing.
Final values were taken from averaging multiple measured points. In the case of imbibition and initial drainage, rolling averages were used, while for presentation of data from the drainage curves over the full 1207-day duration, samples were taken by averaging data over 600 s (60 sample) windows in one-day intervals. Calibration of the sensors was conducted at the completion of the experiment based on an assumed linear response between maximum and minimum signals observed for wet and dry states at each sensor; in-place calibration is necessary because sensor response is influenced slightly by differences in the amount of silicone used for installation, so individual sensor ranges vary slightly. Separate calibration tests with similar sensors in the past found near-linear behavior. Note that drainage PcS curves measured with the sensors in this work are in close agreement with previously measured curves for the same material, supporting the calibration approach.
The water table for all experiments was controlled by maintaining the water level in an external 4 L high-density polyethylene container with a 15 cm inner diameter connected to the bottom of the column by 0.25 inch outer diameter nylon tubing with a ball valve. The reservoir had a screw-cap to minimize evaporation from the reservoir; a 1.5 mm diameter hole in the cap allowed air to enter and leave so the air in the reservoir would remain at atmospheric pressure.
An experiment measuring imbibition from a dry state was conducted by opening the valve and physically raising the reservoir so the water level was below the bottom sensor. After approximately 20 h, the reservoir was raised above the top of the column and held overnight to saturate the column. On October 22, 2018, the reservoir was then again lowered to below the bottom sensor, and the column was allowed to drain; this reservoir position was maintained for 1207 days (3.31 years), until February 10, 2022. Because of the large diameter of the reservoir, the full pore volume of the column above the water table would be equivalent to a 2.4 cm elevation difference on the reservoir itself. As such, even large changes in saturation over much of the column would only be expected to produce changes in reservoir water level on the order of mm. After the first 24 h, a time period that included a large fraction of the observed drainage, the measured water level in the reservoir had increased by less than 1 mm. For the duration of the 1207-day experiment, the water level was maintained within 0.25 cm of its initial position 5.75 cm below the bottom sensor through the addition of water as needed to adjust for evaporative losses through the column. The four column vents were closed for the imbibition experiment and for the first 128 days of the drainage experiment, and were then opened for the remaining 1079 days on day 128 (28 February 2019). Details of water additions are discussed in the Results and Discussion section in the context of saturation results.
Daily outdoor temperature and relative humidity data for Norman, Oklahoma were collected from the Oklahoma Mesonet [19,20]. Absolute humidity values were calculated based on reported relative humidity and calculated saturated water vapor pressure as a function of temperature. Indoor temperatures were monitored at one-hour intervals for the duration of the experiment using an Omega (Norwalk, CT) OM-63 temperature data logger.
The experiments described in this work were conducted using US Silica (Ottawa, IL) F45 sand, a high-purity, foundry-grade medium sand with an average grain size (d50) of 323 µm. Previous work [18] was based on a coarser sand (US Silica F-30, d50 = 590 µm). The use of a finer sand for this work required the use of a taller column, as a result of the thicker capillary fringe for the F45 sand (an air-entry-pressure of ~30 cm for F45 vs. ~12 cm for F30).

3. Results and Discussion

Figure 2 shows elevation vs. saturation for eight different times within the first 12 h following the start of imbibition. It is apparent from the figure that imbibition proceeds rapidly from the dry state, with the water saturation profile approaching its near final configuration after only four hours, and only subtle changes observed after eight hours. The rapid influx of water is not unexpected for imbibition to initially dry sand. Because the column was initially dry, the top of the column (more than about 45 cm above the water table) remains nearly dry after imbibition.
Figure 3 shows the measured saturations at the 12 sensors in the column, illustrating the almost instantaneous saturation increase low in the column, with a slowing saturation increase high in the column as the wetting front moves farther up the column; this is expected behavior. Note that the noise in Figure 3 is due to a low frequency sampling that was used from 45 min to the end of the imbibition experiment, as described previously.
Following the completion of the imbibition experiment, the column was saturated by raising the reservoir above the top of the column, and then a drainage experiment was started by closing the valve, lowering the reservoir to the desired water table level and opening the valve. Figure 4 shows the measured saturations at the 12 sensors in the column over two time scales. Figure 4a corresponds to the first 20 h of the experiment—the same time scale as the imbibition curves in Figure 3—while Figure 4b corresponds to the first 900 h (37.5 days) of the experiment.
From Figure 4a, it is apparent that the onset of drainage is extremely rapid, with the most substantial changes in saturation occurring in the first hour—faster than for imbibition—after which saturation seemingly stabilizes throughout much of the column. It should be noted, however, that sensors near the bottom of the column show values below full saturation at the start of the experiment, possibly due to packing defects near the bottom sensors that allowed water to bypass pockets of porous media, trapping air. It is also possible that the water level was increased too rapidly when the reservoir was raised prior to the drainage experiment, trapping air in the column. Regardless, the fact that portions of the column were below full saturation (so would exhibit low relative permeabilities) may in part explain the fact that the top of the column was at a relatively high saturation after an elapsed time of 20 h. Note that at the end of 20 h, the water content at the top sensor (S~0.55) was actually higher than the water content at the three sensors below it (S~0.4 for the second from the top).
From Figure 4b, it is apparent that the seemingly stabilized saturation values at the end of the 20 h period in Figure 4a were, in fact, still changing, but over a much longer time scale, as redistribution occurs. Saturations corresponding to the bottom sensors continue to rise over approximately 200 h (~8.3 days), although in the time scale of this figure, the portion of the column near the water table remains below full saturation. Similarly, the saturations corresponding to the top two sensors gradually decrease over approximately 200 h. Of particular interest, however, is the change in slopes of the saturation vs. time curves for the top two sensors, notable around 200 h for the top sensor (orange arrow) and just after 400 h for the second from the top (blue arrow). These slope changes appear to correspond to an evaporative drying front moving down from the top of the column. The top sensor was approximately 2.5 cm from the top of the packed sand, while the next one was approximately 7.5 cm from the top of the packed sand. Once the evaporative front reached the top sensor, the saturation at that sensor dropped rapidly, reaching zero (i.e., completely dry) around 300 h (~12 days from the start of the experiment). The saturation at the second sensor from the top of the column dropped to approximately 0.1 around 800 h (~33 days). Saturation at the top sensor remained essentially zero for the duration of the 1207-day experiment, while saturation at the second sensor did not continue to drop, although it fluctuated over time. The evaporative drying front did not appear to reach the third sensor.
Figure 5 shows elevation vs. saturation for sixteen different times within the first 436 days following the start of drainage. Also shown in the figure is a measured PcS relationship for the same US-Silica F45 sand, determined in earlier work [21,22] using methods based on a small (2.54 cm inner diameter × 1.27 cm height) soil cell described previously [10]. Notable in the figure is the fact that the porous medium near the water table began below full saturation, but gradually approached full saturation over hundreds of days as water redistributed in the column. This process was likely accelerated by the opening of column vents at day 128, a step that allowed trapped air to leave the column. The similarity of the late time elevation-saturation curves to the previously measured PcS relationships is also apparent from the figure, except for the two sensors near the top of the column. Previously measured residual saturation in the F45 sand was found to be approximately 0.2, yet the saturations measured by the top two sensors fell below that within the first month of the experiment and remained below it for the full duration of the experiment, an indication of the role of evaporation at the top of the column.
Figure 6 shows saturation measured by the 12 sensors over the full 1207-day duration of the experiment (Figure 6a), as well as details on the environmental conditions during the experiment, including outdoor temperature (Figure 6b), outdoor relative and absolute humidity (Figure 6c) and indoor temperature (Figure 6d). In addition, cumulative water additions over the course of the experiment are shown in Figure 6e. Note that the sensors monitoring saturations experienced several interruptions, likely due to brief power outages that were not immediately detected; most occurred later in the experiment, when the experiment was only monitored intermittently.
Looking at Figure 6a, one of the most interesting features of the experiment is the apparent evidence of seasonal fluctuations in the column saturations. While the water table (the level in the external reservoir) was kept within 0.25 cm of its initial position for the duration of the experiment, the saturation values at many levels in the column were observed to fluctuate by as much as 0.1, and the fluctuations appear to have a seasonal pattern. Except for the measured saturation at the top sensor, which remained nearly dry for most of the experiment, the greatest seasonal effects were observed higher in the column. In general, saturations tended to increase in the spring and summer, but decrease in the fall and winter. This is likely the result of changing indoor humidity.
Note the strong similarity between the seasonal trends in outdoor temperature (Figure 6b) and outdoor absolute humidity (Figure 6c) and column saturations. While outdoor relative humidity varies erratically over the course of a year (Figure 6c, dotted lines), with no strong seasonal trends, the absolute humidity (Figure 6c, solid lines) is strongly tied to temperature, because the saturated vapor pressure increases with increasing temperature. The total moisture content in the air in the summer exhibited peak summer values of ~250 g/m3, more than eight times the typical winter value of ~30 g/m3. That means that substantially more moisture is likely to be in indoor air exchanged with the outdoors during times when outdoor absolute humidity is high. Since indoor temperatures are relatively constant, that means that indoor relative humidity is likely much higher in the spring and summer than the fall and winter. It is relative humidity that controls evaporation rate, with higher relative humidities corresponding to less driving force for evaporation, and thus, slower evaporation. As such, slower evaporation would be expected indoors in the spring and summer, changing the balance between rate of evaporation and rate of capillary driven flow, allowing saturation to increase. Furthermore, note that the saturation change was reversible, with saturation increasing at times when evaporation would be expected to be slow, and decreasing at times when evaporation would be expected to be fast. This result highlights how the saturation profile results from a dynamic equilibrium between capillary-driven flow and evaporation.
Figure 6d shows the indoor temperature monitored throughout the experiment. Overall, temperature only varied within a narrow window over the 1207 days, typically staying within the 22–24 °C range for most of the 3.31 years, with a few exceptions, including one period in summer 2020, when the indoor temperature rose to ~26 °C for an extended period when the laboratories were closed for the pandemic, possibly due to a HVAC system problem. It is interesting to note that the indoor temperature was on the lower end of the range for the final year of the experiment, a time period that corresponded to some of the highest saturations at some levels, particularly at the level of the third sensor from the top (50.75 cm above the water table). This sustained lower indoor temperature may be responsible for the higher saturation late in the experiment, because a lower indoor temperature with the same absolute humidity would correspond to a higher indoor relative humidity, and thus, less driving force for evaporation.
Figure 6e shows the cumulative water addition in the column over the course of the experiment. Water was added intermittently as needed to maintain the level in the covered reservoir, particularly after the initial phase of the experiment. The large diameter of the reservoir relative to the column meant that even large volumes of water added to the reservoir produced extremely small changes in reservoir level. As such, the intermittent addition amounts are not a good measure of the water flux in the column over a short time scale. However, the cumulative amounts are of interest. A total of 533.5 mL of water was added to the reservoir over the duration of the experiment, more than eleven times the volume of water present in the column after the column water profile had stabilized early in the experiment. This means that there was significant upward flow over the course of experiment, with the water in the column continuously evaporating and being replaced by capillary-driven flow. Despite this, once the top of the column had dried, it remained dry for the duration of the experiment.
No water was added until 11 days into the experiment, due to rapid initial drainage. Subsequent water additions were initially high, as evidenced by the steep slope of the curve in Figure 6e, but then slowed, likely as a result of the growing dry layer of sand at the top of the column effectively slowing exchange with the atmosphere. Required water additions were high initially (213 mL in the first six months), but slowed to an average of 127 mL/year later in the experiment. It should be noted that, although the column was vented, the areas of the vents (0.3 cm2 each × 4) were small relative to the open top of the column (4.9 cm2), a fact that may explain why no apparent drying was observed near the vents either through the sensor response or visually. This may also explain why the drying front did not proceed down into the capillary fringe as was observed in the highly-vented sand experiment in earlier work [18]. Finally, although the cumulative upward flow was substantial over the duration of the experiment, the porewater velocity later in the experiment was negligibly small (127 mL/year corresponds to a porewater velocity of ~2 × 10−8 m/s in this system); as such, pressure gradients from flow are essentially negligible after the initial portion of the experiment (on the order of 10−5 cm water/cm sand), meaning that elevation in the column was a good measure of capillary pressure over much of the experiment. That means that the seasonal effects observed in Figure 6 result from the impact of environmental conditions on the PcS relationship itself.

4. Conclusions

Results of this work support the assertion that there is a link between the observed Pc–S relationship in a porous medium and the rate of evaporation. Strong seasonal effects were observed, with saturations increasing during times when the driving forces for evaporation would be expected to be low, and decreasing at times when the driving forces for evaporation would be expected to be high. Column saturation change occurred over different time scales, with initial imbibition and drainage occurring relatively rapidly within hours to tens of hours, but then evaporative drying effects becoming apparent at much longer time scales in the hundreds of hours. Seasonal changes occurred at still longer time scales, and corresponded to reversible increases and decreases in saturation based on environmental conditions, despite the stable water table level. While the existence of upward flow due to evaporation has long been recognized and plays an important role in the movement of contaminants and plant nutrients in the subsurface, the temporal saturation increases and decreases observed in this work are interesting. This is particularly true given that the PcS relationship is hysteretic, with imbibition at a given Pc corresponding to a lower S than drainage at that same Pc. That fact might create the expectation that something that causes a decrease in saturation at a given elevation above the water table might not be followed by a reversible increase in saturation, simply because an increase in saturation due to net flow of water is typically associated with imbibition, which predicts a lower saturation at a given elevation above the water table. In this work, however, it is apparent that the saturations following drainage decreased and increased reversibly around the values associated with drainage.
Many contaminants must pass through the unsaturated zone in order to reach groundwater, and the presence of both air and water in the unsaturated zone means that adsorption to air–water interfaces can significantly retard transport for some contaminants. This means that anything that impacts water content in the unsaturated zone has the potential to influence fate and transport. As such, seasonal effects on the PcS relationship might be expected to contribute to the vertical smearing of environmental contaminants in the unsaturated zone, and may play an important role in determining ultimate fate. Furthermore, recognition of the distinct evaporation-associated environmental effects on the Pc–S relationship may be needed for accurate long-term predictive simulations of transport of contaminants in the unsaturated zone, because use of the hysteretic Pc–S relationship in computations may yield incorrect predictions of water content when evaporation is the major driver of saturation change.
Evaporation at the top of the column was observed, with water in the top of the column evaporating to dryness. Despite substantial water flux upward in the column, the top of the column remained dry for the duration of the experiment. However, unlike previous work in a highly vented column [18], the drying front in this experiment did not continue to progress downward in the column, but instead was likely limited by a dynamic equilibrium between capillary-driven flow and evaporation.

Author Contributions

Conceptualization, A.O.A. and T.C.G.K.; methodology, A.O.A. and T.C.G.K.; investigation, A.O.A.; writing, A.O.A. and T.C.G.K.; software, T.C.G.K.; funding acquisition, T.C.G.K.; project administration, T.C.G.K. All authors have read and agreed to the published version of the manuscript.

Funding

This material is based upon work supported by the National Science Foundation under grant no. 1446264.

Data Availability Statement

Data are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elmaloglou, S.; Diamantopoulos, E. Effects of hysteresis on redistribution of soil moisture and deep percolation at continuous and pulse drip irrigation. Agric. Water Manag. 2009, 96, 533–538. [Google Scholar] [CrossRef]
  2. Levin, I.; Van Rooyen, P.C.; Van Rooyen, F.C. The Effect of Discharge Rate and Intermittent Water Application by Point-source Irrigation on the Soil Moisture Distribution Pattern. Soil Sci. Soc. Am. J. 1979, 43, 8–16. [Google Scholar] [CrossRef]
  3. Chen, L.X.; Kibbey, T.C.G. Transport of nanomaterials in the unsaturated zone. In Nanoscience and Nanotechnology: Environ-Mental and Health Impacts; Grassian, V.H., Ed.; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 107–132. [Google Scholar]
  4. Hoggan, J.L.; Sabatini, D.A.; Kibbey, T.C. Transport and retention of TiO2 and polystyrene nanoparticles during drainage from tall heterogeneous layered columns. J. Contam. Hydrol. 2016, 194, 30–35. [Google Scholar] [CrossRef]
  5. Yan, S.; Kibbey, T.C. The effect of evaporation-induced flow at the pore scale on nanoparticle transport and deposition in drying unsaturated porous media. J. Contam. Hydrol. 2019, 226, 103524. [Google Scholar] [CrossRef]
  6. Seetha, N.; Kumar, M.M.; Hassanizadeh, S.M. Modeling the co-transport of viruses and colloids in unsaturated porous media. J. Contam. Hydrol. 2015, 181, 82–101. [Google Scholar] [CrossRef] [PubMed]
  7. Marshall, J.D.; Shimada, B.W.; Jaffe, P.R. Effect of temporal variability in infiltration on contaminant transport in the unsaturated zone. J. Contam. Hydrol. 2000, 46, 151–161. [Google Scholar] [CrossRef]
  8. French, H.K.; Hardbattle, C.; Binley, A.; Winship, P.; Jakobsen, L. Monitoring snowmelt induced unsaturated flow and transport using electrical resistivity tomography. J. Hydrol. 2002, 267, 273–284. [Google Scholar] [CrossRef]
  9. Le, S.-T.; Gao, Y.; Kibbey, T.C.; Glamore, W.C.; O’Carroll, D.M. Predicting the impact of salt mixtures on the air-water interfacial behavior of PFAS. Sci. Total Environ. 2022, 819, 151987. [Google Scholar] [CrossRef]
  10. Suits, L.D.; Sheahan, T.C.; Chen, L.; Miller, G.A.; Kibbey, T.C.G. Rapid Pseudo-Static Measurement of Hysteretic Capillary Pressure-Saturation Relationships in Unconsolidated Porous Media. Geotech. Test. J. 2007, 30, 474–483. [Google Scholar] [CrossRef]
  11. Arya, L.M.; Paris, J.F. A Physicoempirical Model to Predict the Soil Moisture Characteristic from Particle-Size Distribution and Bulk Density Data. Soil Sci. Soc. Am. J. 1981, 45, 1023–1030. [Google Scholar] [CrossRef]
  12. Arya, L.M.; Leij, F.J.; van Genuchten, M.; Shouse, P.J. Scaling Parameter to Predict the Soil Water Characteristic from Particle-Size Distribution Data. Soil Sci. Soc. Am. J. 1999, 63, 510–519. [Google Scholar] [CrossRef] [Green Version]
  13. Daamen, C.C.; Simmonds, L.P. Measurement of Evaporation from Bare Soil and its Estimation Using Surface Resistance. Water Resour. Res. 1996, 32, 1393–1402. [Google Scholar] [CrossRef]
  14. Brutsaert, W. Daily evaporation from drying soil: Universal parameterization with similarity. Water Resour. Res. 2014, 50, 3206–3215. [Google Scholar] [CrossRef]
  15. Mosthaf, K.; Helmig, R.; Or, D. Modeling and analysis of evaporation processes from porous media on the REV scale. Water Resour. Res. 2014, 50, 1059–1079. [Google Scholar] [CrossRef]
  16. Haghighi, E.; Or, D. Linking evaporative fluxes from bare soil across surface viscous sublayer with the Monin–Obukhov atmospheric flux-profile estimates. J. Hydrol. 2015, 525, 684–693. [Google Scholar] [CrossRef]
  17. Tugwell-Wootton, T.; Skrzypek, G.; Dogramaci, S.; McCallum, J.; Grierson, P.F. Soil moisture evaporative losses in response to wet-dry cycles in a semiarid climate. J. Hydrol. 2020, 590, 125533. [Google Scholar] [CrossRef]
  18. Adegbule, A.O.; Kibbey, T.C. Exploring the impact of long-term evaporation on the relationship between capillary pressure and water saturation in unsaturated porous media. J. Hydrol. 2020, 582, 124557. [Google Scholar] [CrossRef]
  19. Brock, F.V.; Crawford, K.C.; Elliott, R.L.; Cuperus, G.W.; Stadler, S.J.; Johnson, H.L.; Eilts, M.D. The Oklahoma Mesonet: A Technical Overview. J. Atmos. Ocean. Technol. 1995, 12, 5–19. [Google Scholar] [CrossRef]
  20. McPherson, R.A.; Fiebrich, C.A.; Crawford, K.C.; Kilby, J.R.; Grimsley, D.L.; Martinez, J.E.; Basara, J.; Illston, B.G.; Morris, D.A.; Kloesel, K.A.; et al. Statewide Monitoring of the Mesoscale Environment: A Technical Update on the Oklahoma Mesonet. J. Atmos. Ocean. Technol. 2007, 24, 301–321. [Google Scholar] [CrossRef] [Green Version]
  21. Kibbey, T.C.G.; Chen, L. A pore network model study of the fluid-fluid interfacial areas measured by dynamic-interface tracer depletion and miscible displacement water phase advective tracer methods. Water Resour. Res. 2012, 48, W10519. [Google Scholar] [CrossRef]
  22. Kibbey, T.C.G.; Naghavi, B.; Hou, L. Understanding the Role of Water Phase Continuity on the Relationship Between Drainage Rate and Apparent Residual Saturation: A Dynamic Network Model Study. Transp. Porous Media 2016, 112, 689–705. [Google Scholar] [CrossRef]
Figure 1. Column configuration used in this study. (a) The column; (b) Side view illustrating sensor and vent positioning; (c) Top view illustrating sensor positioning; (d) A custom interdigitated capacitive sensor of the kind used for this work.
Figure 1. Column configuration used in this study. (a) The column; (b) Side view illustrating sensor and vent positioning; (c) Top view illustrating sensor positioning; (d) A custom interdigitated capacitive sensor of the kind used for this work.
Water 14 03294 g001
Figure 2. Elevation above the water table (z) vs saturation (S) measured over the first 12 h of the imbibition experiment.
Figure 2. Elevation above the water table (z) vs saturation (S) measured over the first 12 h of the imbibition experiment.
Water 14 03294 g002
Figure 3. Saturation at the twelve sensors measured over the first 20 h of the imbibition experiment. Sensor elevation (z) is reported from the water table.
Figure 3. Saturation at the twelve sensors measured over the first 20 h of the imbibition experiment. Sensor elevation (z) is reported from the water table.
Water 14 03294 g003
Figure 4. Saturation at the twelve sensors over short (a) and intermediate (b) timeframes at the start of the drainage experiment. Sensor elevation (z) is reported from the water table. Arrows indicate when the evaporative drying front from the top of the column reached the top two sensors.
Figure 4. Saturation at the twelve sensors over short (a) and intermediate (b) timeframes at the start of the drainage experiment. Sensor elevation (z) is reported from the water table. Arrows indicate when the evaporative drying front from the top of the column reached the top two sensors.
Water 14 03294 g004
Figure 5. Elevation above the water table (z) vs. saturation (S) measured over 436 days from the start of the imbibition experiment. The dotted line corresponds to a previously measured PcS relationship for the F45 sand.
Figure 5. Elevation above the water table (z) vs. saturation (S) measured over 436 days from the start of the imbibition experiment. The dotted line corresponds to a previously measured PcS relationship for the F45 sand.
Water 14 03294 g005
Figure 6. Saturation at the twelve sensors over the full 1207 days of the drainage experiment (a), and information on conditions over the duration of the experiment: (b) Outdoor temperature; (c) Outdoor absolute and relative humidity (note: relative humidity scale is not shown, but varies from 0 to 100% over vertical axis); (d) Indoor temperature; (e) Cumulative water added. Periods where sensor data collection experienced interruptions are indicated with dotted lines. Sensor elevation (z) in (a) is reported from the water table.
Figure 6. Saturation at the twelve sensors over the full 1207 days of the drainage experiment (a), and information on conditions over the duration of the experiment: (b) Outdoor temperature; (c) Outdoor absolute and relative humidity (note: relative humidity scale is not shown, but varies from 0 to 100% over vertical axis); (d) Indoor temperature; (e) Cumulative water added. Periods where sensor data collection experienced interruptions are indicated with dotted lines. Sensor elevation (z) in (a) is reported from the water table.
Water 14 03294 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adegbule, A.O.; Kibbey, T.C.G. The Influence of Evaporation and Seasonal Effects on the Water Content in the Unsaturated Zone: A Multi-Year Laboratory Study. Water 2022, 14, 3294. https://doi.org/10.3390/w14203294

AMA Style

Adegbule AO, Kibbey TCG. The Influence of Evaporation and Seasonal Effects on the Water Content in the Unsaturated Zone: A Multi-Year Laboratory Study. Water. 2022; 14(20):3294. https://doi.org/10.3390/w14203294

Chicago/Turabian Style

Adegbule, Aderonke O., and Tohren C. G. Kibbey. 2022. "The Influence of Evaporation and Seasonal Effects on the Water Content in the Unsaturated Zone: A Multi-Year Laboratory Study" Water 14, no. 20: 3294. https://doi.org/10.3390/w14203294

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop