Green-Ampt Infiltration Parameters: Difference between revisions

Content deleted Content added
 
(8 intermediate revisions by 3 users not shown)
Line 2:
TUFLOW provides several methods for modelling infiltration from the 2D surface into the sub-surface, including Green-Ampt, Horton, and Initial Loss/Continuing Loss. These methods are used to simulate hydrological losses, particularly when rainfall is applied directly to the 2D surface and runoff is generated.
 
The choice of infiltration method and its parameters areis an important calibration factorsfactor and should be adjusted to match observed flow data. This is especially relevant for whole of catchment modelling, where infiltration is the main way hydrological losses are represented. This page describes the Green-Ampt infiltration parameters and their sensitivity.
 
== Green-Ampt Infiltration ==
Line 82:
 
== Green-Ampt Infiltration: User Parameters ==
Where the inbuilt USDA soil types are not used, the user can specify their own values for the Suction, Hydraulic Conductivity, Porosity and Initial Soil Moisture. What follows is a description of each parameter and the sensitivity to a low, medium and high value based on the USDA soil type summary values.<br>
 
=== Capillary Suction Head ===
Line 95:
<br>
 
[[File:Fig1Suction Head.jpegjpg|border|600px760x760px|Figure 4: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the Suction head parameter in the Green-Ampt infiltration model.|border]]<br>
 
'''Figure 4: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the Suction head parameter in the Green-Ampt infiltration model.'''<br>
<br>
It can also be seen that the higher the suction head value that the longer it takes the hydrograph to start rising, with the high suction head scenario less responsive to the rainfall.
<br>
 
=== Saturated Hydraulic Conductivity ===
Line 110 ⟶ 109:
<br>
 
[[File:Fig6Hydraulic hydroconductConductivity.pngjpg|border|600px760x760px|Figure 5: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the Saturated Hydraulic Conductivity parameter in the Green-Ampt infiltration model when using a clay soil type.|border]]<br>
 
'''Figure 5: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the Saturated Hydraulic Conductivity parameter in the Green-Ampt infiltration model when using a clay soil type.'''<br>
Line 117 ⟶ 116:
[[File:Fig7 hydroconduct.png|600px|Figure 6: Sensitivity of cumulative infiltration in the Plynlimon Gwy catchment to the Saturated Hydraulic Conductivity parameter in the Green-Ampt infiltration model.|border]]<br>
 
'''Figure 6: Sensitivity of cumulative infiltration in the Plynlimon Gwy catchment to the Saturated Hydraulic Conductivity parameter in the Green-Ampt infiltration model.'''<br>
 
=== Porosity ===
Line 123 ⟶ 122:
<br>
 
[[File:Fig7 porosity sensPorosity.pngjpg|border|600px760x760px|Figure 7: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the porosity parameter in the Green-Ampt infiltration model.|border]]<br>
 
'''Figure 7: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the porosity parameter in the Green-Ampt infiltration model.'''<br>
Line 131 ⟶ 130:
 
'''Figure 8: Sensitivity of cumulative infiltration in the Plynlimon Gwy catchment to the porosity parameter in the Green-Ampt infiltration model.'''<br>
<br>
 
=== Initial Moisture ===
The initial moisture value represents the fraction of the soil that is initially wet. As both initial moisture and porosity are expressed as fractions, the soil capacity is defined as the difference between them both. As such, the initial moisture should not exceed the porosity otherwise soil capacity will be set to zero with no infiltration occurring for that soil type. A [[TUFLOW Message 2508 | 2508 WARNING]] is issued if this is the case.<br>
 
In pre-2023 releases of TUFLOW, a single variable storage capacity was calculated by subtracting the initial moisture fraction from the porosity, in order to reduce memory requirements. However, in TUFLOW releases 2023 and onwards, the soil porosity and initial moisture must be stored separately to allow the soil to drain correctly when using the interflow functionality. This updated approach requires that a soil thickness be specified to calculate the soil depth. If a soil thickness is not specified when using the updated approach, an infinite soil depth is assumed for each layer and therefore different initial moisture fractions no longer have an effect on modelled results. These two approaches can generate different results when using the Green-Ampt method.
As you increase initial moisture at the beginning of your simulation, you experience less infiltration (as you are closer to the soil capacity), therefore have more run-off and a quicker response. Figure 9 shows the degree of change to cumulative infiltration with varying initial moisture and the effect on the catchment can be seen in Figure 10. As the event progresses, soils become more saturated and the influence of the initial moisture parameter becomes less significant.
 
The three initial moisture sensitivity tests have been undertaken with the Green-Ampt method using both a pre-2023 release of TUFLOW and a post-2023 release of TUFLOW. Figure 9 shows how variations in the initial moisture affect the simulated cumulative infiltration, whereas Figures 10 and 11 show the effects of varying the initial moisture on flows at the catchment outlet when using the pre-2023 and post-2023 releases of TUFLOW. As the initial moisture is increased at the beginning of your simulation, there is less infiltration (as you are closer to soil capacity) and more runoff, causing the catchment outflows to exhibit a faster response to rainfall upstream. As the event progresses, soils become more saturated and the influence of the initial moisture parameter becomes less significant. In the examples shown here, the catchment outflows, as visible in Figures 10 and 12, show a higher responsiveness to variations in initial moisture at the beginning of the simulations, and attain higher peak values, when using the post-2023 TUFLOW releases.
 
[[File:Init moisture F10.png|600px|Figure 9: Sensitivity of cumulative infiltration in the Plynlimon Gwy catchment to the initial moisture parameter in the Green-Ampt infiltration model.|border]]<br>
'''Figure 9: Sensitivity of cumulative infiltration in the Plynlimon Gwy catchment to the initial moisture parameter in the Green-Ampt infiltration model.'''<br>
[[File:InitInitial moisture F11 catchment.pngjpg|border|600px760x760px|Figure 10: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the initial moisture parameter in the Green-Ampt infiltration model in pre-2023 release of TUFLOW.|border]]<br>
'''Figure 10: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the initial moisture parameter in the Green-Ampt infiltration model in pre-2023 release of TUFLOW.'''<br>
[[File:Fig10Init GAmoisture soilspost 2023 v2.pngjpg.jpg|border|600px760x760px|Figure 11: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDAinitial soil typemoisture parameter in the Green-Ampt infiltration model in post-2023 release of TUFLOW.|border]]<br>
'''Figure 11: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDAinitial soil typemoisture parameter in the Green-Ampt infiltration model in post-2023 release of TUFLOW.'''
 
=== Max Ponding Depth ===
Line 149 ⟶ 151:
 
== In built USDA soil type ==
The model was also run with the default in-build USDA soil types. Figure 1112 shows the outputs. As expected the higher the soil type, then typically the more the infiltration and the lower the produced runoff. Soils 8-11, which represent sandy soils do not show any runoff in this example as the rainfall applied directly to the mesh is all infiltrated.<br>
<br>
[[File:Fig10 GA soils.png|600px|Figure 11: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDA soil type parameter in the Green-Ampt infiltration model.|border]]<br>
'''Figure 11: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDA soil type parameter in the Green-Ampt infiltration model.'''
<br>
[[File:Soil Type.jpg|border|760x760px|Figure 12: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDA soil type parameter in the Green-Ampt infiltration model.]]<br>
'''Figure 12: Sensitivity of simulated flow at the Cefn-Brwn gauge location in the Plynlimon Gwy catchment to the USDA soil type parameter in the Green-Ampt infiltration model.'''
 
== Summary ==
The Green-Ampt infiltration model is one of the infiltration methods available within TUFLOW. There is extensive literature on its application, including suggested parameter values for various soil types, though these are mostly based on soils in the United States.
The Green-Ampt infiltration model is one of the available infiltration models within TUFLOW. There is a wide range of literature on Green-Ampt applications and some suggested parameter values for particular soil types, albeit mostly soil types from the US. The 3 main Green-Ampt parameters have been tested to show the sensitivity of model outputs to the values as well as the variation in initial moisture. The results show that the model is relatively insensitive to the porosity value and suction head parameter. However, the outputs do show significant variations in runoff volume due to variations in both hydraulic conductivity. As part of any calibration exercise it is suggested that it is the hydraulic conductivity, in conjunction with the initial moisture content, which would be the parameters that are most focused on. The hydraulic conductivity appears to affect the runoff volume throughout the event whereas the initial soil moisture has a limited impact at the beginning of the event before soils become saturated and results converge.<br>
 
Three main Green-Ampt parameters have been tested to assess the sensitivity of model outputs to parameter values and variations in initial soil moisture. The results show that the model is relatively insensitive to the porosity and suction head parameters. However, outputs show significant variations in runoff volume in response to changes in hydraulic conductivity.
 
As part of any calibration process, it is recommended that hydraulic conductivity and initial moisture content be prioritised during calibration. Hydraulic conductivity influences runoff volume throughout the event, while initial soil moisture mainly affects the early part of the simulation until soils become saturated and results converge.
 
==Acknowledgements==