Journal of Environmental Protection
Vol.08 No.09(2017), Article ID:78818,11 pages
10.4236/jep.2017.89065

Locally Sourced Iron and Aluminum Byproducts Decrease Phosphorus Leached from Broiler House Dust Deposited near Ventilation Fans

Sheri L. Herron1, Andrew N. Sharpley2*, Kristofor R. Brye2, David M. Miller2

1Ecosustain, LLC, Dripping Springs, TX, USA

2Department of Crop, Soil, and Environmental Sciences, University of Arkansas, Fayetteville, AR, USA

Copyright © 2017 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: July 12, 2017; Accepted: August 27, 2017; Published: August 30, 2017

ABSTRACT

Freshwater impairment by eutrophication, as a result of excessive phosphorus (P) inputs from runoff in particular, remains a ubiquitous environmental concern. A common issue with systems designed to remove P and nitrogen (N) from runoff is their reduced effectiveness under high-flow conditions. To overcome this, P removal from broiler-house fan dust would be more effective if removal occurred at the nutrient source, where the water volume is limited to direct rainfall. The P removal efficiencies of different thicknesses of locally sourced, iron-rich red mud (RM) generated during the manufacture of steel belts for tires and alum-based drinking water treatment residual (WT) byproducts were investigated. Byproduct thicknesses of 4, 8, and 12 cm were tested using 57-L leaching columns. The columns were filled with the specified byproduct thickness and a 3-mm thickness of poultry house dust was surface applied prior to receiving six, 30-min simulated rainfalls (at 7 cm∙hr−1) at 1-day intervals. The 8-cm thickness of both RM and WT outperformed the other thicknesses in terms of sorbing P released from the added broiler house dust, removing 99 and 96% of the added P, respectively, over the six simulated rainfall events. The 12-cm thickness of both RM and WT showed no additional benefit for P removal over the 8-cm thickness. As the 4-cm-thick WT treatment was less effective (89% of added P removed), the 8-cm thickness was the optimal thickness for field testing. Locally sourced materials with large P-sorbing capacities can offer a convenient, relatively inexpensive alternative for P removal from areas around poultry houses impacted by P-containing, exhausted broiler house dust.

Keywords:

Phosphorus Runoff, Poultry Production, Red Muds, Water Quality, Water Treatment Residual

1. Introduction

Freshwater-use impairment by accelerated eutrophication from increased nutrient inputs, particularly phosphorus (P), remains despite widespread conservation efforts [1] [2] [3] . For instance, harmful algal blooms have been linked to excess P in western Lake Erie [4] [5] and Florida [6] [7] , as well as to hypoxia in the Northern Gulf of Mexico [8] [9] [10] .

The potential for P and nitrogen (N) runoff from land around poultry production facilities has been noted with two major sources being dust exhausted from the production houses by fans [11] and poultry litter material tracked from houses during bird removal and clean-out [7] . While dust and broiler litter have similar total P concentrations, the impact of P release from ventilation fan dust is exacerbated due to the water solubility of dust P being, on average, three times greater than that in the litter [12] . Removing nutrients from runoff water prior to water leaving the poultry production area is likely the most cost-effective on-farm mitigation practice. Similarly, in an effort to most cost-effectively decrease impairment, strategies to reduce nutrient runoff at a watershed scale are more successful when targeting nutrient sources rather than treating receiving waters [13] [14] [15] .

Substantial research has been conducted evaluating material that can sorb large amounts of P and structures to contain such materials [16] [17] . To be cost-effective, the P-retention material needs to be locally sourced, readily available, and inexpensive. Two such byproduct materials are a group referred to as red muds (RM) or iron filter cakes, which consist primarily of iron oxides and oxyhydroxides, and an alum-based water treatment (WT) residual.

The byproducts used in this study were an iron filter cake RM from a company in northwest Arkansas, generated during the manufacture of steel belts for tires, and an alum-based WT generated by the Beaver Water District in northwest Arkansas during flocculation of suspended solids from raw municipal water supplies. Currently, a major portion of these materials is land filled. The P-sorption characteristics and hydraulic properties of these materials have been shown to provide potential as on-site sinks for P [11] . For instance, Herron et al. [11] reported RM and WT could sorb 25 and 10 g P kg−1 of byproduct, respectively, with hydraulic conductivities of 8.0 and 15.4 cm∙min−1 [18] , prior to being used to sorb P from broiler house ventilation dust. However, the capacity of these materials to sorb and retain P released from ventilation fan dust must be assessed under conditions approximating those encountered in the field under environmental exposure.

The objective of this study was to evaluate the capacity of different thicknesses (i.e., 4, 8, and 12 cm) of RM and WT byproducts to sorb and retain P and N released from broiler house fan dust as a result of sequential simulated rainfall events (7 cm∙hr−1) applied to leachate columns. It was anticipated that the two byproducts would behave similarly, but that an optimum byproduct thickness would be identified.

2. Materials and Methods

2.1. Experimental Design

It was specifically hypothesized that an on-farm nutrient removal system could be designed to reduce nutrients in stormwater runoff to within acceptable limits by trapping the nutrients at the source. To test this hypothesis, a rainfall simulation and a small-scale containment structure were used to test various byproduct thicknesses for nutrient removal from broiler house dust (BHD). The containment structure (i.e., leachate column) was designed by adapting 56.8-L buckets to contain the RM and WT byproducts. Dust was applied to the surface of the byproducts in each leachate column followed by simulated rainfall application. The leachate produced was then collected for chemical analysis.

2.2. Column Configuration

Prior testing of the RM and WT byproducts [11] showed the hydraulic conductivity of 6-mm-diameter particles of RM was 123 cm 30 min−1, while that of the 6-mm-diameter particles of WT was 138 cm 30 min−1. These hydraulic conductivities were determined to be more than adequate to transmit the calculated volume of rainfall that would fall to the ground adjacent to broiler house fan outlets during a 25-yr, 30 min storm event (i.e., 33 cm per 30 min).

Three byproduct thicknesses (i.e., 4, 8, and 12 cm) were tested using the 56.8-L leachate columns. The columns were constructed by removing the bottoms of 27, 56.8-L buckets, leaving the remainder of the bucket for the specified byproduct thickness plus 7.6 cm of open space at the top. Welded wire mesh, with 1.27 cm openings, was attached to the bottom of the buckets and landscape fabric was placed on top of the wire mesh. The byproduct was placed on the landscape fabric. The cut buckets were placed inside uncut 56.8-L buckets, perforated on the sides to allow airflow, to capture leachate during the experiment (Figure 1). Leaching experiments were conducted in triplicate for each byproduct-thickness treatment combination, including controls for each byproduct that received no dust and dust only treatments (i.e., no byproduct).

Dust collected off fan shutters from four broiler facilities in northwest Arkansas was combined and mixed thoroughly as described by Herron et al. [12] . The dust application rate for the experiment was calculated based on direct measurements of dust thicknesses that had settled on the ground adjacent to fan outlets during one grow-out cycle (flock) of birds (approximately 6 weeks), which averaged 3.2 mm per flock. Based on the surface area of the leachate

Figure 1. Schematic of stacked columns used in rainfall simulation experiment.

column (638 cm2) and the estimated bulk density of the dust (0.23 g∙cm−3), 45.43 g of dust were used for each byproduct-thickness treatment combination. The dust was analyzed by the University of Arkansas Agricultural Diagnostic Laboratory and contained 4312 mg water extractable P kg−1, 10,400 mg total P kg−1, 3,785 mg NH4-N kg−1, 277 mg NO3-N kg−1, and 97,600 mg total N kg−1. The dust was applied by hand uniformly to the surface of byproduct in each column. For the no-byproduct treatment, dust was applied by hand uniformly directly to the landscape fabric.

2.3. Rainfall Simulation

The rainfall simulator used for this experiment was based on the design of Miller (1987). The dimensions for the area under the simulator were 2.0 meters wide by 1.5 meters long. The frame was built of hollow metal tubing, on which a 25-mm-diameter, polyvinyl chloride (PVC), water-supply pipe was mounted. Polyvinyl chloride windscreens were attached to four sides to limit wind effects on rainfall distribution. A water hose supplying water from the city water supply was attached to the water-supply pipe. A single, fixed nozzle (TeeJetTM 1/2HH-SS-5OWSQ, TeeJet, Glendale Heights, IL), designed for a flow rate of 210 mL∙sec−1, was installed on the water-supply pipe and centered at the top of the frame, 3 m above the surface of the leachate columns [19] . A low-pressure regulator and an oil-filled pressure gauge were used to maintain the water pressure at 37.9 kPa (5.5 PSI). An in-line filter was placed in the supply pipe to prevent particles clogging the line and nozzle during water delivery.

Leachate columns were randomly placed under the rainfall simulator in three zones, one zone for each replicate set of treatment combinations (Figure 2, Figure 3). Rainfall intensity was verified and distribution uniformity was determined by

Figure 2. Assignment of columns into three groups (white, gray, and black) based on randomization plan.

Figure 3. Actual layout of columns under rainfall simulator based on rainfall uniformity test.

simulating rainfall three times onto the area where the leachate columns were to be placed and measuring the resulting rainfall volume. The exact locations of the three zones were then selected based on the results of the rainfall uniformity tests. The coefficient of uniformity was determined using the method described by Christiansen [20] , which expresses uniformity as a percentage according to the formula coefficient of uniformity = 100 * (1.0 ? Σd/m * n), where d is the standard deviation of individual observations from the mean value m, and n is the number of observations [21] .

Rainfall simulations were conducted at an intensity of 7 cm∙hr−1 based on a 2-year, 24-hour storm rate and leachate was collected for a total of 30 minutes according to the National Phosphorus Research Project [19] . Leachate columns received 30 minutes of simulated rainfall prior to the first treatment with dust to equalize their water contents. After the first rainfall simulation, dust was surface applied to the byproducts in each column, which then received simulated rainfall. Six consecutive dust applications were followed by 30 minutes of simulated rainfall at a rate of 7 cm∙hr−1, the combination of which occurred at one-day intervals. Dust applications approximated the amount of dust deposited on the soil outside a broiler house fan during one year of typical operation (i.e., the production of six flocks of birds; [12] ).

2.4. Chemical Analyses

Following each rainfall simulation, leachate was collected from the bottom of the columns. Leachate volume was measured and a 1-L sub-sample was stored at 4˚C and analyzed within 24 hours of collection. A sample of source water used for the simulated rainfall was collected and similarly retained for analysis. Forty milliliters of leachate were filtered (0.45-µm membrane) immediately after collection and stored at 4˚C until analyzed for dissolved-reactive P (DRP) by the colorimetric molybdenum-blue method [22] . Ammonium-N (NH4-N) and nitrate-N (NO3-N) were analyzed colorimetrically by flow-injection analyses (Lachat Instruments QuikChem 8500, Loveland, CO). One hundred twenty-five milliliters of leachate were acidified to a pH ≤ 2 with 12 drops of concentrated sulfuric acid for sample preservation and analyzed for total P and N following persulfate-autoclave digestion [23] . Total P (TP) was determined by spectrophotometry (Beckman Coulter, Pasadena, CA) and total N (TN) was determined by flow-injection analysis (Lachat Instruments QuickChem 8500, Loveland, CO). One hundred twenty-five milliliters were analyzed gravimetrically for total solids after oven drying at 105˚C for 12 h.

2.5. Statistical Analyses

The TP results were tested for and had homogeneous variances. Therefore, a two-factor analysis of variance (ANOVA) was conducted on the concentrations of TP to evaluate the effects of byproduct, byproduct thickness, and their interaction. Due to non-homogeneous variances, separate one-factor ANOVAs were conducted on the leachate concentrations of DRP, NH4-N, NO3-N, and TN to separately evaluate the effect of byproduct for each thickness and thickness for each byproduct. All relationships were reported as statistically significant at the p < 0.05 level unless noted otherwise. Wolfram Mathematica (online beta version 2014, Wolfram, Champaign, IL) was used to conduct all statistical analyses.

3. Results and Discussion

3.1. Simulated Rainfall Distribution

The average rainfall uniformity was 95%, 94%, and 87% for zones 1, 2, and 3, respectively, as determined by the coefficient of uniformity (Table 1). Additionally, the coefficient of uniformity among the three consecutive 30 min, 7 cm∙hr−1

Table 1. Coefficient of uniformity of rainfall on the columns for three rainfall events.

rainfalls was 93%. Thus, the configuration of leachate columns in the rainfall simulator produced uniform and repeatable rainfall both within the simulator footprint and among individual simulated rainfall events. For instance, coefficients of uniformity of 90% to 93% were considered excellent by Humphry et al. [21] and Miller [24] , and the National Phosphorus Research Project protocol [19] recommends a coefficient of uniformity of 85% or greater for standardized rainfall simulations using the same rainfall simulation equipment as used in this study.

3.2. Rainfall Simulation on Byproducts Treated with Broiler House Dust

Leachate concentrations for all thicknesses of RM ranged from 0.001 to 1.8 mg DRP L−1, 0.28 to 16.7 mg TP L−1, 0.6 to 442.5 mg NH4-N L−1, 0.0 to 4.8 mg NO3-N L−1, and 5.5 to 450.2 mg TN L−1 (Table 2). Leachate concentrations for all thicknesses of WT ranged from 0.03 to 9.8 mg DRP L−1, 0.43 to 23.3 mg TP L−1, 51.4 to 307.3 mg NH4-N L−1, 0 to 114.0 mg NO3-N L−1, and 0.4 to 346.6 mg TN L−1 (Table 2). Neither thickness, byproduct type, nor their interaction affected (p > 0.05) the concentrations of DRP, TP, NH4-N, NO3-N, or TN (Table 2), when the thicknesses and byproducts were compared together. However, when the thicknesses and type of byproduct were compared separately, the DRP concentration was lower (p < 0.05) for RM than for WT for the 4- and 12-cm byproduct thicknesses (Table 2). Similarly, the NH4-N concentration was also lower (p = 0.02) for RM than for WT for the 8-cm byproduct thickness (Table 2). In addition, the NO3-N concentration was lower (p < 0.002) for RM than for

Table 2. Byproduct column study―average dissolved P, total P, ammonium-N (NH4-N), nitrate-N (NO3-N), and total N concentrations in leachate..

#Dissolved reactive phosphorus (DRP), total phosphorus (TP), ammonium nitrogen (NH4-N), nitrate nitrogen (NO3-N), and total nitrogen (TN). *Two-way ANOVA conducted on data with homogeneous variances. One-way ANOVA conducted on data without homogeneous variances.

WT for the 8- and 12-cm byproduct thicknesses. Comparing the leachate analyses of dust without byproduct to the dust with byproduct treatments, the 8-cm-thick treatments of both byproducts were effective at reducing leachate TP, with RM sorbing up to 94% and WT up to 87% of P leached from dust (Table 3). In the case of N, the RM was more effective than WT in reducing leachate NH4-N concentrations. The 8 cm-depth treatment of RM reduced leachate NH4-N concentrations by 75%.

4. Summary and Conclusions

Based on the observations of the physical characteristics of the byproducts and the average leachate concentrations for DRP and TP, the 8-cm thickness appears to be the desired thickness for subsequent field-testing. Additionally, the intent of this project was to design a cost-effective treatment system. While there was no cost to obtain the byproducts, there was significant cost to pretreat byproducts,

Table 3. Mean dissolved P concentrations of column leachate and P removal by byproducts treated with broiler house dust and subject to six simulated rainfalls.

in terms of drying and creating an appropriate particle-size range. Therefore, it was determined that using a thinner layer (i.e., 8- rather than 12-cm thickness) of byproduct in containment structures would require less byproduct and result in beneficial cost savings to the end user, especially as little difference in P removal was observed between the 8- and 12-cm thicknesses.

Field-testing of RM and WT byproducts for their effectiveness in sorbing nutrients from broiler house dust is needed for future design criteria for on-farm P reduction practices. Preliminary findings from the analyses of the leachate from rainfall simulations indicate both RM and WT byproducts are effective at removing P from broiler house dust prior to the P entering stormwater runoff. Additional opportunities for use of these byproducts in agricultural applications, such as containment in ponds, containment to intercept runoff in drainage swales, and land application at nutrient sources should be explored.

Acknowledgements

This research was supported by a grant from the USGS 104B program funding.

Cite this paper

Herron, S.L., Sharpley, A.N., Brye, K.R. and Miller, D.M. (2017) Locally Sourced Iron and Aluminum Byproducts Decrease Phosphorus Leached from Broiler House Dust Deposited near Ventilation Fans. Journal of Environmental Protection, 8, 1026-1036. https://doi.org/10.4236/jep.2017.89065

References

  1. 1. Dubrovsky, N.M., Burow, K.R., Clark, G.M., Gronberg, J.M., Hamilton, P.A., Hitt, K.J., Mueller, D.K. and Munn, M.D. (2010) The Quality of our Nation’s Waters— Nutrients in the Nation’s Streams and Groundwater, 1992-2004: U.S. Geological Survey Circular 1350. 174 p. http://water.usgs.gov/nawqa/nutrients/pubs/circ1350

  2. 2. Sharpley, A.N., Jarvie, H.P., Buda, A., May, L., Spears, B. and Kleinman, P. (2013) Phosphorus Legacy: Overcoming the Effects of Past Management Practices to Mitigate Future Water Quality Impairment. Journal of Environmental Quality, 42, 1308-1326. https://doi.org/10.2134/jeq2013.03.0098

  3. 3. U.S. Environmental Protection Agency (2015) Watershed Assessment, Tracking & Environmental Results. National Probable Sources Contribution to Impairments. http://iaspub.epa.gov/waters10/attains_nation_cy.control#prob_source

  4. 4. Michalak, A.M., Anderson, E. J., Beletsky, D., Boland, S., Bosch, N.S., Bridgeman, T.B., Chaffin, J.D., Cho, K., Confesor, R., Daloglu, I., DePinto, J.V., Evans, M.A., Fahnenstiel, G.L., He, L., Ho, J.C., Jenkins, L., Johengen, T.H., Kuo, K.C., Laporte, E., Liu, X., McWilliams, M.R., Moore, M.R., Posselt, D.J., Richards, R.P., Scavia, D., Steiner, A.L., Verhamme, E., Wright, D.M. and Zagorski, M.A. (2013) Record-Setting Algal Bloom in Lake Erie Caused by Agricultural and Meteorological Trends Consistent with Expected Future Conditions. Proceedings of the National Academy of Sciences of the United States of America, 110, 6448-6452. https://doi.org/10.1073/pnas.1216006110

  5. 5. Scavia, D., Allan, J.D., Arend, K.K., Bartell, S., Beletsky, D., Bosch, N.S., Brandt, S.B., Briland, R.D., Daloglu, I., DePinto, J.V., Dolan, D.M., Evans, M.A., Farmer, T.M, Goto, D., Han, H., Hoeoek, T.O., Knight, R., Ludsin, S.A., Mason, D., Michalak, A.M., Richards, R.P., Roberts, J.J., Rucinski, D.K., Rutherford, E., Schwab, D.J., Sesterhenn, T.M., Zhang, H. and Zhou, Y. (2014) Assessing and Addressing the Re-Eutrophication of Lake Erie: Central Basin Hypoxia. Journal of Great Lakes Research, 40, 226-246. https://doi.org/10.1016/j.jglr.2014.02.004

  6. 6. Reddy, K.R., Newman, S., Osborne, T.Z., White, J.R. and Fitz, H.C. (2011) Phosphorus Cycling in the Greater Everglades Ecosystem: Legacy Phosphorus Implications for Management and Restoration. Critical Reviews in Environmental Science and Technology, 41, 149-186. https://doi.org/10.1080/10643389.2010.530932

  7. 7. U.S. Environmental Protection Agency—Scientific Advisory Board (2011) Review of EPA’s Draft Approaches for deriving numeric nutrient criteria for Florida’s Estuaries, Coastal Waters, and Southern Inland Flowing Waters. EPA Scientific Advisory Board, Washington DC, 67 p. http://yosemite.epa.gov/sab/sabproduct.nsf/DCC3488B67473BDA852578D20058F3C9/$File/EPA-SAB-11-010-unsigned.pdf

  8. 8. Alexander, R.B., Smith, R.A., Schwarz, G.E., Boyer, E.W., Nolan, J.V. and Brakebill, J.W. (2008) Differences in Phosphorus and Nitrogen Delivery to the Gulf of Mexico from the Mississippi River. Environmental Science & Technology, 42, 822-830. https://doi.org/10.1021/es0716103

  9. 9. Dale, V.H., Kling, C.L., Meyer, J.L., Sanders, J., Stallworth, H., Armitage, T., Wangsness, D., Bianchi, T.S., Blumberg, A., Boynton, W., Conley, D.J., Crumpton, W., David, M.B., Gilbert, D., Howarth, R.W., Lowrance, R., Mankin, K.R., Opaluch, J., Paerl, H.W., Reckhow, K., Sharpley, A.N., Simpson, T.W., Snyder, C. and Wright, D. (2010) Hypoxia in the Northern Gulf of Mexico. Springer Series on Environmental Management. Springer Science, New York. https://doi.org/10.1007/978-0-387-89686-1

  10. 10. Rebich, R.A., Houston, N.A., Mize, S.V., Pearson, D.K., Ging, P.B. and Hornig, C.E. (2011) Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico from Streams in the South-Central United States. JAWRA Journal of the American Water Resources Association, 47, 1062-1086. https://doi.org/10.1111/j.1752-1688.2011.00583.x

  11. 11. Herron, S.L., Sharpley, A.N., Brye, K.R., Miller, D.M., Watkins, S., McCreery, D. and Daniels, M.B. (2016) Determination of Nutrient Concentrations in Simulated Rainfall-Runoff from Poultry House Dust Deposited Adjacent to Exhaust Fans. Journal of Environmental Protection, 7, 27-40. https://doi.org/10.4236/jep.2016.71003

  12. 12. Herron, S.L., Brye, K.R., Sharpley, A.N., Miller, D.M. and Daniels, M.B. (2015) Nutrient Composition of Dust Emitted from Poultry Broiler Houses in Northwest Arkansas. Journal of Environmental Protection, 6, 1257-1267. https://doi.org/10.4236/jep.2015.611110

  13. 13. Herron, S.L., Sharpley, A.N., Brye, K.R. and Miller, D.M. (2016) Optimizing Hydraulic and Chemical Properties of Iron and Aluminum Byproducts for Use in On-Farm Containment Structures for Phosphorus Removal. Journal of Environmental Protection, 7, 1835-1849. https://doi.org/10.4236/jep.2016.712146

  14. 14. Stout, W.J., Sharpley, A.N., Gburek, W.J. and Pionke, H.B. (1999) Reducing Phosphorus Export from Croplands with FBC Flyash and FGD Gypsum. Fuel, 78, 175-178. https://doi.org/10.1016/S0016-2361(98)00141-0

  15. 15. US Department of Agriculture (USDA) and US Environmental Protection Agency (USEPA) (1999) Unified National Strategy for Animal Feeding Operations. US Government Printing Office, Washington DC. http://www.epa.gov/npdes/pubs/finafost.pdf

  16. 16. Oliver, I.W., Grant, C.D. and Murray, R.S. (2011) Assessing Effects of Aerobic and Anaerobic Conditions on Phosphorus Sorption and Retention Capacity of Water Treatment Residuals. Journal of Environmental Economics and Management, 92, 960-966. https://doi.org/10.1016/j.jenvman.2010.11.016

  17. 17. Vohla, C., Koiv, M., Bavor, H.J., Chazarenc, F. and Mander, U. (2011) Filter Materials for Phosphorus Removal from Wastewater in Treatment Wetlands—A Review. Journal of Ecological Engineering, 37, 70-89. https://doi.org/10.1016/j.ecoleng.2009.08.003

  18. 18. Sharpley, A.N., Kleinman, P.J.A., Flaten, D.N. and Buda, A.R. (2011) Critical Source Area Management of Agricultural Phosphorus: Experiences, Challenges and Opportunities. Water Science and Technology, 64, 945-952. https://doi.org/10.2166/wst.2011.712

  19. 19. Southern Extension Research Activities-17 (2017) National Phosphorus Research Project for Simulated Rainfall-Surface Runoff Southern Extension-Research Activity-17. National Phosphorus Research Project for Simulated Rainfall-Surface Runoff Studies. Southern Extension-Research Activity-Information Exchange Group 17, Virginia Tech. University, Blacksburg, VA. https://sera17dotorg.files.wordpress.com/2015/02/national-p-protocol.pdf

  20. 20. Christiansen, J.E. (1942) Irrigation by sprinkling. University of California Agriculture Experiment Station Bulletin 670, Berkeley, CA.

  21. 21. Humphry, B.J., Daniel, T.C., Edwards, D.R. and Sharpley, A.N. (2002) A Portable Rainfall Simulator for Plot-Scale Runoff Studies. Transactions of Applied Engineering in Agriculture, 18, 199-204.

  22. 22. Murphy, J. and Riley, J.P. (1962) A Modified Single Solution Method for Determination of Phosphate in Natural Waters. Analytica Chimica Acta, 27, 31-36. https://doi.org/10.1016/S0003-2670(00)88444-5

  23. 23. Patton, C.J. and Kryskalla, J.R. (2003) Methods of Analysis by the U.S. Geological Survey National Water Quality Laboratory: Evaluation of Alkaline Persulfate Digestion as an Alternative to Kjeldahl Digestion for the Determination of Total and Dissolved Nitrogen and Phosphorus in Water. U.S. Geological Survey, Water Resources Investigations Report 03-4174. U.S. Geological Survey, Branch of Information Services, Federal Center, Denver, CO, 33 p. http://nwql.usgs.gov/pubs/WRIR/WRIR-03-4174.pdf

  24. 24. Miller, W.P. (1987) A Solenoid-Operated, Variable Intensity Rainfall Simulator. Soil Science Society of America Journal, 51, 832-834. https://doi.org/10.2136/sssaj1987.03615995005100030048x