JEQ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (15)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
Related Collections
Right arrow Soil Conservation
Right arrow Soil Fertility and Productivity
Right arrow Sustainable Agriculture
Right arrow Best Management Practices
Right arrow Soil Analysis
Journal of Environmental Quality 31:1848-1857 (2002)
© 2002 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America

TECHNICAL REPORTS
Ecological Risk Assessment

Soil Quality at a National Scale in New Zealand

G. P. Sparling* and L. A. Schipper

Landcare Research, Private Bag 3127, Hamilton, New Zealand

* Corresponding author (SparlingG{at}LandcareResearch.co.nz)

Received for publication May 22, 2001.

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
New Zealand is a signatory to international conventions on environmental performance, and soil quality information is needed for reporting both at a national and regional level. Soil quality was measured at 222 sites in five regions of New Zealand (12 soil orders and 9 land-use categories). Topsoil (0–100 mm) properties measured were total carbon and nitrogen, potentially mineralizable N, pH, Olsen P, cation exchange capacity, bulk density, total porosity, macroporosity, and total available and readily available water. Our objectives were to gauge the representativeness of the sample, determine the contribution from land use or soil order to variability, rationalize the data set, and identify concerns for long-term sustainable land use. Soil and land use combinations were both under- or overrepresented in the data set compared with national distribution. Soil order and land-use categories explained 55 to 76% of the variance in soil properties. Total C contents of pastures were comparable with indigenous forest soils, but pastures were less acidic and with higher N and P contents. Plantation forests had characteristics similar to indigenous forests on comparable soils. Cropland soils comprised <1% of the national land cover and generally had high inorganic fertility and low organic matter, with evidence of compaction. Seven characteristics (total C, total N, mineralizable N, pH, Olsen P, bulk density, and macroporosity) explained 87% of the total variability. The findings are being used by monitoring agencies to raise awareness about soil quality in the wider community, set land management guidelines, and develop policies.

Abbreviations: CEC, cation exchange capacity


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
HUMAN ACTIVITIES can cause accelerated changes in dynamic soil characteristics (Carter et al., 1997; Soane, 1990). In New Zealand, human-induced changes to soil properties have been occurring for about 800 yr, with evidence of forest clearance and erosion following the first human settlement when Maori arrived from Polynesia (Molloy, 1988). Since around 1850, human influences on soils have accelerated following the arrival of European migrants, with mechanization, extensive clearance of native forests, wide-scale introduction of exotic species, and agrochemical use (Taylor and Pohlen, 1968; New Zealand Soil Bureau, 1968). Within the last 15 yr, further change in land-use pattern has occurred with greater intensification of cropland and dairy farming on better-class land. This contrasts with decreased intensity of use on hill country by conversion of pastures to plantation forest, reversion to scrub, and reestablishment of indigenous vegetation. New Zealand is a signatory to international conventions on environmental performance, and soil quality information is needed for reporting both at a national and regional level (e.g., Organisation for Economic Co-Operation and Development, 1996; Ministry for the Environment, 1997).

The principle legislative tool for environmental regulation in New Zealand is the Resource Management Act (New Zealand Government, 1991), which is administered by 17 autonomous Regional Authorities. The Resource Management Act requires effects-based criteria to determine whether a particular activity has a detrimental effect on the environment, but does not specify properties to be monitored or interpretative framework (Ministry for the Environment, 1997). To provide incentive for a unified approach to soil quality monitoring, the New Zealand Ministry for the Environment offered Regional Authorities a 60% subsidy toward monitoring costs during 1998–2001. The subsidy required that the participants use common methods and interpretation. Despite this financial incentive, only five of the 17 Regional Authorities participated in 1998–2000, with a commitment from a further five in 2000–2001. This paper reports on the data collected in the first two years of the project, which is popularly known as the 500 Soils Project. The data were supplemented with 56 sites from three other studies using comparable sampling and analytical methods (Schipper and Sparling, 2000; Sparling et al., 2000a, b).

Our data were collected primarily for state of the environment reporting at a regional scale, rather than for production purposes at a farm scale. Thus, our project differs from many in the USA, where much soil quality monitoring and reporting has a production focus (Karlen et al., 1999, 2001). Several large-scale studies in North America included a high proportion of arable farms and intensive land uses. Examples are Wander and Bollero (1999), who compared soil quality of 36 arable farms in Illinois, and Boehm and Anderson (1997), who measured soil quality on 20 farms under three crop management systems in Saskatchewan. Hellkamp et al. (2000) surveyed 293 sites in six U.S. Mid-Atlantic states with 12% of the sites on cropped land. In another large-scale study in the U.S. central and southern High Plains by Brejda et al. (2000b), 49% of the land supported intensive livestock production, and 31% was under some form of arable cropland.

Our survey is on a similar scale to some of the North American studies, but the 222 sites covered a wider range of soil orders and land uses, including arable, horticulture, pastoral, plantation forestry, and indigenous ecosystems. All sites other than indigenous ecosystems were on commercial farms or plantations. Twelve soil chemical and physical properties previously tested for variability and response to soil management (Schipper and Sparling, 2000) were used to characterize aspects of soil quality response to management.

Our objectives for this paper were to compare the soil orders and land-use frequency in the 500 Soils data set with the national distribution of soil orders and land uses, determine the contribution from land use or soil order to the variability in soil properties, provide justification for reducing the number of measured soil properties, and identify any items of concern for long-term sustainable land use.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The methods used in this study were essentially those described by Schipper and Sparling (2000). Only a brief summary and changes to the previous methodology are presented here.

Soil and Site Selection
Soil and land-use combinations were selected by the Regional Councils based on a combination of the area extent of the soils and land uses within their regions, and we used "targeted" sampling of soil and land use combinations perceived to be a risk to the environment. Soils were classified by soil order, using the New Zealand soil classification system (Hewitt, 1993). Where possible, matched sites on the same soil order but under different land uses were selected to produce a matrix of soil types and land uses (Schipper and Sparling, 2000). Soils were collected from 222 locations mainly from the North Island. The soils included examples from 13 different soil orders (Hewitt, 1993). Land use was divided into nine categories including indigenous forest; plantation forest (mainly Monterey pine, Pinus radiata D. Don); native tussock grassland used for sheep and beef farming; introduced European pasture species for drystock beef, sheep, or deer; intensively managed ryegrass–white clover pasture for dairy cows; orchards–horticulture; mixed arable cropland; and long-term use for intensive vegetable production. Not all land uses occurred on all soil orders (Table 1).


View this table:
[in this window]
[in a new window]
 
Table 1. Soil orders and land use matrix showing distribution of the 222 sites in New Zealand that are included in the 500 Soils data set.

 
Site information included location, map reference, soil series and soil classification, current and previous land use, vegetation, slope, elevation, landform, annual precipitation, parent material, and soil drainage class (Milne et al., 1995). Only soil order and land-use information are reported here. The representativeness of soil orders and land uses in the 500 Soils data set was estimated by comparing the frequency of sampling against the mapped area of each soil order and land cover from the National Soils database and Land Cover databases held by Landcare Research (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Relative proportions (percent by area) of soil orders and land cover category on a national basis (estimated from mapping and aerial photography) and the frequency (percent number of sites in that land use or soil order) within the 500 Soils data set.

 
A 50-m transect was laid out at each site. Soil cores 25 mm in diameter to a depth of 100 mm for chemical and biochemical analyses were taken every 2 m along the transect. The 25 individual cores were bulked, mixed, and sieved to <6 mm while moist. Subsamples were used to determine potentially mineralizable N. Air-dry, sieved (<2 mm) subsamples were used for the other chemical analyses. For physical analyses three undisturbed soil samples were obtained from each plot at 15-, 25-, and 35-m positions along the transect by pressing steel liners 75 mm in depth by 100 mm in diameter into the top 75 mm of soil. The liner and soil cores were removed as a unit by careful excavation around the liner, then bagged, loaded into padded crates, and transported to the Landcare Research Laboratory in Hamilton. Subsamples of the cores were then taken for physical analyses.

We took core samples to a depth of 100 mm, but subsurface and subsoil condition below this depth can affect the overall soil quality characteristics of the site. Hence, at the midpoint of each transect a pit was dug and the soil profile fully described down to at least a 1-m depth (Milne et al., 1995). The full depth of A horizon, total potential rooting depth, and the nature and depth of any limiting horizons were specifically recorded. The profile descriptions are not reported here, nor the site characterizations, but we consider that both are necessary to combine with the topsoil chemical and physical data for meaningful soil quality assessment.

Soil Quality Data Set
The number of soil properties measured for soil quality monitoring was smaller than used in our earlier work (Schipper and Sparling, 2000). Chemical properties were total C, total N, cation exchange capacity (exchangeable Ca2+, Mg2+, K+, and Na+), base saturation, Olsen phosphate, and soil pH. Soil physical properties were water retention (at {psi}s = -5, -10, -100, and -1500 kPa), dry bulk density, and particle density, which allowed calculation of readily and total available water, macroporosity, and total porosity. The only biochemical property measured was anaerobically mineralizable N. Properties previously measured but not included in the present data set were unsaturated hydraulic conductivity, soil respiration, microbial biomass C, and particle size distribution. The rationale for discarding some soil properties from the soil quality measurements was given by Schipper and Sparling (2000). In brief, unsaturated hydraulic conductivity was dropped because high variability meant impractical levels of replication were required to detect significant changes. Particle-size distribution was dropped because, although this is a useful property to identify matched sites, it is not a dynamic soil property responsive to land use. The biological measures of soil microbial biomass and soil respiration were not included because we did not feel able to interpret the data (Carter et al., 1999), and because in many cases there is a reasonable correlation with anaerobically mineralizable N (Hart et al., 1986).

Analyses
Total C and N were determined by dry combustion. For New Zealand soils total C is an acceptable measure of total organic C, as the contribution from carbonates to C content is normally negligible (Metson et al., 1979). The cation exchange capacity was determined by flame spectrometer after leaching of <2 mm air-dry soil with 1 M CH3COONH4 at pH 7.0. Soil pH and concentrations of Ca, Mg, K, Na, NH4, and PO4 were determined with standard methods (Blakemore et al., 1987). Potentially mineralizable N was estimated by the anaerobic (waterlogged) incubation method with a 7-d incubation at 40°C, extraction in 2 M KCl, and measurement by Technicon (Tarrytown, NY) AutoAnalyzer.

Water release was calculated from drainage on pressure plates (Gradwell and Birrell, 1979). Dry bulk density ({rho}b) was measured on a subsample core dried at 105°C (Gradwell and Birrell, 1979), and the remaining soil was analyzed for particle density ({rho}p) as described by Claydon (1989). Readily available water and total available water were measured as the difference in volumetric water ({theta}) between -10 kPa and -100 kPa {psi}s, and -10 kPa and -1500 kPa {psi}s, respectively. The total porosity (St) was calculated from the formula St = 100[1 - ({rho}b/{rho}p)], where {rho}p is the particle density, and macroporosity was calculated from St - {theta} at -5 kPa {psi}s (Klute, 1986).

Statistics
To allow valid land-use comparisons, gravimetric data were recalculated to a soil volume basis using the dry bulk density (Reganold and Palmer, 1995; Doran and Parkin, 1994). The length of each box shows the range within which the central 50% of the values fall, with the box hinge (borders) at the first and third quantiles; the midpoint line shows the median (Fig. 13) . The whiskers show the range of values that fall within the inner fences. Values between the inner and outer fences are plotted with asterisks. Limits for the inner upper or lower fences were set as hinge minus (1.5 times the upper or lower hinge median, respectively). Those for the upper inner or lower fence were set as the hinge plus (1.5 times the upper or lower hinge median, respectively). Values outside the outer fences are plotted with empty circles and used a similar calculation, but were 3 times rather than 1.5 times the upper or lower hinge median.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 1. Box plots showing medians and quartiles for soil organic resources (total C, total N, C to N ratio, and potentially mineralizable N), arranged by land use category.

 


View larger version (43K):
[in this window]
[in a new window]
 
Fig. 3. Box plots showing medians and quartiles for soil physical properties (bulk density, macroporosity, and readily and total available water), arranged by land use category.

 


View larger version (45K):
[in this window]
[in a new window]
 
Fig. 2. Box plots showing medians and quartiles for soil chemical properties (pH, Olsen P, cation exchange capacity, and base saturation), arranged by land use category.

 
Separate analyses of variance (ANOVA) were conducted to determine the proportion of variation in soil properties attributable to land use, soil order, or the combination of both. Significant differences between selected items of interest were determined using one-way ANOVA followed by Bonferoni post-hoc pairwise comparisons at P < 0.05. Principal components analysis using an Oblimin rotation (SYSTAT 7.0; Systat, 1992) was applied to 12 variables in the data set to identify those properties that best explained the variation between sites.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Representativeness of the 500 Soils Data
Compared with the national distribution of soil orders and land uses, the frequency of sampling in the current 500 Soils dataset corresponded to the geographical distribution of the participating Regional Authorities and was biased to those soil and land-use combinations of concern (Table 2). The most common soil orders in terms of land area (Brown, Podzol, and Pallic) were underrepresented, whereas others (Granular, Melanic, and Organic) were overrepresented. Three soil orders were not represented at all (e.g., Semiarid), consequently, the 12 orders in the 500 Soils data set shown in Table 1 are less than the possible maximum of 15 (Table 2). The frequency of the Allophanic (22%), Recent (17%), and Pumice (18%) soil orders in the 500 Soils data set reflected the predominance of these soil types in the participating North Island regions. The underrepresentation of Brown soils reflected the nonparticipation of Regional Councils in areas where these soils are common (Molloy, 1988; Hewitt, 1993). Podzols were underrepresented perhaps because these soils occur in higher rainfall and altitude zones (Molloy, 1988), many of which remain under indigenous vegetation and plantation forest (Table 1), and the perception is that soil quality in these areas is generally not at risk. The Granular, Melanic, and Organic soils were overrepresented and reflect regional concerns that the more intensive land uses on these soil orders may place them in at-risk categories. The main land use of concern on Granular and Melanic soils was arable cropland (including intensive vegetable production), where inorganic nutrient levels can be very high and there has been a marked decline in organic matter (Haynes and Tregurtha, 1999; Schipper and Sparling, 2000; Sparling et al., 2000a,b; Shepherd et al., 2001). Such soils are of concern because of the increased risk to water quality from erosion and nutrient leaching, coupled with deterioration in soil physical and biological condition. Organic soils were of concern to the Regional Councils because peat loss and shrinkage under agriculture infers that this is a nonsustainable use. Pumice soils were the most frequently sampled and the Pumice and Recent soils had the greatest range of land uses (Table 1).

The most frequently sampled land use was pasture for dairy production, followed by pasture for drystock sheep and beef (Table 1). Satellite and aerial photography were not able to distinguish between different uses of pastures for dairy, sheep, and beef, or deer production, and had to be combined to a single category for Table 2. Crop production and horticulture could also not be separated by photography and are combined to a single class. Sixteen percent of the samples in the 500 Soils dataset were from soils used for crop production and horticulture, whereas from aerial and satellite photography this land use was <1% of the current land cover (even when combined). Conversely, low-risk soils are underrepresented; indigenous forest comprised 10% of the samples in the 500 Soils data set despite being 23% of the national vegetation cover (Table 2).

Within the 500 Soils dataset some soil order and land-use combinations are not represented at all (Table 1). This pattern reflects the pragmatic reality of land use within New Zealand, and that the sampling sites were on commercially operated farms rather than research farms with balanced experimental designs. It also reflects the geographic distribution of the participating regions. The strategy of targeted sampling to prioritize at-risk soils is efficient use of monitoring resources, but has led to the over- and underrepresentation of various soil and land use classes when considered on a wider basis. The sites selected provided a good range of soil and land-use comparisons, but caution is needed in aggregating such data into a national overview because of the sampling bias. Such bias can be partially overcome by applying weighting factors based on the area extent of the soil or land use relative to the sampling frequency (Hellkamp et al., 2000).

Land Use and Soil Order
Organic Resources
Total C content was best explained by the soil order (50% of the variability) with a comparatively small (27%) contribution from the land use (Table 3). Total C content (to a 100-mm depth) ranged from 7.8 to 154 Mg ha-1. There were few significant differences in total C between land uses because of high variability. Overall, there was a tendency for total C to be greater under pastures and indigenous vegetation than cropland and plantation forests (Fig. 1a). Land use had a greater effect on total N (43%) with an equivalent contribution from soil order (43%). Total N contents ranged from 0.39 to 10.0 Mg ha-1. Overall, the N content of pastures was greater than of other land uses (Fig. 1b,c). This pattern reflected the substantial organic matter content of permanent pastures, and the build up of organic N through the introduction of exotic legume species, principally white clover (During, 1984). The pattern of organic N accumulation was even more distinct when mineralizable N was considered; the highest values occurred under pastures for dairy and drystock production (Fig. 1d). Values ranged between 0.18 and 380 µg N cm-3, with 45% of the variability explained by land use and 22% by soil order (Table 3). The enhanced N status of the more intensively used soils was shown by the much lower C to N ratio of between 10:1 and 15:1 in managed pastures compared with between 15:1 and 20:1 for the plantation forests, tussock grassland, and indigenous forest (Fig. 1c). The low C to N ratio of managed pasture soils in New Zealand is a consequence of the widespread use of long-term productive ryegrass–white clover leys that maintain high organic matter levels with a high N content (Ledgard and Steele, 1992).


View this table:
[in this window]
[in a new window]
 
Table 3. Percentages of total variance explained by the categorical variables land use, soil order, and their interaction.{dagger}

 
Chemical Properties
Soil pH was influenced more by land use (31%) than soil order (20%), values ranging from 4.1 to 7.4 (Fig. 2a). In their pristine state, many New Zealand soils tend to be acidic (New Zealand Soil Bureau, 1968), and the occurrence of agricultural topsoils with pH > 5.5 reflects the widespread use of lime following the establishment of European farming methods (During, 1984; Cornforth and Sinclair, 1984). The less modified land uses (forests and tussock grassland) have more acid pH, but indigenous plants and introduced plantation forests are tolerant of pH < 5.5, and the soil pH is suitable for that use (Gibson et al., 2000; Parfitt et al., 1997; Schipper and Sparling, 2000).

The Olsen P status of New Zealand soils showed a large range in values, ranging from undetectable on some forest sites, through to >200 µg cm-3 in a soil used for intensive vegetable production (Fig. 2b). However, a surprisingly small proportion of variability in Olsen P was attributed to land use (28%), with a small effect of soil order (13%). A combination of the two categorical variables explained 57% of the variability (see Table 3). There was a clear trend for Olsen P to be substantially greater under cropland and intensive pastures. Olsen P contents under forestry and tussock grassland were lower, although there were numerous outlier points in the range (Fig. 2b). Many unmodified New Zealand soils have a low available P status (Cornforth and Sinclair, 1984). The occurrence of high Olsen P contents in indigenous ecosystems predominantly occurred at sites adjacent to agricultural land without fencing; stock transfer and fertilizer drift could have contributed to the high P status.

The cation exchange capacity ranged from 4 to 69 cmol cm-3 (Fig. 2c). The greater part of the variation (56%) was explained by the soil order, and only 21% by land use (Table 3). Generally, the cation exchange capacity (CEC) in New Zealand topsoils is highly correlated with the organic matter content (Parfitt et al., 1995). However, amount and types of clay minerals also influence CEC. In the 500 Soils data set, the soils with the highest CEC (Fig. 2c) were not those that had the greatest total C contents. The data set included several soils on the Waiareka clay soil, much used for arable crops, which has exceptionally high CEC because of a large proportion of smectite clay (Sparling et al., 2000a). Inclusion of this particular soil means that the correlation between total C and CEC was not as strong as in the larger data set used by Parfitt et al. (1995). Base saturation ranged from 5 to 100%, being markedly lower under tussock grassland, pines, and indigenous forest than other land uses (Fig. 2d).

Physical Characteristics
Bulk density ranged between 0.27 and 1.44 Mg m-3, with 55% of the variability explained by soil order and 30% by land use (Fig. 3a, Table 3). Mixed cropland soils had the highest bulk density. The cropped soils were all sampled at the end of the harvesting operations, but before any further tillage, and probably represent the maximum bulk density for nontrafficked soil. In the majority of cases, bulk densities were well below values that are of concern for root growth, or for adequate drainage and aeration (Drewry et al., 2000; Pabin et al., 1998), although Houlbrooke et al. (1997) caution than on volcanic ash soils, adverse effects on ryegrass (Lolium perenne L.) root and shoot growth can occur at bulk densities < 1.3 Mg m-3.

Total porosity was generally a reverse of the bulk density pattern, with the lowest porosity on the mixed cropland soils (data not shown). Soil order explained 45% of the variability, and land use 29% (Table 3). The maximum porosity was 86%, and the lowest 46%. Macroporosity in these soils ranged from a minimum of 0.95% v/v to a maximum of 46.3% (Fig. 3b). Variability was partly explained by soil order (32%) and partly by land use (36%) (Table 3). Macroporosity is generally considered a more sensitive measure than total porosity because macropores are the dominant pores responsible for aeration and drainage (Dasilva and Kay, 1997), and are preferentially destroyed by compaction (Drewry et al., 2000). In our data set, the soils under mixed cropland, horticulture, and pastures for dairy and drystock production all had a substantial number of examples with macroporosity < 10%. While critical values for macroporosity are not well defined, and differ for different plant species, a precautionary value of >10% macroporosity has been suggested as advisable in New Zealand to maintain pasture production near optimum (Drewry et al., 1999, 2000). Several soils in the data set were below this value, raising concern about compaction by stock treading on dairy pastures and by tillage and vehicle traffic under cropland.

The total available water (Fig. 3d) showed a wide range in values (4–49% v/v), with little differentiation between land uses (21% of the variation) and soil order (25% of the variation, Table 3). Readily available water (Fig. 3c) also showed a wide range (2.6–23.1% v/v) but tended to be lower (median < 8%) under mixed and arable cropland than under other land uses. Even so, land use explained only 19% of the variability (Table 3), soil order explained 26%, and there was strong interaction (60%).

Interpretation
Most authors have found that organizing data by land use category has been a useful approach to show differences in soil quality properties despite spatial variability. Boehm and Anderson (1997) and Wander and Bollero (1999) found differences in soil properties due to tillage management on farms under different crop management systems, despite a range of soil classes and geographic distribution. Hellkamp et al. (2000) grouped sites by land use and management classes, and reported proportions of sites meeting target criteria. The proportion of cropland (12–31%) in the North American states (Brejda et al., 2000a,b; Hellkamp et al., 2000) was very much greater than in New Zealand (<1%). The latter authors weighted their data to adjust for differing sampling frequency and intensity across states and land uses. We also favor reporting soil quality using land use categories, but have also found soil order to be important for interpretation, particularly of organic resources, bulk density, and total porosity. The proportion of sites meeting target criteria (with weighting for area if necessary) is also our intended way to summarize data once the three-year data collection is completed. Established and scientifically justifiable target criteria are a prerequisite for this approach, but setting of such targets has proved contentious (Sojka and Upchurch, 1999; Karlen et al., 2001).

A Minimum Data Set
Principal components analysis using all 12 soil properties in the data set extracted four factors with eigenvalues > 1 that explained 75% of the variability (Table 4). Factor 1 (18.5% explanation) was primarily physical properties (bulk density and porosity), Factor 2 (24.4%) was organic resources (total C and N and mineralizable N), Factor 3 (16.6%) was soil water characteristics (readily and total available water), and Factor 4 (16.5%) was soil chemical characteristics (pH, base saturation, and Olsen P). We appraised the data set to determine whether all soil properties warranted retention. Some items such as pH and base saturation were highly correlated (r = 0.64, P < 0.001, n = 221) and there was no advantage in retaining both. We retained pH in favor of base saturation, as it was easier to measure. The CEC was also discarded because within a soil order much topsoil CEC was explained by total C content (for all soil orders r = 0.45, P < 0.001, n = 221), as was also reported by Parfitt et al. (1995). A further reason for not retaining CEC was that it was not highly responsive to soil management within a soil order (it is not a highly dynamic soil characteristic). Total porosity was dropped because, within soil orders, it was generally inversely related to total C, and was much less responsive than macroporosity. We also did not retain readily and total available water because it is the water storage in the whole soil profile, and not just the 0- to 100-mm depth, that is of importance when assessing water availability for plants (Hargrove, 1988; Rhoton and Lindbo, 1997). To assist with estimating water storage, we recommend recording the profile characteristics of the subsoil and subsurface horizons.


View this table:
[in this window]
[in a new window]
 
Table 4. Rotated pattern loadings after principal component analysis of all properties in the soil quality data set, and after reduction to seven key properties.{dagger}

 
After discarding these items, the minimum data set contains seven key properties: soil pH, total C and N, mineralizable N, Olsen P, bulk density, and macroporosity. Principal components analysis using those seven properties extracted four factors that explained 87% of the variability (Table 4). Factor 1 (37.1%) was primarily organic resources (total C and N, mineralizable N), Factor 2 (18.9%) was dominated by physical properties (bulk density and macroporosity), Factor 3 (15%) was soil P fertility (Olsen P), and Factor 4 (15.7%) was soil acidity (pH). These groupings are similar to those obtained using the full data set, and the four factors split the data in a way that is logical to scientists and land managers. The reduction of complex data sets to four factors using principal components analysis has the potential for simple graphical display of soil quality characteristics in a way that is more easily interpretable by land managers (Romig et al., 1995). Multivariate analyses using principal components analysis (Wander and Bollero, 1999) or discriminant analysis (Brejda et al., 2000a,b) have also been used to group soil properties in the North American studies, and identify which properties were contributing to the separation of the sites. Brejda et al. (2000b) concluded that total organic C and total N were the most sensitive indicators of soil quality, but that total organic C was the only indicator that distinguished between land uses across regions and the two soil types they examined. We also found total C and N to be important contributors to site separation in New Zealand but note that total C was strongly dependent on soil order. Other soil properties, such as potentially mineralizable N, also contributed significantly to the separation, in contrast to the findings of Wander and Bollero (1999), who reported that potentially mineralizable N was not strongly responsive to main effects. It may be that the generally high levels of total C and N and mineralizable N in New Zealand pastoral soils gives these characteristics greater influence during multivariate analysis. Differences between the current study and these others may be due to the wider range of soil orders we sampled. The usefulness of a soil property to distinguish between land uses within a soil order may not apply when examining the impact of different land uses across a number of soil orders. We have subsequently adopted these seven soil properties as our core minimum data set for soil quality monitoring.

Soil Quality and Sustainability
The soil quality properties of concern were those associated with the more intensive land uses of dairy production and arable cropland. Soils under pasture for dairy production and arable cropland showed evidence of soil compaction (decreased macroporosity) relative to other land uses. In more than half the cases, macroporosity was <10% v/v, suggested as a limit below which some deterioration in pasture growth could be expected (Drewry and Paton, 2000; Drewry et al., 1999, 2000; Singleton et al., 2000). We have no information on the consequences of the reduced macroporosity in our study. High levels of inorganic nutrients and substantial organic matter loss were often found under arable cropland, in common with other studies (Boehm and Anderson, 1997; Wander and Bollero, 1999). While organic matter loss is often associated with a decline in desirable soil physical attributes (Soane, 1990; Reeves 1997), high inorganic fertility is undesirable because of the increased risk of leaching losses to the wider environment. The risk is much modified by the low proportion (<1%) of New Zealand soils under cropland, but will be of local concern. There was no evidence in our study that pine plantations had resulted in excessive soil acidification (a perceived risk by land managers in New Zealand). Some pine sites under the third rotation of pine had similar acidity to those under the first rotation, and indigenous vegetation frequently had a lower pH than forest plantations.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
The project is now two-thirds of the way through its three-year duration. We have obtained soil quality samples from 222 sites from five regions in New Zealand. Five additional regions have indicated they will participate in the project in the final year, with an increased intensity of sampling to that already undertaken, and we anticipate the target of 500 sites will be achieved. At this stage of the project, our preliminary conclusions are that the original data set can be simplified to seven key properties meeting the criteria of being readily understood, responsive, of adequate precision, interpretable, cost effective, and internationally acceptable. We acknowledge that there is a bias in the data set in that soils and land uses may be over- or underrepresented on a national basis. However, we believe that the approach of targeted or selective sampling has been very effective in providing soil quality information of interest to the regional users, and increasing their enthusiasm and participation in the project.

The pragmatic, simple, and low-cost approach of the project has been successful in obtaining the participation of the majority of the regions, and has provided quantitative information where none previously existed. We found similar patterns in soil quality characteristics within the same land use across all regions despite differences in soils and climate. Examples were the high total C and N contents under virtually all pastures, loss of macroporosity under pastures, and high chemical fertility under cropped soils. In no region did we find evidence that soil acidity under plantation pine was any lower than that under indigenous forest. The participating regions are waiting until the end of the project to fully assess the implications for future monitoring policy. At present, the data are being used to raise awareness of soil quality issues among regional council staff, scientists, and the general public, and to justify future monitoring. Aspects of the 500 Soils project findings on soil quality under pasture for dairy production have already been incorporated into guidelines for best management practices for the dairy industry (Roberts et al., 2001).


    ACKNOWLEDGMENTS
 
This project was partially funded from the New Zealand Ministry for the Environment Sustainable Management Fund (Contract 5089). We thank the Auckland, Canterbury, Bay of Plenty, Taranaki, and Waikato Regional Councils for their financial support and active participation in site selection and sampling. The cooperation of the landowners for access and sampling on their property and soil management information is gratefully acknowledged. Wim Rijkse, Landcare Research, Hamilton, collected soil samples in the Waikato, Bay of Plenty, and Auckland regions. Hugh Wilde, Landcare Research, Palmerston North, collected samples in Taranaki, and Tony van der Weerden, Crop and Food Research, Lincoln, collected samples in Canterbury. Analyses were completed by Landcare Research staff at Hamilton and Palmerston North.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
J. Environ. Qual.Home page
A. C. R. de Lima, W. Hoogmoed, and L. Brussaard
Soil quality assessment in rice production systems: establishing a minimum data set.
J. Environ. Qual., March 1, 2008; 37(2): 623 - 630.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
G. P. Sparling, D. Wheeler, E.-T. Vesely, and L. A. Schipper
What is Soil Organic Matter Worth?
J. Environ. Qual., March 1, 2006; 35(2): 548 - 557.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
G. Sparling, R. L. Parfitt, A. E. Hewitt, and L. A. Schipper
Three Approaches to Define Desired Soil Organic Matter Contents
J. Environ. Qual., May 1, 2003; 32(3): 760 - 766.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
L. R. Lilburne, A. E. Hewitt, G. P. Sparling, and N. Selvarajah
Soil Quality in New Zealand: Policy and the Science Response
J. Environ. Qual., November 1, 2002; 31(6): 1768 - 1773.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (15)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Sparling, G. P.
Right arrow Articles by Schipper, L. A.
Related Collections
Right arrow Soil Conservation
Right arrow Soil Fertility and Productivity
Right arrow Sustainable Agriculture
Right arrow Best Management Practices
Right arrow Soil Analysis


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
The SCI Journals Agronomy Journal Crop Science
Vadose Zone Journal Journal of Plant Registrations
Journal of Natural Resources
and Life Sciences Education
Soil Science Society of America Journal