Garden with Insight v1.0 Help: Nitrogen - Volatilization

Volatilization, the loss of ammonia to the atmosphere, is estimated simulataneously with nitrification. Volatilization of surface-applied ammonia is estimated as a function of temperature and wind speed using the equation [Equation 194] where AKV is the volatilization regulator from equation 185 and WNF is the wind speed factor for surface application (soil layer 1). (What this means is that a different set of equations is used for the first soil layer than for the other soil layers.)

Equation 194

AKV = TF * WNF
Code:
same
Variables:
AKV = VolatilRegulatorForSurfaceLayer
TF = nitrifAndVolatilTempFactor
WNF = volatilWindSpeedFactorForSurface

The wind speed factor for surface volatilization is estimated with this equation [Equation 195] where V is the mean wind speed in m/sec for the day.

Equation 195

WNF = 0.335 + 0.16 * ln(V)
Code:
same
Variables:
WNF = VolatilWindSpeedFactorForSurface
V = meanWindSpeedForDay_mPsec

Depth of ammonia within the soil, cation exchange capacity of the soil, and soil temperature are used in estimating below surface volatilization [Equation 196] where CECF is the cation exchange capacity factor.

Equation 196

AKV = TF * CECF * DPF
Code:
same
Variables:
AKV = VolatilRegulatorForNonSurfaceLayer
TF = nitrifAndVolatilTempFactor
CECF = volatilCationExchangeCapFactor
DPF = volatilDepthFactor

The cation exchange capacity factor, CECF, is estimated with the equation [Equation 197] where CEC is the cation exchange capacity for soil layer l.

Equation 197

CECF = 1.0 - 0.038 * CEC, CECF >= 0
Code:
same
Variables:
CECF = VolatilCationExchangeCapFactorForLayer
CEC = cationExchangeCapacity_cmolPkg

The depth factor, DPF, is computed with the equation [Equation 198] where ZZ is the depth to the middle of soil layer l in mm. The parameters in equation 198 are set to give DPF = 0.95 at ZZ = 5 mm and DPF = 0.05 at ZZ = 100 mm.

Equation 198

DPF = 1.0 - ZZ / (ZZ + exp(4.706 - 0.0305 * ZZ))
Code:
same
Variables:
DPF = VolatilDepthFactorForLayer
ZZ = depthToMiddleOfSoilLayer_mm

The depth to the middle of the soil layer is computed with [Equation 199] where Z is the depth to the bottom of layer l in mm.

Equation 199

ZZ = 1000 * (Z(l) + Z(l-1)) / 2
Code:
same
Variables:
ZZ = DepthToMiddleOfSoilLayer_mm
Z(l) = depth_m
Z(l-1) = depthUpperLayer_m

To partition nitrification and volatilization rates appropriately, equation 185 is solved for each process separately [Equation 200] and [Equation 201] where RNIT and RVOL are the separate estimates of nitrification and volatilization. However, to satisfy equation 185, equations 200 and 201 must be adjusted. Thus, the final equations [Equation 202] and [Equation 203].

Equation 200

RNIT = 1.0 - exp(-AKN)
Code:
same
Variables:
RNIT = UnadjNitrifForLayer_kgPha
AKN = nitrifRegulatorForLayer

Equation 201

RVOL = 1.0 - exp(-AKV)
Code:
same
Variables:
RVOL = UnadjVolatilForLayer_kgPha
AKV = volatilRegulator

Equation 185

RNV = WNH3 * (1.0 - exp(-AKN - AKV))
Code:
same
Variables:
RNV = CombinedNitrifAndVolatilForLayer_kgPha
WNH3 = ammonia_kgPha
AKN = nitrifRegulator
AKV = volatilRegulator

Equation 202

RNIT* = RNV - RVOL*
Code:
RNIT* = RNV * RNIT / (RVOL + RNIT)
(nitrif and volatil are switched, but it comes to the same thing)
Variables:
RNIT* = NitrificationForLayer_kgPha
RNV = combinedNitrifAndVolatil_kgPha
RNIT = unadjNitrif_kgPha
RVOL = unadjVolatil_kgPha

Equation 203

RVOL* = RNV * RVOL / (RVOL + RNIT)
Code:
RVOL* = RNV - RNIT*
(nitrif and volatil are switched, but it comes to the same thing)
Variables:
RVOL* = VolatilizationForLayer_kgPha
RNV = combinedNitrifAndVolatil_kgPha
RNIT* = nitrification_kgPha