2. Materials and methods

#### 2.1. Experimental soils

was presented in 1856 [1]. In 1907, Edgar Buckingham applied "capillary potential" to soil water for the first time, which showed the energy state of soil water. Green and Ampt proposed an infiltration model based on capillary theory [2]. Richards introduced Darcy's law to describe soil water flow [3]. Philip presented a basic equation to describe the water movement in a one-dimensional vertical soil column [4]. In addition, the formula for infiltration of the Kostiakov infiltration model, Horton infiltration model, and Holtan infiltration model is also used [5–7]. Shu puts forward the model of capillary bundle infiltration [8]. Ghosh combines the one-dimensional vertical infiltration formula and Kostiakov empirical formula to obtain the new infiltration formula [9]. Many experts and scholars have proposed some new methods because it is difficult to find out the results in accordance with the actual results. Parlange presented an approximate solution for Richards equation [10], and Hogarth and Parlange improved the solution [11]. The finite difference method and finite element method were also used to solve the solution of the water eqaution [12]. Yang and Lei established a numerical model for the one-dimensional saturated water flow in the FORTRAN language and verified it in laboratory [13]. There are a large number of basic formulas and empirical formulas to describe the process of one-dimensional soil water movement. However, both the classical infiltration model and empirical model have different parameters which make them difficult and time-consuming. In this study, we want to find a simple and feasible method to determine

Soil thermal conductivity is not only one of the important indexes of soil thermal properties but also an important parameter for simulating the soil hydro-thermal-solute-coupled model. How to estimate soil thermal conductivity quickly and accurately is one of main contents of studying soil thermal properties [14]. At present, a number of indirect estimation models to describe the relationship between thermal conductivity and soil texture, bulk density, water content, and organic matter were proposed by domestic and foreign scholars [15–23]. There are two types of indirect estimation models: empirical models [15, 16] and semi-theoretical models [17–20]. The empirical model mainly established the relationship between thermal conductivity and soil moisture content, such as the Chung-Horton model and Campbell model. These models are simple to calculate, but the model parameters are uncertain which will lead to large errors between the calculated data and measured values due to the difference of soil qualities in different regions [21]. The semi-theoretical model showed the relationship between thermal conductivity and soil saturation, such as the Johansen model, Côté-Konrad model, and Lu-Ren model [22]. These models have a theoretical basis and have given the model parameters for different soil textures. However, the model parameter values varied greatly with different soil particles and organic matter content, which limited the application of this model. In general, different models have their own advantages and disadvantages, but the effect of particle composition and organic matter content on the parameters of different types of soils needs further study. In this chapter, the thermal conductivity of undisturbed soil was measured by heat pulse methods. So, in this study, a new method based on analyzing the influence of soil particle composition on thermal conductivity, the relationship between thermal conductivity, saturation, bulk density, soil particle composition, and organic matter, was established. The improved Côté-Konrad and Lu-Ren models were also proposed to provide a reference method for obtaining soil thermal conductivity in a simple and rapid manner.

soil hydraulic parameters.

128 Hydrology of Artificial and Controlled Experiments

In this study, four soils were collected, and the soil's physical characteristics were listed in Table 1.

#### 2.2. Experimental measurement

#### 2.2.1. Horizonal infiltration experiment

Four soils in Table 1 were collected for the infiltration experiments, and the negative hydraulic heads were designed as 2.5, 6, 9, 12, 15, and 18 cm. Soil samples were filled in the column, every 10 cm, and the bulk density of soil was designed as 1.4 g/cm<sup>3</sup> . The length of h2 and h3 were measured, and the values of h1 were calculated by the formula p = h3-h1 + h2. The


Table 1. The soil's physical characteristics.

water head in pressure regulator pipe values was adjusted by the values of h1 (Figure 1). The standpipe was filled with distilled water before the experiment. Opening the right valve, the experiment continued until the bubble was emptied. For each infiltration measurement, cumulative infiltration was recorded every minute until it reached a steady state.

temperature in the process with time after the middle probe sent the heat pulse (Figure 3). The diameter, length, and space distance of the three probes were 1.3, 40, and 6 mm, respectively (as shown in Figure 1). The 5–6 gL<sup>1</sup> agar solution was used to demarcate in advance in actual, which was to prevent natural convection of water when heated. The Data Collector (US CR1000 Data Collector) controls the heated input via a relay, and the electric current was determined by a precise resistance (10 Ω) of the assigned voltage. The data collector also recorded the temperature change of the sensing probe at intervals of 1 s. The volumetric heat capacity of the agar solution is 4.18MJm<sup>3</sup>C1. The distance r was obtained by a non-linear fitting temperature–time curve and averaged by repeating the calibration process 10 times.

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

131

To study the variation characteristics of the undisturbed soil thermal conductivity, the ring knife was used to take samples in the experimental ground. Each measurement point was arranged as follows: 10 measuring points per column, step length of 3 m between 2 points, and setting in 2 columns; 4 kinds of water contents were given to measure soil thermal conductivity in each measurement point, and the actual water contents were determined by the measured value at the end of the measurement (that means that the actual moisture content of the soil sample is supposed to be equal to that in the ring knife after finishing the measurement).

The basic component of the experiments was a rain simulator, which could generate a variable intensity of rainfall. The nozzles used to simulate rainfall were 15 m from the soil surface. We used six steel soil flumes with the following dimensions: 1 m in length 0.40 m in width 0.50 m in height. The flumes were filled with soil to a depth of 0.35 m; this depth allowed

2.2.3. Nutrient runoff experiment

Figure 3. Schematic diagram of experimental apparatus.

#### 2.2.2. Soil thermal experiment

The test equipment uses a three-probe heat pulse probe (Figure 2) which was connected to the data collector, and sensor probes were used on two sides to observe and monitor the changing

Figure 1. The sketch map of experimental equipment for horizonal soil column.

Figure 2. Schematic diagram of the heat pulse probe.

temperature in the process with time after the middle probe sent the heat pulse (Figure 3). The diameter, length, and space distance of the three probes were 1.3, 40, and 6 mm, respectively (as shown in Figure 1). The 5–6 gL<sup>1</sup> agar solution was used to demarcate in advance in actual, which was to prevent natural convection of water when heated. The Data Collector (US CR1000 Data Collector) controls the heated input via a relay, and the electric current was determined by a precise resistance (10 Ω) of the assigned voltage. The data collector also recorded the temperature change of the sensing probe at intervals of 1 s. The volumetric heat capacity of the agar solution is 4.18MJm<sup>3</sup>C1. The distance r was obtained by a non-linear fitting temperature–time curve and averaged by repeating the calibration process 10 times.

To study the variation characteristics of the undisturbed soil thermal conductivity, the ring knife was used to take samples in the experimental ground. Each measurement point was arranged as follows: 10 measuring points per column, step length of 3 m between 2 points, and setting in 2 columns; 4 kinds of water contents were given to measure soil thermal conductivity in each measurement point, and the actual water contents were determined by the measured value at the end of the measurement (that means that the actual moisture content of the soil sample is supposed to be equal to that in the ring knife after finishing the measurement).

#### 2.2.3. Nutrient runoff experiment

water head in pressure regulator pipe values was adjusted by the values of h1 (Figure 1). The standpipe was filled with distilled water before the experiment. Opening the right valve, the experiment continued until the bubble was emptied. For each infiltration measurement, cumu-

The test equipment uses a three-probe heat pulse probe (Figure 2) which was connected to the data collector, and sensor probes were used on two sides to observe and monitor the changing

lative infiltration was recorded every minute until it reached a steady state.

Figure 1. The sketch map of experimental equipment for horizonal soil column.

Figure 2. Schematic diagram of the heat pulse probe.

2.2.2. Soil thermal experiment

130 Hydrology of Artificial and Controlled Experiments

The basic component of the experiments was a rain simulator, which could generate a variable intensity of rainfall. The nozzles used to simulate rainfall were 15 m from the soil surface. We used six steel soil flumes with the following dimensions: 1 m in length 0.40 m in width 0.50 m in height. The flumes were filled with soil to a depth of 0.35 m; this depth allowed

Figure 3. Schematic diagram of experimental apparatus.

infiltration without causing the bottoms of the flumes to become dank and left a 0.15 cm "lip" above the soil level to prevent water losses from splashing. The flumes' angle of inclination could be varied between 0 and 30 (Figure 4). The experiments were performed from April to September 2010 in a laboratory for simulating artificial rainfall at the Institute for Soil and Water Conservation, Chinese Academy of Sciences, Shaanxi Province, China.

2.3. Theory

as:

2.3.1. Models of infiltration process

The Kostiakov model [5] was presented by large amount of experiments and can be expressed

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

where I is cumulative infiltration (cm), t is the infiltration time (min), a and b are empirical constants. The empirical constants a and b have no physical meaning and are obtained by

Philip [4] proposed an infiltration equation based on Boltzmann transformation as power

where S is the soil sorptivity (cm/min1/2) and A is stable infiltration rate (cm/min). The result is

A horizontal one-dimensional infiltration equation, neglecting gravity action, can be expressed

A vertical one-dimensional infiltration equation, neglecting gravity action, can be expressed as:

Wang et al. [28] use k hð Þ¼ <sup>k</sup>sð Þ <sup>h</sup>d=<sup>h</sup> <sup>m</sup> and ð Þ <sup>θ</sup><sup>0</sup> � <sup>θ</sup><sup>r</sup> <sup>=</sup>ð Þ¼ <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>r</sup> ð Þ <sup>h</sup>d=<sup>h</sup> <sup>n</sup> to calculate unsaturated hydraulic conductivity and soil moisture characteristic curves. The relationship between

m þ n � 1

t 1

<sup>I</sup> <sup>¼</sup> xfð Þ <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>i</sup> <sup>1</sup> � <sup>n</sup>

<sup>q</sup> <sup>¼</sup> <sup>1</sup> xf hd

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2hdksð Þ m þ n � 1 ð Þ m<sup>2</sup> � 1 ð Þ θ<sup>s</sup> � θ<sup>i</sup>

wetting front and cumulative infiltration time can be expressed as:

xf ¼

s

<sup>I</sup> <sup>¼</sup> atb (1)

http://dx.doi.org/10.5772/intechopen.76280

133

<sup>I</sup> <sup>¼</sup> St<sup>0</sup>:<sup>5</sup> <sup>þ</sup> At <sup>þ</sup> Bt<sup>1</sup>:<sup>5</sup> <sup>þ</sup> <sup>⋯</sup> (2)

<sup>I</sup> <sup>¼</sup> St<sup>0</sup>:<sup>5</sup> (3)

<sup>I</sup> <sup>¼</sup> St<sup>0</sup>:<sup>5</sup> <sup>þ</sup> At (4)

� � (5)

<sup>m</sup> � <sup>1</sup> (6)

<sup>2</sup> (7)

2.3.1.1. Kostiakov model

experimental data fitting.

exactly enough with two terms.

2.3.1.2. Philip model

2.3.1.3. Wang's model

series:

as:

Three treatments were established to test our model. In treatment 1, three initial levels of soil moisture (5, 10, and 20%, measured gravimetrically) were used to study the influence of the soil's initial water content on our model. The rainfall rate was 90 mm/h and the slope gradient was 5. Treatment 2 was designed to investigate the influence of variation in the rainfall intensity on our model. Three different rainfall intensities (60, 96, and 129 mm/h) were examined, with an initial soil moisture content of 10% and a slope gradient of 5. Treatment 3 was designed to assess the influence of the slope gradient. Slopes of 5, 15, and 25 were investigated, with a rainfall intensity of 90 mm/h and an initial (gravimetric) soil moisture content of 10%. All treatments were run three times.

The soil samples were sieved (0.004 m in aperture) to remove coarse rock and debris and then air dried (to about 2%, gravimetrically). Potassium, used as a tracer, was dissolved in water and added to the test soils based on their designed soil water contents and potassium concentrations; the soil was then thoroughly mixed. The soil flume was filled with the prepared soil sample layer by layer to achieve a dry bulk density of 1.35 g/cm<sup>3</sup> . To obtain a flat surface, a sharp-edged straight blade was used to remove excess soil. The soil surface was covered with plastic for approximately 24 h before the beginning of the experiments. During the experiments, the outflow from one of the holes in the flume was collected into plastic containers every minute to measure the amount of runoff and its sediment and soil concentrations. We directly measured the depth of the exchange layers along a vertical section. The potassium content in the runoff was measured with an atomic absorption spectrophotometer (Perkin-Elmer 5100ZL). The soil water content was measured by drying, and the sediment was isolated by filtration on filter paper and weighed after drying.

Figure 4. Experimental setup of artificial rainfall.

#### 2.3. Theory

infiltration without causing the bottoms of the flumes to become dank and left a 0.15 cm "lip" above the soil level to prevent water losses from splashing. The flumes' angle of inclination could be varied between 0 and 30 (Figure 4). The experiments were performed from April to September 2010 in a laboratory for simulating artificial rainfall at the Institute for Soil and

Three treatments were established to test our model. In treatment 1, three initial levels of soil moisture (5, 10, and 20%, measured gravimetrically) were used to study the influence of the soil's initial water content on our model. The rainfall rate was 90 mm/h and the slope gradient was 5. Treatment 2 was designed to investigate the influence of variation in the rainfall intensity on our model. Three different rainfall intensities (60, 96, and 129 mm/h) were examined, with an initial soil moisture content of 10% and a slope gradient of 5. Treatment 3 was designed to assess the influence of the slope gradient. Slopes of 5, 15, and 25 were investigated, with a rainfall intensity of 90 mm/h and an initial (gravimetric) soil moisture content of

The soil samples were sieved (0.004 m in aperture) to remove coarse rock and debris and then air dried (to about 2%, gravimetrically). Potassium, used as a tracer, was dissolved in water and added to the test soils based on their designed soil water contents and potassium concentrations; the soil was then thoroughly mixed. The soil flume was filled with the prepared soil

sharp-edged straight blade was used to remove excess soil. The soil surface was covered with plastic for approximately 24 h before the beginning of the experiments. During the experiments, the outflow from one of the holes in the flume was collected into plastic containers every minute to measure the amount of runoff and its sediment and soil concentrations. We directly measured the depth of the exchange layers along a vertical section. The potassium content in the runoff was measured with an atomic absorption spectrophotometer (Perkin-Elmer 5100ZL). The soil water content was measured by drying, and the sediment was isolated

. To obtain a flat surface, a

Water Conservation, Chinese Academy of Sciences, Shaanxi Province, China.

sample layer by layer to achieve a dry bulk density of 1.35 g/cm<sup>3</sup>

by filtration on filter paper and weighed after drying.

Figure 4. Experimental setup of artificial rainfall.

10%. All treatments were run three times.

132 Hydrology of Artificial and Controlled Experiments

#### 2.3.1. Models of infiltration process

#### 2.3.1.1. Kostiakov model

The Kostiakov model [5] was presented by large amount of experiments and can be expressed as:

$$I = at^b\tag{1}$$

where I is cumulative infiltration (cm), t is the infiltration time (min), a and b are empirical constants. The empirical constants a and b have no physical meaning and are obtained by experimental data fitting.

#### 2.3.1.2. Philip model

Philip [4] proposed an infiltration equation based on Boltzmann transformation as power series:

$$I = St^{0.5} + At + Bt^{1.5} + \cdots \tag{2}$$

where S is the soil sorptivity (cm/min1/2) and A is stable infiltration rate (cm/min). The result is exactly enough with two terms.

A horizontal one-dimensional infiltration equation, neglecting gravity action, can be expressed as:

$$I = St^{0.5} \tag{3}$$

A vertical one-dimensional infiltration equation, neglecting gravity action, can be expressed as:

$$I = \mathcal{S}t^{0.5} + At \tag{4}$$

#### 2.3.1.3. Wang's model

Wang et al. [28] use k hð Þ¼ <sup>k</sup>sð Þ <sup>h</sup>d=<sup>h</sup> <sup>m</sup> and ð Þ <sup>θ</sup><sup>0</sup> � <sup>θ</sup><sup>r</sup> <sup>=</sup>ð Þ¼ <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>r</sup> ð Þ <sup>h</sup>d=<sup>h</sup> <sup>n</sup> to calculate unsaturated hydraulic conductivity and soil moisture characteristic curves. The relationship between wetting front and cumulative infiltration time can be expressed as:

$$I = x\_f(\theta\_s - \theta\_i) \left( 1 - \frac{n}{m + n - 1} \right) \tag{5}$$

$$q = \frac{1}{\chi\_f} \frac{h\_d}{m - 1} \tag{6}$$

$$\mathbf{x}\_f = \sqrt{\frac{2h\_dk\_s(m+n-1)}{(m^2-1)(\theta\_s-\theta\_i)}}\mathbf{f}^\ddagger \tag{7}$$

Eq. (5), (6), and (7) also can be expressed as:

$$I = A\_1 \mathbf{x}\_f, \qquad \mathbf{i} = A\_2/\mathbf{x}\_f, \qquad \mathbf{x}\_f = A\_3 t^{1/2} \tag{8}$$

where θ is the volume of water content (cm<sup>3</sup>

8 < :

as follows:

saturated soil.

/cm<sup>3</sup>

calculated according to soil bulk density, clay content, quartz, and other mineral volume ratios

For the unsaturated soil, the relationship between λ and K<sup>e</sup> (Kersten) is established based on thermal conductivity λdry (W/(m�K)) of dry soil and thermal conductivity λsat (W/(m�K)) of

And the relationship between K<sup>e</sup> and conventional soil moisture content or saturation Sr

Ke <sup>¼</sup> <sup>0</sup>:7logSr <sup>þ</sup> <sup>1</sup>:0 0:<sup>05</sup> <sup>&</sup>lt; Sr <sup>≤</sup> <sup>0</sup>:<sup>1</sup> logSr þ 1:0 Sr > 0:1

<sup>λ</sup>sat <sup>¼</sup> <sup>λ</sup><sup>1</sup>�<sup>n</sup>

where λ<sup>w</sup> = 0.594 W/(m�K) under the condition 20�C, n is the soil porosity, and λ<sup>s</sup> is obtained by the quartz content (q) of the whole solid, its thermal conductivity being λ<sup>q</sup> = 7.7 W/(m�K),

> <sup>λ</sup>dry <sup>¼</sup> <sup>0</sup>:135r<sup>b</sup> <sup>þ</sup> <sup>64</sup>:<sup>7</sup> 2700 � 0:947r<sup>b</sup>

In order to simplify the calculation of the logarithmic function formula in the Johansen model, Côté and Konrad proposed a new relationship between K<sup>e</sup> and Sr based on the parameter k:

where k is an independent parameter related to the soil texture and its values for coarse sand, small sand, clay, and higher organic matter content are 4.60, 3.25, 1.40, and 1.20, respectively.

1 þ ð Þ k � 1 Sr

Ke <sup>¼</sup> kSr

<sup>s</sup> λ<sup>n</sup>

<sup>D</sup> <sup>¼</sup> <sup>0</sup>:<sup>03</sup> <sup>þ</sup> <sup>0</sup>:1r<sup>2</sup>

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

<sup>A</sup> <sup>¼</sup> <sup>0</sup>:<sup>65</sup> � <sup>0</sup>:78r<sup>b</sup> <sup>þ</sup> <sup>0</sup>:60r<sup>2</sup>

2:6 m<sup>0</sup>:<sup>5</sup> c

C ¼ 1 þ

where mc is the clay content and r<sup>b</sup> is the soil bulk density.

(Sr = θ/θs, where θ<sup>s</sup> is saturated water content) is established:

and thermal conductivity (λ0) of the other minerals is <sup>λ</sup><sup>s</sup> <sup>¼</sup> <sup>λ</sup><sup>q</sup>

3. Improved Johansen model by Côté and Konrad [19].

And a new formula to estimate λdry is given as follows:

(q > 0.2), λ<sup>0</sup> = 3.0 W/(m�K) (q ≤ 0.2).

�

2. Semi-theoretical model of thermal conductivity by Johansen [18].

), and the parameters A, B, C, D, and E can be

http://dx.doi.org/10.5772/intechopen.76280

(14)

135

(16)

(18)

(19)

<sup>b</sup> B ¼ 1:06r<sup>b</sup>

<sup>b</sup> E ¼ 4

<sup>λ</sup> <sup>¼</sup> <sup>λ</sup>sat � <sup>λ</sup>dry � �Ke <sup>þ</sup> <sup>λ</sup>dry (15)

qλ<sup>1</sup>�<sup>q</sup>

<sup>λ</sup>dry <sup>¼</sup> <sup>χ</sup>10�η<sup>n</sup> (20)

<sup>w</sup> (17)

<sup>0</sup> , where, λ<sup>0</sup> = 2.0 W/(m�K)

$$\begin{aligned} A\_1 &= \frac{k\_s h\_d^m h\_s^{1-m}}{m-1} \\ A\_2 &= (\theta\_s - \theta\_i) \left( 1 - \frac{nh\_d^{m-1} h\_s^{m-1}}{m+n-1} \right) \\ A\_3 &= \frac{2k\_s h\_d^m h\_s^{1-m} (m+n-1)}{(m-1)(\theta\_s - \theta) \left( m+n-1-nh\_d^{m-1} h\_s^{m-1} \right)} \end{aligned} \tag{9}$$

Therefore, parameters n, m, and h<sup>d</sup> can be expressed as:

$$m = \sqrt{\frac{\Theta\_s - \Theta\_r}{A\_1 + \Theta\_i - \Theta\_r} - 1}$$

$$h\_d = \frac{A\_2}{anK\_s}$$

$$m = \frac{A\_3(\Theta\_s - \Theta\_r)}{aK\_s h\_d} + \frac{1+n}{n} \tag{10}$$

where θ<sup>i</sup> is the soil water content (cm<sup>3</sup> /cm<sup>3</sup> ), θ<sup>s</sup> is the saturated soil water content (cm<sup>3</sup> /cm<sup>3</sup> ), θ<sup>r</sup> is the residual water content (cm3 /cm<sup>3</sup> ), h<sup>d</sup> is air entry suction (cm), Ks is saturated hydraulic conductivity (cm/min), and n, m is the parameter.

#### 2.3.1.4. Green-Ampt model

Green and Ampt [2] proposed the model. The equation is expressed as:

$$\dot{\mathbf{u}} = \mathbf{K}\_s \left( \mathbf{1} + \frac{h\_0 + h\_f}{z\_f} \right) \tag{11}$$

$$I = (\theta\_s - \theta\_i)z\_f \tag{12}$$

where i is infiltration rate (cm/min), Ks is saturate hydraulic conductivity (cm/min), h0 is ponder depth (cm), h<sup>f</sup> is wetting front suction (cm), I is cumulative infiltration (cm), and θs, θ<sup>i</sup> are saturated water content (cm<sup>3</sup> /cm<sup>3</sup> ) and initial water content (cm<sup>3</sup> /cm<sup>3</sup> ), respectively.

#### 2.3.2. Models of soil heat conductivity

#### 1. Thermal conductivity empirical model by Campbell [16].

Campbell proposed an empirical formula for calculating soil thermal conductivity based on soil texture, bulk density, and volume moisture content, which can be specifically expressed as

$$\lambda = A + B\theta - (A - D)\exp\left[-\left(\mathbb{C}\theta\right)^{\mathbb{E}}\right] \tag{13}$$

where θ is the volume of water content (cm<sup>3</sup> /cm<sup>3</sup> ), and the parameters A, B, C, D, and E can be calculated according to soil bulk density, clay content, quartz, and other mineral volume ratios as follows:

$$\begin{cases} A = 0.65 - 0.78\rho\_b + 0.60\rho\_b^2 & B = 1.06\rho\_b \\ C = 1 + \frac{2.6}{m\_c^{0.5}} & D = 0.03 + 0.1\rho\_b^2 \qquad E = 4 \end{cases} \tag{14}$$

where mc is the clay content and r<sup>b</sup> is the soil bulk density.

Eq. (5), (6), and (7) also can be expressed as:

134 Hydrology of Artificial and Controlled Experiments

<sup>A</sup><sup>1</sup> <sup>¼</sup> kshm

Therefore, parameters n, m, and h<sup>d</sup> can be expressed as:

where θ<sup>i</sup> is the soil water content (cm<sup>3</sup>

conductivity (cm/min), and n, m is the parameter.

is the residual water content (cm3

are saturated water content (cm<sup>3</sup>

2.3.2. Models of soil heat conductivity

2.3.1.4. Green-Ampt model

<sup>d</sup> <sup>h</sup><sup>1</sup>�<sup>m</sup> s m � 1 <sup>A</sup><sup>2</sup> <sup>¼</sup> ð Þ <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>i</sup> <sup>1</sup> � nhm�<sup>1</sup>

<sup>A</sup><sup>3</sup> <sup>¼</sup> <sup>2</sup>kshm

n ¼

r

hd<sup>¼</sup> <sup>d</sup> <sup>A</sup><sup>2</sup> anKs

/cm<sup>3</sup>

i ¼ Ks 1 þ

where i is infiltration rate (cm/min), Ks is saturate hydraulic conductivity (cm/min), h0 is ponder depth (cm), h<sup>f</sup> is wetting front suction (cm), I is cumulative infiltration (cm), and θs, θ<sup>i</sup>

Campbell proposed an empirical formula for calculating soil thermal conductivity based on soil texture, bulk density, and volume moisture content, which can be specifically expressed as

<sup>λ</sup> <sup>¼</sup> <sup>A</sup> <sup>þ</sup> <sup>B</sup><sup>θ</sup> � ð Þ <sup>A</sup> � <sup>D</sup> exp �ð Þ <sup>C</sup><sup>θ</sup> <sup>E</sup> h i

/cm<sup>3</sup>

Green and Ampt [2] proposed the model. The equation is expressed as:

/cm<sup>3</sup>

1. Thermal conductivity empirical model by Campbell [16].

<sup>m</sup> <sup>¼</sup> <sup>A</sup>3ð Þ <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>r</sup> aKshd

I ¼ A1xf , i ¼ A2=xf , xf ¼ A3t

<sup>d</sup> <sup>h</sup><sup>1</sup>�<sup>m</sup>

ð Þ <sup>m</sup> � <sup>1</sup> ð Þ <sup>θ</sup><sup>s</sup> � <sup>θ</sup> <sup>m</sup> <sup>þ</sup> <sup>n</sup> � <sup>1</sup> � nhm�<sup>1</sup>

<sup>d</sup> hm�<sup>1</sup> s m þ n � 1 � �

> ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>θ</sup><sup>s</sup> � <sup>θ</sup><sup>r</sup> A<sup>1</sup> þ θ<sup>i</sup> � θ<sup>r</sup>

<sup>s</sup> ð Þ m þ n � 1

� �

� 1

þ 1 þ n

h<sup>0</sup> þ hf zf � �

) and initial water content (cm<sup>3</sup>

<sup>d</sup> hm�<sup>1</sup> s

), θ<sup>s</sup> is the saturated soil water content (cm<sup>3</sup>

), h<sup>d</sup> is air entry suction (cm), Ks is saturated hydraulic

I ¼ ð Þ θ<sup>s</sup> � θ<sup>i</sup> zf (12)

/cm<sup>3</sup>

), respectively.

<sup>1</sup>=<sup>2</sup> (8)

<sup>n</sup> (10)

(9)

/cm<sup>3</sup> ), θ<sup>r</sup>

(11)

(13)

2. Semi-theoretical model of thermal conductivity by Johansen [18].

For the unsaturated soil, the relationship between λ and K<sup>e</sup> (Kersten) is established based on thermal conductivity λdry (W/(m�K)) of dry soil and thermal conductivity λsat (W/(m�K)) of saturated soil.

$$
\lambda = (\lambda\_{\rm sat} - \lambda\_{\rm dry})\mathbf{K}\_e + \lambda\_{\rm dry} \tag{15}
$$

And the relationship between K<sup>e</sup> and conventional soil moisture content or saturation Sr (Sr = θ/θs, where θ<sup>s</sup> is saturated water content) is established:

$$K\_{\ell} = \begin{cases} 0.7 \text{log} \text{S}\_{r} + 1.0 & 0.05 < \text{S}\_{r} \le 0.1\\ \log \text{S}\_{r} + 1.0 & \text{S}\_{r} > 0.1 \end{cases} \tag{16}$$

$$
\lambda\_{\rm sat} = \lambda\_s^{1-n} \lambda\_w^n \tag{17}
$$

where λ<sup>w</sup> = 0.594 W/(m�K) under the condition 20�C, n is the soil porosity, and λ<sup>s</sup> is obtained by the quartz content (q) of the whole solid, its thermal conductivity being λ<sup>q</sup> = 7.7 W/(m�K), and thermal conductivity (λ0) of the other minerals is <sup>λ</sup><sup>s</sup> <sup>¼</sup> <sup>λ</sup><sup>q</sup> qλ<sup>1</sup>�<sup>q</sup> <sup>0</sup> , where, λ<sup>0</sup> = 2.0 W/(m�K) (q > 0.2), λ<sup>0</sup> = 3.0 W/(m�K) (q ≤ 0.2).

$$
\lambda\_{dry} = \frac{0.135\rho\_b + 64.7}{2700 - 0.947\rho\_b} \tag{18}
$$

3. Improved Johansen model by Côté and Konrad [19].

In order to simplify the calculation of the logarithmic function formula in the Johansen model, Côté and Konrad proposed a new relationship between K<sup>e</sup> and Sr based on the parameter k:

$$K\_{\ell} = \frac{kS\_r}{1 + (k-1)S\_r} \tag{19}$$

where k is an independent parameter related to the soil texture and its values for coarse sand, small sand, clay, and higher organic matter content are 4.60, 3.25, 1.40, and 1.20, respectively. And a new formula to estimate λdry is given as follows:

$$
\lambda\_{dry} = \chi 10^{-\eta n} \tag{20}
$$

where χ (W/(m�K)) and η are parameters that are affected by particle traits. The χ and η values for crushed rock, mineral soil, and soil with high organic matter were 1.70 and 1.80, 0.75 and 1.2, and 0.30 and 0.87, respectively.

#### 4. Improved Johansen model by Lu and Ren [20].

In order to make the Johansen model more suitable for calculating the thermal conductivity under the condition of low soil water content, Lu and Ren proposed a new exponential function expression of K<sup>e</sup> about Sr:

$$K\_t = \exp\{a\left[1 - S\_r^{a-1.33}\right]\}\tag{21}$$

er <sup>¼</sup> ap rb

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

<sup>a</sup> <sup>¼</sup> <sup>e</sup>

where e<sup>r</sup> is the raindrop-induced water transfer rate, a is the detachability of the bare soil (g/cm) [29], e is the rainfall-induced soil detachment per unit soil area (g/cm), r is a constant parameter, r<sup>b</sup> is the bulk density of the dry soil (g/cm), p is the rainfall intensity (mm/min), and

r tð Þð Þ Ptp þ rbθoHo

In the horizontal one-dimensional suction process, the soil water content increases with times, along with cumulative infiltration. However, the cumulative infiltration amount is different at different negative hydraulic heads and at various soil textures in the same infiltration time. From Figure 5, it can be seen that under different negative hydraulic head conditions, the change of cumulative is the largest in Loessal soil, followed by red glue soil and black loessial soil, with Lou soil showing the smallest cumulative change in the same infiltration time.

It can be seen from Figure 6 that the cumulative infiltration capacity with infiltration time is almost the same for each soil. The regulation in infiltration gradually decreases with a negative hydraulic head increase.. Among them, the most significant change was observed in Loessial soil, and there is a large difference between �9 and � 12 cm; the reason is the soil porosity ratio. Lou soil also shows a large difference. From �2.5 to �18 cm, red glue soil, black Loessial soil, and Lou soil show no significant difference, because no difference in porosity was observed between these two hydraulic heads. As for various soil textures, there are great differences in soil moisture absorption characteristics at different negative hydraulic head conditions, and the

In order to obtain a negative pressure suction effect on soil infiltration characteristics of quantitative analysis, we use the Kostiakov infiltration equation to fit the measured data. The Kostiakov model could fit the cumulative infiltration and infiltration time very well, which were shown

From Table 2, it could be found that the correlation coefficient R<sup>2</sup> are all lager than 0.99 which indicated that the relationship between cumulative infiltration and time all have followed a power function under different negative hydraulic heads. For different soil textures, the order of coefficient a is Loessial soil >red glue soil>dark loessial soil>Lou soil, and index b had no

t

θ is the soil water content (%) [30, 31]. Eqs. (23) and (24) can be combined to give:

3. Results and discussions

3.1.1. Cumulative infiltration with times

3.1. The infiltration in the horizonal soil column

lighter the soil texture, the more is the difference.

Table 2.

c tðÞ¼ apθCoHo

θ (24)

http://dx.doi.org/10.5772/intechopen.76280

137

<sup>p</sup><sup>r</sup> (25)

<sup>b</sup> (26)

where α is the parameter determined by the soil texture, and for coarse soils with sand content larger than 40% and fine soils with sand content of less than 40%, α = 0.96, 0.27. 1.33 refers to the shape parameters. A new formula is given for mineral soil as follows:

$$
\lambda\_{dry} = -an + b \tag{22}
$$

where a and b are the empirical coefficients; when 0.2 < n < 0.6, the value is: 0.56, 0.51.

#### 2.3.3. Models of nutrient runoff on the slope

To better understand the factors affecting the loss of solutes to the runoff, we applied our experimental data to the model developed by Wang et al. [25]. This model is described and justified in full detail in the publication cited above and is only briefly outlined here. The model is based on a soil water system that is divided into three vertically distributed horizontal layers: runoff or water ponding on the surface; an exchange layer below that; and the underlying soil. The variation in solute mass in the exchange layer changes over time and can be modeled using a power function. The transport of solutes from the exchange layer to the surface runoff is assumed to be dependent on the mass exchange rate. The model can be expressed as:

$$\sigma(t) = k\_m \frac{\mathbb{C}\_o \rho\_b H\_o}{r(t)(Ptp + \rho\_b \Theta\_o H\_o)} t^b \tag{23}$$

where c(t) represents the solute concentration (mg/L) in the runoff, km is the exchange rate, co denotes the initial solute concentration in the surface soil (g/g), r<sup>b</sup> is the soil's dry bulk density (g/ cm), H<sup>o</sup> is the depth of the exchange layer (cm), r(t) is the runoff volume (L), p is the rainfall intensity (cm/min), t<sup>p</sup> is the time between the initiation of rainfall and the formation of the runoff (min), θ<sup>o</sup> is the initial soil moisture content (%), t is time (min), and b is an empirical parameter.

We adopted a new model, in which the presumed exchange layer is replaced by a mixing zone, which can be regarded as an extension of the deposited layer or "shield" concept presented by Hairsine and Rose [26]. Assuming that the exchange rate is controlled by raindrop splash and that the effects of diffusion can be neglected, we replaced the exchange rate km with the variable er, developed by Gao et al. [27]. This substitution obviates the need to calibrate km. This new variable is the rate at which soil water is ejected from the soil during rainfall:

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 137

$$
\sigma\_r = \frac{ap}{\rho\_b} \Theta \tag{24}
$$

$$a = \frac{e}{p^{\rho}}\tag{25}$$

where e<sup>r</sup> is the raindrop-induced water transfer rate, a is the detachability of the bare soil (g/cm) [29], e is the rainfall-induced soil detachment per unit soil area (g/cm), r is a constant parameter, r<sup>b</sup> is the bulk density of the dry soil (g/cm), p is the rainfall intensity (mm/min), and θ is the soil water content (%) [30, 31]. Eqs. (23) and (24) can be combined to give:

$$\mathcal{L}(t) = \frac{ap\theta C\_oH\_o}{r(t)(Ptp + \rho\_b\theta\_oH\_o)}t^b \tag{26}$$

#### 3. Results and discussions

where χ (W/(m�K)) and η are parameters that are affected by particle traits. The χ and η values for crushed rock, mineral soil, and soil with high organic matter were 1.70 and 1.80, 0.75 and

In order to make the Johansen model more suitable for calculating the thermal conductivity under the condition of low soil water content, Lu and Ren proposed a new exponential

Ke <sup>¼</sup> exp <sup>α</sup> <sup>1</sup> � <sup>S</sup><sup>α</sup>�1:<sup>33</sup>

where α is the parameter determined by the soil texture, and for coarse soils with sand content larger than 40% and fine soils with sand content of less than 40%, α = 0.96, 0.27. 1.33 refers to

To better understand the factors affecting the loss of solutes to the runoff, we applied our experimental data to the model developed by Wang et al. [25]. This model is described and justified in full detail in the publication cited above and is only briefly outlined here. The model is based on a soil water system that is divided into three vertically distributed horizontal layers: runoff or water ponding on the surface; an exchange layer below that; and the underlying soil. The variation in solute mass in the exchange layer changes over time and can be modeled using a power function. The transport of solutes from the exchange layer to the surface runoff is

> CorbHo r tð Þð Þ Ptp þ rbθoHo

where c(t) represents the solute concentration (mg/L) in the runoff, km is the exchange rate, co denotes the initial solute concentration in the surface soil (g/g), r<sup>b</sup> is the soil's dry bulk density (g/ cm), H<sup>o</sup> is the depth of the exchange layer (cm), r(t) is the runoff volume (L), p is the rainfall intensity (cm/min), t<sup>p</sup> is the time between the initiation of rainfall and the formation of the runoff (min), θ<sup>o</sup> is the initial soil moisture content (%), t is time (min), and b is an empirical parameter. We adopted a new model, in which the presumed exchange layer is replaced by a mixing zone, which can be regarded as an extension of the deposited layer or "shield" concept presented by Hairsine and Rose [26]. Assuming that the exchange rate is controlled by raindrop splash and that the effects of diffusion can be neglected, we replaced the exchange rate km with the variable er, developed by Gao et al. [27]. This substitution obviates the need to calibrate km.

This new variable is the rate at which soil water is ejected from the soil during rainfall:

t

<sup>b</sup> (23)

where a and b are the empirical coefficients; when 0.2 < n < 0.6, the value is: 0.56, 0.51.

assumed to be dependent on the mass exchange rate. The model can be expressed as:

c tðÞ¼ km

the shape parameters. A new formula is given for mineral soil as follows:

r

(21)

λdry ¼ �an þ b (22)

1.2, and 0.30 and 0.87, respectively.

136 Hydrology of Artificial and Controlled Experiments

function expression of K<sup>e</sup> about Sr:

2.3.3. Models of nutrient runoff on the slope

4. Improved Johansen model by Lu and Ren [20].

#### 3.1. The infiltration in the horizonal soil column

#### 3.1.1. Cumulative infiltration with times

In the horizontal one-dimensional suction process, the soil water content increases with times, along with cumulative infiltration. However, the cumulative infiltration amount is different at different negative hydraulic heads and at various soil textures in the same infiltration time. From Figure 5, it can be seen that under different negative hydraulic head conditions, the change of cumulative is the largest in Loessal soil, followed by red glue soil and black loessial soil, with Lou soil showing the smallest cumulative change in the same infiltration time.

It can be seen from Figure 6 that the cumulative infiltration capacity with infiltration time is almost the same for each soil. The regulation in infiltration gradually decreases with a negative hydraulic head increase.. Among them, the most significant change was observed in Loessial soil, and there is a large difference between �9 and � 12 cm; the reason is the soil porosity ratio. Lou soil also shows a large difference. From �2.5 to �18 cm, red glue soil, black Loessial soil, and Lou soil show no significant difference, because no difference in porosity was observed between these two hydraulic heads. As for various soil textures, there are great differences in soil moisture absorption characteristics at different negative hydraulic head conditions, and the lighter the soil texture, the more is the difference.

In order to obtain a negative pressure suction effect on soil infiltration characteristics of quantitative analysis, we use the Kostiakov infiltration equation to fit the measured data. The Kostiakov model could fit the cumulative infiltration and infiltration time very well, which were shown Table 2.

From Table 2, it could be found that the correlation coefficient R<sup>2</sup> are all lager than 0.99 which indicated that the relationship between cumulative infiltration and time all have followed a power function under different negative hydraulic heads. For different soil textures, the order of coefficient a is Loessial soil >red glue soil>dark loessial soil>Lou soil, and index b had no

Figure 5. Relationship between cumulated infiltration and measured infiltration time at different hydraulic heads with various soils. (a) �2.5 cm; (b) �6 cm; (c) �9 cm; (d) �12 cm; (e) �15 cm; (f) �18 cm.

significant changing tendency. For each soil, the parameter increased with the increase of

Figure 6. Relationship between cumulated infiltration and measured infiltration time at different negative hydraulic

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

139

According to horizontal one-dimensional experiments, we can easily obtain the soil sorptivity with the analysis to cumulative infiltration change with t1/2. In Figure 7, it shows the change processes of cumulative infiltration with t1/2 under four negative hydraulic heads. Sub-graphs a–f are Loessal soils, sub-graphs g–l are red glue soils, sub-graphs m–r are dark Loessial soils,

As shown in Figure 7, the cumulative infiltrations of four kinds of soils had a linear relationship with t1/2. We use a linear function to describe the curves, and the results were shown in Table 3.

negative hydraulic heads, while for parameter b, the opposite is true.

heads. (a) Loessal soil; (b) red glue soil; (c) dark loessial soil; (d) Lou soil.

sub-graph s–x are Lou soils.

3.1.2. Determining the soil sorptivity based on horizonal one-dimensional experiments

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 139

Figure 6. Relationship between cumulated infiltration and measured infiltration time at different negative hydraulic heads. (a) Loessal soil; (b) red glue soil; (c) dark loessial soil; (d) Lou soil.

significant changing tendency. For each soil, the parameter increased with the increase of negative hydraulic heads, while for parameter b, the opposite is true.

#### 3.1.2. Determining the soil sorptivity based on horizonal one-dimensional experiments

Figure 5. Relationship between cumulated infiltration and measured infiltration time at different hydraulic heads with

various soils. (a) �2.5 cm; (b) �6 cm; (c) �9 cm; (d) �12 cm; (e) �15 cm; (f) �18 cm.

138 Hydrology of Artificial and Controlled Experiments

According to horizontal one-dimensional experiments, we can easily obtain the soil sorptivity with the analysis to cumulative infiltration change with t1/2. In Figure 7, it shows the change processes of cumulative infiltration with t1/2 under four negative hydraulic heads. Sub-graphs a–f are Loessal soils, sub-graphs g–l are red glue soils, sub-graphs m–r are dark Loessial soils, sub-graph s–x are Lou soils.

As shown in Figure 7, the cumulative infiltrations of four kinds of soils had a linear relationship with t1/2. We use a linear function to describe the curves, and the results were shown in Table 3.


http://dx.doi.org/10.5772/intechopen.76280

141

Table 2. Fitted value for parameters of I-t at different negative hydraulic heads.

The regression coefficients R<sup>2</sup> in Table 3 were all above 0.95, which indicated that Philip equation can describe the infiltration rule very well under different negative hydraulic heads. Meanwhile, the soil sorptivity decreased with increasing soil viscidity (Loessal soil > red glue soil >dark Loessal soil > Lou soil). And beyond that, soil sorptivity decreases with the increasing negative hydraulic head.

#### 3.1.3. Determining the parameters using Wang's proposed equation

We used three kinds of textured soils (red glue soil, dark Loessial soil, and Lou soil) for horizontal one-dimensional infiltration experiments. The length of the soil column is 50 cm. The upper boundary was a constant hydraulic head (i.e., when x = 0, the hydraulic head was designed as a different negative hydraulic head). The hydraulic heads of red glue soil were � 21 cm and � 30 cm. The hydraulic heads of dark Loessial soil were � 18 cm and � 24 cm. The hydraulic heads of Lou soil were � 21 cm and � 34 cm. The lower boundary condition was free discharge. The duration of the experiment was 810 min. The saturated water content, the retention water content, and the initial water content were measured, noting down the changes of the cumulative infiltration with time. The parameters in the Brook-Corey model can be determined by MATLAB programming based on the experimental data. Figure 8 shows the relationship between cumulative infiltration and wetting front under different negative hydraulic head conditions.

As shown in Table 4, there is a good linear relationship between cumulative infiltration and wetting front which is in agreement with the theoretical derivation.

140 Hydrology of Artificial and Controlled Experiments Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 141

The regression coefficients R<sup>2</sup> in Table 3 were all above 0.95, which indicated that Philip equation can describe the infiltration rule very well under different negative hydraulic heads. Meanwhile, the soil sorptivity decreased with increasing soil viscidity (Loessal soil > red glue soil >dark Loessal soil > Lou soil). And beyond that, soil sorptivity decreases with the increas-

Loessal soil A 1.3548 1.1979 0.9380 0.5566 0.4585 0.7439

Red glue soil A 0.353 0.3167 0.2844 0.3820 0.2095 0.2906

Dark loessial soil A 0.2343 0.1942 0.1559 0.1631 0.1421 0.1017

Lou soil A 0.1718 0.1272 0.1244 0.1149 0.073 0.0934

�2.5 �6 �9 �12 �15 �18

B 0.5949 0.6086 0.6502 0.7239 0.7507 0.6453 R<sup>2</sup> 0.9939 0.9937 0.9951 0.9899 0.9847 0.9922

B 0.6025 0.6105 0.6255 0.5787 0.6727 0.6129 R<sup>2</sup> 0.9954 0.9956 0.9956 0.9921 0.9857 0.9925

B 0.6106 0.6418 0.6551 0.6529 0.6666 0.7125 R<sup>2</sup> 0.9941 0.9926 0.9964 0.9889 0.9938 0.9923

B 0.6198 0.6672 0.6561 0.6579 0.744 0.6903 R<sup>2</sup> 0.9976 0.9966 0.9951 0.9985 0.9947 0.9976

We used three kinds of textured soils (red glue soil, dark Loessial soil, and Lou soil) for horizontal one-dimensional infiltration experiments. The length of the soil column is 50 cm. The upper boundary was a constant hydraulic head (i.e., when x = 0, the hydraulic head was designed as a different negative hydraulic head). The hydraulic heads of red glue soil were � 21 cm and � 30 cm. The hydraulic heads of dark Loessial soil were � 18 cm and � 24 cm. The hydraulic heads of Lou soil were � 21 cm and � 34 cm. The lower boundary condition was free discharge. The duration of the experiment was 810 min. The saturated water content, the retention water content, and the initial water content were measured, noting down the changes of the cumulative infiltration with time. The parameters in the Brook-Corey model can be determined by MATLAB programming based on the experimental data. Figure 8 shows the relationship between cumu-

lative infiltration and wetting front under different negative hydraulic head conditions.

wetting front which is in agreement with the theoretical derivation.

As shown in Table 4, there is a good linear relationship between cumulative infiltration and

ing negative hydraulic head.

3.1.3. Determining the parameters using Wang's proposed equation

Table 2. Fitted value for parameters of I-t at different negative hydraulic heads.

Soil textural Parameters Hydraulic head (cm)

http://dx.doi.org/10.5772/intechopen.76280

 Hydrology of Artificial and Controlled Experiments Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 

Substitution of A1 and A2 (listed in Table 5) into Eq. (10) yields the hydrodynamic parameters

Hydraulic head (cm) Loessal soil Red glue soil Dark loessial soil Lou soil

�2.5 1.9904 0.994 0.6448 0.9903 0.4436 0.993 0.3466 0.9882 �6 1.8747 0.9903 0.6083 0.9878 0.4423 0.9893 0.344 0.9759 �9 1.7482 0.9836 0.5968 0.9856 0.3902 0.9814 0.3112 0.9843 �12 1.4331 0.9706 0.6008 0.9959 0.4012 0.9838 0.2947 0.9777 �15 1.3145 0.9686 0.5695 0.986 0.3776 0.9832 0.2875 0.9688 �18 1.3654 0.9866 0.5597 0.9918 0.3552 0.9746 0.2918 0.9684

S R<sup>2</sup> S R<sup>2</sup> S R<sup>2</sup> S R<sup>2</sup>

http://dx.doi.org/10.5772/intechopen.76280

145

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

The soil water characteristic curves can be easily obtained by the values in Table 5. Comparing the calculated results and the experimental results (determined by centrifuge), the results were listed in Figure 9. As shown in Figure 9, the calculated data concur with experimental data. The results indicated that the parameters in the Brooks-Core model can be accurately and easy

Soil thermal conductivity reflects the size of soil thermal conductivity, and soil texture has a certain impact on thermal conductivity. According to the principle of thermal pulse probe, soil thermal parameters were measured and soil thermal conductivity was calculated. Figure 10 shows the curves of soil thermal conductivity with soil moisture content in four experimental sites of Shenmu (sand), Ansai, Yichuan, and Changwu. It can be seen from Figure 10 that soil thermal conductivity increases rapidly with the increase of water content when the soil water

increasing trend of soil thermal conductivity is relatively reduced. Under the same moisture content, the trend of soil thermal conductivity is as follows: Shenmu sand soil>Ansai sandy loam soil>Yichuan clay loam soil> Changwu silty loam soil. So we can see that the higher the sand

A study by Lu and Ren et al. showed that the soil can be divided into two categories according to sand content of the soil: It is coarse soil when the sand content is more than 40%, and Sr = 0.3, Ke-S<sup>r</sup> relationship curve of the coarse soil is divided into two linear ranges; when the sand content is less than 40%, it is fine soil, it is composed of S<sup>r</sup> = 0.13 and S<sup>r</sup> = 0.30, and the Ke-S<sup>r</sup> relationship curve of the fine soil is divided into three linear intervals; Figure 11 shows the Ke-S<sup>r</sup> curves of the normal form of soil thermal conductivity. The soil samples of Changwu and Ankang belong to fine soil, Shenmu, Mizhi, Ansai, Shangnan, Yichuan, Luochuan, and

content, the lower the silt content, and the greater the soil thermal conductivity [22].

. When the soil water content is higher than 0.13 cm3

/cm3 , the

(in the Brooks-Core model). The results are listed in Table 5.

Table 3. The suction of four soil at different negative hydraulic heads.

3.2. Analysis of soil thermal conductivity characteristics

/cm3

computed by the new method.

content is lower than 0.13 cm3

Figure 7. Relationship between cumulated infiltration and square root of measured infiltration time at different negative heads. A. Loessal soil; B. Red glue soil; C. Dark loessial soil; D. Lou soil.

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 145


Table 3. The suction of four soil at different negative hydraulic heads.

Substitution of A1 and A2 (listed in Table 5) into Eq. (10) yields the hydrodynamic parameters (in the Brooks-Core model). The results are listed in Table 5.

The soil water characteristic curves can be easily obtained by the values in Table 5. Comparing the calculated results and the experimental results (determined by centrifuge), the results were listed in Figure 9. As shown in Figure 9, the calculated data concur with experimental data. The results indicated that the parameters in the Brooks-Core model can be accurately and easy computed by the new method.

#### 3.2. Analysis of soil thermal conductivity characteristics

Figure 7. Relationship between cumulated infiltration and square root of measured infiltration time at different negative

heads. A. Loessal soil; B. Red glue soil; C. Dark loessial soil; D. Lou soil.

144 Hydrology of Artificial and Controlled Experiments

Soil thermal conductivity reflects the size of soil thermal conductivity, and soil texture has a certain impact on thermal conductivity. According to the principle of thermal pulse probe, soil thermal parameters were measured and soil thermal conductivity was calculated. Figure 10 shows the curves of soil thermal conductivity with soil moisture content in four experimental sites of Shenmu (sand), Ansai, Yichuan, and Changwu. It can be seen from Figure 10 that soil thermal conductivity increases rapidly with the increase of water content when the soil water content is lower than 0.13 cm3 /cm3 . When the soil water content is higher than 0.13 cm3 /cm3 , the increasing trend of soil thermal conductivity is relatively reduced. Under the same moisture content, the trend of soil thermal conductivity is as follows: Shenmu sand soil>Ansai sandy loam soil>Yichuan clay loam soil> Changwu silty loam soil. So we can see that the higher the sand content, the lower the silt content, and the greater the soil thermal conductivity [22].

A study by Lu and Ren et al. showed that the soil can be divided into two categories according to sand content of the soil: It is coarse soil when the sand content is more than 40%, and Sr = 0.3, Ke-S<sup>r</sup> relationship curve of the coarse soil is divided into two linear ranges; when the sand content is less than 40%, it is fine soil, it is composed of S<sup>r</sup> = 0.13 and S<sup>r</sup> = 0.30, and the Ke-S<sup>r</sup> relationship curve of the fine soil is divided into three linear intervals; Figure 11 shows the Ke-S<sup>r</sup> curves of the normal form of soil thermal conductivity. The soil samples of Changwu and Ankang belong to fine soil, Shenmu, Mizhi, Ansai, Shangnan, Yichuan, Luochuan, and

Soil textural Hydraulic head (cm) A1 R<sup>2</sup> Pressure head (cm) A1 R<sup>2</sup> Red glue soil �21 0.552 0.9997 �30 0.5192 0.9995 Dark loessial soil �18 0.5539 0.9997 �24 0.5059 0.9994 Lou soil �21 0.5356 0.9995 �24 0.5195 0.9966

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

147

Soil textural n Mh<sup>d</sup> m/n Red glue soil 0.17 2.51 26 14.65 Dark loessial soil 0.32 2.94 40.11 9.18 Lou soil 0.23 2.69 54.06 11.69

Figure 9. Compared between observed and calculated soil water characteristic curve. a. Red glue soil; b. Dark loessial

Table 4. Relation fitting values of cumulated infiltration and wetting front.

Table 5. Parameters calculation of three soils.

Figure 10. Trend of soil thermal conductivity with water content.

soil; c. Lou soil.

Figure 8. Relation curve between cumulated infiltration and wetting front. A. Red glue soil; B. Dark loessial soil; C. Lou soil.

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 147


Table 4. Relation fitting values of cumulated infiltration and wetting front.


Table 5. Parameters calculation of three soils.

Figure 9. Compared between observed and calculated soil water characteristic curve. a. Red glue soil; b. Dark loessial soil; c. Lou soil.

Figure 10. Trend of soil thermal conductivity with water content.

Figure 8. Relation curve between cumulated infiltration and wetting front. A. Red glue soil; B. Dark loessial soil; C. Lou

soil.

146 Hydrology of Artificial and Controlled Experiments

Figure 11. Thermal conductivity as a function of water content as indicated by the normalized form (the Kersten number K<sup>e</sup> - the degree of saturation Sr).

Zhangye kind of matter is coarse soil. It can be seen from the figure that the K<sup>e</sup> value of fine soil is obviously smaller than the K<sup>e</sup> value of coarse soil when 0.2 < S<sup>r</sup> < 0.6.

> the following four soils: Shenmu sand soil, Ansai sandy loam soil, Changwu silty loam soil, and Ankang silty clay loam soil. The model's calculated values and measured values are shown in Figure 13. It can be seen from the figure that the calculated values of Johansen model are significantly smaller than the measured values, the calculation error is larger, the coefficient of determination R<sup>2</sup> is in the range of 0.656–0.827, the root mean square error (RMSE) is in the range of 0.0848–0.2548, and the relative error R<sup>e</sup> is in the range of 10.32–20.41%. For fine soil, the Côté-Konrad model and Lu-Ren model have a good fitting effect on soil thermal conductivity and the precision is high. Where the variation coefficient of R<sup>2</sup> is in the range from 0.842 to 0.940, the variation range of RMSE is from 0.0810 to 0.1208, the relative error R<sup>e</sup> is in the range from 9.67 to 10.57%. The coefficient of determination R2 of the Lu-Ren model is in the range from 0.874 to 0.937, RMSE varied from 0.0725 to 0.1238, and the relative error R<sup>e</sup> varied from 8.28% to 9.91%. For coarse soil (sand content greater than 40%), the Côté-Konrad model and Lu-Ren model can still well fit soil thermal conductivity when the saturation Sr < 50%, but the prediction accuracy of the model is poor, and the calculated value is obviously smaller than the measured value when the saturation S<sup>r</sup> > 50%. This phenomenon may be due to the large soil voids and the weak water-holding capacity, resulting in the measured value of water content being lower.

> Figure 12. Comparison of soil thermal conductivity values calculated by Campbell model with measured. (a) Shenmu

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

149

sand soil, (b) Ansai sandy loam soil, (c) Changwu silty loam soil and (d) Ankang silty clay loam soil.

The comparison between the calculated values of Campbell model, Johansen model, Côté-Konrad model, and Lu-Ren model and the measured values of the thermal pulses show that the soil thermal conductivity is closely related to the soil particle composition, organic matter

3.2.2. The improved Côté-Konrad model and the improved Lu-Ren model

#### 3.2.1. Accuracy analysis of the soil thermal conductivity model

#### 3.2.1.1. Campbell model

The sandy soil and sandy loam soil of coarse soil in Shenmu and Ansai, silty loam soil and silty clay loam soil of fine soil in ChangWu and Ankang, respectively, were selected. The thermal conductivities of these four soils were calculated by the Campbell model, and the results were shown in Figure 12. According to the statistical analysis, it can be seen that the difference between the calculated value and the measured value of the heat pulse is small when the water content of the soil is less than 0.20 cm<sup>3</sup> /cm<sup>3</sup> , and the relative error (Re) of Shenmu sand soil and Ansai sandy loam soil are 13.51 and 9.56%, respectively; When the soil water content is higher than 0.20 cm<sup>3</sup> /cm<sup>3</sup> , the measured value of the heat pulse is larger than the calculated value, and the R<sup>e</sup> of Shenmu sand soil and Ansai sandy loam soil are 19.40 and 13.38%, respectively; the larger the volume of moisture content, the greater the difference; relative to the thermal pulse's measured value, the calculation of the coarse soil model is too small for the coarse soil. For the fine soil, when the soil moisture content is less than 0.25 cm3 /cm<sup>3</sup> , the R<sup>e</sup> are 26.29 and 21.19%, respectively, and the measured value of the heat pulse is larger than the calculated value of the model. When the soil moisture content is higher than 0.25 cm<sup>3</sup> /cm<sup>3</sup> , and the R<sup>e</sup> of Shenmu sand soil and Ansai sandy loam soil are 14.15 and 6.60%, respectively, the difference between the calculated value and the measured value of the heat pulse is small. Therefore, the model needs to be improved when calculating the thermal conductivity using the Campbell model [8, 9].

#### 3.2.1.2. Johansen model, Côté-Konrad model and Lu-Ren model

Côté-Konrad model and Lu-Ren model are all semi-theoretical models of thermal conductivity based on Johansen model. The three models are used to calculate soil thermal conductivity for

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 149

Figure 12. Comparison of soil thermal conductivity values calculated by Campbell model with measured. (a) Shenmu sand soil, (b) Ansai sandy loam soil, (c) Changwu silty loam soil and (d) Ankang silty clay loam soil.

the following four soils: Shenmu sand soil, Ansai sandy loam soil, Changwu silty loam soil, and Ankang silty clay loam soil. The model's calculated values and measured values are shown in Figure 13. It can be seen from the figure that the calculated values of Johansen model are significantly smaller than the measured values, the calculation error is larger, the coefficient of determination R<sup>2</sup> is in the range of 0.656–0.827, the root mean square error (RMSE) is in the range of 0.0848–0.2548, and the relative error R<sup>e</sup> is in the range of 10.32–20.41%. For fine soil, the Côté-Konrad model and Lu-Ren model have a good fitting effect on soil thermal conductivity and the precision is high. Where the variation coefficient of R<sup>2</sup> is in the range from 0.842 to 0.940, the variation range of RMSE is from 0.0810 to 0.1208, the relative error R<sup>e</sup> is in the range from 9.67 to 10.57%. The coefficient of determination R2 of the Lu-Ren model is in the range from 0.874 to 0.937, RMSE varied from 0.0725 to 0.1238, and the relative error R<sup>e</sup> varied from 8.28% to 9.91%. For coarse soil (sand content greater than 40%), the Côté-Konrad model and Lu-Ren model can still well fit soil thermal conductivity when the saturation Sr < 50%, but the prediction accuracy of the model is poor, and the calculated value is obviously smaller than the measured value when the saturation S<sup>r</sup> > 50%. This phenomenon may be due to the large soil voids and the weak water-holding capacity, resulting in the measured value of water content being lower.

#### 3.2.2. The improved Côté-Konrad model and the improved Lu-Ren model

Zhangye kind of matter is coarse soil. It can be seen from the figure that the K<sup>e</sup> value of fine soil

Figure 11. Thermal conductivity as a function of water content as indicated by the normalized form (the Kersten number

The sandy soil and sandy loam soil of coarse soil in Shenmu and Ansai, silty loam soil and silty clay loam soil of fine soil in ChangWu and Ankang, respectively, were selected. The thermal conductivities of these four soils were calculated by the Campbell model, and the results were shown in Figure 12. According to the statistical analysis, it can be seen that the difference between the calculated value and the measured value of the heat pulse is small when the water

Ansai sandy loam soil are 13.51 and 9.56%, respectively; When the soil water content is higher

the R<sup>e</sup> of Shenmu sand soil and Ansai sandy loam soil are 19.40 and 13.38%, respectively; the larger the volume of moisture content, the greater the difference; relative to the thermal pulse's measured value, the calculation of the coarse soil model is too small for the coarse soil. For the

respectively, and the measured value of the heat pulse is larger than the calculated value of the

soil and Ansai sandy loam soil are 14.15 and 6.60%, respectively, the difference between the calculated value and the measured value of the heat pulse is small. Therefore, the model needs to be improved when calculating the thermal conductivity using the Campbell model [8, 9].

Côté-Konrad model and Lu-Ren model are all semi-theoretical models of thermal conductivity based on Johansen model. The three models are used to calculate soil thermal conductivity for

, the measured value of the heat pulse is larger than the calculated value, and

/cm<sup>3</sup>

/cm<sup>3</sup>

, and the relative error (Re) of Shenmu sand soil and

, the R<sup>e</sup> are 26.29 and 21.19%,

, and the R<sup>e</sup> of Shenmu sand

/cm<sup>3</sup>

is obviously smaller than the K<sup>e</sup> value of coarse soil when 0.2 < S<sup>r</sup> < 0.6.

3.2.1. Accuracy analysis of the soil thermal conductivity model

fine soil, when the soil moisture content is less than 0.25 cm3

model. When the soil moisture content is higher than 0.25 cm<sup>3</sup>

3.2.1.2. Johansen model, Côté-Konrad model and Lu-Ren model

3.2.1.1. Campbell model

K<sup>e</sup> - the degree of saturation Sr).

148 Hydrology of Artificial and Controlled Experiments

than 0.20 cm<sup>3</sup>

content of the soil is less than 0.20 cm<sup>3</sup>

/cm<sup>3</sup>

The comparison between the calculated values of Campbell model, Johansen model, Côté-Konrad model, and Lu-Ren model and the measured values of the thermal pulses show that the soil thermal conductivity is closely related to the soil particle composition, organic matter

of the parameter fitting are shown in Table 6. The comparison between soil thermal conductivity and the measured values is shown in Figure 14. From the fitting error, it can be seen that the accuracy of the two improved models is not very different, and they have high accuracy. However, it can be seen from Figure 14 that the fitting value is larger than the measured value,

0.9274, and 9.62% for the improved Côté-Konrad model when soil thermal conductivity is less than 0.6 W/(m�K), respectively. For the improved Lu-Ren model, although the fitting value and the measured value are also different, the discrete points in the figure are evenly distributed

Luochuan, and Ankang combined with the model parameter fitting values in Table 6 was predicted. The comparison between the predicted values and the measured values about the different models is shown in Figure 15a–d and Table 7, where, the sand content in four experimental sites was as follows: Shenmu> Shangluo> Luochuan> Ankang; clay content: Shenmu <Shangluo <Luochuan <Ankang; silt content: Shenmu <Shangluo <Luochuan <Ankang. Analysis of the simulation error shows that two improved models can be used to simulate soil thermal conductivity of different soils. For Shenmu sand soil and Ankang silty clay loam soil, the RMSE of the improved Côté-Konrad model is less than 0.1183, the R<sup>2</sup> is greater than 0.9259, and the R<sup>e</sup> is less than 9.47%, which is better than the Côté-Konrad model,

Model Model parameters RMSE R<sup>2</sup> Re/% a<sup>1</sup> / b<sup>1</sup> a<sup>2</sup> / b<sup>2</sup> a<sup>3</sup> / b<sup>3</sup> a<sup>4</sup> / b<sup>4</sup> Improved Côté-Konrad model 4.1381 �0.8413 4.1506 �0.2200 0.0964 0.9274 9.62 Improved Lu-Ren model �0.5863 0.9451 0.1080 0.0567 0.0961 0.9278 9.59 Note: a<sup>i</sup> is the parameter in the improved Côté-Konrad model; b<sup>i</sup> is the parameter in the improved Lu-Ren model; i = 1, 2,

Table 6. Parameters fitted values and errors by improved Côté-Konrad model and improved Lu-Ren model.

Figure 14. Fitted values of soil thermal conductivity by improved Côté-Konrad model and improved Lu-Ren model.

(a) Improved Côté-Konrad model and (b) Improved Lu-Ren model.

,and Re are 0.0961, 0.9278, and 9.59%, respectively.

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

, Re, soil thermal conductivity of four samples of Shenmu (sand), Shangluo,

, and Re are 0.0964,

151

http://dx.doi.org/10.5772/intechopen.76280

and the soil thermal conductivity is more than 1.1 W/(m�K), the RMSE, R2

near the 1:1 line; RMSE, R<sup>2</sup>

According to R<sup>2</sup>

3, 4.

Figure 13. Comparison of soil thermal conductivity values calculated by different models (Johansen model, Côté-Konrad model and Lu-Ren model) with measured. (a) Shenmu sand soil, (b) Ansai sandy loam soil, (c) Changwu silty loam soil and (d) Ankang silty clay loam soil.

content, and bulk density. For different soils with different textures, model parameters are also different. In Johansen model, the parameter λ<sup>s</sup> is related to the quartz content of the whole solid, and the thermal conductivity λdry of the dry soil is related to soil bulk density. In Côté-Konrad model, the parameter k is related to the content of coarse sand, small sand, clay, and organic matter. In Ren model, the parameter α is related to the soil sand content, and the soil is divided into coarse soil and fine soil according to the sand content. Under certain conditions, these models can calculate the soil thermal conductivity more accurately but cannot reflect the effect of soil particle composition and organic matter content on soil thermal conductivity. The improved Côté-Konrad model and the improved Lu-Ren model established the relationship between the model parameters, the composition of the particles, and the content of organic matter, respectively, and can describe the relationship between soil texture and soil thermal conductivity in detail.

In this chapter, the data of the five sites (466 sample points) of Mizhi, Shenmu (sandy loam), Ansai, Yichuan, and Changwu combined with R<sup>2</sup> , Re, and Figure 13, the relationship between the soil texture and the thermal conductivity, is fitting for the improved model, and the results of the parameter fitting are shown in Table 6. The comparison between soil thermal conductivity and the measured values is shown in Figure 14. From the fitting error, it can be seen that the accuracy of the two improved models is not very different, and they have high accuracy. However, it can be seen from Figure 14 that the fitting value is larger than the measured value, and the soil thermal conductivity is more than 1.1 W/(m�K), the RMSE, R2 , and Re are 0.0964, 0.9274, and 9.62% for the improved Côté-Konrad model when soil thermal conductivity is less than 0.6 W/(m�K), respectively. For the improved Lu-Ren model, although the fitting value and the measured value are also different, the discrete points in the figure are evenly distributed near the 1:1 line; RMSE, R<sup>2</sup> ,and Re are 0.0961, 0.9278, and 9.59%, respectively.

According to R<sup>2</sup> , Re, soil thermal conductivity of four samples of Shenmu (sand), Shangluo, Luochuan, and Ankang combined with the model parameter fitting values in Table 6 was predicted. The comparison between the predicted values and the measured values about the different models is shown in Figure 15a–d and Table 7, where, the sand content in four experimental sites was as follows: Shenmu> Shangluo> Luochuan> Ankang; clay content: Shenmu <Shangluo <Luochuan <Ankang; silt content: Shenmu <Shangluo <Luochuan <Ankang. Analysis of the simulation error shows that two improved models can be used to simulate soil thermal conductivity of different soils. For Shenmu sand soil and Ankang silty clay loam soil, the RMSE of the improved Côté-Konrad model is less than 0.1183, the R<sup>2</sup> is greater than 0.9259, and the R<sup>e</sup> is less than 9.47%, which is better than the Côté-Konrad model,


Note: a<sup>i</sup> is the parameter in the improved Côté-Konrad model; b<sup>i</sup> is the parameter in the improved Lu-Ren model; i = 1, 2, 3, 4.

Table 6. Parameters fitted values and errors by improved Côté-Konrad model and improved Lu-Ren model.

content, and bulk density. For different soils with different textures, model parameters are also different. In Johansen model, the parameter λ<sup>s</sup> is related to the quartz content of the whole solid, and the thermal conductivity λdry of the dry soil is related to soil bulk density. In Côté-Konrad model, the parameter k is related to the content of coarse sand, small sand, clay, and organic matter. In Ren model, the parameter α is related to the soil sand content, and the soil is divided into coarse soil and fine soil according to the sand content. Under certain conditions, these models can calculate the soil thermal conductivity more accurately but cannot reflect the effect of soil particle composition and organic matter content on soil thermal conductivity. The improved Côté-Konrad model and the improved Lu-Ren model established the relationship between the model parameters, the composition of the particles, and the content of organic matter, respectively, and can describe the relationship between soil texture and soil thermal

Figure 13. Comparison of soil thermal conductivity values calculated by different models (Johansen model, Côté-Konrad model and Lu-Ren model) with measured. (a) Shenmu sand soil, (b) Ansai sandy loam soil, (c) Changwu silty loam soil

In this chapter, the data of the five sites (466 sample points) of Mizhi, Shenmu (sandy loam),

the soil texture and the thermal conductivity, is fitting for the improved model, and the results

, Re, and Figure 13, the relationship between

conductivity in detail.

and (d) Ankang silty clay loam soil.

150 Hydrology of Artificial and Controlled Experiments

Ansai, Yichuan, and Changwu combined with R<sup>2</sup>

Figure 14. Fitted values of soil thermal conductivity by improved Côté-Konrad model and improved Lu-Ren model. (a) Improved Côté-Konrad model and (b) Improved Lu-Ren model.

In order to further verify whether the improved model can be extended to other soils, soil thermal conductivity of Zhangye samples in Gansu Province is predicted by the improved model. As the soil samples are sandy clay loam soil, the sand content is 60.13%. From the above model comparison analysis, we can see that the improved Côté-Konrad model is better for soil thermal conductivity with higher sand content. Figure 15e and Table 7 show the prediction results of thermal conductivity and the measured values and the simulation error, respectively. Through the error analysis, we can see that the results show that the improved Côté-Konrad model is slightly higher than other three models where the RMSE and R2 of the improved Côté-Konrad model are 0.1026 and 0.9069, respectively, which is slightly higher than other three models, R<sup>e</sup> is 8.15%, slightly lower than the other three models. Therefore, by selecting the appropriate improved model, soil thermal conductivity for different soil textures

Table 7. Soil thermal conductivity simulated values and errors by different soil thermal conductivity models in sampling

Shenmu 0.1208 0.9401 9.67 0.1238 0.937 9.91 0.1183 0.9425 9.47 0.1366 0.9234 10.94 Ankang 0.1088 0.9062 10.57 0.1014 0.9185 9.85 0.0951 0.9259 9.55 0.0986 0.8775 10.94 Shngluo 0.081 0.8422 9.87 0.0725 0.8736 8.83 0.1243 0.8451 13.17 0.0766 0.9412 8.11 Luochuan 0.0946 0.8872 10.5 0.0747 0.9298 8.28 0.1063 0.8514 10.97 0.0815 0.9326 8.21 Zhangye 0.1216 0.8985 8.68 0.1349 0.8911 8.81 0.1026 0.9069 8.15 0.1034 0.9053 8.22

model

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

RMSE R<sup>2</sup> Re/% RMSE R<sup>2</sup> Re/% RMSE R<sup>2</sup> Re/% RMSE R<sup>2</sup> Re/%

Sampling area Côté-Konrad model Lu-Ren model Improved Côté-Konrad

Most of the parameters in our model were measured directly. The depths of the exchange layers, Ho, were measured directly in the soil profiles and are shown in Table 8. Assuming that the exchange layer was saturated when the runoff was generated, the water content θ was

content, θo, was 0.1 g/g except in the experiments where the initial soil moisture content was varied. The rainfall intensity, p, was controlled by a computer. In all experiments except those investigating the influence of varying this parameter, p was fixed at 90 mm/h. The ponding

To determine the relationship between the rainfall-induced soil detachment per unit area, e, and the rainfall intensity, p, we performed a series of rainfall experiments at rainfall intensities of 36 and 90 mm/h with a slope gradient of 5 and a second series of experiments at rainfall

The measured and simulated (assuming r = 2) relationships between the rainfall-induced soil detachment per unit area, e, and the rainfall intensity, p, for a slope gradient of 5 are shown in Figure 16. The model's predictions agree well with the experimental results (R2 > 0.90); this

/cm<sup>3</sup>

. The initial soil moisture

Improved Lu-Ren

153

model

http://dx.doi.org/10.5772/intechopen.76280

, and r(t) was measured directly.

can be calculated accurately.

area.

3.3. Solute transport in runoff by raindrops

assumed to be equal to that in saturated soil, that is, 0.42 cm<sup>3</sup>

times, tp, are shown in Table 8. r<sup>b</sup> has a value of 1.35 g/cm3

intensities of 60 and 84 mm/h with a slope gradient of 15.

Figure 15. Comparison of soil thermal conductivity values predicted by improved Côté-Konrad model and improved Lu-Ren model with measured. (a) Shenmu sand soil, (b) Shangluo loam soil, (c) Luochuan clay loam soil, (d) Ankang silty clay loam soil and (e) Zhangye sandy clay loam soil.

Lu-Ren model, and improved Lu-Ren model. In other words, the improved Côté-Konrad model can be used to simulate the soil thermal conductivity for the soil with high sand content or high silt content. For the Shangnan loam soil and Luochuan clay loam soil, the RMSE of improved Lu-Ren model is less than 0.0815, R2 is greater than 0.9326, and R<sup>e</sup> is less than 8.11%, which are obviously better than the other three models. In other words, the improved Lu-Ren model can be used to simulate soil thermal conductivity.


Table 7. Soil thermal conductivity simulated values and errors by different soil thermal conductivity models in sampling area.

In order to further verify whether the improved model can be extended to other soils, soil thermal conductivity of Zhangye samples in Gansu Province is predicted by the improved model. As the soil samples are sandy clay loam soil, the sand content is 60.13%. From the above model comparison analysis, we can see that the improved Côté-Konrad model is better for soil thermal conductivity with higher sand content. Figure 15e and Table 7 show the prediction results of thermal conductivity and the measured values and the simulation error, respectively. Through the error analysis, we can see that the results show that the improved Côté-Konrad model is slightly higher than other three models where the RMSE and R2 of the improved Côté-Konrad model are 0.1026 and 0.9069, respectively, which is slightly higher than other three models, R<sup>e</sup> is 8.15%, slightly lower than the other three models. Therefore, by selecting the appropriate improved model, soil thermal conductivity for different soil textures can be calculated accurately.

#### 3.3. Solute transport in runoff by raindrops

Lu-Ren model, and improved Lu-Ren model. In other words, the improved Côté-Konrad model can be used to simulate the soil thermal conductivity for the soil with high sand content or high silt content. For the Shangnan loam soil and Luochuan clay loam soil, the RMSE of improved Lu-Ren model is less than 0.0815, R2 is greater than 0.9326, and R<sup>e</sup> is less than 8.11%, which are obviously better than the other three models. In other words, the improved Lu-Ren

Figure 15. Comparison of soil thermal conductivity values predicted by improved Côté-Konrad model and improved Lu-Ren model with measured. (a) Shenmu sand soil, (b) Shangluo loam soil, (c) Luochuan clay loam soil, (d) Ankang silty

model can be used to simulate soil thermal conductivity.

clay loam soil and (e) Zhangye sandy clay loam soil.

152 Hydrology of Artificial and Controlled Experiments

Most of the parameters in our model were measured directly. The depths of the exchange layers, Ho, were measured directly in the soil profiles and are shown in Table 8. Assuming that the exchange layer was saturated when the runoff was generated, the water content θ was assumed to be equal to that in saturated soil, that is, 0.42 cm<sup>3</sup> /cm<sup>3</sup> . The initial soil moisture content, θo, was 0.1 g/g except in the experiments where the initial soil moisture content was varied. The rainfall intensity, p, was controlled by a computer. In all experiments except those investigating the influence of varying this parameter, p was fixed at 90 mm/h. The ponding times, tp, are shown in Table 8. r<sup>b</sup> has a value of 1.35 g/cm3 , and r(t) was measured directly.

To determine the relationship between the rainfall-induced soil detachment per unit area, e, and the rainfall intensity, p, we performed a series of rainfall experiments at rainfall intensities of 36 and 90 mm/h with a slope gradient of 5 and a second series of experiments at rainfall intensities of 60 and 84 mm/h with a slope gradient of 15.

The measured and simulated (assuming r = 2) relationships between the rainfall-induced soil detachment per unit area, e, and the rainfall intensity, p, for a slope gradient of 5 are shown in Figure 16. The model's predictions agree well with the experimental results (R2 > 0.90); this


\* Values shown represent the average of three runs.

Table 8. Experimental parameters used in the numerical model.

Figure 16. The relationship between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p.

result is consistent with the findings of Meyer [32], Foster [33], and Liebenow et al. [34], whose results were incorporated into the USDA WEPP model [35]. When the slope gradient is 15 or greater, r = 1, which is consistent with the results of Gao et al. [29] and corroborates the conclusions drawn by Sharma et al. [36, 37] and Jayawardena and Bhuiyan [38]. The relationships under these conditions between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p, at slope gradients of 15 and 25, are shown in Figure 17. The values of the bare-soil detachability parameter, a, were calculated using Eq. (25) and are shown in Table 9. The values of the bare-soil detachability parameter, a, shown in Table 9, were quite similar, which may be attributed to the fact that the soil detachability is constant when the ponding depth is below a critical or breakpoint depth [30, 31, 38, 39]. Parameter b was estimated using

the best fit to the experimental data listed in Table 9. The predictions made using Eq. (26) are

Figure 17. The relationship between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p, at slope

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

Initial soil moisture content (%) 5 2.18 0.15 2 0.4 0.094

Slope gradient () 5 1.85 0.15 2 0.35 0.082

) p (cm min<sup>1</sup>

) 60 1.16 0.1 2 0.2 0.116

15 1.9 0.15 2 0.4 0.083 20 1.68 0.15 2 0.4 0.075

96 1.86 0.16 2 0.2 0.082 129 3.93 0.215 2 0.2 0.085

15 10.11 0.15 1 0.35 0.072 25 11.93 0.15 1 0.35 0.082

) r a (g cm<sup>3</sup>

http://dx.doi.org/10.5772/intechopen.76280

155

) b\*

The simulated data agreed well with the experimental results for all three treatments except for the experiment where the rainfall intensity was 129 mm/h, suggesting that the use of the raindrop-induced water transfer rate, er, in place of the exchange rate, km, is reasonable for conditions involving relatively natural rainfall. The model, however, did not appear to accurately predict the solute concentrations in the runoff observed under conditions that give rise

compared to our experimental results in Figure 18.

Table 9. Data used in calculating the bare-soil detachability, a.

Calibrated to best fit the runoff solute data.

Treatment <sup>e</sup> (<sup>10</sup><sup>3</sup> g cm<sup>2</sup>

to severe soil erosion.

gradients of 15 and 25.

Rainfall intensity (mm h<sup>1</sup>

\*

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control http://dx.doi.org/10.5772/intechopen.76280 155

Figure 17. The relationship between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p, at slope gradients of 15 and 25.


Calibrated to best fit the runoff solute data.

result is consistent with the findings of Meyer [32], Foster [33], and Liebenow et al. [34], whose results were incorporated into the USDA WEPP model [35]. When the slope gradient is 15 or greater, r = 1, which is consistent with the results of Gao et al. [29] and corroborates the conclusions drawn by Sharma et al. [36, 37] and Jayawardena and Bhuiyan [38]. The relationships under these conditions between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p, at slope gradients of 15 and 25, are shown in Figure 17. The values of the bare-soil detachability parameter, a, were calculated using Eq. (25) and are shown in Table 9. The values of the bare-soil detachability parameter, a, shown in Table 9, were quite similar, which may be attributed to the fact that the soil detachability is constant when the ponding depth is below a critical or breakpoint depth [30, 31, 38, 39]. Parameter b was estimated using

Treatment t<sup>p</sup> (min) Ho (cm) θ (%) Initial soil moisture content 5% 2.83 0.2 5

Rainfall intensity 60 mm h<sup>1</sup> 4.00 0.18 10

Slope gradient 5 3.9 0.16 10

\*

Values shown represent the average of three runs.

154 Hydrology of Artificial and Controlled Experiments

Table 8. Experimental parameters used in the numerical model.

15% 1.83 0.23 15 20% 1.21 0.24 20

96 mm h<sup>1</sup> 2.40 0.23 10 129 mm h<sup>1</sup> 1.70 0.27 10

15 3.3 0.18 10 25 2.3 0.3 10

Figure 16. The relationship between the rainfall-induced soil detachment per unit area, e, and rainfall intensity, p.

Table 9. Data used in calculating the bare-soil detachability, a.

the best fit to the experimental data listed in Table 9. The predictions made using Eq. (26) are compared to our experimental results in Figure 18.

The simulated data agreed well with the experimental results for all three treatments except for the experiment where the rainfall intensity was 129 mm/h, suggesting that the use of the raindrop-induced water transfer rate, er, in place of the exchange rate, km, is reasonable for conditions involving relatively natural rainfall. The model, however, did not appear to accurately predict the solute concentrations in the runoff observed under conditions that give rise to severe soil erosion.

The agreement of the simulated results with the measured data was quantified by calculating

P<sup>N</sup>

where N is the total number of data points, p<sup>i</sup> is a given simulated data point, and o<sup>i</sup> is the corresponding experimental data point. The RMSEs are shown in Table 10. Table 10 and Figure 18 illustrate that the model has correctly captured the temporal behavior of the solute

Figure 18A shows that the measured and simulated solute concentrations for different initial soil moisture contents changed with time. The results indicated that the refined model [25] could predict the movement of solutes in the overland flow under different initial soil moisture contents. Also, the higher initial soil moisture contents were associated with higher solute

Initial soil moisture content (%) 5 0.227 0.90

Rainfall intensity (mm/h) 60 0.295 0.80

Slope gradient (�) 5 0.081 0.94

Table 10. The root means square errors (RMSEs) and R<sup>2</sup> between the measured and simulated data.

Figure 19. Variation of the sediment concentrations over time at different rainfall intensities.

s

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

<sup>i</sup>¼<sup>1</sup> pi � oi � �<sup>2</sup>

(27)

157

) R2

http://dx.doi.org/10.5772/intechopen.76280

15 0.245 0.93 20 0.308 0.91

96 0.229 0.88 129 0.508 0.46

15 0.127 0.88 25 0.336 0.86

N

the root mean square error (RMSE) [31]. RMSE can be expressed as:

concentration in the runoff under all conditions investigated.

Treatment RMSE (mg l�<sup>1</sup>

RMSE ¼

Figure 18. Simulated and measured runoff concentrations of potassium under different initial soil moisture contents (A), rainfall intensities (B), and slope gradients (C).

Our results support the conclusion drawn by Walter et al. [40], who argued that the depth of the exchange layer decreases as the rate of infiltration increases. The initial soil moisture content, rainfall intensity, and slope gradient influence the solute concentration of the runoff solution by virtue of their effects on the depth of the exchange layer, the infiltration rate, and the length of time between the initiation of rainfall and the formation of the runoff.

The agreement of the simulated results with the measured data was quantified by calculating the root mean square error (RMSE) [31]. RMSE can be expressed as:

$$RMSE = \sqrt{\frac{\sum\_{i=1}^{N} \left(p\_i - o\_i\right)^2}{N}} \tag{27}$$

where N is the total number of data points, p<sup>i</sup> is a given simulated data point, and o<sup>i</sup> is the corresponding experimental data point. The RMSEs are shown in Table 10. Table 10 and Figure 18 illustrate that the model has correctly captured the temporal behavior of the solute concentration in the runoff under all conditions investigated.

Figure 18A shows that the measured and simulated solute concentrations for different initial soil moisture contents changed with time. The results indicated that the refined model [25] could predict the movement of solutes in the overland flow under different initial soil moisture contents. Also, the higher initial soil moisture contents were associated with higher solute


Table 10. The root means square errors (RMSEs) and R<sup>2</sup> between the measured and simulated data.

Figure 19. Variation of the sediment concentrations over time at different rainfall intensities.

Our results support the conclusion drawn by Walter et al. [40], who argued that the depth of the exchange layer decreases as the rate of infiltration increases. The initial soil moisture content, rainfall intensity, and slope gradient influence the solute concentration of the runoff solution by virtue of their effects on the depth of the exchange layer, the infiltration rate, and

Figure 18. Simulated and measured runoff concentrations of potassium under different initial soil moisture contents

the length of time between the initiation of rainfall and the formation of the runoff.

(A), rainfall intensities (B), and slope gradients (C).

156 Hydrology of Artificial and Controlled Experiments

concentrations per unit time. Figure 18A and Table 10 indicate that the differences between the measured and simulated solute concentrations under an initial soil moisture content of 20% were more distinct than those under the other two initial soil moisture contents, which implied

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

159

Figure 21. Graph of potassium concentration calculated by Eq. (26) against potassium concentration measured (under

different slope gradient conditions).

Figure 20. Graph of potassium concentration calculated by Eq. (26) against potassium concentration measured (under different initial soil moisture content, rainfall intensity conditions).

concentrations per unit time. Figure 18A and Table 10 indicate that the differences between the measured and simulated solute concentrations under an initial soil moisture content of 20% were more distinct than those under the other two initial soil moisture contents, which implied

Figure 21. Graph of potassium concentration calculated by Eq. (26) against potassium concentration measured (under different slope gradient conditions).

Figure 20. Graph of potassium concentration calculated by Eq. (26) against potassium concentration measured (under

different initial soil moisture content, rainfall intensity conditions).

158 Hydrology of Artificial and Controlled Experiments

that the model did not accurately predict the solute concentrations of the runoff in conditions of severe soil erosion.

content, the improved Côté-Konrad model is superior to Côté-Konrad model, Lu-Ren model, and the improved Lu-Ren model. For the soils with low sand content and silt content, the Lu-Ren model is obviously better than the other three models. The relationship between the parameters of the model, particle composition, and organic matter content can be predicted by two improved models. These models can describe the relationship between the soil's basic physical parameters and thermal conductivity in detail. Thus, soil thermal conductivity can be predicted more accurately by choosing the appro-

Experimental Study on the Mechanisms of Soil Water-Solute-Heat Transport and Nutrient Loss Control

http://dx.doi.org/10.5772/intechopen.76280

161

3. The refined power functions of a model of solute transport were illustrated and tested using simple experiments. The model fit the experimental data very well. Our results also indicated that the constant parameter, r, was equal to 1 when the slope gradient was 15 or larger and equal to 2 when the slope gradient was less than 15. The soil detachability was confirmed to be independent of the rain intensity and was a constant in all treatments. The model, however, could not accurately predict the solute concentrations in the runoff under conditions of severe soil erosion. The initial soil moisture content, rainfall intensity, and slope gradient influenced the solute concentration in the runoff, depth of the exchange layer, infiltration rate, and length of time between the initiation of rainfall and the genera-

This study was financially supported by the National Natural Science Foundation of China (grant nos. 51239009, 41371239), Science and Technology Planning Project of Shaanxi Province (2013kjxx-38), and Doctoral fund of Xi'an University of Technology (106-211301). We also

State Key Laboratory of Eco-Hydraulic Engineering in Arid Area, Xi'an University of

[1] Darcy H. Less Fontaines Publiques de la Ville de Dijon. Paris, France: Dalmont; 1956

[2] Green WH, Ampt CA. Studies on soil physics. I. Flow of air and water through soils.

priate improved model based on the different soil texture.

tion of the runoff.

Acknowledgements

Author details

Technology, Xi'an, China

References

thank Xiaopeng Chen for his helpful comments.

Quanjiu Wang, Beibei Zhou\*, Lijun Su and Yuyang Shan

Journal of Agricultural Science. 1911;4:1-24

\*Address all correspondence to: happyangle222@aliyun.com

Comparisons between the simulated and the experimental solute concentrations for the different rainfall intensities over 60 min are shown in Figure 18B. At a rainfall intensity of 129 mm/h, the solute concentration of the runoff increased substantially between 37 and 49 min after the initial generation of the runoff (Figure 18B). The mass of sediment in the runoff between 37 and 43–49 min showed a corresponding spike (Figure 19), which indicated that solute loss is closely related to sediment loss [41–45]. These results indicated that significant erosion of the surface soil occurred at the bottom of the slope during the experiments. Deeper soil layers were exposed to water in which the solute concentrations were higher than in those washed away. Consequently, the solute concentration of the runoff increased as these solutes were transferred from the soil under the influence of the runoff and the splashing caused by raindrops. Soil erosion thus promoted increased solute concentrations in the runoff.

Figure 18C shows that the measured and simulated solute concentrations for different slope gradients also changed with time. The simulated data were highly correlated with the measured data for solute concentration in the runoff. This degree of correlation demonstrated that the model captured the temporal behavior of the solute transport in the runoff. Increasing the gradient of the slope increased the erosion capacity of rain drops and water flow. Increasing the slope gradient also led to increases in the RMSE (Table 10) and R2 . Figures 20 and 21 show the relationships between potassium concentrations observed in the runoff and predicted using Eq. (26). The graph indicates the model accurately predict the solute transport in the runoff with the solute concentration being at a much lower level.
