Published online 27 August 2007
Published in Soil Sci Soc Am J 71:1629-1635 (2007)
DOI: 10.2136/sssaj2006.0083
© 2007 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
PEDOLOGY
Statistical Methods for Evaluating Soil Salinity Spatial and Temporal Variability
Ahmed Douaika,*,
Marc Van Meirvenneb and
Tibor Tóthc
a Research Unit on Environ. and Conserv. of Natural Resources, National Institute of Agric. Res. (INRA), Avenue Hassan II BP 415, 10000 Rabat, Morocco
b Dep. of Soil Management and Soil Care, Ghent Univ., Coupure 653, 9000 Gent, Belgium
c Research Institute for Soil Science and Agricultural Chemistry, Hungarian Academy of Sciences, Herman Otto ùt 15, PO Box 35, 1525 II Budapest, Hungary
* Corresponding author (ahmdk{at}hotmail.com).
 |
ABSTRACT
|
|---|
Monitoring soil salinity requires knowledge of its magnitude and its spatial and temporal variability. To characterize the spatiotemporal variability of soil salinity in a native sodic grassland in the east of Hungary, we applied several statistical methods. Within a 25-ha study area, soil samples were taken repeatedly from 13 to 20 locations at 19 dates between November 1994 and June 2001 (with intervals between 2 and 9 mo). Electrical conductivity was measured both in the laboratory in a 1:2.5 soil/water suspension (EC2.5) and in situ using a four-electrode probe (ECa). These measurements were converted, via calibration regressions, into predicted EC2.5*, which were compared with EC2.5 in their ability to characterize the spatiotemporal variability of soil salinity. The temporal change in the mean soil salinity level between each subsequent two dates was evaluated using a paired t-test, a test of significance of the regression parameters based on the concept of temporal stability, and a temporal mean shift test. The static–dynamic (uniform–nonuniform) nature of the temporal change in the spatial pattern of soil salinity between two dates was evaluated using the same concept of temporal stability and a spatial shift test. For either type of temporal change (mean level or spatial pattern), the methods agreed for some pairs of dates and did not for others, but these differences were partly due to differences in data input. The method to use depends on the data availability and the aim of the study. The joint use of temporal stability and temporal mean shift and spatial shift tests could result in a drastically reduced sampling effort.
Abbreviations: EC2.5, soil electrical conductivity determined from a 1:2.5 soil/water suspension EC2.5*, electrical conductivity in a 1:2.5 soil/water suspension predicted from apparent electrical conductivity using calibration regression equation ECa, apparent soil electrical conductivity ECe, electrical conductivity of a water-saturated soil paste MD, mean difference
 |
INTRODUCTION
|
|---|
To control its harmful effects, soil salinity needs to be monitored in space and time. This requires knowledge of its magnitude, temporal dynamics, and spatial variability. The magnitude of soil salinity is conventionally assessed in the laboratory by determining the electrical conductivity of a water-saturated soil paste (ECe), either in an extracted suspension or in the paste directly. This procedure, however, is time consuming and expensive.
Alternatively, soil salinity can be evaluated by measuring the apparent electrical conductivity (ECa) in the field using electrode probes or a soil sensor like the EM38 (Geonics Ltd, Mississauga, ON, Canada). This approach is faster and cheaper, and allows a more intensive surveying. Nevertheless, it still requires the collection of soil samples for analysis in the laboratory to establish calibration equations linking ECa to ECe. Different approaches are used but the more statistically rigorous is the stochastic field-calibration technique (Rhoades et al., 1999). It is based on regression models linking ECa to ECe for areas uniformly managed and homogeneous for all but soil salinity factors. For example, Halvorson and Rhoades (1974) found a correlation coefficient between both variables of 0.98 when the soil profile was at or near field capacity, and of 0.96 when it was drier. More details about describing and interpreting soil salinity from ECa can be found elsewhere (e.g., Rhoades et al., 1999). Clay et al. (2001) evaluated the spatial variability of ECa and found that it was positively correlated to ECe, soil water, and clay content. Corwin et al. (2006), studying the space–time variation of soil quality based on an ECa directed sampling, found that salinity levels were reduced by 13% in the 0- to 60-cm depth from 1999 to 2002, mainly due to leaching. When dynamic soil properties, such as salt concentration or water content, dominate ECa measurements, their spatial patterns are less time stable than ECa measurements that are mainly influenced by static properties like soil texture (Corwin, 2005).
As soil salinity can be variable in space and time, monitoring requires a lot of measurements. The spatial pattern of soil salinity can persist with time (Castrignanò et al., 1994), however, so that locations can be identified where deviations of soil salinity from the field mean value are nearly constant at all times. This means that a reduced number of selected locations would suffice to characterize the mean soil salinity status of the field. In such a situation, the sampling effort could be reduced substantially (Vachaud et al., 1985). Temporal persistence, also known as temporal stability, is defined as the time-invariant association between spatial location and classic statistical parameters (Vachaud et al., 1985). To check for its presence, Vachaud et al. (1985) introduced two tests. The first one is based on the Spearman rank correlation. It consists of the determination of the rank correlation coefficient between sampling dates. The second test is based on the relative differences between sampling dates. Kachanoski and de Jong (1988) applied the second approach and showed that a good test for the temporal persistence was the significance and magnitude of the Pearson correlation coefficient and the parameters of a simple linear regression fitted to values observed at two consecutive dates. The significance test of the regression intercept is a check on whether the mean soil salinity level has changed (either increased or decreased) or not between two dates, while the test of significance of the regression slope is a check on whether the spatial pattern of soil salinity between two dates was static (uniform, i.e., changed by the same amount for all locations) or dynamic (nonuniform, i.e., the amount of change was different from one location to another). Although the temporal stability has been used by many researchers, most of their studies have been applied only to soil water properties (Van Wesenbeeck et al., 1988; Da Silva et al., 2001; Martínez-Fernández and Ceballos, 2003; Petrone et al., 2004). Castrignanò et al. (1994) evaluated the temporal stability of ECe, the Na content, and the Na adsorption ratio from a field of 2.8 ha (at 28 locations) during 2 yr (at eight dates). They were able to identify a small number of sampling locations to estimate the mean soil salinity and sodicity at any time during the 2-yr period. Also, they found that the spatial pattern of these soil properties persisted during the time period considered.
The paired t-test, comparing the mean differences, is another way to check if the mean soil salinity level has significantly changed between two dates. This method, however, does not allow a check of the static–dynamic nature of the spatial pattern of soil salinity. This approach was used by Kenny et al. (2002) to check the temporal trend of the mean thickness of the Ap horizon.
There is a third way to check the significance of a temporal change in the mean soil salinity level as well as the static–dynamic nature of its spatial pattern. For this method, Lesch et al. (1998) introduced a different procedure by combining ECa and ECe. Their approach is based on two statistical tests and allows a check of two kinds of change. The first test, called the temporal mean shift, is used to check if the observed mean value at the second date and the mean of the estimated values (obtained from the calibration equation of the first date and ECa measurements at the second date) are significantly different, i.e., if the mean soil salinity level has changed. The second test, called the spatial shift, allows detection of the existence of a dynamic (nonuniform) spatial variation, i.e., a change in soil salinity that is different from one location to another. Otherwise, the spatial variation was static, meaning either no change or a constant change that occurred for all locations between the two sampling dates.
The main objective of this research was to evaluate several statistical methods to characterize the spatiotemporal variability of soil salinity either by direct electrical conductivity measurements or from predictions obtained from ECa observations. The specific objectives were to evaluate: (i) the temporal change in the mean soil salinity level between two sampling dates using a paired t-test, a significance test of the regression parameters based on the concept of temporal stability, and a temporal mean shift test; (ii) the static–dynamic (uniform–nonuniform) nature of the temporal change of the spatial pattern of soil salinity between two dates using the same concept of temporal stability and a spatial shift test; and (iii) the opportunity to reduce the sampling effort for monitoring soil salinity with time depending on whether the mean soil salinity level has changed or not and whether its spatial pattern is static or dynamic.
 |
MATERIALS AND METHODS
|
|---|
Study Site
Our study site was a native sodic grassland of about 25 ha located in the Hortobagy National Park in eastern Hungary with central coordinates 47°30' N and 21°30' E (Fig. 1
). This area has been described in several studies, providing details on soil properties (Tóth and Jozefaciuk, 2002) and the spatial and temporal variation of soil salinity (Douaik et al., 2004; Douaik, 2005; Tóth et al., 2002). Although it was not our aim to characterize the pedological or chemical aspects of the soil salinity of this area, some of its characteristics will be discussed briefly.
The study area is a part of the largest contiguous salt-affected landscape in central Europe. The topography is nearly flat, with a mean elevation of 89 m above sea level. The typical soils of the study site are classified as Typic Natrustolls and Typic Natraquolls according to Soil Taxonomy and they have a characteristic morphology. The A/E (eluvial) horizon is mostly salt free, with a pH of about 7. Its texture is loam with a weakly developed platy structure with many roots. In contrast, the B (illuvial) horizon is characterized by a strong alkalinity, a finer texture, typically clay loam, a columnar structure, and a very low hydraulic conductivity. The C horizon is less affected by salt accumulation, resembling the original loess parent material, but the Natraquolls show many redoximorphic features in this horizon.
The region is characterized by a wide range in monthly average temperatures (23.6°C) and a mean annual precipitation of 524 mm. The wettest month is June (67 mm), while January is the driest (26 mm). The drought index, i.e., the ratio of potential evaporation to precipitation, exceeds unity for most of the year (March–October), hence the salinization susceptibility of the region. Since the region is a recharge area of saline groundwater originating from rock weathering in the surrounding mountains, the groundwater is the main source of salts.
Data Description
Field work was conducted at 19 dates between November 1994 and June 2001, with temporal intervals ranging from 2 to 9 mo. The ECa was measured at up to 413 locations according to a pseudo-regular grid of 25 by 25 m, using the four-electrode probe apparatus of Rhoades and Miyamoto (1990). Based on the principles of the effect of inserting the electrodes (Rhoades et al., 1999) in the sodic soil, two insertion depths were initially used to characterize the 0- to 20- and 0- to 40-cm depths. Due to the temporally changing soil moisture content (ranging from puddles to cracking surface), however, instrumental data collected with 8-cm-deep insertion (measured from the air–litter contact of the grassland) showed great scatter and were ignored, whereas data obtained with 13-cm-deeper insertion showed significant correlation with the soil salinity of the 0- to 40-cm layer. An electric current with intensity Y is induced between the outer electrodes and the potential, E, is measured between the inner electrodes. The soil resistance Rs is obtained from
 | [1] |
which was converted to ECa using a cell constant determined experimentally.
Simultaneously, a limited number of soil samples, between 13 and 20 depending on the date, were collected following a spatial sampling algorithm (Lesch et al., 1995) (Fig. 2
). Normally, 20 calibration locations were sampled but at seven dates, some of the calibration locations (1–7) could not be sampled due to water logging, instrument malfunction, or vandalism of the reference marks since the area is an open pasture. On these samples, EC2.5 was measured in the laboratory in a 1:2.5 soil/water suspension and converted to a reference temperature of 25°C, which we considered as a simple proxy for ECe (Soil and Plant Analysis Council, 1992). The sampling locations were selected following a spatial response surface sampling design (based on the 413 ECa observations) and spanned the whole extent of the study area, which optimized the estimation of the calibration model. Soil samples were collected at four depths from 0 to 40 cm with an increment of 10 cm. In this study, we considered only the ECa of the depth interval 0 to 40 cm and the mean EC2.5 for the four depths. The average EC2.5 values ranged between 1.6 mS m–1, with a standard deviation of 0.63 mS m–1 (March 1997), and 3.3 mS m–1, with a standard deviation of 2.17 mS m–1 (September 1998) (Douaik, 2005).

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 2. Study area with locations where apparent electrical conductivity was measured and where samples were taken to determine electrical conductivity in a 1:2.5 soil/water suspension (in this figure all 20 locations are indicated, but due to practical limitations, locations 1 to 7 could not be sampled at some of the 19 dates).
|
|
The EC2.5 is reported at a standard temperature of 25°C, while ECa was reported at the prevailing temperature in the field. To avoid a proportional shift due to differences in temperature and moisture content, however, ECa was rescaled (using calibration equations) into a predicted EC2.5* separately for each date (Table 1). Thus we could evaluate the characterization of the temporal stability of spatial patterns in the observed (EC2.5) and predicted (EC2.5*) soil salinity at the same 13 to 20 locations. This allowed us to evaluate how well the fast field method using electrodes can replace laboratory analyses in characterizing the temporal stability of soil salinity.
View this table:
[in this window]
[in a new window]
|
Table 1. Linear regression parameters of calibration equations relating apparent soil electrical conductivity (x) to soil electrical conductivity determined from a 1:2.5 soil/water suspension (y) for every sampling date to determine the predicted soil salinity.
|
|
Maps of EC2.5* at different dates have been reported elsewhere (Douaik, 2005).
Statistical Methods
The difference in mean soil salinity level between two dates can be checked using three methods: the paired t-test, the significance test of the regression parameters of the relative differences based on the concept of temporal stability, and the temporal mean shift test. It is also interesting, however, to evaluate the static–dynamic nature of the temporal change of the spatial pattern of soil salinity. Two evaluation methods are available for this: the same concept of temporal stability and the spatial shift test.
Temporal Change of the Mean Soil Salinity Level
Paired t-Test.
The paired t-test (McClave and Sincich, 2000) tests if the mean difference (MD) is significant. Therefore the differences between the soil salinity (generally identified as EC) on two consecutive dates are calculated:
 | [2] |
with
j=(1/n)
i=1nECij and ECij the EC measured at location i (i = 1, ..., n) and date j (j = 1, ..., l). Then, a t-test statistic is computed and compared with tabulated t values or the corresponding probability of the t-test statistic is compared with a given level of significance (generally 0.05) to check if the difference of the means is significant.
Temporal Stability Test.
The concept of temporal stability was introduced by Vachaud et al. (1985). It refers to the tendency of measurements of a soil property at different locations in space to maintain their relative ranking with time. One way to evaluate this temporal stability is to calculate the relative difference (
ij):
 | [3] |
Kachanoski and De Jong (1988) refined the definition of temporal stability by showing that it exists if the relative differences remain constant between two dates. This implies that
 | [4] |
Thus, to evaluate the temporal stability of soil salinity, a linear regression between EC at two different dates can be modeled with an intercept I, and a slope S given by
 | [5] |
The temporal change in the mean soil salinity level can be checked by testing if I = 0. If this hypothesis is rejected, we can assume that the mean soil salinity level has increased (I > 0) or decreased (I < 0) between the two dates. On the other hand, if the hypothesis can't be rejected, the mean soil salinity is considered unchanged between these dates.
Temporal Mean Shift Test.
The temporal stability test investigates the temporal changes of EC2.5 and EC2.5* (and thus indirectly ECa) separately. Lesch et al. (1998) discussed another approach that can test the change in the mean soil salinity level, as well as the static–dynamic nature of the change in the soil salinity spatial pattern between two sampling dates, by integrating both soil properties. Their approach is based on the estimation of the regression parameters (intercept Ij and slope Sj) of the model
 | [6] |
using the data from date j (j = 1, ..., l) observed at n locations. Next, the EC2.5 of date j+1 is predicted from ECa measurements at m locations (i = 1, ..., m), with m
n, using the regression parameters of the previous date:
 | [7] |
Then, the MD is obtained between these predictions and the observed EC2.5 at the m locations of the second date:
 | [8] |
with
j+1=(1/m)
i=1mEC2.5i,j+1, and similarly for the predicted value.
Together with the regression model design matrix, an approximate t-test (with n – 2 degrees of freedom) can be derived to test if MD = 0, i.e., if there is a change in the mean soil salinity level (Lesch et al., 1998).
Static–Dynamic Nature of the Temporal Change of the Soil Salinity Spatial Pattern
Temporal Stability Test.
If there is no, or uniform, change in the EC of all the locations between two sampling dates, we expect that S = 1 (Eq. [5]; Kachanoski and de Jong, 1988). Then, the magnitude of the change is given by I, so when I = 0 there is no change. In this case, the spatial pattern of the soil salinity is said to be static. On the other hand, if S is significantly different from 1, the spatial pattern of the soil salinity can be considered to be dynamic, i.e., the magnitude of change in the soil salinity between two dates is nonuniform and is different from one location to another.
Spatial Shift Test.
The nature of the spatial variation of the temporal change (static–dynamic) is tested by computing a statistic
, which is a function of the normalized squared differences, the design matrix, and the regression model mean squared error estimate. The statistic
is compared to an F distribution with m – 1 and n – 2 degrees of freedom (Lesch et al., 1998).
Combining Tests for Mean Level and Spatial Pattern of Soil Salinity
Temporal Stability Test.
Based on whether I is significantly different from 0, S is significantly different from 1, or both, four possible scenarios can be distinguished (Grant et al., 2004):
- Scenario 1 (I = 0 and S = 1): there is no net change between the two sampling dates
- Scenario 2 (I
0 and S = 1): the mean soil salinity level has changed and there is a static (uniform) change in the soil salinity spatial pattern
- Scenario 3 (I = 0 and S
1): the mean soil salinity level has not changed and the change in the soil salinity spatial pattern is dynamic (nonuniform) but, when averaged, this change is not significant
- Scenario 4 (I
0 and S
1): the mean soil salinity level has changed and the change in the soil salinity spatial pattern is dynamic (nonuniform)
Temporal Mean Shift and Spatial Shift Tests.
Again, four scenarios are possible:
- Scenario 1 (MD = 0 and
= 0): there is no net change with time in the soil salinity
- Scenario 2 (MD
0 and
= 0): the mean soil salinity level has changed and there is a static (uniform) change in the soil salinity spatial pattern
- Scenario 3 (MD = 0 and
0): the mean soil salinity level has not changed and the change in the soil salinity spatial pattern is dynamic (nonuniform); again, when averaged, this change is not significant
- Scenario 4 (MD
0 and
0): the mean soil salinity level has changed and the change in the spatial pattern of the soil salinity is dynamic (nonuniform)
The usefulness of the temporal mean shift and spatial shift tests lies in the possibility of using only the pairs of ECa–EC2.5 values from the first sampling date, along with the EC2.5 values from the second date. If it is found that there was no or a static (uniform) change between the two dates, then the regression equation computed using the data from the first date can be used to predict EC2.5* for the second date using a larger ECa survey data set. If a dynamic spatial variation was detected, however, then a new regression equation should be estimated using the ECa–EC2.5 pairs of values from the second date.
The paired t-test and the temporal stability test were performed for both EC2.5 and EC2.5*, while the temporal mean shift (significance of MD) and the spatial shift tests (significance of
) were applied only to EC2.5, as these tests integrate the ECa.
 |
RESULTS AND DISCUSSION
|
|---|
Calibration Equations
All the coefficients of determination (r2) of the calibration equations, EC2.5 = I + SECa, for each of the observed dates (Table 1) were highly significant and ranged between 0.69 and 0.94, indicating that the regression explained between 69 and 94% of the variance within the data. The intercept (I) values ranged between –0.10 and 1.29 dS m–1, while the slope (S) values varied between 0.41 and 2.56. These fluctuations reflect the differences in field conditions (temperature and moisture content) when ECa was observed, but no clear seasonal effect could be found. Also, the RMSE, reflecting the deviation from the regression line, was generally small. It varied between 0.26 and 1.17 dS m–1. Therefore, it could be concluded that, generally, ECa was found to be closely and linearly linked with EC2.5.
It was assumed that EC2.5 could be used as a simple proxy of ECe, reflecting soil salinity. Since EC2.5 was determined in the laboratory for each date, we considered this soil property to be the reference. The predicted soil salinity (EC2.5*) was obtained from field measurements of ECa using an electrode configuration and a fitted calibration curve. The goodness-of-fit of the calibration equations (Table 1) indicate that generally there was a strong relationship between both variables and that EC2.5 could be reliably predicted from ECa.
Correlation Coefficients
The Pearson correlation coefficients (r) between consecutive sampling dates for both EC2.5 and EC2.5* were generally quite strong and all values were significant (Table 2). They ranged from 0.58 to 0.94 for EC2.5 and from 0.81 to 0.98 for EC2.5*. These values indicate the presence of moderate to strong temporal stability across the observed time intervals. The correlation coefficients of both variables were very similar and the differences were negligible except for the first pair of dates, where the difference was moderate (0.27).
View this table:
[in this window]
[in a new window]
|
Table 2. Pearson correlation coefficients (r) and results of the paired t-test of the observed (EC2.5) and predicted (EC2.5* derived from calibration equations relating apparent electrical conductivity to EC2.5) electrical conductivity in a 1:2.5 soil/water suspension for 18 pairs of observation dates (significant values are in italic and the level of significance was 0.05).
|
|
Temporal Change of the Mean Soil Salinity Level
Paired t-Test.
The paired t-test was applied to both EC2.5 and EC2.5* (Table 2). The MD in EC2.5 between two consecutive dates (Eq. [2]) was significantly different from zero for only four out of the 18 pairs of dates. Of these four pairs, the mean value decreased only once (by 1.46 dS m–1), while it increased for the other three pairs (by 0.27–1.70 dS m–1). The results were strongly different for EC2.5*. For this variable, the hypothesis that the mean difference was equal to zero was rejected for 11 pairs of dates. The significant mean differences decreased for seven pairs and increased for four pairs. The decrease ranged from 0.18 to 1.45 dS m–1, while the increase varied between 0.15 and 1.70 dS m–1. The four pairs of dates where the mean difference in EC2.5 was significant are included in the case of EC2.5* and the magnitudes of the mean differences of both variables were very similar.
Temporal Stability.
From the four possible scenarios of the temporal stability test, the second and the fourth represent the cases where the mean change is statistically significant, whereas for the first and third scenarios, it is not. Regarding EC2.5 (Table 3), the mean level changed significantly for seven pairs of dates, all showing an increase between 0.37 and 1.05 dS m–1, except one pair for which there was a decrease of 0.57 dS m–1. The mean EC2.5* level changed significantly for nine pairs of dates, with an increase between 0.44 and 0.77 dS m–1 for six pairs and a decrease between 0.57 and 0.88 dS m–1 for the remaining three pairs of dates. Again, as with the paired t-test, all the pairs of dates for which the mean EC2.5 level has changed are found for EC2.5*.
View this table:
[in this window]
[in a new window]
|
Table 3. Regression parameters between consecutive values of the observed (EC2.5) and predicted (EC2.5* derived from calibration equations relating apparent electrical conductivity to EC2.5) electrical conductivity in a 1:2.5 soil/water suspension for 18 pairs of observation dates (significant values are in italic and the level of significance was 0.05).
|
|
Temporal Mean Shift Test.
The temporal mean shift test compared the mean value for the observed EC2.5 of one date to the mean of the predicted values based on the regression modeling between EC2.5 and ECa from the previous date (Table 1). The obtained MDs (Eq. [8]) were computed, together with their probability of significance (second and third columns of Table 4). The MD was significantly different for only three pairs of dates, with a decrease of 1.45 dS m–1 for one pair and an increase of 0.37 and 1.70 dS m–1 for the other two pairs of dates. These correspond to the same pairs of dates (no. 1, 11, and 12) for which the MD, based on the paired t-test, was found to be significant (Table 2), but this test indicated one more pair (no. 8) with a significant change in the mean soil salinity level.
View this table:
[in this window]
[in a new window]
|
Table 4. Test of significance of the temporal mean shift and the spatial shift in the observed electrical conductivity in a 1:2.5 soil/water suspension for 18 pairs of observation dates (significant values are in italic and the level of significance was 0.05).
|
|
Static–Dynamic Nature of the Temporal Change of the Spatial Pattern of Soil Salinity
Temporal Stability Test.
The first and second scenarios of the temporal stability test indicate the presence of a static (uniform) change in the soil salinity spatial pattern between two sampling dates, while the third and fourth scenarios indicate a dynamic (nonuniform) change. It can be observed (Table 3) that the spatial pattern of EC2.5 was static for eight pairs of dates, whereas it was dynamic for the remaining 10 pairs of dates. The spatial pattern of EC2.5* was static and dynamic for the same number of pairs of dates (nine). The spatial pattern of both variables agreed for 15 pairs of dates, with a static spatial change for seven pairs and a dynamic spatial change for eight pairs.
Spatial Shift Test.
As for temporal stability, in the first and second scenarios of the spatial shift test, the change in the soil salinity spatial pattern can be considered static, whereas the third and fourth scenarios show a dynamic soil salinity spatial pattern. The
statistic was found to be significant, thus the spatial pattern of EC2.5 was dynamic (nonuniform) for only three pairs of dates (the last two columns of Table 4, no. 8, 11, and 18). For the remaining pairs (15), the spatial pattern can be considered static, i.e., had not changed or changed uniformly by the same amount for all locations.
There was a clear difference in the ability of the different methods to characterize the spatiotemporal variability of the observed (EC2.5) and predicted (EC2.5*) soil salinity. These differences are mainly due to different data input and data modeling. Regarding the data input, the observed soil salinity was measured directly in the laboratory following standard methods, whereas the predicted soil salinity was obtained from rescaling of ECa into soil salinity using calibration equations. Since the ECa is a function of many soil properties such as clay content, water content, and total soluble salts, the predicted electrical conductivity (EC2.5*) is only an indirect measurement of soil salinity and may fluctuate depending on the level of the clay and water content of the study site. This fluctuation is rather limited within our study area, so it can be expected that salinity dominates the ECa measurements.
Regarding the data modeling, the paired t-test compared the means of either the observed or the predicted soil salinity data between two sampling dates. The temporal stability test determined the slope and the intercept of the measured data at two consecutive dates and tested if these parameters were significantly different from one and zero, respectively. The temporal mean shift and the spatial shift methods integrated soil salinity and ECa data measured at a first date by determining the calibration equation linking these two soil properties, then used this regression to convert ECa into predicted soil salinity, and finally computed the differences between the predicted values and those observed at the second date. Although the data input and modeling are different for the methods, there was some agreement, since some pairs of dates were identified by all methods as dates with a significant temporal change in the mean level and in the spatial pattern.
 |
CONCLUSIONS
|
|---|
A strong relationship was found between EC2.5 and ECa for all 19 measurement dates, with large and significant r2 values and small RMSEs. This confirms the potential of the four-electrode probe to monitor the temporal change in soil salinity in native sodic grassland in every season. When ECa was converted into EC2.5*, however, some differences with EC2.5 were found.
The concept of temporal stability requires a moderate number of measurements, but once it is checked, the number of future measurements can be reduced to locations with measurements representing important statistics such as the average value. Likewise, based on the spatial shift test, if it was found that there was a static (uniform) spatial variation between two dates, the sampling effort could be limited to the ECa survey for the first date, while the laboratory analysis could be done for two dates but at a reduced number of locations. So the joint use of the concept of temporal stability and the temporal mean shift and spatial shift tests could result in a drastically reduced sampling effort. In another study using the same data, Douaik et al. (2006) found that it was sufficient to measure soil salinity at just two locations and nine dates to characterize reliably the average field soil salinity of this area.
The statistical method to use will depend on the data availability and the aim of the study. If only soil salinity or ECa data are available, there is only one possibility for checking the spatial pattern of the soil salinity: the concept of temporal stability. For checking the temporal change in the average soil salinity, there is a choice between the paired t-test and the concept of temporal stability. If, on the other hand, both soil salinity and ECa data are available, there are more options: the paired t-test and the concept of temporal stability for checking the temporal change in the average level separately for both soil properties or the temporal mean shift test for checking the temporal change in the average level combining both soil properties. Regarding the check of the temporal stability of the spatial patterns, the same choices are available except for the paired t-test. Finally, the results of this work can be extended to evaluate the temporal change (in the mean level as well as in the static–dynamic nature of the spatial pattern) of any other soil property related to ECa, such as water content or organic matter content.
 |
ACKNOWLEDGMENTS
|
|---|
A. Douaik was financially supported by the Ministry of the Flemish Community (Belgium) through his Ph.D. scholarship. T. Tóth acknowledges the material assistance of the Bilateral Flemish–Hungarian Research Foundation and the grants HNSF T37731 and OM-00124/2001.-NKPF/4. The linguistic improvement of David Simpson, Ghent University, is acknowledged.
 |
NOTES
|
|---|
All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permission for printing and for reprinting the material contained herein has been obtained by the publisher.
Received for publication February 22, 2006.
 |
REFERENCES
|
|---|
- Castrignanò, A., G. Lopez, and M. Stelluti. 1994. Temporal and spatial variability of electrical conductivity, Na content and sodium adsorption ratio of saturation extract measurements. Eur. J. Agron. 3:221–226.
- Clay, D.E., J. Chang, D.D. Malo, C.G. Carlson, C. Reese, S.A. Clay, M. Ellsbury, and B. Berg. 2001. Factors influencing spatial variability of soil apparent electrical conductivity. Commun. Soil Sci. Plant Anal. 32:2993–3008.[CrossRef][Web of Science]
- Corwin, D.L. 2005. Geospatial measurements of apparent soil electrical conductivity for characterizing soil spatial variability. p. 585–615. In J. Alvarez Benedi and R. Munoz-Carpena (ed.) Soil-water-solute process characterization: An integrated approach. CRC Press, Boca Raton, FL.
- Corwin, D.L., S.M. Lesch, J.D. Oster, and S.R. Kaffka. 2006. Monitoring management-induced spatio-temporal changes in soil quality through soil sampling directed by apparent electrical conductivity. Geoderma 131:369–387.[CrossRef][Web of Science]
- Da Silva, A.P., A. Nadler, and B.D. Kay. 2001. Factors contributing to temporal stability in spatial patterns of water content in the tillage zone. Soil Tillage Res. 58:207–218.[CrossRef]
- Douaik, A. 2005. Evaluation of the space–time variability of soil salinity by statistical, geostatistical and Bayesian maximum entropy methods. Ph.D. diss. Ghent Univ., Belgium.
- Douaik, A., M. Van Meirvenne, and T. Tóth. 2006. Temporal stability of spatial patterns of soil salinity determined from laboratory and field electrolytic conductivity. Arid Land Res. Manage. 20:1–13.[CrossRef]
- Douaik, A., M. Van Meirvenne, T. Tóth, and M.L. Serre. 2004. Space–time mapping of soil salinity using probabilistic Bayesian maximum entropy. Stoch. Environ. Res. Risk Assess. 18:219–227.
- Grant, L., M. Seyfried, and J. McNamara. 2004. Spatial variation and temporal stability of soil water in a snow-dominated, mountain catchment. Hydrol. Process. 18:3493–3511.[CrossRef]
- Halvorson, A.D., and J.D. Rhoades. 1974. Assessing soil salinity and identifying potential saline-seep areas with field soil resistance measurements. Soil Sci. Soc. Am. Proc. 38:576–581.
- Kachanoski, R.G., and E. de Jong. 1988. Scale dependence and the temporal persistence of spatial patterns of soil water storage. Water Resour. Res. 24:85–91.[CrossRef]
- Kenny, E.A., J.W. Hall, and C. Wang. 2002. Temporal trends in soil physical properties at a soil quality benchmark site in British Columbia, Canada. In Proc.World Congr. Soil Sci., 17th, Bangkok, Thailand [CD-ROM]. 14–20 Aug. 2002. Soil and Fert. Soc. of Thailand, Bangkok.
- Lesch, S.M., J. Herrero, and J.D. Rhoades. 1998. Monitoring for temporal changes in soil salinity using electromagnetic induction techniques. Soil Sci. Soc. Am. J. 62:232–242.[Abstract/Free Full Text]
- Lesch, S.M., D.J. Strauss, and J.D. Rhoades. 1995. Spatial prediction of soil salinity using electromagnetic induction techniques: 2. An efficient spatial sampling algorithm suitable for multiple linear regression model identification and estimation. Water Resour. Res. 31:387–398.[CrossRef]
- Martínez-Fernández, J., and A. Ceballos. 2003. Temporal stability of soil moisture in a large-field experiment in Spain. Soil Sci. Soc. Am. J. 67:1647–1656.[Abstract/Free Full Text]
- McClave, J.T., and T. Sincich. 2000. Statistics. 8th ed. Prentice Hall, Upper Saddle River, NJ.
- Petrone, R.M., J.S. Price, S.K. Carey, and J.M. Waddington. 2004. Statistical characterization of the spatial variability of soil moisture in a cutover peatland. Hydrol. Process. 18:41–52.[CrossRef]
- Rhoades, J.D., F. Chanduvi, and S.M. Lesch. 1999. Soil salinity assessment: Methods and interpretation of electrical conductivity measurements. Irrig. Drain. Pap. 57. FAO, Rome.
- Rhoades, J.D., and S. Miyamoto. 1990. Testing soils for salinity and sodicity. p. 428–433. In R.L. Westerman (ed.) Soil testing and plant analysis. 3rd ed. SSSA Book Ser. 3. SSSA, Madison, WI.
- Soil and Plant Analysis Council. 1992. Handbook of reference methods for soil analysis. Univ. Stn., Athens, GA.
- Tóth, T., and G. Jozefaciuk. 2002. Physicochemical properties of a solonetzic toposequence. Geoderma 106:137–159.[CrossRef][Web of Science]
- Tóth, T., L. Kuti, I. Forizs, S. Kabos, and A. Douaik. 2002. Spatial and temporal aspects of soil salinization in a sodic grassland. p. 276–288. In A. Fez Cano, R. Ortiz Silla, and A.R. Mermut (ed.) Sustainable use and management of soils in arid and semiarid regions. Int. Symp., Cartegena, Spain. 22–26 Sept. 2002. Vol. 1. Univ. of Murcia, Spain.
- Vachaud, G., A. Passerat De Silans, P. Balabanis, and M. Vauclin. 1985. Temporal stability of spatially measured soil water probability density function. Soil Sci. Soc. Am. J. 49:822–828.[Abstract/Free Full Text]
- Van Wesenbeeck, I.J., R.G. Kachanoski, and D.E. Rolston. 1988. Temporal persistence of spatial patterns of soil water content in the tilled layer under a corn crop. Soil Sci. Soc. Am. J. 52:934–941.[Abstract/Free Full Text]