Journal of Environmental Quality 31:671-681 (2002)
© 2002 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
TECHNICAL REPORT
Vadose Zone Processes and Chemical Transport
Fate of Nitrate and Bromide in an Unsaturated Zone of a Sandy Soil under Citrus Production
S. Paramasivam*,a,
A. K. Alvab,
A. Faresc and
K. S. Sajwand
a Center for Marine, Environmental Sciences, and Biotechnology Research, P.O. Box 20600, Savannah State University, Savannah, GA 30404
b USDA-ARS, Vegetable and Forage Crops Research Unit, 24106 N. Bunn Road, Prosser, WA 99350
c University of Florida, Institute of Food and Agricultural Sciences, Citrus Research and Education Center, 700 Experiment Station Road, Lake Alfred, FL 33850
d CMESBR, Savannah State University Savannah, GA 30404
* Corresponding author (siva{at}savstate.edu)
Received for publication October 10, 2000.
 |
ABSTRACT
|
|---|
Understanding water and nutrient transport through the soil profile is important for efficient irrigation and nutrient management to minimize excess nutrient leaching below the rootzone. We applied four rates of N (28, 56, 84, and 112 kg N ha-1; equivalent to one-fourth of annual N rates being evaluated in this study for bearing citrus trees), and 80 kg Br- ha-1 to a sandy Entisol with >25-yr-old citrus trees to (i) determine the temporal changes in NO3N and Br- distribution down the soil profile (2.4 m), and (ii) evaluate the measured concentrations of NO3N and Br- at various depths with those predicted by the Leaching Estimation and Chemistry Model (LEACHM). Nitrate N and Br- concentrations approached the background levels by 42 and 214 d, respectively. Model-predicted volumetric water content and concentrations of NO3N and Br- at various depths within the entire soil profile were very close to measured values. The LEACHM data showed that 21 to 36% of applied fertilizer N leached below the root zone, while tree uptake accounted for 40 to 53%. Results of this study enhance our understanding of N dynamics in these sandy soils, and provide better evaluation of N and irrigation management to improve uptake efficiency, reduce N losses, and minimize the risk of ground water nitrate contamination from soils highly vulnerable to nutrient leaching.
Abbreviations: ETP, potential evapotranspiration LEACHM, Leaching Estimation and Chemistry Model LEACHN, one-dimensional water flow and nitrogen transport model
 |
INTRODUCTION
|
|---|
WATER FLOW IN the unsaturated zone of agricultural soils is of great importance in determining the fate of surface-applied soluble nutrients and pesticides (Nightingale, 1972; Hallberg, 1989; Hubbard et al., 1986, 1989). Nitrate (NO-3) leaching associated with routine use of nitrogen (N) fertilizers and the potential for subsequent movement below the root zone and ground water contamination is influenced by soil properties, water management, and climatic factors.
In major citrus growing regions in central Florida, soils are typically Entisols, without confining soil horizons. A combination of homogeneous sandy soils (>95% sand), shallow citrus root systems, high (approximately 1300 mm per year) and unevenly distributed rainfall (two-thirds of annual rainfall during May through September) creates favorable conditions for leaching of surface-applied soluble nutrients and agrichemicals with drainage water below the rootzone. Therefore, efficient nutrient and irrigation management are critical to optimize nutrient uptake efficiency and to minimize losses of surface-applied fertilizer nutrients. Maintaining nutrients within the effective rootzone is essential for the crop to get maximum benefit from the applied fertilizer.
Previous citrus N management studies were usually conducted to evaluate the effects of rates and sources of N on the tree nutritional status, yield, and fruit quality, with very little attention to the fate of N below the rootzone and potential nonpoint-source NO3N contamination of ground water (Davies, 1996). Mansell et al. (1986) studied the temporal changes in distributions of NO3N and NH4N concentrations in a profile of Oldsmar fine sand (sandy siliceous, hyperthermic Alfic Arenic Haplaquod) from a citrus grove with 16-yr-old trees near Fort Pierce, Florida. This study was conducted under three different soil tillage practices (i.e., shallow-tilled [015 cm], deep-tilled [0105 cm], and deep-tilled [0105 cm]) with 56 Mg ha-1 dolomitic limestone. They applied 15N-labeled (NH4)2SO4 at a rate of 115 kg N ha-1 and sampled the soil 12, 42, 75, and 134 d after the fertilizer application. Results of their study confirmed earlier observations (Mansell et al., 1977) that deep tillage plus lime incorporation into the profile of a Spodosol enhanced citrus root absorption of fertilizer N, minimized N leaching losses, and thus minimized ground water NO3N contamination. Their study also showed that to minimize leaching, N fertilizer application should not be immediately followed by excessive irrigation. The Spodosol used in this study had a slightly heavier texture compared with Entisols and an impermeable organic layer within 100 cm from the soil surface also not found in Entisols. Therefore, the results of the Mansell et al. (1986) study are not applicable to deep fine sandy Entisols. The current study was designed to provide data to increase understanding of the movement and distribution of fertilizer nitrogen in an Entisol under mature citrus production.
In recent years, solute transport models, such as the Leaching Estimation and Chemistry Model (LEACHM; Wagenet and Hutson, 1989), have been used to evaluate the fate of contaminants in agricultural soils as basis to develop best management practices (Pennell et al., 1990; Soulsby and Reynolds, 1992; Jabro et al., 1993). Jabro et al. (1993) compared the transport of NO-3 measured in a Hagerstown silt loam (fine, mixed, semiactive, mesic Typic Hapludalf) soil with that predicted by LEACHM (Version 3.0) (Hutson and Wagenet, 1992) and NCSWAP (Nitrogen, Carbon, Soil, Water And Plant) (Molina and Richards, 1984) models and concluded that these models did not accurately predict NO3N leaching to a depth of 1.2 m. They concluded that NCSWAP resulted in significant overestimation of the amount of NO-3 leached from the soil profile due to an inability to simulate macropore flow effects. It was speculated that NO-3 was stored in the micropore system and that percolating water bypassed this NO-3. However, LEACHM had the same overestimation problem as a result of its inability to estimate macropore flow effects only for the soil profile of the nontreated plot. Jabro et al. (1993) further cited the use of an inappropriate water retention function fitted by Campbell's equation (Campbell, 1974) as a reason for an overestimation of soil water content, which subsequently resulted in nonsignificant differences between LEACHM simulation and measured data of NO-3 leaching in the profile. On the basis of overall performance of these models, Jabro et al. (1993) concluded that NO-3 leaching was significantly better simulated by the LEACHM than by the NCSWAP model. Our parallel study (Paramasivam et al., 2000a) showed a good agreement between the measured concentrations of NH4N and NO3N following a heavy loading of NH4NO3 liquid to an uncropped Entisol (Candler fine sand) through 270 cm depth and the respective concentrations simulated by LEACHM. Soulsby and Reynolds (1992) used LEACHM to model the soil water flux in an aluminum (Al) leaching study. They calibrated the model using in situ tensiometer data, then compared model predictions against measured tensiometer data for the remainder of the year and found good agreement.
Bromide has been widely used as a tracer to monitor and/or predict water and nutrient transport in soils (Smith and Davis, 1974; Bowman, 1984; Owens et al., 1985) because it is free from chemical or biological transformations in the soil (Bowman, 1984; Owens et al., 1985). Several field studies with bromide tracer were conducted in central Florida citrus groves to determine intrinsic spatial and temporal variability, travel time of bromide through the vadose zone, and lateral movement of bromide in the surficial aquifer (Hornsby et al., 1990; Foussereau and Graham, 1997; Paramasivam et al., 1999). Comfort et al. (1993) reported a good agreement between the measured Br- concentrations in the soil profile to a 1-m depth of an uncropped Brocko silt loam (coarse-silty, mixed, Borollic Calciorthid) soil with application of 56 kg ha-1 Br- and those predicted by LEACHM.
The objectives of this study were to (i) determine temporal changes in distribution of NO3N following a broadcast application of water-soluble dry granular fertilizer mixture containing NH4NO3 as the N source, (ii) evaluate the effect of water input (irrigation and rainfall) on transport and distribution of NO3N, (iii) compare the transport and distribution of NO3N and a nonreactive tracer (Br-) applied with the fertilizer, and (iv) evaluate LEACHM predictions of NO3N and Br- transport in the soil profile under standard citrus management practices and compare them with measured concentrations.
 |
MATERIALS AND METHODS
|
|---|
A field study was conducted in a highly productive commercial citrus grove with 25-yr-old Hamlin orange trees [Citrus sinensis (L.) Osbeck] on Cleopatra mandarin (Citrus reticulata Blanco) rootstock planted (7.62 x 6.52 m) in a moderately well-drained Tavares fine sand (hyperthermic, uncoated, Typic Quartzipsamment). Physicochemical properties of the soil at various depths of the experimental site are given in Table 1. This experiment was conducted as part of an ongoing long-term experiment to evaluate the effect of various N rates and sources of N on fruit yield and quality responses in a vulnerable soil under optimal irrigation conditions. The tree response in terms of leaf mineral concentrations, fruit yield, and quality was reported earlier (Alva and Paramasivam, 1998) and a brief description is given here. Each plot covering an area of 198 m2 consisted of four uniform trees, with the middle two trees being used for fruit yield and quality evaluation N treatments. Trees were irrigated with under-tree microsprinklers (one emitter per tree, with a delivery rate of 0.083 m3 h-1 at 0.276 MPa, which wetted an area of 18.7 m2, approximately equivalent to the ground area under the canopy). Irrigation scheduling was based on matric potential values recorded by tensiometers at the 15- and 30-cm depths along the tree drilling. Irrigation was scheduled when the matric potential at the 15-cm depth attained -10 kPa for January through June and -15 kPa for July through December (Smajstrala et al., 1987; Parsons, 1989). Thus, irrigation water was applied to replenish soil water in the depth of maximum root activity (i.e., 90 cm) to field capacity. Total water delivery at each irrigation was recorded using a flow meter. Rainfall was recorded using an automatic rain gauge.
Experimental plots receiving 0, 112, 224, 336, and 448 kg N ha-1 yr-1 as water-soluble dry granular fertilizer applied in four equal split doses (February, April, May, and October) were used to monitor the NO3N and nonreactive tracer (Br-) movement. A water-soluble granular blend containing N, P, K, Mg, and S at a ratio of 13:2.2:13.5:3.25:6.5 derived from ammonium nitrate, triple super phosphate, potassium chloride (muriate of potash), and potassium magnesium sulfate (sul-po-mag) was broadcasted using commercial broadcast equipment calibrated for single-sided application (Conibear Equipment Co., Lakeland, FL). Since this monitoring started soon after the first split fertilizer application in February, the N application rates for the purpose of this study were 0, 28, 56, 84, and 112 kg N ha-1. The February application was nearly four months after the previous years' last fertilizer application. Therefore, the residual fertilizer N effects of the previous year were negligible.
Bromide was applied at a rate of 80 kg ha-1 (as KBr) on the treated area based on the same day as the first dose of fertilizer. The KBr was applied in solution using a herbicide boom with a band width of 2.4 m on either side of the tree trunk for four trees in each plot. The area of Br- application per plot was 30.5 x 4.8 m2. Since this part of the study involved intensive soil sampling, bromide was applied only to three replicated plots that received water-soluble granular fertilizer out of five replicates of the original randomized complete block experiment. This experiment was conducted as a randomized complete block design, with three replications.
Field Sampling and Soil Analysis
Soil samples were taken 90 cm inside the drip line under the canopy using a 2.5-cm-diameter bucket auger at 15- and 30-cm depths, and thereafter at 30-cm increments to the soilground water interface. The average depth to ground water was about 300 cm. Triplicate cores of soil samples per replicate per treatment were taken from the 15 (five treatments x three repititions) plots before fertilizer and bromide application. Soil sampling continued at weekly intervals after the fertilizer application. Most of the soil sampling was done within 24 to 48 h of irrigation or rainfall events. Soil cores were taken within the wetted zone of fertilizer and Br- application under the tree canopy. Soil samples from each depth (cores and replicates) were stored separately in an ice chest for transportation to the laboratory and were extracted with water at 1:1 (soil to water) on a field-moist basis. At each sampling, nine samples per depth per treatment were analyzed individually for concentrations of NO3N and Br- and for water content.
Twenty grams of field moist soil were weighed into a 50-mL centrifuge tube and 20 mL deionized water were added. The solution was shaken for 30 min, centrifuged at 3270 x g for 10 min, and filtered through Whatman (Maidstone, UK) 42 filter paper. The concentrations of NO3N and Br- were determined within 24 h of extraction using an ion chromatograph (Dionex 300 [Sunnyvale, CA]) using the procedure outlined by the USEPA (1991). Gravimetric soil moisture was determined in a subsample for each depth to calculate the soil moisture percentage on an oven-dry soil weight basis. The concentrations of NO3N and Br- in the soil were expressed as mg kg-1 dry soil basis. Soil sampling continued until the Br- concentrations in the soil profile approached nondetectable concentrations. This simulated the time taken for the complete disappearance of Br- from this soil profile.
Leaching Estimation and Chemistry Model Simulations
The Leaching Estimation and Chemistry Model (LEACHM) is a one-dimensional water flow, solute transport, and plant uptake model. It is composed of different submodels including LEACHN, a one-dimensional water flow and N transport model (Hutson and Wagenet, 1992). The choice of LEACHM was based on our experience (Paramasivam et al., 2000a) and that of other researchers (Pennel et al., 1990; Lamb, 1996; Graham and Wheaton, 1999; Lamb et al., 1999) under Florida sandy conditions in nearby locations. The current modeling exercise benefitted substantially from these previous works with LEACHM; thus, there was no need for an extensive calibration exercise since most of the input parameters used were taken from these simulation works specific to our current conditions.
The input parameters for LEACHN include: (i) daily water input (rainfall or irrigation) (Fig. 1)
, (ii) weekly potential evapotranspiration (ETP) and minimum and maximum temperatures (Fig. 2)
, and (iii) soil physical properties including water retention curve parameters and hydraulic conductivity for each soil profile (Tables 1 and 2). Soil moisture content and concentrations of NO3N and Br- for the soil samples taken prior to fertilizer and Br- application were the initial conditions for the model simulation. Nitrogen transformation rates and partitioning coefficients (Table 3) were the input parameters.

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 1. Water balance components including rainfall, irrigation, simulated drainage below the rootzone, and evapotranspiration for a Tavares fine sand, with 25-yr-old Hamlin orange trees on Cleopatra mandarin rootstock. Arrows along the x axis show the soil sampling times during the study period.
|
|

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 2. Weekly minimum and maximum air temperature and total potential evapotranspiration used as input parameters for the Leaching Estimation and Chemistry Model (LEACHM) model.
|
|
View this table:
[in this window]
[in a new window]
|
Table 2. Best-fit retention parameters for Campbell's water retention equations for the different soil layers as determined by RETFIT and used by the Leaching Estimation and Chemistry Model (LEACHM).
|
|
View this table:
[in this window]
[in a new window]
|
Table 3. Nitrogen transformation rates and transport parameters as input for the Leaching Estimation and Chemistry Model (LEACHM).
|
|
Soil Physical Properties
At the beginning of the study, undisturbed soil cores were taken at 0- to 0.15-, 0.15- to 0.30-, 0.30- to 0.60-, 0.60- to 0.90-, 0.90- to 1.20-, and 1.20- to 1.50-m depths, at five locations from the experimental site. These soil cores were used to determine soil water release curves (Klute, 1986) and saturated hydraulic conductivity (Klute and Dirksen, 1986) at each of the five depths. Details of the soil water release curve were reported earlier by Paramasivam et al. (2000b). Subsamples of these cores were used to determine selected additional soil properties and reported by Paramasivam et al. (2001). These values were in agreement with those values reported by Carlisle et al. (1989) for the profile samples collected from the same site. Data from the deepest-measured soil cores were used for the remainder of the simulated soil profile. A stand-alone program, RETFIT, which accompanies the LEACHN model, was used to determine best-fit retention parameters for Campbell's (1974) water retention equation from these data (Table 2).
Crop Factor
Citrus root distribution in each soil layer was estimated from the literature. Based on Bowman (1984), as well as studies conducted by Castle (1980) and Alva and Syvertsen (1991), about 90% of the roots were found in the top 0.3 m and the remaining 10% of the roots were evenly distributed over the rest of the simulated profile. A constant mature citrus crop was assumed to be present throughout the simulation period, and with a crop cover fraction of 1 (fraction of land covered by the crop). Additional soil- and crop-related parameters used in the simulation are summarized in Table 3.
Weather Data
Daily potential evapotranspiration, ETP, was calculated for the study period using daily minimum and maximum air temperatures (Fig. 2) recorded at the nearby Archbold Biological Station in central Florida as input to the potential evapotranspiration calculation model, ETM, (Fares, 1996). The ETM uses the PriestleyTaylor equation (Priestley and Taylor, 1972) to calculate ETP. Input data for ETM include: daily minimum and maximum temperatures, altitude and latitude of the location, albedo, and sunshine ratio. More detailed information can be found about ETM in Fares (1996) and Fares and Mansell (1996). Daily ETP data were summed to calculate a weekly total ETP as required by LEACHN.
The model assumes no surface runoff. This is justified for the study site due to its flat terrain and its high saturated hydraulic conductivity (Ks). The LEACHN model simulated a transient or a steady state water flow. The model solves the one-dimensional Richards equation for transient water flow using a block-centered finite-difference approximation technique. The upper boundary condition of the simulated system is a flux-controlled surface boundary during the period of evaporation, nonponded infiltration, and zero flux. However, during periods of ponded infiltration, the pressure potential of the first node is set to zero. Freely draining profile boundary conditions (a unit hydraulic potential gradient = 1) were used in the lower boundary condition of the simulated system given the absence of a water table at the simulated depth. The simulated soil profile was equally divided into 16 horizontal segments, each 0.15 m thick. The predicted concentrations of NO-3 and Br- and volumetric water content were reported for the center of each depth segment (0.075, 0.225, up to 2.35 m). A value of 20 mm was used for the soil dispersivity for all layers, as recommended by Hutson and Wagenet (1992). More information about model calibration and sensitivity to Florida sandy soil under similar weather conditions for citrus crop was previously given by several researchers (Pennel et al., 1990; Lamb, 1996; Graham and Wheaton, 1999; Lamb et al., 1999; Paramasivam et al., 2000a).
 |
RESULTS AND DISCUSSION
|
|---|
Water Balance Components
Water balance components, including measured rainfall and irrigation, simulated drainage below the rootzone, and calculated evapotranspiration values for experimental plots, are shown in Fig. 1 and Table 2 for an entire growing season. During the first 21 d of this monitoring period, the experimental area received about 45 and 182 mm of cumulative irrigation and rainfall, respectively. For the same period, calculated evapotranspiration loss was 47 mm. The LEACHM simulation predicted no water drainage losses below the 2.4-m depth profile during this period. Drainage losses were predicted to start 28 d after fertilizer and Br- application and continued through the study period (Fig. 1). During the first 42 d of intensive monitoring of nutrient movement, cumulative irrigation and rainfall accounted for 95 and 213 mm, respectively. During the same period, calculated evapotranspiration and simulated drainage loss below the depth of 2.4 m accounted for 95 and 105 mm, respectively (Fig. 1). Water mass balance for the entire crop year is presented in Table 4.
Amount of available nutrients in the soil profile, amount of available water for leaching, and soil profile characteristics are important factors to determine the nutrient leaching losses from the applied fertilizers (Staver and Brinsfield, 1990). The amount of precipitation and/or irrigation in excess of evapotranspiration along with soil profile characteristics determine the amount of drainage water below the rooting depth. These in turn control the movement of wetting front in the soil profile and thereby determine the water content and nutrient distribution at various depths in a soil profile at each sampling time. Nutrient contents measured at various depths in a soil profile at a given time represent only brief glimpses of nutrient movement in the soil profile with time and do not attempt to represent the fate of the total nutrient applied to soil in the form of a fertilizer mixture. Campbell (1978) reported that sandy soils, along with high rainfall, create conditions conducive to rapid nutrient movement.
Measured and Predicted Soil Water Contents
Mean and standard error of measured soil water contents of samples collected at various depths in the soil profile at various sampling dates are shown with filled circle symbols in Fig. 3
. Corresponding simulated water contents are shown in the same figures as continuous lines. Generally, LEACHM-predicted water content reasonably agreed with measured soil water content in the soil profile to a depth of 0.9 to 1.2 m (Fig. 3) in all treatments. Under cropped conditions the depth of root distribution was an important input for LEACHM to predict plant water uptake and subsequent volumetric water content as a function of soil depth. Since the majority of roots are in the top 90-cm depth (Alva and Syvertsen, 1991; Zhang et al., 1996), the effect of root uptake on soil water and plant nutrients would be substantial between 0 and 90 cm depth compared with the rest of the soil profile. The good agreement between predicted and measured volumetric water contents in soils indicates that the input of 90 cm as rooting depth to the LEACHM root functions accurately predicted water uptake by citrus trees.

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 3. Distribution of measured and Leaching Estimation and Chemistry Model (LEACHM)-simulated volumetric water content in the soil profile of a Tavares fine sand at various sampling times following an irrigation and rainfall events. Measured values are the mean and standard error of 45 samples per depth profile (five treatments x three repititions x three cores).
|
|
Temporal Changes in Bromide Distribution in the Soil Profile
Variability in Br- concentration in samples collected from different N treatments and replicates was nonsignificant (P
0.05). Data were pooled and the mean Br- concentration was calculated across all N rates and replicates (Fig. 4)
. Soil samples taken 7 d after Br- application showed 20 mg kg-1 Br- in the top 15-cm depth, which decreased to 10 mg kg-1 by 14 d following a 15.3-mm irrigation water application (Fig. 4). The Br- concentration in the 0- to 15-cm layer of soil further decreased to about 2 mg kg-1 28 d after surface application of Br-. For the same period, the cumulative water input was 57.6 mm cumulative (Fig. 1). A trace amount (0.04 mg kg-1) of Br- appeared at the bottom of this soil profile for the same period. However, continuous monitoring of Br- indicated that the Br- concentration increased up to 1 mg kg-1 at the bottom of the sampling depth (2.4 m) 98 d after Br- application. This is an indication of downward movement of Br- beyond the rooting depth with drainage water. During the period in between 28 and 98 d, the cumulative water input (rainfall + irrigation) was 410 mm, of which 293 mm of water drained below the 2.4-m depth (Fig. 1).

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 4. Distribution of measured and Leaching Estimation and Chemistry Model (LEACHM)-simulated Br- concentrations in the soil profile of a Tavares fine sand at various sampling times following the application of 80 kg Br- ha-1 as KBr under the citrus tree canopy. Measured values are the mean and standard error of 45 samples per depth profile (five treatments x three repititions x three cores).
|
|
Simulation of water movement through the soil profile using LEACHM indicated that drainage water beyond the soil profile depth of 2.4 m began 28 d after Br- and N application (Fig. 1). Presence of Br- was continuously observed at each sampling event at the bottom of the soil profile (2.4 m) up to 98 d after surface application. However, Br- became nondetectable through the entire depth of sampling by 214 d, indicating that all of the surface-applied Br- was taken up by the trees and leached through the 2.4-m-deep soil profile. Cumulative water input from 98 to 214 d was 872 mm, of which about 434 mm water drained below 2.4 m. (Fig. 1). The cumulative water input during the 98- to 214-d period represented 37 and 50% of the annual irrigation and rainfall, respectively, coinciding with the rainy months (MaySeptember) in the study area. The quantity of water leached during this period was 42% of the average total annual drainage. Water and bromide mass balance for the entire period of this study is given Table 4.
Measured Temporal Changes in Nitrate Distribution in the Soil Profile
Nitrification occurs quite rapidly under these prevailing conditions (<7 d) as demonstrated by Khakural and Alva (1995). Therefore, although N was applied as NH4NO3, transport of NO3N in the soil profile was evaluated in this study. Soil samples were taken to quantify the background levels of NO3N and Br- prior to application of the respective sources 120 d after the last fertilizer application of the previous cropping year (1996). Variability in background concentrations of soil profile NO3N in samples collected at different depths across different N treatments and replicates was nonsignificant (P
0.05). Therefore, the data were pooled and the mean NO3N concentration was calculated across all N rates and replicates and reported here. The background concentrations of soil profile NO3N were varied within a range of 0.3 to 0.8 mg kg-1 in the soil profile (data not presented). The concentrations of NO3N in the soil samples within the top 0.3 m were reasonably related to the rate of N application during the first 28 d of fertilizer application. Thereafter, the differences were not significant (P
0.05) and approached the background concentrations (0.2 to 0.9 mg kg-1) by 42 d (Fig. 5)
. These plots received about 213 mm rainfall and about 58 mm irrigation during 28 d after fertilizer application, which could have resulted in significant downward transport of nitrate. However, the water balance calculations (net water input; i.e., cumulative water input minus losses) and the soil water content measurements from various depth samples (with the exception of the soil samples taken on 28 d and thereafter) showed no water leaching during the first 28 d after the fertilizer application (Fig. 1 and 3). The nutrient requirement of a citrus tree is maximum from early spring through early summer (Tucker et al., 1995). Since the majority of active fibrous roots are in the top 15- to 30-cm depth (Alva and Syvertsen, 1991; Zhang et al., 1996), the effect of root uptake on soil moisture and other plant nutrients would be greatest at this depth. Nutrient uptake is rapid and substantial during the initial 30 d after the fertilizer application. This explains the rapid and significant decrease in NO3N concentrations in the top 15- to 30-cm depth during the 28 d after N application. The results also show some redistribution of N from the surface into the deeper soil with the downward-percolating water front.

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 5. Distribution of measured and Leaching Estimation and Chemistry Model (LEACHM)-simulated NO3N in the soil profile of a Tavares fine sand at various sampling times following the broadcast application of water-soluble dry granular fertilizer containing NH4NO3. Measured values are the mean and standard error of nine samples per depth profile per treatment (three repititions x three cores).
|
|
Nitrate N concentrations through the 2.4-m-deep soil profile decreased to almost background concentrations regardless of N rates by 42 d after fertilizer application. For the same period, the cumulative water input (irrigation and rainfall) was 308 mm. During 28 to 42 d after the fertilizer application, LEACHM simulation indicated a loss of about 122 mm water through drainage beyond the 2.4-m depth in this soil profile (Fig. 1). Disappearance of most N from the soil profile during the 28- to 42-d period could be due to NO3N leaching below the rootzone. Leaching loss of N during this period could be from residual NO3N in the soil profile, decomposition and transformation of plant residue, or residual fertilizer N. Plant residue mineralization, including root turnover, has been reported to contribute up to 156 kg N ha-1 yr-1 for bearing citrus trees under field conditions (Dou et al., 1997). Mineralization parameters used for the current work were based on the work of Clark (1994) and Dou et al. (1997) conducted under Florida sandy soils and used for LEACHM. Further, mineralization rate used in this LEACHM simulation was verified with similar simulation studies conducted in nearby citrus-growing sandy Entisols of central Florida (Pennel et al., 1990; Lamb, 1996; Graham and Wheaton, 1999; Lamb et al., 1999; Paramasivam et al., 2000a). More details about mineralization can be obtained from Clark (1994) and Dou et al. (1997). Mineralized N would affect the distribution pattern of NO3N in the soil profile and leaching, although the contribution may vary with time since the mineralization process is affected by various factors in the soil environment (i.e., soil temperature and water content).
Comparison of Measured and Predicted Nitrate and Bromide Values
Under mature citrus tree growth conditions, concentrations of NO3N at various depths within the entire profile were predicted reasonably by the LEACHM and compared favorably with the measured values at various sampling events across all N rates (Fig. 5). Variations between the measured and predicted NO3N values were within 1 to 2 mg kg-1 over most of the sampling period. Similarly, predicted Br- concentrations were very close to measured Br- values for all the sampling events at various sampling points. Weekly profile soil sampling indicated that the concentrations of NO3N and Br- in the top 15 cm were almost reduced by half of their respective concentrations at the previous sampling at least for the first 28 d of their application. Since the sampling was carried out following every irrigation or rainfall event, this could be considered as a similar transport response of NO3N and Br- to downward movement of water. The LEACHM model predicted that about 40% of the applied Br- was taken up by the trees and the remaining 60% leached below 2.4 m of the simulated profile (Table 5). Similar high Br- uptake by different crops has been reported in the literature. Owens et al. (1985) found that plant uptake of Br- accounted for about 32% of that applied, which was about the same amount that was recovered in drainage water. Even high Br- uptake by plants was reported by Kung (1990) and Schnabel et al. (1995), who found that as much as 55 and 38% of the applied Br- was taken up by potato (Solanum tuberosum L.) and ryegrass (Lolium perenne L.), respectively. This high Br- plant uptake simulated by LEACHM could be due to citrus tree uptake, since the application time coincided with a period of very high plant nutrient uptake demand (i.e., N; Tucker et al., 1995). However, we do not have analytical information to ascertain the Br- plant uptake component other than Br- plant uptake simulated by LEACHM. Therefore, based on the available information, leaching and plant uptake would have been the major processes responsible for the nondetectable amount of Br- in the soil profile at the end of the experiment.
View this table:
[in this window]
[in a new window]
|
Table 5. Nitrogen mass balance for different fertilization rates as simulated by the Leaching Estimation and Chemistry Model (LEACHM).
|
|
Predicted Cumulative Losses from Applied Fertilizer Nitrogen and Nitrogen Mass Balance
Nitrogen mass balance simulated by LEACHM is presented in Table 5. A one-fourth split of the annual application of various N rates (28, 56, 84 and 112 kg N ha-1) resulted in loss of 15, 15, 28, and 30 kg N ha-1, respectively, with 122 mm cumulative water leached during the 42-d period. Various factors could have contributed to N loss from the soil profile, probably a significant contribution from applied fertilizer N and mineralized N under the tree canopy. Even though it is difficult to separate the contribution of mineralized N for leaching, we can estimate the leaching losses of N from applied fertilizer N on a proportional basis. For the calculation of amount of N leached from applied fertilizer, it was assumed that the contribution of mineralized N for leaching had been consistent across all N treatments. Based on these calculations, about 10 (i.e., {[15/(28 + 14)] x 28}), 12 (i.e., {[15/(56 + 14)] x 56}), 24 (i.e., {[28/(84 + 14)] x 84}), and 27 (i.e., {[30/(112 + 14)] x 112}) kg N ha-1 loss came from applied fertilizer N during this 42-d period. Thus, 36, 21, 29, and 24% of the applied fertilizer N from low to high N fertilizer rates, respectively, would have been lost in the leachate. As per LEACHM prediction, other components of N mass balance such as plant uptake, volatilization losses, and variation in soil N storage increased with increasing rate of mineral N fertilizer input.
Analyses of soil cores from plots receiving various amounts of N as dry granular fertilizer revealed that there was very little variation in movement and distribution of NO3N in the soil profile across all N rates. This observation clearly suggests a negligible preferential flow effect under field conditions of this study and further justifies the use of LEACHM. Preferential flow would be an important phenomenon under which soil has layering, cracks, and fissures as a result of the presence of shrinking and swelling clay, natural soil pipes, and dead plant roots and residues (Hillel, 1998). These preferential flow factors are unlikely to be present under a structureless homogeneous sandy soil (>96% sand). Thus, it is not expected that preferential flow could play a dominant role in water and solute movement in extremely sandy soils.
 |
SUMMARY AND CONCLUSION
|
|---|
Transport and distribution of water, NO3N, and Br- in a sandy soil profile were studied under mature citrus tree production conditions. Nitrate N concentrations through the 2.4-m-deep soil profile decreased to levels very close to the background concentrations by about 42 d, irrespective of N rates. Rapid decrease of NO3N concentrations in the top 15- to 30-cm depth occurred during the 28 d after N application, possibly attributable to high plant uptake demand during the early spring season. The results also showed some downward movement of N with wetting front following water application from the surface layer deeper into the soil profile, which could have contributed to a decrease in the NO3N concentration in the topsoil. Bromide, used in this study as a tracer, indicated a somewhat similar trend of distribution and disappearance pattern as fertilizer NO3N. However, about 214 d was taken for the surface-applied Br- to become nondetectable within the 2.4-m soil profile. The LEACHM simulations and mass balance of simulations indicated an increasing trend of leaching losses of applied N (15 to 30 kg ha-1) with increasing N rate of application used in this study. This study also indicated the usefulness of a numerical modeling approach to estimate the N leaching losses and predict other mass balance components of N and water simultaneously and accurately for the entire crop year. Results of this study further suggest that, under extremely sandy nature of soils, leaching losses and ground water contamination of NO3N are driven by excess water leaching and thus could be minimized by maintaining irrigation practices to supply water only to replenish the water deficit within the rooting depth of citrus and applying N fertilizers in small quantities but with more frequent applications. Results also suggest that repeated intensive monitoring of nutrient movement in the soil profile following each split fertilizer application coupled with a numerical modeling exercise might help us to enhance our understanding of N dynamics more thoroughly and to fine tune N and irrigation management practices in similar types of sandy soils. However, the authors are not advocating repeated intensive sampling following each split fertilizer application by every farmer because it is very laborious and not feasible.
 |
ACKNOWLEDGMENTS
|
|---|
This study was made possible, in part, by funding from the Florida Department of Agriculture and Consumer Services (DACS), Tallahassee, FL. We appreciate the support of Richard Budell and Marlene Czerniak of the Florida Department of Agriculture and Consumer Services. We also appreciate assistance from K.H. Hostler and J.T. Kelley. We are grateful to W. Dubberly, J. Whitakar, and A.Tucks of Kahn Grove Incorporated for their cooperation and assistance.
 |
NOTES
|
|---|
Contribution of the Citrus Research and Education Center, Florida Agricultural Experiment Station Journal Series no. R-07149.
 |
REFERENCES
|
|---|
- Alva, A.K., and S. Paramasivam. 1998. Nitrogen management for high yield and quality of citrus in sandy soils. Soil Sci. Soc. Am. J. 62:13351342.[Abstract/Free Full Text]
- Alva, A.K., and J.P. Syvertsen. 1991. Irrigation water salinity affects soil nutrient distribution, root density, and leaf nutrient levels of citrus under drip fertigation. J. Plant Nutr. 14:715727.
- Bowman, R.S. 1984. Evaluation of some new tracers for soil water studies. Soil Sci. Soc. Am. J. 48:987993.[Abstract/Free Full Text]
- Campbell, G.S. 1974. A simple method for determining unsaturated conductivity from moisture retention data. Soil Sci. 117:311314.
- Campbell, K.L. 1978. Pollution in runoff from nonpoint sources. Publ. 42. Research project technical completion report. OWRT Project no. B-023-FLA. Florida Water Resour. Res. Center, Gainesville.
- Carlisle, V.W., F. Sodek, III, M.E. Collins, L.C. Hammond, and W.G. Harris. 1989. Characterization data for selected Florida soils. Soil Sci. Res. Rep. 89-1. Univ. of. Florida, Soil Conserv. Soc. and USDA, Gainesville.
- Castle, W.S. 1980. Fibrous root distribution of pineapple orange trees on rough lemon root stock at three tree spacings. J. Am. Soc. Hortic. Sci. 105:478480.
- Clark, J. 1994. Evaluation of nitrogen movement beneath a leatherleaf fernery. Master's thesis. Univ. of Florida, Gainesville.
- Clothier, B.E., T.J. Sauer, and S.R. Green. 1988. The movement of ammonium nitrate into unsaturated soil during unsteady absorption. Soil Sci. Soc. Am. J. 52:340345.
- Comfort, S.D., W.P. Inskeep, and R.H. Lockerman. 1993. Observed and simulated transport of a conservative tracer under line-source irrigation. J. Environ. Qual. 22:543561.
- Davies, F.S. 1996. Literature review of research related to nitrogen nutrition, fertilization and groundwater pollution of citrus. Florida Dep. of Agric. and Consumer Serv., Tallahassee.
- Dou, H., A.K. Alva, and B.R. Khakural. 1997. Nitrogen mineralization from citrus tree residues under different production conditions. Soil Sci. Soc. Am. J. 61:12261232.[Abstract/Free Full Text]
- Fares, A. 1996. Environmental impact of unharvested forest buffer zones upon cypresspond systems in coastal plains regions: Modeling analyses. Ph.D. diss. Univ. of Florida, Gainesville.
- Fares, A., and R.S. Mansell. 1996. A potential evapotranspiration prediction submodel for multi-dimensional water flow in soils. Soil Crop Sci. Soc. Florida Proc. 55:3645.
- Foussereau, X., and W.D. Graham. 1997. Field-scale subsurface transport of a surface applied tracer at a south-west Florida citrus grove. Soil Crop Sci. Soc. Florida Proc. 56:7882.
- Graham, W.D., and T.A. Wheaton. 1999. Ridge citrus water quality project. Progress report submitted to Florida Dep. of Agric. and Consumer Serv., Tallahassee. Florida Dep. of Environ. Protection, Tallahassee, and Southwest Florida Water Management Dist., Brooksville.
- Hallberg, G.R. 1989. Nitrate in groundwater in the United States. p. 3574. In R.F. Follet (ed.) Nitrogen management and groundwater protection. Elsevier Sci. Publ., Amsterdam.
- Hillel, D. 1998. Environmental soil physics. Academic Press, San Diego, CA.
- Hornsby, A.G., P.S.C. Rao, and R.L. Jones. 1990. Fate of aldicarb in the unsaturated zone beneath a citrus grove. Water Resour. Res. 26:22872302.
- Hubbard, R.K., G.J. Gascho, J.E. Hook, and W.G. Knisel. 1986. Nitrate movement into shallow groundwater through a Coastal Plain sand. Trans. ASAE 29:15641571.
- Hubbard, R.K., R.G. Williams, and M.D. Erdman. 1989. Chemical transport from Coastal Plain soils under simulated rainfall: I. Surface runoff, percolation, nitrate, and phosphate movement. Trans. ASAE 32:12391249.
- Hutson, J.L., and R.J. Wagenet. 1992. LEACHM: Leaching Estimation and Chemistry Model: A process-based model of water and solute movement, transformation, plant uptake and chemical reactions in the unsaturated zone. Water Resour. Inst., Cornell Univ., Ithaca, NY.
- Jabro, J.D., J.M. Jemison, Jr., L.L. Lengnick, R.H. Fox, and D.D. Fritton. 1993. Field validation and comparison of LEACHM and NCSWAP models for predicting nitrate leaching. Trans. ASAE 36:16511657.
- Johnson, H., L. Bergstrom, P.E. Jansson, and K. Paustian. 1987. Simulated nitrogen dynamics and classes in a layered agricultural soil. Agric. Ecosyst. Environ. 18:333356.
- Khakural, B.R., and A.K. Alva. 1995. Hydrolysis of urea in two sandy soils under citrus production as influenced by rate and depth of placement. Commun. Soil Sci. Plant. Anal. 26:21432156.
- Klute, A. 1986. Water potential: Tensiometry. p. 635662. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Klute, A., and C. Dirksen. 1986. Hydraulic conductivity and diffusivity: Laboratory methods. p. 687734. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Kung, K.J.S. 1990. Influence of plant uptake on the performance of bromide tracer. Soil Sci. Soc. Am. J. 54:975979.[Abstract/Free Full Text]
- Lamb, S.T. 1996. Monitoring and modeling the impacts of alternative nutrient and water management practices on groundwater quality in the central Florida Ridge. Master's thesis. Univ. of Florida, Gainesville.
- Lamb, S.T., W.D. Graham, C.B. Harrison, and A.K. Alva. 1999. Impact of alternative citrus management practices on groundwater nitrate in the central Florida Ridge, Part I. Field investigation. Trans. ASAE 42:16531668.
- Mansell, R.S., D.V. Calvert, E.H. Stewart, W.B. Wheeler, J.S. Rogers, D.A. Greatz, L.H. Allen, A.R. Overman, and E.B. Knipling. 1977. Fertilizer and pesticide movement from citrus groves in Florida flatwood soils. Completion report for Project R-800517. USEPA Environ. Res. Lab., Athens, GA.
- Mansell, R.S., J.G.A. Fiskell, D.V. Calvert, and J.S. Rogers. 1986. Distributions of labeled nitrogen in the profile of a fertilized sandy soil. Soil Sci. 141:120126.
- Molina, J.A.E., and K. Richards. 1984. Simulation models of the nitrogen and carbon cycle in the soilwaterplant system, NCSWAP. Guide for the preparation of input data files and execution of NCSWAP. Soil Series 116. Dep. of Soil Sci., Univ. of Minnesota, St. Paul.
- Nightingale, H.I. 1972. Nitrates in soil and groundwater beneath irrigated and fertilized crops. Soil Sci. 114:300311.
- Owens, L.B., R.W. Van Keuren, and W.M. Edwards. 1985. Groundwater quality changes resulting from a surface bromide application to a pasture. J. Environ. Qual. 14:543548.[Abstract/Free Full Text]
- Paramasivam, S., A.K. Alva, and A. Fares. 2000a. Transformation and transport of nitrogen forms in a sandy Entisol following a heavy loading of ammonium nitrate solution: Field measurements and model simulations. J. Soil Contam. 9:6586.
- Paramasivam, S., A.K. Alva, and A. Fares. 2000b. An evaluation of soil water status using tensiometers in a sandy soil profile under citrus production. Soil Sci. 165:343353.
- Paramasivam, S., A.K. Alva, A. Fares, and K.S. Sajwan. 2001. Estimation of nitrate leaching in an entisol under optimum citrus production. Soil Sci. Soc. Am. J. 65:914921.[Abstract/Free Full Text]
- Paramasivam, S., A.K. Alva, and W.D. Graham. 1999. Transport of bromide in an entisol and its dissipation in a surficial aquifer. J. Environ. Sci. Health A34:585604.
- Parsons, L.R. 1989. Management of micro-irrigation systems for Florida citrus. Fruit Crops Fact Sheet. Univ. of Florida Coop. Ext. Serv. FC-81. Univ. of Florida, Gainesville.
- Pennell, K.D., A.G. Hornsby, R.E Jessup, and P.S.C. Rao. 1990. Evaluation of five simulation models for predicting aldicarb and bromide behavior under field conditions. Water Resour. Res. 26: 26792693.
- Priestley, C.H.B., and R.J. Taylor. 1972. On the assessment of surface heat flux and evaporation using large scale parameters. Monthly Weather Rev. 100:8192.
- Schnabel, R.R., W.L. Stout, and J.A. Shaffer. 1995. Uptake of a hydrologic tracer (bromide) by ryegrass from well and poorly-drained soils. J. Environ. Qual. 24:888892.[Abstract/Free Full Text]
- Smajstrala, A.G., D.S. Harrison, F.S. Zazutta, L.R. Parsons, and K.C. Stone. 1987. Trickle irrigation scheduling for Florida citrus. Bull. 208. Florida Coop. Ext. Serv., Inst. Food and Agric. Sci., Univ. of Florida, Gainesville.
- Smith, S.J., and R.J. Davis. 1974. Relative movement of bromide and nitrate through soils. J. Environ. Qual. 3:152155.[Abstract/Free Full Text]
- Soulsby, C., and B. Reynolds. 1992. Modeling hydrological processes and aluminum leaching in an acid soil at Llyn Brianne, Mid-Wales. J. Hydrol. 138:409429.
- Staver, K.W., and R.B. Brinsfield. 1990. Patterns of soil nitrate availability in corn production systems: Implications for reducing groundwater contamination. J. Soil Water Conserv. 45:318322.
- Stevenson, F.J. 1982. Origin and distribution of nitrogen in soils. p. 139. In F.J. Stevenson (ed.) Nitrogen in agricultural soils. Agron. Monogr. 22. ASA, Madison, WI.
- Tucker, D.P.H., A.K. Alva, L.K. Jackson, and T.A. Wheaton. 1995. Nutrition of Florida citrus trees. Univ. of Florida Coop. Ext. Serv. SP169. Univ. of Florida, Gainesville.
- USEPA. 1991. Test method. The determination of inorganic anions in water by ion chromatography. Environ. Monitoring Systems Lab., Cincinnati, OH.
- Wagenet, R.J., and J.L. Hutson. 1989. LEACHM: Leaching Estimation and Chemistry Model: A process based model: A process based model of water and solute movement transformations, plant uptake and chemical reactions in the unsaturated zone. Continuum Vol. 2. Water Resour. Inst., Cornell Univ., Ithaca, NY.
- Zhang, M., A.K. Alva, Y.C. Li, and D.V. Calvert. 1996. Root distribution of grapefruit trees under dry granular dry broadcast vs. fertigation method. Plant Soil 183:7984.