# Journal of Soil and Water Science

ISSN: 2643-5799

### Article Outline

RESEARCH ARTICLE | VOLUME 1 | ISSUE 1 | DOI: 10.36959/624/426 OPEN ACCESS

# Characterizing Sorption and Modeling Phosphorus Movement on Candler and Immokalee Fine Sand

Davie M Kadyampakeni, Peter Nkedi-Kizza, Kelly T Morgan and Arnold W Schumann

• Peter Nkedi-Kizza 2
• Kelly T Morgan 3
• Arnold W Schumann 1
• Citrus Research and Education Center, University of Florida, USA
• Soil and Water Science Department, University of Florida, USA
• Southwest Florida Research and Education Center, University of Florida, USA

Kadyampakeni DM, Kizza PN, Morgan KT, et al. (2017) Characterizing Sorption and Modeling Phosphorus Movement on Candler and Immokalee Fine Sand. J Soil Water Sci 1(1):8-14.

Accepted: February 18, 2017 | Published Online: February 22, 2017

# Characterizing Sorption and Modeling Phosphorus Movement on Candler and Immokalee Fine Sand

## Abstract

Laboratory and computer experiments were conducted to 1) determine the effect of supporting electrolyte on mass distribution coefficient (KD) for predicting P movement at 30- and 60-cm soil depths using HYDRUS-1D and 2) compare the temporal P concentrations as a function of soil type. The results on Candler fine sand at Lake Alfred showed that P contents for the KD estimated with 0.01 M KCl and 0.005 M CaCl2 were 10 to 15% higher than those predicted with a KD value measured with fertilizer mixture. The predictions on Immokalee fine sand showed that P contents for the KD estimated with fertilizer mixture and 0.005 M CaCl2 were 12 to 20% higher than those predicted with a KD value measured with 0.01 M KCl. The outputs with KD measured with 0.005 M CaCl2 appear to be close to those predicted with a KD measured with fertilizer mixture on Immokalee fine sand. However, the analysis of the KD values across all electrolytes on the two soils studied revealed that 0.01 M KCl is the electrolyte that yields KD values fairly close to fertilizer mixture on Candler fine sand and would be appropriate to use for the coated sand while 0.005 M CaCl2 was comparable with fertilizer mixture on Immokalee fine sand and would be an appropriate electrolyte to use for such uncoated sands. The results of study will help in improving estimation of KDs and description of P transport in citrus production systems of Florida sandy soils. This would provide necessary information for sustainable environmental management and reduce problems of eutrophication and prevent groundwater contamination.

## Keywords

Phosphorus adsorption, Simulation modeling, Entisol, Spodosol

## Introduction

Adsorption is the mechanism that soil cations and anions (including various phosphorus (P) compounds) are retained by soil particles. The adsorption process tends to restrict compound's mobility and bioavailability [1]. Thus, the procedure for determining the P sorption isotherms could then provide information on their mobility in the soil. The supporting electrolyte concentration is chosen to mimic that of soil solution. Most commonly 0.01 M CaCl2 [2,3], 0.01 N CaCl2 [4], 0.005 M CaCl2 [1], 0.1 to 3 mM KCl [5], 0.05 M KCl [6,7], and 0.01 M KCl [8,9] have been used as electrolytes in studies on P and K sorption. Nair, et al. (1984) [10] reported that P sorption varies with ionic strength and cation species of the supporting electrolyte. For example, [10] showed that P adsorption was generally lower with K+ as the supporting electrolyte cation compared with Ca2+. These studies and others have not explained the rationale behind use of a particular electrolyte other than equilibrating the solutions in deionized or tap water.

The chemical characteristics of soils dominating the Flatwoods and Ridge regions of Florida are well described in [11] and some were also determined in this study. The soil at the Flatwoods consists of nearly level and poorly drained on the Flatwoods classified as Immokalee fine sand (sandy, siliceous, hyperthemic Arenic Haplaquods) with the spodic horizon lying within 1 m from the soil surface. The soil at the Ridge, classified as Candler fine sand (hyperthermic, coated Typic Quartzipsamments) is a well-drained sandy soil with no continuous layer limiting vertical water movement [11]. The Immokalee and Candler fine sand are moderately acidic (pH ranging from 4.9 to 5.6) with > 94% sand textural composition, have low organic matter content (ranging from 0.41 to 0.61% on Immokalee fine sand and from 1.6 to 1.96% on Candler fine sand) and low cation exchange capacity (CEC) (ranging from 2 to 6 cmol (+) kg-1), have inorganic N in the range of 8.2 and 11.2 mg kg-1, moderate to very high P (in the range of 28.7 to 46.5 mg kg-1 for Immokalee sand and 112.8 to 115.8 mg kg-1 for Candler fine sand) and K in the range of 11.8 to 15.2 mg kg-1 for Immokalee fine sand and 23.0 to 29.7 mg kg-1 for Candler fine sand (Table 1) [12]. It has been accepted that soil CEC, particularly the CEC contribution from organic matter content (OMC) has a significant influence on the soil P movement [13]. This study sought to 1) determine the effect of supporting electrolyte on the mass distribution coefficient (KD) for predicting P movement at 30- and 60-cm soil depths using HYDRUS-1D; and 2) compare the P concentrations with time as a function of soil type.

## Materials and Methods

The baseline soil chemical and physical properties for soils used in the study are described in (Tables 1 and Table 2). Sorption isotherms on the disturbed soil samples (0-15 cm, 15-30 cm) were determined using the batch equilibration procedure [14]. The initial solution concentrations from potassium dihydrogen phosphate (KH2PO4) for P in 0.005 M CaCl2 were 10, 25, 50 mg PL-1. To determine P sorption isotherms, soil samples were obtained from 5 random positions per site at two depths giving a total of 10 samples. Each sample was weighed in triplicates plus a blank control. A 10 g air-dried, subsample of soil < 2 mm particle size was placed in a centrifuge tube and equilibrated with 20 ml (soil solution ratio 1:2) of 3 initial concentrations of P solutions. The centrifuge tubes were shaken for 24 h, centrifuged for 20 min, and filtered through a filter paper (Whatman, #42). All these procedures were done at room temperature ~25 ± 1 °C as recommended by [14] but the filtrate was later stored at < 4 °C until analyzed for P.

### Phosphorus sorption isotherm determination

The solutions from the adsorption study were analyzed for P using ICP-AES and calibration standards of 10, 30 and 50 mg PL-1 preparation in fertilizer mixture (with additions of 0.027 mM NH4NO3 and 0.013 mM KCl), 0.01 M KCl and 0.005 M CaCl2. The amount of chemical sorbed to the soil was calculated from the difference between the initial and equilibrium solution concentration:

Where S is the adsorbed concentration (mg kg-1); Vo is the volume of initial solution (L); m is the soil mass (kg); Co is the initial concentration of the standard solution (mg L-1), and, Ce is the soil solution concentration at equilibrium (mg L-1).

Sorption isotherms for P were calculated using the Freundlich equation:

Where Kf is the Freundlich sorption coefficient (mg1-N kg-1 LN) and N is an empirical constant related to adsorption phenomena [15]. The linearized form of the Freundlich equation was used to calculate Kf and N:

Where S is the adsorbed equilibrium concentration (mg kg-1); Ce is the equilibrium concentration (mg L-1) and Cmax is the estimated maximum concentration (mg L-1) and Kf and N are calculated from the intercept and slope of Equation 3. To find average linearized KD for the Freundlich isotherm, the integrated form of the equation was used:

### Concepts and governing equations for the simulations

The governing flow equations for water flow and nutrient transport are given by the [16] and convection-dispersion equations (CDE) [17-19]:

Where θ is the volumetric water content [L3L-3], h is the pressure head [L], xi (I = 1, 2) are the spatial coordinates [L], t is time [T], ${K}_{ij}^{A}$ are components of a dimensionless anisotropy tensor ${K}^{A}$ (which reduces to the unit matrix when the medium is isotropic), K is the unsaturated hydraulic conductivity function (LT-1), and s is a sink/source term [L3L-3T-1], accounting for root water uptake (transpiration). The sink/source represents the volume of water removed per unit time from a unit volume of soil due to compensated citrus water uptake.

The equation (CDE) governing transport of independent solutes i.e. single-ion transport is given as:

Where c1 and c2 are solute concentrations in the solid (MM-1) and liquid (ML-3) phases, respectively; qi is the ith component of volumetric flux density (LT-1), Ф is the rate of change of mass per unit volume by chemical or biological reactions or other sources (negative) or sinks (positive) (ML-3T-1), respectively, providing connections between individual chain species, ρb is the soil bulk density (ML-3), Dij is the dispersion coefficient tensor for the liquid phase [L2T-1]. The term ra represents the root nutrient uptake (ML-3T-1) which is the sum of actual active and passive nutrient uptake. The solid phase concentration, c1, accounts for nutrient either sorbed in the solid phase or precipitated in various minerals. This is usually quantified by the adsorption isotherm relating c1 and c2 described by the linear equation of the form:

Where KD (L3 M-1) is the mass distribution coefficient of species 1. A tracer (e.g. bromide) are assumed to have a KD = 0 cm3 g-1. The first order decay constant ranges from 0.36-0.56 d-1 [20]. For P, KD is reportedly in the range of 19 to 185 cm3 g-1 [21,22]. Bulk density for the soil is in the range 1.59-1.72 g cm-3 (Immokalee) and 1.55-1.93 g cm-3 (Lake Alfred) (T.A. Obreza, unpublished).

The sink term, s, for the Richards equation represents the volume of water removed per unit time from a unit volume of soil due to plant water uptake. Thus, s is defined as:

Where the water stress response function $\alpha \left(h\right)$ is a prescribed dimensionless function of the soil water pressure head, b is the normalized water uptake distribution, Lt is the width of the soil surface associated with the transpiration process and Tp is the potential transpiration rate (LT-1) and w is the water stress index.

The nonlinear, predictive equations for the unsaturated hydraulic function in terms of soil water retention parameters are given by [23] as:

Where

Where θr, θs, Ks and l are residual water content (L3L-3), saturated water content (L3L-3), saturated hydraulic conductivity (LT-1), and pore connectivity parameter (estimated to be an average of 0.5 for many soils respectively). Parameters α (L-1) and n are empirical coefficients affecting the shape of the hydraulic functions [23]. We estimated the hydraulic functions α and n after fitting the water content and matric potential data using the van Genuchten model in Community Analyses System (CAS) 2007 [24] developed for determination of soil hydraulic functions (Table 3).

### Cumulative flux and concentration

The P concentration and flux were predicted at 30 cm soil depth. Assuming isotropy and a homogeneous profile, the concentration and fluxes at 60 cm soil depth were also predicted.

## Results and Discussion

### Phosphorus movement with microsprinkler irrigation as function of time

Phosphorus movement was predicted using three different KDs estimated with fertilizer mixture, 0.01 M KCl and 0.005 M CaCl2 for of 21 days, assuming no rainfall events (Figure 3). The assumption is that a KD value obtained using fertilizer mixture typifies that of field conditions with regard to chemical processes. The results on Candler fine sand at Lake Alfred showed that P contents for the KD estimated with 0.01 M KCl and 0.005 M CaCl2 were 10 to 15% higher than those predicted with a KD value measured with fertilizer mixture. The predictions on Immokalee fine sand showed that P contents for the KD estimated with fertilizer mixture and 0.005 M CaCl2 were 12 to 20% higher than those predicted with a KD value measured with 0.01 M KCl. The outputs with KD measured with 0.005 M CaCl2 appear to be close to those predicted with a KD measured with fertilizer mixture on Immokalee fine sand. Thus, the analysis of the KD values across all electrolytes on the two soils studied revealed that 0.01 M KCl is the electrolyte that yields KD values (0.42 to 0.69 L kg-1 on Immokalee fine sand and 0.42 to 1.06 L kg-1 on Candler fine sand) fairly close to fertilizer mixture on coated sands while 0.005 M CaCl2 tends to give KD values (2.44 to 4.66 L kg-1 on Immokalee fine sand and 2.81 to 4.93 L kg-1 on Candler fine sand) two to threefold in magnitude to those determined with fertilizer mixture on sands with Fe coatings. This suggests that 0.005 M CaCl2 would tend to overestimate P sorption and retardation during unsaturated or saturated flow on coated sand while estimating the adsorption process satisfactorily on uncoated sands than 0.01 M KCl. It appears the addition of a supporting electrolyte with a divalent or monovalent cation, unlike fertilizer mixture, increases the surface charge for adsorption of orthophosphate anions on sands with coatings. The use of 0.01 M KCl appears not to be influenced by presence of sand coatings compared with CaCl2 and, thus, would present an appropriate supporting electrolyte for Candler fine sand while CaCl2 would be appropriate for the Immokalee fine sand.

## Conclusion

The results show that P adsorption in the top 0-15 cm was greater for Candler than Immokalee sand using the fertilizer mixture, 0.005 M CaCl2 and 0.01 M KCl. The mass distribution coefficients (KD) for P estimated using 0.01 M KCl were similar to values determined using fertilizer mixture for Candler fine sand. The KD values determined using 0.005 M CaCl2 as the supporting electrolyte were two- to threefold greater than the KD of the fertilizer mixture on Candler fine sand suggesting that divalent Ca might result in overestimation of P sorption on Candler sandy soils. On Immokalee fine sand, fertilizer mixture and 0.005 M CaCl2 were comparable and resulted in greater P adsorption than 0.01 M KCl. Thus, it would be appropriate to use 0.01 M KCl as supporting electrolyte for Florida's Candler fine sand and 0.005 M CaCl2 on Immokalee fine sand.

## Abstract

Laboratory and computer experiments were conducted to 1) determine the effect of supporting electrolyte on mass distribution coefficient (KD) for predicting P movement at 30- and 60-cm soil depths using HYDRUS-1D and 2) compare the temporal P concentrations as a function of soil type. The results on Candler fine sand at Lake Alfred showed that P contents for the KD estimated with 0.01 M KCl and 0.005 M CaCl2 were 10 to 15% higher than those predicted with a KD value measured with fertilizer mixture. The predictions on Immokalee fine sand showed that P contents for the KD estimated with fertilizer mixture and 0.005 M CaCl2 were 12 to 20% higher than those predicted with a KD value measured with 0.01 M KCl. The outputs with KD measured with 0.005 M CaCl2 appear to be close to those predicted with a KD measured with fertilizer mixture on Immokalee fine sand. However, the analysis of the KD values across all electrolytes on the two soils studied revealed that 0.01 M KCl is the electrolyte that yields KD values fairly close to fertilizer mixture on Candler fine sand and would be appropriate to use for the coated sand while 0.005 M CaCl2 was comparable with fertilizer mixture on Immokalee fine sand and would be an appropriate electrolyte to use for such uncoated sands. The results of study will help in improving estimation of KDs and description of P transport in citrus production systems of Florida sandy soils. This would provide necessary information for sustainable environmental management and reduce problems of eutrophication and prevent groundwater contamination.