Based on legacy soil data from a soil survey conducted recently in the traditional manner in Hong Kong of China, a digital soil mapping method was applied to produce soil order information for mountain areas of Hong Kong. Two modeling methods (decision tree analysis and linear discriminant analysis) were used, and their applications were compared. Much more eflort was put on selecting soil covariates for modeling. First, analysis of variance (ANOVA) was used to test the variance of terrain attributes between soil orders. Then, a stepwise procedure was used to select soil covariates for linear discriminant analysis, and a backward removing procedure was developed to select soil covariates for tree modeling. At the same time, ANOVA results, as well as our knowledge and experience on soil mapping, were also taken into account for selecting soil covariates for tree modeling. Two linear discriminant models and four tree models were established finally, and their prediction performances were validated using a multiple jackknifing approach. Results showed that the discriminant model built on ANOVA results performed best, followed by the discriminant model built by stepwise, the tree model built by the backward removing procedure, the tree model built according to knowledge and experience on soil mapping, and the tree model built automatically. The results highlighted the importance of selecting soil covariates in modeling for soil mapping, and suggested the usefulness of methods used in this study for selecting soil covariates. The best discriminant model was finally selected to map soil orders for this area, and validation results showed that thus produced soil order map had a high accuracy. Based on legacy soil data from a soil survey conducted recently in the traditional manner in Hong Kong of China, a digital soil mapping method was applied to produce soil order information for mountain areas of Hong Kong. Two modeling methods (decision tree analysis and linear discriminant first, analysis of variance (ANOVA) was used to test the variance of terrain attributes between soil orders. Then, a stepwise procedure was used to select soil covariates for linear discriminant analysis, and a backward removal procedure was developed to select soil covariates for tree modeling. At the same time, ANOVA results, as well as our knowledge and experience on soil mapping, were also taken into account for selecting soil covariates for tree modeling. Two linear discriminant models and four tree models were established finally, and their prediction performances were validated using a multiple jackknifing approach. Results that the discriminant model built on ANOVA results performed best, followed by the discriminant model built by stepwise, the tree model built by the backward filtering procedure, the tree model built according to knowledge and experience on soil mapping, and the tree model built automatically. The results highlighted the importance of selecting soil covariates in modeling for soil mapping, and suggested the usefulness of of used used in this study for selecting soil covariates. The best discriminant model was finally selected to map soil orders for this area, and validation results showed that that produced produced soil order map had a high accuracy.