Soil parameterization
To develop soil parameterization file, first the depth of the different soil horizons and texture (combination of soil, silt and clay) in percentage is required. This information is general available in soil dataset. The texture information is useful to describe the characteristics of soil water retention curve. The soil texture information (in percentage) is provided as input data to the software component ROsetta inside 'HYDRUS 1D' to understand pedotransfer functions of soils in three different hypothetical pressure scenarios. They are: 0 mbar, 60 mbar and 15,000 mbar. The soil texture is classified according to the specific pore volumes of soil (Scheffer and Schachtschabel 1984). The soil texture information provides additional
For example, a soil type A has following depth and texture combination:
| Horizon | Depth (cm) | Texture (Sand,Silt,Clay %) |
|---|---|---|
| A Horizon | 15 | texture (47, 39, 14 ) |
| B Horizon | 15 | texture (40, 42, 19 ) |
| C Horizon | 25 | texture (33, 40, 27 ) |
| D Horizon | 35 | texture (49, 33, 18 ) |
The texture information needs to be supplied to HYDRUS 1D software to get the water holding capacity of the specific texture combination:
- The HYDRUS-1D can be download from the following link: http://www.pc-progress.com/en/Default.aspx?hydrus-1d
- Click Hydrus-1D Downloads-->
- Download the latest version (In this tutorial, version Hydrus1D_4.15.0110 is downloaded)
- Open the software:
Click File--> New-->>
Give the name for the new project. For example: Soil type-A
Click OK
A new window with new project name will appear:
On pre-processing windows select the “Water Flow – Soil Hydraulic Parameters”
The new window will appear:
On this window, select the “Neural Network Prediction”.
A new window will appear:
Select model of “% Sand, Silt and Clay (SSC)”.
Input the percentage of SSC in the window. In this example, the SSC of A-horizon from the table above is provided.
Click ‘Predict’ and the Output will be resulted.
Click ‘Accept’ and a new window will appear.
Note down the value in the volume Qr which is 'Residual Water Content' equivalent to wilting point. The value represents the water content in 15,000 mbar pressure.
Click ‘Next’ a new windows appear.
In the new window, make sure that the Upper and Lower Boundary Condition is set to ‘Constant Pressure Head’ and Initial Condition is ‘In Pressure Heads’.
Then click Next
Click next and a new window will appear:
Click Next
Soil profile summary table will appear
Click Next
A new window with the following message will appear: 'Do you want to run HYDRUS-1D application?' Click OK
Several outputs will appear in Post-processing window
Press Enter to continue.
Now your data is saved and you will be taken to initial windows.
Double click the Soil Hydraulic Properties and a new graph will appear
Keep the horizontal variable 'pressure head' and verticle 'water content'.
Right click on the graph and click 'Edit Chart Data'.
A new window 'data grid editor' will appear.
The C1 is pressure (mbar) and C2 is water holding capacity (%)
Take a note of the pressure of 0, 59.9 and related water holding capacity. Since, the water content equivalent to 15,000 mbar is not available, the vaule is taken from Qr (Residual Water Content) which has been described earlier. You will get the values 0.39, 0.35 and 0.05 respectively.
By substracting water content value from 60mbar to 0 mbar (i.e. 0.39-0.35 = 0.04) will give the water available in large pore storage. The soil can hold water available in LPS against the gravity. By substracting water content value from 15000 mbar to 60 mbar (i.e. 0.35-0.05 = 0.30) will provide water available in middle pore storage. This state of water availability is also known in field capacity.
The total volume of MPS and LPS can be derived by multiplying MPS and LPS value with respective soil depth (mm) of A horizon. i.e. 150 * 0.30 = 45 (MPS) and 150 * 0.04 = 6 (LPS) as shown in Figure below for A horizon. By considering all the available horizon as shown in the first place, the following information can be derived.
The following MPS value (268.5) has to be equally divided into each 1 dm zone of respective soil horizone in the parameter file as shown in the figure below.
For example in the table below, the value for fc_1 (field capacity:MPS) for first 1 dm is 30 because MPS in A horizon (15 cm) is 30. The next fc_2 (2 dm) layer is located in between A and B horizon. Therefore, the average has to be taken from the MPS value of A and B horizon (i.e. 30+31 = 30.5). Similarly, the MPS value for each dm is provided until 22 dm.
The value for LPS does not need to be divided into the different soil depth.
Users has to provide the numerical id for each soil type. For example, the soil type A from the example is provided the id 11 in the parameter file.
