JEQ Grow Your Career With ASA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (18)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Fisher, M.M.
Right arrow Articles by Reddy, K.R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Fisher, M.M.
Right arrow Articles by Reddy, K.R.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Fisher, M.M.
Right arrow Articles by Reddy, K.R.
Related Collections
Right arrow Wetland Soils
Right arrow Wetlands and Aquatic Processes
Right arrow Nutrient Management
Journal of Environmental Quality 30:261-271 (2001)
© 2001 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America

TECHNICAL REPORT
WETLANDS AND AQUATIC PROCESSES

Phosphorus Flux from Wetland Soils Affected by Long-Term Nutrient Loading

M.M. Fisher and K.R. Reddy

Wetland Biogeochemistry Lab., Soil and Water Science Dep., Univ. of Florida, 106 Newell Hall, P.O. Box 110510, Gainesville, FL 32611

Corresponding author (Krr{at}gnv.ifas.ufl.edu)

Received for publication January 4, 2000.

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Wetland soils play a key role in the cycling of nutrients within an ecosystem. Since soils are potentially a source or a sink for inorganic nutrients, it is important to quantify their influence on overlying water quality in order to understand their importance in overall ecosystem nutrient budgets. Laboratory and field studies were performed in the northern Everglades (WCA-2A) to determine the magnitude of phosphorus (P) flux between the soil and the overlying water column, under various redox conditions. The P flux was estimated using three techniques: intact soil cores, in situ benthic chambers, and porewater equilibrators. There was reasonable agreement between the P flux estimated using intact soil cores and benthic chambers; however, P flux estimates using the porewater equilibrators were considerably lower than the other two techniques. Models of solute flux, based solely on soil physico–chemical characteristics, may substantially underestimate soil–water nutrient exchange processes. Phosphorus flux measured with the intact soil cores varied from 6.5 mg m-2 d-1 near nutrient inflow areas to undetectable flux 4 km away from the inflow. Oxygen consumption varied from 4 mg m-2 d-1 near the inflow to a constant 1 to 2 mg m-2 d-1 at a distance of 4 km from the inflow. Rate of consumption of NO-3–N and SO2-4 showed no significant trend with respect to distance from inflow. Nitrate N and SO4 consumption rates averaged 120 and 130 mg m-2 d-1, respectively. Consumption of O2 was correlated with P flux, whereas NO-3–N and SO2-4 consumption were not.

Abbreviations: BOD, biochemical oxygen demand • DRP, dissolved reactive phosphorus • SOD, soil oxygen demand


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
WETLAND soils can function as sources or sinks for nitrogen (N) and P to the overlying water column. This has implications for environmental managers in that the soils underlying many natural systems can potentially supply nutrients in excess of anthropogenic loading. Therefore, the quantification of the magnitude of nutrient flux, especially N and P, from soil can suggest the relative importance of controlling either internal or external loading of nutrients to an ecosystem.

Wetland soils are characterized by slow turnover of organic material due to a limited supply of terminal electron acceptors. This results in accumulation of organic matter and is one of the principal, long-term removal mechanisms of P and N in treatment wetlands (Kadlec, 1999). The mineralization of organic matter and subsequent nutrient release is predominately controlled by quality of organic material (Godshalk and Wetzel, 1978) and the supply of electron acceptors (D'Angelo and Reddy, 1994a,b). Thus, the supply of electron acceptors is one of the key determinants of organic matter decomposition and release of organically bound nutrients (McLatchey and Reddy, 1998). Thus, the rate of utilization of electron acceptors can indicate the rate of cycling of key nutrients, such as N and P.

Nutrient loading increases net accumulation of organic matter and associated nutrients in wetlands through accelerated primary productivity (Craft and Richardson, 1993; Reddy et al., 1993). However, increased loading not only increases total nutrient content of soils but also increases soluble forms, which can potentially be released into the water column (Reddy et al., 1998). This release of nutrients from soils becomes very critical and can maintain eutrophic conditions in wetlands, even when external loads are curtailed.

Several methods have been used to determine the magnitude and direction of nutrient flux from soils to the overlying water column. Deterministic models of solute flux, such as Fick's First Law of Diffusion, rely on close-interval estimation of porewater concentration gradients (Berner, 1980). These gradients are typically measured using finely sliced soil and sediment samples or by using porewater equilibration (dialysis) devices (D'Angelo and Reddy, 1994a; Moore et al., 1991, 1998). Phosphorus flux can also be determined using more direct approaches, such as by observing changes in water column concentrations inside in situ benthic chambers (Gomez-Parra and Forja 1993; Callender and Hammond, 1982; McCaffrey et al., 1980), or by measuring changes in floodwater concentrations in intact, incubated soil and sediment cores (Moore et al., 1998). The results from these approaches to flux measurement are often compared to derive some estimate of the variability between the techniques, and presumably attempt to determine a "best" measurement technique for a given system. Fickian diffusion has been compared with flux determined with incubated cores and benthic chambers in lakes (Sinke et al., 1990), open ocean (Devol, 1987), nearshore marine (Hopkinson, 1987), and estuarine (Callender and Hammond, 1982) ecosystems, but not in wetlands.

There has been a great deal of research on nutrient exchange dynamics between lake sediments and the overlying water column, but comparatively little research on the influence of wetland soils on water quality. Accordingly, we investigated several factors that are known to influence nutrient exchange processes in order to test the hypothesis that the supply of electron acceptors exerts a control on the rate of flux of P from soils. Additionally, we examined several techniques commonly used to estimate exchange rates. This study will have broad-ranging implications in other wetlands where nutrient exchange rates are of interest.

The objectives of this study were to (i) determine P flux from soil to overlying water column as measured by porewater equilibrators, in situ benthic chambers, and intact soil cores; (ii) determine the role electron acceptors play in governing P flux; and (iii) estimate internal P load to the water column.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Site Description
The study was conducted in the Water Conservation Area 2A (WCA-2A) in the northern Florida Everglades (Fig. 1) . The WCA-2A is a vegetated peat marsh that receives agricultural drainage water from extensive farmland located to the north. Hydrologic and nutrient effects have caused trophic changes: mostly a shift from a sawgrass (Cladium jamaicense Crantz)- to a cattail (Typha domingensis Pers.)-dominated vegetative community (Davis, 1991) and a soil P enrichment gradient (DeBusk et al., 1994) (Table 1). The study locations used in the following experiments were located south of a water control structure (S10-C) that discharges agricultural wastewater into WCA-2A. The P gradient caused by nutrient loading from this structure provided a wide range of soil P levels for comparison of different approaches to estimate P flux.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 1. Location of study site and experimental locations. Porewater equilibrators were installed and intact soil cores were taken at Stations 1 through 8. Benthic flux chambers were used at Stations 1, 3, and 8. Field-dried soil cores were taken at Stations F1, F4, and U3

 

View this table:
[in this window]
[in a new window]
 
Table 1. Physico–chemical properties of surface soils (0–10 cm) collected from stations referenced in this study. Total P, total N, and total C data are mean values (n = 3) from samples collected in February and August 1996 and March 1997 (White, 1999). Porewater pH, NH4–N, and dissolved reactive phosphorus (DRP) from February 1996 replicate equilibrators (n = 2)

 
Greenhouse Experiments
Twenty-four intact soil cores were retrieved from WCA-2A on 12 Sept. 1996. Three cores were taken from each of the eight stations established along a soil P enrichment gradient (Fig. 1). Clear acrylic tubes (14.6 cm i.d. x 51 cm long) were used to retrieve intact soil cores. The tubes were driven into the soil by cutting around the base of the core to a depth of approximately 15 cm with a serrated knife to sever all roots. A large sharpened and serrated shovel was then used to cut the peat and advance the tube into the soil. This procedure resulted in minimal soil compaction. The soil cores were removed, stoppered, and sealed in the bottom, and the floodwater depth was adjusted to 10 cm. The cores were returned to a temperature-controlled greenhouse and covered with an opaque foil shroud to limit primary productivity. Four sequential experiments were performed as detailed below.

Experiment I: Phosphorus Flux under Ambient Floodwater Conditions
Floodwater samples of approximately 10 mL were removed from each core at 0, 1, 3, 5, 7, and 10 d, immediately filtered, acidified, and analyzed for dissolved reactive phosphorus (DRP). This resulted in the removal of approximately 60 mL of floodwater over the course of the experiment. In this experiment and all of the following, evaporative and sampling losses were measured over the 10-d experimental period and were accounted for in the determination of mass flux. The average incubation temperature was 28.5°C (±2.6).

Experiment II: Consumption of Electron Acceptors and Phosphorus Release
Following the conclusion of Experiment I, the floodwater was mechanically aerated to obtain an initial dissolved O2 content of approximately 6 mg L-1. The cores were then sealed with polyethylene piston-type closures. Dissolved O2 was then measured in each replicate core by inserting a dissolved O2 bottle electrode through a port in the piston. The volume to surface area ratio of 100 (L to m2) was sufficient to allow the O2 measurements to be made over a short time interval (within 4 to 12 h) (Bowman and Delfino, 1980). Oxygen consumption measurements were made in groups of three due to equipment limitations; however, the whole experiment was completed within 3 d. Even though there was a slight offset in time between the dissolved O2 determinations on the first and last sets of cores, this effect was minimal, since all environmental conditions were held constant.

Immediately following O2 consumption measurements, consumption of O2, NO3–N, and SO4–S were measured in the same cores. The cores were opened, aerated, and spiked with 900 mg of KNO3 to obtain a final NO3–N concentration of 9 mg L-1. This is well above ambient concentrations and thus the measured rates represent a potential maximum estimate of NO3 reduction under NO3 nonlimiting conditions. Ambient water column SO4–S in the cores was 34 (±16) mg L-1 at the initiation of the experiment, so spiking was not necessary. Oxygen consumption rates determined in this experiment are not reported, as only two measurements were obtained and the data from the continuous measurements described previously were believed more accurate. Floodwater samples (27 mL) were withdrawn for analysis at 0, 0.3, 1, 2, 4, 7, and 10 d. During the experimental period, floodwater dissolved O2 concentrations were reduced to <0.1 mg L-1, thus measured P flux represents conditions of O2 depletion and influence of alternate electron acceptors. The average incubation temperature was 25.5°C (±1.5).

In a separate experiment, soil oxygen demand (SOD) was measured on surface soil (0–10 cm) collected from Field Stations 1 through 8 along a nutrient enrichment gradient in WCA-2A (DeBusk et al., 1994) (Fig. 1). Approximately 5 g (wet weight) of soil was placed into a 250-mL biochemical oxygen demand (BOD) bottle placed on a magnetic stirrer. The BOD bottle was filled with deionized water and aerated for 5 min. A dissolved O2 probe was inserted into the bottle (YSI 5905 BOD probe, Yellow Springs Instruments, Yellow Springs, OH); the contents of the BOD bottle were continuously stirred during the 8 h incubation. Dissolved O2 measurements were made at 0, 0.5, 1, 2, 4, 6, and 8 h after initiation of the experiment. Incubations were made in the dark at 25°C. All measurements were made using three-replicate samples. The SOD rate was calculated as follows:

where KSOD is maximum linear slope of the plot between the dissolved oxygen (DO) concentration (mg L-1) and time (h), V is the volume of water in the BOD bottle (L), m is the mass of dry weight of soil (kg), and SOD is soil oxygen demand (mg kg-1 h-1).

Experiment III: Phosphorus Flux under Aerobic Floodwater Conditions
At the end of the Experiment II, the overlying floodwater from each core was siphoned off and immediately replaced with water collected from a nutrient-unaffected area of WCA-2A. This reflood water was coarse-filtered through Whatman qualitative #4 filters. The DRP concentration of the reflood water was 6 µg L-1. A 10-mL sample was withdrawn for analysis at 0, 7, 14, 21, 30, and 36 d. The approximate incubation temperature was 20°C. The conditions of this experiment were similar to Experiment I, except that the period of anoxia imposed in Experiment II presumably lead to the death of invertebrate populations. Thus, flux in this experiment represents nutrient exchange without the augmentation of solute transport due to invertebrates.

Experiment IV: Influence of Water-Level Drawdown and Reflooding on Phosphorus Flux
The soil cores from the previous experiment were then uncovered, drained of all floodwater, and exposed to air-drying for approximately 2 mo. Since the soil cores were relatively short (approx. 30 cm), all drying that occurred was a result of evaporation from the soil surface (i.e., no subsurface drawdown was used to expedite the drying process). Platinum-tipped redox electrodes were installed at a depth of 5 cm in one replication from each of the eight stations to approximately determine the extent of drying and oxidation. On average, the redox potential (Eh) increased from -29 (±237) to 215 (±81) mV during the drying period. After drying, soil cores were reflooded with filtered (Whatman #42) water from the nutrient-unaffected area within WCA-2A. The floodwater volume was adjusted such that each core contained 1 L (6 cm) of overlying water. The average incubation temperature was 26.0°C (±2.2).

In addition to the above described study, several intact soil cores measuring 10 cm in diameter were collected from Stations F1, F4, and U3 on 13 Mar. 1997 after a prolonged dry period in the field (Fig. 1). At the time of collection, the cores were devoid of floodwater. The cores were reflooded with 0.95 L of water collected from the field site on 20 Mar. 1997, resulting in a floodwater depth of 12 cm. The cores were incubated in the dark at 22°C. A 10-mL sample was withdrawn at 0, 1, 2, 4, 8, 15, 22, 29, and 39 d and analyzed for DRP, as described previously. Details concerning the sequence of experiments performed on the intact soil cores and the following field experiments are given in Fig. 2 .



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 2. Flowchart depicting the sequence of intact soil core experiments in the greenhouse and the field studies

 
Field Experiments
Benthic Chambers: Field Measurement of the Flux of Nutrients and Electron Acceptors
Benthic chambers placed at the soil–floodwater interface were used to determine in situ fluxes of dissolved O2, NO3–N, SO4–S, and P. The chambers were constructed of 6.35-mm (0.25 in)-thick acrylic and enclosed a soil surface area of 0.5 m2. Each of the chambers was equipped with a recirculation pump and a port for a dissolved O2 electrode. The approximate floodwater volume enclosed by the chamber was 80 L.

Three replicate chambers were placed in WCA-2A at Stations 1 (1.4 km from inflow), 3 (3.3 km from inflow), and 8 (10.2 km from the inflow) on 16 July 1996 and 4 Aug. 1997 (Fig. 1). The floodwater inside the chambers was spiked with 80 mg of N as KNO3 to obtain a floodwater concentration of approximately 1 mg L-1 NO3–N for the 1996 event. For the 1997 benthic chamber experiment, the chambers were spiked with 8 mg of N as KNO3 for a final water column concentration of 0.1 mg L-1 NO3–N. This concentration more closely approximated ambient NO3–N levels. Sulfate spiking was not needed, as the levels in the water column during the chamber experiment were in the range of 17 (±2) mg L-1. Dissolved O2, NO3–N, SO4–S, temperature, and pH were measured for a period of 24 h. The water surrounding the chambers was monitored for the same analytes as mentioned above. A Hydrolab (Austin, TX) Datasonde 3 was installed at each of the three stations to record ambient temperature, dissolved O2, and pH concurrent with the benthic chamber measurements. The average in situ temperature during the July 1996 deployment was 30.5°C (±1.1). During August 1997, the water column inside each chamber was mechanically aerated to obtain an initial dissolved O2 concentration of approximately 1.5 mg L-1 at Stations 1 and 3. Dissolved O2 levels at Station 8 were initially high, therefore aeration was not necessary. The average temperature during the August 1997 experiment was 28.7°C (±1.6). The recirculation pump was operated for approximately 1 min prior to each O2 determination in order to have sufficient water velocity for the oxygen electrode. Visual inspections were made to insure that sediment resuspension did not occur during pump operation.

Electron acceptor consumption and P flux measured in the soil cores and benthic chambers were calculated by determining the slope of the concentration vs. time curve through linear regression, then multiplying by the floodwater volume to soil surface area ratio of the soil core:

where JI = flux of component i (mg m-2 d-1), C = component concentration in floodwater (mg L-1), V = floodwater volume (L), A = soil surface area (m2), and t = time (d).

Porewater Equilibrators
The porewater equilibrators used in this study were patterned after a device described by Hesslein (1976). Porewater equilibrators consist of 2- x 10- x 50-cm blocks of acrylic into which are milled 8-cm3 cells that are spaced vertically 1 cm apart. The sample cells were filled with deionized water and overlain with 0.2-µm pore diameter polyethersulfone membranes (Gelman Sciences, Ann Arbor, MI, Product no. XE 22061). The equilibrators were placed into acrylic containers, sealed, and purged of O2 with N2 gas. In the field, the equilibrators were removed from the containers and pushed 30 to 40 cm into the soil. They were then left in situ for 2 wk, allowing time for dissolved constituents in the sediment porewater to equilibrate with the deionized water inside the cells. Equilibrium studies have shown that 2 wk is sufficient time for porewater nutrients to equilibrate with the solution inside the equilibrator cells in sediments of high porosity (Carignan, 1984; Carignan et al., 1985). Diffusion modeling has also indicated that the time necessary to achieve 90% equilibration with porewater dissolved P, through molecular diffusion alone, is approximately 3 d (Webster et al., 1998). If the assumption is made that some resupply of soluble P occurs through exchange with the solid phase, 99% equilibration occurs in approximately 2 d (Harper et al., 1997). Since the effects of macrobenthos would likely shorten the equilibration period, 2 wk was deemed adequate for equilibration. After the equilibration period, the equilibrators were withdrawn from the soil and the cells were sampled by withdrawing the contents with a syringe. Samples were withdrawn from the cells at 1-cm intervals to a depth of 10 cm below the soil–water interface, and thereafter every 2 cm. The samples were stored at 4°C until analysis.

To determine temporal variability in porewater concentration gradients, porewater equilibrators were installed in WCA-2A at Stations 1 through 8 on 7 Feb., 8 May, and 11 Sept. 1996.

The concentration gradients were used to estimate the flux of P across the soil–water interface using Fick's First Law:

where J = diffusive flux of component i (mg m-2 d-1), {phi} = soil porosity (cm3 cm-3), Ds = bulk soil diffusion coefficient (cm2 s-1), {partial}C/{partial}Z = concentration gradient of component i with respect to depth, Z (cm), and 8.64 x 105 is a units conversion.

The average pH in the upper 10 cm of soil during the February 1996 equilibrator sampling at Stations 1 through 8 was 7.05 (±0.08). This was used to determine the diffusion coefficient of P in water. Since the speciation of soluble P is pH dependent, an interpolated value of Ds was used. Li and Gregory (1974) reported the diffusion coefficients of HPO2-4 and H2PO-4 in pure water as 7.34 x 10-6 and 8.46 x 10-6 cm2 s-1, respectively. An average of these two values was used in the calculation of diffusive flux, 7.9 x 10-6 cm2 s-1. The diffusion coefficient was modified for the restrictive effect of soil structure by dividing it by the square of the soil tortuosity. Soil tortuosity ({Theta}) was calculated from a relationship developed by Sweerts et al. (1991), or:

where {phi} is the soil porosity.

Analytical Methods
Dissolved O2 was measured with a YSI model 58 O2 meter (Yellow Springs Instrument Company) equipped with a YSI model 5730 stirring electrode using USEPA Method 360.1 (USEPA, 1979). Dissolved O2 and temperature were continuously recorded in the greenhouse experiments with a Campbell Scientific (Logan, UT) Model CR10 datalogger. Nitrate N and SO4–S were determined with a Dionex (Sunnyvale, CA) Series 4500i ion chromatograph, using EPA Method 300.0 (USEPA, 1979). Floodwater samples were analyzed for DRP using a Technicon (Tarrytown, NY) AutoAnalyzer, and EPA Method 365.1 (USEPA, 1979). The pH and Eh were measured with a Fisher Scientific (Pittsburgh, PA) Accumet Model 1002 pH meter.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Greenhouse Experiments
Experiment I: Phosphorus Flux under Ambient Floodwater Conditions
Soils closer to the inflow source of agricultural wastewater (<3.25 km) showed appreciable flux of P to the overlying floodwater. The floodwater DRP concentration in cores from Station 1 increased from approximately 25 to 660 µg L-1 during the 10-d incubation period (Fig. 3) . The concentration in cores from Station 2 were similar; starting out at 25 µg L-1 and increasing to an average of 729 µg P L-1, while the cores from Station 3 increased from 27 to 213 µg L-1. The floodwater concentration in the cores from the remaining stations tended to remain between 5 and 10 µg L-1. This resulted in a P flux of 6.5 mg P m-2 d-1 from soil cores at Stations 1 and 2 to 1.5 mg m-2 d-1 at Station 3 (Table 1). Net P flux from soil cores from areas >4 km from inflow was negligible. The average pH and dissolved O2 in the floodwater were 7.6 (±0.2) and 1.35 mg L-1 (±0.78), respectively. Dissolved O2 was low because the floodwater was not mechanically aerated.



View larger version (36K):
[in this window]
[in a new window]
 
Fig. 3. Phosphorus flux observed in intact soils cores from WCA-2A under aerobic and anaerobic water column conditions (n = 3). Error bars represent one standard deviation. Time for Experiments I, II, and III: 0 to 10, 10 to 20, and 20 to 56 d, respectively. Note different axis scaling in panels (c) and (d)

 
Experiment II: Effects of Consumption of Electron Acceptors on Phosphorus Release from Soils
In addition to P flux measurements, consumption rates of O2, NO3–N, and SO4–S were also determined in this experiment. Oxygen consumption rates were higher in soil cores representative of stations close to the S10-C inflow structure, and decreased exponentially with distance from the inflow (Fig. 4) (Table 2). The O2 consumption rate at Station 1 was approximately 3.6 g m-2 d-1 and decreased to a constant 1.3 g m-2 d-1 at a distance of 4 km from the inflow. In order to distinguish respiration by organisms in the water column from soil O2 consumption, a separate 300-mL floodwater sample was incubated in the dark during the incubation period. Dissolved O2 consumption in the water column was found to be minimal (data not shown).



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 4. Dissolved oxygen, nitrate, and sulfate consumption rates as a function of distance from the nutrient inflow in WCA-2A as measured during Experiment II. Error bars represent one standard deviation

 

View this table:
[in this window]
[in a new window]
 
Table 2. A comparison of the P flux observed in three sequential intact soil core experiments. ND = not detectable. Rates reported represent the mean of three replicates and values in parentheses represent one standard deviation

 
In a separate batch of soil samples, SOD rates were measured under nonlimiting diffusion conditions (Fig. 5) . The SOD rates decreased with distance from inflow, with high O2 demand in nutrient enriched (affected) soils, as compared with unaffected sites. These SOD values represent 1.9 to 3.1 g O2 m-2 d-1 for a soil depth of 10 cm and bulk density of 0.1 g cm-3. These values are approximately in the same range as those observed in intact cores, suggesting O2 diffusion from water to soil is not a limiting factor in controlling overall SOD.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 5. Soil oxygen demand measured on surface soils collected at WCA-2A

 
All of the added NO3–N was consumed within 24 h. The NO3–N concentration decreased to undetectable levels (<5 µg L-1) within 24 h. Rates of NO3–N and SO4–S consumption were similar for soil cores obtained from all eight stations (Fig. 4b and 4c). There was no apparent effect of the WCA-2A nutrient enrichment gradient on soil NO3–N or SO4–S consumption rates. The average NO3–N and SO4–S consumption rates were 120 (±81) and 130 (±52) mg m-2 d-1, respectively.

Experiment III: Phosphorus Flux Under Aerobic Floodwater Conditions
Experiment III (aerobic P flux following anoxic electron acceptor experiment) was performed to determine if the initial DRP flux observed in Experiment I could be maintained at the same intensity and if asphyxiation of aerobic soil biota had any effect on P flux. The concentration of DRP in the floodwater of the Station 1 cores increased from approximately 172 to 332 µg L-1 over the 36-d experimental period, while cores taken from Station 2 increased from an average of 524 to 1196 µg L-1. The remaining stations showed very little or no increase in DRP. In general, all stations showed much lower P flux after reflooding, compared with the two previous experiments (Table 2).

Experiment IV: Influence of Water-Level Drawdown and Reflooding on Phosphorus Flux
At the end of the 60-d flooding period (Experiments I, II, and III), floodwater was removed and soil cores were allowed to dry for approximately 2 mo. This was followed by flooding for an additional period of 60 d. Dissolved P flux from these cores was lower in soil cores obtained from the affected area (2.3 km from inflow), as compared with flux measured in cores immediately after removal from the field (Experiments I and II) (Table 3).


View this table:
[in this window]
[in a new window]
 
Table 3. Influence of prolonged water-level drawdown and reflooding of soil cores obtained from WCA-2A. Soil cores from stations 1 through 8 were allowed to undergo drying when greenhouse conditions for a period of 60 d. Soil cores from Stations F-1, F-4, and U-3 were obtained from the field after prolonged drought. Rates reported represent the mean of three replicates and values in parentheses represent one standard deviation

 
After a prolonged dry period, additional soil cores were obtained from Stations F1, F4, and U3 (Fig. 1) and reflooded in the lab. The DRP flux was rapid at Station F1 (1.8 km from inflow S10-C), as compared with the other two stations sampled. The DRP concentration increased from 5 to 1200 µg L-1 within 2 wk after flooding (Table 3). The calculated P flux for Station F1 was in the same range as the flux observed for Station 1, with average flux of 5.2 (±4.3) mg P m-2 d-1. The flux from the soil cores retrieved F4 and U3 was near the limit of detection.

Field Experiments
Benthic Chambers
The mean initial water column DRP concentration at the initiation of the 1996 chamber experiment was 60 µg L-1, while the water column DRP concentration at the beginning of the 1997 experiment was considerably lower at approximately 15 µg L-1. The DRP concentration at Station 8 (10.1 km from inflow) generally remained below the limit of detection (ca. 2 µg L-1) throughout the July 1996 and the August 1997 experimental periods.

The average P flux from Stations 1 and 3 for the July 1996 experiment was 10 and 9 mg m-2 d-1, respectively, and was calculated using the first 4 h of data for Station 1 and the first 8 h of data for Station 3. Station 1 was the only station that showed any measurable P flux during the August 1997 experiment. Flux from this station was very similar to that measured in the 1996 experiment, 9.8 (±3.7) mg P m-2 d-1 (Table 4).


View this table:
[in this window]
[in a new window]
 
Table 4. Oxygen; nitrate N and dissolved reactive phosphorus (DRP) flux in benthic chambers in WCA-2A on 16 July 1996 and 4 Aug. 1997. Rates reported represent the mean of three replicates and values in parentheses represent one standard deviation. ND = not detectable; NA = not measured

 
The ambient dissolved O2 concentrations at Stations 1 and 3 were very low (ca. 0.5 mg L-1) at the initiation of the July 1996 experiment, probably due to extensive coverage of the water surface with Lemna spp. and dense growth of Typha spp. This precluded obtaining accurate estimations of O2 flux, as the initial concentration inside the chambers was only slightly greater than the sensitivity of the polarographic O2 electrode technique. An estimate of O2 consumption was successfully made at Station 8, as average initial dissolved O2 inside the chambers was 3 mg L-1. The average estimated O2 flux to the soil at Station 8 was 818 (±165) mg m-2 d-1 (Table 4).

The ambient dissolved O2 concentration was also low during the August 1997 experiment. However, an attempt was made to aerate the water column that succeeded in elevating the water column dissolved O2 to approximately 1.5 mg L-1. The consumption rates measured at all three stations were on the order of 100 mg m-2 d-1 using this approach.

Nitrate removal was rapid at Station 1 with almost complete disappearance within 8 h. The calculated NO3–N consumption rate was 274 mg m-2 d-1. Nitrate consumption at Station 3 was slower than at Station 1, with nearly all the NO3–N consumed in 24 h (Table 4). However, in the nutrient unaffected site (Station 8) approximately 39% of NO3–N added was still in the water column after the 24 h. Nitrate consumption in August 1997 was highest at Station 1 (124 mg m-2 d-1). Rates observed in the 1997 chamber experiment were approximately half those observed in 1996, probably due to the lower initial concentration used in the 1997 experiment.

There was very little change in SO4–S concentration inside the chambers during the 24 h period (data not shown). Ambient dissolved O2, pH, and temperature at each of the three stations were continuously recorded during the 24-h period that the benthic chambers were in place (Fig. 4c). Diel cycling in dissolved O2 and pH indicates increased primary productivity in the water column at Station 8, compared with Stations 1 and 3 (Fig. 6) .



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 6. Ambient diel changes in pH, dissolved O2, and temperature at Stations 1, 3, and 8 on 16 July 1996 in WCA-2A

 
Porewater Equilibrators
Phosphorus flux calculated from porewater DRP concentration decreased during the February through September 1996 sampling period (Table 5). This may be related to the hydraulic loading from inflow structure S-10C, a major supplier of P and water to WCA-2A. There was 57% reduction in water delivery rates from S-10C for the period January 1994 through December 1996 (South Florida Water Management District, 1997). Declining porewater P gradients may therefore be related to decreasing hydraulic and nutrient loading in this region of WCA-2A. There was little or no floodwater at Stations 1 through 4 for all but the February sampling.


View this table:
[in this window]
[in a new window]
 
Table 5. Seasonal P flux from selected stations along a nutrient enrichment gradient in WCA-2A, as measured with porewater equilibrators. Rates reported represent the mean of three replicates and values in parentheses represent one standard deviation. If absent, flux determined from one equilibrator. N/A = floating flocculent sediment obscured interface preventing flux calculation

 
Generally, P flux was highest at stations nearest the inflow and lowest at the unaffected stations (Table 5). Phosphorus flux calculated using the February equilibrators (floodwater present) was not correlated to dissolved O2 consumption, as measured in Experiment II. Very little or no flux of DRP was suggested by the porewater concentration gradients from soils located at >5 km from inflow.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Phosphorus flux from soils to the overlying water column was measured under various redox conditions, including water-table drawdown and supply of alternate electron acceptors. Phosphorus flux observed in the intact cores was correlated to soil O2 consumption rate (r = 0.75; p < 0.005), possibly indicating that greater microbial activity near the inflow region of WCA-2A was responsible for a higher rate of mineralization of newly accreted organic material. DeBusk and Reddy (1998) also observed greater carbon mineralization rates in regions closer to nutrient inflows in WCA-2A. Based on the results from the intact soil cores, there was an approximate stoichiometric ratio of P mineralized to O2 consumed of 0.002 (Fig. 7) . In the anaerobic core flux experiment (Experiment II), P flux was actually lower at Stations 1 and 2 and higher at Stations 3, 4, and 5 than in the previous aerobic experiment. There was no clear trend in the uptake of the other electron acceptors (NO3–N and SO4–S), either with respect to the distance from the source of nutrient loading or with respect to P flux. In both the benthic chamber experiment and soil cores experiments, NO3–N was added in excess of ambient concentrations. An anticipated result was the increased mineralization of P and therefore increased P flux. However, either the time scale of mineralization was too long to observe using this experimental approach or the existing soil organic matter was too humified for any further decomposition to take place (DeBusk and Reddy, 1998).



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 7. Stoichiometric relationship between oxygen consumption and P mineralization. Data are from intact cores collected from Stations 1 through 4 along the WCA-2A nutrient gradient

 
Oxygen uptake was on the order of 1.5 to 3.5 g O2 m-2 d-1 in the soil cores, yet only 0.1 to 0.8 g O2 m-2 d-1 in the in situ benthic chambers. Soil oxygen demand measurements are dependent on temperature, overlying water velocity, and initial dissolved O2 concentration (Bowman and Delfino, 1980). Temperature differences between the two experiments were approximately 2°C and therefore could not explain the variability between the results of the two methods. Rate of stirring can also have an effect on SOD due to its influence on the diffusive boundary layer (Sweerts et al., 1991). In the intact soil core experiment, the overlying water column was continuously stirred at a slow rate with the BOD bottle electrodes used. The water column in the in situ chambers was only intermittently stirred for several minutes during each O2 measurement. Though both experiments were stirred, actual water velocities and duration of stirring were probably different. Part of the difference in results may be due to differences in the initial dissolved O2 concentration (Edburg and Hofsten, 1973). It is difficult to control the dissolved O2 content under in situ conditions, though this was attempted with partial success in the August 1997 deployment. The closest agreement between the two techniques was at Station 8. Oxygen flux to the sediment in the soil cores was 1.5 (±0.1) g O2 m-2 d-1 versus 0.8 (±0.2) g O2 m-2 d-1 in the chambers, as measured during the July 1996 deployment. The closer agreement was probably due to the higher ambient initial dissolved O2 content (approximately 3 mg L-1) at Station 8 during the July 1996 chamber experiment. It was possible to elevate the initial dissolved O2 content in the soil cores to approximately 6 mg L-1 at the initiation of the experiment; therefore, the more realistic estimates of SOD are probably those estimated in the intact core experiment. Belanger et al. (1989) estimated rates of sediment O2 consumption at affected and unaffected regions of WCA-2A using a similar intact soil core technique and found results similar to those reported here. Dissolved O2 uptake at nutrient-affected and unaffected stations in 1985 were reported as 2.34 and 2.09 g O2 m-2 d-1, respectively.

Dissolved P flux measured using porewater P concentration gradients was approximately an order of magnitude lower than P flux measured using benthic chambers and intact soil cores. Many researchers have performed comparative studies similar to this one in lakes and estuaries and have also noted disparities between calculated and observed flux (Callender and Hammond, 1982; Devol, 1987; Hopkinson, 1987); however, prior to this study, these comparisons had not been made in wetlands. For instance, Gomez-Parra and Forja (1993) compared flux based on sediment porewater concentrations gradients to benthic chamber P flux in the coastal waters of the southwest of Cadiz, Spain and found that in all cases the benthic chamber flux exceeded the gradient-derived flux, sometimes by as much as 29 times.

Though not directly measured in our study, the advection of porefluids out of the sediment due to the activity of invertebrates can result in a flux at least equal in magnitude to molecular diffusion alone (McCaffrey et al., 1980; Van Rees et al., 1996). In experiments conducted with labeled tracers, McCaffrey et al. (1980) were able to distinguish flux due to concentration gradients from flux caused by bioturbation. The combined flux compared favorably with flux measured with chambers alone. Callender and Hammond (1982) compared the calculated P flux, based on concentration gradients, with flux measurements in benthic chambers and found a significant "flux enhancement" or ratio of in situ to calculated flux. They concluded that nutrient fluxes measured in situ with benthic chambers may be 1 to 10 times higher than flux calculated from porewater nutrient profiles and the increased fluxes were attributed to irrigation of sediments by macrofauna. The authors found that field locations that had the lowest concentrations of macrofauna also had the best agreement between calculated and measured flux. Van Rees et al. (1996) investigated the effect of benthic invertebrates on diffusion coefficients in Lake Okeechobee, Florida and found that the diffusion coefficient varied by as much as 470% within a single sediment type. They attributed this range to the presence of benthic fauna. Several other factors, including temperature, dissolved O2 content of floodwater, rate of stirring, and presence of biota have also been shown to regulate P flux from sediments to the overlying water column (Holdren and Armstrong, 1980; Moore et al., 1998; Khalid and Patrick, 1974; Olila and Reddy, 1997).

Other factors that are difficult to estimate in deterministic models of P exchange with overlying floodwater include microgradients at the sediment–water interface (Berner, 1980; Krom et al., 1994; Davison et al., 1991) and variability in sediment physical characteristics. However, Sweerts et al. (1991) found that the ratio of the diffusion coefficient in water to the diffusion coefficient in sediment varied little between sediments of low and high porosity, indicating that this is a minor concern.

In Experiment II, the floodwater dissolved O2 was completely depleted, which probably decreased the activity of benthic invertebrates, thereby decreasing solute transport due to bioturbation. Rutgers et al. (1984) described a technique to discriminate between the flux of silica caused by molecular diffusion and the flux augmented by bioturbation. They found a dramatic decrease in silica flux as the dissolved oxygen content within the benthic chambers approached zero, suggesting a reduction in biotic activity. If the differences in the calculated and the observed P flux in the WCA-2A study were due to activity of the benthic organisms, then it should be possible to modify the diffusion coefficient to account for their activity.

This can be expressed as:

where Ja = flux observed in soil cores under aerobic conditions, Jx = flux observed in soil cores under anoxic conditions, De = effective diffusion coefficient, and Ds = diffusion coefficient corrected for soil porosity only.

Knowing the values for the aerobic and anaerobic flux and Ds allows the calculation of De for use in subsequent diffusion studies. Since the average ratio between the aerobic and anaerobic (asphyxiated) flux of P was approximately 7, a more realistic value for the soil diffusion coefficient (Ds) in these sediments may be closer to 3.78 x 10-5 cm2 s-1. This value considers both the tortuosity effects of soil structure and the increased solute transport presumably caused by benthic organisms. Further evidence of the contribution of benthic organisms to nutrient transport was found in Experiment IV. The anticipated result of soil drying was the increased mineralization of stored nutrients, resulting in greater P flux upon reflooding. That this was not observed possibly indicates a negative effect of soil drying on the local benthic community.

Porewater concentration of DRP, as measured with the equilibrators, was generally highest during February and lowest in September, with the exception of Station 5 (5.12 km from S-10C), which showed a reverse trend. This is in contrast to porewater profiles observed in WCA-2A in 1990 (Koch-Rose et al., 1994). In the 1990 study, porewater nutrient profiles were generally higher in July and lowest in May. The authors attributed this seasonal sequence to high plant uptake of porewater nutrients (and therefore depressed concentration gradients) in spring and high mineralization rates of organic matter in late summer, leading to elevated porewater nutrient profiles. Sinke et al. (1990) found a high correlation between porewater CH4, DRP, and NH4–N and also concluded that this is evidence that mineralization of native Netherlands peat soils determined porewater concentrations of these nutrients. Stations 6 through 8 showed no distinct seasonal trend in porewater P and remained at or near the detection limit.

Results presented in this study showed that the water column redox conditions and water-level drawdown can significantly influence P flux to the water column. This internal load of P becomes a critical factor in regulating eutrophication status of the wetland, once external loads are curtailed. Thus, identification of internal P loading is essential to determine the time required for recovery. High P flux (expressed as an average of all conditions) from soils affected by P loading suggests that these soils function as a source of P to the overlying water column (Fig. 8) . If we assume that approximately 25% of the total P in the top 30 cm of soil (Reddy et al., 1998) is potentially mobile and can diffuse at a rate of approximately 2 mg P m-2 d-1 into the overlying water column, the measured P flux would be sustained for a period of approximately 5 yr.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 8. Average soluble phosphorus flux from soil to the overlying water column. Error bars represent standard error

 

    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Benthic chambers, intact soil cores, and porewater equilibrators were used in WCA-2A in order to estimate the flux of P. Phosphorus flux from the soil to the overlying water column determined in intact soil cores (Experiment I) was approximately equal to 6.5 mg m-2 d-1 at Stations 1 and 2 and 1.9 mg m-2 d-1 at Station 3. Phosphorus flux from the remaining stations was below the limit of detection. Phosphorus flux measured with the intact soil cores and the in situ benthic chambers gave similar results. Flux measured with the porewater equilibrators was approximately an order of magnitude lower than results observed from the other two techniques, possibly due to effects not generally included in deterministic models of solute flux.

There was a gradient in O2 consumption with respect to distance from the source of P loading in the intact soil cores. Oxygen consumption at Station 1 (1.4 km from inflow) was 3.6 g m-2 d-1 and decreased to a constant 1.5 g m-2 d-1 approximately 4 km from the inflow. There was no apparent trend in NO3–N and SO4–S consumption as measured in intact soil cores with respect to distance from the inflow. Average NO3–N and SO4–S consumption rates were 120 (±81) and 130 (±52) mg m-2 d-1, respectively. Soil oxygen demand was correlated to P flux, while NO3 and SO4 consumption rates were not, indicating the dominant role that aerobic processes have on P mineralization in wetland soils. Phosphorus flux from intact soil cores after the electron acceptor experiment (Experiment II) was considerably lower than P flux in soil cores maintained aerobically, possibly due to the death of invertebrate organisms. Draining, drying, and reflooding of intact soil cores resulted in P flux that was intermediate between the flux observed in anaerobically (Experiment II) and aerobically (Experiment III) incubated soil cores. Intact soil cores retrieved from Station F1 (1.8 km) after field-drying showed P flux that was similar to P flux measured in intact soil cores from Station 1, or 5.2 mg m-2 d-1, on average. The remaining soil cores retrieved from Stations F4 (6.8 km) and U3 (11.0 km) showed little or no P flux. Results of this study suggest that the soils underlying WCA-2A will function as a source of P loading to the water column, even after external loads are reduced.


    ACKNOWLEDGMENTS
 
Florida Agricultural Experiment Station Journal Series no. R-07819. This study was funded by the South Florida Water Management District. We would like to thank Yu Wang for laboratory assistance, Quentin Clark for field assistance, and three anonymous reviewers for their insightful comments.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Soil Sci.Home page
S. Grunwald and K. R. Reddy
Spatial Behavior of Phosphorus and Nitrogen in a Subtropical Wetland
Soil Sci. Soc. Am. J., June 18, 2008; 72(4): 1174 - 1183.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
G. C. Sigua, W.-J. Kang, and S. W. Coleman
Soil Profile Distribution of Phosphorus and Other Nutrients following Wetland Conversion to Beef Cattle Pasture
J. Environ. Qual., October 27, 2006; 35(6): 2374 - 2382.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
S. Grunwald, R. Corstanje, B. E. Weinrich, and K. R. Reddy
Spatial Patterns of Labile Forms of Phosphorus in a Subtropical Wetland
J. Environ. Qual., January 5, 2006; 35(1): 378 - 389.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
R. Corstanje and K. R. Reddy
Response of Biogeochemical Indicators to a Drawdown and Subsequent Reflood
J. Environ. Qual., November 1, 2004; 33(6): 2357 - 2366.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
L. M. Malecki, J. R. White, and K. R. Reddy
Nitrogen and Phosphorus Flux Rates from Sediment in the Lower St. Johns River Estuary
J. Environ. Qual., July 1, 2004; 33(4): 1545 - 1555.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
N. Vaillant, F. Monnet, H. Sallanon, A. Coudret, and A. Hitmi
Use of Commercial Plant Species in a Hydroponic System to Treat Domestic Wastewaters
J. Environ. Qual., March 1, 2004; 33(2): 695 - 702.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow