Sep 30, 2010

Agricultural Lands Assessment of Nitrogen (ALAN) is a one-dimensional process-oriented model designed to describe in reasonable detail the fate of water and nitrogen (N) in croplands at short time-steps (i.e. minutes to hours) using easy-to-obtain input data.  ALAN integrates calculations made by eight sub-models as depicted in Figure 1. The Canopy-Transpiration-Photosynthesis (CTP) sub-model originally simulated crop transpiration, net leaf photosynthesis and stomatal conductance based on a two-leaf canopy model (Kremer et al., 2008), but was extended to include soil sensible and latent heat fluxes, crop senescence, and water interception by the crop. The crop growth sub-model simulates leaf area index and biomass accumulation based on photosynthesis and respiration rates, but restricted by shortages of soil N. Photosynthates are partitioned into leaves, stems, roots and storage organs, based on the crop development stage, which is driven by the accumulation of growing degree days. Crop N uptake is computed by a supply/demand scheme. The soil organic carbon (C) sub-model simulates the turnover of soil organic matter, manure amendments and crop residue additions, CO2 production, and mineralization and immobilization of inorganic N using a one-pool model of soil organic matter as proposed by Kemanian and Stöckle (2010). The urea hydrolysis sub-model computes hydrolysis of urea based on soil organic carbon, soil pH and soil moisture as proposed by Goodwin and Jones (1991). The nitrification sub-model simulates oxidation of soil ammonium (NH4+) with production of nitric oxide (NO) and (nitrous oxide) N2O, based on a first-order process controlled by soil pH, soil temperature, and soil moisture. The denitrification sub-model computes reduction of soil nitrate (NO3-) with production of N2O and dinitrogen (N2). The approach is similar to that of Parton et al. (1996) and del Grosso et al. (2000), in which the degree of reduction is controlled by the soil moisture as an indicator of anoxic conditions, soil NO3- (electron acceptors) and soil respiration rate as an indicator of labile C levels (electron donors). 
Figure 1.  Simplified flow diagram of ALAN model
The sub-model for the transport of water, solutes, gases and heat in the soil profile solves governing partial differential equations by a fully-implicit, cell-centered finite difference formulation. Water, solutes and gases are solved together in a nested, variable-time procedure. Soil water transport is simulated by solving the Richards’ equation in a mass-conserving procedure that couples the modified Picard iteration scheme developed by Celia et al. (1990) with the internodal hydraulic conductivity algorithm developed by Romano et al. (1998). Soil runoff is computed by setting a soil surface ponding capacity, while deep percolation is calculated by tracking the water flux at the lower boundary of the soil profile. Soil NO3- is considered to be transported in the soil only by convection, assuming that numerical dispersion substitute, to some extent, for physical dispersion. Transport of NH4+-N is simulated by a common formulation of the convection equation for aqueous NH4+ coupled with a diffusion equation for gaseous ammonia (NH3). Solution-phase and adsorbed NH4+ are related by a Langmuir isotherm. Convective transport of both NO3- and NH4+ make use of a backward-time, upwind-type scheme. Transport of gaseous NO and N2O is simulated by a diffusion equation, using a source term for the production of NO or N2O due to nitrification and/or denitrification processes. Gas diffusion coefficients are computed by three optional models developed and evaluated by Moldrup et al. (1997, 1999, 2004). Soil heat transport is simulated by the Fourier’s law of heat conduction according the procedure developed by Campbell (1985).

The NH3 exchange sub-model uses information from other sub-models to calculate the network of NH3 fluxes in the soil-plant-atmosphere system as shown in Figure 2. Sub-stomatal NH3 concentration is calculated according to Farquhar et al. (1980) by relating the ratio of apoplastic [NH4+]/[H+] to the plant N concentration, based on studies by Husted et al. (1996), Schjoerring et al. (1998), and Mattsson and Schjoerrinig (2002). NH3 is allowed to be deposited on the leaf cuticle based on the model of Sutton and Fowler (1993), but is washed off by precipitation and irrigation events. Finally, ALAN is driven by meteorological variables and management practices. Daily meteorological variables are distributed into short-term values according to the procedure described by Kremer (2006). Daily rainfall, however, is disaggregated into short-time intensity values using a modification of the model developed by Arnold and Williams (1989). Management practices considered are tillage operations, manure amendments, and applications of water and N-fertilizers. Inorganic N can be applied as urea, NH4+- and NO3--based fertilizers, on the soil surface, incorporated into the soil, or distributed by the irrigation water.  
Figure 2.  Electrical analog of NH3 exchange in the soil-plant-atmosphere system as described by ALAN model


References 
Arnold, J.G. and Williams, J.R. (1989). “Stochastic generation of internal storm structure at a point.” Trans. ASAE 32, 161–167.
Campbell, G.S. (1985). “Soil physics with BASIC, transport models for soil-plant systems.” Developments in Science 14. Elsevier. New York.
Celia, M.A., Bouloutas, E.T. and Zarba, R.L. (1990). “A general mass-conservative numerical solution for the unsaturated flow equation.” Water Resour. Res. 26, 1483–1496.
del Grosso, S.J., Parton, W.J., Mosier, A.R., Ojima, D.S., Kulmala, A.E. and Phongpan, S. (2000). “General model for N2O and N2 gas emissions from soils due to denitrification.” Global Biochem. Cycles 14, 1045–1060.
Farquhar, G.D., Firth, P.M., Wetselaar, R. and Weir, B. (1980). “On the gaseous exchange of ammonia between leaves and the environment: determination of ammonia compensation point.” Plant Physiol. 66, 710–714.
Godwin, D.C. and Jones, C.A. (1991). “ Nitrogen dynamics in soil-plant systems.” In:  Hanks, J. and Ritchie, J.T. (eds.). Modeling plant and soil systems. Agronomy Series 31. American Society of Agronomy, Crop Science Society of America & Soil Science Society of America. Madison, Wisconsin. pp, 287–339
Husted, S., Mattson, M. and Schjoerring, J.K. (1996). “Ammonia compensation points in two cultivars of Hordeum vulgare L. during vegetative and generative growth.” Plant Cell Environ. 19, 1299–1306.
Kemanian, A.R. and Stöckle, C.O. (2010). “C-FARM: a simple model to evaluate the carbon balance of soil profiles.” Eur. J. Agron. 32, 22–29.
Kremer, C. (2006). “Evaluating simple transpiration-based models of crop productivity.” Ph.D. Dissertation. Department of Biological Systems Engineering. Washington State University. Pullman.
Kremer, C., Stöckle, C.O., Kemanian, A.R. and Howell, T. (2008). “A canopy transpiration and photosynthesis model for evaluating simple productivity models.” In:  Ahuja, L (ed.).   Responses to crops to limited water: understanding and modeling water stress on plant growth processes. Advances in Agricultural Systems Modeling Series 1. American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America. pp, 165–189.
Mattson, M. and Schjoerring, J.K. (2002). “Dynamic and steady-state responses of inorganic nitrogen pools and NH3 exchange in leaves of Lolium perenne and Bromus erectus to changes in root nitrogen supply.” Plant Physiol.128, 742–750.
Moldrup, P., Olesen, T., Rolston, D.E. and Yamaguchi, T. (1997). “Modeling diffusion and reaction in soils: VII. Predicting gas and ion diffusivity in undisturbed and sieved soils.” Soil Sci. 162, 632–640.
Moldrup, P,., Olesen, T., Yamaguchi, T., Schjønning, P. and Rolston DE (1999). “ Modeling diffusion and reaction in soils: IX. The Buckingham-Burdine-Campbell equation for gas diffusivity in undisturbed soils.” Soil Sci. 164, 542–551.
Moldrup, P., Olesen, T., Toshikawa, S., Komatsu, T. and Rolston, D.E. (2004). “Three-porosity model for predicting the gas diffusion coefficient in undisturbed soil.” Soil Sci. Soc. Am. J. 68, 750–759.
Parton, W.J., Mosier, A.R., Ojima, D.S., Valentine, D.W., Schimel, D.S., Weier, K. and Kumala, A.E. (1996). “Generalized model for N2 and N2O production from nitrification and denitrification.” Global Biochem. Cycles 10, 401–412.
Romano, N., Brunone, B. and Santini, A. (1998). “Numerical analysis of one-dimensional unsaturated flow in layered soils.” Adv. Water Resour. 21, 315–324.
Schjoerring, J.K., Husted, S. and Mattson, M. (1998). “Physiological parameters controlling plant-atmosphere ammonia exchange.” Atmos. Environ. 32, 491–498.
Sutton, M.A. and Fowler, D. (1993). “A model for inferring bidirectional fluxes of ammonia over plant canopies.” In: Proceedings of the WMO conference on the measurement and modelling of atmospheric composition changes including pollution transport. Sofia, Bulgaria. October 4–8. World Meteorological Organization, Global Atmosphere Watch Report 91. Geneva, Switzerland. pp, 179–182.