|
|
||||||||
a Smithsonian Environmental Research Center, P.O. Box 28, Edgewater, MD 21037
b Duke Univ. Dep. of Biology, Durham, NC 27708
c Univ. of Maryland, Dep. of Biological Resources Engineering, College Park, MD 20742
* Corresponding author (jordanth{at}si.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: TN, total nitrogen TNH+4, total ammonium TOC, total organic carbon TON, total organic nitrogen TOP, total organic phosphorus TPO3-4, total phosphate TSS, total suspended solids
| INTRODUCTION |
|---|
|
|
|---|
Preserving or restoring wetlands may help reduce nonpoint-source pollution. Wetlands can act as filters removing particulate material, as sinks accumulating nutrients, or as transformers converting nutrients to different forms, such as gaseous compounds of nitrogen (N) and carbon (C) (Richardson, 1989). Recent research has shown that constructed or restored wetlands can remove sediments and nutrients from nonpoint sources, including agricultural discharges (e.g., Fleischer et al., 1994; Mitsch, 1994; Raisin and Mitchell, 1995; Whigham, 1995; Jordan et al., 1999). Widespread restoration of wetlands has been suggested as part of a plan for reducing nitrogen releases from the Mississippi River basin (Mitsch et al., 2001).
Nutrient removal by constructed wetlands has been extensively studied for their use in wastewater treatment (Hammer, 1989; Kadlec and Knight, 1996). However, wetlands constructed for wastewater treatment usually receive measured and controlled inflows of wastewater. Also, the outflows from wastewater treatment wetlands are usually monitored to check the wetland's performance. Therefore, much is known about the capabilities and design criteria of such wetlands (Hammer, 1989; Kadlec and Knight, 1996). Much less is known about the nutrient and sediment removal capabilities of natural and restored wetlands that receive unregulated inflows. Unregulated flows are more difficult to measure than the regulated flows. Automated sampling is generally required to quantify unregulated event-driven fluxes (e.g., Kovacic et al., 2000; Braskerud, 2002).
In general, the ability of a wetland to trap or transform nutrients increases as the water retention time increases. Models incorporating the effects of water retention time are used in designing treatment wetlands (Kadlec and Knight, 1996). Similar effects of water retention in natural and restored wetlands have been suggested by several studies (e.g., Mitsch et al., 1995; Carleton et al., 2001). Water retention time may vary widely with weather and season in wetlands with unregulated inflows. Variability of water flow may diminish the ability of wetlands to remove nutrients and sediments, as removal capacities may be temporarily overwhelmed during short-lived high flow events (e.g., Kovacic et al., 2000).
Wetlands are being restored in agricultural watersheds to provide wildlife habitat as well as improve water quality (Whigham, 1995). Some restorations involve minimal alterations of drainage, which produce wetlands with highly variable inflow rates (e.g., Magner et al., 1995). The objective of this study was to quantify the removal of nutrients and sediments of one such wetland under a highly variable flow regime. This wetland, on the Chesapeake Bay shore, receives cropland runoff and may serve as a model for systems mitigating nonpoint-source nutrient discharges toward the goal of 40% reduction set by the Chesapeake Bay Program (1997). We hypothesized that the wetland would remove nutrients and sediments although the variable inflow rate would reduce the removal efficiencies compared with similar wetlands with more constant inflow rates. To test this hypothesis, we monitored fluxes of water, nutrients, and sediments into and out of the wetland for two years using an automated sampling system to permit observation of unpredictable episodes of high flow.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Water enters the wetland via drainage leads carrying surface runoff from the surrounding watershed and via precipitation directly on the wetland surface. Water leaves the wetland via the standpipe drain installed in the dike and via evapotranspiration. When the water was deep enough to flow out the drain, the entire 1.3-ha area of the wetland was submerged and lacked well-defined flow channels. Ground water exchanges are negligible due to the impermeable layer of clay within 0.5 m of the soil surface. We concluded that the clay layer blocks water infiltration, because clay sampled from beneath inundated areas was dry.
Measuring Water Flow
We used automated instruments to measure water flow and to sample water entering and leaving the wetland from 8 May 1995 through 12 May 1997. The instruments included a depth sensor consisting of a float and counter weight suspended in a stilling well that was connected to the impounded water near the wetland drain. A CR10 datalogger (Campbell Scientific, Logan, UT), housed atop the stilling well, recorded the position of the float to monitor water depth. Outflowing water passed over a 120° V-notch weir at the drainpipe. The outflow rate was calculated from the depth of water in the V-notch.
The total rate of water input to the wetland from runoff and direct precipitation combined was calculated by summing the rate of outflow and the rate of increase in water volume held in the wetland, with decrease in volume treated as negative increase. Thus, if water volume in the wetland remained constant, then the total water input rate was assumed to equal the outflow rate; if water volume increased, then the total input rate was assumed to equal the outflow rate plus the rate of increase in volume; and if the water volume decreased, then the total input rate was assumed to equal the outflow rate minus the decrease in volume. This method of calculating the total input of water yielded negative values for input when there was no surface water flow and evapotranspiration decreased the water volume in the wetland. Therefore, we interpreted negative values of total input as indicative of zero input from runoff and precipitation. We also assumed that evapotranspiration was negligible during periods of precipitation and runoff input. The volume of water in the wetland was calculated from water depth and the areas enclosed in 10-cm elevation contours within the wetland basin, which was surveyed with a Total Station CTS-2/2B (Topcon, Tokyo, Japan).
The water input from precipitation directly on the wetland surface was calculated from the surveyed wetland area (1.3 ha) and the precipitation volume measured with standard rain gauges at the wetland and at the Wye Research Center (WRC), 13 km from the wetland. The WRC precipitation data were obtained from the Maryland State Climatologist. The amount of water the wetland received from watershed runoff was calculated by subtracting the direct precipitation input from the measured total water input from runoff and precipitation combined. Evapotranspiration from the wetland was estimated using data from a standard weather-bureau evaporation pan at the Smithsonian Environmental Research Center (SERC), 25 km from the wetland.
Sampling Water
The datalogger controlled pumps that collected separate samples of inflowing runoff and outflowing water in volumes proportional to the respective flow rates. This produced volume-weighted composite samples that represented the water quality of the inflowing runoff and outflowing water. The logger calculated the rates of outflow and total inflow (the sum of runoff plus precipitation) every 15 min. For the purposes of controlling sample pumping, total inflow was assumed to be proportional to runoff. When the amount of outflow or total inflow since the last pumping exceeded a certain threshold, the logger activated the appropriate pump to collect a volume of water proportional to the amount of flow since the last pumping. Thus, the frequency of pumping and the amount of water pumped each time could both vary. During high flow events, the pumps could be activated as often as every 15 min. Increasing the frequency of pumping during high flow events was important because concentrations can change rapidly during runoff events. However, it was also necessary to vary the volume pumped each time because the flow threshold for triggering pumping could be exceeded by different amounts during different 15-min measurement cycles. The pumps sampling inflowing runoff and outflow were controlled independently since inflow and outflow usually differed in timing and rate. The logger recorded each time when each pump was activated and how long it was pumping for each sampling. The signal to activate the inflow pump was transmitted via wire from the logger to the inflow pump located near the drainage ditch that collects the runoff from about 70% of the wetland's watershed. Inflowing runoff was sampled from about 5 cm above the bottom of this ditch. Previous comparisons showed that runoff carried by the main ditch had similar chemical composition to runoff entering from two other points (Jordan et al., 1999). Outflowing water was pumped from the water column near the V-notch weir at the wetland drain (Fig. 1). A submersible impeller pump (Model 1P811A; Teel) was used to sample outflow but a self-priming peristaltic pump (Model LG100; Little Giant Pump Co., Oklahoma City, OK) was needed to sample inflow because the inlet ditch usually dried up between runoff events. At both the inflow and outflow sampling points, the samples were pumped through plastic tubing, which was first rinsed with stream water. The pumped sample stream was split between two carboys, one with about 3 mL of sulfuric acid per liter of sample added as a preservative and one without preservative. In addition to the automatically collected composite samples, grab samples were collected whenever there was water flowing at the inlet or outlet during the weekly visits to retrieve the composite samples.
Each week the composite samples that accumulated during the week were brought into the laboratory for analysis. The acid-preserved samples were analyzed for nutrient concentrations. The measured nutrient concentrations represent the total of dissolved nutrients plus particulate nutrients that were dissolved by the acid preservative. The unpreserved samples collected by the automated instruments were analyzed for total suspended solids (TSS). Unpreserved grab samples were analyzed for pH and conductivity immediately after return to the laboratory.
The content of N, P, and organic C in bulk precipitation was measured in an ongoing monitoring program at the Smithsonian Environmental Research Center (e.g., Correll et al., 1994; Jordan et al., 1995). After each event of more than 5 mm of precipitation, samples of bulk precipitation were collected with a 28-cm-diameter polyethylene funnel and bottle. These samples were analyzed by the same methods as for the wetland water samples.
Chemical Analyses
Standard techniques were used for analysis of nitrogen (N) and phosphorus (P) compounds. Samples for total inorganic nutrients (including originally dissolved species and those dissolved by the acid preservative) were filtered before analysis with prewashed 0.45-µm Millipore (Bedford, MA) filters. Total Kjeldahl nitrogen (TKN), total P, and total organic C were measured on unfiltered samples. Total P was digested to phosphate with perchloric acid (King, 1932). Phosphate in the digestate and total phosphate (TPO3-4) in undigested samples were analyzed by reaction with stannous chloride and ammonium molybdate (American Public Health Association, 1995). The TKN was digested with sulfuric acid, Hengar granules, and hydrogen peroxide (Martin, 1972). The resultant ammonia was distilled and analyzed with a Dionex (Sunnyvale, CA) ion chromatograph. In undigested aliquots, total ammonium (TNH+4) was oxidized to nitrite by alkaline hypochlorite (Strickland and Parsons, 1972), dissolved nitrate was reduced to nitrite by cadmium amalgam, and then the nitrite was analyzed by reaction with sulfanilamide (American Public Health Association, 1995). We present data on the sum of nitrite and nitrate concentrations, which we refer to as NO-3. From results of the above analyses we calculated total nitrogen (TN) by adding NO-3 to TKN, total organic nitrogen (TON) by subtracting TNH+4 from TKN, and total organic phosphorus (TOP) by subtracting TPO3-4 from total P.
Total organic carbon (TOC) was analyzed as chemical oxygen demand by drying samples at 60°C, followed by reaction with potassium dichromate in 67% sulfuric acid at 100°C for 3 h (Maciolek, 1962; American Public Health Association, 1995). Organic carbon was calculated from the amount of unreacted dichromate measured colorimetrically (Maciolek, 1962; Gaudy and Ramanathan, 1964).
Total suspended solids (TSS) were measured by filtering the nonacidified samples through prewashed, preweighed 0.40-µm Nuclepore filters (Whatman, Maidstone, UK), after which the filters were dried and reweighed.
Measurements of pH were made on air-equilibrated samples with an expanded range pH meter and a ROSS electrode (Thermo Orion, Beverly, MA). Conductivity was measured with a Model 32m conductivity meter (YSI, Yellow Springs, OH).
Statistical Analyses
To assess whether annual net fluxes were statistically significant, we calculated the 95% confidence limits around the annual net fluxes using the bootstrap technique (Efron, 1982). Differences among annual fluxes may arise due to the variability among weekly fluxes, because a few weeks with high flux can dominate the calculation of the annual flux. The bootstrap technique measures the consequences of randomly including or excluding certain weekly fluxes from the calculation of the annual flux. In other words, it accounts for the chance occurrence of weeks with differing water and nutrient flow within a given year. The bootstrap procedure begins by creating 1000 sets of data by randomly selecting data points from the original data set, replacing the selected points so they can be chosen again. In this case, the data points are the weekly net fluxes within the one-year or combined two-year periods. We included only weekly net fluxes for weeks when measurements were available for both influx and outflux. We did not include weeks for which net flux was estimated due to missing measurements. Each of the data sets created by the bootstrap procedure has the same number of samples as the original data set. The means of the created sets are calculated, and the 2.5 and 97.5 percentiles of these means represent the 95% confidence limits of the original mean (Efron, 1982). If the confidence limits of the net flux do not overlap zero, then the net flux is significantly different from zero at the p < 0.05 level. The bootstrap analysis and all other statistical analyses were performed using the Statistical Analysis System (SAS Institute, 1989).
| RESULTS |
|---|
|
|
|---|
|
|
Inflow to the wetland was episodic, depending almost entirely on rain events. Half of the total water inflow occurred in only 24 d scattered throughout the two-year study. During eight of those days, outflow exceeded 2500 m3, the water holding capacity of the wetland (Fig. 4) . Even on a weekly basis, total water inflow to the wetland was very uneven (Fig. 5) . More than half of the total annual inflow occurred during only 12 weeks scattered throughout the study.
|
|
Our estimates of annual water gains and losses in the wetland come within 5% of balancing (Table 1) . This close agreement supports our assumption that the underlying clay layer prevented ground water exchanges. Runoff from the watershed was the main source of water input. The total input from runoff and rainfall over the two years (150200 m3) was about 60 times the water holding capacity of the wetland (2500 m3). Surface outflow was the main water loss. Surface flows were especially dominant in the second year of the study, which had a wetter summer than the first year. Annual net change in standing water volume was orders of magnitude smaller than the surface flow. Our estimate of evapotranspiration, assumed equal to pan evaporation, is probably the least certain component of the water budget. However, pan evaporation (Table 1) was very similar to potential evapotranspiration estimated from the Thornthwaite equation (Veihmeyer, 1964), which was 770 mm yr-1 or 9995 m3 yr-1 for the whole wetland.
|
|
|
|
The concentrations of some materials were serially correlated (i.e., there were multiweek trends in concentrations). We assessed the serial correlation for each material from the correlation between its concentration in each week with the average of its concentrations in the weeks immediately before and after. The concentrations of all of the materials in outflow except TSS had significant serial correlation (Pearson, p < 0.05). Among outflow concentrations, NO-3 had the strongest serial correlation with an r2 of 0.76, TPO-34 and TOC had the next highest serial correlations with r2 values of 0.37 and 0.33, respectively, and correlations for other materials ranged down to r2 = 0.10 for TON. Materials in inflowing runoff showed less serial correlation than materials in outflowing water. Only two materials in inflowing runoff had significant serial correlation: TPO-34 and TOC with r2 values of only 0.18 and 0.29, respectively.
In some cases the concentration of a material in outflow was correlated with its concentration in inflowing runoff. Such correlations are more likely as the volume of outflow increases, indicating more rapid passage of water through the wetland. We analyzed correlations of runoff and outflow concentrations for weeks with outflow volumes of more than 500 m3 (i.e., more than 20% of the water holding capacity of the wetland). These represented 47% of the weeks studied and accounted for 97% of the outflow that occurred. During these weeks, concentrations in inflowing runoff and in outflowing water were significantly correlated (Pearson, p < 0.05) for NO-3, TPO-34, TSS, TNH+4, and TOC, with r2 values of 0.46, 0.25, 0.24, 0.18, and 0.18, respectively.
Correlations between concentrations and water flow rates might be expected due to the effects of flow on erosion, resuspension, or dilution. However, the only significant correlations (Pearson, p < 0.05) we observed were for TNH+4 and TOC in inflowing runoff, which were negatively correlated with inflow rate (r2 = 0.11 and 0.12, respectively). The weak negative correlations suggest a slight tendency for high water flows to dilute the materials. Relationships between concentration and flow rate may be difficult to demonstrate with our weekly data because weeks differ not only in the total amount of flow but also in how the flow is distributed during the week. For example, a week with only moderate total flow may include a short period of very high flow that could affect concentrations in the weekly composite samples. Therefore, event-based sampling might be needed to reveal correlations between concentrations and flow rates. However, conditions antecedent to the flow event, such as soil saturation, may also influence the effect of water flow rate and thereby obscure correlations between flow and concentration.
Concentrations of materials other than TNH+4 and NO-3 were usually much higher in inflowing runoff than in precipitation (Fig. 6 and 7). Therefore, precipitation falling directly onto the wetland surface dilutes most materials entering the wetland via runoff. This dilution is important, although only about one-fifth of the water entering the wetland enters via direct precipitation (Table 1). Rather than diluting TNH+4 and NO-3, direct precipitation represents a considerable source of those materials to the wetland because their concentrations were similar in precipitation and runoff (Fig. 8).
Fluxes of Nutrients and Suspended Solids
We calculated fluxes of materials based on weekly concentrations and water flows. Our concentration measurements represented most of the water flow volume but were more complete for outflow than for inflow. The percentage of outflow volume for which we measured concentrations was 90% for TSS, 96% for TOC, and 99% for forms of N and P. The percentage of inflowing runoff for which we measured concentrations was 69% for TSS and 74% for other materials.
When concentration data were missing, we substituted annual flow-weighted mean concentrations, which were calculated from all the available measurements. To avoid bias in weekly net flux calculations, we also substituted flow-weighted means for measured concentrations when concentrations in the opposing flow were not available. This usually applied to weeks when concentrations were measured in outflow but not in inflowing runoff. In such cases we would substitute mean concentrations for outflow as well as for runoff. We followed this procedure because runoff concentrations of NO-3, TPO-34, TSS, TNH+4, and TOC correlated with their respective outflow concentrations. Thus, substituting a mean concentration for only one of the flow directions could create an artificial imbalance of concentrations leading to a less accurate estimate of net flux.
The protocol for filling missing data was selected because it could be applied consistently for all the missing data and because alternate protocols would require predictions based on weak correlations. In some cases, missing concentrations of NO-3, TPO-34, TSS, TNH+4, and TOC could have been predicted from the correlations between their concentrations in runoff and outflow. However, those correlations had r2 values of 0.46 to 0.18 and would therefore provide a poor basis for making predictions. Moreover, a measured concentration was not always available to use for predicting a missing concentration in the opposing flow. Similarly, the weak serial correlations for outflow concentrations would yield only imprecise predictions of missing values and could be applied in only a limited number of cases.
The total annual influxes of nutrients and TSS differed between the two years of our study, but generally not as much as would be expected from the difference between the annual water flows. Influxes of TOP and NO-3 were actually greater in the first year than in the second year (Fig. 9) . Influxes of other nutrients were greater in the second year than in the first year (Fig. 9), but not by a factor of two, which was the difference between the annual water inflows from runoff (Table 1). Clearly, the increases in inflows of water from runoff in the second year were offset by decreases in concentrations of materials in the inflowing runoff. Direct atmospheric deposition provided <4% of the total annual influx of most materials but provided 28 to 33% of the influx of TNH+4 and 18 to 36% of the influx of NO-3 (Fig. 9). Unlike influxes, the annual outfluxes of TOP, TPO3-4, TON, and TOC were elevated in the second year to about the same extent as the outflow of water (Fig. 9, Table 1). However, this was not true for NO-3 and TSS, which flowed out of the wetland in lesser amounts during the second year than during the first (Fig. 9). In the first year, influxes of TOC and all forms of P and N exceeded outfluxes, suggesting a net removal of these materials by the wetland. In the second year, there were apparent net releases of TOP, TPO3-4, and TON from the wetland, while other materials appeared to be removed by the wetland.
|
|
We assessed the statistical significance of annual net fluxes by using the bootstrap technique (see Materials and Methods section) to measure the consequences of randomly including or excluding different weekly net fluxes from the calculation of annual net flux. The bootstrap analysis highlights the differences between the two years. In the first year, the net influxes of all forms of P, N, and organic C were statistically different from zero (Table 2) . In the second year, only the net influxes of NO3, TNH+4, and TOC, and the net release of TPO3-4, were statistically different from zero (Table 2). For the combined two-year period, only the net influxes of NO3, TNH+4, and TOC were statistically different from zero (Table 2).
|
When calculating net flux for weeks when concentration was measured in outflow but not in inflowing runoff, we substituted annual flow-weighted mean concentrations for both outflow and inflow. As mentioned, we think this produces the best estimate of net flux because it avoids possible biases that may arise due to correlations between concentrations in outflow and inflowing runoff. Another way to deal with missing concentration would be to only substitute mean concentrations when no concentration measurement was available. Annual net fluxes calculated that way (with minimal substitutions) generally agreed well with net fluxes calculated by our preferred method. For annual mass per area net fluxes that were statistically different from zero (p < 0.05, bootstrap) the differences between the alternate estimates were <5% except for three materials. Differences greater than 5% were as follows: compared with using minimal substitutions, our preferred method predicted 16% less net influx of TOC in the first year, and 13% less net influx of TPO-34 and 74% more net influx of TNH+4 in the second year. Calculated with minimal substitutions, the net influx of TNH+4 in the second year is 1.1 kg N ha-1 yr-1. By comparison, our preferred estimate is 2.0 kg N ha-1 yr-1 (Table 2). This is a small absolute difference compared with the net influxes of other forms of N (Table 2). Thus, the major conclusions of our study would not be changed by estimating net influxes with minimal substitutions.
| DISCUSSION |
|---|
|
|
|---|
Comparing Removal Rates
Different wetlands remove materials at widely differing rates (e.g., see reviews by Verhoeven and van der Toorn, 1990; Mitsch and Gosselink, 1993; Whigham, 1995; Mitsch et al., 2000). To put our results in the broadest possible context, we compared removal rates for our wetland with average rates for a wide variety of wetlands reviewed by Kadlec and Knight (1996). By this comparison, our wetland seems to remove nutrients and suspended sediments at below average rates (in mass per area), even in the year when nutrient removal was highest (Table 2). For example, in that year our wetland removed total P at about one-third the average rate and total N at about one-ninth the average rate (Table 2). For some forms of nutrients the differences were even greater (Table 2). However, if we compare the percentages of inflowing nutrients removed in the same year, our wetland does not seem very different from average (Table 2). This may reflect the fact that influxes of materials to our wetland are lower than average due to the large area of our wetland (1.3 ha) relative to the area of its watershed (14 ha).
The concentrations of materials in inflowing water may influence their rates of removal. Often, removal rates are modeled according to first-order kinetics, with removal rates proportional to concentration (Kadlec and Knight, 1996). This could explain why our wetland, which usually had <1 mg NO-3N L-1 in inflowing water (Fig. 8), removed NO-3 at a lower rate than did the wetlands studied by Kovacic et al. (2000) and Hunt et al. (1999), with 9 to 13 mg NO-3N L-1 and 3 to 9 mg NO-3N L-1 in inflowing water, respectively. The highest NO-3 concentrations entering our wetland followed the extended dry period from JuneOctober 1995 (Fig. 8). The highest concentration observed (5.3 mg NO-3N L-1) was in the first runoff event after the dry period. After that, inflowing NO-3 concentrations remained elevated for about 5 mo (Fig. 8). The antecedent dry conditions may have promoted NO-3 accumulation in the watershed soil by enhancing nitrification while preventing NO-3 removal via runoff or denitrification. The differences between concentrations of NO-3 in inflowing runoff and outflowing water tended to be greatest during weeks with higher concentrations in runoff (Fig. 8). This suggests that our wetland may have the capacity to remove NO-3 at higher rates if runoff entering the wetland had higher NO-3 concentrations. Jordan et al. (1997a) found that Delmarva watersheds with 80% cropland (the proportion in our wetland's watershed) typically discharge water with about 3 mg NO-3N L-1, 1.2 mg TON L-1, and 4 to 13 mg TOC L-1. By comparison, discharges from our wetland's watershed generally had <1 mg NO-3N L-1, 1 to 5 mg TON L-1, and 15 to 80 mg TOC L-1. These concentrations probably reflect the lack of ground water flow from the watershed draining into the wetland because NO-3 concentrations decrease and total organic N and C concentrations increase as the proportion of ground water in watershed discharge decreases (Jordan et al., 1997b).
Fluxes of particulate and dissolved materials are likely to differ, but, because of our acid preservative, we could only measure the combined total fluxes of particulate and dissolved materials except for NO-3, which is essentially all dissolved. A previous study analyzed dissolved and particulate materials separately in grab-sampled water flowing in and out of our study wetland (Jordan et al., 1999). That study found that inflowing TOC, TON, and TOP were 85, 35, and 15% dissolved matter, respectively. The differences in the proportions of dissolved matter may account for the differences among the temporal variations of TOC, TON, and TOP concentrations (Fig. 6).
Effects of Flow Variability
The low absolute rates of nutrient removal by our wetland (Table 2) may reflect the unregulated inflow. The removal rates reviewed by Kadlec and Knight (1996) are based on a diversity of surface flow wetlands, including constructed wetlands and natural wetlands, receiving water from a variety of municipal and agricultural sources. However, most of the wetlands reviewed had regulated inflows. The few published studies that give absolute rates of TN or TP removal from unregulated inflows (e.g., Table 3)
represent a wide variety of wetland types with a wide range of hydraulic loading rates. For example, the wetlands studied by Kovacic et al. (2000) received tile drain effluent leached from cropland soils. Thus, nitrate was the main form of N input and there was little input of particulate matter or P. These wetlands were effective at removing nitrate, probably via denitrification, but were less effective at removing organic N and ineffective at removing P (Kovacic et al., 2000). In contrast, the wetlands studied by Braskerud (2000)(2002) received stream water carrying agricultural runoff with high particulate loads. These wetlands were inefficient at removing nitrate but effective at removing organic N, particulate matter, and P. Similarly, a restored prairie pothole wetland in an agricultural watershed was effective at removing particulate matter and P (Magner et al., 1995). Usually P is associated with particulate matter but one of the wetlands studied by Reinelt and Horner (1995) received and removed relatively high amounts of dissolved phosphate from ground water flowing through P-rich deposits.
|
The effects of inflow variability have not been addressed by many studies. However, Raisin and Mitchell (1995) used automated samplers to measure mass balances of N and P during high flow events in three wetlands that receive agricultural runoff. They found that mass balances differed greatly among different high flow events with events in winter causing net releases due to flushing (Raisin and Mitchell, 1995). We could not resolve the effects of individual flow events because our samples were weekly composites. We found large differences in net flux among weeks, with net removal in some weeks and net export in others (Fig. 10). However, we could find no correlations between concentrations or net fluxes of materials and either the total weekly water inflow or the weekly maximum water inflow rate. Event-based automated sampling would probably be better than weekly automated sampling for revealing correlations between water flow and net fluxes of materials. Due to the importance of unpredictable high-flow events, automated sampling is essential for quantifying mass balances for wetlands with variable unregulated inflow. A few studies have used flow-proportional automated sampling (e.g., Reinelt and Horner, 1995; Kovacic et al., 2000; Braskerud, 2002), as we did, while others have used automated sampling at fixed time intervals (e.g., Magner et al., 1995; Hunt et al., 1999), or at selected inflow rates (Raisin and Mitchell, 1995; Raisin et al., 1997).
Carleton et al. (2001) reviewed studies of 49 wetlands receiving unregulated inputs of urban or agricultural runoff and concluded that the wetlands performed similarly to wetlands with regulated flow in removing pollutants. However, they noted that there was high variability of performance that could be related to the temporal variability of inflows. They summarized wetland performance by regressing the percentage of material influx removed versus the ratio of wetland area to watershed area. Such regressions should be interpreted cautiously because the regressed variables are both correlated with the amount of water inflow, which may cause spurious correlations as defined by Kenney (1982) and Garsd (1984). Comparing a few studies that report absolute as well as relative removal rates, we found that the percentage of N and P influx that is removed tends to increase as the hydraulic loading rate decreases and the detention time increases (Table 3). However, net exports of N or P sometimes occurred from wetlands with the lowest hydraulic loading rates (e.g., Table 3: our study and Kovacic et al., 2000). Moreover, the highest absolute rates of nutrient removal were reported for wetlands receiving the highest hydraulic loading rates (e.g., Table 3: Braskerud 2000, 2002; Fleischer et al., 1994).
One factor that can affect nutrient removal is hydraulic efficiency, the degree to which inflowing water is dispersed over the wetland area (Persson et al., 1999). The even dispersion of inflowing water over the wetland surface maximizes hydraulic efficiency and nutrient removal. In contrast, channeled flow may limit the exposure of inflowing water to the wetland surface and thereby limit nutrient removal. For example, the wetlands studied by Reinelt and Horner (1995)(Table 3) carry water in channels that are only 11 to 25% of the total wetland areas. In contrast, our wetland is completely submerged during periods when the water was deep enough to flow over the weir. Thus, there is the potential for the inflowing water to interact with the entire wetland surface before flowing out. However, it is likely that flow was not evenly dispersed over our wetland and, lacking measurements of dispersion, we do not know whether uneven flow may have limited nutrient removal. Persson et al. (1999) discuss design features that can maximize hydraulic efficiency of constructed wetlands.
For wetlands with unregulated variable inflow, nutrient removal may be improved by reducing the variability of the outflow. This could be achieved by designing outflow control structures, such as dikes and drains, to maximize the residence time of water within the wetland. For example, a drain that allows the wetland to slowly empty after a storm inflow prolongs water residence in the wetland by providing holding capacity for later storm inflows. In contrast, a flat-topped standpipe drain, which maintains nearly constant water volume in the wetland, makes the wetland unable to hold additional water from storm inflows.
Seasonal and Interannual Changes
Besides differing among high flow events, nutrient removal may vary at seasonal and interannual time scales. Several studies have observed seasonal changes in NO-3 removal presumably linked to the temperature dependence of denitrification (e.g., Hunt et al., 1999; Spieles and Mitsch, 2000). For our wetland, variability among high flow events may have obscured seasonal patterns. However, we did observe striking differences between the two years of the study, with the greatest N and P removal in the first year (Table 2). The two years also differed hydrologically, with drying period from JuneOctober in the first year only (Fig. 2). We cannot be certain that differences in nutrient removal are linked to the drying period, but the elevated concentrations of NO-3 in inflow after the dry period (Fig. 8) may have enhanced NO-3 removal, as we have discussed. Interannual differences in hydrology and nutrient removal were also noted for a restored prairie pothole wetland (Magner et al., 1995).
Besides interannual differences linked to rainfall, there may be long-term trends in nutrient trapping as the wetland ages. Our study was conducted about a decade after restoration of the wetland when emergent vegetation was well established. After wetland restoration, nutrient removal efficiency could improve as growing vegetation helps trap and hold sediment and produces organic matter to support denitrification (e.g., Mitsch and Carmichael, 1996; Mitsch et al., 2000; Spieles and Mitsch, 2000). During our study, plant biomass in our wetland was not increasing monotonically but varied from year to year in response to variations in precipitation (Whigham et al., 2002). Aboveground biomass and plant N and P were lower in the first year of our study (Whigham et al., 2002) when nutrient retention was highest. This suggests that plant biomass did not limit nutrient removal in our wetland.
Removal of water-borne N and organic C may continue indefinitely if these elements are converted to gaseous forms in the wetland and released to the atmosphere. However, removal of suspended sediment and P may cease sometime after wetland restoration as the wetland fills with sediment and the sediment becomes saturated with P (Richardson, 1989). One multiyear study of constructed wetlands documented an increase in sediment retention in the first four years as vegetation coverage increased (Braskerud, 2001), followed by a gradual decline in organic N removal from 3 to 10 yr after construction (Braskerud, 2002). In a recent review, Mitsch et al. (2000) suggest that sustainable removal rates range from about 5 to 50 kg ha-1 yr-1 for P and 100 to 400 kg ha-1 yr-1 for N (with 1000 to 2000 kg ha-1 yr-1 N removal possible in warm climates). By comparison, the two-year average removal rates for our wetland are near the low end of the sustainable range for P (7.6 kg P ha-1 yr-1, Table 2) and below the sustainable range for N (17 kg N ha-1 yr-1, Table 2). In contrast, some other constructed wetlands with unregulated inflow (Table 3) removed P (Braskerud, 2000) or N at rates above the sustainable range (Fleischer et al., 1994; Hunt et al., 1999; Braskerud, 2002).
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|