Groundwater Modelling Advice Draft

Page Under Construction

Common Questions Answered (FAQ)

Why is the groundwater model affecting areas beyond the expected flow path?

Groundwater movement in TUFLOW is influenced by topography, soil properties, and groundwater parameters. If groundwater is introduced, infiltrated water may reappear on the surface depending on these factors.

Key considerations:

• Without horizontal groundwater movement, infiltrated water is lost from the system. When groundwater movement is enabled, water can accumulate in low-lying areas or flow towards the downstream boundary.
• Unexpected groundwater presence in areas without assigned soil depths may be due to groundwater behaviour at boundaries or model-wide settings.
• The 2d_po regional outputs introduced in version AF can assist with analysing groundwater movement and verifying model behaviour.

Should soil layers be assigned across the entire model?

If groundwater is being simulated, defining soil layers across the full model domain can provide more control over groundwater behaviour.

Options include:

• Assigning a depth of zero in areas where groundwater should not be present.
• Setting a large depth in areas where groundwater storage is needed but horizontal transmission is not desired.

Model validation using observed data is recommended to confirm that groundwater interactions align with real-world conditions.

How can a lined filtration trench be modelled in TUFLOW?

A lined filtration trench can be represented in TUFLOW using one of the following methods:

1D Channel Approach

• Model the trench as a 1d_nwk “Q” type structure, which uses a depth-discharge relationship.
• Connect the trench to multiple 2D cells using a 2d_bc SX line.
• If the trench is long, divide it into multiple sections to improve accuracy.

1D Pit and Culvert Approach

• Use multiple 1d_nwk Q pits connected via a 1d_nwk culvert.
• The connected 2D cells are automatically selected using the sag or on-grade method (see section 5.12.3.3 of the TUFLOW manual).
• To manually control cell connections, set the 1d_nwk Conn_No attribute to a negative value. For example, Conn_No = -1 ensures each pit connects to only one 2D cell.
• This method allows the trench to be represented as multiple pit points, each selecting one 2D cell and linking to the 1D node using X connectors.

Additional Considerations

• Pit Inlet Discharge Curve: If using a 1D network, pre-compute the discharge for various depths to define a suitable pit inlet discharge curve.
• Interflow Functionality: The use of interflow depends on the cell size relative to the trench feature.
• Two-Stage Modelling Approach: One method involves running a 2D infiltration simulation first to determine infiltration rates. The infiltration can then be applied as a 1D boundary condition in a second simulation to represent flow into the drainage system.

TUFLOW is currently developing functionality to dynamically link interflow to a 1D network, which will allow infiltrated subsurface flows to connect directly to 1D nodes or pipes in future releases.

What methods and result outputs can be used to quantify infiltration losses over a given area in TUFLOW?

To assess infiltration losses, the following outputs can be used:

Map Output Data Types

• CI (Cumulative Infiltration): Displays the total infiltration over time.

• IR (Infiltration Rate): Shows the infiltration rate at each timestep.

Point Output (2d_po) Approach

• If infiltration data is needed for a specific area, a 2d_po region can be set up.

Both methods help understand how infiltration happens across different areas.

Why does changing the initial soil moisture in the Green-Ampt (GA) infiltration method not affect infiltration rates as expected?

In the 2023 release of TUFLOW, the treatment of initial soil moisture in the GA method was updated to accommodate horizontal soil water movement. This change allows soil moisture to flow between sub-surface cells and be removed via evapotranspiration, making it more dynamic over time. Instead of remaining static in the infiltration equation, the initial soil moisture is now used as the cumulative infiltration (F) at the first timestep of the simulation. This means:

• The infiltration rate at the start of the simulation is influenced by the initial soil moisture, but subsequent infiltration behaviour depends on the balance of infiltration, drainage, and evapotranspiration.

• If no soil thickness is defined, the model assumes dZ = 0, meaning the initial moisture value is ignored altogether.

Suggested Workarounds

1. Reverting to the pre-2023 method.

Use the backward compatibility switch:

Defaults == Pre 2023

Note: This changes multiple default settings, not just the GA method.

2. Adjusting the soil porosity in the .tsoil file.

Instead of defining initial soil moisture separately set: Adjusted Porosity = Porosity - Initial Moisture

Example: If soil porosity is 0.385 and initial soil moisture is 0.200, set the porosity to 0.185.

TUFLOW is updating documentation to clarify this change and is working on a dedicated command to allow users to revert to the original GA method without affecting other default settings. Further refinements are also being considered for long-term simulations involving multiple wet/dry cycles.

If you have feedback or specific cases where this affects your modelling, please let contact TUFLOW Support via support@tuflow.com.

Can TUFLOW trap groundwater beneath impervious areas to prevent unrealistic exfiltration?

Currently, TUFLOW does not directly simulate pressurised groundwater flow beneath impervious surfaces. The existing groundwater model calculates exfiltration based on mass balance, meaning that groundwater can still migrate upward even in areas where the surface is defined as impervious.

Recognising the importance of this feature for improving groundwater representation in urban modelling, consideration is being given to incorporating a method to account for trapped groundwater as part of future TUFLOW release.

In the meantime, users may mitigate this issue by:

• Removing impervious locations from the soil layer entirely.

• Setting horizontal hydraulic conductivity to zero for soil polygons beneath impervious areas.

• Adjusting the soil thickness and properties to reduce unrealistic seepage effects.

If this issue is recurring or if there are specific use cases where this feature would be beneficial, feedback is encouraged via support@tuflow.com.

What is the recommended method for representing a railway ballast area in TUFLOW?

A railway ballast area in TUFLOW can be represented as a soil layer (specific soil ID) with high infiltration, suitable porosity, and high hydraulic conductivity in both horizontal and vertical directions.

A global soil thickness value can define the depth of this layer relative to the surface elevation. Similarly, a global soil base elevation value beneath the railway ballast area can set the absolute elevation of its bottom. Both parameters can also be varied spatially with GIS and/or grid layers.