Next Article in Journal
Change Point Enhanced Anomaly Detection for IoT Time Series Data
Previous Article in Journal
Land-Use/-Cover Changes and Their Effect on Soil Erosion and River Suspended Sediment Load in Different Landscape Zones of European Russia during 1970–2017
Previous Article in Special Issue
The Past, Present, and Future of Phosphorus Removal Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Field Scale Demonstration of Fly Ash Amended Bioretention Cells for Stormwater Phosphorus Removal: A Review of 12 Years of Work

1
Oklahoma Water Survey, University of Oklahoma, Norman, OK 73019, USA
2
Bureau of Overseas Building Operations, US Department of State, Arlington, VA 22209, USA
3
Parsons Corp., Chicago, IL 60606, USA
4
Biosystems and Agricultural Engineering, Oklahoma State University, Stillwater, OK 74078, USA
*
Author to whom correspondence should be addressed.
Water 2021, 13(12), 1632; https://doi.org/10.3390/w13121632
Submission received: 16 April 2021 / Revised: 4 June 2021 / Accepted: 8 June 2021 / Published: 10 June 2021

Abstract

:
In 2007, ten bioretention cells were constructed in Oklahoma as part of a full-scale technology project to demonstrate stormwater phosphorus reduction. The filter media used was amended with 5%, Class C fly ash by weight to increase phosphorus and heavy metal retention. In 2014, core samples were collected from four of the cells, and three were instrumented for continuous water monitoring for the following year. This paper will review the design, construction, computer modeling of phosphorus retention, and measured phosphorus removal after seven years of operation. Total phosphorus retained in the sampled cells showed reductions in effluent water concentrations of 68 to 75%, while total effluent mass reductions of 51 to 93% were achieved. Total phosphorus accumulation in the cells measured in cores ranged from 0.33 to 0.60 kg/year, which was somewhat greater than the annual calculated effluent reduction of 0.27 to 0.41 kg/year. While good, phosphorus retention was not as high as computer modeling predicted. Other research on the cells, including hydraulics, heavy metal adsorption, and microbial transport, is summarized. Experimental challenges with phosphorus extraction from samples are also discussed. All experience and results suggest that fly ash amendments are an effective option for phosphorus removal in bioretention cells.

1. Introduction

Ten bioretention cells (BRC) were constructed in 2007 as part of a US EPA 319h demonstration project administered by the Oklahoma Conservation Commission with the goal to reduce phosphorus levels in Grand Lake, Oklahoma (OK) [1]. Like other water in the US, Grand Lake experiences algae blooms and other water quality impairments due largely to excess nutrients [2]. Eight of the constructed cells are located in Grove, OK, as shown in Figure 1. Two additional BRC were constructed in Stillwater, OK in 2008. The watershed drainage area ranged from 0.07 to 0.77 ha and is consistent with potential applications in urban areas [3]. Land use on the cells’ drainage varied. Two are residential, six are public or educational, and two are commercial. These cells were unique in three aspects. First, they incorporated Class C fly ash into the filter media to greatly increase phosphorus and metal retention. Secondly, they utilized surface sand plugs to increase infiltration and ensure that the cells would not pool water for more than two days. Finally, and perhaps most importantly, there has been extensive laboratory testing and computer modeling carried out both before and after their construction to better understand the BRC’s long-term removal of multiple pollutants, including phosphorus, heavy metals, and microbes.
In 2014, under US EPA 319h funding administered by the Oklahoma Secretary of the Environment, efforts began to quantify the performance of the BRC and compare it to previous modeling predictions [4]. Core samples were taken from four of the Grove cells, and three of those were instrumented for continuous water quality sampling to quantify effluent phosphorus reductions and phosphorus retention. This paper provides an overview of the fly ash used, cell design, construction, testing, and overall phosphorus retention. Likewise, modeling predictions are contrasted with the cell sampling. Additional research on the cells is also briefly described, and a comparison of core sample phosphorus extraction methods is made.
This is a “big picture” project review intended for both researchers and practitioners. While somewhat unusual, its focus is only to present relevant information on fly ash and the most important project procedures and results for those considering similar work. Some, but not all, of the following has been published previously. In an effort to keep the paper a reasonable length and straightforward for practitioners, many details are omitted. Researchers should consult the citations for relevant related literature, experimental procedures, data and other factors. Marvin et al. [5] have recently presented a systematic literature review of phosphorous sorption in bioretention media. They note that a number of other researchers have more recently conducted laboratory or mesocosm scale tests of phosphorus retention with fly ash amended filter media [6,7,8,9,10,11,12]. Those tests showed that fly ash is generally superior to other additives for phosphorus reduction. However, Marvin et al. [5] also note that only four field scale bioretention studies have been published on any phosphorus adsorbing amendment [13,14,15,16]. Finally they note that the Oklahoma BRC are the only cells to have undergone testing after several years of service.

2. Fly Ash

Fly ash is a combustion product of coal fired power plants that exits with the flue gases and is collected by electrostatic precipitators, gas centrifuges, or filter bags. United States production in 2019 was 29 million tons, of which 18 million tons were used in construction and consumer products [17]. Fly ash should not be confused with other coal combustion products, such as bottom ash, boiler slag, and flue gas desulfurization gypsum.
Since fly ash is a byproduct generated in large quantities, most coal plants give it away to responsible companies to haul it off. Thus, the wholesale price is the cost of transportation, a few dollars a ton. The authors’ experience is that power companies are supportive of an environmental application of their waste, and the fly ash used here was obtained free. In addition, since fly ash is a commonly used concrete additive, many cement suppliers sell it in bulk. The primary safety concern with the use of fly ash is the fine dust particles with large amounts of silica. Workers should use respirators and the public should not be exposed to the dust during storage and handling. Fly ash is not classified as a hazardous waste [18].
Zhang et al. [19,20] carried out pollutant retention screening tests on several materials including fly ash, limestone, and expanded shale. Unweathered Class C fly ash [21] was obtained from the Grand River Energy Center, near Chouteau, OK, and the Sooner Power Plant at Red Rock, OK. The fuel source for both plants was sub-bituminous coal from the Powder River Basin, Wyoming. Table 1 lists properties for the Sooner fly ash. The Grand River fly ash was nearly identical. Chemical composition, was determined by X-ray fluorescence (Activation Laboratories Ltd., Ancaster, ON, Canada). Metal leaching potential was determined by the TCLP Method 1311 [22]. All tested metals concentrations were well below the regulatory level [23].
Fly ash is strongly pozzolanic. When wetted, pure fly ash will form a relatively impermeable mass [24]. Thus, pure fly ash is inappropriate for filter media, and only mixtures of fly ash and clean sand with less than 5% fines were considered. With a final hydraulic conductivity of 9.1 mm/h after 28 days of saturation, construction sand mixed with 5% by weight fly ash possesses an adequate hydraulic conductivity to drain 0.3 m of ponded water within 33 h.
Zhang et al. [20] measured the phosphorus adsorption isotherms of several potential filter media by ASTM D 4646-03 with a 20:1 solution to solid ratio [25]. Five initial concentrations of 1, 3, 10, 20, and 30 mg/L P, as sodium phosphate in 0.01 mol/L KCl solution, were used. The 5% fly ash-construction sand mix adsorbed twice the phosphorus of the next best material, an expanded shale. Construction sand used in the mix originated from a sand and gravel pit excavating the aeolian Dougherty sand that was later processed to contain less than 5% fines (<0.05 mm). Table 2 lists laboratory measurements fits for Langmuir and Freundlich isotherms and field samples described in following sections. The mechanics of the adsorption are ill-defined. However, the mixture released negligible phosphorus in batch desorption tests, indicating that the adsorption is relatively irreversible. After these promising results, it was decided to use fly ash in the demonstration BRC.

3. Bioretention Cells

3.1. Cells Sampled in 2014

Four of the BRC were sampled in 2014 and are listed in Table 3 along with their physical dimensions. Elm Creek Plaza (ECP) is at a busy strip mall and captures runoff from the parking lot, which would otherwise flow directly into a tributary of Grand Lake. The Grand Lake Association (GLA) is a visitor center that has the largest drainage area of the sites chosen for this project. Runoff comes from both a parking lot and the lawn. Grove High School (GHS) readily lent itself to cell construction. There was an existing swale transporting runoff from the staff parking lot to the storm water drain in the city right-of-way. The cell was designed to fit inside the swale. A private residence (SR) is a single family home with a property extending to the lake. The cell intercepts runoff from the property before entering the lake.
USDA soil mapping plot Eldorado silt loam (thermic Typic Paleudolls) at GLA and GHS, Shidler silty clay loam (thermic Lithic Haplustolls) at SR, and Razort silt loam (mesic Mollic Hapludalfs) at ECP [26]. However, all sites exhibited construction related to earth moving. Thus, natural soil profiles were not present.

3.2. Cell Design

In general, the designs for this project are consistent with the BRC guidelines set by Prince George’s County [27], Hunt and White [28], and the LID Center [29]. Figure 2 depicts a typical section of the BRC design for this project. Variations from the references cited include a 1:1 side slope for improved safety during construction, sand plugs for adequate infiltration through the top soil layer, fly ash-sand blend for the filter medium, and cover vegetation suitable for eastern Oklahoma. Construction for this project was formally bid through the Oklahoma Department of Central Services, requiring a complete Plans, Specifications, and Estimate package [1]. The funding agency’s highest priority for these cells was to provide a technology demonstration for the public. Sites were selected that had good public access, cooperative landowners, a variety of land uses, and an attractive installation. As such, cell sizing could not be held constant but had to be adjusted due to site limitations and landowner requests. Recommended cell surface areas for BRC range from 3 to 8% of the runoff area [28]. The four cells ranged from 2.3 to 6.3%. Three of the BRC had a ponding depth of 0.3 m, which produced a pond volume ranging from 7 to 19 mm of watershed first-flush runoff. The ponding depth for ECP was increased by 50% since its watershed is entirely pavement, which corresponded to 20 mm of watershed runoff. Available data on how much of the total contamination are contained in the first-flush is highly variable [30,31]. Nevertheless, monitoring in 2014–2015 presented later shows that the BRC ponds were able to catch almost all the normal rainfall events for the year.
The topsoil thickness in the BRC was 0.3 m, based on the advice of the landscape professional who did the vegetation design for the Grove Cells [32]. The available topsoil in Grove, Oklahoma, is a silty loam, with a hydraulic conductivity much less than 25 mm/h. Even blending equal parts sand and topsoil did not provide adequate hydraulic conductivity to ensure infiltration into the filter layer. Sand plugs that created a channel from the surface to the filter media were utilized. Small areas of clean sand were placed above the surface of the filter media to form the sand plugs. Topsoil was then backfilled around them, rendering forms unnecessary. Figure 3 shows the construction of drains and the filter material, and sand plug placement. Plugs totaling 25% of the cell surface were randomly placed with no two plugs touching. This area of sand plugs was sized to drain the cell within two days, based on one half of sand’s hydraulic conductivity, which would allow for clogging over time. The two-day limit on ponding prevents mosquito reproduction, which is important since several of the cells were at schools and public locations.
Total BRC depth was adjusted due to site topography and the distance the drain had to run to daylight. Average filter depth varied due to the 1:1 side slope and aesthetics adjustments to the cell shape. Filter volumes for the irregular cell shapes were calculated with a 3D model [33]. Average BRC filter depths ranged from 0.9 to 1.2 m, which is consistent with the transport calculations by Zhang et al. [20].
Surface inlet channel dimensions were calculated using a combination of the rational method and Manning’s equation. Since the drainage area size is limited to less than one hectare, the rational method was an acceptable means of calculating the expected surface runoff. A rainfall intensity of 91 mm/h was used that corresponds to a 50-year 1-h storm at Grove [34], which will provide a reasonable safeguard for these non-critical structures on small watersheds. The dimensionless runoff coefficient was assumed to be 0.95 for impermeable surfaces such as pavement and roofs. Manning’s equation was employed to size inlet and outlet channels. Slopes were estimated from the site topography. Depth of flow was set to 0.15 m and channel width calculated. An overflow bypass weir was placed on the back end of the cell. Values for weir height and total head above the weir were held constant at 0.3 m and 0.15 m, respectively, and the length of the weir sized using the broad-crested Weir equation. Chavez et al. [35] present additional design details and calculations.
BRC bottom drains are nominal 0.1 m (4″) diameter slotted corrugated polyethylene pipe with a filter fabric wrap. Within the cells, the drain pipe was spaced on a 1.5 to 2.7 m grid, with one or two pipes exiting the cell. Christianson et al. [36] utilizing agricultural field drainage relationships showed that this number of drain pipes was excessive. Even for the largest BRC, a single pipe down the middle would be adequate. Similarly, it was showed that a 50 mm (2″) pipe would be more than adequate for the flow. However, it was decided to over design this aspect of the BRC to ensure uniform flow through the filter media. Cell bottoms were not sealed since the native soils had low permeablity and no significant leakage or inflow was expected. Stormwater monitoring, as presented later, showed that was not the case for at least two of the Grove BRC.

3.3. Cell Construction

The BRC construction, with one exception, followed normal procedures. The cell volume was excavated with frontend loaders. Evacuated spoil was used in grading berms and for the BRC top soil. Drain pipes were laid on the cell bottom and effluent exit pipes laid in level trenches until they day lighted to the natural grade. Filter material was backfilled, sand plugs and top soil placed, and the cells planted. The one unusual aspect of construction was mixing the fly ash on site. An unnamed alluvial construction sand with less than 5% fines and fly ash were combined in the field by two means shown in Figure 4. The first method involved mixing with heavy equipment before placement in the cell. A load of sand was deposited near the cell site and the required amount of fly ash was then mixed into the sand by repeatedly filling the bucket of a front end loader and pouring it back over the pile containing the sand and fly ash, until an even blend was achieved. The second method used a roto-tiller to mix media in 150 mm lifts inside the cell. While considerable effort was expended on both methods, sampling of the materials during placement showed substantial and nearly equal variation in fly ash content throughout the filter layer [37]. Either method may be used in the future, but other options should be considered.
BRC were planted from a list of native and nonnative species that were wet and dry tolerant, not N fixers, noninvasive, low maintenance, aesthetically pleasing, easily attainable, and replaceable. The list included trees, shrubs, flowering perennials, ornamental grasses, and rock accents. Property owners had the discretion to select final plantings from the approved list, and were responsible for maintenance. Thus, there soon developed a large variation between plant cover on the BRC [38].
To maintain landscaping aesthetics and to help retain moisture for the vegetation, a 50 mm hardwood mulch layer was applied over the topsoil/sand plug layer after planting the vegetation. Hardwood was used to minimize mulch loss due to floating when the cell pond has runoff. No compost was used in the cell due to potential phosphorus leaching. Figure 5 shows two completed BRC, soon after initial planting, while Figure 6 shows two cells in 2014 undergoing hydraulic testing [39].
Figure 7 presents a plot of the BTC construction cost for all 10 cells as a function of the total cell volume. Costs include excavation, drains, filter sand, filter and topsoil placment, site grading, and construction erosion control [37]. Costs do not include design, contract supervision, flyash, or planting. The grove’s cell cost represents full outside contractor construction. The two OSU cells’ cost was lower due to the use of inhouse construction labor. It is expected that large scale projects would be more economical. Chavez et al. [37] reported planting cost ranging from $526 for the smallest cell to $4481 for the largest, GLA.
After construction BRC maintenance was the responsibility of the landowners. GLA, GHS, and SP were well maintained. GLA and GHS sustained a landscaped appearance while SP was cultivated into a butterfly garden. ECP had no maintenance and quickly became overgrown with the planted and volunteer species. Nevertheless, plant cover did not appear to affect the BCE hydraulics.
Initial testing was performed on the BRC to support simulations of the flow hydraulics in the BRC [36,40]. With the simulated hydraulics, Chavez [41] performed stocastic, 3D, finite element modeling of phosphorus transport. A BRC similar in size to ECP was simulated for 144 years, with an influent concentration of 1 mg/L of phosphorus. That relatively high influent concentration was selected as a worst case scenario. Simulations using 75,000 elements examined three ranges of cell variability with 20 random realizations each. Figure 8 shows one simulation concentration distribution at 20 and 144 years, while Figure 9 presents the mean drain effluent concentration with time. Note that the phosphorus water quality standard of 0.035 mg/L would not be exceeded for 34 years. The earliest any simulation exceeded the standard was 17 years. The BRC were predicted to continue to reduce stormwater phosphorus for more than 100 years.

4. 2014 Sampling and Monitoring

Quantification of the BRC performance is not a simple issue. There are three questions to address: (1) how much phosphorus is being retained now, (2) how much phosphorus has been retained over the seven years, and (3) how much more phosphorus will the BRC hold? The first question was addressed by stormwater sampling, the second was quantified by both stormwater sampling and field core sample extraction, while the third was predicted by laboratory adsorption measurements of field core samples.

4.1. Stormwater Flow and Phosphorus Sampling

Stormwater was monitored and sampled at three cells: ECP, GLA, and GHS in 2014 and 2015 [42]. Fourteen-bottle autosamplers with rain gauges collected BRC influent and drain samples. Continuous flow volumes were measured by Palmer–Bowlus flumes [43]. Storm events were defined by the successful collection of water samples using the automated samplers. The flumes triggered their respective samplers programed on an increasing time interval schedule. Rainfall events that did not produce runoff were ignored. Samplers were serviced within one day. Thus, sampled storms were separated by at least 24 h. The collected sampler bottles were proportionally combined based on the flume flow hydrograph to create a flow weighted composite sample. Unfiltered samples were analyzed by the Soil Water and Forage Analytical Laboratory, Oklahoma State University by inductively coupled plasma atomic emission spectroscopy (ICP-AES) for total phosphorus. The difference between the influent and drain sample provides the phosphorus concentration reduction. Mass fluxes for both influent and drains were calculated as a product of the event mean concentration and the total runoff volume measured during a runoff event. Table 4 presents the sampling results. The Wilcoxon Rank-Sum Test, a non-parametric analysis, was used to determine if the BRC treatment made a significant improvement on phosphorus reduction [44].
At the GHS cell, ten sampled storm events occurred from September 2014 to September 2015. The mean storm size was 23 mm, while the median storm size was 19 mm. Event watershed runoff depths ranged from 2 mm to 80 mm. Three BRC overflow events occurred but were not recorded due to autosampler failures. A mean flow volume reduction of 13% was achieved. Total phosphorus effluent concentration was 67% lower than the influent, while T-P drain effluent mass was 84% lower than the influent mass.
At the ECP cell, 20 storm events were sampled storm events occurred from May 2014 to October 2015. Event watershed runoff depth ranged from 6 mm to 97 mm, the mean storm size was 26 mm, and the median storm size was 19 mm. No overflow occurred during the events. Drain flow for three storm events were not recorded due to autosampler failure. Both flow volume and phosphorus concentration were reduced. A flow volume reduction of 73% was achieved. This reduction is attributed to both evapotranspiration in the cell and seepage out of the unlined cell bottom. Total phosphorus effluent concentration was significantly reduced 75% from the influent, while T-P effluent mass was 93% lower than the influent mass.
At the GLA cell, 14 measureable storm events occurred from June 2014 to September 2015. Mean storm size was 37 mm, and the median storm size was 32 mm. Event watershed runoff depths ranged from 11 mm to 92 mm. One overflow event occurred during the event. Flow in the GLA underdrain was higher than the inlet for most of the storms monitored due to groundwater seepage into the cell. The study area has low relief, so the source of the inflow was likely shallow saturated soil water from a large adjacent grassed field, which would have normally seeped into the adjacent lake, but was intercepted by the BRC. Flow from the underdrain was noticed at times with no precipitation or when there was no influent. Overall, underdrain flow volume was 12 times greater than influent during the monitoring period. However, even with this large seepage, T-P drain effluent concentration was 68% lower than the influent. Likewise, T-P drain effluent mass was 51% lower than the inlet influent mass, with a statistical significance of p = 0.09. While a complicated issue for this study, the cell is effectively treating water from other sources before it reaches Grand Lake. Thus, the unlined cell bottom is beneficial to reducing the lake’s phosphorus loading.
Of final note, Table 4 shows that the average drain effluent T-P concentrations ranged from 0.03 to 0.07 mg/L. This compares well with the Scenic Rivers Total Phosphorus Criterion of 0.037 mg/L [45]. However, Chavez’s [41] modeling with a much higher influent phosphorus load predicted that these concentrations could not be reached in the effluent for at least another 10 years.

4.2. Core Sampling Phosphorus Concentrations

Four cells, ECP, GLA, GHS, and SR, were selected for core sampling in 2014 [42]. Six, 38 mm cores were collected along the center line of each cell to a depth of 0.6 m, and it was deemed impractical to core the full depth of the BRC due to potential harm to the cell surface from the increased loads of the truck mounted core tool, and damage to the drains if pierced by the core tool. Visual examination of the cores did not show the sharp layering of construction. The top soil simply faded into the filter media. Apparently, planting, vegetation growth, maintenance, soil settlement, and the high infiltration blended the boundary of the top soil and filter material. Likewise, due to high organic matter, it was not possible to visually identify cores that sampled a surface sand plug versus the top soil. Thus, the 24 cores were simply sectioned into four, 0.15 m lengths. The upper section is assumed to reflect the original top soil, while the lower three are more representative of filter material.
Those 96 samples were subjected to three increasing vigorous extractions: water soluble by 20:1 water extraction (WS-P), weak acid by Mehlich-3 (M3-P), and total by EPA Method 3050 (T-P). M3-P is considered plant available phosphorus and an indicator of environmental availability, while Method 3050 is a total acid digestion which provides total elemental concentrations. Along with those subsamples, 32 samples of the original top soil and filter media collected and stored after the 2007 construction were subjected to identical procedures.
Table 5 and Figure 10, Figure 11 and Figure 12 present the core phosphorus extraction results. Note that the analysis does not differentiate depth within the filter media since the initial samples were hand collected during construction and do not correspond to a specific position within that BRC. Not too surprisingly, a cursory examination of the graph scales shows clear differences between the extraction methods. The WS-P concentrations were less than 2 mg/kg, indicating that most phosphorus was strongly bound and would not be easily leached by rainwater. The M3-P concentrations were an order of magnitude higher at 10 to 40 mg/kg and showed the clearest increase from the near zero initial samples. Total phosphorus had by far the greatest magnitude, ranging from 170 to 440 mg/kg. It also had the smallest relative increase and significant variation between samples. However, it had an average increase of 17 to 83 mg/kg. Taken as a whole, these data show that the fly ash media is strongly retaining a significant mass of phosphorus.
Analysis of variance was performed using a Tukey test for the comparison of phosphorus accumulated over the seven years after construction. Both WS-P and M3-P showed statistically significant increases in phosphorus. Total phosphorus increased on all samples, but at a low significance level (p > 0.05). This is not surprising since the initial samples had high variance relative to the magnitude.
Table 6 presents the phosphorus concentrations by the three extraction methods as a function of depth in the 2014 cores. General transport theory would predict a near exponential reduction in the adsorbed solute concentrations with depth, which was not the case. However, there was a decreasing downward trend. An ANOVA using the General Linear Model [44] shows a linear trend for both WS-P and M3-P. Again, due to the large variance of the T-P samples, statistical significance was low, but the raw data show a trend. The lack of a stronger trend with depth is curious. One-dimensional modeling by Zhang et al. [20] and three-dimensional modeling by Chavez et al. [41] predicted that most of the retained phosphorus would still be near the top of the BRC at the seven-year mark. However, both of those modeling efforts did not consider solute redistribution between infiltration events, evapotranspiration, or vegetation sequestration and removal. It appears that accurate prediction of the transport and the ultimate phosphorus retention within the BRC will require more complex modeling and/or longer term monitoring.

4.3. Annual Phosphorus Retention

The total annual phosphorus reduction was estimated by two methods. First, the measured stormwater concentrations were used for input to the Load Estimator (LOADEST) model [46]. LOADEST is a USGS program for estimating constituent loads in streams and rivers when only partial measurements are available. Given a time series of streamflow, and parameter concentration, LOADEST creates a linear regression model to predict the instantaneous load based on one or more sample input variables including discharge and concentration. LOADEST automatically creates several multiple regression models and selects the best model from those based on the lowest Akaike Information Criteria statistic. Table 7 presents the LOADEST predicted annual T-P reduction of 0.3, 0.27, and 0.41 kg/year at ECP, GHS, and GLA, respectively. These are relatively significant quantities when compared to the cells’ watershed areas.
The second method used the difference between the initial and final core phosphorus concentrations multiplied by the cell filter volume and dry density. Dry densities of 1400 kg/m3 and 1500 kg/m3 were used for the top soil and filter media, respectively. Phosphorus retention rates over the seven years of operation are presented in Table 7. The T-P trapped in the 0.6 m of the media sampled on the four BRCs were 0.4, 0.33, 0.51, and 0.6 kg/year at ECP, GHS, GLA, and SR, respectively. These values compare favorably to the LOADEST calculations. It is the opinion of the authors that, within the context of stormwater, they are the same.

4.4. Remaining Adsorption Potential of 2014 BRC Media

The ability of the fly ash-sand media to continue to retain phosphorus was quantified by measuring adsorption isotherms on samples of the 2014 core filter media and the unweathered initial material. Procedures were identical to those presented in Section 2. Figure 13 displays the actual measurements, while Table 2 presents the fitted Langmuir and Freundlich isotherm parameters. These data reflect adsorption above the phosphorus already retained on the media. Isotherm fits were good but substantially reduced in magnitude from the initial. The differences between the initial and 2014 filter media is a reflection of the phosphorus adsorption and any mineral changes during the seven years of BRC operation. The Langmuir fits produced roughly the same maximum sorption of 400 g/kg. However, the B parameter, which is related to the binding energy, is reduced by two orders of magnitude from roughly 2.5 to 0.03 L/kg. This difference is interpreted to reflect the minerals with the highest binding energy becoming saturated with phosphorus or other ions.
The Freundlich KF parameter provides a convenient comparison of the potential adsorption. KF will equal the mass adsorption potential when the phosphorus solute concentration in the aqueous phase is 1 mg/L. Potential adsorption has been decreased an order of magnitude from roughly 200 to 20 mg/kg, or a 180 mg/kg reduction. This mass is greater than the 17 to 143 mg/kg of T-P extracted from the cores samples listed in Table 5.
There are one clear-cut and two potential causes for the difference between the initial and the 2014 core isotherms. First, obviously, the cores have adsorbed significate phosphorus that can be expected to lower the isotherms proportionally. Next, the initial and construction measurements were carried out quickly with samples that had never been previously wetted, while the cores had been wetted and aged for years. The pozzolanic mineralization with this fly ash is significant for at least 28 days after wetting [47]. Thus, the initial isotherm may not be representative of aged material. Finally, the initial and construction isotherms were carried out with simple KCl solutions. Stormwater contains other ions, which, over the seven years of cell operation, may have taken adsorption sites in addition to those occupied by phosphorus. Nevertheless, the filter media continues to retain significant phosphorus as demonstrated by the water sampling. The confounding differences between the initial and aged samples make prediction of long-term filter phosphorus removal efficiency problematic without another round of sampling.

4.5. Observations on the Extraction Methods

Before ending, it is appropriate to comment on the different soil extraction methods used. It was anticipated and confirmed that WS-P would be low and only a fraction of the phosphorus in the filter. The main purpose of performing WS-P was to show that the phosphorus was strongly bound. Basic mass balance considerations indicate that T-P would be the best measure of phosphorus retention in the BRC. However, the high magnitude and variability of the initial and final sample clouds the estimate of performance. It was hoped that M3-P would provide a clearer picture of adsorption, as it had been used by others [48]. It did show an increase from near zero with time. Unfortunately, the magnitude of the increase was 2 to 18 times lower than the increase of T-P. It appears that, in these cells, most of the phosphorus is retained in forms not released by the M3-P extraction. Therefore, even with its high variability, the T-P extraction is the best for quantifying phosphorus retention.
While greatly dissimilar, the three extractions should not be considered conflicting. The correlation between T-P, M3-P and WS-P was determined using a Pearson correlation coefficient test with a significance level of α = 0.05 chosen for the comparison [49]. Total phosphorus correlated with WS-P (r = −0.01), and M3-P (r = −0.28). Similarly, a correlation between M3-P and WS-P (r = 0.49) was obtained. Thus, the three extractions should be considered complementary and essentially a measure of the spectrum of phosphorus adsorption.
Cursory mineralogy of the phosphorus in three of the initial 2007 samples was determined at the Brookhaven National Synchrotron Light Source II by X-ray absorption near edge structure analysis (XANES) [50] and listed in Table 8. Phosphorus was held as calcium, aluminum, and iron phosphates. The important point to note is the mineral variability between samples, which may explain some of the wide range of measured concentrations of the WS-P and M3-P extractions. These extraction procedures likely have different efficiency on the various minerals.
The ultimate inference from these results is that accurately quantifying phosphorus retention using field filter media samples is challenging. Core sampling aged cells has a great advantage in labor and expense compared to continuous influent and effluent monitoring. However, sample concentration variability weakens most statistical tests. It is recommended that future research use all three methods used here to address variation. However, it may be appropriate to explore the use of an extraction between M3-P and T-P in strength in the hope of reducing sample variability.

5. Other Research on These BRC

Several other aspects of these cells have been examined. Zhang et al. [19] performed batch adsorption of heavy metal retention on fly ash, fly ash and sand mixtures, and other materials. They also carried out bench scale column transport experiments with copper, lead, and zinc. These results were used as inputs to a one-dimensional leaching model. This model predicted that the fly ash-sand mix would retain moderate influent concentrations of those metals for hundreds of years. Chavez et al. [40] performed three-dimensional modeling of a BRC with sand plugs and the measured fly ash sample variability measured during the construction. They found that the spatial variability in the flow and adsorption present in the cells would not seriously impact the cell performance. Christianson et al. [36] performed hydraulic testing on the GLA and GHS cells, using fire hydrant flow to simulate storms. These results were compared to a one-dimensional unsaturated flow model for both wet and dry conditions. After calibration, model flow volumes matched measurements, but the model hydrographs did not display the tailing of the measurements. Chavez et al. [37] provide details of the design, construction, and lessons learned during these projects. After examining one cell failure, they recommend specific guidance for future construction. Coffman et al. [38] quantified plant survival and performance in the Grove BRC, and recommended two plant palettes for future consideration. Kandel [49], in parallel to the phosphorus measurements of Kandel et al. [49], quantified copper, lead, and zinc retention in the BRC. He found that metal concentrations significantly increase in the BRC. No conclusions were made on effluent reductions in copper and lead since most measurements were below detection. However, total zinc mass was shown to be reduced 43 to 88%. McLemore [39] repeated the hydraulic experiments of Christianson et al. [36]. He found that, in the continuous flooding test, both BRC experienced a reduction in the drain flow of about one-half from their earlier rates. However, both cells were still meeting hydraulic design standards and providing pollution retention. Youngblood et al. [51] in concert with Vogel et al. [4] quantified E. coli, enterococci and coliphage effluent removal in three of these BRC. They found that, while removal was highly variable, fly ash amended bioretention in these cells performed 49% better than those with a pure sand filter media layer reported in the literature.
Two of the original fly ash BRC are part of an environmental research and education program in partnership with the Oklahoma State University Botanic Gardens in Stillwater, OK. One cell receives runoff from a short stretch of roadway serving as an entrance into the botanic gardens and a nearby gravel parking lot, respectively, and is experimentally paired with a section of permeable pavement. The second cell is paired with an equivalent cell that does not have the fly ash amendment. In an exploratory study, over six storm events, T-P was low in both influent and drain effulent, probably due to adsorption in the parking lot gravel [16]. However, T-P in the fly ash drain was reduced from the influent by 69%, while the pure sand removed 54%.

6. Summary and Conclusions

This research has demonstrated the potential for using fly ash as an additive in BRC. Average drain effluent T-P concentrations on the three cells sampled ranged from 0.03 to 0.07 mg/L. This compares well with the Scenic Rivers Total Phosphorus Criterion of 0.037 mg/L. Drain effluent T-P mass was reduced 51 to 93% compared to the stormwater influent over the year of monitoring. Performance is made even more noteworthy when it was observed in testing and monitoring that the cell design was able to filter almost all surface runoff. Core sampling on four cells showed T-P retention of 0.33 to 0.60 kg/year over their seven years of operation, which correlates to the 0.27 to 0.41 kg/year estimated from water sampling in the seventh year. These BRC are clearly performing at a good level of phosphorus removal. In addition, other research has shown improved heavy metal and microbe retention. All lab, modeling, and field results justify expanded use of fly ash amendments in stormwater systems.
The design of these cells with the sand plugs to improve infiltration and prevent prolonged surface ponding has proven successful. After seven years of vegetation growth and filtering, the cells still maintain at least one-half of their initial filter rate. Their hydraulic performance should be adequate for the foreseeable future. The only recommendation for construction changes is to use a better method to mix the fly ash and sand.
Quantifying the long-term performance of these BRC will need additional research. Specifically, the measure phosphorus adsorption in the media and concentrations in the effluent was not well predicted by laboratory and computer modeling. Initial isotherm measurements used in both simple 1D modeling and complex 3D modeling predicted that, at seven years, the effluent T-P would be well below the 0.037 mg/L criterion. Likewise, the modeling also predicted that phosphorus would still be concentrated near the cell top. However, coring showed it more uniformly distributed with only a small decreasing trend to the bottom. It is suspected that there are two causes for the inability of the models to better fit the actual cell mechanics. First, phosphorus isotherms of the aged core material indicate loss of adsorption capacity exceeding the mass retained over the seven years of operations. The most likely cause of the decrease is a change in the fly ash-sand filter mineralogy as it aged. However, this decrease could also be the result of adsorption site competition from other stormwater solutes. Second, modeling may require simulation of solute redistribution between infiltration events, evapotranspiration, and/or vegetation sequestration and removal. These added processes could increase modeling complexity by an order of magnitude and necessitate additional field measurement of input parameters.
Core sampling is in the authors’ opinion the gold standard for measuring long-term pollutant retention. It is not easy to question its foundation in simple mass conservation. Likewise, it is much simpler to take a core than to install and maintain multiple stormwater samplers for a year or more. However, when applied to filter materials, the ability to distinguish between initial and adsorbed phosphorus makes the selection of the extraction process problematic. None of the three extractions used here were ideal, but each provided a useful measure. WS-P showed that most of the adsorbed phosphorus was strongly bound—a practical, if not a somewhat obvious, result. M3-P gave the clearest indication of the difference between the initial and seven year samples, but it appeared to measure one-half or less of the phosphorus adsorbed from the stormwater. Average T-P values compared well with the estimate of annual stormwater retention; however, the sample variability greatly weakens the interpretation of the results. A better extraction method is needed for any future core sampling program. Until that occurs, it is recommended that all three of these extractions be carried out.
It has been seven years since the 2014 sampling. A new round of laboratory measurements, modeling, and field monitoring should provide a stronger assessment of the fly-ash BRC long-term performance. If conducted, the procedures used here, incorporating the recommendations for improvement, should be adequate. However, the phosphorus adsorption processes on the weathered fly ash-sand media should be given careful consideration. Extending the fly ash research to other applications is reasonable. In particular, granulized fly ash filters may produce another effective design option.
With 12 years of design, construction, modeling, and testing, the fly ash BRC have one of the best, if not the best, documentation and quantification of any stormwater filter demonstration. However, the nature and limitations of piecing together a long-term research project from a range of agencies prevented better integration. In addition, of course, the researcher’s most common hindsight lament is magnified by the project duration, “If only we knew the results at the start of the project…”.

Author Contributions

Conceptualization, J.R.V. and G.O.B.; methodology, all authors; software, J.R.V., R.A.C. and S.K.; validation, R.A.C. and S.K.; formal analysis, R.A.C. and S.K.; investigation, all authors; resources, J.R.V. and G.O.B.; data curation, R.A.C. and S.K.; writing—original draft preparation, G.O.B., R.A.C. and S.K.; writing—review and editing, all authors; visualization, all authors; supervision, project administration and funding acquisition, J.R.V. and G.O.B. All authors have read and agreed to the published version of the manuscript.

Funding

The design: construction and initial testing of the BRC were funded under the US EPA 319h contract C9-996100-12 administered by the Oklahoma Conservation Commission. The 2014 sampling and monitoring was funded under the US EPA 319h contract C9-00F56701 administered through the Oklahoma Office of the Secretary of the Environment. Additional support was provided by the Oklahoma Agricultural Experiment Station Project OKL02779 and the Oklahoma Cooperative Extension Service. Fly ash was provided by the Grand River Dam Authority and Oklahoma Gas and Electric.

Data Availability Statement

Acknowledgments

The many contributions of colleagues and students over the 12 years of this research are sincerely appreciated. The support of Dan E. Storm, Chad J. Penn, Michael Smolen, Billy J. Barfield, Hailin Zhang, Shanon Phillips, Reid Coffman, Wei Zhang, Reid D. Christianson, Alex J. McLemore, Sheila Youngblood, Grant Graves, and Taber L. Midgley were notable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brown, G.O.; Storm, D.E.; Chavez, R.A. Grand lake, Oklahoma watershed implementation project, task 5.3 bioretention cell design, evaluation and technology demonstration. EPA grant C9-996100-12. In Report to the Oklahoma Conservation Commission, Biosystems & Ag Eng; Oklahoma State University: Stillwater, OK, USA, 2008; p. 64. [Google Scholar]
  2. Grand Lake O’ the Cherokees Watershed Alliance Foundation. Grand Lake Watershed Plan. 2008. Available online: www.ok.gov (accessed on 26 March 2021).
  3. Chavez, R.A.; Brown, G.O.; Storm, D.E. Bioretention cell design and construction specifications. In Proceedings of the ASABE Annual International Meeting, Minneapolis, MN, USA, 17–20 June 2007; p. 072268. [Google Scholar]
  4. Vogel, J.R.; Brown, G.O.; Storm, D.E.; Coffman, R.; Phillips, S. Bioretention cells for mass load reduction of phosphorus and sediment in urban watersheds in Oklahoma, US-EPA grant C9-00F56701. In Report to the Oklahoma Secretary of the Environment, Biosystems & Ag Eng; Oklahoma State University: Stillwater, OK, USA, 2016; p. 98. [Google Scholar]
  5. Marvin, J.T.; Passeport, E.; Drake, J. State-of-the-art review of phosphorus sorption amendments in bioretention media: A systematic literature review. J. Sustain. Water Built Environ. 2020, 6, 03119001. [Google Scholar] [CrossRef]
  6. Wendling, L.A.; Douglas, G.B.; Coleman, S.; Yuan, Z. Nutrient and dissolved organic carbon removal from natural waters using industrial by-products. Sci. Total Environ. 2013, 442, 63–72. [Google Scholar] [CrossRef] [PubMed]
  7. Li, J.; Jiang, C.; Lei, T.; Li, Y. Experimental study and simulation of water quality purification of urban surface runoff using nonvegetated bioswales. Ecol. Eng. 2016, 95, 706–713. [Google Scholar] [CrossRef]
  8. Li, J.; Li, Y.; Zhang, J.; Li, Y. Bioswale column experiments and simulation of pollutant removal from urban road stormwater runoff. Desalin. Water Treat. 2016, 57, 24894–24912. [Google Scholar] [CrossRef]
  9. Houle, J.J.; Ballestero, T.P.; Puls, T.A. The performance analysis of two relatively small capacity urban retrofit stormwater controls. J. Water Manag. Model. 2017, 25, 417. [Google Scholar] [CrossRef]
  10. Li, J.; Li, L.; Dong, W.; Li, H. Purification effects of amended bioretention columns on phosphorus in urban rainfall runoff. Water Sci. Technol. 2018, 78, 1937–1945. [Google Scholar] [CrossRef] [PubMed]
  11. Li, J.; Liang, Z.; Li, Y.; Li, P.; Jiang, C. Experimental study and simulation of phosphorus purification effects of bioretention systems on urban surface runoff. PLoS ONE 2018, 13, e0196339. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, B.; Li, J.; Li, Y.; Li, H. Adsorption characteristics of several bioretention-modified fillers for phosphorus. Water 2018, 10, 831. [Google Scholar] [CrossRef] [Green Version]
  13. Cording, A.; Hurley, S.; Adair, C. Influence of critical bioretention design factors and projected increases in precipitation due to climate change on roadside bioretention performance. J. Environ. Eng. 2018, 144, 04018082. [Google Scholar] [CrossRef]
  14. Liu, J.; Davis, A.P. Phosphorus speciation and treatment using enhanced phosphorus removal bioretention. Environ. Sci. Technol. 2014, 48, 607–614. [Google Scholar] [CrossRef]
  15. Shrestha, P.; Hurley, S.E.; Wemple, B.C. Effects of different soil media, vegetation, and hydrologic treatments on nutrient and sediment removal in roadside bioretention systems. Ecol. Eng. 2018, 112, 116–131. [Google Scholar] [CrossRef]
  16. Brown, G.O.; Chavez, R.A.; Storm, D.E.; Smolen, M.D. Construction and performance of bioretention cells. In Proceedings of the World Environmental and Water Resources Congress 2009, Kansas City, MO, USA, 17–21 May 2009; ASCE: Reston, VA, USA, 2009; pp. 981–990. [Google Scholar]
  17. American Coal Council. Coal-Ash-Production-and-Use Survey Report. 2021. Available online: www.americancoalcouncil.org (accessed on 26 March 2021).
  18. US Environmental Protection Agency. Disposal of coal combustion residuals from electric utilities. Fed. Regist. 2015, 80, 37988. [Google Scholar]
  19. Zhang, W.; Brown, G.O.; Storm, D.E. Enhancement of heavy metals retention in sandy soil by amendment with fly ash. Trans. Am. Soc. Agric. Biol. Eng. 2021, 51, 1247–1254. [Google Scholar] [CrossRef]
  20. Zhang, W.; Brown, G.O.; Storm, D.E.; Zhang, H. Fly ash-amended sand as filter media in bioretention cells to improve phosphorus removal. Water Environ. Res. 2008, 80, 507–516. [Google Scholar] [CrossRef]
  21. ASTM. Standard Specification for Coal Fly Ash and Raw or Calcined Natural Pozzolan for Use in Concrete, C618-19; ASTM International: West Conshohocken, PA, USA, 2019. [Google Scholar]
  22. US Environmental Protection Agency. National Water Quality Inventory: 2000 Report, EPA-841/R-02-001; US Environmental Protection Agency: Washington, DC, USA, 2002.
  23. U.S. Environmental Protection Agency. Toxicity Characteristic. 40 Code of Federal Regulations, Section 261.24; US Environmental Protection Agency: Washington, DC, USA, 2006.
  24. Zhang, W.; Brown, G.O.; Storm, D.E. Enhancement of phosphorus removal in bioretention cells by soil amendment. In Proceedings of the ASABE Annual International Meeting, Portland, OR, USA, 9–12 July 2006; p. 062193. [Google Scholar]
  25. ASTM. Standard Test Method for 24-h Batch Type Experiments of Contaminant Sorption by Soils and Sediments, D4646-03; ASTM International: West Conshohocken, PA, USA, 2004. [Google Scholar]
  26. US Department of Agriculture. Web Soil Survey. Available online: Websoilsurvey.nrcs.usda.gov (accessed on 1 June 2021).
  27. Prince George’s County Department of Environmental Resources. Bioretention Manual. Program & Planning Division; Department of Environmental Resources: Prince George’s County, MD, USA, 2002.
  28. Hunt, W.F.; White, N. Designing Rain Gardens (Bio-Retention Areas). Urban Waterways. AG-588-3; North Carolina Cooperative Extension Service: Raleigh, NC, USA, 2001.
  29. Low Impact Development (LID) Center. Urban Design Tools; Low impact Development Center, Inc.: Beltsville, MD, USA, 1999; Available online: www.lid-stormwater.net (accessed on 5 October 2005).
  30. Lee, J.H.; Bang, K.W.; Ketchum, L.H.; Choe, J.S.; Yu, M.J. First flush analysis of urban storm runoff. Sci. Total Environ. 2002, 293, 163–175. [Google Scholar] [CrossRef]
  31. Hathaway, J.M.; Tucker, R.S.; Spooner, J.M.; Hunt, W.F. A traditional analysis of the first flush effect for nutrients in stormwater runoff from two small urban catchments. Water Air Soil Pollut. 2012, 223, 5903–5915. [Google Scholar] [CrossRef]
  32. Megan, P. (Oklahoma State University, Stillwater, OK, USA). Personal communication, 2006. [Google Scholar]
  33. Autodesk Inc. Autodesk Inventor 10; Autodesk, Inc.: San Rafael, CA, USA, 2006. [Google Scholar]
  34. Hershfield, D.M. Technical Paper No. 40: Rainfall Frequency Atlas of the United States for Durations from 30 Minutes to 24 Hours and Return Periods from 1 to 100 Years; Engineering Division, Soil Conservation Commission Service, USDA: Washington, DC, USA, 1961.
  35. Chavez, R.A.; Brown, G.O.; Storm, D.E. Bioretention cell design for full scale project in Grove, Oklahoma. In Proceedings of the ASABE Annual International Meeting, Portland, OR, USA, 9–12 July 2006; p. 062305. [Google Scholar]
  36. Christianson, R.D.; Brown, G.O.; Chavez, R.A.; Storm, D.E. Modeling field-scale bioretention cells with heterogeneous infiltration media. Trans. Am. Soc. Agric. Biol. Eng. 2012, 55, 1193–1201. [Google Scholar] [CrossRef]
  37. Chavez, R.A.; Brown, G.O.; Coffman, R.R.; Storm, D.E. Design, construction and lessons learned from Oklahoma bioretention cell demonstration project. Appl. Eng. Agric. 2015, 31, 63–71. [Google Scholar]
  38. Coffman, R.R.; Graves, D.; Vogel, J.R.; Brown, G.O. Vegetation in dryland bioretention systems. Landsc. Res. Rec. 2015, 4, 58–71. [Google Scholar]
  39. McLemore, A.J. Evaluation of Established Low Impact Development Techniques: Assessing Aged Bioretention Cells and Clogging Pervious Concrete. Ph.D. Dissertation, Oklahoma State University, Stillwater, OK, USA, 2017. [Google Scholar]
  40. Chavez, R.A.; Brown, G.O.; Storm, D.E. Impact of variable hydraulic conductivity on bioretention cell performance and implications for construction standards. J. Hydraul. Eng. 2013, 139, 707–715. [Google Scholar] [CrossRef]
  41. Chavez, R.A. Oklahoma Bioretention Technology Demonstration and the Effects of Filter Media Heterogeneity on Performance. Ph.D. Dissertation, Oklahoma State University, Stillwater, OK, USA, 2015. [Google Scholar]
  42. Kandel, S.; Vogel, J.R.; Penn, C.J.; Brown, G.O. Phosphorus retention by fly ash amended filter media in aged bioretention cells. Water 2017, 9, 746. [Google Scholar] [CrossRef] [Green Version]
  43. Todeschini, S.; Manenti, S.; Volponi, F.; Ciaponi, C. An analytical methodology for the discharge-stage relation of flexible shape “Palmer-Bowlus” flumes. J. Irrig. Drain. Eng. 2020, 146, 04020017. [Google Scholar] [CrossRef]
  44. SAS Institute Inc. Interface to ADABAS: Reference; SAS Institute Inc.: Cary, NC, USA, 2013. [Google Scholar]
  45. Oklahoma Water Resources Board (OWRB). Oklahoma Scenic Rivers Phosphorus Criteria Review—Majority Report. 2002. Available online: www.owrb.ok.gov/quality/standards/scenicrivers.php (accessed on 26 March 2021).
  46. Runkel, R.L.; Crawford, C.G.; Cohn, T.A. Load Estimator (LOADEST): A FORTRAN Program for Estimating Constituent Loads in Streams and Rivers; US Geological Survey Techniques and Methods Book 4; USGS: Reston, VA, USA, 2004; Chapter A5.
  47. Brown, G.O.; The Application of Fly Ash to Treat Stormwater around Poultry Houses. Oklahoma Water Resources Center. 2017. Available online: Ojs.library.okstate.edu/osu/index.php/OWRC/article/view/7866 (accessed on 3 April 2021).
  48. Komlos, J.; Traver, R.G. Long-term orthophosphate removal in a field-scale stormwater bioinfiltration rain garden. J. Environ. Eng. 2012, 138, 991–998. [Google Scholar] [CrossRef]
  49. Kandel, S. Phosphorus and Metal Sorption in Aged Bioretention Cells with Fly ash Amended Filter Media. Ph.D. Dissertation, Oklahoma State University, Stillwater, OK, USA, 2016. [Google Scholar]
  50. Penn, C.J. (Oklahoma State University, Stillwater, OK, USA). Personal communication, 2014. [Google Scholar]
  51. Youngblood, S.; Vogel, J.R.; Brown, G.O.; Storm, D.E.; McLemore, A.; Kandel, S. Field studies of microbial removal from stormwater by bioretention cells with fly ash amendment. Water 2017, 9, 526. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Locations of Grove BRC.
Figure 1. Locations of Grove BRC.
Water 13 01632 g001
Figure 2. Typical BRC cross section (not to scale).
Figure 2. Typical BRC cross section (not to scale).
Water 13 01632 g002
Figure 3. BRC during construction; (left) Grand Lake Association, (right) Grove High School.
Figure 3. BRC during construction; (left) Grand Lake Association, (right) Grove High School.
Water 13 01632 g003
Figure 4. Mixing fly ash; (left) by front end loaders, (right) by rototiller.
Figure 4. Mixing fly ash; (left) by front end loaders, (right) by rototiller.
Water 13 01632 g004
Figure 5. Completed BRC; (left) GLA, (right) SR.
Figure 5. Completed BRC; (left) GLA, (right) SR.
Water 13 01632 g005
Figure 6. Mature BRC in 2014 ponded during hydraulic testing; (left) GLA, (right) GHS.
Figure 6. Mature BRC in 2014 ponded during hydraulic testing; (left) GLA, (right) GHS.
Water 13 01632 g006
Figure 7. BRC construction cost as a function of volume [37].
Figure 7. BRC construction cost as a function of volume [37].
Water 13 01632 g007
Figure 8. Simulated adsorbed phosphorus concentrations in fly ash-sand BRC; (left) at 20 years, (right) at 144 years [41].
Figure 8. Simulated adsorbed phosphorus concentrations in fly ash-sand BRC; (left) at 20 years, (right) at 144 years [41].
Water 13 01632 g008
Figure 9. Mean effluent concentrations from 20 simulations in 3-D model, for six sand plug cell with influent concentration of 1 mg/L [41].
Figure 9. Mean effluent concentrations from 20 simulations in 3-D model, for six sand plug cell with influent concentration of 1 mg/L [41].
Water 13 01632 g009
Figure 10. WS-P mass concentrations in top soil and filter media in 2007 and 2014.
Figure 10. WS-P mass concentrations in top soil and filter media in 2007 and 2014.
Water 13 01632 g010
Figure 11. M3-P mass concentrations in top soil and filter media in 2007 and 2014.
Figure 11. M3-P mass concentrations in top soil and filter media in 2007 and 2014.
Water 13 01632 g011
Figure 12. T-P mass concentrations in top soil and filter media in 2007 and 2014.
Figure 12. T-P mass concentrations in top soil and filter media in 2007 and 2014.
Water 13 01632 g012
Figure 13. Phosphorus adsorption isotherms in aged BRC media and unweathered fly ash-sand.
Figure 13. Phosphorus adsorption isotherms in aged BRC media and unweathered fly ash-sand.
Water 13 01632 g013
Table 1. Fly ash mineral composition, metal leaching, and bulk properties [20].
Table 1. Fly ash mineral composition, metal leaching, and bulk properties [20].
MineralContent
(%)
MetalConcentration in Leachate (mg/L)
Acetic AcidDeionized WaterRegulatory Level
SiO238.1Arsenic 0.070.025.0
Al2O318.4Cadmium 0.000.001.0
Fe2O35.93Lead 0.000.005.0
MnO0.02Zinc 0.260.01
MgO5.43Chromium 0.330.035.0
CaO22.9Manganese 0.130.00
Na2O1.82Copper 0.020.01
K2O0.56Selenium 0.280.021.0
Ti2O1.39Iron 0.130.01
P2O51.37
BaO0.69Bulk Properties
Cr2O30.01pH11.5
SrO0.30CEC (meg/100 g)78
Loss on ignition0.69Exchangeable Ca (mg/kg)14,300
Total97.6Extractable P (mg/kg)13
Table 2. Langmuir and Freundlich phosphorus isotherm parameters obtained for initial screening, 2007 construction, and 2014 cores from four BRC [20].
Table 2. Langmuir and Freundlich phosphorus isotherm parameters obtained for initial screening, 2007 construction, and 2014 cores from four BRC [20].
MaterialLangmuir Freundlich
qmBr2KFnr2
(mg/kg)(L/kg) (L/kg)
5% Fly Ash-SandLaboratory Screening Sample
3762.460.9492050.2680.994
5% Fly Ash-Sand2007 Construction Sample
3852.890.9982030.2950.981
BRC2014 Core Sample
ECP1570.2250.98831.22.020.945
GLA3830.0310.97119.41.490.986
GHS4310.0210.92117.01.480.976
SR4910.0280.98719.21.350.992
Table 3. BRC locations and physical dimensions for four cells sampled in 2014.
Table 3. BRC locations and physical dimensions for four cells sampled in 2014.
BRCProperty Type-LocationImpervious Cover (%)Drainage Area (ha)Cell Area (m2)Area Ratio (%)Filter Volume (m3)Total Volume (m3)Average Depth (m)Ponding Depth
(m)
Ponded Runoff Depth (mm)
ECPCommercial 36°34′44.8″ N 94°46′02.9″ W1000.251074.397.11281.200.4520
GHSEducational 36°35′55.6″ N 94°44′52.5″ W900.651492.31121611.080.307
GLAPublic 36°37′27.4″ N 94°48′59.3″ W360.773234.23404351.350.3013
SRResidential 36°38′58.2″ N 94°49′34.9″ W130.161016.363.5930.920.3019
Table 4. Stormwater and drain mean T-P concentration, mean mass loading, number of storms sampled (n), and Wilcoxon Rank-Sum Test significance of difference between influent and drain (p).
Table 4. Stormwater and drain mean T-P concentration, mean mass loading, number of storms sampled (n), and Wilcoxon Rank-Sum Test significance of difference between influent and drain (p).
BRCnT-P (mg/L)Mean Mass Loading (g)
InfluentDrain EffluentReduction %pInfluentDrain EffluentReduction %p
ECP200.120.0375<0.053.250.2293<0.05
GHS90.150.0567<0.055.130.8384<0.05
GLA120.230.0768<0.0512.656.2251>0.05
Table 5. Mean and standard deviation of initial and final phosphorus concentrations by extraction method, and Tukey test significance of phosphorus increase in BRCs before and after seven years of operation (p).
Table 5. Mean and standard deviation of initial and final phosphorus concentrations by extraction method, and Tukey test significance of phosphorus increase in BRCs before and after seven years of operation (p).
BRCMaterialExtractionInitial
(mg/kg)
n = 8
Final
(mg/kg)
n = 20
Increase (mg/kg)Increase Significance
p
ECPTop SoilWS-P0.11 ± 0.11.4 ± 0.21.30.001
M3-P1.7 ± 0.127.6 ± 4.625.90.001
T-P225 ± 14308 ± 8783>0.05
Filter MediaWS-P0.1 ± 0.061 ± 0.30.90.001
M3-P3.2 ± 0.27.6 ± 3.74.40.05
T-P361 ± 110440 ± 14079>0.05
GHSTop SoilWS-P0.2 ± 0.011.5 ± 0.31.30.01
M3-P8 ± 0.0734.7 ± 7.726.70.01
T-P265 ± 3.5331 ± 11566>0.05
Filter MediaWS-P0.4 ± 0.050.8 ± 0.30.40.05
M3-P5 ± 1.419 ± 7.5140.01
T-P243 ± 3281 ± 33.738>0.05
GLATop SoilWS-P0.2 ± 0.012.7 ± 1.22.50.05
M3-P10 ± 0.130 ± 5200.01
T-P276 ± 19290 ± 60.814>0.05
Filter MediaWS-P0.3 ± 0.181 ± 0.80.70.05
M3-P13 ± 322.7 ± 8.69.70.05
T-P196 ± 28223 ± 5327>0.05
SRTop SoilWS-P0.1 ± 0.015.7 ± 1.25.60.001
M3-P5 ± 0.0839.8 ± 18.834.80.05
T-P170 ± 17.7312 ± 70142>0.05
Filter MediaWS-P0.2 ± 0.041 ± 0.70.80.05
M3-P3.8 ± 0.312.7 ± 88.90.05
T-P355 ± 129372 ± 3517>0.05
Table 6. Phosphorus concentration and standard deviation by extraction method and significance level of linear reduction with depth after seven years of operation (p).
Table 6. Phosphorus concentration and standard deviation by extraction method and significance level of linear reduction with depth after seven years of operation (p).
BRCExtraction Concentration (mg/kg)Depth (m)Linear Trend
0–0.150.15–0.300.30–0.450.45–0.60Significance
p
WS-P1.4 ± 0.20.97 ± 0.21 ± 0.31.2 ± 0.30.05
ECPM3-P27.5 ± 4.66.5 ± 0.86 ± 0.610 ± 70.001
T-P308 ± 87447 ± 99469 ± 153406 ± 196>0.05
WS-P1.5 ± 0.30.67 ± 0.10.9 ± 0.30.8 ± 0.40.001
GHSM3-P34.6 ± 7.615.4 ± 2.816.6 ± 4.631.2 ± 6.30.001
T-P331 ± 114279 ± 50282 ± 23285 ± 9>0.05
WS-P2.6 ± 1.20.7 ± 0.11.4 ± 11.4 ± 0.60.01
GLAM3-P29 ± 515 ± 1.826 ± 825 ± 90.01
T-P290 ± 60207 ± 33255 ± 75205 ± 300.05
WS-P5 ± 1.21 ± 0.41 ± 0.61.2 ± 0.90.001
SRM3-P39.7 ± 189.6 ± 3.57.5 ± 1.720.6 ± 100.001
T-P312 ± 170414 ± 267387 ± 238315 ± 230>0.05
Table 7. Phosphorus mass retention rate within the BRC by core sampling and LOADEST modeling with water sampling results.
Table 7. Phosphorus mass retention rate within the BRC by core sampling and LOADEST modeling with water sampling results.
BRCMediaDepth
(m)
Phosphorus Retained by Extraction Method (kg/year)Phosphorus Retained by
LOADEST
T-P
(kg/year)
WS-PM3-PT-P
ECPTop Soil0.150.0020.050.160.30
Filter media0.450.0060.030.24
Total0.600.0080.080.4
GHSTop Soil0.150.0040.080.290.27
Filter media0.450.0040.130.04
Total0.600.0080.210.33
GLATop Soil0.150.0120.090.070.41
Filter media0.450.010.140.44
Total0.60.0220.230.51
SRTop Soil0.150.0290.180.43Not Monitored
Filter media0.450.0130.150.17
Total0.600.0420.330.60
Table 8. XANES analysis of three initial samples taken at construction in SR cells; (a) 0 to 150 mm, (b,c) 300 to 450 mm separated by 2 m [50].
Table 8. XANES analysis of three initial samples taken at construction in SR cells; (a) 0 to 150 mm, (b,c) 300 to 450 mm separated by 2 m [50].
GroupMineral (%)Sample
abc
Calcium PhosphatesBrushite1320
Monetite152112
Hydroxyapatite192128
b-tricalcium phosphate7106
Octacalcium phosphate81113
Aluminum PhosphatesVariscite99
Al Oxide Sorbed PO49 12
Iron PhosphatesStrengite 11
Amorphous Fe Phospate88
Fe Oxide Sorbed P0411 13
Phytic Acid 5
R-Factor 0.02040.06660.0207
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vogel, J.R.; Chavez, R.A.; Kandel, S.; Brown, G.O. Field Scale Demonstration of Fly Ash Amended Bioretention Cells for Stormwater Phosphorus Removal: A Review of 12 Years of Work. Water 2021, 13, 1632. https://doi.org/10.3390/w13121632

AMA Style

Vogel JR, Chavez RA, Kandel S, Brown GO. Field Scale Demonstration of Fly Ash Amended Bioretention Cells for Stormwater Phosphorus Removal: A Review of 12 Years of Work. Water. 2021; 13(12):1632. https://doi.org/10.3390/w13121632

Chicago/Turabian Style

Vogel, Jason R., Rebecca A. Chavez, Saroj Kandel, and Glenn O. Brown. 2021. "Field Scale Demonstration of Fly Ash Amended Bioretention Cells for Stormwater Phosphorus Removal: A Review of 12 Years of Work" Water 13, no. 12: 1632. https://doi.org/10.3390/w13121632

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop