# County-Scale Spatial Variability of Macronutrient Availability Ratios in Paddy Soils

Authors: Mingkai Qu; Weidong Li; Chuanrong Zhang
Journal: Applied and Environmental Soil Science
 (2014)

---

## Abstract

Macronutrients (N, P, and K) are essential to plants but also can be harmful to the environment when their available concentrations in soil are excessive. Availability ratios (available concentration/total concentration) of macronutrients may reflect their transforming potential between fixed and available forms in soil. Understanding their spatial distributions and impact factors can be, therefore, helpful to applying specific measures to modify the availability of macronutrients for agricultural and environmental management purposes. In this study, 636 topsoil samples (0–15 cm) were collected from paddy fields in Shayang County, Central China, for measuring soil properties. Factors influencing macronutrient availability ratios were investigated, and total and available concentrations of macronutrients were mapped using geostatistical method. Spatial distribution maps of macronutrient availability ratios were further derived. Results show that (1) availability of macronutrients is controlled by multiple factors, and (2) macronutrient availability ratios are spatially varied and may not always have spatial patterns identical to those of their corresponding total and available concentrations. These results are more useful than traditional soil macronutrient average content data for guiding site-specific field management for agricultural production and environmental protection.


---

## Body

## 1. Introduction

Soil macronutrients (i.e., Nitrogen (N), phosphorus (P), and potassium (K)) are essential to plants [1]. They provide nutrients necessary for plant growth, which are important to maintain ecosystems and high crop yields. However, macronutrients, particularly N, and P can be potentially hazardous to water resources when their available components in soils are excessive, because available macronutrients can be transported off site in runoff due to rain or irrigation [2–4]. Improper or excessive fertilization has increasingly become a serious problem and the eutrophication problem caused by the losses of N and P from farmland to water bodies has caught people’s attentions [5–7]. Therefore, proper management of soil N, P, and K is necessary to avoid deteriorating the environment while meeting the requirement of high crop productivity. In addition, reducing the losses of macronutrients from farmland also can save the costs spent on fertilizers.Most of macronutrient contents exist in fixed forms in soils (e.g., contained in organic matter and minerals) and thus cannot be directly utilized by plants or transported to water bodies. Part of fertilizers applied to soil also can be fixed by soil and thus become unavailable to plants. This means that the total content of a macronutrient in soil is only a potentially available content in a long term, rather than its currently available content. Apparently the total content and the available content of a macronutrient are two different measures for the macronutrient in soil, and it is the availability ratio (i.e., available concentration/total concentration) that represents the potential effectiveness of a specific macronutrient in soil. That is to say, all the three indices may be necessary to understand the general situation of a macronutrient in soil. It is, therefore, important to investigate the spatial variability of availability ratios of soil macronutrients and corresponding controlling factors so that proper measures may be taken to modify the availability of the macronutrients for site-specific management.In addition to macronutrients, micronutrients and trace metals in soil also have similar problems, and their availability ratios in soil are also important indices for agricultural and environmental management. Recently, Wu et al. [8] studied the variability of Cu availability in soil and its relation to some soil properties in a 731 km2 region in southeast China. Their results showed that Cu availability ratio has a distributional trend similar to those of total Cu concentration and available Cu concentration, which are similar to each other. Zhang et al. [9] modeled the availability ratio of Cu in soil in Fuyang County, China, using the general regression neural network method. Although many studies on the spatial variability of soil N, P, and K were conducted in last several decades [10–20], most of them, however, were focused on their total contents or available contents in agricultural fields, and studies on the spatial distributions of their availability ratios in soils and corresponding impact factors at regional scales were rarely seen. Therefore, a study on the spatial variability of availability ratios of soil macronutrients and their impact factors in soil at a regional scale should be interesting and may provide valuable information to soil scientists and agricultural managers.The objectives of this study were to(
1
) analyze the relationships between some soil properties and the availability of soil N, P, and K, (
2
) determine the controlling factors on the availability ratios of macronutrients, (
3
) map the spatial distributions of the indices of macronutrients in paddy soils at a regional scale, and (
4
) understand the characteristics of spatial patterns of macronutrient availability ratios. Shayang County in central China was selected as our case study area. This paper will provide an example for the availability ratio study of soil macronutrients, and results of such kind of studies will be helpful to decision making in applying concrete measures to modify the availability of soil macronutrients for site-specific agricultural management purposes.


## 2. Materials and Methods

### 2.1. Study Area

Shayang County, Hubei Province in central China was selected as the study area. The county is among the most developed regions for cropping in Hubei Province and is located in the Jianghan Plain, bounded by the longitudes 112°03′ and 112°43′ east, and the latitudes 30°23′ and 30°57′ north, with an area of 2044 km2. It belongs to the northern subtropical monsoonal climate zone, with a temperate-humid climate throughout the year and four distinct seasons. The average annual temperature is 16.1°C and the mean annual precipitation is 1025.6 mm. In this region, paddy fields constitute the dominant land use type of arable lands.


### 2.2. Soil Sampling and Chemical Analysis

636 topsoil samples were collected from paddy fields in 2007. The sampling scheme considered the distributional uniformity of samples and also ensured that all samples were located in paddy fields (Figure1). All samples were taken in fall after harvest and before next cropping season to avoid the effect of fertilization during crop cultivation. For each soil sample, soils at 6–8 points within a small area of approximately 0.01 ha were collected from surface layers (0–15 cm) and then mixed and delivered to a laboratory for analysis. Exact sample locations were recorded using a hand-held global position system. All samples were air-dried at room temperature (20–22°C), crushed, and sieved into soil particles less than 2 mm. Portions of each sample (about 100 g) were further ground and sieved through 0.149 mm mesh. The prepared soil samples were stored for chemical analysis.Figure 1
The study area and sampling sites.Soil organic matter (SOM) was measured using the wet oxidation method at 180°C with a mixture of potassium dichromate and sulfuric acid [21]. Soil cation exchange capacity (CEC) was measured using 1.0 mol/L ammonium acetate solution [21]. Soil pH was measured using a pH meter (Sartorius Basic pH meter PB-10) with a soil/water ratio of 1 : 2.5 [21]. Total nitrogen (TN) was determined by the Kjeldahl method with H2SO4 + H2O2 digestion [22]. Available nitrogen (AN) was extracted using 2.0 mol/L KCl [22]. Total phosphorus (TP) was extracted using melt sodium hydroxide and measured using the molybdenum stibium anticolor method [22]. Available phosphorus (AP) was extracted using 0.03 mol/L NH4F and 0.025 mol/L HCl or 0.5 mol/L NaHCO3 (based on pH values), and analyzed using the molybdenum-blue method [22]. Total potassium (TK) was extracted using melt sodium hydroxide and then measured by flame emission spectrometry [22]. Available potassium (AK) was extracted using 1 mol/L NH4OAc and then measured by flame emission spectrometry [22]. Available S (AS) was extracted using 500 mg P/L Monocalcium phosphate [Ca(H2PO4)2] and then determined by the turbidimetric method [23]. Available Fe (AFe), Available Mn (AMn), Available Cu (ACu), and Available Zn (AZn) were extracted with diethylenetriaminepentaacetic acid (DTPA), and the extracted AFe, AMn, ACu, and AZn were determined with flame atomic absorption spectrometry [23]. Available B (AB) in soil was extracted with near-boiling deionized water and analyzed using a curcuma colorimetric procedure [23]. Quality control was based on the use of certified samples (GBW 07413) and analysis duplicates.


### 2.3. Availability Ratio

The availability ratio of a soil macronutrient (N, P, or K) at a given location is calculated using the following equation:(1)
R
=
A
(
u
)
T
(
u
)
×
1000
×
100
,
where R represents the availability ratio (%); u refers to a specific location; A
(
u
) is the available concentration of a soil macronutrient at location u with a unit of mg/kg; T
(
u
) is the total concentration of the soil macronutrient at location u with a unit of g/kg.


### 2.4. Correlation Analysis

Correlation analysis was conducted to reveal relationships between availability ratios of N, P, and K and selected soil properties (i.e., TN, AN, TP, AP, TK, AK, AFe, AMn, AS, ACu, AZn, AB, pH, SOM, and CEC). The significant level reported (i.e.,P
<
0.01) was based on Pearson’s correlation coefficients. The software used is SPSS R13.0 for Windows.We also performed stepwise regression analysis to select the main factors that affect the availability ratios of soil N, P, and K using the SPSS R13.0 software. Usually many possible explanatory variables may exist in a data set. Through a stepwise regression process, explanatory variables can be examined one by one, and the one that explains the most variation in the dependent variable will be added to the model at each step. When adding an extra variable will make no significant improvement in the amount of explained variation, the regression process will stop [24].


### 2.5. Geostatistical Analysis

Geostatistics provides methods to predict values at unsampled locations from values at sampled locations by taking into account the spatial correlation of sampled points. The variogram, the spatial measure for kriging, is an effective tool for evaluating spatial variability. The variogram may describe the spatial autocorrelation structure of a continuous variable and provide some insight into possible factors that affect data distribution. Spatial patterns of soil attributes following intrinsic stationarity can be described using variogram models [25].In this study, no apparent anisotropy was found for any studied variable through experimental variograms. So, all experimental variograms were in isotropic form and were fitted using basic math models, such as the spherical and exponential models, for kriging interpolation. Ordinary kriging was chosen to create the spatial distribution maps of soil N, P, and K contents, with the maximum search radius being set to the autocorrelation range of the corresponding variable. The interpolated grid had a resolution of 200.79 m × 200.79 m. Readers may see Chiles and Delfiner [26] for detailed description of kriging and the variogram. The Software GSLIB was used to conduct all the geostatistical work.


## 2.1. Study Area

Shayang County, Hubei Province in central China was selected as the study area. The county is among the most developed regions for cropping in Hubei Province and is located in the Jianghan Plain, bounded by the longitudes 112°03′ and 112°43′ east, and the latitudes 30°23′ and 30°57′ north, with an area of 2044 km2. It belongs to the northern subtropical monsoonal climate zone, with a temperate-humid climate throughout the year and four distinct seasons. The average annual temperature is 16.1°C and the mean annual precipitation is 1025.6 mm. In this region, paddy fields constitute the dominant land use type of arable lands.


## 2.2. Soil Sampling and Chemical Analysis

636 topsoil samples were collected from paddy fields in 2007. The sampling scheme considered the distributional uniformity of samples and also ensured that all samples were located in paddy fields (Figure1). All samples were taken in fall after harvest and before next cropping season to avoid the effect of fertilization during crop cultivation. For each soil sample, soils at 6–8 points within a small area of approximately 0.01 ha were collected from surface layers (0–15 cm) and then mixed and delivered to a laboratory for analysis. Exact sample locations were recorded using a hand-held global position system. All samples were air-dried at room temperature (20–22°C), crushed, and sieved into soil particles less than 2 mm. Portions of each sample (about 100 g) were further ground and sieved through 0.149 mm mesh. The prepared soil samples were stored for chemical analysis.Figure 1
The study area and sampling sites.Soil organic matter (SOM) was measured using the wet oxidation method at 180°C with a mixture of potassium dichromate and sulfuric acid [21]. Soil cation exchange capacity (CEC) was measured using 1.0 mol/L ammonium acetate solution [21]. Soil pH was measured using a pH meter (Sartorius Basic pH meter PB-10) with a soil/water ratio of 1 : 2.5 [21]. Total nitrogen (TN) was determined by the Kjeldahl method with H2SO4 + H2O2 digestion [22]. Available nitrogen (AN) was extracted using 2.0 mol/L KCl [22]. Total phosphorus (TP) was extracted using melt sodium hydroxide and measured using the molybdenum stibium anticolor method [22]. Available phosphorus (AP) was extracted using 0.03 mol/L NH4F and 0.025 mol/L HCl or 0.5 mol/L NaHCO3 (based on pH values), and analyzed using the molybdenum-blue method [22]. Total potassium (TK) was extracted using melt sodium hydroxide and then measured by flame emission spectrometry [22]. Available potassium (AK) was extracted using 1 mol/L NH4OAc and then measured by flame emission spectrometry [22]. Available S (AS) was extracted using 500 mg P/L Monocalcium phosphate [Ca(H2PO4)2] and then determined by the turbidimetric method [23]. Available Fe (AFe), Available Mn (AMn), Available Cu (ACu), and Available Zn (AZn) were extracted with diethylenetriaminepentaacetic acid (DTPA), and the extracted AFe, AMn, ACu, and AZn were determined with flame atomic absorption spectrometry [23]. Available B (AB) in soil was extracted with near-boiling deionized water and analyzed using a curcuma colorimetric procedure [23]. Quality control was based on the use of certified samples (GBW 07413) and analysis duplicates.


## 2.3. Availability Ratio

The availability ratio of a soil macronutrient (N, P, or K) at a given location is calculated using the following equation:(1)
R
=
A
(
u
)
T
(
u
)
×
1000
×
100
,
where R represents the availability ratio (%); u refers to a specific location; A
(
u
) is the available concentration of a soil macronutrient at location u with a unit of mg/kg; T
(
u
) is the total concentration of the soil macronutrient at location u with a unit of g/kg.


## 2.4. Correlation Analysis

Correlation analysis was conducted to reveal relationships between availability ratios of N, P, and K and selected soil properties (i.e., TN, AN, TP, AP, TK, AK, AFe, AMn, AS, ACu, AZn, AB, pH, SOM, and CEC). The significant level reported (i.e.,P
<
0.01) was based on Pearson’s correlation coefficients. The software used is SPSS R13.0 for Windows.We also performed stepwise regression analysis to select the main factors that affect the availability ratios of soil N, P, and K using the SPSS R13.0 software. Usually many possible explanatory variables may exist in a data set. Through a stepwise regression process, explanatory variables can be examined one by one, and the one that explains the most variation in the dependent variable will be added to the model at each step. When adding an extra variable will make no significant improvement in the amount of explained variation, the regression process will stop [24].


## 2.5. Geostatistical Analysis

Geostatistics provides methods to predict values at unsampled locations from values at sampled locations by taking into account the spatial correlation of sampled points. The variogram, the spatial measure for kriging, is an effective tool for evaluating spatial variability. The variogram may describe the spatial autocorrelation structure of a continuous variable and provide some insight into possible factors that affect data distribution. Spatial patterns of soil attributes following intrinsic stationarity can be described using variogram models [25].In this study, no apparent anisotropy was found for any studied variable through experimental variograms. So, all experimental variograms were in isotropic form and were fitted using basic math models, such as the spherical and exponential models, for kriging interpolation. Ordinary kriging was chosen to create the spatial distribution maps of soil N, P, and K contents, with the maximum search radius being set to the autocorrelation range of the corresponding variable. The interpolated grid had a resolution of 200.79 m × 200.79 m. Readers may see Chiles and Delfiner [26] for detailed description of kriging and the variogram. The Software GSLIB was used to conduct all the geostatistical work.


## 3. Results and Discussions

### 3.1. Correlation Analysis and Controlling Factors

The summary statistics of soil N, P, and K and selected soil properties (pH, SOM, CEC, AFe, AMn, AS, ACu, AZn, and AB) are listed in Table1. The coefficients of variation (CVs) for availability ratios of N, P, and K are 14.07%, 125.36%, and 61.8%, respectively. This means that the availability ratios of P and K have much larger variability than those of N in the study area. This situation may be caused by some extrinsic factors such as different fertilizer applications and seasonal cropping, which have stronger impacts on available P and K. In fact, TP, AP, and AK indeed have large CV values. Table 1 also shows that CEC, AB, ACu, AZn, AMn, AFe, and AS have relatively high CV values (>25%), reflecting the complex heterogeneity of paddy soils in the study area.Table 1
Summary statistics for macronutrient elements and other selected properties of paddy soils(
n
=
636
).
Range
Minimum
Maximum
Mean
S.D.
Skewness
Kurtosis
C.V. (%)
N availability ratio (%)
14.13
3.38
17.50
9.43
1.33
1.09
5.45
14.07
P availability ratio (%)
31.88
0.12
32.00
2.45
3.08
4.76
31.11
125.36
K availability ratio (%)
4.51
0.13
4.64
0.56
0.35
3.97
33.08
61.81
TN (g/kg)
1.87
0.48
2.35
1.44
0.31
-
0.37
0.14
21.69
AN (mg/kg)
181.00
45.00
226.00
134.64
28.15
-
0.45
0.18
20.91
TP (g/kg)
15.65
0.15
15.80
0.74
0.66
18.75
426.86
88.73
ln TP
4.66
-
1.90
2.76
-
0.41
0.46
-
0.29
4.07
-
112.04
AP (mg/kg)
93.00
2.00
95.00
14.03
11.21
2.55
9.73
79.89
ln AP
3.86
0.69
4.55
2.39
0.71
-
0.02
-
0.01
29.54
TK (g/kg)
28.36
7.33
35.69
19.96
3.41
0.91
1.87
17.07
ln TK
1.58
1.99
3.57
2.98
0.17
-
0.05
3.00
5.61
AK (mg/kg)
400.00
25.00
425.00
111.11
64.66
1.87
4.43
58.20
ln AK
2.83
3.22
6.05
4.57
0.53
0.01
0.34
11.60
SOM (g/kg)
35.65
4.89
40.54
23.11
5.31
-
0.60
0.07
22.97
pH
3.00
5.10
8.10
6.61
0.67
0.76
-
0.19
10.08
CEC (cmol (+)/kg)
30.42
1.08
31.50
12.78
4.19
0.04
0.83
32.79
AB (mg/kg)
5.04
0.10
5.14
0.64
0.55
3.59
18.16
86.73
ACu (mg/kg)
5.28
0.61
5.89
2.70
0.89
0.45
-
0.11
33.08
AZn (mg/kg)
8.39
0.12
8.51
1.38
0.73
2.40
15.20
52.58
AMn (mg/kg)
545.41
5.79
551.20
33.08
25.79
12.92
256.27
77.96
AFe (mg/kg)
349.74
6.03
355.77
118.78
89.18
0.67
-
0.99
75.08
AS (mg/kg)
182.65
7.85
190.50
41.25
21.99
1.42
4.58
53.30Pearson’s correlation coefficients were calculated for each variable to reveal the relationships between the availability ratios of macronutrients and the selected soil properties (see Table2). Apparently, the availability ratio of every macronutrient is negatively correlated with its total content and positively correlated with its available content, as shown in Table 2. However, the correlation situations between availability ratio and total and available contents are different for different macronutrients. The availability ratio of K has a very weak negative correlation with TK and a very strong positive correlation with AK. This implies that the TK content in soil is relatively stable over the study area, and a considerable amount of AK may come from outside input. The strong positive correlation between P availability ratio and AP also indicates that outside input may contribute a lot to soil AP. However, the availability ratio of N does not have such characteristic, and it is relatively more strongly impacted by TN. This means that most AN in paddy soils may come from SOM, and this point is supported by the strong positive correlation between AN and SOM. In general, the available contents of K and P are more decisive to their availability ratios; on the contrary, the total content of N has more impact on its availability ratio.Table 2
Pearson’s correlation coefficients between macronutrient elements and other selected soil properties.
TN
AN
TP
AP
TK
AK
pH
SOM
CEC
AB
AZn
ACu
AFe
AMn
AS
N availability ratio
-
0.374
*
*
0.186
*
*
-
0.036
-
0.038
-
0.032
-
0.117
*
*
0.022
-
0.135
*
*
-
0.025
-
0.014
-
0.088
*
-
0.102
*
*
-
0.049
-
0.058
-
0.099
*
P availability ratio
0.196
*
*
0.180
*
*
-
0.208
*
*
0.755
*
*
-
0.096
*
0.242
*
*
-
0.098
*
0.219
*
*
0.122
*
*
-
0.098
*
0.218
*
*
0.261
*
*
0.266
*
*
-
0.023
0.193
*
*
K availability ratio
0.255
*
*
0.182
*
*
0.065
0.332
*
*
-
0.039
0.901
*
*
-
0.047
0.236
*
*
0.252
*
*
-
0.010
0.367
*
*
0.226
*
*
0.123
*
*
0.053
0.064
TN
0.826
*
*
0.063
0.298
*
*
-
0.297
*
*
0.199
*
*
-
0.578
*
*
0.877
*
*
0.437
*
*
-
0.269
*
*
0.305
*
*
0.519
*
*
0.473
*
*
-
0.090
*
0.240
*
*
AN
0.042
0.272
*
*
-
0.345
*
*
0.118
*
*
-
0.629
*
*
0.853
*
*
0.459
*
*
-
0.320
*
*
0.246
*
*
0.484
*
*
0.469
*
*
-
0.142
*
*
0.200
*
*
TP
0.004
0.017
0.080
*
-
0.024
0.055
0.025
0.016
0.043
0.023
0.029
0.017
0.030
AP
-
0.122
*
*
0.306
*
*
-
0.180
*
*
0.295
*
*
0.205
*
*
-
0.109
*
*
0.200
*
*
0.310
*
*
0.244
*
*
-
0.044
0.188
*
*
TK
0.268
*
*
0.583
*
*
-
0.353
*
*
-
0.246
*
*
0.487
*
*
-
0.154
*
*
-
0.303
*
*
-
0.260
*
*
0.166
*
*
-
0.235
*
*
AK
0.089
*
0.169
*
*
0.197
*
*
0.097
*
0.327
*
*
0.161
*
*
0.066
0.098
*
0.002
*
P
<
0.05; *
*
P
<
0.01.Checking the correlation situations of other soil properties with N, P, and K, one can find that these selected properties each have similar positive or negative influences on AN and TN (Table2). This may explain why the availability ratio of N is less correlated with other soil properties (the largest correlation coefficient is with SOM, which is −0.135, P
<
0
.
01). However, the selected soil properties have much stronger correlations with AP and TK than with TP and AK, although these correlations are mostly weaker than those with TN and AN. For example, both pH and AB show strong correlations with TK but weak correlations with AK. The correlations of availability ratios of N, P, and K with the selected soil properties are generally not strong. The availability ratio of P is moderately correlated with ACu and AFe (r
s
=
0.261 and 0.266, respectively, P
<
0
.
01). The availability ratio of K is moderately correlated with CEC and AZn (r
s
=
0.252 and 0.367, respectively, P
<
0
.
01).Stepwise regression analysis was performed for the availability ratios of N, P, and K, by using all of other soil properties as independent variables. The results are given as follows:(2)
Y
(
N
availability
ratio
%
)
=
8.669
-
7.063
X
(
TN
)
+
0.075
X
(
AN
)
+
0.106
X
(
PH
)
+
0.004
X
(
AP
)
+
0.061
X
(
AZn
)
+
0.069
X
(
AB
)
,
000000000000000000000000000000000
(
R
2
=
0.91
)
,
Y
(
P
availability
ratio
%
)
=
-
2.533
+
0.203
X
(
AP
)
-
0.992
X
(
TP
)
+
0.004
X
(
AFe
)
+
0.359
X
(
PH
)
,
(
R
2
=
0.86
)
,
Y
(
K
availability
ratio
%
)
=
0.497
+
0.005
X
(
AK
)
-
0.036
X
(
TK
)
+
0.026
X
(
PH
)
+
0.027
X
(
AB
)
,
(
R
2
=
0.93
)
.
The coefficients of determination (R
2) provided after each of above regression models indicate that all of three multiple linear regression models are well fitted by sample data. It can be seen that the availability ratios of N, P, and K are all influenced by multiple factors which are also interrelated with each other. Stepwise regression could find major controlling factors. Above equations indicate that the availability ratio of N is mainly controlled by TN, AN, pH, AP, AZn, and AB, the availability ratio of P mainly depends on AP, TP, AFe, and pH, and that of K is mainly contributed by AK, TK, pH, and AB. These results are not fully consistent with the Pearson’s correlation coefficients in Table 2. For example, pH appears to be an important positive contributing factor to the availability ratios of all macronutrients here; however, its correlation coefficients with these availability ratios are quite small in absolute values and are not all positive. The reason may be that a Pearson’s correlation coefficient only reflects the linear correlation between two variables, but a regression coefficient in a multiple linear regression model is also related with other predictor variables and impacted by the multicollinearity problem.


### 3.2. Spatial Structure and Distributional Patterns

Geostatistical methods were used to analyze the spatial correlation structures of the total and available contents of N, P, and K in soil and spatially estimate their values at unsampled locations. Because kriging assumes the normal distribution for each estimated variable, it is necessary to check whether the total and available contents of N, P, and K in soil samples are approximately normally distributed or not. From Table1, one can see that the skewness and kurtosis indices of TN and AN are close to the standard value of 0, but those of TP, AP, TK, and AK are not. In this study, natural logarithmic transformation was used to reduce the skewness of the data distributions of P and K. The histograms of N, P, and K data used for geostatistical interpolation are shown in Figure 2. One can see that they tend to be normally distributed.Histograms of TN (a),ln
⁡
⁡
(
TP
) (b), ln
⁡
⁡
(
TK
) (c), AN (d), ln
⁡
⁡
(
AP
) (e), and ln
⁡
⁡
(
AK
) (f).
(a)
(b)
(c)
(d)
(e)
(f)Experimental variograms and fitted models for TN, AN,ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) are presented in Figure 3. Spherical models are chosen for TN and ln
⁡
⁡
(
AK
), Gaussian models for ln
⁡
⁡
(
TP
), and ln
⁡
⁡
(
AP
), and exponential models for AN and ln
⁡
⁡
(
TK
). The nugget/sill ratio (
C
0
/
(
C
+
C
0
)
) may be used as a criterion to reflect the spatial autodependency of a variable. Ratio values lower than 25%, between 25% and 75%, and higher than 75% correspond to strong, moderate, and weak spatial dependencies, respectively [11]. The nugget/sill ratios of variogram models for TN, AN, ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) all fall between 25% and 75%, which indicates that all of the six properties exhibit moderate spatial dependency. According to Cambardella et al. [11], these spatial autodependencies may be attributed to both intrinsic factors such as other soil properties and extrinsic factors such as human activities. The correlation range measures the spatial separate distance within which data are autocorrelated [27]. The approximate correlation ranges for TN, AN, ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) are 36.15 km, 59.46 km, 2.04 km, 5.00 km, 25.29 km, and 3.00 km, respectively. Apparently, AN, TN, and TK are autocorrelated in longer ranges than TP, AP, and AK are. This result is consistent with their CV values. This may imply that TP, AP, and AK contents are more sensitive to extrinsic factors such as fertilization.Experimental variograms of TN (a),ln
⁡
⁡
(
TP
) (b), ln
⁡
⁡
(
TK
) (c), AN (d), ln
⁡
⁡
(
AP
) (e), and ln
⁡
⁡
(
AK
) (f), and their fitted models.
(a)
(b)
(c)
(d)
(e)
(f)To map the spatial distributions of TN, AN, TP, AP, TK, and AK, we used ordinary kriging to interpolate their respective sample data (Figure4). After the spatial distribution maps of the total and available contents of a macronutrient were obtained, its availability ratio map could be further derived from the total content map and available content map by using (1) (i.e., at each spatial location, availability ratio = available content/total content). The spatial maps of TN and AN (Figures 4(a) and 4(b)) show similar patterns. This may be why TN and AN are strongly correlated (r
s
=
0.826, P
<
0
.
01). However, the spatial map of the availability ratio of N (Figure 5(a)) shows a different spatial pattern. This means that the contents of AN and TN are not spatially proportional and they are impacted by other factors with different degrees. AP and the availability ratio of P exhibit similar spatial trends (see Figures 4(d) and 5(b)). This should be caused by the strong correlation (r
s
=
0.755, P
<
0.01) between these two variables. Because AP and TP have almost no correlation (see Table 2), consequently the spatial pattern of P availability ratio is quite different from that of TP. The spatial distribution maps of AK and K availability ratio also show similar spatial patterns, which should be caused by their strong (nearly linear) correlation (r
s
=
0.901, P
<
0
.
01) (see Figures 4(f) and 5(c)). The correlation between AK and TK is not strong, which may be why they have different spatial patterns and the spatial pattern of K availability ratio is very different from that of TK.Spatial distribution maps of TN (a), AN (b), TP (c), AP (d), TK (e), and AK (f), interpolated by ordinary kriging.
(a)
(b)
(c)
(d)
(e)
(f)Spatial distribution maps of availability ratios of N (a), P (b), and K (c). Note that (a) is derived from Figures4(a) and 4(b), (b) is derived from Figures 4(c) and 4(d), and (c) is derived from Figures 4(e) and 4(f).
(a)
(b)
(c)The short autocorrelation ranges of TP, AP, and AK are reflected on their spatial distribution maps with fragmentary patterns and this characteristic is further propagated to the spatial distribution maps of the availability ratios of P and K. Consequently, one can see in Figure5 that the spatial distributional characteristics of the availability ratios of N, P, and K are very different: (i) the availability ratio of N changes smoothly with large patches while that of K and P vary abruptly as small pieces or spots in the study area; (ii) high availability ratio values of N mainly appear in several subareas in the east and north of the study area, and high availability ratio values of P mainly occur as a series of spots along a strip in the southwest, while high availability ratio values of K mainly emerge as isolated spots, often bordered with very low values. These availability ratio maps clearly show where a macronutrient has low availability and where its availability is high. Thus, concrete measures might be taken to increase the low availability or decrease the high availability of a macronutrient in a specific area to meet the requirements of agricultural production and environmental protection. For example, changing soil pH may help change the efficiency of soil nitrogen to be absorbed by plants.It should be noted that the data value ranges (from minimum value to maximum value in each legend) in the legends of Figures4 and 5 are not consistent with the corresponding value ranges of sample data provided in Table 2. The major reason may be that the locations of interpolated grid values by kriging for mapping do not overlap those of sample data, thus the legends show only the value ranges of interpolated data and sample data values are ignored in the maps.


## 3.1. Correlation Analysis and Controlling Factors

The summary statistics of soil N, P, and K and selected soil properties (pH, SOM, CEC, AFe, AMn, AS, ACu, AZn, and AB) are listed in Table1. The coefficients of variation (CVs) for availability ratios of N, P, and K are 14.07%, 125.36%, and 61.8%, respectively. This means that the availability ratios of P and K have much larger variability than those of N in the study area. This situation may be caused by some extrinsic factors such as different fertilizer applications and seasonal cropping, which have stronger impacts on available P and K. In fact, TP, AP, and AK indeed have large CV values. Table 1 also shows that CEC, AB, ACu, AZn, AMn, AFe, and AS have relatively high CV values (>25%), reflecting the complex heterogeneity of paddy soils in the study area.Table 1
Summary statistics for macronutrient elements and other selected properties of paddy soils(
n
=
636
).
Range
Minimum
Maximum
Mean
S.D.
Skewness
Kurtosis
C.V. (%)
N availability ratio (%)
14.13
3.38
17.50
9.43
1.33
1.09
5.45
14.07
P availability ratio (%)
31.88
0.12
32.00
2.45
3.08
4.76
31.11
125.36
K availability ratio (%)
4.51
0.13
4.64
0.56
0.35
3.97
33.08
61.81
TN (g/kg)
1.87
0.48
2.35
1.44
0.31
-
0.37
0.14
21.69
AN (mg/kg)
181.00
45.00
226.00
134.64
28.15
-
0.45
0.18
20.91
TP (g/kg)
15.65
0.15
15.80
0.74
0.66
18.75
426.86
88.73
ln TP
4.66
-
1.90
2.76
-
0.41
0.46
-
0.29
4.07
-
112.04
AP (mg/kg)
93.00
2.00
95.00
14.03
11.21
2.55
9.73
79.89
ln AP
3.86
0.69
4.55
2.39
0.71
-
0.02
-
0.01
29.54
TK (g/kg)
28.36
7.33
35.69
19.96
3.41
0.91
1.87
17.07
ln TK
1.58
1.99
3.57
2.98
0.17
-
0.05
3.00
5.61
AK (mg/kg)
400.00
25.00
425.00
111.11
64.66
1.87
4.43
58.20
ln AK
2.83
3.22
6.05
4.57
0.53
0.01
0.34
11.60
SOM (g/kg)
35.65
4.89
40.54
23.11
5.31
-
0.60
0.07
22.97
pH
3.00
5.10
8.10
6.61
0.67
0.76
-
0.19
10.08
CEC (cmol (+)/kg)
30.42
1.08
31.50
12.78
4.19
0.04
0.83
32.79
AB (mg/kg)
5.04
0.10
5.14
0.64
0.55
3.59
18.16
86.73
ACu (mg/kg)
5.28
0.61
5.89
2.70
0.89
0.45
-
0.11
33.08
AZn (mg/kg)
8.39
0.12
8.51
1.38
0.73
2.40
15.20
52.58
AMn (mg/kg)
545.41
5.79
551.20
33.08
25.79
12.92
256.27
77.96
AFe (mg/kg)
349.74
6.03
355.77
118.78
89.18
0.67
-
0.99
75.08
AS (mg/kg)
182.65
7.85
190.50
41.25
21.99
1.42
4.58
53.30Pearson’s correlation coefficients were calculated for each variable to reveal the relationships between the availability ratios of macronutrients and the selected soil properties (see Table2). Apparently, the availability ratio of every macronutrient is negatively correlated with its total content and positively correlated with its available content, as shown in Table 2. However, the correlation situations between availability ratio and total and available contents are different for different macronutrients. The availability ratio of K has a very weak negative correlation with TK and a very strong positive correlation with AK. This implies that the TK content in soil is relatively stable over the study area, and a considerable amount of AK may come from outside input. The strong positive correlation between P availability ratio and AP also indicates that outside input may contribute a lot to soil AP. However, the availability ratio of N does not have such characteristic, and it is relatively more strongly impacted by TN. This means that most AN in paddy soils may come from SOM, and this point is supported by the strong positive correlation between AN and SOM. In general, the available contents of K and P are more decisive to their availability ratios; on the contrary, the total content of N has more impact on its availability ratio.Table 2
Pearson’s correlation coefficients between macronutrient elements and other selected soil properties.
TN
AN
TP
AP
TK
AK
pH
SOM
CEC
AB
AZn
ACu
AFe
AMn
AS
N availability ratio
-
0.374
*
*
0.186
*
*
-
0.036
-
0.038
-
0.032
-
0.117
*
*
0.022
-
0.135
*
*
-
0.025
-
0.014
-
0.088
*
-
0.102
*
*
-
0.049
-
0.058
-
0.099
*
P availability ratio
0.196
*
*
0.180
*
*
-
0.208
*
*
0.755
*
*
-
0.096
*
0.242
*
*
-
0.098
*
0.219
*
*
0.122
*
*
-
0.098
*
0.218
*
*
0.261
*
*
0.266
*
*
-
0.023
0.193
*
*
K availability ratio
0.255
*
*
0.182
*
*
0.065
0.332
*
*
-
0.039
0.901
*
*
-
0.047
0.236
*
*
0.252
*
*
-
0.010
0.367
*
*
0.226
*
*
0.123
*
*
0.053
0.064
TN
0.826
*
*
0.063
0.298
*
*
-
0.297
*
*
0.199
*
*
-
0.578
*
*
0.877
*
*
0.437
*
*
-
0.269
*
*
0.305
*
*
0.519
*
*
0.473
*
*
-
0.090
*
0.240
*
*
AN
0.042
0.272
*
*
-
0.345
*
*
0.118
*
*
-
0.629
*
*
0.853
*
*
0.459
*
*
-
0.320
*
*
0.246
*
*
0.484
*
*
0.469
*
*
-
0.142
*
*
0.200
*
*
TP
0.004
0.017
0.080
*
-
0.024
0.055
0.025
0.016
0.043
0.023
0.029
0.017
0.030
AP
-
0.122
*
*
0.306
*
*
-
0.180
*
*
0.295
*
*
0.205
*
*
-
0.109
*
*
0.200
*
*
0.310
*
*
0.244
*
*
-
0.044
0.188
*
*
TK
0.268
*
*
0.583
*
*
-
0.353
*
*
-
0.246
*
*
0.487
*
*
-
0.154
*
*
-
0.303
*
*
-
0.260
*
*
0.166
*
*
-
0.235
*
*
AK
0.089
*
0.169
*
*
0.197
*
*
0.097
*
0.327
*
*
0.161
*
*
0.066
0.098
*
0.002
*
P
<
0.05; *
*
P
<
0.01.Checking the correlation situations of other soil properties with N, P, and K, one can find that these selected properties each have similar positive or negative influences on AN and TN (Table2). This may explain why the availability ratio of N is less correlated with other soil properties (the largest correlation coefficient is with SOM, which is −0.135, P
<
0
.
01). However, the selected soil properties have much stronger correlations with AP and TK than with TP and AK, although these correlations are mostly weaker than those with TN and AN. For example, both pH and AB show strong correlations with TK but weak correlations with AK. The correlations of availability ratios of N, P, and K with the selected soil properties are generally not strong. The availability ratio of P is moderately correlated with ACu and AFe (r
s
=
0.261 and 0.266, respectively, P
<
0
.
01). The availability ratio of K is moderately correlated with CEC and AZn (r
s
=
0.252 and 0.367, respectively, P
<
0
.
01).Stepwise regression analysis was performed for the availability ratios of N, P, and K, by using all of other soil properties as independent variables. The results are given as follows:(2)
Y
(
N
availability
ratio
%
)
=
8.669
-
7.063
X
(
TN
)
+
0.075
X
(
AN
)
+
0.106
X
(
PH
)
+
0.004
X
(
AP
)
+
0.061
X
(
AZn
)
+
0.069
X
(
AB
)
,
000000000000000000000000000000000
(
R
2
=
0.91
)
,
Y
(
P
availability
ratio
%
)
=
-
2.533
+
0.203
X
(
AP
)
-
0.992
X
(
TP
)
+
0.004
X
(
AFe
)
+
0.359
X
(
PH
)
,
(
R
2
=
0.86
)
,
Y
(
K
availability
ratio
%
)
=
0.497
+
0.005
X
(
AK
)
-
0.036
X
(
TK
)
+
0.026
X
(
PH
)
+
0.027
X
(
AB
)
,
(
R
2
=
0.93
)
.
The coefficients of determination (R
2) provided after each of above regression models indicate that all of three multiple linear regression models are well fitted by sample data. It can be seen that the availability ratios of N, P, and K are all influenced by multiple factors which are also interrelated with each other. Stepwise regression could find major controlling factors. Above equations indicate that the availability ratio of N is mainly controlled by TN, AN, pH, AP, AZn, and AB, the availability ratio of P mainly depends on AP, TP, AFe, and pH, and that of K is mainly contributed by AK, TK, pH, and AB. These results are not fully consistent with the Pearson’s correlation coefficients in Table 2. For example, pH appears to be an important positive contributing factor to the availability ratios of all macronutrients here; however, its correlation coefficients with these availability ratios are quite small in absolute values and are not all positive. The reason may be that a Pearson’s correlation coefficient only reflects the linear correlation between two variables, but a regression coefficient in a multiple linear regression model is also related with other predictor variables and impacted by the multicollinearity problem.


## 3.2. Spatial Structure and Distributional Patterns

Geostatistical methods were used to analyze the spatial correlation structures of the total and available contents of N, P, and K in soil and spatially estimate their values at unsampled locations. Because kriging assumes the normal distribution for each estimated variable, it is necessary to check whether the total and available contents of N, P, and K in soil samples are approximately normally distributed or not. From Table1, one can see that the skewness and kurtosis indices of TN and AN are close to the standard value of 0, but those of TP, AP, TK, and AK are not. In this study, natural logarithmic transformation was used to reduce the skewness of the data distributions of P and K. The histograms of N, P, and K data used for geostatistical interpolation are shown in Figure 2. One can see that they tend to be normally distributed.Histograms of TN (a),ln
⁡
⁡
(
TP
) (b), ln
⁡
⁡
(
TK
) (c), AN (d), ln
⁡
⁡
(
AP
) (e), and ln
⁡
⁡
(
AK
) (f).
(a)
(b)
(c)
(d)
(e)
(f)Experimental variograms and fitted models for TN, AN,ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) are presented in Figure 3. Spherical models are chosen for TN and ln
⁡
⁡
(
AK
), Gaussian models for ln
⁡
⁡
(
TP
), and ln
⁡
⁡
(
AP
), and exponential models for AN and ln
⁡
⁡
(
TK
). The nugget/sill ratio (
C
0
/
(
C
+
C
0
)
) may be used as a criterion to reflect the spatial autodependency of a variable. Ratio values lower than 25%, between 25% and 75%, and higher than 75% correspond to strong, moderate, and weak spatial dependencies, respectively [11]. The nugget/sill ratios of variogram models for TN, AN, ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) all fall between 25% and 75%, which indicates that all of the six properties exhibit moderate spatial dependency. According to Cambardella et al. [11], these spatial autodependencies may be attributed to both intrinsic factors such as other soil properties and extrinsic factors such as human activities. The correlation range measures the spatial separate distance within which data are autocorrelated [27]. The approximate correlation ranges for TN, AN, ln
⁡
⁡
(
TP
), ln
⁡
⁡
(
AP
), ln
⁡
⁡
(
TK
), and ln
⁡
⁡
(
AK
) are 36.15 km, 59.46 km, 2.04 km, 5.00 km, 25.29 km, and 3.00 km, respectively. Apparently, AN, TN, and TK are autocorrelated in longer ranges than TP, AP, and AK are. This result is consistent with their CV values. This may imply that TP, AP, and AK contents are more sensitive to extrinsic factors such as fertilization.Experimental variograms of TN (a),ln
⁡
⁡
(
TP
) (b), ln
⁡
⁡
(
TK
) (c), AN (d), ln
⁡
⁡
(
AP
) (e), and ln
⁡
⁡
(
AK
) (f), and their fitted models.
(a)
(b)
(c)
(d)
(e)
(f)To map the spatial distributions of TN, AN, TP, AP, TK, and AK, we used ordinary kriging to interpolate their respective sample data (Figure4). After the spatial distribution maps of the total and available contents of a macronutrient were obtained, its availability ratio map could be further derived from the total content map and available content map by using (1) (i.e., at each spatial location, availability ratio = available content/total content). The spatial maps of TN and AN (Figures 4(a) and 4(b)) show similar patterns. This may be why TN and AN are strongly correlated (r
s
=
0.826, P
<
0
.
01). However, the spatial map of the availability ratio of N (Figure 5(a)) shows a different spatial pattern. This means that the contents of AN and TN are not spatially proportional and they are impacted by other factors with different degrees. AP and the availability ratio of P exhibit similar spatial trends (see Figures 4(d) and 5(b)). This should be caused by the strong correlation (r
s
=
0.755, P
<
0.01) between these two variables. Because AP and TP have almost no correlation (see Table 2), consequently the spatial pattern of P availability ratio is quite different from that of TP. The spatial distribution maps of AK and K availability ratio also show similar spatial patterns, which should be caused by their strong (nearly linear) correlation (r
s
=
0.901, P
<
0
.
01) (see Figures 4(f) and 5(c)). The correlation between AK and TK is not strong, which may be why they have different spatial patterns and the spatial pattern of K availability ratio is very different from that of TK.Spatial distribution maps of TN (a), AN (b), TP (c), AP (d), TK (e), and AK (f), interpolated by ordinary kriging.
(a)
(b)
(c)
(d)
(e)
(f)Spatial distribution maps of availability ratios of N (a), P (b), and K (c). Note that (a) is derived from Figures4(a) and 4(b), (b) is derived from Figures 4(c) and 4(d), and (c) is derived from Figures 4(e) and 4(f).
(a)
(b)
(c)The short autocorrelation ranges of TP, AP, and AK are reflected on their spatial distribution maps with fragmentary patterns and this characteristic is further propagated to the spatial distribution maps of the availability ratios of P and K. Consequently, one can see in Figure5 that the spatial distributional characteristics of the availability ratios of N, P, and K are very different: (i) the availability ratio of N changes smoothly with large patches while that of K and P vary abruptly as small pieces or spots in the study area; (ii) high availability ratio values of N mainly appear in several subareas in the east and north of the study area, and high availability ratio values of P mainly occur as a series of spots along a strip in the southwest, while high availability ratio values of K mainly emerge as isolated spots, often bordered with very low values. These availability ratio maps clearly show where a macronutrient has low availability and where its availability is high. Thus, concrete measures might be taken to increase the low availability or decrease the high availability of a macronutrient in a specific area to meet the requirements of agricultural production and environmental protection. For example, changing soil pH may help change the efficiency of soil nitrogen to be absorbed by plants.It should be noted that the data value ranges (from minimum value to maximum value in each legend) in the legends of Figures4 and 5 are not consistent with the corresponding value ranges of sample data provided in Table 2. The major reason may be that the locations of interpolated grid values by kriging for mapping do not overlap those of sample data, thus the legends show only the value ranges of interpolated data and sample data values are ignored in the maps.


## 4. Conclusion

The spatial variability of soil N, P, and K and their availability ratios in paddy soils in a county were explored. After a long history of various land management, the spatial variability of N, P, and K and their availability ratios in soils are not only dependent on soil parent materials but also impacted by anthropogenic activities. Statistical analysis shows that the availability ratios of P and K have larger variation than that of N in the study area, probably due to the heterogeneous land use and management. The contents of TN, AN, and TK in paddy fields have relatively longer correlation ranges than those of TP, AP, and AK. The availability ratios of N, P, and K are all affected by multiple factors that are also interrelated. The main factors controlling the availability ratios of N, P, and K were determined using stepwise regression analysis. N availability ratio is mainly related to TN, AN, pH, AP, AZn, and AB; P availability ratio is mainly controlled by AP, TP, AFe, and pH; and K availability ratio is mainly affected by AK, TK, pH, and AB.TN and AN contents have similar spatial patterns, but N availability ratio has a different spatial pattern; however, the spatial patterns of availability ratios of P and K are only similar to those of their available components (i.e., AP and AK), respectively. Availability ratios of nutrients in soil may reflect soil’s nutrient holding (or fixing) ability. The different spatial variability characteristics among the total content, the available content, and the availability ratio of each macronutrient indicate that these three indices are all useful to reflect the general situation of the macronutrient in soil and their spatial maps are all useful data for agricultural and environmental management.Because the availability ratio of a macronutrient is related with its availability (or efficiency) to plants and its risk of loss into deep soil or water body, proper adjustment of macronutrient availability ratios in soil through some measures (e.g., fertilization, liming to change soil pH) may be helpful to improving agricultural productivity and reducing environmental hazards. Therefore, spatial distribution maps of macronutrient availability ratios may provide useful information to guide site-specific soil management.



---

*Source: 689482-2014-03-11.xml*