Abstract:
This research was conducted in Khanikan forests located in lowland of Mazandaran province (North of Iran). Eighteen profiles were dug and several chemical, physical and biological soil properties were investigated. The soil properties evaluated were soil pH, bulk density, saturation moisture content, electrical conductivity, organic carbon, total nitrogen, cation exchangeable capacity, available phosphorous, soil texture, calcium carbonate content, number and biomass of earthworms, litter carbon and litter nitrogen. Principal Component Analysis (PCA) was used to identify the variation of soil properties. PCA, a technique which reduces the dimensionality of multivariate data by removing Interco relations among variables, has a number of useful applications in forest researches. The results showed significant relationships between some soil factors with PC1 and PC2 axes, also, among different soil factors, the distribution of forest types was most strongly controlled with some soil characteristics such as acidity, bulk density, texture, phosphorous, organic carbon, total nitrogen and cation exchangeable capacity.
INTRODUCTION
The accepting or rejecting the null hypothesis, judgments must be made about the extent of differences between dependent variables or whether the strength of any relationship between the X and Y is sufficient to merit the rejection of the null hypothesis. These decisions are more easily made when quantitative data is available. Statistical analysis allows such judgments to be made in a quantifiable way (Wilson et al., 2001).
Statistics provides a set of objective procedures to efficiently summarize the data and compare the values. Where the results of analysis suggest that the null hypothesis should be rejected, the use of statistics allows the reliability of this decision to assess. This is the great value of statistics. In an ideal world, predictions would be tested and a hypothesis accepted or rejected with complete certainly (Austin and Noy-Meir, 1971). Unfortunately, in biology, ecology and environmental sciences it is rarely possible to formulate the null hypothesis so it can be accepted or rejected in a simple yes or no manner. A decision has to be made on the balance of probability that a decision is correct, allowing choices to be made in an informed manner (Mesdaghi, 2001).
Since variations of soil properties are high and soil characteristics relate to several variables, multivariate analysis should be more useful than other univariate statistic methods. Multivariate analysis methods assume that several random variables have interred relation between each other and that they are important for researcher (Moghadam, 2000). Classification and ordination are two main branches of multivariate analysis and in this study as ordination method were used. In general, ordination is the ordering of objects along axes according to their resemblance (Mc Cune and Mefford, 1999). The major objective is to achieve an effective data reduction, expressing many dimensions. One of the most famous ordination methods is Principal Component Analysis (PCA) (Mesdaghi, 2001, 2005).
PCA is the basic Eigen-analysis technique and was first introduced by Pearson used in ecology (Mesdaghi, 2001). It maximize the variance explained by each successive axis. Although it has severe faults with many community data sets, it is probably the best technique to use when a data set approximates multivariate normality and variables have linear relationships. This method was gradually introduced as an applied multivariate analysis method in ecology after 1966 and is one widespread method to ordinate data in ecology (Pourbabaei, 2004). Although using this method for ordination of floristic data meets some limitations (Mesdaghi, 2001), results of PCA on site data collected with floristic data can be valuable (Goldsmith, 1973).
Several studies investigating the relationship between environmental factors and vegetation have been conducted using PCA in Iran (Zahedi Amiri, 1998; Zahedi Amiri and Mohammadi Limayee, 2002; Hassanzad Navroodi et al., 2004, Salehi, 2004; Taleshi, 2004; Mahmoodi et al., 2005; Salehi and Zahedi Amiri, 2005; Salehi et al., 2005).
Also, using of principal component analysis are well studied and results for models have been reported by Schoenholtz et al. (2000), Ewald (2000), Hedman et al. (2000), Wilson et al. (2001), Stendhal et al. (2002), Phillips et al. (2002), Crowley et al. (2003), Lyon and Gross (2005), Ajbilou et al. (2006) and Corney et al. (2006). In this study we used PCA as valuable method to define the most significant and important soil properties in differentiation of forest types in lowland forest of khanikan chaloose (North of Iran).
MATERIALS AND METHODS
Study site location and description: Khanikan forests are located in the lowland and midland of Mazandaran province in North of Iran with the area of 2807 ha-1 (Between 36° 33´ 15", 36° 37` 45" latitude and between 51° 23` 45", 51° 27` 45" longitude). The maximum elevation is 1400 m and the minimum elevation is 50 m. Minimum temperature in December (7.5°C) and the highest temperature in June (24.6°C) are recorded, respectively. Mean annual precipitation of the study area were from 237.6 to 47.5 mm at the Noushahr city metrological station, which is 10 km far from the study area. This research was performed in mid-summer of 2006.
Soil sampling and analysis method: In lowland region 268.7 ha-1 of this forest was selected. For investigation of tree and shrub covers sixty quadrate (20x20 m AR.) (Hedman et al., 2000; Grant and Loneragan, 2001; Mesdaghi, 2001) and sub quadrate (1 m2 AR.) in each quadrate for investigation of herbaceous covers (Mesdaghi, 2001), were taken by randomized-systematic method. Considering variation of vegetation and environmental factors, floristic list and canopy cover percentage were determined in each quadrate. Vegetation cover data were recorded using ordinal scale of Van-der-Marel (1979). Soil samples were selected from organic horizon (litter layer) and mineral layers (0-10, 10-20 and 20-30 cm).
Soil pH, bulk density, saturation moisture content, electrical conductivity, organic carbon, total nitrogen, cation exchangeable capacity, available phosphorous, soil texture, calcium carbonate content, number and biomass of earthworms, litter carbon and litter nitrogen were determined. Slope, elevation and aspect were recorded for each plot using a clinometers, altimeter and compass, respectively.
Data vegetation for each species in all of the sampling plots and the total vegetation for all of the species within each plot were calculated. Two-way indicator species analysis method (TWINSPAN) was used for classification of vegetation groups to determine forest types (Hill et al., 1976). For above mentioned multivariate analysis PC-ORD program version 3.17 was used. Comparing of means of environmental factors amongst forest types and also study of inter-relationships between these variables was evaluated by one way ANOVA (Analysis of ANOVA) method in SAS of statistical program.
PCA is the ordination technique that constructs the theoretical variable that minimizes the total residual sum of squares after fitting straight lines to the species data. PCA does so by choosing the best values for the sites. The apply PCA; data standardization is necessary if we are analyzing variables that are measured in different units. Also, species with high variance, often the abundant ones, therefore dominate the PCA solution, whereas species with low variance, often the rate ones, have only minor influence on the solution. These may be reasons for applying the standardized PCA, in which all species receive equal weight (Jongman et al., 1995). Therefore, data was centered and standardized by standard deviation.
RESULTS
Distribution of soil properties and forest types by PCA axes: The Eigen values and proportion of variance explained by the axes are listed in Table 1. The proportion of the variance is simply the Eigen value for that axis divided by the total variance i.e., the sum of the diagonal of the cross-products matrix.
The soil properties and forest types occupied different regions of the
diagram (Fig. 1). According to the Pearson and Kendall
correlation coefficient among soil variables and between them and main
axes of PCA. The correlation structure is used to position objects in
the ordination space and objects close in the ordination space are generally
more similar than objects distant in the ordination space (Mc Cune and
Mefford, 1999).
Table 1: | Percentage of variance and Eigen values of PCA axes |
Fig. 1: | Distribution of soil factors and forest types in the two first PCA axes |
As shown in Table 1, the first axis of PCA is the most important to explain variance across the variables (soil factors). The percentage of Eigen value for the first axis and second axis of PCA are about 90 and 89%, respectively. As mentioned before, the location of soil properties and forest types in different regions of PCA axes is on the basis of correlation coefficient between variables, therefore the location of the variables in the diagram is significant and important. For example, forest types that occupied a specific region and are close to each other have several similar characteristics between each other and consequently the correlation coefficient between those characteristics should be high. On the other hand, the characteristics of profiles occupying opposite places in the diagram, for example opposite direction of an axis show distinct variations.
Eigen vectors and correlation coefficients between soil variables and the PCA axes: The values of the Eigen vector for soil variables and the axes of PCA are shown in Table 2. Eigen vectors contain set of scores that show the weight of each variable (soil variables) on each axis of PCA. Eigen vectors very between -1 to +1 and if the value of the Eigen vector for a specific variable is close to absolute of 1, it is more important to weight on the axes (Mesdaghi, 2001). As shown in Table 2, the Eigen vectors for content of clay, nitrogen, bulk density and acidity in second layer and acidity in one and third layers for the first axis and the value of clay, silt, electrical conductivity, nitrogen in one layer, phosphorus, carbon to nitrogen ratio in second layer and carbon to nitrogen ratio in third layer for second axis, are higher compared to the Eigen vectors of other variables. The last mentioned subject showed that these variables play a more important role to explain the first and second axes of variation.
These coefficients express the linear (Pearson`s r) and rank (Kendall`s
tau) relationships between the PCA scores and the individual variables
used to construct the axes. The table of correlation coefficient can be
quite helpful in providing a quick interpretation of the ordination. As
shown in Table 3, acidity of one layer (r = 0.999),
acidity of second layer (r = 0.977), acidity of third layer (r = 0.970),
bulk density of third layer (r = 0.915), nitrogen of litter (r = 0.940),
content sand of second layer (r = 0.911), lime of one layer (r = 0.920),
lime of second layer (r = 0.935) have a negative relationship to the first
Table 2: | Eigen vectors of soil variables for the first six PCA axes |
a, b, c cods = 0-10, 10-20, 20-30 layers, respectively |
Table 3: | Pearson and Kendall correlation between soil variables and PCA axes |
r = Linear (Pearson`s r) correlation coefficient tau = Rank (Kendall`s tau) correlation coefficient |
axis of PCA. Saturation moisture (r = 0.972), percentage of organic carbon (r = 0.973), content of nitrogen (r = 0.956), carbon to nitrogen ratio of litter (r = 0.910), cation exchangeable capacity of third layer (r = 0.926), content of clay of second layer (r = 0.956), content of clay of third layer (r = 0.991) are positively correlated to the first axis of PCA. In relation to the second axis of PCA, it is negatively correlated to the carbon to nitrogen ratio in third layer (r = 0.949) and positively correlated to content of clay in one layer (r = 0.909) (Table 3).
Forest types and relationships between them and soil properties: The
result of TWINSPAN method determined five main tree groups. We used these
groups as forest
Table 4: | Mean of soil chemical and physical properties for forest types in study area |
types to study relationships between them and soil properties. Relationships between some important soil properties and forest types have been displayed in Table 4.
DISCUSSION
Forests, as complex ecosystems, usually have ranges of soils with varying properties. Some of these soil properties are more important than the others and play an obvious role to distinguish the soils covering specific regions. Understanding of the segregated factors among several different soil properties is useful for forest management (Yau, 1965; Salehi and Zahedi Amiri, 2005). PC analysis is a useful data reduction technique which works by removing Interco relations among variables (components). By using PC analysis, not only is the number of comparisons between treatment means reduced, but the meaningfulness of these comparisons is enhanced. In this research, using of multivariate analysis obviously set apart differences of soil properties among soil profiles (Mesdaghi, 2001).
The results showed that in the study area, among different soil factors, the distribution of forest types was most strongly controlled with some soil characteristics such as acidity, bulk density, texture, phosphorous, organic carbon, total nitrogen and cation exchangeable capacity (Table 2). Result of principal component analysis showed, the first two principal components together accounted for 78.55% of the total variance in data set. Therefore, 59.11 and 19.43% variance were accounted for by the first and second principal components, respectively (Table 1).
The obtained result showed that the first axis has the most correlation with productively factors and the second axis has the most correlation with physical factors of soil. This result has been reported by many investigations (Zahedi Amiri and Mohammadi Limayee, 2002; Stendhal et al., 2002; Hassanzad Navroodi et al., 2004; Mahmoodi et al., 2005; Salehi et al., 2005). In humid and sub humid regions, the relation between species distribution and acidity gradient has been reported by many investigators (Oland et al., 1998; Ewald, 1999; Coker, 2000; Ewald, 2000;Phillips et al., 2002; Crowley, 2003; Taleshi, 2004).
Also, soil texture and bulk density controls distribution of plant species by affecting moisture availability, ventilation and distribution of plant roots. Soil texture is the most fundamental soil physical property controlling water, nutrient and oxygen exchange an uptake (Schoenholtz et al., 2000) and influences the growth and distribution of vegetation (Fisher and Binkley, 2000). Organic carbon and nitrogen are the effective factors in the differentiation of vegetation types (Zahedi Amiri and Mohammadi Limayee, 2002; Salehi et al., 2005). The role CEC and available phosphorous, as key elements in the distribution of plant species, is described by Zahedi Amiri and Mohammadi Limayee (2002).
On the other hand, although several environmental factors affect the rate of litter decomposition, the rate of litter fall is remarkably uniform among tree species growing under similar soil and climatic conditions (Fisher and Binkley, 2000). Finzi et al. (1998) showed that there were differences in the size and distribution of carbon and nitrogen pools at varying soil depth beneath six different tree species in temperature forests. Edwards et al. (1970) reported that the decomposition rate of beech leaves is lower than oak, elm, birch and ash leaves. Comparison of leaf decomposition rates among different plant species revealed that litter of Carpinus orientalis was decomposed more rapidly than Quercus castanifolia (Cornelissen, 1996). In the one forest type that dominant tree species was Carpinus betulus, thus, collection of climatic conditions and activity of earthworm (especially in the third layer) were effective in low and high content of carbon. Therefore, in the one forest type with Carpinus betulus dominant species, content of carbon was lower, presence and activity of earthworms were more visible, thus, decomposition of organic matter were higher (Table 4).
In this study, PCA showed noticeable variations of soil properties in the study site. Since these methods are of high accuracy and have different abilities, they could be used for habitat analysis and determination of effective ecological factors. Analyzing ecological data using ordination methods makes simpler understanding of the complex relationship between plants and environmental gradients. In addition, these methods prevent presence of ineffective factors and data complexity from affecting ecological models.