The Scientific World Journal

Volume 2014 (2014), Article ID 109310, 11 pages

http://dx.doi.org/10.1155/2014/109310

Research Article

Particle-Size Distribution Models for the Conversion of Chinese Data to FAO/USDA System

Wei Shangguan,¹ YongJiu Dai,¹ Carlos García-Gutiérrez,² and Hua Yuan¹

¹College of Global Change and Earth System Science, Beijing Normal University, No. 19, Xinjiekouwai Street, Beijing 100875, China

²Department of Applied Mathematics, School of Agricultural Engineering Technical, University Polytechnic of Madrid, 28040 Madrid, Spain

Received 9 April 2014; Accepted 23 May 2014; Published 7 July 2014

Academic Editor: Cezary Sławiński

Copyright © 2014 Wei Shangguan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigated eleven particle-size distribution (PSD) models to determine the appropriate models for describing the PSDs of 16349 Chinese soil samples. These data are based on three soil texture classification schemes, including one ISSS (International Society of Soil Science) scheme with four data points and two Katschinski’s schemes with five and six data points, respectively. The adjusted coefficient of determination , Akaike’s information criterion (AIC), and geometric mean error ratio (GMER) were used to evaluate the model performance. The soil data were converted to the USDA (United States Department of Agriculture) standard using PSD models and the fractal concept. The performance of PSD models was affected by soil texture and classification of fraction schemes. The performance of PSD models also varied with clay content of soils. The Anderson, Fredlund, modified logistic growth, Skaggs, and Weilbull models were the best.

1. Introduction

Particle-size distribution (PSD) is a basic physical property of soils that affects many important soil properties. PSD has been widely used for estimating various soil hydraulic properties [1, 2]. The PSD prediction has been used for comparing and converting texture measurements from different classification systems [3–6]. Many soil databases do not contain detailed PSD data but contain only several mass fractions and use different cutoff points to separate the size fractions. To obtain a more complete description of texture, various parametric PSD models are used [7–12]. In the Second National Soil Surveys of China, soil textures were measured by ISSS (International Society of Soil Science) and Katschinski’s systems. The conversion from ISSS and Katschinski’s to the popular USDA (United States Department of Agriculture) system is required to achieve compatibility of soil data and to use pedotransfer functions as estimators of soil hydraulic properties [13, 14].

PSD models can be categorized as several classes: regression [15], spline [13], fractal method [16, 17], method based on PSD statistics [12], similarity procedure [13], gray model [18], and so on. In this study, we investigated regression models for their simplicity and effectiveness when there are a small number of fractions.

A few studies have been performed with the purpose of determining the best model for fitting particle-size distribution curves of soils [9, 15, 19–23]. Hwang et al. [15] compared seven PSD models to fit PSD data sets of Korean soils and found that the Fredlund model with four fitting parameters [11] showed the best performance for the majority of texture types. Hwang [19] examined nine PSD models and found that the performance of the PSD models was affected by soil texture and most of them improved with an increase in clay content. He also reported that the Fredlund model and two Skaggs models [12] showed better performance for several specific soil textures. Liu et al. [20] evaluated the suitability of five models to fit particle-size analysis data of Chinese soils of three texture types. The results indicated that the four-parameter Fredlund model showed the best representation of particle-size distributions, and the three-parameter Fredlund model and modified logistic growth model produced comparative results for silty clay loam and silt loam soils but yielded worse results for sandy loam soils. da Silva et al. [21] tested 14 different models with feasibility to fit the cumulative soil particle-size distribution curve based on four measured points and experienced that the most recommendable models were the Skaggs et al. [12], Lima and da Silva [24], Weibull [25], and Morgan et al. [26] models, which are all three-parameter models. Botula et al. [23] evaluated ten PSD models for soils of humid tropics and found that the performance was dependent on the bimodal character of the soil PSD. However, these studies have not fully investigated the performance of PSD models and used data measured only according to USDA standard.

Several studies have proposed conversion methods between texture classification systems. Nemes et al. [13] evaluated the accuracy of four interpolation procedures using a data set of European countries, and these procedures were used to accomplish the standardization of particle-size description by the FAO (Food and Agriculture Organization)/USDA system. Minasny and McBratney [4] used an empirical multiple linear regression to accomplish the conversion between the ISSS and USDA system. Shirazi et al. [27] developed a conversion table between the USDA and ISSS system based on a lognormal particle-size distribution assumption within each size fraction (sand, silt, and clay). Shirazi et al. [5] expanded the development of the PSD statistics (i.e., the geometric mean particle diameter and its standard deviation) to compare USDA and ISSS systems including rock fragments, which can be applied to the conversion between the two systems. Most PSD models are incapable of supporting the conversion from Katschinski’s to USDA systems. These PSD models are not capable of extrapolating the soil fraction above the upper limit (i.e., 1 mm) of Katschinski’s system, but the fraction between 2 mm and 1 mm is included by the USDA system. Few studies have been devoted to the conversion between USDA and Katschinski’s systems [3]. The conventional method is plotting the particle-size curve by semilogarithm regression [28], but it is time-consuming and of low accuracy. Rousseva [3] used the fraction between 1 mm and 3 mm to interpolate the cumulative fraction at 2 mm limit by exponential and power-law distribution models during the transformation from Katchinski’s to ISSS and USDA systems. However, the gravel contents are usually not measured (e.g., the Second National Soil Surveys of China). The following question arises: is it possible to extrapolate the soil fraction between 2 mm and 1 mm based on the data below 1 mm limit? The fractal method [29] provides an optional solution, which was used in this work.

The objective of this study is to evaluate the performances of several PSD models to find the best describing the soil PSD, to convert the data from ISSS and Katschinski’s schemes to the USDA standard, and to investigate the effects of soil texture and classification schemes on the efficiency of PSD models.

2. Material and Methods

2.1. The Chinese Soil Profile Database

The Chinese soil profile database was established using the results of the Second National Soil Surveys of China conducted in the 1980s. It contains data from 33010 soil horizons representing 8979 profiles. The data have been published by the National Soil Survey Office [30], soil survey offices of provinces, and some Tibetan counties’ soil survey offices. This work collected and digitized these data, which are available in Microsoft Excel worksheet file format. Not all soil horizons have PSD data. Fine size fractions were determined by the hydrometer or pipette method, whereas the coarse fractions were obtained by sieving [28]. Particle-size fraction data were classified by ISSS or Katschinski’s schemes. Katschinski’s [31] scheme was not only used in China but also widely used in countries of former USSR and Eastern Europe. Because some schemes have small number of soil samples and some schemes have too few particle-size limits (e.g., 2, 0.02, and 0.002) for PSD models to fit the curve, we used only part of the data (i.e., 16349 samples) in this study. Four to six particle-size fractions are given in one ISSS scheme and in two Katschinski’s schemes (Table 1).

Table 1: Soil particle-size data used in this study.

Soil fractions schemes
The particle-size limits (mm)
Number

T1 (ISSS)
2, 0.2, 0.02, 0.002
15262

T2 (Katschinski)
1, 0.25, 0.05, 0.01, 0.005, 0.001
671

T3 (Katschinski)
1, 0.05, 0.01, 0.005, 0.001
1090

ISSS: international system of International Soil Science Society.

Katschinski: Katschinski classification system of Russia.

T3 contains the data of T2 after summing up the 1–0.25 mm and 0.25–0.05 mm fraction.

2.2. Particle-Size Distribution Models

Eleven PSD models were tested in this study (Table 2). Ten models are shown in Table 2 and a self-similar model (SELF) is described here because this model cannot be given in a simple equation [32] and it has not yet been compared with other models in previous studies.

Table 2: Particle-size distribution models.

Name
Model
Parameters

Anderson

(AD [8])

b, c, ,

Fredlund4P

(F4P [11])

, , ,  

()

Fredlund3P

(F3P [11])

, ,  

()

Modified logistic growth

(ML [43])

a, b, c

Offset-nonrenormalized

lognormal

(ONL [9])
, 

where  

, ,

Offset-renormalized lognormal

(ORL [9])

defined in ONL model]
, ,

Skaggs

(S [12])
, 

where  

= 0.002 mm for T1; = 0.001 mm for T2, 

T3; : fraction <
u, c

van Genuchten type

(VG [7])

where
, n

van Genuchten type modified

(VGM [34])
, 

where ; 

, fraction of minimum particle size
a, n

Weibull

(W [10])
, 

where  

, for T1;  

, for T2, T3
a, b, c,

d: particle diameter in mm.

: error function.

The self-similar model uses the soil textural data to define an iterated function system, which determines how the self-similar distribution reproduces its fractal structure at different length scales [33]. A set of proportions of soil mass are selected (for simplification, assume that ). The fractal model for PSD derived from the self-similarity hypothesis is constructed as follows (see Martín and Taguas [32] for the details): let , , and be the subfractions of sizes corresponding to the three size classes and let , , and be the proportions of mass for the fractions , , and , respectively. Notice that and . Associated with these definitions, we consider the following functions: , , and , where , , and and is any point (or particle size value) in the fraction . , , and are the linear mappings (similarities) that transform the points of the fraction in the points of the fractions , , and , respectively. The set , is called an iterated function system (IFS) and it defines a self-similar distribution μ supported on , which satisfies for every . The set of textural data, together with the self-similarity hypothesis mentioned above, determine a self-similar mass distribution for particle size determined in , which fits exactly the input data; that is, . From the self-similarity, the contribution to soil mass of particles with sizes in the fraction is given by . In general, for any size fraction , the IFS model allows one to determine the mass of soil corresponding to the soil particles with sizes in the fraction . In this study, the values of the fraction limits chosen are , , and (in μm); and the model is seen as a biparametric model ( and ) which varies until the best simulation, in terms of minimizing the distance between simulated and real data, is reached.

Six models which showed relatively good performance in the work of Hwang et al. [15, 19] were tested, including the Fredlund4P (F4P), Skaggs (S), van Genuchten (VG), Weibull (W), offset-nonrenormalized lognormal (ONL), and offset-renormalized lognormal (ORL) models. In addition, four other models with different number of fitting parameters were tested. These models are selected among scores of PSD models in the literature according to their performance in previous studies. The Anderson model (AD) [8] assumes that the log mass of the particles is arcus-tangent distributed (i.e., Cauchy distributed). Fredlund et al. [11] found that a value of 0.001 for the parameter in F4P model (Table 2) provided a reasonable fit in most cases, so the Fredlund3P (F3P) model with three parameters was also tested. The modified logistic growth model (ML) [20] is modified from the logistic growth model whose graph is in a sigmoid shape. The van Genuchten type model is modified by Zhuang et al. [34] (VGM) which is also based on an assumption that the shape of particle-size distribution is sigmoidal. Spline is not included because Liu et al. [20] found that it performs worse than modified logistic model.

2.3. Model Comparison and Fitting Techniques

We used three statistical indices to determine the performance of the PSD models. The adjusted coefficient of determination () was used as a relative measure of goodness-of-fit between predicted data and observed data, which is a better criterion than the coefficient of determination [23]. is defined as

with

where is the number of observed data points, is the number of model parameters, and and are the observed and predicted cumulative mass fractions, respectively. A further index was Akaike’s information criterion (AIC) [35] which examines the complexity of the model together with its goodness of fit to the sample data, balances between the two, and discourages overfitting. AIC is defined as

with

where RSS is the residual sum of squares. Geometric mean error ratio (GMER) was used to test whether the underestimation (values bigger than 1) or overestimation (values smaller than 1) in model’s prediction occurs:

The “soiltexture” package of language was developed to fit the parametric models to the raw data [36]. The optimization methods in the “optim” function of are used, including Nelder-Mead, quasi-Newton and conjugate-gradient, box-constrained, and simulated annealing algorithms [37]. We ran all methods and chose the best results with minimum RSS. The nonlinear optimization procedures were carried out using at least five random initial parameter estimates for all soils. When the final solution for each soil converged to different parameter values, the parameter values with the best fitting statistics (i.e., RSS) were kept. In most cases, the parameters were similar.

2.4. Conversion to the USDA Standard

The soil data were converted to the USDA standard, including soil fractions (i.e., clay fraction ( mm), silt fraction ( mm), and sand fraction ( mm)).

After estimating the parameters of each model for every soil sample, we can easily obtain the predicted values of the cumulative mass fractions at a specific point between the upper and lower limits of the three soil fraction schemes, so the fractions in USDA standard can be obtained. The result of the model with the biggest value was adopted for a sample in the conversion to the USDA standard. For T1 scheme, the clay fraction was already known, the silt fraction was the predicted cumulative fraction at the observed 0.05 mm limit minus clay fraction, and the sand fraction was calculated by subtracting silt fraction and clay fraction from 100%. For T2 and T3 scheme, the clay fraction was predicted by PSD models, and the silt fraction was calculated by subtracting clay fraction from the cumulative fraction at the observed 0.05 mm limit. To predict the cumulative mass fraction at the 2 mm limit, the fractal method was adopted. In recent years, PSD of soil samples has been successfully studied by fractal methods. Turcotte [16] and Tyler and Wheatcraft [17] developed a relationship relating mass measurements and diameters for the analysis of PSD, which is expressed as

where is the mass of soil particles with a radius smaller than a prespecified size ; is the total mass of particles; is the upper particle size limit; and is the fractal dimension of fragmentation. The fractal dimension of fragmentation can be found by estimating the slope coefficient of the versus plot. Many studies using detailed experimental data have shown multiple fractal dimensions on log-transformed PSD of soil samples. Wu et al. [38] found three domains within PSDs determined over six orders of magnitude in the particle size. Bittelli et al. [39] reported that three main power-law domains could characterize the PSD across the whole range of measurements, including the sand domain between and 2 mm. Millán et al. [40] divided the whole PSD domain into two fractal domains (i.e.,  mm and  mm), where the  mm domain spanned approximately one order of magnitude with fractal concepts. Prosperini and Perugini [41, 42] found that all soils in their analysis displayed two scaling domains in the range from 0.0014 to 50.8 mm. According to these studies, the slope of the log-log plot between the cumulative fractions and particle sizes has a linear trend in the domain from 0.5 to 2 mm. The relationship between the cumulative fraction at 2 mm and the cumulative fractions at 1 and 0.5 mm can be derived from (6). The equation is given by

where , , and are cumulative fractions with a diameter smaller than 2, 1, and 0.5 mm, respectively. and were predicted by PSD models. All three fractions were divided by the predicted value at 2 mm so that the cumulative final value was 100 percent. We did not apply the fractal model alone to do the conversion, because there are not enough observed data points to determine the subdomains and fractal dimensions to accomplish the interpolation for each soil sample in this study.

3. Results and Discussion

3.1. Model Performance

Among all of the three soil fraction schemes and all of the models, values of ranged from 0.873 to 1 (Figure 1). In T1 scheme of ISSS, the S and W models showed the highest values; F3P, ML, and ORL models gave similar performance; and the other four models (ONL, VG, VGM, and SELF) had lower values. In the two Katschinski’s schemes (i.e., T2 and T3), the performances of PSD models were quite similar to each other. The S model with two parameters yielded the highest values. The AD, F3P, F4P, ML, and W models indicated similar performance. The other five models had lower values. The F4P model with four parameters performed slightly better than the F3P model with three parameters in both Katschinski’s schemes, which indicates that F3P is enough for most soils with the advantage of one less parameter [11]. The ONL model had lower values in the T1 scheme of ISSS than in Katschinski’s scheme, while the ORL had better performance in the T1 scheme. The performance of AD model became better from the T2 scheme to the T3 scheme as one limit (0.25 mm) of PSD is omitted in T3, while the opposite happened to the SELF model. The results of the AIC test were generally consistent with those of assessment (Figures 1 and 2). The performance of AD and F4P model with four parameters in the T1 scheme can be evaluated by AIC but not by , because the denominator (i.e., -) will be zero in (3). The AD model had the lowest AIC, while the F4P had similar AIC with the F3P model. Though AIC has penalty on parameter number, the AD model with four parameters and the W model with three parameters were still among the best PSD models.

(a)

(b)

(c)

Figure 1: Box plot for percentiles as the goodness-to-fit of 11 models, (a) T1 scheme soils, (b) T2 scheme soils, and (c) T3 scheme soils. AD = Anderson model, F3P = Fredlund3P model, F4P = Fredlund4P model, ML = modified logistic growth model, ONL = offset-nonrenormalized lognormal model, ORL = offset-renormalized lognormal model, S = Skaggs model, VG = van Genuchten type model, VGM = van Genuchten type modified model, W = Weibull model, and SELF = self-similar model. There is no box plot for AD and F4P model because cannot be calculated for them.

(a)

(b)

(c)

Figure 2: Box plot for Akaike’s information criterion (AIC) percentiles as a criterion for assessing the performance of 11 models, (a) T1 scheme soils, (b) T2 scheme soils, and (c) T3 scheme soils.

Though the performances evaluated by and AIC were quite similar on the whole (Figures 1 and 2), the best model according to different criteria can be different for a specific soil. Figure 3 shows the percentage of cases where a model was the best according to and AIC for soils of the T3 scheme. The percentage of best cases according to for the AD and F4P models was smaller than that according to AIC, while the opposite happened to the S and ONL models.

Figure 3: Percentage of soils of T3 scheme for which or AIC are the best for a given model.

Figure 4 shows GMER for all of the three soil fraction schemes and all of the models. Most models overestimated the fractions while the S and VG model underestimated them in general. The VGM model showed an overestimation in T1 scheme but it showed an underestimation in T2 and T3, while the opposite happened to the ONL model. The AD and W models showed the best results without apparent underestimation or overestimation.

(a)

(b)

(c)

Figure 4: Box plot for geometric mean error ratio (GMER) percentiles as a criterion for assessing the performance of 11 models, (a) T1 scheme soils, (b) T2 scheme soils, and (c) T3 scheme soils.

Over all, the AD, S, and W models provided the best performance for soils in ISSS and Katschinski’s schemes according to the analysis, AIC tests, and GMER. The F4P, F3P, and ML models also displayed relative good performances.

3.2. The Validation of Combined Use of PSD Model and Fractal Method

PSD models to interpolate and fractal method to extrapolate have different statistical assumptions. We need to know whether the combined use of two distinct models results in a real representing soil sample. For this, 20 soil samples were used to validate the combined use, which has the observed cumulative particle-size percentages at 2 mm, 0.05 mm, and 0.002 mm of USDA standard and at the limits of T2 scheme, respectively. These samples were only available with the given specifics. We used each PSD model combined with the fractal method to predict PSD with data points of T2 scheme. Figure 5 illustrates the test of suitability of the combined method to predict percentages of particles finer than 0.05 and 0.002 mm. The value was 0.960, which was slightly lower than those of PSD models. Thus, the combined method with interpolation and extrapolation is suitable for transferring data from Katschinski’s to USDA scheme. However, it cannot be determined whether this method is suitable for a specific soil textural class, because the sample size is too small.

Figure 5: Percentages predicted by one of the eleven PSD models and fractal method versus measured percentages of particles finer than 0.05 mm and 0.002 mm of 20 soils.

3.3. Effect of Soil Textural Classes and Fraction Schemes

The AIC analysis had revealed that the performance of a PSD model can be affected by soil textural classes and soil fraction schemes (Tables 3 and 4). The dataset used in this study represents all the 12 textural classes of ISSS system and 10 textural classes of Katschinski system but in various degrees (Tables 3 and 4). The dominant textural classes in the ISSS scheme are light clay (25.7%), clay loam (21.7%), sandy loam (15%), and loam (9%). The dominant textural classes in the T2 scheme are light clay (40.4%), moderate clay (27.6%), heavy loam (14.8%), and heavy clay (12%). The dominant textural classes in the T3 scheme are light clay (33.6%), moderate clay (21%), heavy loam (19%), and moderate loam (11.3%). Table 3 shows the number of cases as the best model due to the smallest AIC value for each soil textural class of ISSS in the T1 scheme. All the models performed best in some cases. The AD model had the largest number of cases for most texture, while the S model did so for sand and silt clay. The AD model had the largest number of cases as the best model, followed by the S model. Although the S and W models performed similarly in the and AIC analysis (Figures 1 and 2), the S model performed better with a bigger number of cases which resulted in the smallest AIC values for all soil textural classes. Although the SELF model did not perform well in the overall analysis (Figures 1 and 2), it had similar number of best cases with the W model. These results suggest that the PSD models that show good performance for overall soil PSDs are not always good models for each soil textural class. Table 4 shows the number of cases of all the tested PSD models having the smallest AIC value for each soil textural class in the T2 and T3 schemes. In the T2 scheme (Table 4(a)), most soils are clay soils in Katschinski system (i.e., heavy clay, moderate clay, and light clay). The AD model had the largest number of cases which resulted in the smallest AIC values for these clay soils. However, the S model performed better for heavy loam. In the T3 scheme (Table 4(b)), the AD model also had the largest number of cases having the smallest AIC values for almost all soil textures, though the F4P, S, and W models also have a large number of cases. From the comparisons among the T1, T2, and T3 schemes, it appears that the performance of PSD models vary not only with the soil textural classes but also with the soil fraction schemes. However, the dataset used in this work does not seem to be sufficient to draw a well-founded conclusion for some textural classes with a small number of samples, especially in the T2 and T3 schemes.

Table 3: Number of cases as the best model with the smallest AIC value for each soil textural class in T1 schemes.

Texture
AD
F3P
F4P
ML
ONL
ORL
S
VG
VGM
W
SELF
Sum

Sa
64
33
17
9
0
32
73
0
0
45
47
320

LoSa
163
39
53
46
3
37
153
0
0
78
54
626

Lo
799
16
3
24
1
70
380
0
1
36
47
1377

SaLo
979
56
55
36
22
35
589
0
4
189
312
2277

SaClLo
464
25
56
23
7
12
229
0
2
31
97
946

SiClLo
244
22
1
6
0
14
196
1
0
27
0
511

SiLo
137
9
0
4
0
17
95
0
0
24
0
286

ClLo
1929
31
0
62
9
72
1024
0
2
80
110
3319

SaCl
108
4
16
3
0
1
36
0
0
7
6
181

SiCl
236
16
1
55
2
18
264
3
0
39
0
634

LCl
2086
33
34
138
1
66
1346
0
1
97
115
3917

HCl
385
10
21
39
1
23
335
0
0
47
7
868

Sum
7594
294
257
445
46
397
4720
4
10
700
795
15262

For example, if the best model according to AIC value of a soil sample is the AD model, the number of the AD model in the table will increase by one. If a model has a big number in a soil textural class, it has advantages in fitting the PSD curve compared to models with small number. The compared eleven models are AD = Anderson model, F3P = Fredlund3P model, F4P = Fredlund4P model, ML = modified logistic growth model, ONL = offset-nonrenormalized lognormal model, ORL = offset-renormalized lognormal model, S = Skaggs model, VG = van Genuchten type model, VGM = van Genuchten type modified model, W = Weibull model, and SELF = self-similar model.

Textural classes of ISSS system: Sa = sand, LoSa = loamy sand, Lo = loam, SaLo = sandy loam, SaClLo = sandy clay loam, SiClLo = silty clay loam, SiLo = silt loam, ClLo = clay loam, SaCl = sandy clay, SiCl = silty clay, LCl = light clay, and HCl = heavy clay.

The biggest number of cases as the best model for the specific soil texture class.

Table 4: Number of cases as the best model with the smallest AIC value for each soil textural class in T2 (a) and T3 (b) schemes.

(a)

Texture
AD
F3P
F4P
ML
ONL
ORL
S
VG
VGM
W
SELF
Sum

SL
1
1
0
0
0
0
0
0
0
0
0
2

HL
17
5
7
4
8
2
32
5
5
14
0
99

ML
8
4
1
4
1
6
3
0
0
1
0
28

LL
2
1
0
0
0
1
1
1
0
0
0
6

HC
38
6
19
1
1
0
4
0
0
3
8
80

MC
62
8
50
5
9
1
31
3
1
10
5
185

LC
79
0
60
7
12
2
67
8
6
29
1
271

Sum
207
25
137
21
31
12
138
17
12
57
14
671

(b)

Texture
AD
F3P
F4P
ML
ONL
ORL
S
VG
VGM
W
SELF
Sum

TS
2
1
1
0
0
0
0
0
0
0
0
4

LS
0
1
2
0
0
0
0
0
0
0
0
3

SL
9
2
1
0
1
2
5
0
0
0
0
20

HL
74
1
16
10
22
7
45
0
3
23
6
207

ML
53
7
9
4
6
15
14
0
1
13
1
123

LL
21
4
7
1
0
5
4
1
0
4
0
47

HC
40
7
14
0
0
0
12
0
0
7
11
91

MC
90
1
37
2
9
7
44
0
0
30
9
229

LC
143
5
39
0
39
3
85
1
4
39
8
366

Sum
432
29
126
17
77
39
209
2
8
116
35
1090

For example, if the best model according to AIC value of a soil sample is the AD model, the number of the AD model in the table will increase by one. If a model has a big number in a soil textural class, it has advantages in fitting the PSD curve compared to models with small number.

Textural classes of Katschinski system: TS = tight sand, LS = loose sand, SL = sandy loam, HL = heavy loam, ML = moderate loam, LL = light loam, HC = heavy clay, MC = moderate clay, and LC = light clay.

The biggest number of cases as the best model for the specific soil texture class.

The performance of PSD models appears to be affected by clay content of soils in different manners. In the T2 scheme (Figure 6), the clay content is defined as particles smaller than 0.001 mm in Katschinski system. All models performed best with the clay content between 30% and 40%; the AD model had the smallest values with clay content between 10% and 20%, and other models had the smallest values with clay content larger than 40%; the S and W model became better with the increase in the clay content but were very poor with the clay content over 40%. However, the performance of these models with different clay content in the T1 and T3 schemes is quite different from that in the T2 scheme (results not shown here). Hwang et al. [15] found better fitting with the increase in the clay content. Botula et al. [23] observed no clear trend of performance with clay content using a more clayed soil database. Our results generally confirmed the findings of Botula et al. [23].

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Box plot for percentiles of the PSD models on different clay content of soils in the T2 scheme.

4. Conclusions

Eleven models for soil PSD were compared, including the AD, F4P, F3P, ML, ONL, ORL, S, SELF, VG, VGM, W, and SELF models, using one ISSS scheme (i.e., T1) data with four data points and two Katschinski’s schemes (i.e., T2 and T3) with five and six data points. The results from values, AIC, and GMER test indicated that the AD, F4P, F3P, S, and W models had the best performance in three schemes. The combined use of PSD model and fractal method was suitable for transferring soil texture data from Katschinski’s to USDA scheme when there are only sparse observed data points. Soil texture can affect the performance of PSD models. And the performance of PSD models has some difference using experimental data of different soil fraction schemes. The performance of PSD models can be affected by clay content of soils in quite different manners for different classification of soil fraction schemes. This work contributes to the selection of soil PSD models to describe the soil particle-size distribution curve which can be used to estimate the hydraulic properties and to compare and convert the different soil texture classification systems.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was supported by the Chinese NSFC under Grant 41205037, MOST no. 2010CB951802, and the Fundamental Research Funds for the Central Universities. The data collection was helped by a group of students in Beijing Normal University. This research was also partially supported by Plan Nacional de Investigación Científica (I+D+i) under reference AGL2011/25175 and DGUI (Comunidad de Madrid) and UPM under reference QM100245066.

References

S. I. Hwang and S. P. Hong, “Estimating relative hydraulic conductivity from lognormally distributed particle-size data,” Geoderma, vol. 133, no. 3-4, pp. 421–430, 2006.

B. Minasny and A. B. McBratney, “Estimating the water retention shape parameter from sand and clay content,” Soil Science Society of America Journal, vol. 71, no. 4, pp. 1105–1110, 2007.

S. S. Rousseva, “Data transformations between soil texture schemes,” European Journal of Soil Science, vol. 48, no. 4, pp. 749–758, 1997.

B. Minasny and A. B. McBratney, “The Australian soil texture boomerang: a comparison of the Australian and USDA/FAO soil particle-size classification systems,” Australian Journal of Soil Research, vol. 39, no. 6, pp. 1443–1451, 2001.

M. A. Shirazi, L. Boersma, and C. Burch Johnson, “Particle-size distributions: comparing texture systems, adding rock, and predicting soil properties,” Soil Science Society of America Journal, vol. 65, no. 2, pp. 300–310, 2001.

W. Shangguan, Y. Dai, B. Liu et al., “A China data set of soil properties for land surface modeling,” Journal of Advances in Modeling Earth Systems, vol. 5, no. 2, pp. 212–224, 2013.

R. Haverkamp and J. Y. Parlange, “Predicting the water-retention curve from particle-size distribution: 1. Sandy soils without organic matter,” Soil Science, vol. 142, no. 6, pp. 325–339, 1986.

S. Andersson, “Markfysikaliska undersökningar I odlad jord, XXVI. Om mineraljordens och mullens rumsutfyllande egenskaper,” in En Teoretisk Studie, Swedish University of Agricultural Sciences, Uppsala, Sweden, 1990, (Swedish).

G. D. Buchan, K. S. Grewal, and A. B. Robson, “Improved models particle diameter and standard deviation from sand, silt, and clay of particle-size distribution: an illustration of model comparison fractions,” Soil Science, vol. 152, no. 4, pp. 427–431, 1993.

S. Assouline, D. Tessier, and A. Bruand, “A conceptual model of the soil water retention curve,” Water Resources Research, vol. 34, no. 2, pp. 223–231, 1998.

M. D. Fredlund, D. G. Fredlund, and G. W. Wilson, “An equation to represent grain-size distribution,” Canadian Geotechnical Journal, vol. 37, no. 4, pp. 817–827, 2000.

T. H. Skaggs, L. M. Arya, P. J. Shouse, and B. P. Mohanty, “Estimating particle-size distribution from limited soil texture data,” Soil Science Society of America Journal, vol. 65, no. 4, pp. 1038–1044, 2001.

A. Nemes, J. H. M. Wösten, A. Lilly, and J. H. Oude Voshaar, “Evaluation of different procedures to interpolate particle-size distributions to achieve compatibility within soil databases,” Geoderma, vol. 90, no. 3-4, pp. 187–202, 1999.

Y. Dai, W. Shangguan, Q. Duan, B. Liu, S. Fu, and G. Niu, “Development of a china dataset of soil hydraulic parameters using pedotransfer functions for land surface modeling,” Journal of Hydrometeorology, vol. 14, no. 3, pp. 869–887, 2013.

S. I. Hwang, K. P. Leeb, D. S. Leeb, and S. E. Powers, “Models for estimating soil particle-size distributions,” Soil Science Society of America Journal, vol. 66, no. 4, pp. 1143–1150, 2002.

D. L. Turcotte, “Fractals and fragmentation,” Journal of Geophysical Research: Solid Earth, vol. 91, no. 2, pp. 1921–1926, 1986.

S. W. Tyler and S. W. Wheatcraft, “Fractal scaling of soil particle-size distributions: analysis and limitations,” Soil Science Society of America Journal, vol. 56, no. 2, pp. 362–369, 1992.

C. Wu, J. Wen, and K. Chang, “Evaluation of the gray model GM(1,1) applied to soil particle distribution,” Soil Science Society of America Journal, vol. 73, no. 6, pp. 1775–1785, 2009.

S. I. Hwang, “Effect of texture on the performance of soil particle-size distribution models,” Geoderma, vol. 123, no. 3-4, pp. 363–371, 2004.

J. Liu, S. Xu, and H. Liu, “Investigation of different models to describe soil particle-size distribution data,” Advances in water Science, vol. 14, pp. 588–592, 2003 (Chinese).

E. M. da Silva, J. E. F. W. Lima, L. N. Rodrigues, and J. A. de Azevedo, “Comparison of mathematical models for fitting particle-size distribution curves,” Pesquisa Agropecuaria Brasileira, vol. 39, no. 4, pp. 363–370, 2004.

P. Zhao, M. Shao, and R. Horton, “Performance of soil particle-size distribution models for describing deposited soils adjacent to constructed dams in the China Loess Plateau,” Acta Geophysica, vol. 59, no. 1, pp. 124–138, 2011.

Y. Botula, W. M. Cornelis, G. Baert, P. Mafuka, and E. van Ranst, “Particle size distribution models for soils of the humid tropics,” Journal of Soils and Sediments, vol. 13, no. 4, pp. 686–698, 2013.

J. E. F. W. Lima and E. M. da Silva, “Utilização do modelo modificado de Genuchten para o traçado da curva granulométrica,” in Encontro Nacional de Engenharia de Sedimentos, vol. 5, pp. 121–125, ABRH, São Paulo, Brazil, 2002.

W. A. Weibull, “A statistical distribution function of wide applicability,” Journal of Applied Mechanics, vol. 11, pp. 292–297, 1951.

P. H. Morgan, L. P. Mercer, and N. W. Flodin, “General model for nutritional responses of higher organisms,” Proceedings of the National Academy of Sciences of the United States of America, vol. 72, no. 11, pp. 4327–4331, 1975.

M. A. Shirazi, L. Boersma, and J. W. Hart, “A unifying quantitative analysis of soil texture: improvement of precision and extension of scale,” Soil Science Society of America Journal, vol. 52, no. 1, pp. 181–190, 1988.

National Soil Survey Office, Soil Survey Technology of China, Agriculture Press, Beijing, China, 1992.

Z. Jin, Y. S. Dong, Y. C. Qi, W. G. Liu, and Z. S. An, “Characterizing variations in soil particle-size distribution along a grass-desert shrub transition in the ordos plateau of Inner Mongolia, China,” Land Degradation and Development, vol. 24, no. 2, pp. 141–146, 2013.

National Soil Survey Office, China Soil Species Description, vol. 1–6, China Agriculture Press, Beijing, China, 1996, (Chinese).

N. A. Katschinski, “Die mechanische Bodenanalyse und die Klassifikation der Böden nach ihrer mechanischen Zusammensetzung,” in Rapports au Sixiéme Congrés International de la Science du Sol, pp. 321–327, Paris, France, 1956.

M. A. Martín and F. J. Taguas, “Fractal modelling, characterization and simulation of particle-size distributions in soil,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 454, no. 1973, pp. 1457–1468, 1998.

K. Falconer, Fractal Geometry, Mathematical Foundations and Applications, John Wiley & Sons, Chichester, UK, 1990.

J. Zhuang, J. Yan, and T. Miyazaki, “Estimating water retention characteristic from soil particle-size distribution using a non-similar media concept,” Soil Science, vol. 166, no. 5, pp. 308–321, 2001.

H. Akaike, “Information theory and an extension of the maximum likelihood principle,” in Proceedings of the 2nd International Symposium on Information Theory (TSAhkadsor, 1971), B.N. Petrov and F. Csaki, Eds., pp. 267–281, Akademiai Kiado, Budapest, Hungary, 1973.

J. Moeys and W. Shangguan, “Package—Soil texture: Functions for soil texture plot, classification and transformation (Version: 1.2.13),” 2014, http://cran.r-project.org/web/packages/soiltexture/

R Core Team, R: A Language and Environment for Statistical Computing, Foundation for Statistical Computing, Vienna, Austria, 2014, http://www.r-project.org/.

Q. Wu, M. Borkovec, and H. Sticher, “On particle-size distributions in soils,” Soil Science Society of America Journal, vol. 57, no. 4, pp. 883–890, 1993.

M. Bittelli, G. S. Campbell, and M. Flury, “Characterization of particle-size distribution in soils with a fragmentation model,” Soil Science Society of America Journal, vol. 63, no. 4, pp. 782–788, 1999.

H. Millán, M. González-Posada, M. Aguilar, J. Domínguez, and L. Céspedes, “On the fractal scaling of soil data. Particle-size distributions,” Geoderma, vol. 117, no. 1-2, pp. 117–128, 2003.

N. Prosperini and D. Perugini, “Application of a cellular automata model to the study of soil particle size distributions,” Physica A: Statistical Mechanics and its Applications, vol. 383, no. 2, pp. 595–602, 2007.

N. Prosperini and D. Perugini, “Particle size distributions of some soils from the Umbria Region (Italy): fractal analysis and numerical modelling,” Geoderma, vol. 145, no. 3-4, pp. 185–195, 2008.

J. Liu, S. Xu, H. Liu, and F. Guo, “Application of parametric models to description of particle-size distribution in loamy soils,” Acta Pedologica Sinica, vol. 41, pp. 375–379, 2004 (Chinese).