Conceptual Overview of Salinity Modeling
The shallow groundwater is not always well mixed in terms of salinity. In some areas, the upper portion of the aquifer interacts quickly with the surface, but further down, the interaction is slower. This process can be modeled as two fully mixed layers within the groundwater object.
The current implementation of the ground water quality model allows you to model salinity with the Layered Salt method.
Figure 11.1 shows the groundwater object with two layers in terms of flow and volume slots.
Figure 11.2 shows the two layers in terms of salinity slots. All slots and links shown in (see
“solveGWMB_givenPreviousElevations” in Objects and Methods) are still present; the object just has additional slots and links for salinity modeling.
Figure 11.1 Slots on the groundwater object to model flow and volume with two layers
Figure 11.2 Slots on the groundwater object to model salinity with two layers
In this approach, the thickness of the upper layer is constant. Each mass balance component except Evaporation and Evapotranspiration have an equivalent Salt Concentration component. They are shown in red in
Figure 11.1. On the groundwater object, the Flow Upstream/Downstream/Left/Right, Inflow From Surface Water, Pumped Flow and Percolation could be either negative or positive.When the flow is coming into the groundwater object, the salt concentration comes from the linked slot or be specified (input or rules). If the flow is out of the groundwater object, the concentration is computed by the groundwater object. When the groundwater object is defining the concentration, the previous values are used to avoid iteration problems.
Following are assumptions used in this computation:
• Inflow from Surface Water is added to (or removed from) the upper layer.
• Specified Inflow is added to either layer based on a user specified proportion.
• Evaporation is removed from the upper layer and no salt is removed.
• Evapotranspiration is removed from the upper layer and no salt is removed.
• Percolation is removed from (or added to) the lower layer.
• Pumped Flow is removed from (or added to) either layer as specified by the user.
Following is a conceptual overview of the computation to determine storage and salt concentration in each layer. See
“Dispatch Method - Solve 2 Layer Salt” for details on this computation.
1. Determine the layer storages and the proportion based on the capacity of the upper layer.
2. Compute the lateral flow and salt flux from each layer based on the storage proportion on the given object or adjacent objects.
3. Compute pumped flow proportion and salt flux.
4. Compute inflow from surface water salt flux.
5. Compute specified inflow salt flux.
6. Compute percolation salt flux.
7. Compute the intermediate layer storage and concentration based on the computed incoming and outgoing components.
8. By comparing the intermediate and final storage and salt concentration, transfer water and salt between the two layers.
9. Compute the salt concentration of each layer.
10. Set salt mass slots if method is selected.
Note: This conceptual design is for two layers. It can be reduced to one layer by specifying a large Upper Thickness.