SCHRIFTENREIHE Institut für Pflanzenernährung und Bodenkunde ...
SCHRIFTENREIHE Institut für Pflanzenernährung und Bodenkunde ...
SCHRIFTENREIHE Institut für Pflanzenernährung und Bodenkunde ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
cross-semivariograms (factorial kriging analysis). The LMC is a set of auto- and<br />
cross-semivariogram models in which all its semivariograms are linear<br />
combinations of the same set of elementary structures. A LMC with k = 1, . . ., q<br />
structures may be written as:<br />
g (h)<br />
k<br />
1 1 k k<br />
q q<br />
γ h)<br />
= b g ( h)<br />
+ b g ( h)<br />
+ ⋅⋅<br />
⋅ + b g ( h)<br />
(3)<br />
i , j ( i,<br />
j<br />
i,<br />
j<br />
i,<br />
j<br />
k<br />
b i,<br />
j<br />
where is the partial sill for the i,jth semivariogram for structure k, while<br />
represents the type of semivariogram model (i.e. spherical, exponential,<br />
1<br />
etc.) for structure k. The first structure g ( h)<br />
represents the nugget effect model.<br />
The type of semivariogram model used in Eq. 3 was based on the experimental<br />
semivariograms of the variables (standardized to unit variance and zero mean)<br />
and knowledge of the main geological and anthropogenic factors. Modified<br />
kriging programs in ArcGIS were used to fit the LMC, in which semivariograms<br />
were fitted with appropriate model functions using the maximum likelihood<br />
cross-validation method (Samper and Carrera, 1990).<br />
The scale-dependent correlations between water content and other<br />
properties was determined from the structural correlation coefficients, ρ k i.j<br />
k<br />
bi<br />
, j<br />
k k<br />
i , ib<br />
j , j<br />
k<br />
ρ i,<br />
j =<br />
b<br />
(4)<br />
which are calculated from the partial sill value, , of the<br />
cross-semivariogram model between i and j and the two partial sill values,<br />
k<br />
b j,<br />
j<br />
and , of semivariogram models for i and j, respectively (Casa and<br />
Castrignanò, 2007). Furthermore, constructed variation contributions of each<br />
component to the total variation were calculated based on Eq. 3.<br />
3. Results and discussion<br />
Grazing intensity influenced a range of soil properties listed in Table 1. The<br />
heavily grazed plot (HG) was the most dense, had the least carbon, and the<br />
greatest shear strength. In contrast, the plot protected from grazing for 25 yr (UG<br />
79) was the least dense, had the most carbon and the smallest shear strength.<br />
There appears to be consistent trends between grazing intensity and many of<br />
46<br />
k<br />
b i,<br />
j<br />
k<br />
b i,<br />
i