Interrelationships of Land Use / Cover Change and Topography with Soil Acidity and Salinity as Indicators of Land Degradation

As soil is the basis of all terrestrial ecosystems, degraded soil means lower fertility, reduced biodiversity and reduced human welfare. Therefore the focus of this paper is on elucidating the influence of land use and land cover (LULC) change on two important soil quality indicators that are fundamental to effective measures for ameliorating soil degradation; namely soil acidity and soil salinity in the Lower Hunter Valley of New South Wales, Australia. First, Analysis of Variance was used to elucidate the effects of LULC categories on soil acidity and salinity. The results indicate that soils under Vineyard have significantly higher pH. In contrast there is no significant effect of LULC or its change on soil salinity. To further elucidate the complex interactions of these soil quality indicators with landscape attributes over 20 years and other terrain attributes, multivariate ordination techniques (correspondence analysis and canonical correspondence analysis) were used. The results show that elevation exerted a more dominant influence on pH than the LULC types and their dynamics. In comparison, salinity of the soil appears to be higher in subsoil layers under woodland than under other LULC categories. The environmental implications of these interactions, as evidenced by this study, provide some insights for future land use planning in the region.


Introduction
Land use and land cover change (LUCC) has been recognized as a major driving force of global environmental perturbations.Through human history, land surface cover has been altered through land clearing of natural forests, subsistence agriculture, intensification of farmland production, modification of rangeland and urbanisation.The consequence of this change in the natural land cover is the modification of the world's landscape in pervasive ways [1,2], in spite of the positive outcomes of LUCC in terms of increased food and fibre production to meet the increasing needs of world's growing population.As such LUCC has brought about unpleasant consequences such as global and regional climate change, distressed global biogeochemical cycles (e.g., carbon, nitrogen, and water cycles), declining biodiversity and soil degradation.As soil is the basis of all terrestrial ecosystems, degraded soil means lower fertility, reduced biodiversity and reduced human welfare [3].Thus, protecting soils from such degradation is extremely important to providing safe and sound habitat for the future generations.Therefore the focus of this paper is on elucidating the influence of LUCC on two important soil quality indicators that are fundamental to effective measures for ameliorating soil degradation; namely soil acidity and soil salinity.
It is well known that soil acidity is caused by acidification as a natural process via mineral weathering in the presence of a continuous input of very dilute carbonic acid, followed by the leaching of the weathering products [4].However, accelerated acidification occurs when soil is exploited for agriculture, especially under intensive systems where large amounts of N are recycled through the system, and when substantial amounts of organic products and residues are removed [4,5].Acid soils are infertile through a number of mechanisms which restrict plant growth: (i) there is increase in concentration of toxic elements such as aluminium and manganese in the soil solution; (ii) as a consequence essential nutrients, such as Ca, Mg, K, are displaced and leached from the root zone; and (iii) other nutrients like P and Mo are made less available to plant roots.Such soils with high soil reaction are estimated to cover approximately 40% of the land area in Australia [5], where accelerated acidification has been reported to be more widespread than salinity, especially in the southeastern regions of Australia where the study region is located.About 18 million ha of the agriculturally intensive regions of Australia is affected by strong soil acidification (pH 4. 3-4.8), of which more than a quarter of the affected areas (5 to 7 million ha) exists in New South Wales [6].For sustainable management of acidic soils, at least in agricultural production systems and ultimately for economic and social sustainability, the panacea to acidification in farming systems is to apply a neutralizing agent such as lime.The question then is, under the agricultural production systems, how do the opposing effects of soil acidification and liming determine ultimate soil acidity?
Like sodicity, soil salinization is another important land degradation process that limits crop production.Soil salinization, which is as a result of accretion of Na + , Cl − , Mg 2+ and SO 4  2− ions in soils, is recognised as the major cause of downfall of past civilisations [7].In Australia, saline soils with the prevalence of sodium salts constitute about 30% of the land area [7].Out of this, dryland salinity is one of the most serious environmental and resource management problems in the country [8].Dryland salinity is mainly a function of rising groundwater tables, caused by increased recharge following replacement of native vegetation with annual/pasture systems [9].As such salts deposited below ground are brought to the surface causing damage to agricultural production and infrastructure, and contaminating streams.A report published by the Australian Natural Resource Atlas [10] estimated that about 5.7 million ha of Australia have a high potential for developing dryland salinity through shallow water tables, with this is expected to rise to 17 million ha by 2050.Similarly, the National Dryland Salinity Program [11] projected the cost of dryland salinity will increase by 60-70 percent over the following 20 years.
If soil acidity and salinity are to be managed effectively, their causes and behaviour, vis-à-vis management practices, must be understood.The monitoring of soil acidity and salinity through understanding of their properties as influenced by LUCC is important, as these soil quality indicators are good indicators of soil productivity.Additionally, these quality indicators are particularly prone to change in short time spans.Specifically the influence of LUCC on these quality indicators has special relevance in the Lower Hunter Valley which has experienced rapid land use change in the last 20 years.Salinity and acidity problems in the area have been reported by a number of studies [12][13][14].Salinity, in particular, has been a problem in the Hunter Valley since it was first settled in early 1800s as recorded by early explorers who noted high salinity levels in creeks around the Bayswater area [15].Human activities, such as agriculture, mining and power generation, have exacerbated this problem resulting in more than 50,000 hectares of land affected by dryland salinity in the Valley [12].
Based on the foregoing discourse, the aim of this paper is to test the two hypotheses that: LUCC has no significant influence on soil acidity/salinity in the lower Hunter valley.
The soil acidity and salinity interactions with LUCC can be differentiated from interactions with other landscape features such as landscape processes influencing soil acidity/salinity; in other words landform attributes combined with LUCC can better explain the spatial patterns of soil acidity/salinity.
To test the first hypothesis, the widely used analysis of variance (ANOVA) was performed on LULC categories in relation to soil pH and electrical conductivity (EC1:5) as measured using 1:5 soil: water suspension.The second hypothesis was tested by performing multivariate ordination analysis to unravel the complex interrelationships of the landform attributes and LULC derivatives with soil pH and EC1:5.However, before proceeding to these analyses a brief review of the statistical methods is covered.
There exist both simple and multivariate approaches that could be used to decipher the influence of LUCC on soil patterns or soil quality.Simple statistical techniques, such as regression, correlation analysis and ANOVA, are commonly used to study the effect LULC change on soil properties [16,17].For example the one-way ANOVA table could be used to summarise the response information for each level of the nominal or ordinal factor-which could be a land use/land cover category.One-way ANOVA can also be used for multiple comparisons to test differences among group means, whereby student's t test, are required for pair-wise comparisons.
When confronted by more complex situations in which there are two or more target response variables and a set of explanatory or predictor variables, methods such as simple/multiple regression analysis or one-way ANOVA, or even MANOVA, may not be robust in elucidation the interrelationships.These methods cannot provide a broad overview of how soil properties behave in relation to the combined effects of the predictor variables.For example while the transient soil quality indicators, such as acidity and salinity, are greatly influenced by the terrain attributes, namely: digital elevation model (DEM), slope, compound topographic index, etc. [18][19][20], it is important to deconstruct the influence of land use/cover on these indicators from the complex interactions with the terrain attributes.It could be more informative and economical to handle the two sets of data all simultaneously using direct or indirect gradient techniques [21], by including the terrain attributes in the analysis.In such cases, multivariate ordination gradient techniques have been recommended as more appropriate for exploring the complex interrelationships between a set of explanatory attributes on one hand and another set of target variables [22][23][24][25][26][27][28].Moreover, LULC attributes or factors are qualitative in nature, and as such are not amenable to or meets the assumptions of the classical ANOVA or regression analysis.ANOVA requires a number of assumptions to hold true such as independence of observation, linearity, and additivity of effects, homogeneity of variance, normally distributed errors [29,30].As these assumptions are not met in many situations, robust multivariate techniques such as CA (correspondence analysis) and CCA (canonical correspondence analysis) have been recommended as they do not require meeting these stringent assumptions.
CA and CCA techniques are powerful tools for investigating and summarizing underlying trends in complex data structures [31].The techniques have been widely applied in species environment relationship [22,23,28].CA, as one of the popular indirect gradient ordination techniques in ecology [22,31], uses the efficient reciprocal weighted averaging algorithm.The technique is an indirect gradient ordination technique because it is followed by environmental gradient identification [22].In contrast CCA as a direct gradient technique is a canonical form of CA in which the site scores are restricted to be linear combination of measured environmental variables.Although, these methods have potential to be applied in soil-environment studies in order to understand a complex relationship, only a few studies [18,32,33] have applied these techniques to soil studies.No single study to our knowledge has included land use/land cover factor in an environmental gradient analysis.This paper is aimed at exploring this possibility.

Land Use/Land Cover Data
In this paper the historical LULC images, covering the study region for 1985, 1991, 1995, 2000 and 2005, were used to decipher the patterns of LU/LC change and to explore the effect of LUCC on soil acidity/salinity.Detailed description on how LULC images were obtained can be found in Manandhar et al. [34].

Soil Quality Indicators Data
Soil samples were collected using a Latin hypercube sampling scheme [35,36].Thus the scheme used is based on a stratified random procedure which provides an efficient way of sampling from the multivariate distribution of the covariates, as it covers a full range of each covariate by maximally stratifying the marginal distribution.Based on this scheme 246 soil sample locations were selected (Figure 1).At each location soil samples were collected at predetermined depths: 0-10, 10-20, 20-30, 40-50 and 60-100 cm.After oven drying at 105 °C, each sample was ground and fine earth separated through a 2-mm diameter sieve.The electrical conductivity and pH of each sample were measured in a 1:5 soil: water suspension after 30 min of shaking, the pH and Electrical Conductivity (EC) measurements were done using pH and EC meters respectively [36].

Derivation of Topographic Attributes
The topographic attributes that would potentially have influence on the two soil quality indicator values include digital elevation model (DEM) and its derivatives-slope, upslope area (A i ) and compound topographic index (CTI).These attributes have been observed to influence a number of soil properties [18,37,38].Large CTI values, usually found in lower part of a watershed and convergent hollows associated with gentle slopes and/or soils that are characterised by low hydraulic conductivity, are indicative of increased likelihood of saturated conditions [36].The DEM of the study region was obtained from the then New South Wales Department of Natural Resources (DNR).As a generic attribute, the DEM was used to derive the following topographic attributes which are deemed influential on soil acidification and salinization processes, using SAGA GIS tools.
Slope, as an important attribute controlling the flux of materials within the landscape, is an important controlling feature in hydrological and geomorphological processes.It can be defined as the rate of change of elevation down the steepest descent across a point in the landscape [39].Slope can be derived in either degree or as gradient (decimal or in per cent).For the purpose of this study, slope was derived in degree from the DEM, using the slope function within the SAGA GIS environment.
Specific contributing area: Upslope contributing area (A) can be computed as the total area above a point or a certain length of contour that contribute flow to the point or across the contour [39].More appropriate for this study is the specific contributing area, A s , which is the ratio of the contributing area New South Wales and the contour length, A/l.In calculating A s from a gridded DEM, the contour length was assumed to be approximately the size of a grid cell; in the simplest case the contributing area was determined by the number of cells contributing flow to a cell [39].Upslope contributing area is an important landform attribute that measures or approximates the quantity of water that passes through a point or across a contour length.
Compound topographic wetness index (CTI) was derived using the equation: where A s is the specific upslope contributing area as described above and β is slope in degrees.CTI represents the topographic control on soil wetness, based on a number of assumptions listed elsewhere [39].

Derivation of LULC and Terrain Attributes for Soil Sampling Points
Ground co-ordinate points from where soil samples were collected (Figure 1), were overlain on LULC maps using ArcGIS (ESRI, 2009) from which the historical LULC information for each point was obtained.The derived terrain attributes (elevation, slope and CTI) were also overlain on the LULC maps and their values were then obtained for each soil sampling points.

Assessing the Effect of Land Use/Land Cover Change on Soil pH and Salinity
In the first instance, the effects of LUCC on soil acidity and salinity were each tested using ANOVA.The latter is used here as an attribution and tests the part of the total variability in a response which is due to the differences in mean responses among the LULC categories.The ANOVA was performed to test if LULC indices, shown in Table 1, had significant impact on soil acidity/salinity.While untransformed pH values were used for the analysis, EC values were log-transformed prior to ANOVA in order to meet the requirement of normality of distribution.Pair-wise comparison was then performed using student's t test to determine if the mean value of a given category is significantly different from the mean of every other category.

Multivariate Analysis of Interrelationships of the Soil Quality Indicators with Landform Attributes and Land Use Covariates
Out of the various multivariate ordination techniques used in environmental studies, principal component analysis (PCA), correspondence analysis (CA), redundancy analysis (RDA) and canonical correspondence analysis (CCA) were initially tested.Based on percentage variance accounted for in the target variable-environment interactions by the resulting ordination axes, CA and CCA were chosen for further analysis as they have variance distributed more or less uniformly across the first few axes.CA, as a popular ordination technique in ecology [22,31], uses an efficient reciprocal weighed averaging algorithm based on the non-monotonic Gaussian species response models.Thus, CA is an indirect gradient analysis technique while CCA is a canonical form of CA in which the site scores are restricted to be linear combinations of the environmental variables [18,22].In ecology, CA extracts continuous axes of variation from species occurrences or abundance data, typically interpreted with the help of external knowledge and data on environmental variables; this two-step approach (ordination followed by environmental gradient identification), is termed as indirect gradient analysis [22].The eigenvalues in CCA are usually smaller than those in CA because of the restrictions imposed on the site scores in CCA.In order to incorporate the influence of LUCC on the soil quality indicators over the period between 1985 and 2005, several types of new LULC dynamic indicators (Table 1) were derived from the series of LULC maps.For example, Vineyard consistency (Vine_C) = area that remains as vineyard throughout the accounting period (1991-2005); Pasture/scrubland consistency (Past_C) = area that remains as pasture/scrubland throughout the accounting period (1991-2005); F95vine = area that was converted from other LULC types to vineyard between 1991 and 1995 and Aforest = Area converted to woodland between 1991 and 1995 or 2005.These indices were derived as nominal data for every sample location, with 1 representing true and zero not.
The CA and CCA were carried out using a FORTRAN program (CANOCO version 4.5), written by [40].The selected soil properties (pH and EC) are the target variables and LULC categories, along with the indices derived above, and the landform attributes (DEM, Slope, and CTI) were the environmental (explanatory) variables used.The list of LULC and other landform attributes and soil parameters used for derivation CA and CCA are shown in Table 1.The transition indices (conversion of one LULC to another during the accounting period, also shown in Table 1) were not used for CA and CCA due to their high correlation with the derived consistency indices and the resultant inflation factors.Biplots of both CA and CCA were obtained to display the complex interrelationships between the soil quality indicators on one hand, and topographic attributes and LULC types on the other.In the biplots the predictor (LULC and topographic) variables are represented by arrows and soil variables are represented by points.The ordination diagram, as it is sometime referred to, shows the main pattern of variation in soil variables as accounted for by the predictor variables.Thus the complex relationships among the selected soil properties (soil pH and EC) and the predictor variables (LULC information, DEM, slope, CTI) can be elucidated.Topographic attributes-DEM, slope and CTI-were chosen because they are the natural attributes that affect most of the landscape processes such as filtration, runoff, soil reaction and salinization.

Effects of LULC on Soil Acidity
The ANOVA of pH in relation to LULC systems indicates a significant difference (P = 0.05) in pH values.Vineyard consistency (land area continuously under Vineyard since 1985) and F91vine (area converted to Vineyard between 1985 and 1991) are found to have higher pH value, significantly higher than that under Woodland consistency and Pasture/scrubland consistency (Table 2).This is specifically true for the topsoil (0-10 cm). Figure 2 illustrates the diamond plots of the pH versus different LULC consistencies for the different soil layers.The line across each diamond represents the group mean.The vertical span of each diamond represents the 95% confidence interval for each group, and the horizontal extent of each group along the x-axis of diamond is proportional to the sample size for each category of the LULC [41].The diamonds of F95vine (i.e., area converted to Vineyard between 1991 and 1995) and F05vine (area converted to Vineyard between 2000 and 2005) are horizontally narrow but vertically long, indicating they have very few observations with high variability.However, the diamond of Past_C (consistent Pasture/scrubland over the study period) has the largest number of observations as it constitutes the largest component of the landscape, and by chance has more soil samples.Wood_C (consistent Woodland) has similar pattern to the Past_C, as this constitutes the second largest component of the landscape.The circles on the right of the plots are indicative of pair-wise comparison means of each soil property under different LULC types.Circles for means that are significantly different either do not intersect or intersect slightly so that outside angle of intersection is less than 90 degrees (Figure 2a).If the circles intersect by an angle more than 90 degrees or they are nested, the means are said to be "not significantly different" Thus, the highlighted Vine_C (thick red circle) is not significantly different from other red circles but significantly different from those other LUCC types or transition represented by dark grey circles indicating higher pH for Vineyard soils, compared to mean pH under these LULC systems.Overall soil layers under Wood_C (Woodland consistency) had the lowest mean pH value followed by Aforest (other LULC categories converting to Woodland), and then by Wood2past (Area converted from Woodland to Pasture/scrubland) and Past_C (Figure 2).While the pH of land under the area that was converted to Vineyard was generally characterised by the mean pH values; i.e., lower than that under the Vine_C but higher than that under the other LULC categories (i.e., Woodland, Aforest and Wood2Past and Past_C).
The depth profiles of the mean pH under the three major LULC categories (Figure 3) indicate Vine_C is characterised by the highest pH values at all the depths followed by Past_C and then by Wood_C.This result is corroborated by the finding of Minasny and McBratney [42] who found that lands under viticulture have significantly larger than under non-viticulture lands, by about 0.5 pH unit.The higher pH of under Vineyard could be due to either or both of the two following pathways: (i) Neutralizing effect of Vineyard land use-this may be due to application of lime under vineyard cultivation [42].(ii) Landforms which attracted cultivation of Vineyard originally had higher soil pH.
As reported by Faine and others [14] the intensive grape-producing area is located where the acid buffering capacity is considered to be moderate, requiring 5.5 to 11 M of H+ to decline to 4.5.Therefore the second pathway is probably the dominant factor.

Effect of LULC on Soil Salinity
Salinity, as measured by EC, does not vary significantly under different LULC categories (Table 3).It is also interesting to note that land management practices under vineyard did not have significant impact on soil salinity level.However, Woodland consistency was found to have lower topsoil EC compared with the same values under other LULCs whereas in the case of the subsoil, Woodland has slightly higher EC compared to others.This is probably due to increased leaching of salts from the topsoil to the lower soil layers under forest, and/or due to preponderance of soil organic matter and hence higher microbial activity under forest soils compared to cultivated soils [16,43].In the case of the layer just below the surface layer, there was no definite trend in EC across different LULC categories (Table 3).Consistent with the above findings, the depth profiles of EC1:5 shown in Figure 4 indicate that the topsoil salinity is generally low.However, there is a distinctive increase of salt level for all LULC categories in the deeper soil layer (>50 cm).This may be due to provenance of salt from the parent material, and perhaps more importantly due to the leaching of salts from the upper soil layers caused by percolating water from rainfall and/or irrigation.The foregoing discussion (i.e., higher level of salinity of subsoil in comparison to upper layers) is corroborated by the work of Taylor [38] who reported that subsoil EC influenced Vineyard performance in Pokolbin wine-producing region, based on comparison of EC and NDVI maps of the region.The high amount of salt in the soil solution probably impacted on viticultural production systems (overall concentration of soluble salts) or sodicity (the relative concentration of sodium to magnesium and calcium).Similarly, Fitzpatrick et al. [44] and Cass et al. [45] reported that salt-affected soils, especially those associated with sodic soils are a major threat to many viticultural/agriculture production systems in southern Australia.Faine et al. [14] also reported occurrence of moderate to high salinity in a large area in the northwest region of the study area, the soils of which were found to be predominantly sodic.However, in contrast to all these findings, this study shows soil EC as below the hazardous level of salinity (4 dS•m −1 or 4000 µS•cm −1 ) even at the deeper sub-soil depths (>50 cm), which in our study have higher EC than the upper soil layers.Thus our findings here indicate that soil salinity is not a problem under cultivation and/or management practices.In fact, the local farmers group called Polkobin Irrigation District has an ongoing salt monitoring system for irrigation water; it has installed an EC measuring device at many of the river pumping stations, which are programmed to stop pumping river water when the EC reaches a certain critical level [14].This might have prevented or delay the rise of soil EC.

Interrelationships of Soil Quality Indicators with Environmental Covariates
Ecosystems are complex consisting of many interacting biotic and abiotic components.Multivariate ordination techniques such as Principal Component Analysis (PCA), CA and CCA have advantage over ANOVA as they enable the inclusion of numerous target variables in the analysis and, thus highlight their associations with the predictor variables.A useful outcome of these multivariate techniques is the biplot, which is an extension of the ordination diagram, in which the properties of the columns of the matrix (the variables analysed) are overlain on top of the main ordination diagram in order to highlight the relationships between the two sets of data (the target attributes and the predictor variables).The biplots rapidly allow the user to explore which target attributes and/or predictor variables are strongly associated since proximity implies close association.
In order to understand the ordination diagram, each arrow in biplots represent a predictor variable which determines the direction relative to the two axes in the diagram (Figure 5); soil attributes (pH and EC), as represented by points, are projected onto each of the axes.The lengths of the vectors and angles between them are more important rather than their inter-point distances [18,22,46].The predictor variables with long arrows, closer to a given ordination axis, are more strongly correlated with the ordination axes than those with short arrows; thus those represented by long arrows are more closely related to the pattern of landscape variation shown in the ordination diagram.Thus the ordination diagram characteristically shows the main pattern of variation in soil variables as accounted for by the predictor variables, in this case, LULC categories and systems and selected landform attributes.
As unimodal models have been considered as more general than monotonic [18], this study opted for analysing the data using CA and its canonical counterpart, CCA.Both CA and CCA biplots convey

Depth of soil in cm
EC of Soil similar overall information about the complex interrelationships, although the effect of the predictor variables is less prominent in the CA biplot in comparison to that shown by the CCA biplot (Figure 5).The biplot of CCA (Figure 5b) demonstrates that Vine_C (Vineyard consistency) is more closely associated with slope; the arrows (predictor variables) pointing roughly in the same direction, indicating high positive correlation, while arrows pointing in the opposite directions imply high negative correlation and that arrows crossing at right angle indicate zero correlation [21,22,32].This is expected, as lands characterised by sloppy terrain are preferred for vineyards because of their advantage of better drainage and a reduced risk of water logging [38].Unlike the results of ANOVA, the CCA biplots shows no evidence of a close association of pH with Vine_C, rather elevation is evidently more closely association with pH than Vine_C.Thus, the higher the elevation the larger is the pH level, i.e., the more alkaline the soils become, especially in the topsoil layers.In contrast, Wood_C is negatively associated with pH, indicating soil under Woodland to be more acidic compared to soils under other LULC categories (Pasture/scrubland and Vineyard).The influence of LULC types and landform attributes on salinity is also well illustrated by the CCA biplot (Figure 5b), which shows EC values at different layers spread out in the biplot.This is indicative of how EC at different layers were variably associated with the different LULC types.As such, Wood_C has negative association with topsoil EC (i.e., EC10) but is somewhat more closely associated with EC50 and EC30.The large values of soil EC at lower depth in contrast to the lesser values of EC and lower values of pH in the topsoil under Woodland LULC type may be attributable to microclimate of forest flora which may have increased organic matter and soil microflora especially at surface layer [43,47].As demonstrated by the biplot of CCA (Figure 5b) CTI and Past_C have somewhat closer association with ECdeep which is logical in the sense that soluble salts would have leached to the lower depths.Additionally, the larger values of CTI are usually found in the lower parts of watersheds and convergent hollow areas associated with soils with low hydraulic conductivity or areas of low slope [42].Generally this study found a few discrepancies between CA and CCA: In CCA, the first ordination axis is largely due to elevation gradient and the second axis is largely accounted for by CTI, which mainly explains the variation in Pasture/scrubland, while in CA, the first ordination axis is mainly accounted for by Vineyard followed by slope and second ordination axis is largely explained by Woodland.
As the length of the arrow quantifies the rate of change in that direction, these two methods indicate differences in variation in the included variables that differ along the environmental gradient.Again, CCA is constrained to environmental gradient rather than the real value.
Table 4 shows the eigenvalues and variance accounted for by the first four axes of CA and CCA.The variance accounted for by the first two axes of CA (i.e., represented by the biplots) increased from 75.5% in CA to 82.2% in case of CCA.The soil-environment correlation also improved slightly from 0.514 for CA to 0.545 for CCA.Ordination provides information as to whether important environmental variables have been overlooked: an important variable has been definitely been missed if there is no relation between the mutual positions of the sites in the ordination diagram and the measured environmental variables [23].
In our study, the significantly higher pH of Vineyard areas, as observed from the ANOVA, can be due to cultivation of Vineyards in the upper slopes which has close association with higher pH.
Vineyard cultivation is preferred in areas of slightly steep slopes and a higher elevation rather than the lower part of watershed in order to have katabatic drainage of air, in addition to the requirement of good water drainage and reduced risk of water logging [38].Therefore, higher pH of Vineyard soils must have resulted due to cultivation of Vineyard in mid-to upper-slopes apart from lime application under Vineyard paddocks.This outcome was missed by ANOVA alone, as the latter technique could not examine the complex relationship under multiple sets of factors.

Conclusions
A key finding of this study is that soils under Vineyard were found to have significantly higher pH and hence increased alkalinity, which may be attributable to land use practices under Vineyard cultivation, and/or also due to inherent soil properties or landform that attracted Vineyard cultivation.However, there is no significant effect of land use and land cover (LULC) on soil salinity, indicating Vineyard soils were not as degraded in terms of soil salinity as have been published [11].Using a robust multivariate ordination technique of canonical correspondence analysis (CCA), the complex interrelationships revealed that elevation had the dominant influence on pH in comparison to Vineyard cultivation, while Woodland had the closer positive association with subsoil salinity, but negative association with topsoil salinity and acidity.Thus, ordination analysis in addition to analysis of variance (ANOVA) opens up our mind in understanding the complex but minute variation in acidity and salinity caused by LULC change in relation to the effect of other landform attributes.

Figure 1 .
Figure 1.Location of the study area in New South Wales, Australia, showing soil sample sites as red dots.Note: the thin grey lines are the suburb boundaries.
are done across the rows.Levels not connected by same superscript letter are significantly different.Also note: Wood_C = Woodland consistency, Past_C = Pasture/scrubland consistency, Vine_C = Vineyard Consistency, F91vine = Area converted to Vineyard between 1985 and 1991, F95vine = Area converted to Vineyard between 1991 and 1995, F2000vine = Area converted to Vineyard between 1995 and 2000, and F05vine = Area converted to Vineyard between 2000 and 2005.Aforest = Area converted to Woodland, Wood2past = Area converted from Woodland to Pasture/scrubland.Others = other than previously stated.pH10 = pH of soil at a depth of 0-10 cm, pH20 = pH of soil at a depth of 10-20 cm, pH30 = pH of soil at a depth of 20-30 cm, pH50 = pH of soil at a depth of 40-50 cm and pHdeep = pH of soil at a depth > 50 cm.

Figure 2 .
Figure 2. Diamond plots of analysis of variance of soil pH at different depths.Note that the circles on the right indicate pairwise comparison; the highlighted Vine_C (thick red circle) is not significantly different from other red circles, but significantly different from thick gray circles.See Table 1 for explanation of abbreviations.(a) Soil depth (0-10 cm); (b) soil depth (10-20 cm); (c) soil depth (20-30 cm); (d) soil depth (40-50 cm).

− 1 )
of soils under different land use and land cover (LULC) categories; values in brackets are the standard errors.
are done across the rows based on log-transformed data.Levels not connected by same superscript letter are significantly different.Also note: Wood_C = Woodland consistency, Past_C = Pasture/scrubland consistency, Vine_C = Vineyard Consistency, F91vine = converted Vineyard between 1985 and 1991, F95vine = converted Vineyard between 1991 and 1995, F2000vine = converted Vineyard between 1995 and 2000, and F05vine = converted Vineyard between 2000 and 2005.Aforest = converted Woodland, Wood2past = Converted Pasture/scrubland from Woodland.Others = other than previously stated.EC10 = EC of soil at a depth of 0-10 cm, EC20 = EC of soil at a depth of 10-20 cm, EC30 = EC of soil at a depth of 20-30 cm, EC50 = EC of soil at a depth of 40-50 cm and ECdeep = EC of soil at a depth > 50 cm.

Figure 5 .
Figure 5. Biplots from (a) Correspondence analysis and (b) Canonical correspondence analysis.Note that: Wood_C = Woodland consistency, Past_C = Pasture/scrubland consistency, Vine_C = Vineyard Consistency, Elev = Elevation, CTI = Compound topographic index, pH10 = pH of soil at a depth of 0-10 cm, pH20 = pH of soil at a depth of 10-20 cm, pH30 = pH of soil at a depth of 20-30 cm, pH50 = pH of soil at a depth of 40-50 cm and pHdeep = pH of soil at a depth > 50 cm; EC10 = pH of soil at a depth of 0-10 cm, EC20 = pH of soil at a depth of 10-20 cm, EC30 = pH of soil at a depth of 20-30 cm, EC50 = EC of soil at a depth of 40-50 cm and pHdeep = EC of soil at a depth > 50 cm.

Table 1 .
Description of land use and land cover (LULC) abbreviations.
Aforest Area converted to Woodland, i.e., Pasture/scrubland in earlier maps but Woodland in the final map Wood2 past Area converted from Woodland to Pasture/scrubland, i.e., Woodland in earlier maps but Pasture/scrubland in the final map Others Other than above, such as Vineyard converting to Pasture/scrubland, or Built-up area etc.

Table 2 .
Mean pH values of soils under different land use and land cover (LULC) categories; values in brackets are the standard errors.

Table 4 .
Variance accounted for by the first four axes of correspondence analysis (CA) and canonical correspondence analysis (CCA).(Sum total eigen value for CA is 0.4 and for CCA is 0.021).