Achievements

STUDY ON THE NUMERICAL SIMULATION OF MOISTURE BEHAVIOR IN FIXED DUNE UNDER WATER EXTRACTION BY PLANT ROOTS

Updated :10,23,2012

Qian Ju1Long Xunjian1Ma Jinzhu1Zhang Huichang1*Yi Lixing2

1.Key Laboratory of Western China’s Environmental Systems, MOE, Lanzhou University, Lanzhou 730000, China;

2.College of Environmental Science and Engineer, Nankai University, Tianjin 300071,China;*Author for correspondence

 

Abstract: Water uptake by roots in fixed dune was discussed in this paper. The field experiments had been daily carried out on time within the 0cm175cm depth of the fixed dune planted Nitraria tangutorum which is located in the north of Dengmaying basin at the southwestern edge of Tengger desert in Shiyanghe river basin in the northwestern China’s Gansu province to measure moisture content and moisture potential with sand layers by Neutron moisture meter and Tensionmeter. Based on the measured data in the field, the parameters of dune moisture transport, such as unsaturated conductivity, moisture diffusivity etc., were made using mathematical statistics method. Then water uptake model by Nitraria tangutorum roots which was related to the rate of evaporation measured by evaporator and dune moisture contents was established at the macroscopic angle. In addition, the numerical simulations of moisture behavior in fixed dune during the growth of Nitraria tangutorum and rainless period in June and August were conducted. The relative errors between simulated data and measured ones were -21.97%5.67%, mainly -0.016%3.56%, which was indicated in the numerical simulation of dune moisture at the depth of 25cm, 55cm and 100cm on June 5th, 10th and 12th. Meantime, the predicted curves of moisture distribution and the measured ones in the dune profile in autumn were compared. The results from the verification and analysis of dune moisture simulated data showed that the precision of root uptake rate was higher, and it was able to indicate the regularity of water uptake during the rainless period, the predicted moisture distribution curve in the dune profile during Nitraria tangutorumgrowth season by the use of mathematical model was basically consistent with the measured one.

Key words: fixed dune; water uptake by plant roots; numerical simulation; moisture prediction

 



Introduction

Limited precipitation is the main source for savageness vegetation existent of no irrigation in desert environment. The plants which grow in desert absorb precipitation mainly depended on dune moisture. Root system is the plants apparatus for absorbing, and it happens every minute. Even though lamina stoma closed, and without transgress of rising moisture during night, the root system are still absorbing little moisture by one certain speed, in order to replenish overfull loss of plants in the daytime, and modulate moisture balance of plants (Kang et al., 1994). It is important part of moisture recycle in groundwater-dune-plant-atmosphere continuum (GSPAC), and is also one of the main contents of plant and moisture connection study (Shi et al., 2003; Zhao et al., 2002). The GSPAC environment factors including dune temperature, moisture content, water solution concentration, aerate status and weather condition would affect the water uptake by plant roots that grow in desert (Zhao et al., 2002; Ji et al., 2006). When the dune containing water in section is different, it affects the distributing of boot system absorbing in section. For example, the main bibulous region appears in shallow dune after precipitation, while sand bed becoming dryer, water uptake zone moves to down direction (Kang et al., 1994).

1 General Situation of Study Area

The study plot lies in the north of Dengmaying Lake103°15′E38°7′Nwhich is the largest lake bottomland at the southwestern edge of Tengeer desert in Shiyanghe river basin of Gansu province in northwestern China, and its weather belongs to typical aridity hungriness climate. Mean annual precipitation is 188.0mm, and fasten from June to August, rainfall occupy 82.4% of the whole year. Mean annual evaporation of water surface is 1880.9mm. Normal annual temperature is 7.8, and frost-free period is 150d. Ordinary year prevail northwester, mean annual sandstorm day is 120d, and mainly from March to May. Endow with groundwater belongs to unconfined water in desert; the lithology of containing water layer is Upper Permian graining sand and silt sand. Groundwater level is 270cm, and submit period change, the changing range is small, average value is 36 cm (Qian et al., 1999, Yi et al., 1997).

Setting observation point on the elliptic dune (research dune) covered with Nitraria tangutorum bush of study plot that is 5m high and wide, 20m long. Using the ZG-1 Neutron moisture meter and WM-1 mercury Tensionmeter regular timing observe moisture content and moist potential at the depth of 0175cm day by day at the same time. And there is atmosphere observation station in open place that is apart from the dune100m away, which observes precipitation, water surface evaporation and temperature.

Research dune is mainly made up of middle sand and graining sand, which occupied about 90.74%, and the particle diameter smaller than 0.05mm is just 0.03%. Main moist physics constants are field moisture capacity4.51%, maximum hygroscopicity0.17%. Capillarity water rise height is 57cm, total porosity is 41.93%. Cover degree of Nitraria tangutorum bush that develop well in nature on dune is upper, and boot system is up growth.

2 General Resume of the Water  Uptake Models by Plant Roots

The intensity of root system water absorption is showed by root system water absorptivity, that means per volume soil of root system absorbed moist volume in unit period (Shao, 1996; Lei et al., 1988). The research of root system absorption rule and water absorptivity shows, the studies are most based on root system absorption models. And totally including three types: (1) experience model, just as simple root absorption models put forward by Gardner (Gardner,1960), wheatcorn root system absorption models established by Yao Jianwen (Yao,1989). (2) semi-theory and semi-experience model, Molz-RemsonNimah-Hanks models etc.(Molz F. J. et al., 1970; Nimah M. N. et al., 1973). (3) Theory model, for example, Hillel lodged the model on the basis of simulate water movement by electricity (Kang et al., 1994; Hillel D., 1976; Zhao et al.,1999).And many root system models were established successfully. The representative modelscalculated by root system water absorptioninclude: (1) Assuming root system absorptivity pro rata with soil-water potential or soil containing moisture quantity. (2) Using transpiration consume water capacity in root system to estimating root system absorptivity (Kang et al., 1994). The research of root system that begin with the end of 50 decades and early 60decades in 20 century suffer from develop modeling, rework and perfect, and have gained quite systemic study production as yet, but researching of configuration, absorption mechanism of root system, specialty application and process visual and etc. are still in investigating (Zhao et al., 1999).

The plant root system absorption models in existence mostly include parameters difficult to assured outside such as density, resistance or perviousness and etc., and some difficulties in actually application. But whether found root system absorption model, or analysis distributing rule of root system water absorptivity, both need to take the root system water absorptivity dynamic data. So, we can consider the root system as a whole one, confirm root section moist movement equation from microscope angle, and estimate root system water absorptivity by observed data in fields (Kang et al., 1994).

3 Numerical Simulation of Moisture Behavior in Dune under Water Extraction by Plant Roots

3.1Confirming of dune moisture dynamic parameters

Determination of unsaturated coefficient of conductivity, dune moisture diffusivity in field is important content for dune moisture dynamics research, and it is also crucial precondition of dune moisture dynamic simulation in root system absorption condition.

3.1.1 Dune moisture characteristics curve

It means the relation curve that dune moisture substratum potential jm or suction s change with volume moisture contentq, and shows the relation between dune moisture energy and quality, also important curve that response the basic characteristics of dune moisture (Lei et al., 1988).

This paper chose observed jm-q data in filed, using elements of least-squares method and experience function s=abq, and by regression analysis gets the expression of dune moisture characteristic curve and it can be written as:

 s=61.7946´0.8314100q      

r=0.9228                          1

3.1.2 Unsaturated transmissibility

Basis on Darcy’s law and moisture balance theory corresponding unsaturated transmissibility could be calculated due to q. Through observed jm and q data, regression by least-squares method obtains as following (Ma et al.,1998; Zhang et al.,1998):

Kq=3.5157´10-8´1.6375100q   

r=0.9259                     2

Where Kqis the unsaturated transmissibility (cm/s).

3.1.3 Diffusivity of dune moisture

Calculating from expressions of diffusivity Dqand sq,i.e.:

Dq=Kq/Cq=4.0904´10-4´1.3614100q3

Where Dqis diffusivity of dune moisturecm2/sCqis specific capacitycm-1,that is named dq/djm.

3.2 Establishment of water uptake model  by plant roots

Single root experience model is one of the sublimate model, it is propitious to analyze root system absorption mechanism, relation of root-moist potential and soil-moisture potential, but inapplicable to dynamic dune base moist modeling (Kang et al.,1994; Zhao et al.,1999). Therefore, it should resolve actual problems of dune moisture movement based on macroscopic level.

Study of root system absorption on macroscopic points is adding root absorption item S on the right of dune moisture movement equation, when vertical is one-dimensional flow, the mathematics equation is given as (Kang et al.,1994; Liu et al., 2004):



           4



Where t is time, z is depth of sand bed, and downwards is positive; DT is moist diffusivity effected by temperature, is temperature, Sz, troot absorption item, named root water absorption.

Other denotation meaning is the same as frontal.

When dune moisture content change little, temperature dispatch direct evocable moisture movement is small,  can dismission, and also ignore dune moisture movement parameter infected by temperature challenge, and formula4is transformed into:



                             5

                            6



Formula6is nonlinear, non-homogeneous partial differential equation, adopting finite-difference equation evaluates the closest result.

Change formula6into difference equation by implicit difference equation, that is using the difference coefficient instead partial differential equation of the equation and then gets the result:



   7



Where i is special step serial number, =0, 1, …, n; j is time step serial number, j=01m

Each root system water absorption  of dune section in facultative period of time can calculate by substituting parameter KqDqto formula7. The calculation depth = 175cm, calculation period is during Nitraria tangutorum growth period (from May to September), special step is:

 cmperiod step Dt=1d.

Adopting nonlinear regress analysis method gets research region root system water adsorption changing with moisture content q, and the experience equation is:

(zt)= E(t) AB100qzt    

r = 0.9492                     8

Where E(t) is water surface evaporation in t period, A=1.9728´10-4 B=1.1858qztis average volume water content of dune depth z in period t%.

3.3 Numerical simulation of moisture behavior in dune

3.3.1 Elementary condition

Distributing of dune elementary moisture content has already knew, thus:

                   9

3.3.2 Boundary condition

The precipitation of proving ground is relatively little, so dune surface didn’t form water and run-flow. The simulation period during continuum evaporation after secondary precipitation was selected. The upper boundary of evaporation control condition is (lei et al.,1988):

        10

Where E(t) is dune surface evaporation intensity in t period, it is calculated by regression Stat. relation of formula11 (Yi et al.,1997; Ma et al.,1998; Zhang et al.,1998).

   11

Where qis average volume water content of 25cm depth of dune surface layer%.

The lower boundary condition takes the first kind edition by the button measure point of neutron instrumentz=175cm:

                 12

3.3.3 Numerical simulation

The modeling calculation region is dune aeration zone, the depth is 175cm. Express by equationq, mathematics positive value questions considering vertical one-dimensional moisture movement of plant root system absorption is:



           13



Where the symbols are the same as forward.

Formula13is two- dimensional nonlinear partial differential equation, evaluation by numerical value method. Researching adopt unit horizon area and height for calculation sand column. By two group perpendicular with one another liner divided calculation region into rectangular gridding, deviding depth into nn=10 unit from z direction, and node number is j=012,…,mspecial step is Dzthe same as front, and provide time into m period, node number is j=012,…,mperiod step is Dt=1d. each node including inside node and outside nodeadopt implicit iterative method separately, and get the quantity of water equilibrium equation:



           14

Where b0c0d0dn-1aibicidare calculated by formula15(Kang et al., 1994).

              15



Where r1=Dt/Dz2r=Dt/2Dz), is average evaporation intensity of period,  are sand pillar zi-1/2zi+1/2 section moisture diffuse ratio in j+1 period respectively,  are sand pillar zi-1/2zi+1/2 section unsaturated conductivity in j+1 period respectively,  are Dand K, qji is point i in period j volume water content,  is root system water adsorption of point iin period j+1/2.

Formula 14is nonlinear three diagonal equation, evaluation by pursuit method after iterative linear, the evaluation progress includes elimination and back substitution. While numerical value calculating, Kq,Dqare presented by working experience equation, and other value of parameters are calculated by geometry average method.

4 Results Verification and Analysis of Numerical Simulation of Moisture Behavior in Dune

Modeling the dune moisture dynamic state during dry period of Nitraria tangutorum growth time in June and August by mathematic model, and contrasted with observed data. Total modeling time is 27d, and chosen time of June 5, June 10, and June12 simulation data of dune water content compare with observed datatable 1, relative error is between -21.97%5.67%, and most between -0.016%3.56%. And also present dune water content distributing simulation section of June 22, August 24 comparison progress with observed sectionfig.1



June 22                         August 24

Fig.1 Comparison of dune moisture simulation valuesawith observed valuesbduring Nitraria tangutorum growth period

 

Table 1 Moisture simulation and observed values at different depths in dune in Junecm%

Date

 

25

 

 

55

 

 

100

 

day/month

Observed value

Simulationvalue

Relative error

Observed value

Simulationvalue

Relative error

Observed value

Simulationvalue

Relative error

5/6

2.385

2.470

3.56

4.982

4.720

-5.26

6.089

6.434

5.67

10/6

1.925

1.502

-21.97

4.862

4.975

2.32

5.647

5.756

1.93

12/6

1.856

1.603

-13.63

5.109

5.221

2.19

6.433

6.432

-0.016



In order to verify the reliability of model, moisture behavior in dune was predicted during Nitraria tangutorum growth later period in autumn.



                               

Fig.2 Comparison of dune moisture forecast valuesawith observed valuesbon Sep. 2

 



5 Conclusions and Discussion

The root system absorption model in this paper can correctly reflect root system absorption rule in dry spell, dynamic numerical value modeling result is testing well with observed value, and this model also could forecast dune section water content distributing during plant growing perioddry period. However, the key of exactly build up of full purport root system model is accurate determinate dune water content and moisture potential, improve dune moisture movement parameter precision, check out distributing characteristic of root system structure, extend depth, root system density and etc, perfect dune moisture movement elements researching. Dynamic numerical value modeling and forecast dependability depend on mathematics assure problem correction, modeling period special step and period step confirm reasonably. In a word, these problems still need more discussion, resolving and perfection.

Acknowledgements

This study was funded by the National Natural Science Foundation of China (No. 49272143) and the National Key Development and Programming Project for Basic Research of China (No. G2000048701).



 

Reference

[1]      Kang Shaozhong, Liu Xiaoming, Xiong Yunzhang. Theory and its application of moisture

[2]      movement in soil-plant-atmosphere continuum[M]. Beijing: Water Conservancy and Electric Power Press1994. 5-1451-5685-115 (in Chinese).

[3]      Shi XiaohongLiu TingxiLi Changyou. Analysis of Soil Water Characteristics in Zone of Main

[4]      Plant Roots in Low Sand Dune of Horqin Sandy Land[J]. Journal of Inner Mongolia Agricultural

[5]      University2003243):23-29 (in Chinese).

[6]      Zhao Wenzhi. Impact of Plantation on Spatial Heterogeneity of Soil Moisture in Horqin Sandy

[7]      Land[J]. Acta Pedologica Sinica2002391):113-119 (in Chinese). 

[8]      Ji XibinKang ErsiChen Renshengetc..Research Advances about Water-uptake Models by Plant Roots[J]. Acta Bot Boreal Occident Sin2006265):1079-1086 (in Chinese).

[9]      Qian Ju, Ma Jinzhu, Zhang Huichang, et al. Study on the Characteristics of Moisture Behavior and Moisture Potential Changes in Fixed Dunes on the Southwestern Edge of Tengger Desert[J]. Journal of LanzhouUniversity (Natural Science Edition). 1999351):218-224 (in Chinese).

[10]  Yi Lixin, Zhang Huichang. experimental study on moisture movement in unsaturated zone of arid’s dune[J]. Arid Zone Research1997142):14-18 (in Chinese).

[11]  Shao Aijun. Numerical simulation of soil moisture movement as a case study of water uptake item by crops roots[J]. Journal of Hydrogeology and Engineering Geology1996232):5-9 (in Chinese).

[12]  Lei ZhidongYang ShixiuXie Senchuan. Soil moisture hydrodynamics[M]. Beijing: Tsinghua University Press1988. 18-24205-214 (in Chinese). 

[13]  Gardner W. R.. Dynamic aspects of water availability to plants[J]. Soil Science1960892):63- 73.

[14]  Yao Jianwen. A numerical model of predicting soil moisture content under the condition of crop growing[J]. Journal of Hydraulic Engineering1989,(9):32- 38.

[15]  Molz F. J.Remson I.. Extraction term models of soil moisture use by transpiring plants[J]. Water

[16]  Resources Research1970,(6):1346-1356.

[17]  Nimah M. N.Hanks R. J.. Model for estimating soil water, plant and atmospheric interrelations.1.

[18]  Description and sensitivity[J]. Soil Science Society of American Journal1973,(37):522-527.

[19]  Hillel D.. A macroscopic scale model of water uptake by a non-uniform root system and salt

[20]  movement in the soil profile[J]. Soil Science1976,(121):242-255.

[21]  Zhao ChengyiHuang JunmeiWang Yuchaoet al.. Study on the specificity of water uptake by plant roots[J]. Arid Land Geography1999222):88-96 (in Chinese).

[22]  Ma Jinzhu, Zhang Huichang,Yi Lixin,eta. A coupled model of water-vapor-heat transport and numerical simulation in unsaturated zone of tengeer desert[J]. Journal of Desert Research1998184):340-345 (in Chinese).

[23]  Zhang Huichang, Shi Peize, Wu Xiuzhen. Study on the Mechanism of dune moisture non-equilibrium movement and its exploitation and utilization[M]. Beijing: Earthquake Press1998. 102-108 (in Chinese).

[24]  Liu ZengjinChai HongminXU Jianxin. Soil Water Numerical Calculation of Winter Wheat[J].

[25]  Journal of Irrigation and Drainage2004232):73-76 (in Chinese).