Quick links to often-used pages:
Quick links to often-used pages:
There are two model options for processes deep below ground; regional groundwater flow and aquifers. The regional groundwater flow model simulates groundwater flow between subbasins, but no explicit storage of deep groundwater. The aquifer model simulates aquifer storage and delay before regional groundwater reaches the destination.
Regional groundwater flow from a subbasin is calculated before the subbasin and class loops. This water is added to the downstream subbasin (in the class-loop part of the code) and/or to the outlet lake of the present subbasin (in the routing calculation part of the code). A subbasin dependent variable (grwtolake) governs the allocation between them. The outflow from a subbasin is calculated as a percentage of the water-retention capacity with a recession coefficient (rcgrw) e.g. for the lowest soil layer (k) the groundwater outflow for a class becomes:
Groundwater flows are summarized for all classes in proportion to their areas and converted to volume flow out of the subbasin.
The groundwater flow to be added to the bottom layer (and mixed) is added in the class loop portion of the code. If the layer then becomes full of water, it will be pressed into the layer above and mixed, etc. Please note that if downstream subbasin is much smaller than the upstream subbasin, groundwater inflows can become large. If part of the groundwater flow is to go to the outlet lake in the subbasin this is added to the volume of the lake in the routing part of the code.
With an input variable in GeoData.txt, an outlet lake may be flagged to receive inflow from the regional groundwater. The coefficient (grwtolake) specifies the portion of the produced regional groundwater flow which is led to the subbasins outlet lake.
|regional_groundwater_module (regional_groundwater.f90)||calculate_regional_groundwater_flow||soil regional groundwater flow creation|
|calculate_soillayer_groundwaterflow_removal||soil regional groundwater flow creation|
|add_regional_groundwater_flow_to_soil||soil addition of regional groundwater flow|
|add_regional_groundwater_flow_to_olake||regional groundwater to outlet lake|
Deep percolation is determined by a recession coefficient (). The recession coefficient is set based on the recession parameter (rcgrwst) which is soil type dependent, but if it is not set the general parameter (rcgrw) is used instead. The recession coefficient can also be regionally adjusted by the correction parameter (aqpercorr).
Deep percolation from lowest soil layer in model is determined by the recession coefficient and the available water for runoff in the soil layer (qperc (mm/timestep)).
Deep percolation from main river in model is determined by the recession coefficient and the water in the river (qperc (m3/timestep)).
The water is removed from river water and the queue in proportion to their relative volumes.
Deep percolation from soil and rivers in all recharge subbasins of an aquifer is added together. This total deep percolation () is delayed before it recharges the aquifer ().
The delay property (delayts, in timesteps per day) is specified for each aquifer (delay in days). It can be adjusted by the parameter region dependent model parameter aqdelcor.
Outflow from the aquifer depends on the water in the aquifer (vol) and a rate. The rate is determined by a aquifer property (retrate) and a regional adjustment parameter (aqretcor, parameter region dependent).
The outflow may be divided between several receiving subbasins (retfrac). The water is added to the inflow of the main river of these subbasins.
The return flow is added to the main river of the receiving subbasin.
The organic nitrogen, particulate phosphorus, and organic carbon in water percolating from soil and river stay in the soil/river, while temperature, inorganic nitrogen and SRP follow the water to the aquifer. Inorganic nitrogen is subjected to denitrification in the aquifer (see Chapter Nitrogen and phosphorus in land routines - Denitrification).
|Aquifers||aquid, subid, parreg||AquiferData.txt|
|basedepth, inidepth, porosity, area|
|Deep percolation||calculated from||rcgrw, rcgrwst, aqpercorr||par.txt|
|wp||wcwp, wcwp1, wcwp2, wcwp3|
|fc||wcfc, wcfc1, wcfc2, wcfc3|
|Deep percolation delay and aquifer recharge||recharge||AquiferData.txt|
|Return flow from aquifer||retrate, retfrac||retrate, retfrac||AquiferData.txt|
|Nutrient processes related to aquifers||topdepth, temp, conc_IN, conc_SP||AquiferData.txt|
|denitaq, hsatINs, wcep, wcfc, wcwp||par.txt|
|calculate_percolation_delay||deep percolation delay and aquifer recharge|
|calculate_aquifer_outflow||return flow from aquifer|
|calculate_current_flow_from_aquifer||return flow from aquifer|
|nutrient processes related to aquifers|
|set_percolation_concentration||nutrient processes related to aquifers|