Journal of Integrative Agriculture 2018, 17(0): Available online at ScienceDirect

Transcription

1 Journal of Integrative Agriculture 8, 7(): 6-7 Available online at ScienceDirect RESEARCH ARTICLE Updating conventional soil maps by mining soil environment relationships from individual soil polygons CHEG Wei,, ZHU A Xing,,,, QI Cheng-Zhi,,, QI Feng 6 State Key Lab of Resources and Environmental Information System, Institute of Geographic Sciences and atural Resources Research, Chinese Academy of Sciences, Beijing, P.R.China University of Chinese Academy of Sciences, Beijing 9, P.R.China Key Laboratory of Virtual Geographic Environment, Ministry of Education, anjing ormal University, anjing, P.R.China Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, anjing, P.R.China Department of Geography, University of Wisconsin-Madison, Madison 76, USA 6 School of Environmental and Sustainability Sciences, Kean University, Union 78, USA Abstract Conventional soil maps contain valuable knowledge on soil environment relationships. Such knowledge can be extracted for use when updating conventional soil maps with improved environmental data. Existing methods take all polygons of the same map unit on a map as a whole to extract the soil environment relationship. Such approach ignores the difference in the environmental conditions represented by individual soil polygons of the same map unit. This paper proposes a method of mining soil environment relationships from individual soil polygons to update conventional soil maps. The proposed method consists of three major steps. Firstly, the soil environment relationships represented by each individual polygon on a conventional soil map are extracted in the form of frequency distribution curves for the involved environmental covariates. Secondly, for each environmental covariate, these frequency distribution curves from individual polygons of the same soil map unit are synthesized to form the overall soil environment relationship for that soil map unit across the mapped area. And lastly, the extracted soil environment relationships are applied to updating the conventional soil map with new, improved environmental data by adopting a Soil Land Inference Model (SoLIM) framework. This study applied the proposed method to updating a conventional soil map of the Raffelson watershed in La Crosse County, Wisconsin, United States. The result from the proposed method was compared with that from the previous method of taking all polygons within the same soil map unit on a map as a whole. Evaluation results with independent soil samples showed that the proposed method exhibited better performance and produced higher accuracy. Keywords: update conventional soil map, soil environment relationships, knowledge extraction, individual soil polygons Received January, 8 Accepted 8 March, 8 CHEG Wei, chengw@lreis.ac.cn; Correspondence ZHU A Xing, azhu@wisc.edu 8 CAAS. Publishing services by Elsevier B.V. All rights reserved. doi:.6/s9-9(8)698-. Introduction Information on the spatial distribution of soils is increasingly important for watershed management and ecological modeling applications (McBratney et al. ; Brus et al.

2 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 8). However, polygon-based conventional soil maps with crisp boundaries cannot satisfy the model requirements for more detailed and accurate soil data (Zhu and Mackay ; Quinn et al. ; Li et al. ). Previous studies have shown that conventional soil maps can be used as a source of both data and knowledge in digital soil mapping (DSM, McBratney et al. ), and be updated to higher accuracy (Häring et al. ; Kerry et al. ; Subburayalu and Slater ; Collard et al. ; Heung et al., 6; Subburayalu et al. ). One way of updating conventional soil maps involves extracting the soil environment relationships embedded in these maps (Bui et al. 999; Moran and Bui ; Qi and Zhu ; Yang et al. ). The extracted knowledge is then applied to updating soil maps with new and improved data on environmental factors which co-vary with the soil forming process (Jenny 9, 96). Currently there are two main approaches to extracting such embedded knowledge from soil maps. One is to select training samples from conventional soil maps, and then use a machine learning algorithm (such as decision trees, classification trees, or random forest) to generalize the soil environment relationships represented by these samples. The training samples are often randomly selected using different methods, such as selecting equal number of training samples for each soil map unit (Moran and Bui ), selecting equal number of training samples for each polygon (Subburayalu and Slater ; Heung et al., 6; Odgers et al. ; Subburayalu et al. ; Holmes et al. ), or selecting training samples according to an area-weighted proportion of each map unit extent over entire study area (Bui and Moran, ; Moran and Bui ; Grinand et al. 8; Heung et al., 6; Vincent et al. 8). Randomly sampling from conventional soil maps, however, cannot avoid the inclusion of low quality samples that may result from mapping errors, such as incorrect labels, inclusions and misplacements of boundaries (Ehlschlaeger and Goodchild 99). To guarantee the selection of high quality, representative training samples from conventional soil maps, it has been proposed to select samples only from typical areas for each map unit. The locations of typical areas of each map unit can be represented by the clustering centroids using clustering algorithms (Yang et al. ) or the modes of frequency distribution histograms on individual environmental covariates for each map unit (Qi and Zhu ; Qi ). The other approach to extracting the knowledge on soil environment relationships from soil maps is through constructing frequency distribution curves of individual environmental covariates (or, environmental frequency curve) for each map unit by overlaying soil map with data layers of the soil forming environmental covariates. This is based on a reasonable assumption that the majority of the polygon area on the conventional soil map is correctly categorized (Qi and Zhu ). The resulted environmental frequency curves indicate the degree that the environmental conditions are suitable for the development of the mapped soil types (Qi et al. 8; Du et al. ). Qi and Zhu () and Silva et al. (6) compared the above two approaches to extracting knowledge from existing soil maps and updating these maps. Extracting soil environment relationships as individual environmental frequency curves have achieved higher accuracy than sampling the original map. Previous studies derived environmental frequency curves from existing soil maps by taking all polygons within the same soil map unit on a map as a whole. The resulting soil environment relationships inevitably represent bias toward larger polygons in the map unit. As the basic spatial unit of a conventional soil map, each polygon represents soil surveyors understanding of the local soil and its environment at that particular location. This means that each individual polygon, no matter small or large, represents some specific environmental conditions for the development of that particular soil type. Grouping the smaller polygons with larger ones within the map unit to construct environmental frequency curves leaves the smaller areas underrepresented at the low frequency ends. It introduces bias of the soil environment relationships favoring environmental conditions represented by polygons that occupy larger areas spatially. This bias may influence the performance of DSM when using such extracted soil environment relationships to update a conventional soil map. In this paper we present a method of mining soil environment relationships as environmental frequency curves from individual soil polygons. With a case study we used the method to update a conventional soil map of the Raffelson watershed in La Crosse County, Wisconsin, United States. The proposed method was then compared with previous method that extracts the soil environment relationship only at the map unit level.. Methods.. Basic idea As in previous studies (Qi and Zhu ; Qi ) which assume that normally the majority in polygon area on the conventional soil map is correctly categorized, environmental frequency curves (Qi et al. 8; Du et al. ) are used to represent the soil environment relationships extracted from mapped soil polygons in this study. To obtain the final soil environment relationships, we first calculate the frequency distribution of each environmental covariate for each individual soil polygon and

3 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 construct the corresponding environmental frequency curve. ext, for each environmental covariate, these environmental frequency curves from individual soil polygons within the same map unit are synthesized through combining the categories together for a categorical environmental covariate or selecting the max frequency value associated with the same environmental covariate value for a continuous environmental covariate. ote that a map unit on a conventional soil map sometimes contains more than one soil type (complex map units). The environmental frequency curves from individual soil polygons within same map unit but represent different soil types should not be synthesized. In such case, we introduce a parent material variable to supplement the map unit delineations. Parent material, as a categorical covariate, plays a decisive role in controlling the development of soil types. One soil type is normally developed on only one parent material, while the same parent material usually develops into different soil types. Overlaying the conventional soil map with a parent material or sometimes bedrock geology layer which is often available for the mapped area, we obtain combinations of soil map unit and parent material called the map unit parent material pairs that help break down complex map units. Environmental frequency curves are then extracted from these map unit parent material pairs to differentiate soil types within complex map units developed on different parent material. For the parent material environmental covariate itself, environmental frequency curves will still be synthesized at the map unit level. Whereas for other environmental covariates, the environmental frequency curves from individual soil polygons will be synthesized at the level of map unit parent material pair. The detailed steps are elaborated below... Step : Extracting soil environment relationships implied by each individual soil polygon The soil environment relationships implied by individual soil polygons are first extracted as environmental frequency curves. The form of these environmental frequency curves is dependent on the type of environmental covariate. For a categorical environmental covariate, a step function with a membership value of for the optimal category and for other categories is used for the environmental frequency curve. The dominant category (i.e., the one with the highest area proportion) in a soil polygon is assigned membership for the soil type at that location. As aforementioned, we overlay the soil map with a parent material map to first create a map unit parent material composite map so that each polygon on this map has both a map unit and a parent material attribute. ote that small polygons located at the border of the soil map are often incomplete polygons whose dominant part may not be included in the map. Using such incomplete polygons to extract the soil environment relationships may introduce errors. Therefore we recommend a data preprocessing step in this study to exclude these small polygons at the map borders. For other categorical environmental covariates (such as land use or vegetation) that come in secondary in terms of their impact on soil formation, the categorical map of each covariate is then overlapped with the above-generated soil map unit parent material composite map. The dominant category (i.e., the one with the highest area proportion) for each soil polygon on the composite map is finally identified to be assigned a membership of. For a continuous environmental covariate (such as elevation, slope gradient), a curve is used to fit its frequency distribution within a soil polygon on the soil map unit parent material composite map. There are parametric or non-parametric methods to fit the curve. Du et al. () experimented with both methods to fit frequency distribution curves for continuous covariates and the results exhibited no significant difference when used to predicting soil type or soil property. We adopt a non-parametric method, i.e., kernel density estimation (KDE) (Gatrell et al. 996; Law et al. 9), in this study as it does not need to assume curve shapes in advance. The basic assumption of KDE is that a soil type associated with an optimal environmental covariate value (i.e., x i ) may occur at locations with a different environmental covariate value (i.e., x) at a lower probability. The probability decreases as the distance between x i and x increases. The following equation shows the probability density estimation with KDE: ^ f p n x x j (x)= K( i ) () h= nh nh Where, f ^ p j (x) is the probability density associated with environmental covariate value x among a set of environmental covariate values (x, x,..., x i,..., x n ) in polygon p j, n is the number of cells in that polygon. K( ) is a kernel density function whose value (the probability) is dependent on the distance between x i and x. A commonly used function is the Gaussian kernel (Silverman 986). h is the bandwidth which is calculated with the rule of thumb algorithm (Silverman 986). KDE estimates the probability density associated with a particular environmental covariate value x by summing up the density contributions from all environmental covariate values (Diggle 98; Fotheringham et al. ). The derived frequency f ^p j (x) (the probability density) values are then normalized to f ^ p j (x) ([,]) using the following equation: ^p ^ ^ f p j (x)= f j p (x) f j min(x) () ^p ^ f j p (x) f j (x) max min

4 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 Where, ^ f p j (x) is the normalized frequency value associated with environmental covariate value x within polygon p j, the subscript means normalized. After normalization, each environmental covariate value x will correspond to a normalized frequency value. The highest normalized frequency value is which means the corresponding environmental covariate value is the optimal value for the soil polygon under consideration... Step : Synthesizing soil environment relationships across polygons Synthesizing the relationship between soil and parent material at the map unit level Theoretically, when a soil type is found to have developed only on a certain parent material, it should not distribute outside the spatial range of this parent material. However, due to errors on conventional soil maps (such as incorrect labels, and inaccurate placements of polygon boundaries), the soil map unit parent material composite map resulted from the last step might contain some unrealistic combinations between soil map unit and parent material. These unrealistic combinations should be excluded when synthesizing knowledge derived from each individual polygon. This is realized through selecting only the dominant parent material within each map unit. Specifically, map unit with single soil type should correspond to one parent material only. Under this assumption, we add the polygons with the same parent material together to get the area of each parent material and select the parent material with the highest area proportion as the optimal one for this map unit, the minor ones which may well be results of mapping errors are discarded. Again, a step function with a membership value of for the optimal parent material and for other parent materials is used for the synthesized soil parent material relationships of the map unit. Complex map units contain more than one soil type that could have developed on different parent materials. In this case, the area proportions of different parent material categories are calculated to reflect the dominant ones in the complex map unit (Fig. ). During this process, a threshold (i.e., thr in Fig. ) is set to judge if the inclusion of a parent material is trivial (due to map errors) or significant enough to represent an actual inclusion. An area proportion of %, which has been used as a threshold value in Soil Survey Manual (SSDS 7) to determine whether a component is a minor one, was used in the case study. Sometimes however, the total area of a parent material (e.g., b+c, in Fig. ) is small, so that its area proportion occupied in a specific map unit (e.g., b/(a+b), in Fig. ) may not be significant but shows a large proportion of its own area (e.g., b/(b+c) in Fig. is larger than thr in Fig. ), it should not be neglected and is still added to be one of the parent Complex map unit Parent material map Area proportion of each parent material on this map unit Area proportion of a parent material on this map unit >thr? Yes o Area of a parent material includes in this map unit Total area of this parent material across the study area >thr? Yes Corresponding parent material categories o ot the corresponding parent material category Fig. Flowchart of synthesizing the corresponding parent material group for a complex map unit.

5 6 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 material categories for this map unit. In our case study, a parent material category included in the complex map unit that occupies more than half of its own area is considered to be such case (i.e., thr=%). After all parent material categories are identified for a complex map unit a similar step function as before is then used to represent the synthesized soil parent material relationships of this complex map unit (Fig. -A). This synthesized step function for the complex map unit can also be decomposed to separate step functions (Fig. -B and C), each for one map unit parent material pair and can be conveniently applied to updating the conventional soil map in Step. Synthesizing the relationship between soil and other environmental covariates at the level of map unit parent material pair For each categorical environmental covariate other than parent material, the dominant category derived from each individual polygon with the same map unit parent material pair is combined together to constitute the category combination for each map unit parent material pair. For a continuous environmental covariate, a fuzzy MAXoperator is adopted to synthesize the frequency distribution curves derived from each individual soil polygon within the same map unit parent material pair (Fig. ): ^ p f (x)=max{^ f (x), f ^ p (x),..., f ^ p j (x),..., f ^ p m (x)} () Where, f ^ (x) is the membership value ([,]) to a map unit parent material pair for a particular environmental covariate value x, m is the number of polygons within this map unit parent material pair. The resulted curve depicts both the optimal environmental covariate value for each polygon within this map unit parent material pair and the change of membership when the environmental covariate value differs from the optimal value... Step : Updating the conventional soil map After the relationships between soil and each environmental covariates are extracted and synthesized in the form of a set of environmental frequency curves, the soil land inference a b c Complex map unit A Parent material I a+b: the total area of complex map unit A b+c: the total area of parent material I Fig. Example of a complex map unit overlapped with a small parent material polygon. B A Optimality Optimality P P P P P Parent material P P P P P C Parent material Optimality P P P P P Parent material Fig. Example of a synthesized step function for a complex map unit and the decomposition of it to separate step functions.

6 7 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 model (SoLIM) (Zhu et al. 997; Zhu et al. ) can be adopted to update the soil map. SoLIM is a DSM model based on the similarity of environmental conditions. It combines the knowledge on soil environment relationships with detailed environmental data layers to infer the spatial distribution of soils (Zhu et al. 996; Zhu 997, 999; Zhu et al. ). Under the SoLIM framework, those synthesized environmental frequency curves on soil map unit parent material pairs are used as fuzzy membership functions to compute the corresponding similarity values between each cell in the mapping area and each map unit parent material pair. This process produces a set of fuzzy membership maps, with each of them shows the spatial variation of membership to a certain map unit parent material pair across the area. A crisp map of soil map units could also be created by hardening the set of fuzzy membership maps, which is achieved by assigning each location the soil map unit bearing the highest membership value at that location (Zhu 997).. Case study.. Study area and data Our study site is the Raffelson watershed in La Crosse County, Wisconsin, United States. It covers a small area of approximately km, and the elevation ranges from m to m (Fig. ). The study area is in the so-called driftless area of southwestern Wisconsin, which has remained free of direct impact from late Pleistocene era continental glaciers. As a result, it has a typical ridge and valley terrain with relatively flat, narrow ridges and wide. ^ fp(x) Frequency ^.8 fp(x).6 fp(x) ^ fp(x) ^ ^ f (x) Elevation (m) Fig. An example of the synthetization of continuous environmental frequency curves. Evaluation samples Elevation (m) 6 Fig. Map of the Raffelson watershed in La Crosse County, Wisconsin, United States. m

7 8 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 valley. The slope gradient in the study area ranges from to 6% and slopes with gradients below % accounts for about % of the area. The land use in the ridges and valley region is primarily cultured land while the side slopes is predominantly covered by forests. The conventional soil map at :, scale of the study area (Fig. 6) was obtained from the Soil Survey Geographic (SSURGO) database created by the atural Resources Conservation Service, United States Department of Agriculture. The map includes map units, which are Valton (Val), Lamoile (Lam), Churchtown (Chu), orden (or), Greenridge (Gre), Council (Cou), Kickapoo (Kic), Orion (Ori), Hixton (Hix), Dorerton-Elbaville (Dor-Elb), Gaphill-Rockbluff (Gap-Roc), and Council-Elevasil-orden (Cou-Ela-or). Fig. 6 shows a rasterized format of this soil map to be used in the case study. It shows that some map units (i.e., Dor-Elb, Gap-Roc, and Cou-Ela-or) are complex map units. Five of the map units (i.e., Chu, or, Cou, Dor-Elb, and Cou-Ela-or) cover 79% of the study area. The environmental covariates used for knowledge extraction and soil map updating were selected based on the relevant soil forming factors in the study area. In this study area, only one categorical environmental covariate was considered: parent material. According to the soil survey report, bedrock geology in this area is complex and plays an important role in the soil formation. We thus used a bedrock geology map (Fig. 7) to depict the parent material conditions. Fig. 7 shows five distinct types of parent materials within the study area: Oneota (dolomite), Glauconite (sandstone), Jordan (sandstone), Wonewoc (sandstone), and alluvial materials (Alluvial for short). Among them, Oneota, Glauconite, and Alluvial occupy 9.6% of the total area, while the remaining two (i.e., Jordan and Wonewoc) only take up. and.% respectively. Land use and vegetation were not chosen in this case study, because they are disturbed seriously by human activities (Zhu 8) and the relationships between them and soil are not clear in this area. Topography is a primary factor that determines the soil variation in the study area. This is confirmed by the SSURGO soil map which shows clear correlation between soil and terrain conditions. In the case study, elevation, slope gradient, and topographic wetness index were chosen to characterize topography and were derived from a DEM of -meter resolution produced by the United States Geological Survey (Qin et al. ). Aspect was not chosen because soil polygons on the conventional soil map (Fig. 6) almost wrap around slopes with all aspects, which means that aspect has no significant indication for soil changes in the case study area. In the study area, a total of 99 field samples (Fig. ) were collected as the independent evaluation set. Among them, samples were collected along a few transects to cover different landscape positions and soil types in the study area. The remaining 6 samples were scattered to cover the major landscape units (such as ridge tops, side slopes, and valley bottoms) throughout the watershed. Sixteen soil types were identified from these samples. Fourteen of them were shown on the conventional soil map and the other two (containing seven samples; Table ) were not... Knowledge extraction Small polygons located at the border of the soil map were Map unit Valton Lamoille Churchtown orden Greenridge Council Kickapoo Orion Hixton Dorerton-Elbaville Gaphill-Rockbluff Council-Elevasil-orden m Fig. 6 Rasterized conventional soil map of the Raffelson watershed in La Crosse County, Wisconsin, United States.

8 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): Bedrock type Oneota Jordan Glauconite Wonewoc Alluvial m Fig. 7 Bedrock geology map of the Raffelson watershed in La Crosse County, Wisconsin, United States. Table Parent material(s) corresponding to each map unit, extracted by the proposed method and the reference method which extract soil-environment relationships only at map unit level Parent material(s) Soil map unit By the reference method By the proposed method Val Oneota Oneota Lam Oneota Oneota Chu Alluvial Alluvial or Glauconite Glauconite Gre Glauconite Glauconite Cou Alluvial Alluvial Kic Alluvial Alluvial Ori Alluvial Alluvial Hix Wonewoc Wonewoc Dor-Elb Oneota Oneota, Alluvial Gap-Roc Alluvial Oneota, Jordan, Alluvial Cou-Ela-or Alluvial Glauconite, Wonewoc, Alluvial excluded first during preprocessing as mentioned in Section.. These polygons accounts for.% of the total area. After that the relationships between soil and parent material were extracted as detailed in Section.. Specifically, Dor-Elb is a complex map unit that contains two soil types, Dorerton and Elbaville. The area proportions of two dominant parent material Oneota and Alluvial in this map unit are 67 and 8%, respectively. The corresponding parent material categories were determined to be Oneota and Alluvial. Gap-Roc is a complex map unit that contains two soil types, Gaphill and Rockbluff. Three dominant parent materials (i.e., Oneota, Jordan, and Alluvial) were found in this map unit, with area proportions of,, and 9%, respectively. All three were determined as parent material categories for Gap-Roc. Cou-Ela-or is also a complex map unit that contains three soil types, Council, Elevasil, and orden. Two parent materials Glauconite and Alluvial cover. and % of the map unit respectively. A third parent material, Wonewoc, covers only % of this map unit, but the area included in this map unit accounts for 7% of all Wonewoc in the study area. It thus should not be neglected. All three were identified to be corresponding parent material categories for this map unit. ext, the step functions of the relationships between map units and their corresponding parent materials were constructed as described in Section.. And the decomposed step function between a complex map unit and each of its corresponding parent material was also constructed. The result is 7 map unit parent material pairs (i.e., step functions) (Table, Fig. 8). The area covered by these 7 map unit parent material pairs accounts for 88.9% of the total area. The rest area contains either border polygons or combinations between soils and parent material that is not typical for the soils and thus considered to be unreasonable and excluded from the steps below. For each continuous environmental covariate, the frequency distribution curves implied by individual soil polygons were constructed by the proposed method detailed in Section. and then synthesized for the map unit parent material pair as described in Section.. With the 7 identified map unit parent material pairs and the selected environmental covariates, a total of 68 environmental frequency curves were constructed. These environmental frequency curves, which represent knowledge of soil environment relationships embedded in the conventional soil map, were input to SoLIM to update the soil map with the newly obtained environmental data layers.

9 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 Table umber of samples for each soil type correctly classified in the conventional soil map and in soil maps updated with the reference method, the proposed method, and the proposed method with two possible modifications Soil type Valton Lamoille Dorerton Elbaville Churchtown Gaphill Rockbluff orden Greenridge Urne Council Kickapoo Orion Elevasil Hixton Boone Total umber Conventional soil map 7 6 Reference method Proposed method Proposed method with the Proposed method with the first modification second modification Soil map unit parent material Val Oneota Lam Oneota Chu Alluvial or Glauconite Gre Glauconite Cou Alluvial Kic Alluvial Ori Alluvial Hix Wonewoc Dor-Elb Oneota Dor-Elb Alluvial Gap-Roc Oneota Gap-Roc Jordan Gap-Roc Alluvial Cou-Ele-or Glauconite Cou-Ele-or Wonewoc Cou-Ele-or Alluvial Unreasonable area m Fig. 8 Distribution of 7 map unit parent material pairs of the Raffelson watershed in La Crosse County, Wisconsin, United States... Evaluation The proposed method was compared to the previous method that extracts soil environment relationship without differentiating individual polygons (called the reference method) for evaluation. To be consistent, the reference method was applied to extract environmental frequency curves for each map unit parent material pair as well instead of just map unit. ote that when the reference method was applied to extracting the relationship between soil and parent material in the study area, there was no soil map unit corresponding to Jordan due to its small area in this study area (Table ). For the sake of completeness, the environmental frequency curve between soil and parent material derived with the proposed method was also used for the reference method in our comparison. With the 7 map unit parent material pairs and environmental covariates, a total of 68 environmental frequency curves were constructed with the reference method and were input to SoLIM to update the conventional soil map. The soil map updated with the proposed method was compared to that with the reference method in two aspects. First, the difference in spatial patterns was analyzed to evaluate whether the proposed method was able to capture the detailed spatial change of soils across the study

10 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 area better than the reference method. Second, the 99 independent samples (Fig. ) were used to evaluate map accuracies. When a sample located in a complex map unit matches any of the soil types in this complex map unit, it is considered to be correctly mapped. Under the assumption of independent classification errors, the number of correctly classified samples should satisfy a binomial distribution. When comparing the two methods, a null hypothesis can be that the accuracy of the soil map updated with the proposed method is not higher than that of the reference soil map. When the number of correctly classified samples follows the binomial distribution B(n, p), where, p is the proportion of correctly classified samples from the reference soil map and n is the total number of samples. The probability, P, with k correctly classified samples is calculated using the following formula: n () P(X=k)=( )pk( p)n k k The probability of having k or more correctly classified samples is calculated using the following equation: n P(X k)= X=k P(X), X=,,..., K,..., n () If P(X k) is smaller than. (or.), the null hypothesis should be rejected and the accuracy improvement is statistically significant at the. (or.) level. Otherwise, the accuracy improvement is not significant at this statistical significance level. method. For those map units with multiple polygons on the conventional soil map (such as Greenridge, Kickapoo, orden), the results from the two methods are notable. The proposed method produced more, smaller patches for these map units than the reference method did as it captured more detailed environmental conditions for local soil distribution from individual soil polygons. Table shows that the accuracy of the soil map updated with the proposed method is higher than that with the reference method, based on the 99 evaluation samples. This is an indicator that the soil environment relationships extracted from the proposed method could capture the relationship between soil and its environmental conditions better than the reference method. However, the statistical test results (P-value) of the accuracy improvement from the conventional map to the map updated with the proposed method and that from the map updated with the reference method to the map updated with the proposed method are. and.7 respectively and are both not statistically significant at the. level (Table ). In order to further improve the proposed method, we need to look into the implementation of the method for possible modifications.. Discussion.. Comparing the environmental frequency curves extracted with the proposed and the reference methods.. Evaluation results The soil maps updated with the proposed method and the reference method are shown in Fig. 9. The spatial patterns of map units on the two maps are consistence in general. The main differences are found in the border areas of map units. For map units with very few polygons on the conventional soil map (such as Valton and Lamoille), there is little difference between results from the proposed method and the reference method. This is natural because the environmental frequency curves for map unit parent material pair with only one polygon are the same with the proposed method and the reference A The number of samples correctly classified with the proposed method is six more than that with the reference method (Table ). Most of the accuracy improvement occurred in the Churchtown map unit, which contains only one soil type and accounts for.9% of the total study area. Here, the Churchtown map unit was taken as an example to discuss the performance of the proposed method in more detail. Churchtown was developed on the Alluvial parent B Map unit o data Valton Lamoille Churchtown orden Greenridge Council Kickapoo Orion Hixton Dorerton-Elbaville Gaphill-Rockbluff Council-Elevasil-orden m m Fig. 9 Soil maps updated with the reference method (A) and the proposed method (B).

11 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 material, and out of the 99 evaluation samples were identified as this soil type (Table ). Five and nine samples were correctly classified by the soil maps updated with the reference method and with the proposed method, respectively. Fig. shows that the environmental frequency curve for elevation of the Churchtown Alluvial pair extracted with the reference method has a much narrower peak value region. This curve thus represents a more centralized distribution for the elevation environmental covariate. This shows that the reference method tends to ignore or underestimate the variance among individual polygons within the same map unit parent material pair, which is also the case with other environmental covariates. On the other hand, the peak region of the environmental frequency curve extracted by the proposed method is much broader (Fig. ), which resembles the empirical knowledge obtained from the local soil expert (the green dotted curve in Fig. ; Zhu 8) and the distribution of actual elevation values of the evaluation samples recognized as this soil type... Possible modifications to the proposed method to further improve performance The proposed method is able to extract environmental frequency curves that have wider peak regions to reflect the distribution of soils in the feature space defined by the environmental covariates, but the so-extracted knowledge is still limited to the local situation. The optimal environmental condition for the development of a certain soil type may not be completely represented by the mapped polygons and the distribution of the optimal environmental covariate value range for the soil distribution could be even wider (Qi and Zhu ; also see the dotted curve of expert knowledge obtained from the local soil expert as shown in Fig. ). Inference with a narrow environmental frequency curve could result in low membership values for large area of the soil type (also map unit) and hence possibly a restrained spatial extent for the soil type (Qi and Zhu ). In this section we discuss two possible modifications to the proposed method for continuous environmental covariates to address the limitation mentioned above... The two possible modifications to the proposed method The first possible modification is to simply extend the optimal value range of environmental covariates to make sure that the value ranges of these environmental covariates for each parent material are completely covered. Specifically, for a continuous environmental covariate, the environmental frequency curves extracted within the same parent material are first aligned along the environmental covariate axis, as shown in Fig. -A. It is shown that within the value range, there exist blank area and low frequency regions at both Table P-values of the accuracy improvement of the soil map updated with the proposed method and with two possible modifications of the proposed method compared with the conventional soil map and the soil map updated with the reference method Reference method Proposed Proposed method with the first Proposed method with the second method modification modification Conventional soil map.. *.6 ** Reference method.7..8 * Proposed method.7. * and ** indicate the accuracy improvement is statistically significant at the. and. level, respectively. Frequency Expert knowledge The reference method The proposed method Evaluation samples. Elevation (m) Fig. Environmental frequency curves for elevation of the Churchtown Alluvial pair extracted with the proposed and the reference methods. The dotted curve of expert knowledge is obtained from the local soil expert according to Zhu (8), elevations for the evaluation samples recognized as this soil type are also shown.

12 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 ends. This is unreasonable as besides the three curves showed on the figure, no other soil type developed on this parent material. Therefore the curves located on both ends of the value range should have optimal membership to complete the feature space. Thus we modify the leftmost and rightmost curves to be Z-shaped or S-shaped (Zhu et al. ), respectively (Fig. -B). To avoid the possible situation that when both sides of a curve are modified it might cover other curves unintentionally, the immediate adjacent curve will be modified on one side instead. Furthermore, if there is only one map unit corresponding to a certain parent material, this map unit should have the optimal membership to each environmental covariate s entire value range of that parent material. This means the original curve is modified A Frequency B Frequency Map unit A parent material I Map unit B parent material I Map unit C parent material I The elevation range of parent material I 6 8 Elevation (m) Elevation (m) C. Frequency Elevation (m) Fig. Examples of possible modifications on extracting the environmental frequency curves. A, the environmental frequency curves of elevation for three soil map unit parent material pairs developed on the same parent material extracted with the proposed method. B, the first possible modification. C, the second possible modification. into a straight line. The second possible modification is a slight variant from the first and is used when more than one curve is located very close to either end of the axis, such as the gray and green curves on the right end as illustrated in Fig. -A. In this case, these curves with adjacent peaks are supposed to have the same possibilities to get the optimal membership associated with those environmental covariate values. So these curves are grouped together and modified to be Z-shaped or S-shaped in a similar way as above. An example is shown in Fig. -C. More specifically, for the rightmost curve, if another curve intersects with its right downslope side and the frequency value at the intersection point is no smaller than. (see the red point in Fig. - A), the intersecting curve is modified together. The same applies to the leftmost curve... Improvement with the two modifications We experimented with the two modifications following extraction of environmental frequency curves using the proposed method. Table shows that the accuracy of the updated map was improved by these two modifications, with the second modification improved the proposed method more. As shown in Table, compared to the conventional soil map, the accuracy improvements from the proposed method with the first modification and the second modification are statistically significant at the. and. level respectively. Compared to the soil map updated with the reference method, the accuracy improvements from the proposed method with the second modification is also statistically significant at the. level. The number of samples correctly classified with the proposed method plus the first modification is two more than that with the proposed method alone (Table ). The most accuracy improvement is from the Gaphill-Rockbluff complex map unit. One of the two samples was recognized as Rockbluff and the other Gaphill in the field. Both were inferred as odata on the soil map updated with the proposed method without any further modification. This is because those environmental frequency curves cannot cover the environmental covariate value range of its parent material. For those locations with uncovered environmental covariate values, the membership values to any map unit parent material pair are. Thus SoLIM set odata for these locations on the soil map. With the first modification, these two samples were correctly mapped to the Gaphill- Rockbluff map unit. The second modification further improved accuracy by classifying two more samples correctly (Table ). The two samples are found to be of the orden map unit. This map unit accounts for 6.% of the total study area on

13 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 the conventional soil map and is mostly developed on the Glauconite parent material. A total of evaluation samples were identified as this soil type. The two modifications to the proposed method correctly classified and of these samples respectively. The dotted curve on Fig. shows the responsible environmental frequency curve for the elevation covariate of the orden Glauconite pair. It also shows that the curve shape stayed unchanged with the first modification (see the gray curve in Fig. -B) and modified to S-shaped with the second modification (Fig. -C). The second modification made the curve cover the distribution of actual elevation values of the evaluation samples (Fig. ), which indicates its effectiveness. Frequency The proposed method and the first modification The second modification Evaluation samples 6 8 Elevation (m) Fig. Environmental frequency curves for elevation of the orden Glauconite pair extracted with the proposed method and two possible modifications to the proposed method (elevations for the evaluation samples recognized as this soil type are also shown).. Conclusion This paper proposes a method of mining soil environment relationships as environmental frequency curves from individual soil polygons to update conventional soil maps. Compared with previous method which takes all polygons within the same soil map unit on a map as a whole during knowledge extraction, the proposed method performed better in our case study. The environmental frequency curves extracted by the proposed method approximates the expert knowledge better and can produce a more accurate soil map when used to update the conventional soil map with detailed environmental data layers. Modifications to the proposed method also showed that the performances could be further improved. It s worth noting that the proposed method does not disaggregate a complex map unit which contains more than one soil type. In our future work, for a complex map unit, we could try to disaggregate the synthesized soil environment relationships into individual soil type components. Acknowledgements This study is supported by the ational atural Science Foundation of China (77 and 9) and the Innovation Project of State Key Laboratory of Resources and Environmental Information System of China (O88RACYA). Supports to A-Xing Zhu through the Vilas Associate Award, the Hammel Faculty Fellow Award, the Manasse Chair Professorship from the University of Wisconsin-Madison, and the One-Thousand Talents Program of China are greatly appreciated. References Brus D J, Bogaert P, Heuvelink G B M. 8. Bayesian maximum entropyprediction of soil categories using a traditional soil map as soft information. European Journal of Soil Science, 9, Bui E, Loughhead A, Corner R Extracting soillandscape rules from previous soil surveys. Australian Journal of Soil Science, 7, 9 8. Bui E, Moran C J.. Disaggregation of polygons of surficial geology and soil maps using spatial modeling and legacy data. Geoderma,, Bui E, Moran C J.. A strategy to fill gaps in soil survey over large spatial extents: An example from the Murray- Darling Basin of Australia. Geoderma,,. Collard F, Kempen B, Heuvelink G B M, Saby P A, Richer de Forges A C, Lehmann S, ehlig P, Arrouays D.. Refining a reconnaissance soil map by calibrating regression models with data from the same map (ormandy, France). Geoderma Regional,,. Diggle P. 98. A kernel method for smoothing point process data. Journal of the Royal Statistical Society (Series C: Applied Statistics),, 8 7. Du F, Zhu A X, Band L E, Liu J.. Soil property variation mapping through data mining of soil category maps. Hydrological Processes, 9, 9. Ehlschlaeger C R, Goodchild M F. 99. Dealing with uncertainty in categorical coverage maps: Defining, visualizing, and managing errors. In: Proceedings of the Workshop on Geographical Information Systems at the Conference on Information and Knowledge Management. Gaithersburg, Maryland. pp Fotheringham A S, Brunsdon C, Charlton M.. Quantitative Geography: Perspectives on Spatial Data Analysis. Thousand Oaks, Sage, CA. Gatrell A C, Bailey T C, Diggle P J, Rowlingson B S Spatial point pattern analysis and its application in geographical epidemiology. Transactions of the Institute of British Geographers,, 6 7.

14 CHEG Wei et al. Journal of Integrative Agriculture 8, 7(): 6-7 Grinand C, Arrouays D, Laroche B, Martin M P. 8. Extrapolating regional soil landscapes from an existing soil map: Sampling intensity, validation procedures, and integration of spatial context. Geoderma,, 8 9. Häring T, Dietz E, Osenstetter S, Koschitzki T, Schröder B.. Spatial disaggregation of complex soil map units: A decision-tree based approach in Bavarian forest soils. Geoderma, 8, 7 7. Heung B, Bulmer C E, Schmidt M G.. Predictive soil parent material mapping at a regional-scale: A Random Forest approach. Geoderma,,. Heung B, Ho H C, Zhang J, Knudby A, Bulmer C E, Schmidt M G. 6. An overview and comparison of machine-learning techniques for classification purposes in digital soil mapping. Geoderma, 6, Holmes K W, Griffin E A, Odgers P.. Large-area spatial disaggregation of a mosaic of conventional soil maps: Evaluation over Western Australia. Soil Research,, Jenny H. 9. Factors of Soil Formation: A System of Quantitative Pedology. McGraw-Hill, ew York. Kerry R, Goovaerts P, Rawlins B G, Marchant B P.. Disaggregation of legacy soil data using area to point kriging for mapping soil organic carbon at the regional scale. Geoderma, 7, 7 8. Law R, Illian J, Burslem D F R P, Gratzer G, Gunatilleke C V S, Gunatilleke I A U. 9. Ecological information from spatial patterns of plants: Insights from point process theory. Journal of Ecology, 97, Li R K, Zhu A X, Song X F, Li B L, Pei T, Qin C Z.. Effects of spatial aggregation of soil spatial information on watershed hydrological modelling. Hydrological Processes, 6, 9. McBratney A B, Mendonça Santos M L, Minasny B.. On digital soil mapping. Geoderma, 7,. Moran C J, Bui E.. Spatial data mining for enhanced soil map modelling. International Journal of Geographical Information Science, 6, 9. Odgers P, Sun W, Mcbratney A B, Minasny B, Clifford D.. Disaggregating and harmonising soil map units through resampled classification trees. Geoderma,, 9. Qi F.. Knowledge discovery from area-class resource maps: Sata preprocessing for noise reduction. Transactions in GIS, 8, Qi F, Zhu A X.. Knowledge discovery from soil maps using inductive learning. International Journal of Geographical Information Science, 7, Qi F, Zhu A X, Pei T, Qin C Z, Burt J E. 8. Knowledge discovery from area-class resource maps: Capturing prototype effects. Cartography and Geographic Information Science,, 7. Qi F, Zhu A X.. Comparing three methods for modeling the uncertainty in knowledge discovery from area-class soil maps. Computers & Geosciences, 7, 6. Qin C Z, Zhu A X, Pei T, Li B L, Scholten T, Behrens T, Zhou C H.. An approach to computing topographic wetness index based on maximum downslope gradient. Precision Agriculture,,. Quinn T, Zhu A X, Burt J E.. Effects of detailed soil spatial information on watershed modeling across different model scales. International Journal of Applied Earth Observation and Geoinformation, 7, 8. Silva S H G, de Menezes M D, Owens P R, Curi. 6. Retrieving pedologist s mental model from existing soil map and comparing data mining tools for refining a larger area map under similar environmental conditions in Southeastern Brazil. Geoderma, 67, Silverman B W Density Estimation for Statistics and Data Analysis. Chapman and Hall, London, UK. SSDS (Soil Survey Division Staff). 7. Soil Survey Manual. USDA Handbook 8. Government Printing Office, Washington, D. C. Subburayalu S K, Jenhani I, Slater B K.. Disaggregation of component soil series on an Ohio County soil survey map using possibilistic decision trees. Geoderma,,. Subburayalu S K, Slater B K.. Soil series mapping by knowledge discovery from an Ohio County soil map. Soil Science Society of America Journal, 77, 68. Vincent S, Lemercier B, Berthier L, Walter C. 6. Spatial disaggregation of complex Soil Map Units at the regional scale based on soil-landscape relationships. Geoderma,,. Yang L, Jiao Y, Fahmy S, Zhu A X, Hann S, Burt J E, Qi F.. Updating conventional soil maps through digital soil mapping. Soil Science Society of America Journal, 7,. Zhu A X A similarity model for representing soil spatial information. Geoderma, 77, 7. Zhu A X A personal construct-based knowledge acquisition process for natural resource mapping. International Journal of Geographical Information Science,, 9. Zhu A X, Band L E, Dutton B, imlos T Automated soil inference under fuzzy logic. Ecological Modelling, 9,. Zhu A X, Band L E, Vertessy R, Dutton B Derivation of soil properties using a Soil Land Inference Model (SoLIM). Soil Science Society of America Journal, 6,. Zhu A X, Hudson B, Burt J, Lubich K, Simonson D.. Soil mapping using GIS, expert knowledge, and fuzzy logic. Soil Science Society of America Journal, 6, 6 7. Zhu A X, Mackay D S.. Effects of spatial detail of soil information on watershed modelling. Journal of Hydrology, 8, 77. Zhu A X. 8. Model and Method of Detail Digital Soil Survey. Science Press, Beijing. (in Chinese) Zhu A X, Yang L, Li B L, Qin C Z, Pei T, Liu B Y.. Construction of membership functions for predictive soil mapping under fuzzy logic. Geoderma,, Section editor SHI Zhou Managing editor SU Lu-juan