Temporal moments analysis of preferential solute transport in soils

Temporal moments analysis of solute breakthrough curves is used to investigate the preferential leaching of chloride, nitrate and phosphate through an Australian soil. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and non-equilibrium between soil constituents, water and solutes. Therefore simple models are required to accurately characterise solute transport in natural and agricultural soils under non-equilibrium conditions. A multiple sample percolation system, consisting of 25 individual collection wells was constructed to study the effects of localised soil heterogeneities on the transport of nutrients (NO₃ ^–, Cl^–, PO₄ ^{3 –}) in the vadose zone of an agricultural soil predominantly dominated by clay. Using data collected from the multiple sample percolation experiments, this paper compares and contrasts the performance of temporal moments analysis with two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR) and a two-region preferential flow model (TRM) suitable for modelling preferential transport. The values for solute transport parameters predicted by temporal moments analysis were in excellent agreement with experimental data and results from ADR and TRM. It is concluded that temporal moments analysis when applied with other physical models such as the ADR and TRM, provide an excellent means of obtaining values for important solute transport parameters and gaining insight of preferential flow. These results have significant ramifications for modelling solute transport and predicting nutrient loadings.     

A mathematical model to estimate errors associated with closed flux chambers

Errors associated with the closed flux chamber technique, used to measure surface emissions from landfills, were investigated by using a combination of numerical modeling and laboratory studies. A transient-state, advective–dispersive–reactive model was developed and used in conjunction with its steady-state version to quantify the errors associated with closed flux chambers. In developing the model, all four major gases, CH₄, O₂, CO₂, and N₂, and the oxidation of CH₄ to CO₂ were considered. Laboratory experiments were conducted on a monolayered as well as a two-layered landfill cover system to calibrate and verify the model. The model was used to develop a plot of the percentage errors associated with closed flux chambers of different dimensions and surface flux rates.     

Adsorption–desorption of metolachlor and atrazine in Indian soils: effect of fly ash amendment

The effect of two fly ashes as soil amendment on the adsorption–desorption of metolachlor [2-chloro-N-(2-ethyl-6-methylphenyl)-N-(2-methoxy-1-methylphenyl)] and atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine) was studied in alluvial and laterite soils. The adsorption data for both the herbicides fitted well the Freundlich equation, and Freundlich adsorption coefficient (K _f) increased with an increase of fly ash amount. Both the fly ashes differed in their extent to increase herbicide sorption, and the effect was different in different soils. Atrazine was sorbed more in the soils/soils?+?fly ash mixtures than the metolachlor. The K _f values showed significant correlation with the amount of fly ash amendment (correlation coefficient, R?>?0.982). The desorption isotherms also fitted the Freundlich equation, and desorption showed hysteresis which increased with an increase in the content of fly ash amendment. The free energy change (ΔG) indicated that the sorption process is exothermic, spontaneous, and physical in nature. The study has shown that fly ash as soil amendment significantly increased the sorption of metolachlor and atrazine, but the effect is soil- and fly ash-specific.     

Atrazine Transport Within a Coastal Zone in Southeastern Puerto Rico: a Sensitivity Analysis of an Agricultural Field Model and Riparian Zone Management Model

Agrichemical runoff from farmland may adversely impact coastal water quality. Two models, the Agricultural Policy/Environmental eXtender (APEX) and the Riparian Ecosystem Management Model (REMM), were used to evaluate the movement of the herbicide atrazine to the Jobos Bay National Estuarine Research Reserve from adjacent fields. The reserve is located on Puerto Rico’s southeast coast. Edge-of-field atrazine outputs simulated with the APEX were routed through a grass-forest buffer using the REMM. Atrazine DT₅₀ (half-life) values measured in both field and buffer soils indicated that accelerated degradation conditions had developed in the field soil due to repeated atrazine application. APEX simulations examined both the measured field and buffer soil atrazine DT₅₀ and the model’s default value. The use of the measured field soil atrazine degradation rate in the APEX resulted in 33 % lower atrazine transport from the field. REMM simulations indicated that the buffer system had the potential to reduce dissolved atrazine transport in surface runoff by 77 % during non-tropical storm events by increasing infiltration, slowing transport, and increasing time for pesticide degradation. During a large runoff event due to a tropical storm that occurred close to the time of an atrazine application, the REMM simulated only a 37 % reduction in atrazine transport. The results indicate that large storm events soon after herbicide application likely dominate herbicide transport to coastal waters in the region. These results agree with water quality measurements in the reserve. This study demonstrated the sensitivity of these models to variations in DT₅₀ values in evaluating atrazine fate and transport in the region and emphasizes that the use of measured DT₅₀ values can improve model accuracy.     

Measuring and modeling ammonium adsorption by calcareous soils

The aim of this study was assessment of ammonium (NH ₄ ⁺ ) adsorption isotherms in some agricultural calcareous soils and modeling of that by using the mechanistic exchange model. Ten surface soils (0–30 cm) were collected from areas covered with different land uses in Hamedan, western Iran. Isotherm experiments were carried out by concentrations of NH ₄ ⁺ prepared from NH₄Cl salt (0, 10, 20, 30, 40, 50, 100, and 150 mg NH ₄ ⁺ ?l^?1) in presence of 0.01 M CaCl₂ solution. The empirical models including simple adsorption isotherm and Freundlich equations were fitted well to the experimental data. The average amounts of adsorbed NH ₄ ⁺ in studied soils varied from 8.95 to 35.23 %. Adsorption percentage indicated positive correlation with pH, cation-exchange capacity (CEC), equivalent calcium carbonate, and clay content and had negative correlation with sand content. In order to predict and model NH ₄ ⁺ adsorption, cation-exchange model in PHREEQC program was used. The model could simulate the NH ₄ ⁺ adsorption very well in all studied soils. The values of CEC played the major role in modeling of NH ₄ ⁺ adsorption in this study indicating that cation-exchange process was the major mechanism controlling NH ₄ ⁺ adsorption in studied soils.     

Probability-based nitrate contamination map of groundwater in Kinmen

Groundwater supplies over 50 % of drinking water in Kinmen. Approximately 16.8 % of groundwater samples in Kinmen exceed the drinking water quality standard (DWQS) of NO₃ ^?-N (10 mg/L). The residents drinking high nitrate-polluted groundwater pose a potential risk to health. To formulate effective water quality management plan and assure a safe drinking water in Kinmen, the detailed spatial distribution of nitrate–N in groundwater is a prerequisite. The aim of this study is to develop an efficient scheme for evaluating spatial distribution of nitrate–N in residential well water using logistic regression (LR) model. A probability-based nitrate–N contamination map in Kinmen is constructed. The LR model predicted the binary occurrence probability of groundwater nitrate–N concentrations exceeding DWQS by simple measurement variables as independent variables, including sampling season, soil type, water table depth, pH, EC, DO, and Eh. The analyzed results reveal that three statistically significant explanatory variables, soil type, pH, and EC, are selected for the forward stepwise LR analysis. The total ratio of correct classification reaches 92.7 %. The highest probability of nitrate–N contamination map presents in the central zone, indicating that groundwater in the central zone should not be used for drinking purposes. Furthermore, a handy EC–pH-probability curve of nitrate–N exceeding the threshold of DWQS was developed. This curve can be used for preliminary screening of nitrate–N contamination in Kinmen groundwater. This study recommended that the local agency should implement the best management practice strategies to control nonpoint nitrogen sources and carry out a systematic monitoring of groundwater quality in residential wells of the high nitrate–N contamination zones.     

Numerical Modeling on the Effect of Dissolved Oxygen on Nitrogen Transformation and Transport in Unsaturated Porous System

Nitrogen pollution in groundwater resulting from wastewater application to land is a common problem, and it causes a major threat to groundwater-based drinking water supplies. In this study, a numerical model is developed to study the nitrogen species transport and transformation in unsaturated porous media. Further, a new mass transfer module for dissolved oxygen (DO) is incorporated in the one-dimensional numerical model for nitrogen species transport to describe the fate and transport of nitrogen species, dissolved oxygen, dissolved organic carbon (DOC), and biomass. The spatial and temporal variation of dissolved oxygen is incorporated in the model through the mass transfer from gaseous phase to water phase in an unsaturated porous system. The numerical results of the water flow model and single species and multispecies transport model in an unsaturated zone developed for this purpose have been validated with the available analytical/numerical solution. The developed model is applied in clay loam, silt, and sand soils to analyze the transport behavior of nitrogen species under unsaturated condition. The numerical results suggest that the high rate of oxygen mass transfer from the air phase to the water phase positively increases the dissolved oxygen in the applied wastewater and enhances the nitrification process. Because of this high oxygen mass transfer, the nitrate nitrogen concentration significantly increases in the unsaturated zone and the same is transported to a larger depth at higher simulation period. On the other hand, the low rate of oxygen mass transfer implicitly enhances the denitrification process and finally reduces the nitrate nitrogen concentration in the unsaturated zone. The numerical results also show that the nitrate nitrogen transport is rapid in sandy soil when compared with clay loam and silty soils under high oxygen mass transfer rate. In essence, the high oxygen mass transfer rate significantly increases the nitrate nitrogen in the unsaturated zone, especially at a greater depth at larger time levels and eventually affects the groundwater quality.     

Nitrogen mineralization in O horizon soils during 27 years of nitrogen enrichment at the Bear Brook Watershed in Maine,USA

Chronic elevated nitrogen (N) deposition has altered the N status of temperate forests, with significant implications for ecosystem function. The Bear Brook Watershed in Maine (BBWM) is a whole paired watershed manipulation experiment established to study the effects of N and sulfur (S) deposition on ecosystem function. N was added bimonthly as (NH₄)₂SO₄ to one watershed from 1989 to 2016, and research at the site has studied the evolution of ecosystem response to the treatment through time. Here, we synthesize results from 27 years of research at the site and describe the temporal trend of N availability and N mineralization at BBWM in response to chronic N deposition. Our findings suggest that there was a delayed response in soil N dynamics, since labile soil N concentrations did not show increases in the treated watershed (West Bear, WB) compared to the reference watershed (East Bear, EB) until after the first 4 years of treatment. Labile N became increasingly available in WB through time, and after 25 years of manipulations, treated soils had 10× more extractable ammonium than EB soils. The WB soils had 200× more extractable nitrate than EB soils, driven by both, high nitrate concentrations in WB and low nitrate concentrations in EB. Nitrification rates increased in WB soils and accounted for ~?50% of net N mineralization, compared to ~?5% in EB soils. The study provides evidence of the decadal evolution in soil function at BBWM and illustrates the importance of long-term data to capture ecosystem response to chronic disturbance.     

The Impacts of Different Meteorology Data Sets on Nitrogen Fate and Transport in the SWAT Watershed Model

In this study, we investigated how different meteorology data sets impacts nitrogen fate and transport responses in the Soil and Water Assessment Tool (SWAT) model. We used two meteorology data sets: National Climatic Data Center (observed) and Mesoscale Model 5/Weather Research and Forecasting (simulated). The SWAT model was applied to two 10-digit hydrologic unit code watersheds in the Coastal Plain and Piedmont zones of North Carolina. Nitrogen cycling and loading response to these meteorological data were investigated by exploring 19 SWAT nitrogen outputs relating to landscape delivery, biogeochemical assimilation, and atmospheric deposition. The largest difference in model output using both meteorology data sets was for large loads/fluxes. Landscape delivery outputs (e.g., NO^? ₃ watershed discharge, groundwater NO^? ₃ flux, soil NO^? ₃ percolation) showed the largest difference across all values. Use of the two weather data sources resulted in a nearly twofold difference in NO^? ₃ watershed discharge and groundwater NO^? ₃ flux. Differences for many nitrogen outputs were greater than those for sub-basin flow. Nitrogen outputs showed the greatest difference for agricultural land covers and there was no flow-related pattern in output differences across sub-basins or over time (years). In general, nitrogen parameter models that had a greater number of nitrate concentration, flow, and temperature terms (equation variables) in each transport model showed the greatest difference between both meteorology applications.     

Implications of climate change for evaporation from bare soils in a Mediterranean environment

The purpose of this study was to predict quantitative changes in evaporation from bare soils in the Mediterranean climate region of Turkey in response to the projections of a regional climate model developed in Japan (hereafter RCM). Daily RCM data for the estimation of reference evapotranspiration (ET _r) and soil evaporation were obtained for the periods of 1994–2003 and 2070–2079. Potential evaporation (E _p) from bare soils was calculated using the Penman–Monteith equation with a surface resistance of zero. Simulation of actual soil evaporation (E _a) was carried out using Aydin model (Aydin et al., Ecological Modelling 182:91–105, 2005) combined with Aydin and Uygur (2006, A model for estimating soil water potential of bare fields. In Proceedings of the 18th International Soil Meeting (ISM) on Soils Sustaining Life on Earth, Managing Soil and Technology, Sanliurfa, 477–480pp.) model of predicting soil water potential at the top surface layer of a bare soil, after performances of Aydin model (R ² = 94.0%) and Aydin and Uygur model (R ² = 97.6) were tested. The latter model is based on the relations among potential soil evaporation, hydraulic diffusivity, and soil wetness, with some simplified assumptions. Input parameters of the model are simple and easily obtainable such as climatic parameters used to compute the potential soil evaporation, average diffusivity for the drying soil, and volumetric water content at field capacity. The combination of Aydin and Aydin and Uygur models appeared to be useful in estimating water potential of soils and E _a from bare soils, with only a few parameters. Unlike ET _r and E _p projected to increase by 92 and 69 mm (equivalent to 8.0 and 7.3% increases) due to the elevated evaporative demand of the atmosphere, respectively, E _a from bare soils is projected to reduce by 50 mm (equivalent to a 16.5% decrease) in response to a decrease in rainfall by 46% in the Mediterranean region of Turkey by the 2070s predicted by RCM, and consequently, to decreased soil wetness in the future.     

Uptake of arsenic by alkaline soils near alkaline coal fly ash disposal facilities

The attenuation of arsenic in groundwater near alkaline coal fly ash disposal facilities was evaluated by determining the uptake of arsenic from ash leachates by surrounding alkaline soils. Ten different alkaline soils near a retired coal fly ash impoundment were used in this study with pH ranging from 7.6 to 9.0, while representative coal fly ash samples from two different locations in the coal fly ash impoundment were used to produce two alkaline ash leachates with pH 7.4 and 8.2. The arsenic found in the ash leachates was present as arsenate [As(V)]. Adsorption isotherm experiments were carried out to determine the adsorption parameters required for predicting the uptake of arsenic from the ash leachates. For all soils and leachates, the adsorption of arsenic followed the Langmuir and Freundlich equations, indicative of the favorable adsorption of arsenic from leachates onto all soils. The uptake of arsenic was evaluated as a function of ash leachate characteristics and the soil components. The uptake of arsenic from alkaline ash leachates, which occurred mainly as calcium hydrogen arsenate, increased with increasing clay fraction of soil and with increasing soil organic matter of the alkaline soils. Appreciable uptake of arsenic from alkaline ash leachates with different pH and arsenic concentration was observed for the alkaline soils, thus attenuating the contamination of groundwater downstream of the retired coal fly ash impoundment.     

Adsorption-desorption and leaching potential of glyphosate and aminomethylphosphonic acid in acidic Malaysian soil amended with cow dung and rice husk ash

This study investigates adsorption-desorption and the leaching potential of glyphosate and aminomethylphosphonic acid (AMPA) in control and amended—addition of cow dung or rice husk ash—acidic Malaysian soil with high oxide mineral content. The addition of cow dung or rice husk ash increased the adsorptive removal of AMPA. The isotherm data of glyphosate and AMPA best fitted the Freundlich model. The constant K_f for glyphosate was high in the control soil (544.873 mg g^?1) followed by soil with cow dung (482.451 mg g^?1) then soil with rice husk ash (418.539 mg g^?1). However, for AMPA, soil with cow dung was high (166.636 mg g^?1) followed by soil with rice husk ash (137.570 mg g^?1) then the control soil (48.446 mg g^?1). The 1/n values for both glyphosate and AMPA adsorptions were <?1 indicating their strong affinity for adsorbents. Desorption of both glyphosate and AMPA occurred only in the control soil. The compounds were not detected in soils with added cow dung or rice husk ash. The addition of cow dung or rice husk ash increased glyphosate mobility. However, ground water ubiquity scores for both control and amended soils were <?2.8. This indicated glyphosate is a transitional herbicide; therefore, its leaching potential in the soil is low, despite the addition of cow dung or rice husk ash. Addition of these wastes decreased the mobility and leaching potential of AMPA. The addition of cow dung or rice husk ash could be beneficial in increasing adsorption and enhancing degradation of these compounds.     

Spatial and Seasonal Variations in Chlorophyll-Nutrient Relationships in The Shallow Hypertrophic Lake Manyas, Turkey

Regression and correlation analyses were used to predict responses of phytoplankton biomass (chlorophyll) (μg L⁻¹) to nitrate (NO₃) (mg L⁻¹), phosphate (PO₄) (mg L⁻¹) and ammonium (NH₄) (mg L⁻¹) dynamics in the shallow hypertrophic Lake Manyas, Turkey. Nutrient concentrations showed a descending gradient with distance, while chlorophyll concentrations showed an ascending gradient with the distance from the Sığırcı Inlet to the Karadere Outlet. Higher nutrient concentrations did always not coincide with higher chlorophyll concentrations. The results showed that regression models developed using seasonal data were more accurate in predicting chlorophyll concentrations than those developed using the pooled data from whole year (based on R ² and the difference between the measured and predicted values). The findings also revealed that within a single large shallow lake, chlorophyll-nutrient relationships might show significant variations spatially. The objective of this study was to determine the seasonal and spatial variations in the relationships between chlorophyll, nitrate, phosphate and ammonium in the shallow hypertrophic Lake Manyas, Turkey.     

Transport of NO<Subscript>X</Subscript> Emissions from Sugarcane Fertilisation into the Great Barrier Reef Lagoon

The Great Barrier Reef World Heritage Area contains highly sensitive ecosystems that are threatened by the effects of anthropogenic activity including eutrophication. The nearby sugarcane plantations of tropical north Queensland are fertilised annually and there has been ongoing concern about the magnitude of the loss of applied nitrogen to the environment. Previous studies have considered the potential of rainwater run-off to deposit reactive nitrogen species into rivers and ultimately into the Great Barrier Reef Lagoon, but have neglected the possibility of transport via the atmosphere. This paper reports the results of a modelling study commissioned by Australia’s National Heritage Trust aimed at assessing whether or not atmospheric deposition of reactive nitrogen from Queensland’s sugarcane plantations posed a potential threat to the Great Barrier Reef Lagoon. Atmospheric dispersion modelling was undertaken using The Air Pollution Model, developed by Australia’s Commonwealth Scientific and Industrial Research Organisation. Despite the predominance of onshore southeasterly winds, the dispersion model results indicate that 9% of the time during the sugarcane fertilization season (in the modeled years 2001–2006) the meteorological conditions resulted in emissions from the coastal regions of north Queensland being transported out over the ocean around the Great Barrier Reef. The results suggest that there may be a greater efficiency for transport out over the reef during October than for November and December. For the 2 months that exhibited the greatest potential for transport of coastal pollution to the Great Barrier Reef, the modeled deposition of nitrogen oxides (NO_X) into the Great Barrier Reef lagoon was less than 1% of the total emissions from the sugarcane plantations, but was not zero. Our model has a simple chemical scheme that does not cover the full chemistry of all reactive nitrogen compounds and so the results are only indicative of the potential levels of deposition. Nevertheless, our study shows that small amounts of NO_X that originate from sugarcane fertilization may be transported and dry deposited into the Great Barrier Reef lagoon. Other pathways not included in the modeling scheme may provide a more efficient transport mechanism. Whilst modern practices for the application of fertilizer to sugarcane plantations have drastically reduced emissions, the potential efficiency of transport of pollutants via the atmosphere may be of concern for other more highly polluting agricultural industries.     

Assessing uncertainty in annual nitrogen,phosphorus, and suspended sediment load estimates in three agricultural streams using a 21-year dataset

Accurate estimation of constituent loads is important for studies of ecosystem mass balance or total maximum daily loads. In response, there has been an effort to develop methods to increase both accuracy and precision of constituent load estimates. The relationship between constituent concentration and stream discharge is often complicated, potentially leading to high uncertainty in load estimates for certain constituents, especially at longer-term (annual) scales. We used the loadflex R package to compare uncertainty in annual load estimates from concentration vs. discharge relationships in constituents of interest in agricultural systems, including ammonium as nitrogen (NH₄-N), nitrate as nitrogen (NO₃-N), soluble reactive phosphorus (SRP), and suspended sediments (SS). We predicted that uncertainty would be greatest in NO₃-N and SS due to complex relationships between constituent concentration and discharge. We also predicted lower uncertainty with a composite method compared to regression or interpolation methods. Contrary to predictions, we observed the lowest uncertainty in annual NO₃-N load estimates (relative error 1.5–23%); however, uncertainty was greatest in SS load estimates, consistent with predictions (relative error 19–96%). For all constituents, we also generally observed reductions in uncertainty by up to 34% using the composite method compared to regression and interpolation approaches, as predicted. These results highlight differences in uncertainty among different constituents and will aid in model selection for future studies requiring accurate and precise estimates of constituent load.     

A Novel Technique for Monitoring Contaminant Transport Through Soils

An experimental setup which is capable of simulation as well as monitoring contaminant transport through soil mass has been developed. Efficiency of the experimental setup has been demonstrated by comparing the results obtained from the argentometric method. Reynolds number (R_e) and the Peclet number (P_e) have been found to be less than unity. This indicates that flow of the solute through soil mass is laminar and the dominant contaminant transport mechanism is diffusion.     

A physicochemical model of sorption processes in NO2 passive sampling with air humidity effects

Aqueous triethanolamine (TEA) solutions are widely used as sorption medium for passive sampling of ambient NO₂, with NO₂ trapped and accumulated as nitrite ion. The results of test measurements of ambient NO₂ concentrations using passive sampling method showed that the simple approach commonly used to describe passive sampling process might lead to substantial systematic errors. Presented in the article is a new physicochemical model of the process of passive sampling of gaseous NO₂, with aqueous TEA solution used as a trapping medium. The model is based on the available results of experimental studies of interaction of gaseous NO₂ with TEA/water solutions. The key principles underlying the model are: (1) when absorbed by a trapping solution, NO₂ forms nitrite ion only on the condition that TEA is hydrated; (2) coefficient of conversion of NO₂ to NO ₂ ^? is equal to one when reacting with hydrated TEA; and (3) the fraction of hydrated TEA molecules depends on air humidity at the moment of measurement. Validation of the model was made using the data of the field measurements carried out in the Middle Urals in 2007–2009. The new model was used to calculate average NO₂ concentrations. Concentrations calculated agreed well with the results obtained by reference methods. The difference between the datasets was statistically insignificant.     

上海市秋季典型PM2.5污染过程数值预报分析

基于2012年10月上海出现的一次典型PM2.5污染案例,验证评估上海市空气质量数值预报系统Model-3/CMAQ预报性能,采用过程分析技术,定量评估不同大气物理化学过程对上海代表性点位PM2.5浓度变化的作用规律。结果表明:Model-3/CMAQ模式系统能较好地反映PM2.5的浓度变化趋势与特点。对于上海市区点位(徐汇上师大)和东南部点位(奉贤海湾和浦东惠南),PM2.5浓度上升主要受本地源排放影响,其贡献比例超过40%,其次是区域大气传输作用的影响。对于西北部点位(崇明监测站和青浦淀山湖),区域大气传输是PM2.5浓度上升的主要原因,贡献比例超过70%,其次是源排放。各点位PM2.5浓度的主要去除途径均为大气传输,贡献比例均超过70%,其次是干沉降。气溶胶过程对PM2.5主要起二次颗粒物生成的作用,特别是市区及东南部点位,贡献比例较西北部点位更高。     

MODELING NITROGEN UPTAKE AND POTENTIAL NITRATE LEACHING UNDER DIFFERENT IRRIGATION PROGRAMS IN NITROGEN-FERTILIZED TOMATO USING THE COMPUTER PROGRAM NLEAP

Readily available nitrogen (N) sources such as ammonium nitrate with excessive irrigation present a potential hazard for the environment. The computer program Nitrate Leaching and Economic Analysis Package (NLEAP) is a mechanistic model developed for rapid site-specific estimates of nitrate-nitrogen (NO₃-N) moving below the root zone in agricultural crops and potantial impacts of NO₃-N leaching into groundwater. In this study, the value of NLEAP was tested to simulate N uptake by crops and NO₃-N leaching parameters in large lysimeters under the tomato crop. Three seedlings of tomato variety of H-2274 (Lycopersicum esculentum L.) were transplanted into each lysimeter. N fertilizer at the rate of 140 kg N ha^–1 was sidedressed in two split applications, the first half as ammonium sulphate and the second half as ammonium nitrate. The lysimeters were irrigated based on programs of C 0.75, 1.00, 1.25 and 1.50, C referring to class A-Pan evaporation coefficients. Parameters such as leaching index (LI), annual leaching risk potential (ALRP), N available for leaching (NAL), amount of NO₃-N leached (NL) and amount of N taken up by the crops (NU) were estimated using the NLEAP computer model. To test the ability of model to simulate N uptake and NL, measured values were compared with simulated values. Significant correlations, R² = 0.92 and P < 0.03 for the first year and R² = 0.86 and P < 0.06 for the second year, were found between measured and simulated values for crop N consumption, indicating that the NLEAP model adequately described crop N uptake under the varied irrigation programs using an optimal N fertilization program for the experimental site. Significant correlations, R² = 0.96 and P < 0.01 for the first year and R² = 0.97 and P < 0.01 for the second year, were also found between measured and simulated values of NL, indicating that the NLEAP model also adequately predicted NL under the varied irrigation programs. Therefore, this computer model can be useful to estimate the NO₃-N moving beyond the root zone under conditions in which the present experiment was carried out. Also, the NLEAP-estimated NAL values and other parameters can also be used to improve N management practices and N fertilizer recommendations that will help to decrease the adverse effect of N fertilizer on groundwater quality and farm profitability.