If nitrate levels in the root zone are insufficient to meet the demand of a legume, SWAT+ allows the plant to obtain additional nitrogen through nitrogen fixation. Nitrogen fixation is calculated as a function of soil water, soil nitrate content and growth stage of the plant.
Nfixβ=Ndemandββfgrββmin(fswβ,fno3β,1) 5:2.3.9
where Nfixβ is the amount of nitrogen added to the plant biomass by fixation (kg N/ha), Ndemandβ is the plant nitrogen demand not met by uptake from the soil (kg N/ha), fgrβ is the growth stage factor (0.0-1.0), fswβ is the soil water factor (0.0-1.0), and fno3β is the soil nitrate factor (0.0-1.0). The maximum amount of nitrogen that can be fixed by the plant on a given day is Ndemandβ.
Growth stage exerts the greatest impact on the ability of the plant to fix nitrogen. The growth stage factor is calculated:
when 5:2.3.10
when 5:2.3.11
when 5:2.3.12
when 5:2.3.13
when 5:2.3.14
where is the growth stage factor and is the fraction of potential heat units accumulated for the plant on a given day in the growing season. The growth stage factor is designed to reflect the buildup and decline of nitrogen fixing bacteria in the plant roots during the growing season.
The soil nitrate factor inhibits nitrogen fixation as the presence of nitrate in the soil goes up. The soil nitrate factor is calculated:
when 5:2.3.15
when 5:2.3.16
when 5:2.3.17
where is the soil nitrate factor and is the nitrate content of the soil profile (kg NO-N/ha).
The soil water factor inhibits nitrogen fixation as the soil dries out. The soil water factor is calculated:
5:2.3.18
where is the soil water factor, is the amount of water in soil profile (mm HO), and is the water content of soil profile at field capacity (mm HO).
Plant nitrogen uptake is controlled by the plant nitrogen equation. The plant nitrogen equation calculates the fraction of nitrogen in the plant biomass as a function of growth stage given optimal growing conditions.
5:2.3.1
where is the fraction of nitrogen in the plant biomass on a given day, is the normal fraction of nitrogen in the plant biomass at emergence, is the normal fraction of nitrogen in the plant biomass at maturity, is the fraction of potential heat units accumulated for the plant on a given day in the growing season, and and are shape coefficients.
The shape coefficients are calculated by solving equation 5:2.3.1 using two known points (, ) and (, ):
n2β=frPHU,100%ββfrPHU,50%β(1n[(1β(frN,1ββfrN,3β)(frN,2ββfrN,3β)β)frPHU,50%βββfrPHU,50%β]β1n[(1β(frN,1ββfrN,3β)(frN,βΌ3ββfrN,3β)β)frPHU,100%βββfrPHU,100%β])β 5:2.3.3
where n1β is the first shape coefficient, n2β is the second shape coefficient, frN,1β is the normal fraction of nitrogen in the plant biomass at emergence, frN,2β is the normal fraction of nitrogen in the plant biomass at 50% maturity, frN,3β is the normal fraction of nitrogen in the plant biomass at maturity, frN,βΌ3β is the normal fraction of nitrogen in the plant biomass near maturity, frPHU,50%β is the fraction of potential heat units accumulated for the plant at 50% maturity (frPHU,50%β=0.5), and frPHU,100%β is the fraction of potential heat units accumulated for the plant at maturity (frPHU,100%β=1.0). The normal fraction of nitrogen in the plant biomass near maturity (frN,βΌ3β) is used in equation 5:2.3.3 to ensure that the denominator term (1β(frN,1ββfrN,3β)(frN,βΌ3ββfrN,3β)β) does not equal 1. The model assumes (frN,βΌ3ββfrN,3β)=0.00001
To determine the mass of nitrogen that should be stored in the plant biomass on a given day, the nitrogen fraction is multiplied by the total plant biomass:
bioN,optβ=frNββbio 5:2.3.4
where bioN,optβ is the optimal mass of nitrogen stored in plant material for the current growth stage (kg N/ha), frNβ is the optimal fraction of nitrogen in the plant biomass for the current growth stage, and bio is the total plant biomass on a given day (kg haβ1).
Originally, SWAT+ calculated the plant nitrogen demand for a given day by taking the difference between the nitrogen content of the plant biomass expected for the plantβs growth stage and the actual nitrogen content Nupβ=bioN,optββbioNβ. This method was found to calculate an excessive nitrogen demand immediately after a cutting (i.e. harvest operation). The equation used to calculate plant nitrogen demand is now
Nupβ=Min{bioN,optββbioNβ4βfrN,3ββΞbioβ 5:2.3.5
where Nupβ is the potential nitrogen uptake (kg N/ha), bioN,optβ is the optimal mass of nitrogen stored in plant material for the current growth stage (kg N/ha), bioNβ is the actual mass of nitrogen stored in plant material (kg N/ha), frN,3β is the normal fraction of nitrogen in the plant biomass at maturity, and Ξbio is the potential increase in total plant biomass on a given day (kg/ha).
The depth distribution of nitrogen uptake is calculated with the function:
Nup,zβ=[1βexp(βΞ²nβ)]Nupβββ[1βexp(βΞ²nββzrootβzβ)] 5:2.3.6
where Nup,zβ is the potential nitrogen uptake from the soil surface to depth z (kg N/ha), Nupβ is the potential nitrogen uptake (kg N/ha), Ξ²nβ is the nitrogen uptake distribution parameter, z is the depth from the soil surface (mm), and zrootβ is the depth of root development in the soil (mm). Note that equation 5:2.3.6 is similar in form to the depth distribution for water uptake described by equation 5:2.2.1. The potential nitrogen uptake for a soil layer is calculated by solving equation 5:2.3.6 for the depth at the upper and lower boundaries of the soil layer and taking the difference.
Nup,lyβ=Nup,zlββNup,zuβ 5:2.3.7
where Nup,lyβ is the potential nitrogen uptake for layer ly (kg N/ha), Nup,zlβ is the potential nitrogen uptake from the soil surface to the lower boundary of the soil layer (kg N/ha), and Nup,zuβ is the potential nitrogen uptake from the soil surface to the upper boundary of the soil layer (kg N/ha).
Root density is greatest near the surface, and nitrogen uptake in the upper portion of the soil will be greater than in the lower portion. The depth distribution of nitrogen uptake is controlled by Ξ²nβ, the nitrogen uptake distribution parameter, a variable users are allowed to adjust. Figure 5:2-4 illustrates nitrogen uptake as a function of depth for four different uptake distribution parameter values.
Nitrogen removed from the soil by plants is taken from the nitrate pool. The importance of the nitrogen uptake distribution parameter lies in its control over the maximum amount of nitrate removed from the upper layers. Because the top 10 mm of the soil profile interacts with surface runoff, the nitrogen uptake distribution parameter will influence the amount of nitrate available for transport in surface runoff. The model allows lower layers in the root zone to fully compensate for lack of nitrate in the upper layers, so there should not be significant changes in nitrogen stress with variation in the value used for Ξ²nβ.
The actual amount if nitrogen removed from a soil layer is calculated:
Nactualup,lyβ=minβNup,lyβ+Ndemandβ,NO3lyββ 5:2.3.8
where Nactualup,lyβ is the actual nitrogen uptake for layer ly (kg N/ha), Nup,lyβ is the potential nitrogen uptake for layer ly (kg N/ha), Ndemandβ is the nitrogen uptake demand not met by overlying soil layers (kg N/ha), and NO3lyβ is the nitrate content of soil layer ly (kg NO3β-N/ha).
