**Time Stability of Soil Water Content**

Wei Hu, Lindsay K. Tallon, Asim Biswas and Bing Cheng Si

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/52469

## **1. Introduction**

Soil water is a key variable controlling water and energy fluxes in soils (Vereecken et al., 2007). It is necessary for plant and vegetation growth and development. Research has indi‐ cated that soil water content (*SWC*) varies both in space and time. Variations in both space and time present a substantial challenge for applications such as precision agriculture and soil water management.

Since the contribution of Vachaud et al. (Vachaud et al., 1985), a large body of research has indicated the presence of time stability of *SWC* (Biswas and Si, 2011c; Comegna and Basile, 1994; Grayson and Western, 1998; Hu et al., 2009; Martínez-Fernández and Ceballos; 2005; Mohanty and Skaggs, 2001), which means that the spatial pattern of *SWC* does not change with time at a certain probability. According to this concept, if a field is repeatedly surveyed for *SWC*, there is a high probability that a location with certain wetness characteristics (i. e., wet, dry, intermediate) will maintain those characteristics on subsequent occasions. Time stability has also been extended to describe the characteristics of *SWC* at point scales. A loca‐ tion will be regarded as time stable provided it can estimate the average *SWC* of an area.

Time stability of *SWC* has been observed at a large variety of scales ranging from plot (Pa‐ chepsky et al., 2005) to region (Martínez-Fernández and Ceballos, 2003) and related studies cover a range of investigated areas, sampling schemes, sampling depths, investigation peri‐ ods, and land uses (Biswas and Si; 2011c; Brocca et al., 2009; Cosh et al., 2008; Hu et al., 2010a;Tallon and Si, 2004; Vachaud et al., 1985). As a result, a variety of methods have been developed to evaluate time stability of *SWC*, each with its own advantages and disadvantag‐ es. Time stability is usually used to characterize time persistence of the spatial pattern of *SWC* between measurement occasions, either at the measurement scale or at different scales (Biswas and Si, 2011c; Cosh et al., 2006; Kachanoski and De Jong, 1988; Vachaud et al., 1985).

© 2013 Hu et al.; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 Hu et al.; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

On the other hand, *SWC* at time stable locations can be used to estimate the spatial average *SWC* of a given area (Grayson and Western, 1998). Therefore, quick and accurate estimation of soil water content at a field or large areas may be possible with only one representative location. Areal estimation of *SWC* from point source data has the potential to substantially reduce both the capital and labor costs involved in estimating average *SWC*, making the method appealing to a wide range of disciplines.

pattern of *SWC* is scale dependent because of the interaction between measurement loca‐ tions (Kachanoski and De Jong, 1988). Scale dependence of the time stability was also inves‐ tigated by different methods including spatial coherency analysis (Kachanoski and De Jong, 1988), wavelet coherency analysis (Biswas and Si, 2011c), and multivariate empirical mode decomposition (Rehman and Mandic, 2010). A detailed introduction of these methods is giv‐

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 49

For time stability at point scales, *SWC* at each location is usually scaled in terms of the field average (Vachaud et al., 1985). This indicates that at a given location, the change in scaled *SWC* exhibits little change over time. Many indices are available to examine the degree of time stabil‐ ity at point scales. These indices include standard deviation of relative difference (Mohanty and Skaggs, 2001; Schneider et al., 2008; Vachaud et al., 1985), root mean square error (Bosch et al., 2006; Jacobs et al., 2004), width of the 90% empirical tolerance interval of relative water con‐ tent (Guber et al., 2008), chi-squared statistic (Guber et al., 2008), root-mean-squared differen‐ ces (Guber et al., 2008), temporal coefficient of variability (Starr, 2005), and mean absolute bias error (Hu et al., 2010a, 2010b). One of the applications of time stability at point scales is to iden‐ tify catchment average soil moisture monitoring (CASMM) site (Grayson and Western, 1998) or benchmark location (Tallon and Si, 2004) for average *SWC* evaluation. Detailed introduc‐

Generally, these two types of time stability are correlated. If the spatial pattern is time stable, there is a larger possibility of time stability at point scales. On the other hand, if more points are time stable, there is a larger possibility of time stability of spatial pattern. However, these two types of time stability are also distinguished for the following three aspects: (1) time stability of spatial pattern is used to evaluate the time stability of an area, while time stability of relative *SWC* is used to evaluate the time stability of a point. Therefore, time in‐ stability of spatial pattern does not mean the absence of time stable points (Grayson and Western, 1998; Schneider, 2008). On the other hand, if no points are time stable, it does not mean that spatial pattern is time unstable for any two measurement occasions; (2) they have different evaluation criteria as mentioned above, and (3) they have different applications. Time stability of spatial pattern is used to qualitatively describe the similarity of spatial dis‐ tribution of *SWC* between different measurement occasions for the better understanding of soil water related process and influencing factors (Lin; 2006), while time stable locations can be selected to estimate the average *SWC* for an area (upscaling) by directly using the soil wa‐ ter content of the time stable location (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Martínez-Fernández and Ceballos, 2005) or by considering the offset between the mean value and that of the time stable location (Grayson and Western, 1998; Hu et al.,

There are different methods for analyzing time stability. We have documented the methods

in two sections: (1) time stability of spatial patterns and (2) time stability of points.

tion of these indices and applications is given in Section 3.

en in Section 3.

2010b; Starks et al., 2008).

**3. Methodology**

Time stability of *SWC* is controlled by various factors including soil, topography, vegetation, and climate (Brocca et al., 2009; Gómez-Plaza et al., 2001; Grayson and Western, 1998; Hu et al, 2010a; Tallon and Si, 2004). Information on the controls of time stability provides an es‐ sential insight into the mechanisms of soil water movement and storage. In addition, it is also important to identify the time stable location for average *SWC* estimation *a priori*.

The goal of this chapter is to provide a comprehensive review on the study of time stability of *SWC*. In doing so, we have introduced the concept, methodology, application, and the controlling factors step by step.

## **2. Concept**

Soil water content varies with time and space. However, if a field is repeatedly surveyed for *SWC*, there is a high probability that a location with a certain relative wetness condition (i. e., wet, dry, and intermediate) at a given time will remain at the same relative wetness con‐ dition at other times. This phenomenon was first explained as time stability by Vachaud et al. (1985). This can be defined as the "time invariant association between spatial location and classical statistical parametric values" of *SWC*, most often the mean.

Time stability of *SWC* can be divided into two types, one is to describe the overall similarity of the spatial pattern between measurements, and the other is to describe the time-invari‐ ance of the relative *SWC* of a given location.

Repeated survey indicated that some locations are always wetter or drier than the average *SWC* of an area, resulting in the preservation of their ranks. Hence, the spatial pattern of *SWC* measured at one time will be similar to those measured on subsequent occasions. Stat‐ istically, this kind of time stability indicates that the orders, or ranks, of *SWC* at different lo‐ cations do not change over time at some probability. If all locations were to maintain their rank on subsequent measurement occasions, then the spatial pattern of *SWC* will be identi‐ cal. From this point, Chen (2006) augured that the term of "rank stability" or "order stabili‐ ty'' should be better than "time stability" because "the stability is the order or rank of a soil property at different spatial points that does not change at some probability".

The most widely used method to describe time stability of *SWC* spatial patterns is Spear‐ man's rank correlation analysis (Vachaud et al., 1985), while Pearson correlation analysis has also been widely used (Cosh et al., 2004). This is because Pearson correlation measures strictly linear relationships, while Spearman Rank correlation also measures nonlinear corre‐ lation between SWC values measured at two different occasions. Time stability of the spatial pattern of *SWC* is scale dependent because of the interaction between measurement loca‐ tions (Kachanoski and De Jong, 1988). Scale dependence of the time stability was also inves‐ tigated by different methods including spatial coherency analysis (Kachanoski and De Jong, 1988), wavelet coherency analysis (Biswas and Si, 2011c), and multivariate empirical mode decomposition (Rehman and Mandic, 2010). A detailed introduction of these methods is giv‐ en in Section 3.

For time stability at point scales, *SWC* at each location is usually scaled in terms of the field average (Vachaud et al., 1985). This indicates that at a given location, the change in scaled *SWC* exhibits little change over time. Many indices are available to examine the degree of time stabil‐ ity at point scales. These indices include standard deviation of relative difference (Mohanty and Skaggs, 2001; Schneider et al., 2008; Vachaud et al., 1985), root mean square error (Bosch et al., 2006; Jacobs et al., 2004), width of the 90% empirical tolerance interval of relative water con‐ tent (Guber et al., 2008), chi-squared statistic (Guber et al., 2008), root-mean-squared differen‐ ces (Guber et al., 2008), temporal coefficient of variability (Starr, 2005), and mean absolute bias error (Hu et al., 2010a, 2010b). One of the applications of time stability at point scales is to iden‐ tify catchment average soil moisture monitoring (CASMM) site (Grayson and Western, 1998) or benchmark location (Tallon and Si, 2004) for average *SWC* evaluation. Detailed introduc‐ tion of these indices and applications is given in Section 3.

Generally, these two types of time stability are correlated. If the spatial pattern is time stable, there is a larger possibility of time stability at point scales. On the other hand, if more points are time stable, there is a larger possibility of time stability of spatial pattern. However, these two types of time stability are also distinguished for the following three aspects: (1) time stability of spatial pattern is used to evaluate the time stability of an area, while time stability of relative *SWC* is used to evaluate the time stability of a point. Therefore, time in‐ stability of spatial pattern does not mean the absence of time stable points (Grayson and Western, 1998; Schneider, 2008). On the other hand, if no points are time stable, it does not mean that spatial pattern is time unstable for any two measurement occasions; (2) they have different evaluation criteria as mentioned above, and (3) they have different applications. Time stability of spatial pattern is used to qualitatively describe the similarity of spatial dis‐ tribution of *SWC* between different measurement occasions for the better understanding of soil water related process and influencing factors (Lin; 2006), while time stable locations can be selected to estimate the average *SWC* for an area (upscaling) by directly using the soil wa‐ ter content of the time stable location (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Martínez-Fernández and Ceballos, 2005) or by considering the offset between the mean value and that of the time stable location (Grayson and Western, 1998; Hu et al., 2010b; Starks et al., 2008).

## **3. Methodology**

On the other hand, *SWC* at time stable locations can be used to estimate the spatial average *SWC* of a given area (Grayson and Western, 1998). Therefore, quick and accurate estimation of soil water content at a field or large areas may be possible with only one representative location. Areal estimation of *SWC* from point source data has the potential to substantially reduce both the capital and labor costs involved in estimating average *SWC*, making the

Time stability of *SWC* is controlled by various factors including soil, topography, vegetation, and climate (Brocca et al., 2009; Gómez-Plaza et al., 2001; Grayson and Western, 1998; Hu et al, 2010a; Tallon and Si, 2004). Information on the controls of time stability provides an es‐ sential insight into the mechanisms of soil water movement and storage. In addition, it is

The goal of this chapter is to provide a comprehensive review on the study of time stability of *SWC*. In doing so, we have introduced the concept, methodology, application, and the

Soil water content varies with time and space. However, if a field is repeatedly surveyed for *SWC*, there is a high probability that a location with a certain relative wetness condition (i. e., wet, dry, and intermediate) at a given time will remain at the same relative wetness con‐ dition at other times. This phenomenon was first explained as time stability by Vachaud et al. (1985). This can be defined as the "time invariant association between spatial location and

Time stability of *SWC* can be divided into two types, one is to describe the overall similarity of the spatial pattern between measurements, and the other is to describe the time-invari‐

Repeated survey indicated that some locations are always wetter or drier than the average *SWC* of an area, resulting in the preservation of their ranks. Hence, the spatial pattern of *SWC* measured at one time will be similar to those measured on subsequent occasions. Stat‐ istically, this kind of time stability indicates that the orders, or ranks, of *SWC* at different lo‐ cations do not change over time at some probability. If all locations were to maintain their rank on subsequent measurement occasions, then the spatial pattern of *SWC* will be identi‐ cal. From this point, Chen (2006) augured that the term of "rank stability" or "order stabili‐ ty'' should be better than "time stability" because "the stability is the order or rank of a soil

The most widely used method to describe time stability of *SWC* spatial patterns is Spear‐ man's rank correlation analysis (Vachaud et al., 1985), while Pearson correlation analysis has also been widely used (Cosh et al., 2004). This is because Pearson correlation measures strictly linear relationships, while Spearman Rank correlation also measures nonlinear corre‐ lation between SWC values measured at two different occasions. Time stability of the spatial

property at different spatial points that does not change at some probability".

classical statistical parametric values" of *SWC*, most often the mean.

ance of the relative *SWC* of a given location.

also important to identify the time stable location for average *SWC* estimation *a priori*.

method appealing to a wide range of disciplines.

controlling factors step by step.

48 Advances in Agrophysical Research

**2. Concept**

There are different methods for analyzing time stability. We have documented the methods in two sections: (1) time stability of spatial patterns and (2) time stability of points.

#### **3.1. Time stability of spatial patterns**

#### *3.1.1. Pearson correlation analysis*

Pearson correlation coefficient is sensitive only to a linear relationship between two varia‐ bles. A Pearson correlation coefficient, *rj,j'*, between two spatial series of soil water content measured at time *j* and *j'* can be defined by

$$r\_{j,j} = \frac{\sum\_{i} (\text{SWC}\_{j}(i) - \left\{ \text{SWC}\_{j} \right\}) (\text{SWC}\_{j}(i) - \left\{ \text{SWC}\_{j} \right\})}{\sqrt{\sum\_{i} (\text{SWC}\_{j}(i) - \left\{ \text{SWC}\_{j} \right\})^{2}} \sqrt{\sum\_{i} (\text{SWC}\_{j}(i) - \left\{ \text{SWC}\_{j} \right\})^{2}}} \tag{1}$$

*3.1.3. Average spatial coefficient of determination and temporal coefficient of variability*

be a time independent scaling factor, *r* (*i*), which is expressed as (Starr, 2005)

*SWC i*

needs an extra term to fit a more general situation (Starr, 2005)

*SWC i*

*SWC*

the variance of *sj* (*i*) over space at any time, *γ<sup>i</sup>*

<sup>2</sup> *<sup>r</sup>* (*i*) and *γ<sup>i</sup>*

where *γ<sup>i</sup>*

where *γ<sup>i</sup>*

The *γ<sup>i</sup>*

g

Then, the average spatial coefficient of determination, *Rs*

<sup>2</sup> *sj* (*i*) is the mean value of *γ<sup>i</sup>*

*j*

If perfectly time stable soil water content pattern exists, the ratio of *SWC <sup>j</sup>* (*i*) to *SWC <sup>j</sup>* , can

( ) ( ) *<sup>j</sup>*

Complete time stability is usually disturbed by a series of processes and therefore Eq. (3)

( ) () () () *<sup>j</sup>*

where *sj* (*i*) is a scaled factor depending on both time and location, and *ej* (*i*) is an addi‐ tional term that accounts for random measurement errors, random sampling variability, and any true deviations from time stability. Assuming that *ej* (*i*) and *r* (*i*) are independent, then

2 22 () () () *ij i ij*

the total variance explained by the time stability model of Eq. (3) can be written as

2

*s*

*R*

2

( ) ( )

*r i*

*s i*

2

*i j*

value of zero over time, *r*(*i*) can be approximated by time average of *sj* (*i*), *sj* (*i*) , i. e.,

é ù ë û <sup>=</sup> é ù ë û

*i*

g

g

gg

<sup>2</sup> *sj* (*i*) can be calculated by averaging the calculated variances for each observation.

*j j*

*ri e i s i*

*r i*

*SWC* <sup>=</sup> (3)

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 51

=+ = (4)

<sup>2</sup> *sj* (*i*) , can be written as

<sup>2</sup> *ej* (*i*) are the variances of *<sup>r</sup>* (*i*) and *ej* (*i*), respectively over space.

2

é ù *s i ri e i* = + é ù é ù ë û ë û ë û (5)

, which describes the proportion of

<sup>2</sup> *sj* (*i*) over time. Given *ej* (*i*) with an average

() () *<sup>j</sup> ri s i* = (7)

(6)

*j*

where *SWC <sup>j</sup>* (*i*) and *SWC <sup>j</sup>* are soil water content at location *i* and spatial average at a giv‐ en time *j*, respectively. *SWC <sup>j</sup> '* (*i*) and *SWC j'* are soil water content at location *i* and spatial average at another time *j '* , respectively. The resulting coefficients refer to the correlation of the spatial patterns of *SWC* from one time to another. It is expected that closely correlated patterns have a *rj*, *j'* near one, while uncorrelated patterns are indicated by *rj,j'* values near zero (Cosh et al., 2004). Student's *t* is generally used to test the significance of Pearson corre‐ lation coefficient. It can be implemented by some statistical software, such as Excel, SPSS, MathCad, Matlab, and SAS.

#### *3.1.2. Spearman's rank correlation analysis*

Spearman's rank correlation coefficient is a non-parametric measure of statistical depend‐ ence between two variables. It assesses how well the relationship between two variables can be described using a monotonic function. It is the Pearson correlation between the ranks of one series and the ranks of another series. Because ranking linearize some of the nonlinear relationships, it is sensitive to nonlinear relationships. Let *Rij* be the rank of soil water con‐ tent *SWC <sup>j</sup>* (*i*) and *Ri ' j* the rank of *SWC j'* (*i*). The Spearman rank correlation coefficient, *rs*, is calculated by

$$r\_s = 1 - \frac{6\sum\_{i=1}^{n} (R\_{ij} - R\_{ij})}{n(n^2 - 1)}\tag{2}$$

where *n* is the number of observations. A value of *rs*=1 corresponds to identity of rank for any sites, or perfect time stability between time *j* and *j'*. The closer *rs* is to 1, the more stable the spatial pattern will be. Student's *t* can be used to test the significance of *rs* (Zar, 1972). Software such as Excel, SPSS, MathCad, Matlab, and SAS can implement this analysis.

#### *3.1.3. Average spatial coefficient of determination and temporal coefficient of variability*

**3.1. Time stability of spatial patterns**

measured at time *j* and *j'* can be defined by

å

Pearson correlation coefficient is sensitive only to a linear relationship between two varia‐ bles. A Pearson correlation coefficient, *rj,j'*, between two spatial series of soil water content

( ) ( )

*SWC i SWC SWC i SWC*

*j j <sup>j</sup> <sup>j</sup> <sup>i</sup>*

( )( )


( ) ( )

the spatial patterns of *SWC* from one time to another. It is expected that closely correlated

zero (Cosh et al., 2004). Student's *t* is generally used to test the significance of Pearson corre‐ lation coefficient. It can be implemented by some statistical software, such as Excel, SPSS,

Spearman's rank correlation coefficient is a non-parametric measure of statistical depend‐ ence between two variables. It assesses how well the relationship between two variables can be described using a monotonic function. It is the Pearson correlation between the ranks of one series and the ranks of another series. Because ranking linearize some of the nonlinear relationships, it is sensitive to nonlinear relationships. Let *Rij* be the rank of soil water con‐

( ')


*R R*

( 1)

where *n* is the number of observations. A value of *rs*=1 corresponds to identity of rank for any sites, or perfect time stability between time *j* and *j'*. The closer *rs* is to 1, the more stable the spatial pattern will be. Student's *t* can be used to test the significance of *rs* (Zar, 1972). Software such as Excel, SPSS, MathCad, Matlab, and SAS can implement this analysis.

*ij ij <sup>i</sup>*

1 2

= - -

*n n* =

*n*

6

1

*s*

*r*

2


'

(*i*) and *SWC <sup>j</sup>* are soil water content at location *i* and spatial average at a giv‐

'

å å (1)

'

(*i*) and *SWC j'* are soil water content at location *i* and spatial

, respectively. The resulting coefficients refer to the correlation of

near one, while uncorrelated patterns are indicated by *rj,j'* values near

the rank of *SWC j'* (*i*). The Spearman rank correlation coefficient, *rs*, is

å (2)

'

( ) ( )

*SWC i SWC SWC i SWC*

*j j <sup>j</sup> <sup>j</sup> i i*

, 2 2

*3.1.1. Pearson correlation analysis*

50 Advances in Agrophysical Research

'

=

*j j*

en time *j*, respectively. *SWC <sup>j</sup> '*

average at another time *j '*

MathCad, Matlab, and SAS.

*3.1.2. Spearman's rank correlation analysis*

*' j*

patterns have a *rj*, *j'*

tent *SWC <sup>j</sup>* (*i*) and *Ri*

calculated by

*r*

where *SWC <sup>j</sup>*

If perfectly time stable soil water content pattern exists, the ratio of *SWC <sup>j</sup>* (*i*) to *SWC <sup>j</sup>* , can be a time independent scaling factor, *r* (*i*), which is expressed as (Starr, 2005)

$$\frac{\text{SWC}\_{j}(i)}{\left\langle \text{SWC}\_{j} \right\rangle} = r(i) \tag{3}$$

Complete time stability is usually disturbed by a series of processes and therefore Eq. (3) needs an extra term to fit a more general situation (Starr, 2005)

$$\frac{\text{SWC}\_{j}(i)}{\left\{\text{SWC}\_{j}\right\}} = r(i) + e\_{j}(i) = s\_{j}(i) \tag{4}$$

where *sj* (*i*) is a scaled factor depending on both time and location, and *ej* (*i*) is an addi‐ tional term that accounts for random measurement errors, random sampling variability, and any true deviations from time stability. Assuming that *ej* (*i*) and *r* (*i*) are independent, then the variance of *sj* (*i*) over space at any time, *γ<sup>i</sup>* <sup>2</sup> *sj* (*i*) , can be written as

$$\left[\gamma\_i^2\left[s\_j(i)\right]\right] = \gamma\_i^2\left[r(i)\right] + \gamma\_i^2\left[e\_j(i)\right] \tag{5}$$

where *γ<sup>i</sup>* <sup>2</sup> *<sup>r</sup>* (*i*) and *γ<sup>i</sup>* <sup>2</sup> *ej* (*i*) are the variances of *<sup>r</sup>* (*i*) and *ej* (*i*), respectively over space. The *γ<sup>i</sup>* <sup>2</sup> *sj* (*i*) can be calculated by averaging the calculated variances for each observation. Then, the average spatial coefficient of determination, *Rs* 2 , which describes the proportion of the total variance explained by the time stability model of Eq. (3) can be written as

$$\left| R\_s^2 = \frac{\chi\_i^2 \left[ r(i) \right]}{\left\langle \chi\_i^2 \left[ s\_j(i) \right] \right\rangle} \right. \tag{6}$$

where *γ<sup>i</sup>* <sup>2</sup> *sj* (*i*) is the mean value of *γ<sup>i</sup>* <sup>2</sup> *sj* (*i*) over time. Given *ej* (*i*) with an average value of zero over time, *r*(*i*) can be approximated by time average of *sj* (*i*), *sj* (*i*) , i. e.,

$$r(i) = \left\langle s\_j(i) \right\rangle \tag{7}$$

The value of *Rs* <sup>2</sup> approaches one for a complete spatial pattern described by Eq. (3) and zero for a situation where no variances can be explained by the time stability model. Therefore, greater *Rs* 2 means stronger time stability of spatial pattern. However, Starr (2005) pointed out it is not suitable for the situation where no spatial variability of *SWC* and completely time stability exists.

To deal with this issue, temporal coefficient of variability, *CV <sup>t</sup>*, was developed to describe the degree of time stability (Starr, 2005) and can be written as

$$CV\_t = 100 \frac{1}{n} \sum\_{i=1}^{n} \frac{\sigma\_t \left[ s\_f(i) \right]}{\left\langle s\_f(i) \right\rangle} \tag{8}$$

where *z j'* (*i*) is the spatial anomaly at the same location *i* but time *j'*.

*Z*

and eigenvalues, *L*. Mathematically, they should satisfy the equation

where *E* is an *m*×*m* matrix that contains the eigenvectors as columns

*E*

respectively, where *m* is the number of sampling times.

variance*, V*, can be written as

and

ros at off-diagonals

To consider all the measurement times, the matrix of spatial anomalies, *Z*, and its spatial co‐

11 1

*z z*

*z z* æ ö ç ÷ <sup>=</sup> ç ÷ ç ÷ è ø K MO M L

<sup>1</sup> *<sup>T</sup> V ZZ*

The next step is to diagonalize the spatial covariance matrix *V* by finding its eigenvectors, *E*,

11 1

*e e*

*e e* æ ö ç ÷ <sup>=</sup> ç ÷ ç ÷ è ø L MO M L

*m mm*

and *L* is an *m*×*m* matrix that contains the associated eigenvalues along the diagonal and ze‐

<sup>11</sup> 0

æ ö ç ÷ <sup>=</sup> ç ÷ ç ÷ è ø L MO M L

0 *mm*

*l*

The eigenvectors in *E* represent the weights applied to each component in *V* to diagonalized *V*. This transformation is a rotation of the original axes in multi-dimensional space, where each dimension corresponds to a sampling time. Axes are orthogonal to each other and explain dif‐ ferent amount of covariance in the spatial anomaly dataset. The eigenvalues contained in *L* represent the variance that occurs in the direction of each newly identified axis. After diagonal‐

*m*

1

*l*

*L*

*n nm*

*m*

*<sup>n</sup>* <sup>=</sup> (12)

*VE LE* = (13)

(11)

53

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469

(14)

(15)

1

where *σ<sup>t</sup> sj* (*i*) is the standard deviation of *sj* (*i*) over time, *sj* (*i*) is the time average of *sj* (*i*). The smaller *CV <sup>t</sup>* value indicates stronger time stability. The *CV <sup>t</sup>* approaching zero in‐ dicates a perfect time stable pattern.

#### *3.1.4. Empirical orthogonal function method*

Empirical orthogonal function (EOF) analysis is a decomposition of a data set in terms of or‐ thogonal basis functions which are determined from the data. It is the same as performing a principal components analysis on the data, except that the EOF method finds both time ser‐ ies and spatial patterns. The basic principle of EOF method is to partition a series into timeinvariant spatial patterns (EOFs) of *SWC* and coefficients (ECs) which vary temporally but is constant spatially (Perry and Niemann, 2007). The original spatial series of *SWC* can be ob‐ tained by taking the sum of product of EOFs by ECs. A limited number of EOFs that present significant spatial variation of *SWC* can be selected and used to identify the dominant fac‐ tors determining the spatial pattern.

First, the spatial anomalies of *SWC* are computed. The spatial anomaly at location *i* and time *j*, *z <sup>j</sup>* (*i*), can be calculated as

$$z\_{\dot{f}}\left(\dot{i}\right) = \text{SWC}\_{\dot{f}}\left(\dot{i}\right) - \frac{1}{n} \sum\_{i=1}^{n} \text{SWC}\_{\dot{f}}\left(\dot{i}\right) \tag{9}$$

Considering the spatial anomalies, the spatial covariance *v <sup>j</sup>*, *j'* at time *j* and *j'* can be calcu‐ lated as

$$w\_{j,j'} = \frac{1}{n} \sum\_{i=1}^{n} z\_j(i) z\_{j'}(i) \tag{10}$$

where *z j'* (*i*) is the spatial anomaly at the same location *i* but time *j'*.

To consider all the measurement times, the matrix of spatial anomalies, *Z*, and its spatial co‐ variance*, V*, can be written as

$$Z = \begin{pmatrix} z\_{11} & \mathsf{K} & z\_{1m} \\ \mathsf{M} & \mathsf{O} & \mathsf{M} \\ z\_{m1} & \mathsf{L} & z\_{nm} \end{pmatrix} \tag{11}$$

and

The value of *Rs*

2

52 Advances in Agrophysical Research

time stability exists.

dicates a perfect time stable pattern.

tors determining the spatial pattern.

*j*, *z <sup>j</sup>* (*i*), can be calculated as

lated as

*3.1.4. Empirical orthogonal function method*

greater *Rs*

<sup>2</sup> approaches one for a complete spatial pattern described by Eq. (3) and zero

means stronger time stability of spatial pattern. However, Starr (2005) pointed

for a situation where no variances can be explained by the time stability model. Therefore,

out it is not suitable for the situation where no spatial variability of *SWC* and completely

To deal with this issue, temporal coefficient of variability, *CV <sup>t</sup>*, was developed to describe

*<sup>n</sup> t j*

1 ( ) <sup>1</sup> <sup>100</sup>

=

*i j*

where *σ<sup>t</sup> sj* (*i*) is the standard deviation of *sj* (*i*) over time, *sj* (*i*) is the time average of *sj* (*i*). The smaller *CV <sup>t</sup>* value indicates stronger time stability. The *CV <sup>t</sup>* approaching zero in‐

Empirical orthogonal function (EOF) analysis is a decomposition of a data set in terms of or‐ thogonal basis functions which are determined from the data. It is the same as performing a principal components analysis on the data, except that the EOF method finds both time ser‐ ies and spatial patterns. The basic principle of EOF method is to partition a series into timeinvariant spatial patterns (EOFs) of *SWC* and coefficients (ECs) which vary temporally but is constant spatially (Perry and Niemann, 2007). The original spatial series of *SWC* can be ob‐ tained by taking the sum of product of EOFs by ECs. A limited number of EOFs that present significant spatial variation of *SWC* can be selected and used to identify the dominant fac‐

First, the spatial anomalies of *SWC* are computed. The spatial anomaly at location *i* and time

1 *<sup>n</sup>*

*i*

*n* <sup>=</sup>

( ) ( ) ( )

' ' ( ) ( ) , <sup>1</sup> 1 *<sup>n</sup> <sup>j</sup> j j <sup>j</sup> <sup>i</sup> v z iz i n* <sup>=</sup>

*jj j*

Considering the spatial anomalies, the spatial covariance *v <sup>j</sup>*, *j'*

*z i SWC i SWC i*

1

= - å (9)

<sup>=</sup> å (10)

at time *j* and *j'* can be calcu‐

*n s i* s

( )

*s i*

é ù

ë û <sup>=</sup> å (8)

the degree of time stability (Starr, 2005) and can be written as

*t*

*CV*

$$V = \frac{1}{n} Z^T Z$$

respectively, where *m* is the number of sampling times.

The next step is to diagonalize the spatial covariance matrix *V* by finding its eigenvectors, *E*, and eigenvalues, *L*. Mathematically, they should satisfy the equation

$$V\mathbf{E} = LE\tag{13}$$

where *E* is an *m*×*m* matrix that contains the eigenvectors as columns

$$E = \begin{pmatrix} e\_{11} & \mathsf{L} & e\_{1m} \\ \mathsf{M} & \mathsf{O} & \mathsf{M} \\ e\_{m1} & \mathsf{L} & e\_{mm} \end{pmatrix} \tag{14}$$

and *L* is an *m*×*m* matrix that contains the associated eigenvalues along the diagonal and ze‐ ros at off-diagonals

$$L = \begin{pmatrix} l\_{11} & \mathbb{L} & 0 \\ \mathbf{M} & \mathbf{O} & \mathbf{M} \\ 0 & \mathbb{L} & l\_{mm} \end{pmatrix} \tag{15}$$

The eigenvectors in *E* represent the weights applied to each component in *V* to diagonalized *V*. This transformation is a rotation of the original axes in multi-dimensional space, where each dimension corresponds to a sampling time. Axes are orthogonal to each other and explain dif‐ ferent amount of covariance in the spatial anomaly dataset. The eigenvalues contained in *L* represent the variance that occurs in the direction of each newly identified axis. After diagonal‐ ization of *V*, *E* and *L* are arranged accordingly to keep the eigenvalues in *L* sorted in a descend‐ ing order. Therefore, the first axis explains the most covariance in the spatial anomaly dataset, and the second axis explains the second most covariance and so on. Each axis represents a di‐ rection in the multidimensional space. The total variance of the spatial anomaly data is the sum of the diagonal values in *L* and is equal to the total variance of the original spatial anomalies. Therefore, the portion of the variance, *P <sup>j</sup>* , that the *j* th axis can explain is

$$P\_j = \frac{I\_{j\bar{j}}}{\sum\_{k=1}^{m} I\_{kk}} \tag{16}$$

where *f <sup>K</sup>* = *K* / *N* , *K* = 0, 1, 2, …, N/2 cycles *h*-1, and *h* is the smoothing coefficient determin‐ ing the degree of averaging of adjacent independent frequencies and the degrees of freedom for the individual spectral variance estimates. The value of *h* cannot be too large because power spectrum of a variable may incorporate more bias from smoothing of the spectra, al‐ though larger value of *h* will result in smaller variance of the estimate (Gómez-Plaza et al.,

The covariance between *SWC <sup>j</sup>* (*i*) and *SWC j'* (*i*) can be estimated as a function of spatial

*F <sup>j</sup>* ( *f <sup>K</sup>* +1) \*

 p

> p

= (22)

<sup>=</sup> é ù - å ë û (21)

<sup>=</sup> é ù - å ë û (20)

*F <sup>j</sup>* '( *f <sup>K</sup>* +1) (19)

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 55

*N* -1 ∑ *l*=-*h h*

( ) ( ) ( ) ( )

' ' ( ) ( ) ( ) ( )

( ) ( )

*R f V fV f*

, 2 , 22 *<sup>K</sup>*

cos 2 sin 2 *<sup>K</sup>*

*F f SWC i f i i f i*

*j j K K*

p

p

Spectral coherency function *<sup>R</sup> <sup>j</sup>*, *j'* <sup>2</sup> ( *<sup>f</sup> <sup>K</sup>* ) can then be calculated in a similar manner to calculate

'

*V f*

The coherency function *R <sup>j</sup>*, *j'*( *f <sup>K</sup>* ) estimates the proportion of the spatial variance of *SWC <sup>j</sup>* (*i*) which can be explained by the spatial variance of *SWC j'* (*i*), as a function of spa‐ tial scale. Thus, spectral coherency is useful to measure time stability of soil water content as

The significance test method for spectral coherency consist of parametric methods based on an assumed theoretical distribution and nonparametric method such as reshuffling, boot‐ strapping, and bagging method. We encourage readers to see Si (2008) for a more detailed

*j j K*

*jK K j*

( ) ( )

'

2

cos 2 sin 2 *<sup>K</sup>*

*F f SWC i f i i f i*

*j jK K*

scale by calculating sample cross spectrum *V <sup>j</sup>*, *j'*( *f <sup>K</sup>* )

where the asterisk (\*) is a complex conjugate of

*<sup>V</sup> <sup>j</sup>*, *<sup>j</sup>* '( *<sup>f</sup> <sup>K</sup>* ) <sup>=</sup> (2*<sup>h</sup>* <sup>+</sup> 1) -1

1


*N*

*i*

*N*

*i*

=

where *f <sup>K</sup>* = *K* / *N* , *K* = 0, 1, 2, …, N/2 cycles *h*-1, and *i* = (-1)0.5.

coefficient of determination (Kachanoski and De Jong, 1988)

'

*j j*

a function of spatial scale.

account on significance test.

=

0

1


0

2000).

The relative importance of each axis in explaining the variability of spatial anomalies can be defined by the relative magnitude of the associated eigenvalue in *L*.

Each spatial anomaly can be described in terms of the new variable axis. A matrix *F* contain‐ ing the coordinates of the spatial anomalies on the new axis can be obtained by projecting the spatial anomalies onto the rotated axis. Mathematically, this operation is

$$F = ZE \tag{17}$$

Each column in *F* is EOF, and corresponding column in *E* is expansion coefficient (EC). From the aspect of time stability application, usually a limited number of EOFs that explains significant amount of the spatial variability of the original SWC series is selected to present the underlying time-invariable spatial pattern. North et al. (1982) indicated a set of selection criterion. An EOF is considered significant provided the lower confidence limit of its eigen‐ value is larger than the upper confidence limit of the eigenvalue of the next most important EOF (North et al., 1982).

#### *3.1.5. Spectral coherency analysis*

Spectral coherency is used to measure the similarity between two spatial series in a frequen‐ cy domain, which can be converted to a spatial scale. This method assumes the spatial series to be stationary and linear. Therefore, it can be used to find the scale specific information on time stability of spatial pattern between different measurements.

Spectral coherency analysis involves calculation of the power spectrum of a variable *Vt* 2 ( *f <sup>K</sup>* ) which estimates the observed variance as a function of spatial scale.

$$\mathcal{V}\_t^2\left(f\_K\right) = \left(2h+1\right)^{-1}N^{-1}\sum\_{l=-h}^h \left[ \left(\sum\_{j=1}^N \text{SWC}\_t\left(j\right)\cos\left(2\pi f\_{K+1}j\right)\right)^2 + \left(\sum\_{j=1}^N \text{SWC}\_t\left(j\right)\sin\left(2\pi f\_{K+1}j\right)\right)^2 \right] \tag{18}$$

where *f <sup>K</sup>* = *K* / *N* , *K* = 0, 1, 2, …, N/2 cycles *h*-1, and *h* is the smoothing coefficient determin‐ ing the degree of averaging of adjacent independent frequencies and the degrees of freedom for the individual spectral variance estimates. The value of *h* cannot be too large because power spectrum of a variable may incorporate more bias from smoothing of the spectra, al‐ though larger value of *h* will result in smaller variance of the estimate (Gómez-Plaza et al., 2000).

The covariance between *SWC <sup>j</sup>* (*i*) and *SWC j'* (*i*) can be estimated as a function of spatial scale by calculating sample cross spectrum *V <sup>j</sup>*, *j'*( *f <sup>K</sup>* )

$$\left(\boldsymbol{V}\_{\boldsymbol{f}\_{j},\boldsymbol{f}\_{j}}\left(\boldsymbol{f}\_{k}\right)=\left(2h+1\right)^{-1}\boldsymbol{N}^{-1}\sum\_{l=h}^{h}\boldsymbol{F}\_{\boldsymbol{f}\_{j}}\left(\boldsymbol{f}\_{K+1}\right)^{\*}\boldsymbol{F}\_{\boldsymbol{f}\_{j}}\cdot\left(\boldsymbol{f}\_{K+1}\right)\tag{19}$$

where the asterisk (\*) is a complex conjugate of

ization of *V*, *E* and *L* are arranged accordingly to keep the eigenvalues in *L* sorted in a descend‐ ing order. Therefore, the first axis explains the most covariance in the spatial anomaly dataset, and the second axis explains the second most covariance and so on. Each axis represents a di‐ rection in the multidimensional space. The total variance of the spatial anomaly data is the sum of the diagonal values in *L* and is equal to the total variance of the original spatial anomalies.

, that the *j*

*jj j m*

*l*

*P*

the spatial anomalies onto the rotated axis. Mathematically, this operation is

defined by the relative magnitude of the associated eigenvalue in *L*.

time stability of spatial pattern between different measurements.

which estimates the observed variance as a function of spatial scale.

=

1

*kk k*

The relative importance of each axis in explaining the variability of spatial anomalies can be

Each spatial anomaly can be described in terms of the new variable axis. A matrix *F* contain‐ ing the coordinates of the spatial anomalies on the new axis can be obtained by projecting

Each column in *F* is EOF, and corresponding column in *E* is expansion coefficient (EC). From the aspect of time stability application, usually a limited number of EOFs that explains significant amount of the spatial variability of the original SWC series is selected to present the underlying time-invariable spatial pattern. North et al. (1982) indicated a set of selection criterion. An EOF is considered significant provided the lower confidence limit of its eigen‐ value is larger than the upper confidence limit of the eigenvalue of the next most important

Spectral coherency is used to measure the similarity between two spatial series in a frequen‐ cy domain, which can be converted to a spatial scale. This method assumes the spatial series to be stationary and linear. Therefore, it can be used to find the scale specific information on

Spectral coherency analysis involves calculation of the power spectrum of a variable *Vt*

( ) ( ) ( ) ( ) ( ) ( )

=- = =

1 1 2 1 cos 2 sin 2 *h N N t K t K t K lh j j V f h N SWC j f j SWC j f j*

p

é ù <sup>æ</sup> ö æ <sup>ö</sup> ê ú = + <sup>ç</sup> ÷ ç <sup>+</sup> <sup>÷</sup> ê ú <sup>ç</sup> ÷ ç <sup>÷</sup> <sup>è</sup> ø è <sup>ø</sup> ê ú ë û


*l* =

th axis can explain is

å (16)

*F ZE* = (17)

2 2

 p

1 1

+ +

å å å (18)

2 ( *f <sup>K</sup>* )

Therefore, the portion of the variance, *P <sup>j</sup>*

54 Advances in Agrophysical Research

EOF (North et al., 1982).

1 2 1

*3.1.5. Spectral coherency analysis*

$$F\_j(f\_{\chi}) = \sum\_{i=0}^{N-1} \text{SWC}\_j(i) \left[ \cos(2\pi f\_K i) - i \sin(2\pi f\_K i) \right] \tag{20}$$

$$F\_{\dot{f}}\left(f\_{\kappa}\right) = \sum\_{i=0}^{N-1} \text{SWC}\_{\dot{f}}\left(i\right) \left[\cos\left(2\pi f\_{k}i\right) - i\sin\left(2\pi f\_{k}i\right)\right] \tag{21}$$

where *f <sup>K</sup>* = *K* / *N* , *K* = 0, 1, 2, …, N/2 cycles *h*-1, and *i* = (-1)0.5.

Spectral coherency function *<sup>R</sup> <sup>j</sup>*, *j'* <sup>2</sup> ( *<sup>f</sup> <sup>K</sup>* ) can then be calculated in a similar manner to calculate coefficient of determination (Kachanoski and De Jong, 1988)

$$\left| \mathcal{R}\_{j,j} \right|^2 \left( f\_{\times} \right) = \frac{\left| V\_{j,j} \left( f\_{\times} \right) \right|^2}{V\_j^2 \left( f\_{\times} \right) V\_j^{-2} \left( f\_{\times} \right)} \tag{22}$$

The coherency function *R <sup>j</sup>*, *j'*( *f <sup>K</sup>* ) estimates the proportion of the spatial variance of *SWC <sup>j</sup>* (*i*) which can be explained by the spatial variance of *SWC j'* (*i*), as a function of spa‐ tial scale. Thus, spectral coherency is useful to measure time stability of soil water content as a function of spatial scale.

The significance test method for spectral coherency consist of parametric methods based on an assumed theoretical distribution and nonparametric method such as reshuffling, boot‐ strapping, and bagging method. We encourage readers to see Si (2008) for a more detailed account on significance test.

#### *3.1.6. Wavelet coherency analysis*

Wavelet analysis, differing with spectral analysis, can identify localized features of soil proc‐ esses (Si, 2008). It is suitable to reveal the scale and location specific time-persistence of spa‐ tial pattern of *SWC* between sampling occasions assuming the soil water content system is linear (Biswas and Si, 2011c). Wavelet coherency of two spatial series can describe the time stability of spatial pattern. It requires the calculation of wavelet coefficients for each of the two data series and associated cross-wavelet spectrum. Many publications on the introduc‐ tion of wavelet analysis (Farge, 1992; Kumar and Foufoula-Georgiou, 1993, 1997) and wave‐ let coherency (Biswas and Si, 2011c; Grinsted et al., 2004; Si, 2008; Si and Zeleke, 2005) can be found. Here, we will present the basic procedure to calculate wavelet coherency.

First, wavelet coefficient, *Wi <sup>Y</sup>* (*s*), is calculated with the continuous wavelet transform (CWT) for a *SWC* series of length *n (Y*<sup>i</sup> , *i* = 1, 2, …, *n*) with equal incremental distance *δx*. This can be defined as the convolution of *Y*<sup>i</sup> with the scaled and normalized wavelet using fast Fourier transform (Si and Zeleke, 2005; Torrence and Compo, 1998):

$$\mathcal{W}\_i^Y\left(s\right) = \sqrt{\frac{\delta \mathbf{x}}{s}} \sum\_{j=1}^n Y\_j \mathcal{W}\left[ (j-i)\frac{\delta \mathbf{x}}{s} \right] \tag{23}$$

( ) ( ) ( ) *YZ Y Z W s Ws Ws i ii* <sup>=</sup> é ùé ù ë ûë û (25)

sur sur (26)

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469

2

*ω* 2), where *ω* is the frequency. Eq. (29)

2

å (29)

 w*<sup>s</sup>* é ù <sup>=</sup> é ù - ë ûê ú ë û (30)

*x*

t

(27)

57

(28)

and the wavelet coherency of two spatial series can be written as

( )

*s*

2

*i*

r=

where ( <sup>⋅</sup> ) <sup>↔</sup> is a smoothing operator, which can be written as

has a similar footprint as the Morlet wavelet

sur

( )

*YZ i*

sur

( ) *W SM SM W scale space* <sup>=</sup> é ù ë û

1

1

=

*k*

*SM W s W s*

The Fourier transform of Equation (28) is exp(-2*s* <sup>2</sup>

convolution theorem and is written as

The smoothing along scales can be written as:

t

p

exp *s* 2 2*s*

where *τ* denotes location. Therefore, the smoothing along locations can be written as

( ) ( ) ( )

*<sup>m</sup> <sup>k</sup>*

ì ü é ù ï ï - <sup>=</sup> í ý ê ú ê ú ï ï î þ ë û

p

t

*s s*

<sup>1</sup> , , , exp <sup>2</sup> <sup>2</sup>

 t

( ) { ( ) ( ) } 2 2 , , , exp 2 *space SM W s x IFFT FFT W s*

*space s*

can be implemented using Fast Fourier Transform (FFT) and Inverse FFT (IFFT) based on

where *SM space* denotes smoothing along the wavelet scale axis and *SM scale* smoothing in spa‐ tial distance. The following smoothing function is the normalized real Morlet wavelet and

> 2 2

t

æ ö ç ÷ ç ÷ è ø

*W s*

*Y Z i i*

*WsWs*

2

2 2

( ) ( )

where *ψ* is the mother wavelet function and *s* is the scale. The mother wavelet functions include Morlet, Mexican hat, Harr, and others (Si, 2008). Depending on the purpose, differ‐ ent mother wavelet functions can be selected. The Morlet wavelets allow us to detect both location dependent amplitude and phase for different frequencies in the spatial series (Tor‐ rence and Compo, 1998), which can be written as

$$
\sqrt{\left[\eta\right]} = \pi^{-1/4} e^{i\alpha\eta - 0.5\eta^2} \tag{24}
$$

where *ω* is dimensionless frequency and *η* is dimensionlessspace. The wavelet is stretched in space (*x*) by varying its scale (*s*), so that *η* =*s* / *x*. The Morlet wavelet (with *ω* =6) is a good choice for feature extraction purpose like identifying the scales and locations, since it pro‐ vides a good balance between space and frequency localization.

The wavelet coefficients *Wi <sup>Y</sup>* (*s*) can be expressed as *a* + *ib* where *a* and *b* are the real and imaginary components of *Wi <sup>Y</sup>* (*s*). For the polar form of complex numbers, *Wi <sup>Y</sup>* (*s*)= <sup>|</sup>*Wi <sup>Y</sup>* (*s*)|(*cosθ* + *isinθ*), where *θ* =*arctan b <sup>a</sup>* is called the phase or argument of *Wi <sup>Y</sup>* (*s*). The wavelet power spectrum is defined as |*Wi <sup>Y</sup>* (*s*)|<sup>2</sup> and the local phase is defined as the complex argument of *Wi <sup>Y</sup>* (*s*).

After calculating the wavelet spectra *Wi <sup>Y</sup>* (*s*) and *Wi <sup>Z</sup>* (*s*) corresponding to two *SWC* spatial series Y and Z, respectively for two different times, cross wavelet power spectrum *Wi YZ* (*s*) at scale *s* and location *i* can be defined as

$$\mathbf{W}\_{i}^{YZ}\left(\mathbf{s}\right) = \left[\mathbf{W}\_{i}^{Y}\left(\mathbf{s}\right)\right]\overline{\left[\mathbf{W}\_{i}^{Z}\left(\mathbf{s}\right)\right]}\tag{25}$$

and the wavelet coherency of two spatial series can be written as

*3.1.6. Wavelet coherency analysis*

56 Advances in Agrophysical Research

First, wavelet coefficient, *Wi*

The wavelet coefficients *Wi*

complex argument of *Wi*

*Wi*

*<sup>Y</sup>* (*s*)= <sup>|</sup>*Wi*

imaginary components of *Wi*

for a *SWC* series of length *n (Y*<sup>i</sup>

defined as the convolution of *Y*<sup>i</sup>

Wavelet analysis, differing with spectral analysis, can identify localized features of soil proc‐ esses (Si, 2008). It is suitable to reveal the scale and location specific time-persistence of spa‐ tial pattern of *SWC* between sampling occasions assuming the soil water content system is linear (Biswas and Si, 2011c). Wavelet coherency of two spatial series can describe the time stability of spatial pattern. It requires the calculation of wavelet coefficients for each of the two data series and associated cross-wavelet spectrum. Many publications on the introduc‐ tion of wavelet analysis (Farge, 1992; Kumar and Foufoula-Georgiou, 1993, 1997) and wave‐ let coherency (Biswas and Si, 2011c; Grinsted et al., 2004; Si, 2008; Si and Zeleke, 2005) can be

*<sup>Y</sup>* (*s*), is calculated with the continuous wavelet transform (CWT)

, *i* = 1, 2, …, *n*) with equal incremental distance *δx*. This can be

with the scaled and normalized wavelet using fast Fourier

ë û <sup>å</sup> (23)


*<sup>Y</sup>* (*s*) can be expressed as *a* + *ib* where *a* and *b* are the real and

*<sup>Y</sup>* (*s*). For the polar form of complex numbers,

*<sup>a</sup>* is called the phase or argument of *Wi*

and the local phase is defined as the

*<sup>Z</sup>* (*s*) corresponding to two *SWC* spatial

*<sup>Y</sup>* (*s*).

*YZ* (*s*)

found. Here, we will present the basic procedure to calculate wavelet coherency.

( ) ( )

=

*j j x x W s Y ji*

y

where *ψ* is the mother wavelet function and *s* is the scale. The mother wavelet functions include Morlet, Mexican hat, Harr, and others (Si, 2008). Depending on the purpose, differ‐ ent mother wavelet functions can be selected. The Morlet wavelets allow us to detect both location dependent amplitude and phase for different frequencies in the spatial series (Tor‐

> <sup>2</sup> 1/4 0.5 *i e* wh h

where *ω* is dimensionless frequency and *η* is dimensionlessspace. The wavelet is stretched in space (*x*) by varying its scale (*s*), so that *η* =*s* / *x*. The Morlet wavelet (with *ω* =6) is a good choice for feature extraction purpose like identifying the scales and locations, since it pro‐

*b*

*<sup>Y</sup>* (*s*) and *Wi*

series Y and Z, respectively for two different times, cross wavelet power spectrum *Wi*

*<sup>Y</sup>* (*s*)|<sup>2</sup>

é ù <sup>=</sup> - ê ú

*s s*

 d

<sup>1</sup> *<sup>i</sup> <sup>n</sup> <sup>Y</sup>*

yh

vides a good balance between space and frequency localization.

*<sup>Y</sup>* (*s*)|(*cosθ* + *isinθ*), where *θ* =*arctan*

The wavelet power spectrum is defined as |*Wi*

After calculating the wavelet spectra *Wi*

at scale *s* and location *i* can be defined as

*<sup>Y</sup>* (*s*).

 p

d

transform (Si and Zeleke, 2005; Torrence and Compo, 1998):

rence and Compo, 1998), which can be written as

$$\rho\_i^2(s) = \frac{\begin{vmatrix} \mathfrak{su}\_{YZ} \\ W\_i \end{vmatrix}^2(s)}{\begin{vmatrix} \mathfrak{su}\_Y \\ W\_i \end{vmatrix}^2(s) \Big| \begin{vmatrix} \mathfrak{su}\_Z \\ W\_i \end{vmatrix}^2(s)} \tag{26}$$

where ( <sup>⋅</sup> ) <sup>↔</sup> is a smoothing operator, which can be written as

$$\stackrel{\text{SSI}}{\text{W}} = \text{SM}\_{\text{scale}} \left[ \text{SM}\_{\text{space}} \left( \text{W} \right) \right] \tag{27}$$

where *SM space* denotes smoothing along the wavelet scale axis and *SM scale* smoothing in spa‐ tial distance. The following smoothing function is the normalized real Morlet wavelet and has a similar footprint as the Morlet wavelet

$$\frac{1}{s\sqrt{2\pi}}\exp\left(-\frac{\tau^2}{2s^2}\right) \tag{28}$$

where *τ* denotes location. Therefore, the smoothing along locations can be written as

$$\left| SM\_{\text{space}}\left(W, s, \tau\right) = \sum\_{k=1}^{m} \left| \mathcal{W}\left(s, \tau\right) \frac{1}{s\sqrt{2\pi}} \exp\left[-\frac{\left(\tau - \chi\_{k}\right)^{2}}{2s^{2}}\right] \right|\right|\_{s} \tag{29}$$

The Fourier transform of Equation (28) is exp(-2*s* <sup>2</sup> *ω* 2), where *ω* is the frequency. Eq. (29) can be implemented using Fast Fourier Transform (FFT) and Inverse FFT (IFFT) based on convolution theorem and is written as

$$\text{SSM}\_{\text{space}}\left(\mathcal{W}, \mathbf{s}, \mathbf{x}\right) = \text{IFFT}\left(\text{FFT}\left[\mathcal{W}\left(\mathbf{s}, \mathbf{r}\right)\right]\right|\exp\left(-2s^2\alpha^2\right)\right) \tag{30}$$

The smoothing along scales can be written as:

$$\text{SM}\_{\text{scale}}\left(\mathcal{W}, \mathbf{s}\_{k'}, \mathbf{x}\right) = \frac{1}{2p + 1} \sum\_{j=k-p}^{k+p} \left[ \text{SM}\_{\text{space}}\left(\mathcal{W}, \mathbf{s}\_{j'}, \mathbf{x}\right) \Pi\left(\mathbf{0}.\mathbf{6s}\_{j}\right) \right] \mathbf{x} \tag{31}$$

**2.** Calculate a projection, denoted by *p*

**3.** Find the spatial instants *si*

*θk* , v(*si*

reference (Huang et al., 2003).

**3.2. Time stability at points**

*<sup>θ</sup><sup>k</sup>* (*s*)}*<sup>k</sup>* =1 *K* .

sets considered.

sets *p*

**4.** Interpolate *si*

whole set of direction vectors), giving *p*

*θk*

*<sup>θ</sup><sup>k</sup>* (*t*)}*<sup>s</sup>*=1

*<sup>θ</sup><sup>k</sup>* (*s*)}*<sup>k</sup>* =1

*<sup>θ</sup><sup>k</sup>* ) to obtain multivariate envelope curves e

**5.** For a set of *K* direction vectors, the mean *m*(*s*) of the envelope curves is calculated as:

**6.** Extract the "detail" *d*(*s*) using *d*(*s*)=v(*s*) - *m*(*s*). If the "detail" *d*(*s*) fulfills the stoppage criterion for a multivariate IMF, apply the above procedure to v (*s*) - *d*(s), otherwise ap‐ ply it to *d*(*s*). The stoppage criterion for multivariate IMFs is similar to that proposed by

MEMD has the ability to align "common scales" present within multivariate data. Each com‐ mon scale is manifested in the common oscillatory modes in all the variates within an *m*variate IMF. After MEMD analysis, scale and location specific time stability of *SWC* can be easily identified by comparing and calculating spearman's rank correlation coefficient be‐ tween the IMFs with the same numerical numbers for different measurements using Eq. (2). The exact scale for each IMF can be obtained from the instantaneous frequencies by Hilbert transformation with IMF (Huang et al., 1998). The instantaneous frequencies of soil water content can be converted to period (1/frequency), which was further converted to the spatial scale after multiplying the period with the sampling interval. MEMD can be completed us‐ ing Matlab program such as that written by Rehman and Mandic (2009) (http:// www.commsp.ee.ic.ac uk/~mandic/research/emd.htm). This method requires that soil water

Time stability at point scales are usually obtained in terms of average *SWC* estimation quali‐ ty. Direct and indirect methods can be used to estimate average *SWC* by the measurement at time stable locations. The direct method estimates average *SWC* directly by measuring *SWC* at a time-stable location (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Pachepsky et al., 2005). The indirect method estimates average *SWC* by considering the off‐ set between the average *SWC* and the measurement value at a time-stable location (Grayson and Western, 1998; Starks et al., 2006). All the indices are listed in Table 1. Each index has its own advantages and disadvantages, and these are discussed in the application (Section 4.2). For all indices, smaller value of time stability index at a given location indicates stronger

time stability and also higher quality of average *SWC* estimation by that location.

*<sup>m</sup>*(*s*)= <sup>1</sup> *<sup>K</sup>* ∑ *k*=1 *K* e

content measurement should be obtained in intervals with equidistance.

number of sampling point for each time) along the direction vector x

<sup>S</sup> , of the spatial data sets {v(*s*)}

*<sup>K</sup>* as the set of projections.

corresponding to the maxima of the set of projected data

*s*=1

*θk*

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469

*<sup>θ</sup><sup>k</sup>* (*s*)}*<sup>k</sup>* =1

*<sup>θ</sup><sup>k</sup>* (*s*) (32)

<sup>S</sup> (*S* is the

59

, for all *k* ( the

*<sup>K</sup>* for all the data

where *P* is the number of terms on each symmetrical half of the window, and П is the rec‐ tangle function. The factor of 0.6 is the empirically determined scale decorrelation length for the Morlet wavelet (Si and Zeleke, 2005; Torrence and Compo, 1998).

There are many methods available to test the statistical significance of wavelet coherency. The Monte Carlo simulation or reshuffling method is among the suggested ones. A detailed description of these methods can be found in Si (2008). Matlab codes for calculating wavelet coherency are available at URL: http://www.pol.ac.uk/home/research/waveletcoherence/ (Grinsted et al., 2004).

## *3.1.7. Multivariate empirical mode decomposition*

Empirical mode decomposition (EMD), extracts oscillations from the soil water content ser‐ ies into a finite and often small number of intrinsic mode functions (IMFs) according to the energy associated with different space scales (Huang et al., 1998; Huang and Wu, 2008). Un‐ like spectral and wavelet methods, EMD does not call for any assumption of the data and works directly in the spatial domain with the basis completely derived from the data. There‐ fore, it is intuitive, direct, *a posteriori* and adaptive (Huang et al., 1998; Huang and Wu, 2008). The locality and adaptivity of EMD can deal with different types of spatial series including non-stationary and nonlinear (Huang et al., 1998). Therefore, EMD has a great potential to find the underlying scale of spatial series of soil moisture without imposing any mathemati‐ cal assumption on the measured data (Biswas and Si, 2011b).

Multivariate empirical mode decomposition (MEMD) is the multivariate extensions of standard EMD (Huang et al., 1998). An important step for MEMD is the computation of the local mean because local extrema are not well defined directly for multivariate spatial data. Moreover, the notion of "oscillatory model" defining an IMF is rather confusing for multi‐ variate spatial data. To deal with these problems, Rehman and Mandic (2010) produced a multiple *m*-dimensional envelopes by taking projections of multiple inputs along different directions in an *m*-dimensional space.

Assuming {v(*s*)} *s*=1 <sup>S</sup> ={υ<sup>1</sup> (*s*), υ<sup>2</sup> (*s*), …, υm(s)} being the *m* spatial data sets of soil water content as a function of space (*s*), and x *<sup>θ</sup><sup>k</sup>* ={x1 *k* , x2 *k* , …, xm *<sup>k</sup>* } denoting a set of direction vectors along the directions given by angles *<sup>θ</sup> <sup>k</sup>* ={*θ*<sup>1</sup> *k* , *θ*<sup>2</sup> *k* , …, *θ*m-1 *<sup>k</sup>* } on a *m*-1-dimensional sphere (*k* is the number of direction used to calculate the projections and envelope curves). Then, IMFs of the *m* spatial data sets can be obtained by MEMD as follows:

**1.** Choose a suitable point set for sampling on an *m*-1-dimensional sphere. This can be done by sampling unit hyperspheres (*m*-spheres) based on both uniform angular sam‐ pling methods and quasi-Monte Carlo-based low-discrepancy sequences.


$$\mathbf{m(s)} = \frac{1}{K} \sum\_{k=1}^{K} \mathbf{e}^{\theta\_k}(\mathbf{s}) \tag{32}$$

**6.** Extract the "detail" *d*(*s*) using *d*(*s*)=v(*s*) - *m*(*s*). If the "detail" *d*(*s*) fulfills the stoppage criterion for a multivariate IMF, apply the above procedure to v (*s*) - *d*(s), otherwise ap‐ ply it to *d*(*s*). The stoppage criterion for multivariate IMFs is similar to that proposed by reference (Huang et al., 2003).

MEMD has the ability to align "common scales" present within multivariate data. Each com‐ mon scale is manifested in the common oscillatory modes in all the variates within an *m*variate IMF. After MEMD analysis, scale and location specific time stability of *SWC* can be easily identified by comparing and calculating spearman's rank correlation coefficient be‐ tween the IMFs with the same numerical numbers for different measurements using Eq. (2). The exact scale for each IMF can be obtained from the instantaneous frequencies by Hilbert transformation with IMF (Huang et al., 1998). The instantaneous frequencies of soil water content can be converted to period (1/frequency), which was further converted to the spatial scale after multiplying the period with the sampling interval. MEMD can be completed us‐ ing Matlab program such as that written by Rehman and Mandic (2009) (http:// www.commsp.ee.ic.ac uk/~mandic/research/emd.htm). This method requires that soil water content measurement should be obtained in intervals with equidistance.

#### **3.2. Time stability at points**

( ) ( ) ( ) <sup>1</sup> , , , , 0.6 2 1

where *P* is the number of terms on each symmetrical half of the window, and П is the rec‐ tangle function. The factor of 0.6 is the empirically determined scale decorrelation length for

There are many methods available to test the statistical significance of wavelet coherency. The Monte Carlo simulation or reshuffling method is among the suggested ones. A detailed description of these methods can be found in Si (2008). Matlab codes for calculating wavelet coherency are available at URL: http://www.pol.ac.uk/home/research/waveletcoherence/

Empirical mode decomposition (EMD), extracts oscillations from the soil water content ser‐ ies into a finite and often small number of intrinsic mode functions (IMFs) according to the energy associated with different space scales (Huang et al., 1998; Huang and Wu, 2008). Un‐ like spectral and wavelet methods, EMD does not call for any assumption of the data and works directly in the spatial domain with the basis completely derived from the data. There‐ fore, it is intuitive, direct, *a posteriori* and adaptive (Huang et al., 1998; Huang and Wu, 2008). The locality and adaptivity of EMD can deal with different types of spatial series including non-stationary and nonlinear (Huang et al., 1998). Therefore, EMD has a great potential to find the underlying scale of spatial series of soil moisture without imposing any mathemati‐

Multivariate empirical mode decomposition (MEMD) is the multivariate extensions of standard EMD (Huang et al., 1998). An important step for MEMD is the computation of the local mean because local extrema are not well defined directly for multivariate spatial data. Moreover, the notion of "oscillatory model" defining an IMF is rather confusing for multi‐ variate spatial data. To deal with these problems, Rehman and Mandic (2010) produced a multiple *m*-dimensional envelopes by taking projections of multiple inputs along different

, …, xm

number of direction used to calculate the projections and envelope curves). Then, IMFs of

**1.** Choose a suitable point set for sampling on an *m*-1-dimensional sphere. This can be done by sampling unit hyperspheres (*m*-spheres) based on both uniform angular sam‐

pling methods and quasi-Monte Carlo-based low-discrepancy sequences.

, …, *θ*m-1

(*s*), …, υm(s)} being the *m* spatial data sets of soil water content

*<sup>k</sup>* } denoting a set of direction vectors along

*<sup>k</sup>* } on a *m*-1-dimensional sphere (*k* is the

= é ù Õ

<sup>+</sup> å ë û (31)

*scale k space j j x jkp SM W s x SM W s x s*

*k p*

+

= -

*p*

the Morlet wavelet (Si and Zeleke, 2005; Torrence and Compo, 1998).

(Grinsted et al., 2004).

58 Advances in Agrophysical Research

*3.1.7. Multivariate empirical mode decomposition*

directions in an *m*-dimensional space.

the directions given by angles *<sup>θ</sup> <sup>k</sup>* ={*θ*<sup>1</sup>

(*s*), υ<sup>2</sup>

the *m* spatial data sets can be obtained by MEMD as follows:

*<sup>θ</sup><sup>k</sup>* ={x1 *k* , x2 *k*

> *k* , *θ*<sup>2</sup> *k*

*s*=1 <sup>S</sup> ={υ<sup>1</sup>

as a function of space (*s*), and x

Assuming {v(*s*)}

cal assumption on the measured data (Biswas and Si, 2011b).

Time stability at point scales are usually obtained in terms of average *SWC* estimation quali‐ ty. Direct and indirect methods can be used to estimate average *SWC* by the measurement at time stable locations. The direct method estimates average *SWC* directly by measuring *SWC* at a time-stable location (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Pachepsky et al., 2005). The indirect method estimates average *SWC* by considering the off‐ set between the average *SWC* and the measurement value at a time-stable location (Grayson and Western, 1998; Starks et al., 2006). All the indices are listed in Table 1. Each index has its own advantages and disadvantages, and these are discussed in the application (Section 4.2). For all indices, smaller value of time stability index at a given location indicates stronger time stability and also higher quality of average *SWC* estimation by that location.


confirmed the presence of time stability of the spatial pattern (Biswas and Si, 2011c; Hu et al., 2009, 2010b; Martínez-Fernández and Ceballos, 2003; Vachaud et al., 1985; Vivoni et al., 2008). In the first of this kind, Vachaud et al. (1985) measured soil water storage of 0 to 1.0 m depth at 17 neutron access tubes located at 10 m intervals for 24 occasions over a time span of two and a half years. According to the cumulative probability function for the two ex‐ treme conditions, i. e., the wettest on 7 September 1982, and the driest on 17 August 1983, many locations maintained their ranks (Figure 1). For example, location #1 always had the maximum soil water storage, location #4 had the minimum soil water storage irrespective of the environmental conditions (e. g., wet or dry). Spearman's rank correlation coefficients of the soil water storage between measurements were also calculated to examine the degree of time stability of the spatial pattern. As an example, the correlation between seven measure‐ ments representing the range of the water storages recorded over the entire observation pe‐ riod is presented in Table 2. All rank correlation coefficients were highly significant at the 0.1% two-tailed test, indicating strong time stability of the spatial pattern of soil water stor‐ age. On the other hand, time instability of spatial pattern between measurements was also reported (Comegna and Basile, 1994; Mohanty and Skaggs, 2001). Although it is very popu‐ lar among researchers, Spearman's rank correlation analysis should be viewed only as a statistical tool measuring the degree of concordance between two rankings (Vachaud et al., 1985). It may be questionable if differences between measured values are smaller than ex‐ perimental uncertainties. This may be the case in situations where either the probability

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 61

density function is very uniform or the experimental determinations are very crude.

**Figure 1.** Cumulative probability function of soil water storage at 0-1.0 m measured in the most dry, and the most wet

conditions. Numbers refer to measuring locations. (Reproduced from Vachaud et al. (1985)).

†Direct refers to estimating average *SWC* directly by measuring *SWC* at a time-stable location and indirect refers to estimating average *SWC* by considering the offset between the mean value and the measurement value at a timestable location.

‡ σ(β*<sup>i</sup>* )=*SDRD* (*i*).

**Table 1.** Time stability indices at point scales and their formulas

## **4. Application**

#### **4.1. Time stabilityof spatial pattern**

Information on the spatial pattern of soil water content with time is important to understand the changes in hydrological processes and in making predictions with physical and statisti‐ cal models. In this section, we have documented the use of each method introduced above to characterize the similarity of the spatial pattern of *SWC*.

Spearman's rank correlation analysis is the most widely used method to examine the simi‐ larity of the spatial patterns of *SWC*. Numerous studies have used this method and have confirmed the presence of time stability of the spatial pattern (Biswas and Si, 2011c; Hu et al., 2009, 2010b; Martínez-Fernández and Ceballos, 2003; Vachaud et al., 1985; Vivoni et al., 2008). In the first of this kind, Vachaud et al. (1985) measured soil water storage of 0 to 1.0 m depth at 17 neutron access tubes located at 10 m intervals for 24 occasions over a time span of two and a half years. According to the cumulative probability function for the two ex‐ treme conditions, i. e., the wettest on 7 September 1982, and the driest on 17 August 1983, many locations maintained their ranks (Figure 1). For example, location #1 always had the maximum soil water storage, location #4 had the minimum soil water storage irrespective of the environmental conditions (e. g., wet or dry). Spearman's rank correlation coefficients of the soil water storage between measurements were also calculated to examine the degree of time stability of the spatial pattern. As an example, the correlation between seven measure‐ ments representing the range of the water storages recorded over the entire observation pe‐ riod is presented in Table 2. All rank correlation coefficients were highly significant at the 0.1% two-tailed test, indicating strong time stability of the spatial pattern of soil water stor‐ age. On the other hand, time instability of spatial pattern between measurements was also reported (Comegna and Basile, 1994; Mohanty and Skaggs, 2001). Although it is very popu‐ lar among researchers, Spearman's rank correlation analysis should be viewed only as a statistical tool measuring the degree of concordance between two rankings (Vachaud et al., 1985). It may be questionable if differences between measured values are smaller than ex‐ perimental uncertainties. This may be the case in situations where either the probability density function is very uniform or the experimental determinations are very crude.

**Index Formula and explanation**

*<sup>m</sup>*−<sup>1</sup> ∑ *j*=1 *m* (δ*j*

(*i*)− *SW Cj SW Cj*

*<sup>P</sup>*=0.95 and β(*i*)

(*i*)− δ*<sup>i</sup>* )2 and

error, *RMSE*(*i*) Direct Jacobs et al.,

and

*<sup>P</sup>*=0.05 are the relative *SWC* values

observation time *j*. Indirect Guber et al.,

<sup>2</sup> + *SDRD*(*i*)2)1/2

at cumulative probability of 95% and 5%, respectively.

is the standard deviation of *SWC* at

error, *MABE*(*i*) Indirect Hu et al.,

†Direct refers to estimating average *SWC* directly by measuring *SWC* at a time-stable location and indirect refers to estimating average *SWC* by considering the offset between the mean value and the measurement value at a time-

Information on the spatial pattern of soil water content with time is important to understand the changes in hydrological processes and in making predictions with physical and statisti‐ cal models. In this section, we have documented the use of each method introduced above

Spearman's rank correlation analysis is the most widely used method to examine the simi‐ larity of the spatial patterns of *SWC*. Numerous studies have used this method and have

where *SDRD*(*i*)= <sup>1</sup>

where *RMSE*(*i*)=( δ*<sup>i</sup>*

δ*j* (*i*)= *SW Cj*

<sup>δ</sup>*<sup>i</sup>* <sup>=</sup> <sup>1</sup> *<sup>m</sup>* ∑ *j*=1 *m* δ*j* (*i*)

σ(β*<sup>i</sup>* )= <sup>1</sup> *<sup>m</sup>*−<sup>1</sup> ∑ *j*=1 *m* (β*j* (*i*)− β*<sup>i</sup>* )<sup>2</sup>

where β(*i*)

where σ *<sup>j</sup>*

**Table 1.** Time stability indices at point scales and their formulas

to characterize the similarity of the spatial pattern of *SWC*.

**4.1. Time stabilityof spatial pattern**

Standard deviation of relative difference, *SDRD*(*i*)

60 Advances in Agrophysical Research

Root mean square

Standard deviation of relative *SWC*,

Width of the 90% empirical tolerance interval of relative water content, *T* (*i*)

Chi-squared statistic, χ 2(*i*)

Root-meansquared differences, *D*(*i*)

stable location.

‡ σ(β*<sup>i</sup>*

Mean absolute bias

)=*SDRD* (*i*).

**4. Application**

σ(β*<sup>i</sup>* ) ‡ **Average SWC estimation Method†**

Direct and indirect

**Reference**

Vachaud et al., 1985

2004

al., 2005

2008

2008

2008

2010a, 2010b

Indirect Pachepsky et

Indirect Guber et al.,

Indirect Guber et al.,

**Figure 1.** Cumulative probability function of soil water storage at 0-1.0 m measured in the most dry, and the most wet conditions. Numbers refer to measuring locations. (Reproduced from Vachaud et al. (1985)).


EC values (Figures 2c and 2d). There was a large decrease in EC1 weight during the dry sea‐ son, which agrees well with the finding of Western et al. (1999) that soil moisture patterns on moderately wet occasions show a strong topographic influence. The EOF method can also be used to predict the *SWC* distribution for unobserved times if the average *SWC* and the tempo‐ ral varying coefficient (ECs), which reflect the relative importance of different spatial patterns can be obtained. Average *SWC* can be estimated by many methods such as linear regression (Western et al., 1999), a dynamic multiple linear regression of soil moisture against topograph‐ ic attributes (Wilson et al., 2005), and by measurement from representative location (Grayson and Western, 1998). Examples of EOF method can also be found in other study (Ibrahim and Huggins, 2011; Jawson and Niemann, 2007; Joshi and Mohanty, 2010; Korres et al., 2010; Yoo and Kim, 2004). Different numbers, usually one to four, of significant EOFs were obtained to explain most of *SWC* variations (usually more than 70%). The time-invariant spatial patterns were usually related to topography-related factors (Yoo and Kim, 2004), sand content (Jawson and Niemann, 2007), mixed influences of rainfall, topography, and soil texture (Joshi and Mo‐ hanty, 2010), elevation and wetness index under wet conditions and soil properties (ECa and

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 63

**Figure 2.** a) EOF1 and (b) EOF2 computed from the spatial anomalies of soil moisture. The bars in (c) and (d) are EC1 and EC2, respectively. The spatial average and variance of the soil moisture on each sampling date are also shown in

(c) and (d) for comparison. (Reproduced from Perry and Niemann (2007))

bulk density) in dry conditions (Jawson and Niemann, 2007).

**Table 2.** Matrix of Spearman's rank correlation coefficient between series of soil water storage measurements obtained on the 17 locations on seven dates. (Reproduced from Vachaud et al. (1985))

As an alternative to Spearman's rank correlation analysis, Pearson correlation coefficients were also used to describe time stability of the spatial pattern (Cosh et al., 2004). It was ob‐ served that the Pearson correlation analysis revealed a pattern similar to that of the Spear‐ man's rank correlation analysis (Cosh et al., 2004). Note that the tendency of soil moisture sensors to occasionally report erroneous measurements may negatively affect the value of Spearman's rank correlation coefficient, while the Pearson correlation coefficient is less sen‐ sitive to this problem. This is because the total number of measurements will not affect Pear‐ son correlation coefficient as it affects the rank correlation coefficient (Cosh et al., 2006, 2008). However, both Spearman's rank correlation and Pearson correlation analysis are only suitable to describe how much spatial pattern at one measurement time can be preserved at another time (Rolston et al., 1991). If we want to get the information on the similarity of spa‐ tial pattern over multiple measurement occasions, mean value of correlation coefficients can be used (Hu et al., 2009).

With the aim to characterize the degree of time stability of spatial pattern among multiple data series, Starr (2005) suggested two indices, i. e., average spatial coefficient of determina‐ tion (*Rs*<sup>2</sup> ) and temporal coefficient of variability (*CVt* ). The authors reported that 26 to 76% (with mean of 47%) of the observed variability was explained by the time stable pattern. The reported *CVt* value ranged from 11% to 18% (with mean of 16%) for different sampling transects. On the other hand, the *Rs*<sup>2</sup> may not be a good measure of how time stable a *SWC* pattern is at a field that is uniform and completely time stable. This is because the *Rs*<sup>2</sup> will be zero in this case (Starr, 2005).

Perry and Niemann (2007) used the EOF method to characterize the time stable spatial pattern of *SWC*. They identified two dominant spatial patterns of soil water variability at Tarrawarra catchment (Figure 2). The first spatial pattern was more important for medium wet periods and was associated with lateral redistribution of soil water. The second spatial pattern served to modify the first spatial pattern at some days which was associated to evapotranspiration. The importance of each spatial pattern on each measurement occasion can be reflected by the EC values (Figures 2c and 2d). There was a large decrease in EC1 weight during the dry sea‐ son, which agrees well with the finding of Western et al. (1999) that soil moisture patterns on moderately wet occasions show a strong topographic influence. The EOF method can also be used to predict the *SWC* distribution for unobserved times if the average *SWC* and the tempo‐ ral varying coefficient (ECs), which reflect the relative importance of different spatial patterns can be obtained. Average *SWC* can be estimated by many methods such as linear regression (Western et al., 1999), a dynamic multiple linear regression of soil moisture against topograph‐ ic attributes (Wilson et al., 2005), and by measurement from representative location (Grayson and Western, 1998). Examples of EOF method can also be found in other study (Ibrahim and Huggins, 2011; Jawson and Niemann, 2007; Joshi and Mohanty, 2010; Korres et al., 2010; Yoo and Kim, 2004). Different numbers, usually one to four, of significant EOFs were obtained to explain most of *SWC* variations (usually more than 70%). The time-invariant spatial patterns were usually related to topography-related factors (Yoo and Kim, 2004), sand content (Jawson and Niemann, 2007), mixed influences of rainfall, topography, and soil texture (Joshi and Mo‐ hanty, 2010), elevation and wetness index under wet conditions and soil properties (ECa and bulk density) in dry conditions (Jawson and Niemann, 2007).

**Dates 7 Sep 1982**

62 Advances in Agrophysical Research

(1) 1

be used (Hu et al., 2009).

transects. On the other hand, the *Rs*<sup>2</sup>

zero in this case (Starr, 2005).

tion (*Rs*<sup>2</sup>

reported *CVt*

(2) 0.953 1

(3) 0.941 0.953 1

**Average storage (mm)**

28 Oct 1981

(4) 0.988 0.961 0.946 1

(5) 0.953 0.922 0.882 0.968 1

obtained on the 17 locations on seven dates. (Reproduced from Vachaud et al. (1985))

) and temporal coefficient of variability (*CVt*

(6) 0.863 0.789 0.843 0.836 0.882 1

**Table 2.** Matrix of Spearman's rank correlation coefficient between series of soil water storage measurements

(7) 0.860 0.824 0.824 0.794 0.863 0.939 1

As an alternative to Spearman's rank correlation analysis, Pearson correlation coefficients were also used to describe time stability of the spatial pattern (Cosh et al., 2004). It was ob‐ served that the Pearson correlation analysis revealed a pattern similar to that of the Spear‐ man's rank correlation analysis (Cosh et al., 2004). Note that the tendency of soil moisture sensors to occasionally report erroneous measurements may negatively affect the value of Spearman's rank correlation coefficient, while the Pearson correlation coefficient is less sen‐ sitive to this problem. This is because the total number of measurements will not affect Pear‐ son correlation coefficient as it affects the rank correlation coefficient (Cosh et al., 2006, 2008). However, both Spearman's rank correlation and Pearson correlation analysis are only suitable to describe how much spatial pattern at one measurement time can be preserved at another time (Rolston et al., 1991). If we want to get the information on the similarity of spa‐ tial pattern over multiple measurement occasions, mean value of correlation coefficients can

With the aim to characterize the degree of time stability of spatial pattern among multiple data series, Starr (2005) suggested two indices, i. e., average spatial coefficient of determina‐

(with mean of 47%) of the observed variability was explained by the time stable pattern. The

Perry and Niemann (2007) used the EOF method to characterize the time stable spatial pattern of *SWC*. They identified two dominant spatial patterns of soil water variability at Tarrawarra catchment (Figure 2). The first spatial pattern was more important for medium wet periods and was associated with lateral redistribution of soil water. The second spatial pattern served to modify the first spatial pattern at some days which was associated to evapotranspiration. The importance of each spatial pattern on each measurement occasion can be reflected by the

pattern is at a field that is uniform and completely time stable. This is because the *Rs*<sup>2</sup>

value ranged from 11% to 18% (with mean of 16%) for different sampling

**2 Dec 1982**

**29 Jul 1981**

428.1 425.77 424.77 382.06 354.72 347.23 344.46

**18 Aug 1981**

**16 Jul 1982**

). The authors reported that 26 to 76%

will be

may not be a good measure of how time stable a *SWC*

**25 Aug 1981**

**Figure 2.** a) EOF1 and (b) EOF2 computed from the spatial anomalies of soil moisture. The bars in (c) and (d) are EC1 and EC2, respectively. The spatial average and variance of the soil moisture on each sampling date are also shown in (c) and (d) for comparison. (Reproduced from Perry and Niemann (2007))

Time stability of soil water content is usually scale dependent because of the interactions of soil moisture between measurement locations (Kachanoski and De Jong, 1988). The first ef‐ fort to reveal the scale specific time stability was made by Kachanoski and De Jong (1988) who used spectral coherency to explore the time stability of soil water storage in a rolling landscape (Figure 3). They observed that for both wet and drying periods, the spectral co‐ herency was high at large scales (>40 m). For smaller scales (<40 m), however, the spectral coherency was high for the drying period but greatly weakened during the alternating wet‐ ting and drying period. Therefore, their results clearly showed the loss of time stability at small scales (<40 m). Gómez-Plaza et al. (2000) also used the spectral coherency to explore the scale dependence of time stability in burned and unburned areas in a semi-arid environ‐ ment. They clearly observed that soil moisture patterns were maintained in all cases at large scales (transect scale), while soil moisture tended to become time unstable, depending on the range of scale and section of transect considered at medium scales.

**Figure 3.** Coherency spectrum for soil water storage for the drying (June 28-July 12) and recharge (July 12-October 1) period. (Reproduced from Kachanoski and De Jong (1988))

Spectral coherency was verified to be robust to look for the scale dependency of time stabili‐ ty. Its multi-variate extension is also feasible (Si, 2008), although no publications have been found on this point. However, spatial coherency analysis assumes the stationarity in the measured data series over space, hence may not capture localized features of time stability (Biswas and Si, 2011a).

**Figure 4.** The inter-season wavelet coherencies between the soil water storage measured during (a) spring (2 May 2008) and summer (23 August 2008), (b) spring (2 May 2008) and fall (22 October 2008), and (c) summer (23 August 2008) and fall (22 October 2008) for the whole soil profile. Cross sectional view of the transect with landform ele‐ ments is shown at the top. The X-axis indicates distance along the transect; the Y-axis indicates the scale in meter; the solid black line indicates 5% significance level; the color bar indicates strength of correlation, and the direction of ar‐

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 65

row indicates the phase relationship or type of correlation. (Reproduced from Biswas and Si (2011c)).

Time stability of soil water content is usually scale dependent because of the interactions of soil moisture between measurement locations (Kachanoski and De Jong, 1988). The first ef‐ fort to reveal the scale specific time stability was made by Kachanoski and De Jong (1988) who used spectral coherency to explore the time stability of soil water storage in a rolling landscape (Figure 3). They observed that for both wet and drying periods, the spectral co‐ herency was high at large scales (>40 m). For smaller scales (<40 m), however, the spectral coherency was high for the drying period but greatly weakened during the alternating wet‐ ting and drying period. Therefore, their results clearly showed the loss of time stability at small scales (<40 m). Gómez-Plaza et al. (2000) also used the spectral coherency to explore the scale dependence of time stability in burned and unburned areas in a semi-arid environ‐ ment. They clearly observed that soil moisture patterns were maintained in all cases at large scales (transect scale), while soil moisture tended to become time unstable, depending on

**Figure 3.** Coherency spectrum for soil water storage for the drying (June 28-July 12) and recharge (July 12-October 1)

Spectral coherency was verified to be robust to look for the scale dependency of time stabili‐ ty. Its multi-variate extension is also feasible (Si, 2008), although no publications have been found on this point. However, spatial coherency analysis assumes the stationarity in the measured data series over space, hence may not capture localized features of time stability

period. (Reproduced from Kachanoski and De Jong (1988))

(Biswas and Si, 2011a).

64 Advances in Agrophysical Research

the range of scale and section of transect considered at medium scales.

**Figure 4.** The inter-season wavelet coherencies between the soil water storage measured during (a) spring (2 May 2008) and summer (23 August 2008), (b) spring (2 May 2008) and fall (22 October 2008), and (c) summer (23 August 2008) and fall (22 October 2008) for the whole soil profile. Cross sectional view of the transect with landform ele‐ ments is shown at the top. The X-axis indicates distance along the transect; the Y-axis indicates the scale in meter; the solid black line indicates 5% significance level; the color bar indicates strength of correlation, and the direction of ar‐ row indicates the phase relationship or type of correlation. (Reproduced from Biswas and Si (2011c)).

Wavelet coherency can deal with the non-stationary series, it is therefore suitable to identify both scale and location specific time stability (Biswas and Si, 2011a). An excellent application of wavelet coherency in time stability analysis was made by Biswas and Si (2011a). They identified scales and locations of time stability between the spatial patterns of soil water storage in a hummocky landscape in central Saskatchewan, Canada. As Figure 4 shows, strong time stability at all scales and locations existed between summer and fall, when the vegetation was present in field. The spatial pattern in spring was different from that of summer or fall. The large-scale (> 72 m) time stability was present over the whole transect. The medium-scale spatial patterns were time stable in large depressions that had fewer landform elements. The small-scale time stability was mostly random.

Wavelet coherency can also be used to identify the correlation between different variables at different scales and positions (Lakshmi et al., 2004; Si and Zeleke, 2005). One limitation of wavelet coherency analysis is that it is based on the assumption of linearity of systems. Moreover, the application of wavelet coherency on time stability analysis was only confined to two spatial series. Extension of this method to multiple spatial series is expected to under‐ stand the scale and location specific time stability when soil water content from different seasons are considered together.

MEMD is a data-driven method that can deal with both non-stationary and nonlinear series. Recently, Hu et al. (unpublished data) (2012a) applied MEMD to time stability analysis using soil water storage measurements from 20 occasions in a hummocky landscape in central Sas‐ katchewan, Canada. At first, multiple spatial series of soil water storage were separated into different IMFs. The IMFs for the 0-20 cm layer on a selection of measurement occasions is giv‐ en in Figure 5. Each IMF corresponds to a specific scale common to all the spatial series consid‐ ered. Spearman's rank correlation analysis was then conducted for each scale. According to their study, the dominant scale of variations of soil water storage was 104-128 m (IMF5), for 0-20 cm, 20-40 cm, 40-60 cm, and 60-80 cm, and 315-412 m (IMF7) for 80-100 cm, 100-120 cm, and 120-140 cm. Time stability generally increased with scale and it was the strongest at the dominant scale of variation. At different scales, the time stability generally ranked in the order of intra-season>inter-annual>inter-season. Therefore, the study of Hu et al. (2012a) verified that scale specific time stability of soil moisture can be identified using MEMD method.

According to the study reviewed above, time stability analysis is developing from two spa‐ tial series to multiple series, from one scale to multiple scales. With the advance of method‐ ology, a deeper understanding of spatial patterns of soil moisture can be expected. Future work should make full use of various methods, especially those most recently developed, such as EOF, wavelet coherency, and MEMD to gain greater insight into the spatial pattern dynamics at different scales.

#### **4.2. Identify time stable locations for average soil water estimation**

One of the most important applications of time stability is to identify a representative loca‐ tion for average soil water content estimation for a given area. Identification of representa‐ tive locations cover different land uses, different climate zones, and places with various topographical attributes (Gao et al., 2011).

**Figure 5.** Intrinsic mode functions (IMFs) and residues based on multivariate empirical mode decomposition using soil water storage series (top) of the surface layer (0-20 cm). The horizontal axis is the distance between a sampling loca‐ tion and the origin of the transect. The vertical solid bar in Y-axis shows the scale of IMFs. (Reproduced from Hu et al.

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 67

(2012a)).

Wavelet coherency can deal with the non-stationary series, it is therefore suitable to identify both scale and location specific time stability (Biswas and Si, 2011a). An excellent application of wavelet coherency in time stability analysis was made by Biswas and Si (2011a). They identified scales and locations of time stability between the spatial patterns of soil water storage in a hummocky landscape in central Saskatchewan, Canada. As Figure 4 shows, strong time stability at all scales and locations existed between summer and fall, when the vegetation was present in field. The spatial pattern in spring was different from that of summer or fall. The large-scale (> 72 m) time stability was present over the whole transect. The medium-scale spatial patterns were time stable in large depressions that had fewer

Wavelet coherency can also be used to identify the correlation between different variables at different scales and positions (Lakshmi et al., 2004; Si and Zeleke, 2005). One limitation of wavelet coherency analysis is that it is based on the assumption of linearity of systems. Moreover, the application of wavelet coherency on time stability analysis was only confined to two spatial series. Extension of this method to multiple spatial series is expected to under‐ stand the scale and location specific time stability when soil water content from different

MEMD is a data-driven method that can deal with both non-stationary and nonlinear series. Recently, Hu et al. (unpublished data) (2012a) applied MEMD to time stability analysis using soil water storage measurements from 20 occasions in a hummocky landscape in central Sas‐ katchewan, Canada. At first, multiple spatial series of soil water storage were separated into different IMFs. The IMFs for the 0-20 cm layer on a selection of measurement occasions is giv‐ en in Figure 5. Each IMF corresponds to a specific scale common to all the spatial series consid‐ ered. Spearman's rank correlation analysis was then conducted for each scale. According to their study, the dominant scale of variations of soil water storage was 104-128 m (IMF5), for 0-20 cm, 20-40 cm, 40-60 cm, and 60-80 cm, and 315-412 m (IMF7) for 80-100 cm, 100-120 cm, and 120-140 cm. Time stability generally increased with scale and it was the strongest at the dominant scale of variation. At different scales, the time stability generally ranked in the order of intra-season>inter-annual>inter-season. Therefore, the study of Hu et al. (2012a) verified that scale specific time stability of soil moisture can be identified using MEMD method.

According to the study reviewed above, time stability analysis is developing from two spa‐ tial series to multiple series, from one scale to multiple scales. With the advance of method‐ ology, a deeper understanding of spatial patterns of soil moisture can be expected. Future work should make full use of various methods, especially those most recently developed, such as EOF, wavelet coherency, and MEMD to gain greater insight into the spatial pattern

One of the most important applications of time stability is to identify a representative loca‐ tion for average soil water content estimation for a given area. Identification of representa‐ tive locations cover different land uses, different climate zones, and places with various

**4.2. Identify time stable locations for average soil water estimation**

landform elements. The small-scale time stability was mostly random.

seasons are considered together.

66 Advances in Agrophysical Research

dynamics at different scales.

topographical attributes (Gao et al., 2011).

**Figure 5.** Intrinsic mode functions (IMFs) and residues based on multivariate empirical mode decomposition using soil water storage series (top) of the surface layer (0-20 cm). The horizontal axis is the distance between a sampling loca‐ tion and the origin of the transect. The vertical solid bar in Y-axis shows the scale of IMFs. (Reproduced from Hu et al. (2012a)).

Various time stability indices for use at the point scale have been used to identify represen‐ tative locations. Most studies have been aimed at finding locations where soil water content can directly represent the mean of an area (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Martínez-Fernández and Ceballos; 2005). In this way, *SWC* at the time stable location should approximate the average *SWC* of the study area. By extension, point scale time stability representing the areal average also means that the mean relative difference of soil water content, δ<sup>i</sup> , at the stable location *i* should be close to zero at all times. Practically, it is difficult to find a location where the mean relative difference is the close to zero at all times. In this situation, the time stability index *SDRD* is widely used to judge the persistence of the stable location. Therefore, the mean relative difference and associated *SDRD* should be considered together. Successful application of this method can be found in many publica‐ tions (Brocca et al., 2009; Cosh et al., 2006, 2008; Gómez-Plaza et al., 2000; Grayson and West‐ ern, 1998; Tallon and Si, 2004). In the study of Grayson and Western (1998), although the overall spatial pattern of soil moisture was time unstable, some locations were time stable and can represent the average *SWC* of a landscape. In the Tarrawarra watershed (Figure 6), the mean relative differences for locations 1, 4, 12, and 13 were close to zero, and their *SDRD* were also relatively small, so these locations were time stable and could be used as represen‐ tative location for average *SWC* estimation. In the study of Tallon and Si (2004) which was conducted on a chernozemic soil on the Canadian prairies (75 km northeast of Saskatoon, SK, Canada) time stability of spatial patterns of soil moisture for different depths (30, 60, 90, 120 and 160 cm) were observed. With the exception of one location that occurred at two sep‐ arate depths, the locations representative of average *SWC* differed with soil depths.

Although mean relative difference and *SDRD* have been used successfully in the past in identifying representative locations for average *SWC*, the question becomes a matter of which index is more important when the location with mean relative difference closest to zero disagrees with that with the minimum *SDRD*. Most researchers believed that the *SDRD* was more important (Brocca et al., 2009; Cosh et al., 2006, 2008; Gómez-Plaza et al., 2000; Grayson and Western, 1998), while mean relative difference was also emphasized (Martí‐

Instead of considering two indices simultaneously, one index, *RMSE*, combines both indices into one, thus making the identification of time stable location more objective. By this index, drier locations were usually found to have lower *RMSE* values than the wetter points (Cosh et al., 2004, 2006). This suggests that improved consistency will result from the selection of slight‐ ly drier sampling locations (Cosh et al., 2004). Operationally, it facilitates the selection of repre‐ sentative locations. However, it may confuse the concept of time stability, because for locations with large absolute values of relative difference (i. e. extremely dry or wet conditions) but very small *SDRD*, they may have large value of *RMSE* although they are completely time stable.

Another method for estimating average soil water content with time stable location involves a consideration of the constant offset (mean relative difference) for time stable location

( ) ( ) *<sup>j</sup>*

*SWC i*

*SWC*

In this situation, time stable locations are usually defined by the value of *SDRD (*Grayson and Western, 1998). However, there are two limitations in using *SDRD* as a criterion. First, the value of the *SDRD* may be affected to some extent by the magnitude of the mean relative difference; second, the value of the *SDRD* cannot be directly related to the estimation error

Indices that can consider the estimation error are needed for identifying representative time stable locations. One example was conducted by Guber et al. (2008). They applied *χ* 2(*i*) and *D*(*i*) to rank the utility of sampling locations to estimate average *SWC* by considering the estimation error. Another alternative is the mean absolute bias error (*MABE*), which was de‐ veloped directly to relate to the relative error of mean estimation (Hu et al., 2010a, 2010b). Although *MABE* was positively correlated to *SDRD* in some cases (Gao et al., 2011; Hu et al., 2010a, 2010b), the lowest value of *SDRD* did not coincide with the lowest value of *MABE (*Figure 7). Therefore, *MABE* outperformed *SDRD* in terms of representative location identi‐

The degree of representativeness of time stable locations may also change with season or year due to the change of the relative importance of different factors affecting time stability (Guber et al., 2008; Martínez-Fernández and Ceballos, 2003). Schneider et al. (2008) showed that the selected representative locations were appropriate to predict average soil water con‐

*j*

= (33)

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 69

(Grayson and Western, 1998). Then *SWC <sup>j</sup>* can be estimated by

fication for indirect estimation of average *SWC* (Hu et al., 2010b).

*j*

b

*i*

nez-Fernández and Ceballos, 2003).

of the mean value (Hu et al., 2010b).

**Figure 6.** Ranked mean relative differences of soil water content for 0-60 cm measured by neutron probe in the Tarra‐ warra watershed. Also shown is one standard deviation error and site numbers (Reproduced from Grayson and West‐ ern (1998)).

Although mean relative difference and *SDRD* have been used successfully in the past in identifying representative locations for average *SWC*, the question becomes a matter of which index is more important when the location with mean relative difference closest to zero disagrees with that with the minimum *SDRD*. Most researchers believed that the *SDRD* was more important (Brocca et al., 2009; Cosh et al., 2006, 2008; Gómez-Plaza et al., 2000; Grayson and Western, 1998), while mean relative difference was also emphasized (Martí‐ nez-Fernández and Ceballos, 2003).

Various time stability indices for use at the point scale have been used to identify represen‐ tative locations. Most studies have been aimed at finding locations where soil water content can directly represent the mean of an area (Brocca et al., 2009; Cosh et al., 2008; Grayson and Western, 1998; Martínez-Fernández and Ceballos; 2005). In this way, *SWC* at the time stable location should approximate the average *SWC* of the study area. By extension, point scale time stability representing the areal average also means that the mean relative difference of soil water content, δ<sup>i</sup> , at the stable location *i* should be close to zero at all times. Practically, it is difficult to find a location where the mean relative difference is the close to zero at all times. In this situation, the time stability index *SDRD* is widely used to judge the persistence of the stable location. Therefore, the mean relative difference and associated *SDRD* should be considered together. Successful application of this method can be found in many publica‐ tions (Brocca et al., 2009; Cosh et al., 2006, 2008; Gómez-Plaza et al., 2000; Grayson and West‐ ern, 1998; Tallon and Si, 2004). In the study of Grayson and Western (1998), although the overall spatial pattern of soil moisture was time unstable, some locations were time stable and can represent the average *SWC* of a landscape. In the Tarrawarra watershed (Figure 6), the mean relative differences for locations 1, 4, 12, and 13 were close to zero, and their *SDRD* were also relatively small, so these locations were time stable and could be used as represen‐ tative location for average *SWC* estimation. In the study of Tallon and Si (2004) which was conducted on a chernozemic soil on the Canadian prairies (75 km northeast of Saskatoon, SK, Canada) time stability of spatial patterns of soil moisture for different depths (30, 60, 90, 120 and 160 cm) were observed. With the exception of one location that occurred at two sep‐

arate depths, the locations representative of average *SWC* differed with soil depths.

**Figure 6.** Ranked mean relative differences of soil water content for 0-60 cm measured by neutron probe in the Tarra‐ warra watershed. Also shown is one standard deviation error and site numbers (Reproduced from Grayson and West‐

ern (1998)).

68 Advances in Agrophysical Research

Instead of considering two indices simultaneously, one index, *RMSE*, combines both indices into one, thus making the identification of time stable location more objective. By this index, drier locations were usually found to have lower *RMSE* values than the wetter points (Cosh et al., 2004, 2006). This suggests that improved consistency will result from the selection of slight‐ ly drier sampling locations (Cosh et al., 2004). Operationally, it facilitates the selection of repre‐ sentative locations. However, it may confuse the concept of time stability, because for locations with large absolute values of relative difference (i. e. extremely dry or wet conditions) but very small *SDRD*, they may have large value of *RMSE* although they are completely time stable.

Another method for estimating average soil water content with time stable location involves a consideration of the constant offset (mean relative difference) for time stable location (Grayson and Western, 1998). Then *SWC <sup>j</sup>* can be estimated by

$$\mathcal{B}\_{j}\left(i\right) = \frac{\text{SWC}\_{j}\left(i\right)}{\left\langle \text{SWC}\_{j} \right\rangle} \tag{33}$$

In this situation, time stable locations are usually defined by the value of *SDRD (*Grayson and Western, 1998). However, there are two limitations in using *SDRD* as a criterion. First, the value of the *SDRD* may be affected to some extent by the magnitude of the mean relative difference; second, the value of the *SDRD* cannot be directly related to the estimation error of the mean value (Hu et al., 2010b).

Indices that can consider the estimation error are needed for identifying representative time stable locations. One example was conducted by Guber et al. (2008). They applied *χ* 2(*i*) and *D*(*i*) to rank the utility of sampling locations to estimate average *SWC* by considering the estimation error. Another alternative is the mean absolute bias error (*MABE*), which was de‐ veloped directly to relate to the relative error of mean estimation (Hu et al., 2010a, 2010b). Although *MABE* was positively correlated to *SDRD* in some cases (Gao et al., 2011; Hu et al., 2010a, 2010b), the lowest value of *SDRD* did not coincide with the lowest value of *MABE (*Figure 7). Therefore, *MABE* outperformed *SDRD* in terms of representative location identi‐ fication for indirect estimation of average *SWC* (Hu et al., 2010b).

The degree of representativeness of time stable locations may also change with season or year due to the change of the relative importance of different factors affecting time stability (Guber et al., 2008; Martínez-Fernández and Ceballos, 2003). Schneider et al. (2008) showed that the selected representative locations were appropriate to predict average soil water con‐ tent of the sites for multiple years, although the time stability characteristics of some loca‐ tions varied between years. An interesting study conducted by Guber et al. (2008) indicated that the utility of a location to estimate the estimation accuracy of field average *SWC* was not influenced by measurement frequency.

and soils that could be used to determine the representative locations *a priori (*Martínez-Fer‐

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 71

Time stability of soil water content is controlled by factors that consistently influence soil moisture distribution in a similar way at all times. Most studies focus on factors influencing the underlying spatial pattern of soil moisture, and relatively few studies are associated with factors influencing the degree of time stability at point scales. We will review the con‐

Controlling factors of underlying spatial pattern are usually defined by comparing spatial patterns of time stability and various factors (Hu et al., 2010a; Perry and Niemann, 2007; Schneideret al., 2008). Related studies referred to different climate conditions, from semiarid (Hu et al., 2010a) to humid environments (Perry and Niemann, 2007) and different land uses, such as rangeland (Mohanty and Skaggs, 2001), grass (Jacobs et al., 2004), forest (Lin,

Soil properties such as soil texture are usually found to influence the spatial pattern of time stability. For example, Vachaud et al. (1985) found that the scaled factor of soil water content were significantly correlated to the silt plus clay contents on a field in St. Martin d'Heres, France. Mohanty and Skaggs (2001) observed that the sandy loam soil produced the best time-stable soil moisture pattern and silty loam soil produced an intermediate level of time stability in the Little Washita watershed of the southern part of the Great Plains of the US. Jacobs et al. (2004) pointed out that lower sand content may produce more time stable pat‐ terns in the Iowa. Grant et al. (2004) found a significant correlation between clay content, coarse fragment and the rank of mean relative difference of soil water storage in the Rey‐ nolds Mountain East catchment in the Owyhee Mountains. Starr (2005) found that coarser particle size classes were generally drier than adjacent finer particle size class soil on a pota‐ to farm located near Houlton, Maine. Hu et al. (2010a) recognized that sand and silt content can explain about 41.4% - 65.5 % of the spatial pattern of mean relative difference of soil wa‐ ter content. Except for soil texture, other soil properties such as organic matter content (Hu et al., 2009), bulk density (Jacobs et al., 2004), and soil thickness (Zhu and Lin, 2011) have also been recognized as the controlling factors for time stability of spatial pattern. However, soil properties may also have no influences on the spatial pattern of time stability. For exam‐ ple, Tallon and Si (2004) reported that clay content showed the least amount of control of spatial patterns on a chernozemic soil on the Canadian prairies (75 km northeast of Saska‐ toon, SK, Canada). On ungrazed sites, Schneider et al. (2008) found that neither bulk densi‐ ty, organic carbon nor sand and clay content could explain the time-stable characteristics of the sampling points. Comegna and Basile (1994) and Starr (2005) both observed the absence

nándezandCeballos, 2005).

**5. Controlling factors**

trolling factors of the two types of time stability below.

2006), and agricultural land (Starr et al., 2005).

**5.1. Controlling factors of time stability of spatial pattern**

**Figure 7.** Rank ordered mean absolute bias error (*MABE*) and standard deviation of relative differences (*SDRD*) for the soil water storage of 0-1.0 m measured with neutron probe from October, 2004 to September, 2006. Also shown is location number (Reproduced from Hu et al. (2010b))

Many indices have been available to identify the representative location for average *SWC* estimation. However, few works have been conducted to compare their performances with a few exceptions (Gao et al., 2011; Hu et al., 2010a, 2010b, 2012b). *MABE* was reported to be a better index than *SDRD* in terms of average soil water storage estimation by the indirect method (Hu et al., 2010a, 2010b). More recently, with the soil water storage datasets from a transect in a Canadian prairie area and a watershed on the Chinese Loess Plateau, Hu et al. (2012b) evaluated the performance of different time stability indices in terms of the average soil water storage estimation quality judged by different goodness-of-fit indices. Their re‐ sults showed that *MABE*, *χ<sup>2</sup>* , *D*, and *CVt* outperformed *SDRD* and *T*. If root mean squared deviations (*RMSD*) and Nash-Sutcliffe coefficient of efficiency (*NSCE*) were used as good‐ ness-of-fit indices, *D* was the best time stability index. If absolute mean difference (*BIAS*) and absolute bias relative to mean (*RBIAS*) were adopted, *MABE* was the best time stability index. Their results also showed that average soil water storage estimation by the indirect method was more accurate than that by the direct method. It is necessary to compare their performances under other different hydrological backgrounds to obtain a more suitable in‐ dex for identifying representative location. Furthermore, all these methods still need an in‐ tensive measurement of soil water content over time before a representative location can be identified. According to Martínez-Fernández and Ceballos (2005) about one year (a com‐ plete seasonal cycle) of measurements is needed to determine the representative location. Therefore, there is an urgent need to look for the possible definable features of the terrain and soils that could be used to determine the representative locations *a priori (*Martínez-Fer‐ nándezandCeballos, 2005).

## **5. Controlling factors**

tent of the sites for multiple years, although the time stability characteristics of some loca‐ tions varied between years. An interesting study conducted by Guber et al. (2008) indicated that the utility of a location to estimate the estimation accuracy of field average *SWC* was

**Figure 7.** Rank ordered mean absolute bias error (*MABE*) and standard deviation of relative differences (*SDRD*) for the soil water storage of 0-1.0 m measured with neutron probe from October, 2004 to September, 2006. Also shown is

Many indices have been available to identify the representative location for average *SWC* estimation. However, few works have been conducted to compare their performances with a few exceptions (Gao et al., 2011; Hu et al., 2010a, 2010b, 2012b). *MABE* was reported to be a better index than *SDRD* in terms of average soil water storage estimation by the indirect method (Hu et al., 2010a, 2010b). More recently, with the soil water storage datasets from a transect in a Canadian prairie area and a watershed on the Chinese Loess Plateau, Hu et al. (2012b) evaluated the performance of different time stability indices in terms of the average soil water storage estimation quality judged by different goodness-of-fit indices. Their re‐

deviations (*RMSD*) and Nash-Sutcliffe coefficient of efficiency (*NSCE*) were used as good‐ ness-of-fit indices, *D* was the best time stability index. If absolute mean difference (*BIAS*) and absolute bias relative to mean (*RBIAS*) were adopted, *MABE* was the best time stability index. Their results also showed that average soil water storage estimation by the indirect method was more accurate than that by the direct method. It is necessary to compare their performances under other different hydrological backgrounds to obtain a more suitable in‐ dex for identifying representative location. Furthermore, all these methods still need an in‐ tensive measurement of soil water content over time before a representative location can be identified. According to Martínez-Fernández and Ceballos (2005) about one year (a com‐ plete seasonal cycle) of measurements is needed to determine the representative location. Therefore, there is an urgent need to look for the possible definable features of the terrain

outperformed *SDRD* and *T*. If root mean squared

, *D*, and *CVt*

not influenced by measurement frequency.

70 Advances in Agrophysical Research

location number (Reproduced from Hu et al. (2010b))

sults showed that *MABE*, *χ<sup>2</sup>*

Time stability of soil water content is controlled by factors that consistently influence soil moisture distribution in a similar way at all times. Most studies focus on factors influencing the underlying spatial pattern of soil moisture, and relatively few studies are associated with factors influencing the degree of time stability at point scales. We will review the con‐ trolling factors of the two types of time stability below.

#### **5.1. Controlling factors of time stability of spatial pattern**

Controlling factors of underlying spatial pattern are usually defined by comparing spatial patterns of time stability and various factors (Hu et al., 2010a; Perry and Niemann, 2007; Schneideret al., 2008). Related studies referred to different climate conditions, from semiarid (Hu et al., 2010a) to humid environments (Perry and Niemann, 2007) and different land uses, such as rangeland (Mohanty and Skaggs, 2001), grass (Jacobs et al., 2004), forest (Lin, 2006), and agricultural land (Starr et al., 2005).

Soil properties such as soil texture are usually found to influence the spatial pattern of time stability. For example, Vachaud et al. (1985) found that the scaled factor of soil water content were significantly correlated to the silt plus clay contents on a field in St. Martin d'Heres, France. Mohanty and Skaggs (2001) observed that the sandy loam soil produced the best time-stable soil moisture pattern and silty loam soil produced an intermediate level of time stability in the Little Washita watershed of the southern part of the Great Plains of the US. Jacobs et al. (2004) pointed out that lower sand content may produce more time stable pat‐ terns in the Iowa. Grant et al. (2004) found a significant correlation between clay content, coarse fragment and the rank of mean relative difference of soil water storage in the Rey‐ nolds Mountain East catchment in the Owyhee Mountains. Starr (2005) found that coarser particle size classes were generally drier than adjacent finer particle size class soil on a pota‐ to farm located near Houlton, Maine. Hu et al. (2010a) recognized that sand and silt content can explain about 41.4% - 65.5 % of the spatial pattern of mean relative difference of soil wa‐ ter content. Except for soil texture, other soil properties such as organic matter content (Hu et al., 2009), bulk density (Jacobs et al., 2004), and soil thickness (Zhu and Lin, 2011) have also been recognized as the controlling factors for time stability of spatial pattern. However, soil properties may also have no influences on the spatial pattern of time stability. For exam‐ ple, Tallon and Si (2004) reported that clay content showed the least amount of control of spatial patterns on a chernozemic soil on the Canadian prairies (75 km northeast of Saska‐ toon, SK, Canada). On ungrazed sites, Schneider et al. (2008) found that neither bulk densi‐ ty, organic carbon nor sand and clay content could explain the time-stable characteristics of the sampling points. Comegna and Basile (1994) and Starr (2005) both observed the absence of the controlling influence of soil texture and they attribute it to the homogeneous distribu‐ tion of soil texture in their study areas.

will dominate the spatial distribution of soil moisture under very wet conditions (Western et al., 1999). Non-local control can also be found in the study of Biswas and Si (2012). In their study, the coefficient of determination between soil water storage and elevation increased more than eightfold after shifting the spatial series of soil water storage by a length equal to that of the existing slope. Although local and non-local control was widely recognized, it is still difficult to make a very clear limit between them, most likely due to the possible influ‐

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 73

Because of the different controls in different soil water conditions, time stability of soil water content tends to be stronger under the similar soil water conditions than those with quite different soil water conditions. For example, Martinez-Fernández and Ceballos (2003) ob‐ served that the periods with the lowest time stability coincided with situations involving the transition from dry to wet. Gao et al. (2011) revealed that the time stability of root zone soil water (0–60 cm) was higher in either dry or wet season than that including both, and soil

In addition to soil and topography, other factors such as vegetation (Pachepsky et al., 2005; Schneideret al., 2008), soil depth (Huet al., 2009) can also influence the time stability of soil moisture. For example, Gómez- Plaza et al. (2000) found that the existence of vegetation could destroy the time stability of spatial patterns of soil moisture. Mohanty and Skaggs (2001) found varying degrees of time stability depending on vegetation cover as well as top‐ ography. Generally, time stability of spatial pattern is strongest in deep soil layers than in shallow layers (Cassel et al., 2000). However, stronger time stability of spaital pattern at shallower depth (20 cm) than deeper soil depths was also observed (Hu et al., 2009, 2010a). Although many factors have been revealed to influence the time stability of spatial pattern, no single, or even combined, factor can explain all time stable patterns. This suggests that some influencing factors may still remain undiscovered. In addition, random measurement errors also contribute to some variability of the time stable pattern (Starr, 2005). According to the study of Starr (2005), random measurement error can contribute 20% of the variability of underlying spatial pattern. Therefore, to better understand the spatio-temporal variability and time stability of spatial pattern of soil moisture, controlling factors should be made clearer when more accurate measurement of soil water content should also be expected.

Revealing the controlling factors of time stability at the point scale is crucial in identifying time stable locations for average *SWC* estimation. Soil and topography are usually found to influence the time stability of soil water content at point scales. Brocca et al. (2009) observed that the time stable positions were linked to topographic characteristics, primarily the up‐ slope drainage area but also elevation and slope. Cosh et al. (2008) found that almost 60% of the *SDRD* variability can be credited to the soil type parameters. Soil parameters as major controlling factors is a finding that is in agreement with those of Hu et al. (2010a) who found that soil texture can significantly (P<0.05) affect the stability of soil water content. Gao et al. (2011) observed that time-stable locations corresponded to relatively high clay contents, mild slopes and planar surfaces, which agreed with the previous studies, where both soil

water exhibited very low time stability during the transition period from dry to wet.

**5.2. Controlling factors of time stability at the point scale**

ences of topography on soil properties

Topographic characteristics have also been recognized to influence the time stability of spa‐ tial patterns. Gómez-Plaza et al. (2000) showed that when the factors affecting soil moisture were limited to topographical position or local topography at the transect scale, spatial pat‐ terns of soil moisture presented time stability. Joshi et al. (2011) observed that fields with flat topography (LW21) showed the worst time stability features compared to the fields having gently rolling topography (LW03 and LW13) in Oklahoma. The influence of topography on the time stability of spatial pattern may be scale dependent. For example, although the Pear‐ son correlation coefficient between soil water storage and elevation was very small, Biswas and Si (2011a) identified that strong scale-specific correlations existed between soil water storage and elevation at different scales using wavelet coherency, which contributed to the time stability of soil water storage at different scales. Kachanoski and De Jong (1988) attrib‐ uted to the loss of time stability of soil water storage at scale smaller than 40 m to the role of spatial pattern of surface curvature.

Usually, soil properties and topographical properties jointly control the time stability of soil moisture spatial pattern. According to Joshi et al. (2011), soil properties (i. e., percent silt, percent sand, and soil texture) and topography (elevation and slope) jointly affected the spa‐ tiotemporal evolution and time stability of soil moisture at both point and footprint scales in the Little Washita watershed, Oklahoma and in the Walnut Creek watershed, Iowa. Zhu and Lin (2011) also observed that both soil and terrain controlled soil moisture variation at dif‐ ferent seasons and soil depths at farm scale.

The relative role of soil and topographic properties in controlling the spatial variability and time stable pattern of soil moisture is usually related to the dominant hydrological process taking place for different soil water conditions between wet and dry periods. The *SWC* pat‐ terns existing between wet and dry periods are referred to as preferred states (Graysonet al., 1997; Western et al., 1999). When dry conditions dominate, evapotranspiration is greater than precipitation, and local controls dominate mainly by vertical movement of soil water with no connection between adjacent points. In this situation, differences in vegetation and soil properties such as texture would be responsible for the spatial pattern of soil moisture (Graysonet al., 1997). During wet states, evapotranspiration is less than precipitation and non-local controls govern soil moisture distribution. These controls are related to upslope topography and include catchment area, aspect, depth and soil profile curvature (Gómez-Plaza et al., 2001). Western at al. (1999) attributed these differences in local and non-local controls to the high degree of organization during wet periods that consist of connected bands of high soil water content in drainage lines. Soil water in these drainage lines is later‐ ally redistributed by both overland and subsurface flow. Using EOF analysis, Perry and Nie‐ mann (2007) found that underlying spatial patterns of soil moisture were controlled by local soil properties in wet and dry conditions and topographic characteristics during intermedi‐ ate conditions. Their results agreed with the findings of Western et al. (1999) that soil mois‐ ture patterns on moderately wet dates show a strong topographic influence, one that seems to diminish on both dry and very wet days. Note that soil porosity, rather than topography, will dominate the spatial distribution of soil moisture under very wet conditions (Western et al., 1999). Non-local control can also be found in the study of Biswas and Si (2012). In their study, the coefficient of determination between soil water storage and elevation increased more than eightfold after shifting the spatial series of soil water storage by a length equal to that of the existing slope. Although local and non-local control was widely recognized, it is still difficult to make a very clear limit between them, most likely due to the possible influ‐ ences of topography on soil properties

of the controlling influence of soil texture and they attribute it to the homogeneous distribu‐

Topographic characteristics have also been recognized to influence the time stability of spa‐ tial patterns. Gómez-Plaza et al. (2000) showed that when the factors affecting soil moisture were limited to topographical position or local topography at the transect scale, spatial pat‐ terns of soil moisture presented time stability. Joshi et al. (2011) observed that fields with flat topography (LW21) showed the worst time stability features compared to the fields having gently rolling topography (LW03 and LW13) in Oklahoma. The influence of topography on the time stability of spatial pattern may be scale dependent. For example, although the Pear‐ son correlation coefficient between soil water storage and elevation was very small, Biswas and Si (2011a) identified that strong scale-specific correlations existed between soil water storage and elevation at different scales using wavelet coherency, which contributed to the time stability of soil water storage at different scales. Kachanoski and De Jong (1988) attrib‐ uted to the loss of time stability of soil water storage at scale smaller than 40 m to the role of

Usually, soil properties and topographical properties jointly control the time stability of soil moisture spatial pattern. According to Joshi et al. (2011), soil properties (i. e., percent silt, percent sand, and soil texture) and topography (elevation and slope) jointly affected the spa‐ tiotemporal evolution and time stability of soil moisture at both point and footprint scales in the Little Washita watershed, Oklahoma and in the Walnut Creek watershed, Iowa. Zhu and Lin (2011) also observed that both soil and terrain controlled soil moisture variation at dif‐

The relative role of soil and topographic properties in controlling the spatial variability and time stable pattern of soil moisture is usually related to the dominant hydrological process taking place for different soil water conditions between wet and dry periods. The *SWC* pat‐ terns existing between wet and dry periods are referred to as preferred states (Graysonet al., 1997; Western et al., 1999). When dry conditions dominate, evapotranspiration is greater than precipitation, and local controls dominate mainly by vertical movement of soil water with no connection between adjacent points. In this situation, differences in vegetation and soil properties such as texture would be responsible for the spatial pattern of soil moisture (Graysonet al., 1997). During wet states, evapotranspiration is less than precipitation and non-local controls govern soil moisture distribution. These controls are related to upslope topography and include catchment area, aspect, depth and soil profile curvature (Gómez-Plaza et al., 2001). Western at al. (1999) attributed these differences in local and non-local controls to the high degree of organization during wet periods that consist of connected bands of high soil water content in drainage lines. Soil water in these drainage lines is later‐ ally redistributed by both overland and subsurface flow. Using EOF analysis, Perry and Nie‐ mann (2007) found that underlying spatial patterns of soil moisture were controlled by local soil properties in wet and dry conditions and topographic characteristics during intermedi‐ ate conditions. Their results agreed with the findings of Western et al. (1999) that soil mois‐ ture patterns on moderately wet dates show a strong topographic influence, one that seems to diminish on both dry and very wet days. Note that soil porosity, rather than topography,

tion of soil texture in their study areas.

72 Advances in Agrophysical Research

spatial pattern of surface curvature.

ferent seasons and soil depths at farm scale.

Because of the different controls in different soil water conditions, time stability of soil water content tends to be stronger under the similar soil water conditions than those with quite different soil water conditions. For example, Martinez-Fernández and Ceballos (2003) ob‐ served that the periods with the lowest time stability coincided with situations involving the transition from dry to wet. Gao et al. (2011) revealed that the time stability of root zone soil water (0–60 cm) was higher in either dry or wet season than that including both, and soil water exhibited very low time stability during the transition period from dry to wet.

In addition to soil and topography, other factors such as vegetation (Pachepsky et al., 2005; Schneideret al., 2008), soil depth (Huet al., 2009) can also influence the time stability of soil moisture. For example, Gómez- Plaza et al. (2000) found that the existence of vegetation could destroy the time stability of spatial patterns of soil moisture. Mohanty and Skaggs (2001) found varying degrees of time stability depending on vegetation cover as well as top‐ ography. Generally, time stability of spatial pattern is strongest in deep soil layers than in shallow layers (Cassel et al., 2000). However, stronger time stability of spaital pattern at shallower depth (20 cm) than deeper soil depths was also observed (Hu et al., 2009, 2010a).

Although many factors have been revealed to influence the time stability of spatial pattern, no single, or even combined, factor can explain all time stable patterns. This suggests that some influencing factors may still remain undiscovered. In addition, random measurement errors also contribute to some variability of the time stable pattern (Starr, 2005). According to the study of Starr (2005), random measurement error can contribute 20% of the variability of underlying spatial pattern. Therefore, to better understand the spatio-temporal variability and time stability of spatial pattern of soil moisture, controlling factors should be made clearer when more accurate measurement of soil water content should also be expected.

#### **5.2. Controlling factors of time stability at the point scale**

Revealing the controlling factors of time stability at the point scale is crucial in identifying time stable locations for average *SWC* estimation. Soil and topography are usually found to influence the time stability of soil water content at point scales. Brocca et al. (2009) observed that the time stable positions were linked to topographic characteristics, primarily the up‐ slope drainage area but also elevation and slope. Cosh et al. (2008) found that almost 60% of the *SDRD* variability can be credited to the soil type parameters. Soil parameters as major controlling factors is a finding that is in agreement with those of Hu et al. (2010a) who found that soil texture can significantly (P<0.05) affect the stability of soil water content. Gao et al. (2011) observed that time-stable locations corresponded to relatively high clay contents, mild slopes and planar surfaces, which agreed with the previous studies, where both soil type and topography influenced time stability at point scales (Jacobs et al., 2004; Joshi et al., 2011; Mohanty and Skaggs, 2011). However, poor relationships between time stable loca‐ tions and soil and topographic properties were also found (Tallon and Si, 2004). Land use was observed to have no influences on time stability of soil water content at point scales (Hu et al., 2010a; Joshi et al., 2011).

all similarity of soil water content spatial patterns between different occasions, and the other is to describe the time invariance of the relative soil water content with time at point scales; (2) Time stability can be used to analyze the time persistence of spatial pattern of soil water content. Spearman's rank correlation, Pearson correlation analysis, average spatial coeffi‐ cient of determination, and temporal coefficient of variability can be used to characterize the time stability of spatial pattern for two spatial series. Empirical orthogonal function can be used to extract the same spatial patterns of soil moisture for multiple spatial series. Spectral coherency, wavelet coherency, and multivariate empirical mode decomposition can be used to identify the scale specific time stability of a spatial pattern. However, the spectral coher‐ ency method assumes spatial series are stationary and linear, wavelet methods assume the spatial series to be linear, and both of them are only applicable to the case of two spatial ser‐ ies. For multivariate empirical mode decomposition, it can deal with non-stationary and nonlinear systems with multiple spatial series; (3) Time stability can be used to identify the most time-stable locations for average *SWC* evaluation for a given field. Standard deviation of relative difference, root mean square error, standard deviation of soil water content, width of the 90% empirical tolerance interval of relative water content, and mean absolute bias error can be used to identify the most time-stable location. However, the performance of these indices should be compared in different environmental backgrounds to look for the most suitable index; (4) Time stability of soil water content can be influenced by many fac‐ tors, including soil, topography, vegetation, and climate. It also shows non-local control and local control, depending on the soil water conditions. Knowledge of the controlling factors for time stability of spatial pattern and that at point scale is important for understanding the soil water processes, soil water management, and identification of time-stable locations for average *SWC* evaluation. Future work should quantify the influences of different factors on time stability, and make deeper understanding of nonlocal and local control on time stabili‐ ty of soil water content; (5) Time stability of soil water content is also depth-wise and scale specific. Usually, time stability of soil moisture increased with depth and scale. Time stabili‐ ty analysis can not only facilitate our understanding to soil water related processes but also

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 75

greatly reduce the sampling work for soil water content in fields.

, Asim Biswas2

1 University of Saskatchewan, Department of Soil Science, Saskatoon, Canada

The project was funded by Natural Science and Engineering council (NSERC) of Canada.

and Bing Cheng Si1

**Acknowledgements**

**Author details**

, Lindsay K. Tallon1

\*Address all correspondence to: bing.si@usask.ca

Wei Hu1

Time stability at point scales was usually associated with the soil water conditions of that location. For example, Martinez-Fernández and Ceballos (2005) found that time stability as characterized by *SDRD* was usually stronger in dry locations than wet locations for all depths, and they attribute it to the predominance of the sandy fraction with the resultant weak ability to retain water in dry locations. Similar results were obtained in other studies (Bosch et al., 2006; Jacobs et al., 2004; Martínez-Fernández and Ceballos, 2003). However, the relationship of time stability and soil water condition is likely influenced by the time stabili‐ ty index used. In the LYMQ watershed on the Loess Plateau in China, Hu et al. (2010a) also found stronger time stability in drier locations using *SDRD*, while they observed stronger time stability in wetter locations when using *MABE*.

Time stability at point scales depends on different periods and soil depths. Generally, stron‐ ger time stability was found in wetter periods (Guber et al., 2008; Zhao et al., 2010; Zhou et al., 2007). Zhou et al. (2007) believed that the lateral redistribution of soil moisture in the wet/dry transition period and uneven evapotranspiration and root water uptake in drying process was the reasons for the relatively higher time stability during wetting periods. Gub‐ er et al. (2008) related the weaker time stability in the May–June period to the active vegeta‐ tive growth of corn. Zhao et al. (2010) attributed the high temporal stability under wet conditions to an enhanced capillary movement of water from the subsoil to the topsoil. Time stability at point scale was usually found to be stronger at deeper depths observed by differ‐ ent time stability indices (Guberet al., 2008; Hu et al., 2010a, 2010b; Starks et al., 2006). This was generally in agreement with the relationship of time stability of spatial pattern with depth (Cassel et al., 2000).

Due to the complex influences of various factors on the time stability at the point scale, fu‐ ture work should be devoted to quantitatively define the contribution of various factors on time stability of soil water content. Models of the relationship of time stability and environ‐ mental factors should also be developed for the purpose of identification of time stable posi‐ tions with known soil and topographic information. In addition, because of the different relationships between factors and different time stability indices (Hu et al., 2010a), attention should also be paid to the selection of time stability index for average *SWC* estimation in a field.

## **6. Conclusion**

We reviewed the studies on time stability of soil water content, including its concept, meth‐ odology, application, and controlling factors. The following conclusions can be drawn: (1) Two types of time stability of soil water content can be classified; one is to describe the over‐ all similarity of soil water content spatial patterns between different occasions, and the other is to describe the time invariance of the relative soil water content with time at point scales; (2) Time stability can be used to analyze the time persistence of spatial pattern of soil water content. Spearman's rank correlation, Pearson correlation analysis, average spatial coeffi‐ cient of determination, and temporal coefficient of variability can be used to characterize the time stability of spatial pattern for two spatial series. Empirical orthogonal function can be used to extract the same spatial patterns of soil moisture for multiple spatial series. Spectral coherency, wavelet coherency, and multivariate empirical mode decomposition can be used to identify the scale specific time stability of a spatial pattern. However, the spectral coher‐ ency method assumes spatial series are stationary and linear, wavelet methods assume the spatial series to be linear, and both of them are only applicable to the case of two spatial ser‐ ies. For multivariate empirical mode decomposition, it can deal with non-stationary and nonlinear systems with multiple spatial series; (3) Time stability can be used to identify the most time-stable locations for average *SWC* evaluation for a given field. Standard deviation of relative difference, root mean square error, standard deviation of soil water content, width of the 90% empirical tolerance interval of relative water content, and mean absolute bias error can be used to identify the most time-stable location. However, the performance of these indices should be compared in different environmental backgrounds to look for the most suitable index; (4) Time stability of soil water content can be influenced by many fac‐ tors, including soil, topography, vegetation, and climate. It also shows non-local control and local control, depending on the soil water conditions. Knowledge of the controlling factors for time stability of spatial pattern and that at point scale is important for understanding the soil water processes, soil water management, and identification of time-stable locations for average *SWC* evaluation. Future work should quantify the influences of different factors on time stability, and make deeper understanding of nonlocal and local control on time stabili‐ ty of soil water content; (5) Time stability of soil water content is also depth-wise and scale specific. Usually, time stability of soil moisture increased with depth and scale. Time stabili‐ ty analysis can not only facilitate our understanding to soil water related processes but also greatly reduce the sampling work for soil water content in fields.

## **Acknowledgements**

type and topography influenced time stability at point scales (Jacobs et al., 2004; Joshi et al., 2011; Mohanty and Skaggs, 2011). However, poor relationships between time stable loca‐ tions and soil and topographic properties were also found (Tallon and Si, 2004). Land use was observed to have no influences on time stability of soil water content at point scales (Hu

Time stability at point scales was usually associated with the soil water conditions of that location. For example, Martinez-Fernández and Ceballos (2005) found that time stability as characterized by *SDRD* was usually stronger in dry locations than wet locations for all depths, and they attribute it to the predominance of the sandy fraction with the resultant weak ability to retain water in dry locations. Similar results were obtained in other studies (Bosch et al., 2006; Jacobs et al., 2004; Martínez-Fernández and Ceballos, 2003). However, the relationship of time stability and soil water condition is likely influenced by the time stabili‐ ty index used. In the LYMQ watershed on the Loess Plateau in China, Hu et al. (2010a) also found stronger time stability in drier locations using *SDRD*, while they observed stronger

Time stability at point scales depends on different periods and soil depths. Generally, stron‐ ger time stability was found in wetter periods (Guber et al., 2008; Zhao et al., 2010; Zhou et al., 2007). Zhou et al. (2007) believed that the lateral redistribution of soil moisture in the wet/dry transition period and uneven evapotranspiration and root water uptake in drying process was the reasons for the relatively higher time stability during wetting periods. Gub‐ er et al. (2008) related the weaker time stability in the May–June period to the active vegeta‐ tive growth of corn. Zhao et al. (2010) attributed the high temporal stability under wet conditions to an enhanced capillary movement of water from the subsoil to the topsoil. Time stability at point scale was usually found to be stronger at deeper depths observed by differ‐ ent time stability indices (Guberet al., 2008; Hu et al., 2010a, 2010b; Starks et al., 2006). This was generally in agreement with the relationship of time stability of spatial pattern with

Due to the complex influences of various factors on the time stability at the point scale, fu‐ ture work should be devoted to quantitatively define the contribution of various factors on time stability of soil water content. Models of the relationship of time stability and environ‐ mental factors should also be developed for the purpose of identification of time stable posi‐ tions with known soil and topographic information. In addition, because of the different relationships between factors and different time stability indices (Hu et al., 2010a), attention should also be paid to the selection of time stability index for average *SWC* estimation in a

We reviewed the studies on time stability of soil water content, including its concept, meth‐ odology, application, and controlling factors. The following conclusions can be drawn: (1) Two types of time stability of soil water content can be classified; one is to describe the over‐

et al., 2010a; Joshi et al., 2011).

74 Advances in Agrophysical Research

depth (Cassel et al., 2000).

field.

**6. Conclusion**

time stability in wetter locations when using *MABE*.

The project was funded by Natural Science and Engineering council (NSERC) of Canada.

## **Author details**

Wei Hu1 , Lindsay K. Tallon1 , Asim Biswas2 and Bing Cheng Si1

\*Address all correspondence to: bing.si@usask.ca

1 University of Saskatchewan, Department of Soil Science, Saskatoon, Canada

2 Department of Natural Resource Sciences, McGill University, Ste-Anne-de-Bellevue, Que‐ bec, Canada

[14] Gao, X. D., Wu, P. T., Zhao, X. N., Shi, Y. G., & Wang, J. W. (2011). Estimating Spatial Mean Soil Water Contents of Sloping Jujube Orchards Using Temporal Stability. Agr.

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 77

[15] Gómez- Plaza, A., Alvarez-Rogel, J., Albaladejo, J., & Castillo, V. M. (2000). Spatial Patterns and Temporal Stability of Soil Moisture across a Range of Scales in a Semi-

[16] Gómez-Plaza, A., Martínez-Mena, M., Albaladejo, J., & Castillo, V. M. (2001). Factors Regulating Spatial Distribution of Soil Water Content in Small Semi-Arid Catch‐

[17] Grant, L., Seyfried, M., & McNamara, J. (2004). Spatial Variation and Temporal Sta‐ bility of Soil Water in a Snow-Dominated, Mountain Catchment. Hydrol. process., 18,

[18] Grayson, R. B., & Western, A. W. (1998). Towards Areal Estimation of Soil Water Content from Point Measurements: Time and Space Stability of Mean Response. J.

[19] Grayson, R. B., Western, A. W., Chiew, F. H. S., & Blöschl, G. (1997). Preferred States in Spatial Soil Moisture Patterns: Local and Nonlocal Controls. Water resour. res., 33,

[20] Grinsted, A. J., Moore, C., & Jevrejeva, S. (2004). Application of the Cross Wavelet Transform and Wavelet Coherence to Geophysical Time Series. Nonlinear proc. ge‐

[21] Guber, A. K., Gish, T. J., Pachepsky, Y. A., van Genuchten, M. T., Daughtry, C. S. T., Nicholson, T. J., & Cady, R. E. (2008). Temporal Stability in Soil Water Content Pat‐

[22] Hu, W., Biswas, A., & Si, B. C. (2012a). Application of Multivariate Empirical Mode Decomposition for Revealing Scale Specific Time Stability of Soil Water Storage. J.

[23] Hu, W., Shao, M. A., Han, F. P., Reichardt, K., & Tan, J. (2010a). Watershed Scale

[24] Hu, W., Shao, M. A., & Reichardt, K. (2010b). Using a New Criterion to Identify Sites

[25] Hu, W., Shao, M. A., Wang, Q. J., & Reichardt, K. (2009). Time Stability of Soil Water Storage Measured by Neutron Probe and the Effects of Calibration Procedures in a

[26] Hu, W., Tallon, L. K., & Si, B. C. (2012b). Evaluation of Time Stability Indices for Soil

[27] Huang, N. E., Shen, Z., Long, S. R., Wu, M. C. L., Shih, H. H., Zheng, Q., Yen, N. C., Tung, C. C., & Liu, H. H. (1998). The Empirical Mode Decomposition and the Hilbert

Temporal Stability of Soil Water Content. Geoderma, 158, 181-198.

for Mean Soil Water Storage Evaluation. Soil sci. soc. am. j., 74, 762-773.

water manage., 102, 66-73.

ments. J. hydrol., 253, 211-226.

3493-3511.

2897-2908.

hydrol., 207, 68-82.

oph., 11, 561-566.

hydrol., submitted.

Small Watershed. Catena, 79, 72-82.

Water Storage Upscaling. J. hydrol., accepted.

Arid Environment. Hydrol. process., 14, 1261-1277.

terns across Agricultural Fields. Catena, 73, 125-133.

## **References**


[14] Gao, X. D., Wu, P. T., Zhao, X. N., Shi, Y. G., & Wang, J. W. (2011). Estimating Spatial Mean Soil Water Contents of Sloping Jujube Orchards Using Temporal Stability. Agr. water manage., 102, 66-73.

2 Department of Natural Resource Sciences, McGill University, Ste-Anne-de-Bellevue, Que‐

[1] Biswas, A., & Si, B. C. (2011a). Identifying Scale-Specific Controls of Soil Water Stor‐ age in a Hummocky Landscape Using Wavelet Coherency. Geoderma, 165, 50-59. [2] Biswas, A., & Si, B. C. (2011b). Revealing the Controls of Soil Water Storage at Differ‐

[3] Biswas, A., & Si, B. C. (2011c). Scales and Locations of Time Stability of Soil Water

[4] Biswas, A., & Si, B. C. (2012). Identifying Effects of Local and Nonlocal Factors of Soil Water Storage Using Cyclical Correlation Analysis. Hydrol. Process., doi: 10.1002/

[5] Bosch, D. D., Lakshmi, V., Jackson, T. J., Choi. M., & Jacobs, J. M. (2006). Large Scale Measurements of Soil Moisture for Validation of Remotely Sensed Data: Georgia Soil

[6] Brocca, L., Melone, F., Moramarco, T., & Morbidelli, R. (2009). Soil Moisture Tempo‐ ral Stability over Experimental Areas in Central Italy. Geoderma, 148, 364-374.

[7] Cassel, D. K., Wendroth, O., & Nielsen, D. R. (2000). Assessing Spatial Variability in an Agricultural Experiment Station Field: Opportunities Arising from Spatial De‐

[8] Chen, Y. (2006). Letter to the Editor on Rank Stability or Temporal Stability. Soil sci.

[9] Comegna, V., & Basile, A. (1994). Temporal Stability of Spatial Patterns of Soil Water

[10] Cosh, M. H., Jackson, T. J., Bindlish, R., & Prueger, J. H. (2004). Watershed Scale Tem‐ poral and Spatial Stability of Soil Moisture and Its Role in Validating Satellite Esti‐

[11] Cosh, M. H., Jackson, T. J., Moran, S., & Bindlish, R. (2008). Temporal Persistence and Stability of Surface Soil Moisture in a Semi-Arid Watershed. Remote sens. environ.,

[12] Cosh, M. H., Jackson, T. J., Starks, P., & Heathman, G. (2006). Temporal Stability of Surface Soil Moisture in the Little Washita River Watershed and Its Applications in

[13] Farge, M. (1992). Wavelet transforms and their applications to turbulence. Annu. rev.

Satellite Soil Moisture Product Validation. J. hydrol., 323, 168-177.

Storage in a Cultivated Vesuvian Soil. Geoderma, 62, 299-310.

ent Scales in a Hummocky Landscape. Soil sci. soc. am. j., 75, 1295-1306.

Storage in a Hummocky Landscape. J. hydrol., 408, 100-112.

Moisture Experiment of 2003. J. hydrol., 323, 120-137.

pendence. Agron. j., 92, 706-714.

mates. Remote sens. environ., 92, 427-435.

soc. am. j., 70, 306.

112, 304-313.

of fluid mech., 24, 395-457.

bec, Canada

76 Advances in Agrophysical Research

**References**

hyp.8459.


Spectrum for Nonlinear and Non-Stationary Time Series Analysis. Proc. r. soc. a, 454, 903-995.

[41] Martínez-Fernández, J., & Ceballos, A. (2003). Temporal Stability of Soil Moisture in

Time Stability of Soil Water Content http://dx.doi.org/10.5772/52469 79

[42] Martínez-Fernández, J., & Ceballos, A. (2005). Mean Soil Moisture Estimation Using

[43] Mohanty, B. P., & Skaggs, T. H. (2001). Spatio-Temporal Evolution and Time-Stable Characteristics of Soil Moisture within Remote Sensing Footprints with Varying Soil

[44] North, G. R., Bell, T. L., Cahalan, R. F., & Moeng, F. J. (1982). Sampling errors in the estimation of empirical orthogonal functions. Mon. weather rev., 110, 699-706.

[45] Pachepsky, Y. A., Guber, A. K., & Jacques, D. (2005). Temporal Persistence in Vertical

[46] Perry, M. A., & Niemann, J. D. (2007). Analysis and Estimation of Soil Moisture at the

[47] Rehman, N., & Mandic, D. P. (2009). http://www.commsp.ee.ic.ac.uk/~mandic/

[48] Rehman, N., & Mandic, D. P. (2010). Multivariate Empirical Mode Decomposition.

[49] Rolston, D. E., Biggar, J. W., & Nightingale, H. I. (1991). Temporal Persistence of Spa‐

[50] Schneider, K,, Huisman, J. A., Breuer, L., Zhao, Y., & Frede, H. G. (2008). Temporal Stability of Soil Moisture in Various Semi-Arid Steppe Ecosystems and Its Applica‐

[51] Si, B. C. (2008). Spatial Scaling Analyses of Soil Physical Properties: A Review of

[52] Si, B. C., & Zeleke, T. B. (2005). Wavelet Coherency Analysis to Relate Saturated Hy‐ draulic Properties to Soil Physical Properties. Water resour. res., 41, W11424, doi:

[53] Starks, P. J., Heathman, G. C., Jackson, T. J., & Cosh, M. H. (2006). Temporal Stability

[54] Starr, G. C. (2005). Assessing Temporal Stability and Spatial Variability of Soil Water Patterns with Implications for Precision Water Management. Agr. water manage., 72,

[55] Tallon, L. K., & Si, B. C. (2004). Representative Soil Water Benchmarking for Environ‐

[56] Torrence, C., & Compo, G. P. (1998). A Practical Guide to Wavelet Analysis. Bull. am.

tial Soil-Water Patterns under Trickle Irrigation. Irrig. sci., 12, 181-186.

Distributions of Soil Moisture Contents. Soil sci. soc. am. j., 69, 347-352.

a Large-Field Experiment in Spain. Soil sci. soc. am. j., 67, 1647-1656.

Temporal Stability Analysis. J. hydrol., 312, 28-38.

Slope and Vegetation. Adv. water resour., 24, 1051-1067.

Catchment Scale Using Eofs. J. hydrol., 334, 388-404.

tion in Remote Sensing. J. hydrol., 359, 16-29.

of Soil Moisture Profile. J. hydrol., 324, 400-411.

mental Monitoring. J. environ. inform., 4, 581-590.

Spectral and Wavelet Methods, Vadose zone j., 7, 547-562.

research/emd.htm.

Proc. r. soc. a, 466, 1291-1302.

10.1029/2005WR004118.

meteor. soc., 79, 61-78.

223-243.


[41] Martínez-Fernández, J., & Ceballos, A. (2003). Temporal Stability of Soil Moisture in a Large-Field Experiment in Spain. Soil sci. soc. am. j., 67, 1647-1656.

Spectrum for Nonlinear and Non-Stationary Time Series Analysis. Proc. r. soc. a, 454,

[28] Huang, N. E., Wu, M. C. L., Long, S. R., Shen, S. S. P., Qu, W., Gloersen, P., & Fan, K. L. (2003). A Confidence Limit for the Empirical Mode Decomposition and Hilbert

[29] Huang, N. E., & Wu, Z. (2008). A Review on Hilbert-Huang Transform: Method and Its Applications to Geophysical Studies. Rev. geophys, 46, RG2006, doi:

[30] Ibrahim, H. M., & Huggins, D. R. (2011). Spatio-Temporal Patterns of Soil Water Stor‐ age under Dryland Agriculture at the Watershed Scale. J. Hydrol., 404, 186-197.

[31] Jacobs, J. M., Mohanty, B. P., Hsu, E. C., & Miller, D. (2004). SMEX02: Field Scale Var‐ iability Time Stability and Similarity of Soil Moisture. Remote sens. environ., 92,

[32] Jawson, S. D., & Niemann, J. D. (2007). Spatial Patterns from EOF Analysis of Soil Moisture at a Large Scale and Their Dependence on Soil, Land-Use, and Topographic

[33] Joshi. C., & Mohanty, B. P. (2010). Physical Controls of Near-Surface Soil Moisture across Varying Spatial Scales in an Agricultural Landscape During SMEX02. Water

[34] Joshi. C., Mohanty, B. P., Jacobs, J. M., & Ines, A. V. M. (2011). Spatiotemporal Analy‐ ses of Soil Moisture from Point to Footprint Scale in Two Different Hydroclimatic Re‐

[35] Kachanoski, R. G., & De Jong, E. (1988). Scale Dependence and the Temporal Persis‐ tence of Spatial Patterns of Soil Water Storage. Water resour. res., 24, 85-91.

[36] Korres, W., Koyama, C. N., Fiener, P., & Schneider, K. (2010). Analysis of Surface Soil Moisture Patterns in Agricultural Landscapes Using Empirical Orthogonal Func‐

[37] Kumar, P., & Foufoula-Georgiou, E. (1993). A Multicomponent Decomposition of Spatial Rainfall Fields: 1. Segregation of Large- and Small-Scale Features Using

[38] Kumar, P., & Foufoula-Georgiou, E. (1997). Wavelet Analysis of Geophysical Appli‐

[39] Lakshmi, V., Piechota, T., Narayan, U., & Tang, C. L. (2004). Soil Moisture as an Indi‐ cator of Weather Extremes. Geophys. res. lett., 31, L11401, doi:10.1029/2004GL019930.

[40] Lin, H. (2006).Temporal Stability of Soil Moisture Spatial Pattern and Subsurface Preferential Flow Pathways in the Shale Hills Catchment. Vadose zone j., 5, 317-340.

gions. Water resour. res., 47, W01508, doi:10.1029/2009WR009002.

Spectral Analysis. Proc. r. soc. a, 459, 2317-2345.

Properties. Adv. water resour., 30, 366-381.

tions. Hydrol. earth syst. sci., 14, 751-764.

cations. Rev. geoph., 35, 385-412.

Wavelet Transforms. Water resour. res., 29, 2515-2532.

resour. res., 46, W12503, doi: 10.1029/2010WR009152.

903-995.

78 Advances in Agrophysical Research

436-446.

10.1029/2007RG000228.


[57] Vachaud, G., Passerat De Silans, A., Balabanis, P., & Vauclin, M. (1985). Temporal Stability of Spatially Measured Soil Water Probability Density Function. Soil sci. soc. Am. J., 49, 822-828.

**Chapter 4**

**Model Averaging**

Asim Biswas and Bing Cheng Si

http://dx.doi.org/10.5772/52339

**1. Introduction**

**for Semivariogram Model Parameters**

Soil varies considerably from location to location (Nielsen et al., 1973). Knowledge on soil spa‐ tial variability is important in ecological modelling (Burrough, 1983; Corwin et al., 2006), envi‐ ronmental prediction (Trangmar et al., 1985), precision agriculture (Goderya, 1998), soil quality assessment (Heuvelink and Pebesma, 1999; McBratney et al., 2000), and natural resour‐ ces management. Adequate understanding of the variability in soil properties as a function of space and time is necessary for developing logical, empirical and physical models of soil and landscape processes (Burrough, 1993; Foussereau et al., 1993; Wilding et al., 1994). Geostatis‐ tics, a widely used approach, has been used to identify the spatial structure in the variability of soil attributes (Vieira, 2000; Carvalho et al., 2002; Vieira et al., 2002). Semivariance function characterizes the spatial continuity between points. When the semivariance is plotted against the lag distance or separation distance between points, the plot is called semivariogram (Fig. 1; McBratney and Webster, 1986; Isaaks and Srivastava, 1989). The structure of the semivario‐ gram is explained by three properties; the nugget, the sill and the range (Fig. 1). These spatial structures of semivariogram help in identifying autocorrelation and replicating samples, re‐ vealing dominant pattern in data series, identifying major ongoing processes (Si et al., 2007), designing experiments (Fagroud and van Meirvenne, 2002) and monitoring networks (Pra‐ kash and Singh, 2000), selecting proper data analysis method and interpreting data (Lambert et al., 2004), and assessing simulation and uncertainty analysis in a better way (Papritz and Du‐ bois, 1999). The semivariogram structures also help to quantify spatial dependence between observations. Modelling of observed semivariance values helps in predicting the spatial distri‐ bution of attribute values (Goovaerts, 1998). The spatial distribution of attribute values is very important in separating random noise in semivariance, and interpolation and mapping analy‐

sis such as kriging (Deutsch and Journel, 1998; Nielsen and Wendorth, 2003).

© 2013 Biswas and Si; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2013 Biswas and Si; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Additional information is available at the end of the chapter


## **Chapter 4**
