Thursday, June 23, 2016

What processes drive soil moisture dynamics?


The primary processes driving soil moisture dynamics in the Vacant to Vibrant parcels are rainfall events that act as sources (both direct rainfall and runoff) and losses due to infiltration, evaporation, and plant mediated evapotranspiration.  In the previous blog post (here), I presented evidence that the patterns of diurnal variation in soil moisture were driven by variation in soil temperature and did not reflect the effects of the primary drivers.  The regularity of the diurnal variability and its low amplitude allows use of smoothing functions to characterizing longer term declines in soil moisture associated with the loss processes (infiltration, evaporation, and evapotranspiration).  Figure 1 shows the application a simple linear regression to a time series between rainfall events.  Using the linear fit to the time series, Figure 2 shows the that the removal of the longer-term trend emphasizes the diurnal variability, and Figure 3 shows that the variation in the residuals are correlated with observed soil temperature.
Figure 1.  Semi logarithmic plot of decline of soil moisture over the period May 25 to June 6, 2016 for the Gary E1 parcel at 3 cm.  The blue line is a linear regression.





Figure 2.  Plot of the residuals for the regression in Figure 1.

Figure 3.  Plot of the residuals from Figure 2 showing association with measured soil temperature at 3 cm depth.  The blue line is a regression between the two variables (r= 0.80, p <0.0001).


The advantage of the linear regression in Figure 1 is that the slope is a first order (i.e. an exponential decay rate) estimate of the rate of change of soil moisture between rainfall events.  The estimated rate for the Gary E1 parcel at 3 cm soil depth is -1.68e-07 (1/s).  This method can be applied to trends at other depths or for the weighted average soil moisture and provides data for comparison of experimental and control parcels in the different soil types of the neighborhoods of the three cities in the project.

The remaining decline pattern (presented in the previous blog post) is the transient associated with a rainfall event.  Figure 4, presents the observed relation between increment of soil moisture and the subsequent first order rate of decline.



Figure 4.  Relation between soil moisture increment and subsequent rate of decline in soil moisture at 3 cm depth in the Gary E1 rain garden.  

Figure 5.  Relation between soil moisture increment and subsequent rate of decline in soil moisture at 3 cm depth in the Gary E1 rain garden for soil moisture increments greater than 0.012.  The blue line is a regression line and the shaded area represents the standard error bounds.

The slope of the relation in Figure 5 is -1.488e-05 (r=-0.5408624, p=0.1663, NS).  The hint of an inverse relation between the increment of soil moisture following a rainfall event exists, but is not statistically significant.  The average rate of decline, however, is 1.02e-6 (1/s).    An outlier occurs on April 30, 2016.  On 4/29 and 4/30 (see Figure 6), there was a double increase.  The identification of extrema near this interval is a problem and it is reasonable to regard the April 30 point as anomalous.  Figure 7 shows the result of eliminating this point.  The estimate of the decline rate is a statistically significant -1.90e-5, which is nearly two orders of magnitude greater than the rate of decline of the inter-event interval in Figure 1.  Clearly, more data are needed to explore this relationship and its drivers, but it seems reasonable to expect that the rate of decline following a rainfall event is a function of the soil moisture gradient and the permeability of the soil.

Figure 6.  Pattern of variation in soil moisture (m3/m3) at 3 cm soil depth in the Gary E1 rain garden.

Figure 7.  Relation between soil moisture increment and subsequent rate of decline in soil moisture at 3 cm depth in the Gary E1 rain garden for soil moisture increments greater than 0.012 and eliminating.  The blue line is a regression line (r= -0.96, p<0.001) and the shaded area represents the standard error bounds.



Tuesday, June 21, 2016

What Causes Short-term Variability in Soil Moisture?


On June 8, 2017, Ryan, Sandra, and I completed installation of a weather station and four soil moisture monitoring platforms in the PUSH neighborhood of Buffalo for the Vacant to Vibrant Project.  With all of the monitoring platforms now in place for the project, we need to have quantitative methods with which to analyze the data.  The goal of the analysis methods is to measure the storage and processing capacity of the green infrastructure installed in the experimental parcels of the three cities in the project.   As documented in a previous post (here), experimental and control plots display patterns of variability in soil moisture at various time scales. 

Figure 1 illustrates the types of variability patterns observed in soil moisture at four of the sites in Gary’s Aetna neighborhood.   Basically, there are three types of variability in Figure 1.  Rainfall events have an abrupt increase in soil moisture followed by a rapid decline (see selection A in Figure 1).  Between rainfall events there is a slower decline in soil moisture (see selection B in Figure 1) with a pattern of regular short-term variability (see selection C in Figure 1).  This blog post focuses on finding the drivers of short-term variability.


Figure 1.  Variation in soil moisture at 3 cm depth in four study sites in the Aetna neighborhood of Gary.  Parcels are as follows: 1200 Oklahoma (Gary E1), 1035 Oklahoma (Gary E2), 1252 Dakota (Gary E3), and 910 Idaho (Gary C1).  The emphasized patterns include:  A—changes associated with rainfall events, B—steady decline during rain free intervals, and C—short-term variability in the decline time series between rain events.
The regular pattern of short-term variability appears to be diurnal.  Figure 2 shows an expanded view of the decline during the May 28 to June 6 interval in Figure 1, selection C.  The diurnal patterns of soil moisture variability and temperature are synchronous.  An obvious hypothesis is that these short-term changes in soil moisture are driven by temperature changes.

Figure 2.  Comparison of patterns of change in soil moisture and temperature over the interval of May 26 to June 6 at 3 cm in the Gary E1 site.
If temperature changes are driving short-term variability of soil moisture, there should be both theoretical and experimental evidence to support the claim.  From a theoretical perspective, there is ample evidence indicating that soil moisture measurements are temperature sensitive.  The soil moisture probes that we use (Onset Part No. S-SMx-M005) use time domain reflectometry (TDR) method to estimate volumetric water content.  It is well known that temperature dependence of TDR includes both instrumental and physical effects (Or and Wraith, 1999), and there are methods of correcting for instrumental sensitivity to temperature change (e.g. Chanzy et al., 2012).   Because the physical effects of temperature on soil moisture reflects the underlying influence of temperature on the bound to free ratio of soil moisture, soil composition and soil moisture content interact with temperature effects to limit the generality of correction methods (Seyfried and Grant, 2007).  Nevertheless, diurnal variability of soil moisture with temperature changes is an expected pattern.

From an empirical perspective, diurnal variability of soil moisture should reflect diurnal variation in soil temperature.  Using the extrema identification method discussed in the previous post, I estimate the soil moisture increments for the entire time series for Gary E1 in Figure 1.  Figure 3 shows the resulting distribution
Figure 3.  Histogram of soil moisture (SM, m3/m3) for the Gary E1 time series in Figure 1.  The red vertical line is the mean increment, 0.012, and the blue vertical lines are the quartiles.
The distribution in Figure 3 is skewed with the mean near the third quartile.  Using the mean as a demarcation of diurnal soil moisture increments from rainfall event increments, Figure 4 shows a very strong linear dependence of soil moisture increment on temperature change.  Figure 5 provides evidence that the subsequent decline rate also depends on temperature change.
Figure 4.  Association of soil moisture increment and temperature change for soil moisture increments less than 0.012.  In blue, is the regression line (r=0.95, p<0.0001).   The shaded area represents the standard error bounds.



Figure 5.  Association between temperature change and the instantaneous slope of the decline in soil moisture for increments of soil moisture less than 0.012.  In blue, is the regression line (r=0.77, p<0.0001).   The shaded area represents the standard error bounds.

As noted by Or and Wraith (1999), separation temperature-dependent variation of soil moisture from changes associated with evapotranspiration and soil infiltration is necessary to understand the contributions of evapotranspiration and infiltration in individual experimental and control parcels of the V to V project.  Although there are proposed methods for removal of temperature-dependent transients, their lack general reliability across types of soils and range of soil moisture conditions is a limitation (Seyfried and Grant, 2007).  Furthermore, as shown in Figure 6, the diurnal pattern of variability varies with depth, reflecting both attenuation and phase shift of temperature variation.  Even with temperature correction, diurnal fluctuation of soil moisture associated with temperature dependent release of bound soil moisture will remain in the time series.  Given the regularity of the diurnal variation of soil moisture, a smoothing technique may offer a satisfactory method to estimate rates of decline during longer event free intervals.



Figure 6.  Comparison of soil moisture variability in Gary E1 for three soil depths (3 cm, bottom black line; 10 cm, middle red line; and 20 cm, top blue line).

References:

Chanzy, AndrĂ©, Gaudu, Jean-Claude, and Marloie, Olivier.  2012.  Correcting the Temperature Influence on Soil Capacitance Sensors Using Diurnal Temperature and Water Content Cycles.  Sensors, 12, 9773-9790.

Or, Dani, and Wraith, Jon M.  1999.  Temperature effects on soil bulk dielectric permittivity measured by time domain reflectometry: A physical model.  Water Res. Research, 35(2): 371-383.

Seyfried, Mark S.  and Grant, Laura E.  2007.  Temperature Effects on Soil Dielectric Properties Measured at 50 MHz.  Vadose Zone J. 6:759–765.