DETECTING CHANGE IN WATER QUALITY FROM IMPLEMENTATION OF LIMESTONE TREATMENT SYSTEMS IN A COAL-MINED WATERSHED, PENNSYLVANIA

Charles A. Cravotta III and Jeffrey B. Weitzel U.S. Geological Survey, 215 Limekiln Road, New Cumberland, Pennsylvania 17070

During 1996-97, a variety of limestone-based treatment systems were implemented to neutralize acidic mine drainage and reduce the transport of dissolved metals in the northern part of the Swatara Creek watershed, which drains a 43-mi2 (112-km2) area in the Southern Anthracite Field upstream from Ravine, Pa. Since 1996, the current project has monitored water quality upstream and downstream of each treatment and at integrator sites on lower reaches of Swatara Creek. Continuous measurements of pH and specific conductance and periodic sampling for , acidity, sulfate, and metals upstream and downstream of each treatment system show that (1) open limestone channels and limestone-sand dosing generally had negligible effects on water quality and (2) limestone diversion wells and limestone drains generally were effective at producing near-neutral pH and attenuating dissolved metals during baseflow but were less effective during stormflow conditions. Storm runoff in this area commonly is acidic, and, as streamflow volume increases during stormflow conditions, a smaller fraction of total flow is treated and (or) residence time in the treatment system is reduced.

Monitoring on the mainstem of Swatara Creek indicates watershed-scale effects owing primarily to changes in mining practices and secondarily to watershed-wide implementation of treatment systems. Most underground mines in the Swatara Creek Basin were abandoned before 1960 and are presently flooded. Drainage from these mines contributes substantially to baseflow in Swatara Creek. For Swatara Creek at Ravine, Pa., which is immediately downstream of the mined area, long-term data collected since 1959 indicate sulfate concentration declined from about 150 mg/L in 1959 to 75 mg/L in 1999; pH increased sharply from 3.5-4.4 (median ~4) to 4.6-7.0 (median ~6) after 1975. These trends resulted from a decline in pyrite oxidation and the onset of carbonate buffering. Because these long-term attenuation processes have had such a pronounced effect on water quality in Swatara Creek, the effects of recent implementation of limestone treatments are difficult to detect at a watershed scale. Nevertheless, during ecological surveys prior to 1991, no fish were found in Swatara Creek at Ravine. Only six species of fish were found in 1994 and 1996. However, increasing numbers of fish have been found annually since 1996. In 1999, 21 species of fish were documented.

The recent monitoring on the mainstem of Swatara Creek indicates the limestone treatments mitigate extreme fluctuations in pH during storm events; however, additional buffering capacity is needed to maintain near-neutral pH of Swatara Creek during large storm events. Concentration

1 of sulfate, specific conductance, and pH are inversely related to streamflow at Ravine, indicating dilution and acidification during stormflow. Declines in stream-water pH to values approaching 5.0 could result in the remobilization of adsorbed or precipitated metals associated with sediments; declines in pH below 5.0 could cause injury to aquatic organisms. Generally, to maintain stream pH during storms, additional or larger limestone diversion wells could be constructed to begin or increase alkalinity production as the stream stage rises and/or additional or larger limestone drains could be constructed to produce greater amounts of alkalinity and enhance the buffering capacity of baseflow.

Introduction The Pennsylvania Anthracite region consists of four large coalfields within an area of about 3,400 mi2 (8,850 km2) in the Appalachian Mountain Section of the Ridge and Valley Physiographic Province in eastern Pennsylvania (Wood et al. 1986; Berg et al. 1989; Way 2000). The coalfields are the sites of a series of parallel, moderately to deeply downwarped synclinoria. Most mines in the region were developed to access multiple coalbeds of the Llewellyn and Pottsville Formations of Pennsylvanian Age. In the Southern Anthracite Field, a total of 38 coalbeds with average thicknesses ranging from 1 to 11.5 ft (0.3 to 2.5 m) have been identified and mined to depths exceeding 3,280 ft (1,000 m) (Wood et al. 1968, 1986).

More than 150 years of mining in the Southern Anthracite Field has adversely affected surface- water and ground-water supplies (Growitz et al. 1985; Wood 1996; Brady et al. 1998; Way 2000). For example, losses of stream water to and contaminated drainage from abandoned anthracite mines within the upper 43 mi2 (112 km2) of the 576-mi2 (1,492-km2) Swatara Creek Basin degrade the aquatic ecosystem and impair uses of Swatara Creek to its mouth on the 59 mi (95 km) downstream from the mined area (fig. 1). Consequently, the Swatara Creek Basin was designated a “high priority watershed” for reducing nonpoint-source pollution (Pennsylvania Department of Environmental Protection 1998). Current land use in the upper 43-mi2 (112-km2) area, upstream from Ravine, Pa., is classified as 86.6 percent forested, 4.9 percent agricultural, and only 6.4 percent barren, mined (U.S. Geological Survey 2000). The land-use classification for this extensively mined area is misleading, however, because underground mines extend beneath much of the surface and “natural” reforestation conceals large tracts of unreclaimed spoil. Downstream from the mined area, forested and agricultural land uses predominate. For example, land use in the 116-mi2 (300-km2) area of the Swatara Creek Basin upstream from Pine Grove (fig. 1), inclusive of the area above Ravine, is classified as 69.7 percent forested, 25.0 percent agricultural, and 2.4 percent barren, mined.

Although several surface and underground anthracite mines presently are active, most mines in the Swatara Creek Basin were abandoned before 1960 (McCarren et al. 1964; Stuart et al. 1967). Barren, steep banks of spoil and culm and fine coal debris in siltation basins are sources of sediment (suspended solids) and dissolved ions in water that infiltrates or runs off the surface during storms (e.g. Olyphant et al. 1991). The abandoned underground mines have flooded and have collapsed locally causing subsidence. Surface flow is diverted through subsidence pits, fractures, and mine openings to the underground mines where the water becomes contaminated (Skelly & Loy, Inc. 1987; Ladwig et al. 1984; Growitz et al. 1985; Wood 1996). In downstream reaches, the contaminated water resurges as “acidic” or “abandoned” mine drainage (AMD), contaminating Swatara Creek and its while contributing substantially to baseflow (Fishel 1988).

2 76o30’ 76o25’ 76o20’ 76o15’

PENNSYLVANIA Southern Anthracite Coalfield

40o40’ Swatara Creek Watershed C1 C3 (above proposed dam site) Newtown . Cr Tremont ra ata Sw . r. Cr a C ck tar Bla wa E2-244 e S ittl o E2-0 Ravine r L 40 35’ pe E2-1 D1 Up

Pine . S Grove a Cr CH atar L U Sw E Y Little B L wer A KI D2 Lo NO L N L C CO o r. O . 40 30’ C . a ar at w S 0 2 4 6 8 10 MILES Proposed Swatara State 0 2 4 6 8 10 KILOMETERS Park Lake A

Anoxic Open NORTHERN SWATARA CREEK Limestone Swatara Limestone Channel USGS Monitoring Stations: Drain Cr. Oxic B0 Miscellaneous record B1 Limestone A1A2 Drain Continuous record A3 B2 Limestone E5-0 B3 Limestone C0 Sand C4 C1 Diversion E1-0 E5-1 Limestone C6 C2 Wells E1-1 E5-2 C3 E1-2 Diversion Martin E1-3 Well Run n C7 u E1-4 R . C8 l Cr . a e Cr C9 o l pring idd od S C M . Tracy Go Cr Airhole S2 ra M1-0 ata M1-1 Oxic L w o S M2-0 Limestone w S3 e M2-1 Drain r r. R k C a lac u B E3-S0 s Limestone Rowe c Tunnel E3-S1 h Diversion Wells C E2-244 r & Sludge Tanks E3-1 . Constructed E3-2 E2-0 E2-1 Wetlands E2-2 S4 Limestone Lorberry C Sand r. D1 B

Figure 1. Locations of water-quality and streamflow monitoring sites in the Swatara Creek Basin, Lebanon and Schuylkill Counties, Pennsylvania: A, continuous monitoring sites on Swatara Creek above the proposed dam for Swatara State Park Reservoir; B, monitoring sites within the Southern Anthracite Coalfield, above Ravine (area denoted in A).

3 2- Coal-mine drainage commonly has elevated concentrations of sulfate (SO4 ), (Fe), manganese (Mn2+), aluminum (Al3+), and other solutes that result from the oxidation of pyrite (FeS2) and the subsequent dissolution of carbonate, oxide, and aluminosilicate minerals by acidic water (Cravotta 1994; Rose and Cravotta 1998). Pyrite oxidation (reaction 1) takes place primarily in the unsaturated zone and at the land surface, where oxygen (O2) and moisture are available and where (H+) and other oxidation products tend to concentrate in fluids and solids.

2+ 2- + FeS2 + 3.5 O2 + H2O ® Fe + 2 SO4 + 2 H (1)

2+ + 3+ Fe + 0.25 O2 + H ® Fe + 0.5 H2O (2)

3+ + Fe + 3 H2O ® Fe(OH)3 + 3 H (3)

Infiltrating water or surface runoff can dissolve and transport the acidic oxidation products. In 2- contrast with SO4 , which is transported primarily as a dissolved ion, Fe can be transported as ferrous (Fe2+) and ferric (Fe3+) ions and as suspended Fe(III) solids (Cravotta 2000). In the 2+ 3+ presence of O2, Fe tends to oxidize to Fe (reaction 2) (Stumm and Morgan 1996). At pH >3, concentrations of Fe3+ are limited by the formation of Fe(III) oxyhydroxides and related solids (reaction 3) (Bigham et al. 1996).

The transport of dissolved Fe, Al, and trace metals typically is attenuated owing to precipitation and adsorption (Bigham et al. 1996; Smith et al. 1998; Webster et al. 1998; Cravotta and Trahan 1999). Attenuation generally is most effective as pH approaches neutrality (pH 6-7). However, at near-neutral pH and under anoxic conditions, concentrations of Fe2+ and Mn2+ can be elevated owing to relatively high solubility of Fe(II) and Mn(II) oxyhydroxides and carbonates. Aeration of water containing Fe2+ and Mn2+ can promote oxidation and hydrolysis, producing Fe(III) and Mn(III-IV) oxyhydroxides and 2 moles H+ for each mole Fe2+ and Mn2+. The potential for the production of H+ (or consumption of OH-) by hydrolysis reactions involving Fe, Mn, Al, and other metal ions is measured as acidity (Rose and Cravotta 1998).

The acid produced by pyrite oxidation or by hydrolysis can be neutralized by reaction with calcite (CaCO3) and dolomite [CaMg(CO3)2].

+ 2+ - CaCO3 + H = Ca + HCO3 (4)

+ 2+ 2+ - CaMg(CO3)2 + 2 H = Ca + Mg + 2 HCO3 (5)

These calcareous minerals are the dominant components of limestone and can occur in nodules, cementing agents, or fractures in sandstone, siltstone, shale, and associated strata of coal- - bearing rocks. Alkalinity, represented by bicarbonate (HCO3 ), and base cations including (Ca2+) and (Mg2+) are common products of neutralization by these calcareous minerals. Where absent or deficient at a mine site, the addition of calcite, dolomite, or other alkalinity-producing materials to mine spoil or mine drainage can be effective for prevention or neutralization of AMD and the attenuation of metals transport.

Where reclamation of a mine or mining-related surface effects is not possible, treatment of AMD may be necessary to neutralize acidity and remove dissolved and suspended metals from the

4 hydrologic system. Conventional treatment of acidic drainage involves the addition of strong alkaline chemicals to increase pH, which causes Fe, Al, Mn, and other metals to hydrolyze and precipitate (Skousen et al. 1998). Although effective, chemical treatment is expensive and funds to mitigate AMD are limited. Alternative treatment methods for AMD include wetlands and limestone-based systems (Hedin et al. 1994a; Skousen et al. 1998). These “passive systems” generally are limited by slower rates of neutralization and pollutant removal than for conventional treatments but can be cost effective where water chemistry meets suggested criteria and land and component materials are locally available (Skousen et al. 1998). Generally, if acidity exceeds alkalinity, limestone-based treatments could be appropriate to add alkalinity. If alkalinity exceeds acidity, oxidation ponds or aerobic wetlands are useful to remove metals as solids.

A variety of passive to semi-passive limestone treatment systems recently was installed at selected locations to neutralize the AMD and restore the aquatic ecosystem in approximately 25 mi (40 km) within the northern Swatara Creek watershed above Ravine, Pa. (figs. 1 and 2). The treatments, which include limestone-sand dosing, open limestone channels, anoxic and oxic limestone drains, limestone diversion wells, and limestone-based wetlands, were implemented mainly during fall 1996, spring 1997, and winter 1998 by the Schuylkill County Conservation District and the Northern Swatara Creek Watershed Association, with technical assistance from the U.S. Geological Survey (USGS) and the Pennsylvania Department of Environmental Protection (PaDEP). Additional treatments have been implemented since 1998 or are planned as funds can be obtained. A. Anoxic/Oxic Limestone Drain C. Limestone Diversion Well Monitoring pipe Acidic Inflow Diversion Dam

9 ft Neutral Acidic Outflow 5 ft Inflows Cleanout Pipe Diversion Well (All piping within Limestone Limestone, 1.25 - 4 inch drain is perforated) Neutral Outflow Fragments Limestone, < 2 inch (< 1 inch) Landfill liner 1 ton weekly (total 400 tons = 1.5 tons/yd3 x 270 yd3)

B. Open Limestone Channel D. Limestone-Compost Based Wetland

Acidic Inflow Neutral Outflow? 1.5 acre t f 1.5 acre 0 2

20 ft 110 ft (1 ft. deep)

Limestone, 1.25 - 4 inch (66 tons) Limestone sand, < 0.2 inch (44 tons) Mushroom compost Limestone, 8 - 16 inch (R4)

Figure 2. Schematic illustration of limestone treatment system designs: A, open limestone channel; B, limestone diversion well; C, limestone drain; D, limestone-compost wetland.

5 Open limestone channels and limestone-sand dosing generally are simple treatment systems where limestone is added directly to the stream channel (Ziemkiewicz et al. 1997; Skousen et al. 1998). An open limestone channel was constructed within a 110-ft (33.5 m) long segment of Swatara Creek at site B2 (fig. 1) on March 21, 1997. A total of 44 tons (40 Mg) of sand-size fragments (<0.2 inch; <0.5 cm) and 70 tons (63 Mg) of larger fragments (1.25-4 inches; 3-11 cm) were installed as a series of alternating berms extending part way across the 15-ft- (4.6-m) wide channel from opposite sides of the stream (fig. 2). Limestone sand, which can dissolve rapidly because of its small diameter, also was dumped into Coal Run (14 tons; 12.6 Mg) between sites C4 and C6 on September 4, 1996, and into Lorberry Creek (150 tons; 136 Mg) below site E2 on February 13-14, 1997 (fig. 1).

A limestone drain is another relatively simple treatment method that involves the burial of coarse limestone in air-tight trenches that intercept acidic water (Hedin et al. 1994a, b; Skousen et al. 1998; Cravotta and Trahan 1999). Keeping within the drain can enhance limestone dissolution and alkalinity production (Cravotta and Trahan 1999). Keeping O2 out of contact with the discharge water minimizes the potential for oxidation of Fe2+ and the consequent precipitation of Fe(III) encrustation, or armoring, on the limestone surfaces, while allowing O2 into the drain can promote oxidation and hydrolysis reactions and the removal of Fe, Mn, and trace metals (Cravotta and Trahan 1999). Limestone drains designed for varying flow rates and chemistry were constructed on March 15, 1995, at site E3-S0 to treat a small acidic discharge (10-30 gal/min (38-113 L/s), oxic inflow; 44 tons (40 Mg) limestone) along Lower Rausch Creek; on May 21, 1997, at site A1 to treat a large, anoxic discharge (50-200 gal/min (189-756 L/s); 400 tons (363 Mg) limestone) at the headwaters of Swatara Creek; and on June 22, 2000, at site C0 to treat a large oxic discharge (100-500 gal/min (378-1,890 L/s); 800 tons (727 Mg) limestone) near the headwaters of Swatara Creek (figs. 1 and 2). The design of the larger two systems was based on results for the smaller system where pH increased from 3.5 to 6.0-6.5 through the drain during 1-3 hour residence time (Cravotta and Trahan 1999).

In a limestone diversion well (fig. 2), acidic water is diverted from upstream points, and the hydraulic force of the piped flow is deflected upward through limestone fragments inside 4-ft (1.2- m) diameter “wells” (Arnold 1991). Approximately 0.5-1 ton of limestone is consumed weekly by each operating diversion well. Hydraulic churning abrades limestone to fine particles and prevents encrustation of the limestone by Fe(III) or Al oxyhydroxides. Dissolution of limestone within and downstream of the diversion wells promotes pH increases. In addition to pulverized limestone, Fe(III) and Al oxyhydroxides precipitate and accumulate downstream of the diversion wells. On November 14, 1995, a pair of diversion wells was installed to treat water diverted from Swatara Creek at site C2; on July 13, 1997, a single diversion well was installed to treat water from Martin Run at site C8; and on December 15, 1998 a pair of diversion wells was installed to treat water from Lorberry Creek above site E2-0 (fig. 1).

Constructed wetlands for treatment of mine drainage can attenuate the transport of dissolved and suspended pollutants by promoting the production of alkalinity and the precipitation and deposition of iron and other metals (Hedin et al. 1994a; Skousen et al. 1998). For net acidic water (acidity > alkalinity), wetlands that have compost and/or limestone substrates can be appropriate. The organic matter in the compost provides a substrate for plant rooting and for microbial reduction of SO4. On December 1997, near the mouth of Lower Rausch Creek at site E3 (figs. 1, 2), a 3-acre (1.2 ha) compost-limestone based wetland was constructed to remove metals from streamflow that commonly had near-neutral pH but had potential to be net acidic if

6 untreated during stormflow conditions. The Lower Rausch Creek wetlands were constructed downstream from the outflow of the “oxic limestone drain” (OLD) built in 1995 (fig. 1). Although the OLD was effective for neutralizing acid and converting dissolved metals to solid forms as described in another report (Cravotta and Trahan 1999), a settling basin or wetland was needed to attenuate the transport of suspended metals from the OLD and other sites in the Lower Rausch Creek watershed.

In general, passive-treatment systems are designed to be effective for the typical baseflow water-quality conditions. Despite recent documentation of case studies for passive systems in the Northern Appalachian Region (Skousen et al. 1998), hydrological and geochemical factors affecting the performance of passive limestone-treatment systems are poorly characterized, particularly for high-flow conditions. Each treatment has different advantages and disadvantages; however, all suffer from possible complications associated with variability of flow rates and chemistry of the AMD-contaminated water and from uncertainties about efficiency and longevity of the treatment. Furthermore, every site requiring treatment has unique characteristics.

This paper characterizes site conditions and evaluates the effectiveness of limestone treatments for neutralizing acidity and removing metals and other pollutants from AMD and affected stream water in the Swatara Creek watershed. Data collected by USGS primarily during June 1996 through June 2000 for sites within the Swatara Creek watershed above Pine Grove, Pa., are used for this evaluation. Additionally, historical USGS data for Swatara Creek at Ravine, Pa., collected periodically since 1959 are used to indicate long-term trends in water quality.

Methods To characterize untreated AMD, treatment-system performance, and cumulative downstream effects over a range of environmental conditions, the USGS established monitoring sites upstream and downstream of each treatment and along lower reaches of Swatara Creek (fig. 1). During base-flow and high-flow conditions in 1995-2000, data on streamflow rate and water quality were collected at these sites.

Three sites on Swatara Creek, site C3 near Newtown (station 0157155014), site D1 near Ravine (station 01571820), and site D2 near Pine Grove (station 01572025), were equipped with automatic stage-recording, water-quality monitoring, and (or) water-sampling devices in 1996. Automatic water-quality monitoring devices also were installed in 1996 at site C1 (station 0157155010) above site C3, and in 1999 at sites E2-244 (station 403542076263201), E2-0 (station 01571774), and E2-2 (station 01571778) on Lorberry Creek. Where equipped, stream stage was measured continuously with a pressure transducer; temperature, pH, and specific conductance (SC) were measured with a multiparameter sonde. The stage and water-quality data were recorded at 15-minute intervals with an electronic data logger. Continuous streamflow was computed on the basis of a stage-discharge rating developed for each site (Rantz et al. 1982a, b). At these and other sites, instantaneous data for temperature, SC, dissolved O2 (DO) and pH also were measured when samples were collected or processed by use of calibrated instruments (Wilde et al. 1998).

At all sites, baseflow and stormflow samples were collected manually as grab samples from well- mixed zones in the stream or mine discharge. At sites equipped with automatic samplers (sites C3, D1, and D2), discrete stormflow samples were collected using pumping samplers containing

7 24 1-L polyethylene bottles. The automatic samplers were programmed to initiate sampling after a specified rise in stream stage above the current base stage and to proceed at intervals based on rate of change in stage until the stream returned to the base stage. Samples submitted for analysis were selected to cover rising, peak, and falling stages of the storm hydrograph. Stormflow samples were analyzed for more than 20 events, including significant storms in October and December 1996, May 1997, January 1998, January and September 1999, and May 2000 (fig. 3).

2000

1000

500

200 s / 3 t f

,

W 100 O L F M A

E 50 R T S 20

10

5 M J J A S O N D J F M A M J J A S O N D J F M A M J J A S O N D J F M A M J J A S O N D J F M A M J J 1996 1997 1998 1999 2000 Figure 3. Streamflow hydrograph for Swatara Creek near Ravine, Pa., June 1996 - June 2000. Square symbols indicate streamflow at times water-quality samples were collected.

Water samples were split into subsamples in the field or in the USGS laboratory in Lemoyne, Pa., and stored in sample-rinsed polyethylene bottles at 4°C until analyzed. Samples for dissolved (0.45-mm filter) and total recoverable (whole-water; in-bottle nitric and hydrochloric acid digestion) metal analysis were stored in acid-rinsed polyethylene bottles and acidified with nitric acid (HNO3). Samples were analyzed for major ions, metals, and nutrients by inductively coupled plasma optical emission spectroscopy (ICP-OES), ion chromatography (IC), colorimetry, and electrometric titration at the PaDEP Bureau of Laboratories facility in Harrisburg, Pa., following methods of Greenberg et al. (1992), Hoffman et al. (1996), and Fishman and Friedman (1989). Most unknown samples met ionic charge-balance criteria (+10%) and had dissolved concentrations less than, or equal to, total concentrations. Furthermore, duplicates, blanks, and standard reference samples that were submitted periodically with batches of unknown samples typically indicated precise and accurate measurements. The water-quality and streamflow data, which are maintained in the USGS National Water Information System (NWIS) data base, were published annually (Durlin and Schaffstall 1998, 1999, 2000).

8 Fish were collected annually at the Ravine site by electrofishing over a 500-ft (150-m) reach consisting of mixed riffle, run, and pool habitats as described by Bilger et al. (1999). Individual fish were collected with dip nets, identified, and measured before releasing most specimens. Some specimens were sacrificed for analysis of metals in fish tissue as reported by Cravotta and Bilger (2001).

Results and Discussion

Characterization of Stream-Water Quality Prior to Treatment Historical and current data from 1959 to 2000 indicate progressive improvement in water quality for Swatara Creek at Ravine. Baseflow samples were collected for comparable streamflow conditions for the historical and current record. Baseflow SO4 declined from a median of about 150 mg/L in 1959 to 75 mg/L in 1999; pH increased sharply from 3.5-4.4 (median ~4) to 4.6-7.0 (median ~6) after 1975 (fig. 4).The decline in SO4 concentration probably was caused by a decline in pyrite oxidation after flooding of the abandoned mines had minimized inflows of oxygenated air and water. The associated increase in pH was caused by the onset of carbonate buffering which occurred when the rate of alkalinity production equalled or exceeded acid production (Cravotta et al. 1999). Although a variety of environmental factors could affect pH and SO4 concentrations, consistently near-neutral pH values with variable SO4 concentrations at Ravine during 1998-2000 (fig. 4) imply that the recently implemented limestone treatments have neutralized acid, further improving water quality.

300 8 A B

7 ) L /

g 200 m

( 6

E H T p A

F 5 L 100 U S 4

0 3 1957 1968 1979 1990 2001 1957 1968 1979 1990 2001 YEAR YEAR

Figure 4. Long-term baseflow water-quality trends for Swatara Creek at Ravine, Pa.: A, Sulfate; B, pH. Data from McCarren et al. (1964), Stuart et al. (1967), Skelly & Loy, Inc. (1987), Fishel (1988), and Durlin and Schaffstall (1998, 1999, 2000).

As a consequence of the improved water quality, the fish community has rebounded (fig. 5). Fish were nonexistent in Swatara Creek at Ravine during ecological surveys prior to 1991 (fig. 5). However, in 1994 and 1996, six species of fish were found. Increasing numbers of fish species have been documented annually since 1996; 21 species were found in 1999 (fig. 5). A similar, though less dramatic increase in benthic macroinvertebrate species is also occurring. In 1995, no macroinvertebrates were found at Ravine; in 1999, 16 taxa were found (Cravotta and Bilger 2001).

9 Y 24 During the 1996-2000 study period, stream L 9 L

9 water of Swatara Creek at Newtown (site C1 in A 22 9 n 1 U

a fig. 1) and Ravine (site D1 in fig. 1) was mildly N g N 20 e 8 acidic to near-neutral (pH 4.5 - 8.0) with A B

9

t 9 D

n moderate concentrations of dissolved solids 18 1 E e 7 T

m 9 (SC 60 - 400 mS/cm). The pH, SC, and t N 16 9 a U 1

e concentration of SO were correlated and

r 4 O T C 14

inversely related to streamflow (figs. 6, 7A). D S M E

I 12 Higher values of pH, SC, and SO4 were A C

E associated with baseflow conditions sustained 10 P

S by near-neutral AMD in the watershed, such as

4 6 H 8 9 9 the Tracy Airhole, Colket, and Marshfield S I 9 9 1 F 1 6 discharges (Durlin and Schaffstall 1998, 1999, F 0 O 9 2000). Lower values of pH, SC, and SO were 4 4 9 R 1

, E associated with stormflow (Cravotta 2000; 5 B 2 8

M Cravotta and Bilger 2001). Mixing of baseflow 9 1 U 0

N with acidic rainfall and storm runoff having low pH and low dissolved solids explains the major Figure 5. Recent trends in fish species stormflow characteristics of declining pH, SC, abundance for Swatara Creek at Ravine, Pa. Year of survey indicated above bar; data and SO4 with increased streamflow for these from Cravotta and Bilger (2001). reaches of Swatara Creek (Cravotta 2000).I

950 8.5 900 Swatara Creek at Ravine (D1) s / 3

t pH f 850

, 8.0

W 800

O STREAMFLOW L 750 F

M 700 7.5 A E

R 650 T S

600

; 7.0 m

c 550 / S

u 500

, 6.5 H p E 450 C

N 400

A SPECIFIC T CONDUCTANCE 6.0

C 350 U

D 300 N

O 250 5.5 C

C 200 I F I 150 C

E 5.0 100 P S 50 0 4.5 M A M J J A S O N D J F M A M J J 1999 2000

Figure 6. Continuous data for streamflow, specific conductance (SC), and pH, Swatara Creek at Ravine, Pa., March 1999 - July 2000. Black diamonds indicate SC of automatically or manually collected grab samples.

10 420 8.2 s

/ 8.0

3 400 Swatara Creek ab diversion wells (C1) t f

A 7.8 , 380 7.6 W 360 O 7.4 L

F 340 7.2

M SPECIFIC CONDUCTANCE 320 7.0 A

E 300 6.8 R

T 6.6 280 S

pH 6.4 ; 260 m 6.2 c / 240 6.0 S u

220 5.8 , H E 5.6 200 p C

N 180 5.4 A 5.2 T 160 C 5.0 U 140 4.8 D

N 120 4.6 O 4.4 C 100

C 4.2

I 80

F 4.0 I 60 C STREAMFLOW 3.8 E 40 P 3.6 S 20 3.4 0 3.2 M A M J J A S O N D J F M A M J J 1999 2000 8.2 8.0 Swatara Creek bl diversion wells (C3) B 7.8 7.6 7.4 7.2 7.0 6.8 6.6 6.4

H 6.2 p 6.0 5.8 5.6 5.4 5.2 5.0 4.8 4.6 Swatara Creek 4.4 ab diversion wells (C1) 4.2 M A M J J A S O N D J F M A M J J 1999 2000

Figure 7. Continuous data for streamflow, SC, and pH, Swatara Creek at Newtown, Pa., March 1999 - July 2000: A, Streamflow, pH, and SC, upstream of diversion wells (black diamonds indicate SC of manually collected grab samples); B, pH upstream and downstream of diversion wells.

11 950 7.2 Lorberry Creek bl Rowe Drainage Tunnel, ab diversion wells (E2-244) 900 32 A 7.0 pH DISCHARGE 6.8 850 30 SPECIFIC RATE

) CONDUCTANCE 6.6 s

/ 800 28 3 t

f 6.4 ( 750

E 26 6.2 G

R 700 A 24 6.0 H

C 650

S 22 5.8 I D

, 600 ) 20 5.6 m c

/ 550 5.4

S 18 u ( H 500 5.2 p E 16 C

N 450 5.0 A

T 14 4.8 C 400 U

D 12 4.6

N 350

O 4.4

C 10 300 C I 4.2 F I 250 8 C 4.0 E

P 200 6 S 3.8 150 4 3.6 100 2 3.4 50 3.2 M A M J J A S O N D J F M A M J J 1999 2000 7.2 Lorberry Cr bl diversion wells (E2-0) 7.0 B 6.8 6.6 6.4 6.2 6.0 5.8 5.6 5.4 H

p 5.2 5.0 4.8 4.6 4.4 4.2 4.0 3.8 3.6 Lorberry Cr at downstream gage (E2-1) 3.4 3.2 Lorberry Cr ab diversion wells (E2-244) M A M J J A S O N D J F M A M J J 1999 2000 Figure 8. Continuous data for SC and pH, Lorberry Creek above Lorberry Junction, Pa., March 1999 - July 2000: A, pH and SC, upstream of diversion wells (black diamonds indicate SC of manually collected grab samples); B, pH upstream and downstream of diversion wells.

12 In contrast with data for Swatara Creek at Newtown and Ravine, the pH and SC for headwaters of Lorberry Creek, below the Rowe Drainage Tunnel (site E2-244 in fig. 1), were inversely correlated and varied more widely than those for Swatara Creek sites (fig. 8A). Decreased pH and increased SC resulted from additions of acidic water containing elevated concentrations of SO4 and other dissolved ions. Two processes contributed to the temporal fluctuations in pH and SC of Lorberry Creek. Although the Rowe Drainage Tunnel drains the abandoned, flooded New Lincoln Colliery complex, an underground mine that is currently active below the complex periodically pumps untreated, acidic water to the New Lincoln mine pool. On the shorter scale of a few days, periodic fluctuations in pH and SC result from the addition of this untreated AMD to the water drained from the abandoned mine complex. When the pumping was active, pH of Lorberry Creek declined by 0.5 to 1 unit while SC increased 50 to 200 mS/cm. These short-term fluctuations in pH and SC were apparent for the continuous monitoring data but generally were not captured in the data for monthly grab samples (fig. 8A). On the broader time scale, pH declined and SC increased during a drought from May through August 1999 followed by two tropical storms in September 1999. Prior to this investigation, the Rowe Drainage Tunnel was characterized as near-neutral to moderately acidic AMD (pH 5-6; SC <400 mS/cm) requiring only the removal of particulate iron. As illustrated by the continuous monitoring data, however, the mine discharge was extremely acidic for sustained periods lasting hours to weeks (pH 3-4; SC 500-950 mS/cm) during 1999-2000.

Declines in stream-water pH to values approaching 5.0, as for extreme events during 1996-2000, could cause injury to aquatic organisms (e.g. Baker and Schofield 1982; Earle and Callaghan 1998) and could result in the remobilization of adsorbed or precipitated metals associated with streambed sediments (e.g. Francis et al. 1989). Solid forms of the metals, as particulate and particle coatings, can be transported during storms and ingested and accumulated by aquatic organisms (e.g. Elder 1988; Cravotta and Bilger 2001). Hence, pH adjustment as well as removal of metals from the hydrologic system would be needed to meet water-quality goals.

Evaluation of Treatment Effects EXPLANATION OF BOXPLOT The pretreatment and post-treatment data for (2) Number of observations pH, net alkalinity (alkalinity - acidity), dissolved Outlier data value more than 3 times the interquartile range iron, and dissolved aluminum at sites upstream outside the quartile and downstream of six different treatment Outlier data value less than or equal to 3 and more than 1.5 times the interquartile range outside the quartile systems in the Swatara Creek watershed are illustrated as boxplots (fig. 9). The boxplots show Data value less than or equal to 1.5 times the interquartile range outside the quartile data for all samples. Table 1 indicates results of 75th percentile tests for difference between paired samples (same date) from the upstream and downstream Median sites. Differences between the pairs were 25th percentile evaluated for the range of flow and for low-, normal-, and high-flow classes. If streamflow at Ravine on the date of sampling was less than the Figure 9A. Boxplots of pre- and post- 25th percentile for the study period, the sample implementation water-quality data for grab was classified “low-flow”; between the 25th and samples upstream and downstream of 75th percentiles, the sample was classified selected treatment systems in Swatara “normal-flow”; or greater than the 75th percentile, Creek Basin. Explanation, only (boxplots the sample was classified “high-flow” (fig. 10). showing data follow).

13 ANOXIC LIMESTONE DRAIN AT SWATARA CRK HEADWATERS (A2 - A3) 8.5 225 28 0.8 MINE.DISCH DNSTRM MINE.DISCH DNSTRM 27 MINE.DISCH DNSTRM MINE.DISCH DNSTRM )

3 26 200 (33) (33)

O (10) 8.0 25 0.7 C 24 (34) a 23 C 175 7.5 (37) 22 L S /

A 21 0.6

G

150 L L / 20 / M

7.0 G , G 19 (4) M M M

18

125 , U , 17 0.5 N N 6.5 D I I

O 16 C M 100 R 15 I A U

H -

14 L 0.4

6.0 D p A K E 13 75 L

V (10) 12 D A L E (

5.5 11 O V 0.3

Y 50 (32) S 10 L (33) T S I O

I 9 N (12) S D 5.0 I 8

25 S L 7 (33) I 0.2 A D (4) K 6 4.5 L 0 5 (33) A (4) (9) 4 (4) T 3 0.1 4.0 E -25

N 2 1 3.5 -50 0 0 LIMESTONE DIVERSION WELLS NR SWATARA CRK HEADWATERS (C1 - C3) 8.5 30 1.3 6.0 UPSTRM DNSTRM UPSTRM DNSTRM UPSTRM DNSTRM UPSTRM DNSTRM ) (87) 3 1.2 5.5

8.0 (93) O (92) C

a 1.1 5.0 C 7.5 20 (92) L S

1.0 /

A 4.5

G L L

/ (40) / M

7.0 0.9 , G G 4.0 M M M

(43)

, U , 0.8

6.5 10 N N D I I 3.5 O C M

R 0.7 I A U

(38)

H - 6.0 L 3.0 p

D A K E 0.6 L V D A

L 2.5 E 5.5 ( 0

O 0.5 V Y L S T

I S 2.0

O (1) I 0.4 N 5.0 S I D S L I 1.5 A

0.3 D 4.5 K -10 (40) (1) L

A 1.0

0.2 (1) T

4.0 E (1) 0.5

N 0.1 (0) (0) (0) (0) 3.5 -20 0 0 LIMESTONE DIVERSION WELLS AT LORBERRY CRK HEADWATERS (E2-244 - E2-0) 8.5 70 18 2.4 MINE.DISCH DNSTRM MINE.DISCH DNSTRM MINE.DISCH DNSTRM 2.3 MINE.DISCH DNSTRM ) 17 3 60 (14) 2.2

8.0 O 16 (14) 2.1 (13) C

a 15 2.0

C 50 7.5 1.9 L S (7) 14 / A 1.8

G

40 L

L 13 / / M 1.7

7.0 , G (14) (14) G 12 1.6

(17) M M (19) M

30 , U , (15) (12) 11 1.5 N N 6.5 D (13) (14) I I (7) 1.4 O C 20 10 (7) M 1.3 R A I U

H - 6.0 9 L 1.2 p D A K

10 E 1.1

L 8 (7) V D

A 1.0 L E 5.5 ( 7 O V 0.9

Y 0 L S T

I 0.8 S

6 O I N 5.0 S I

D 0.7

-10 5 S L I 0.6 A 4 D 4.5 K 0.5 L -20

A 3 0.4

T 0.3

E 2 4.0 -30 0.2 N 1 0.1 3.5 -40 0 0

Figure 9B. Boxplots of pre- and post-implementation water-quality data for grab samples upstream and downstream of selected treatment systems in Swatara Creek Basin. Pre- implementation data shaded. See explanation of boxplot in figure 10A.

14 OPEN LIMESTONE CHANNEL AT SWATARA CRK HEADWATERS (B1 - B3) 8.5 30 2.8 1.0 DNSTRM UPSTRM DNSTRM 2.7 UPSTRM DNSTRM UPSTRM DNSTRM UPSTRM ) (6) 3 2.6

8.0 O 25 (32) (32) 2.5 0.9 C 2.4

a (6)

(34) C 2.3 7.5 (35) 20 0.8 L

S 2.2 /

A 2.1 (33)

G L L / /

2.0 M

7.0 15 0.7 ,

G (33) G 1.9 M M M

1.8

(33) , U , 1.7

10 N 0.6 N 6.5 D (33) I I

O 1.6 C M

R 1.5 A I U

H - 5 1.4 L 0.5 6.0 p D A K

E 1.3 L V 1.2 D A L E ( 0 0.4 5.5 1.1 O V Y L S 1.0 T I S O

I 0.9 N -5 S 0.3

5.0 I D 0.8 S L 0.7 I A D (6) K -10 (6) 0.6 0.2 4.5 L 0.5 A 0.4 T

4.0 E -15 0.3 0.1

N 0.2 (0) (0) 0.1 (0) (0) 3.5 -20 0 0 LIMESTONE SAND AT CARBON RUN HEADWATERS (C4 - C6) 8.5 35 2.8 0.60 UPSTRM DNSTRM UPSTRM DNSTRM 2.7 UPSTRM DNSTRM UPSTRM DNSTRM )

3 2.6 8.0 0.55 O 30 2.5 (8)

C (8) 2.4 a 0.50

C 2.3 (8) 7.5 25 L S 2.2 (8) /

A 2.1 0.45

(9) G L L /

/ 2.0 (3) 7.0 M (9) 20 (8) (8) , G G 1.9 0.40 M

M (3) M

1.8

(4) , U 6.5 , 1.7

15 N N D I I 0.35

O 1.6 C M

R 1.5 A I U H

-

6.0 L p 10 1.4 0.30

D A K

(3) E 1.3 L V 1.2 D

A 0.25 L E 5.5 ( 5 1.1 O V Y L

S 1.0 T

I 0.20 S O

I 0.9 N 5.0 0 S I D 0.8 S L 0.7 I 0.15 A D

4.5 K -5 0.6 L 0.5

A 0.10 0.4 T

E 0.3 4.0 -10 0.05 (0) N (0) 0.2 (0) 0.1 (0) 3.5 -15 0 0 LIMESTONE-COMPOST WETLANDS AT LOWER RAUSCH CRK (E3-1 - E3-2) 10.0 200 6.0 1.7 UPSTRM DNSTRM UPSTRM DNSTRM UPSTRM DNSTRM UPSTRM DNSTRM ) 1.6 9.5 3 180 (6) 5.5 (14) (15) O (14) 1.5 C

9.0 a (15) 5.0 1.4 C 160

L S 8.5 / 1.3 A 4.5

G

140 L L / / M 1.2 8.0 , G G 4.0 M M

M 1.1

120 , U 7.5 (6) , N N D 1.0 I (1) I 3.5 O C 100 M 7.0 R A I 0.9 U

H - 3.0 L (14) p

D

A 0.8 K

6.5 80 E

L (5) V D

A (13) L 2.5 0.7

(14) E (

O 6.0 V

Y 60 L (13) S 0.6 T I S

2.0 O 5.5 I N S

I 0.5 D

40 S L 1.5 I A 5.0 D 0.4

K (1) L 20 (1) 0.3

A 1.0 4.5 (5) T 0.2

E (1) 0 0.5 4.0 N 0.1 3.5 -20 0 0

Figure 9C. Boxplots of pre- and post-implementation water-quality data for grab samples upstream and downstream of selected treatment systems in Swatara Creek Basin. Pre-implementation data shaded. See explanation of boxplot in figure 10A.

15 16 site, Solids .31 .31 .23 .17 .35 .00 (>) .50 .22 .50 .42 .03 (<) .18 .36 .53 .48 .94 n.p. n.p. n.p. n.p. 1.00 1.00 0.16 n.p. Suspended Total .83 .04 (<) .23 .10 (>) .01 (<) .00 (<) .18 .00 (<) .01 (<) .00 (<) .08 (<) .00 (<) .02 (<) .00 (<) .00 (<) .29 .00 (<) .44 .02 (>) .11 n.p. (<) 1.00 0.00 (<) n.p. Manganese .69 .11 .06 (>) .44 .03 (<) .21 .03 (<) .00 (<) .11 .46 .00 (<) .00 (<) .04 (<) .00 (<) .04 (<) .14 .05 (>) .00 (>) .33 Dissolved n.p. (<) n.p. n.p. (<) 0.00 (<) n.p. Total .00 (<) .07 (<) .02 (<) .39 .20 .65 .72 .00 (<) .03 (<) .00 (<) .46 .00 (<) .29 .07 (<) .21 (>) .20 .07 (>) .20 .14 .65 he same on the basis of the Wilcoxon n.p. (>) 0.61 n.p. (<) n.p. ons for less than 25th, 25th to 75th, and site was subtracted from that for the downstream site was subtracted from that for the downstream Aluminum 0 (see figure 10). .35 .00 (<) .00 (<) .02 (<) .00 (<) .06 (<) .14 .59 .65 .59 .04 (<) Dissolved n.p. (<) n.p. n.p. 1.00 1.00 n.p. 1.00 n.p. 1.00 1.00 1.00 n.p. 1.00 dence for rejection of Ho. If a p-value was not computed, dence for rejection of Ho. If a p-value was not < 2/3 total count and all differences are <0 or >0. < 2/3 total count and all differences are <0 or Total .68 .46 .00 (<) .02 (<) .01 (<) .00 (<) .00 (<) .02 (<) .00 (<) .00 (<) .06 (>) .00 (<) .00 (<) .00 (<) .09 (>) .78 .37 .14 n.p. 1.00 n.p. n.p. 0.00 (<) 1.00 Iron .35 .03 (<) .00 (<) .01 (<) .00 (<) .20 .00 (<) .00 (<) .00 (<) .00 (<) .00 (<) .08 (<) .97 .02 (<) .00 (<) .00 (<) .33 .17 .21 .27 Dissolved n.p. n.p. 0.00 (<) 1.00 Sulfate .11 .22 .22 .40 .04 (>) .67 .01 (<) .00 (<) .13 .40 .00 (>) .12 .00 (<) .03 (<) .22 .50 .22 Dissolved n.p. (<) n.p. n.p. n.p. n.p. 0.24 n.p. (<) .42 .60 .00 (>) .00 (>) .00 (>) .00 (>) .29 .00 (>) .86 .53 .35 .48 .03 (<) .00 (<) .00 (<) .12 .84 .12 .00 (>) .00 (>) Calcium Dissolved n.p. n.p. n.p. 0.00 (<) Net .00 (>) .36 .00 (>) .00 (>) .00 (>) .00 (>) .01 (>) .65 .14 .64 .86 .35 .00 (<) .00 (<) .00 (<) .14 .21 .09 (>) .27 Alkalinity n.p. (>). n.p. n.p. 0.00 (<) 1.00 respectively; n.p., p-value not computed] n.p., p-value respectively; pH .00 (>) .00 (<) .00 (>) .00 (>) .00 (>) .00 (>) .00 (>) .00 (>) .05 (>) .18 .60 .58 .93 .01 (>) .00 (>) .47 .73 .00 (>) .26 .60 n.p. n.p. n.p. 0.00 (>) .17 .13 .65 .04 (>) .02 (>) .03 (>) .01 (<) .00 (>) .06 (<) .98 .00 (>) .35 .00 (>) .06 (>) .00 (>) .00 (<) .01 (<) .00 (<) .02 (>) .50 .05 (>) .08 (>) Specific Conduct. Limestone Sand (LSC) implemented September 4, 1996 near Coal Run Headwaters (C4 - C6) Limestone Sand (LSC) implemented September n.p. 0.00 (<) Anoxic Limestone Drain (ALD) implemented May 21, 1997 near Swatara Creek Headwaters (A2 - A3) (ALD) implemented May 21, 1997 near Swatara Anoxic Limestone Drain ) that concentration or transport of constituents at downflow monitoring site is the same as that from the is the same as monitoring site at downflow of constituents or transport ) that concentration 1 Open Limestone Channel (OLC) implemented March 21, 1997 near Swatara Creek Headwaters (B1 - B3) (OLC) implemented March 21, 1997 near Open Limestone Channel Limestone Diversion Wells (LDZ) implemented November 11, 1995 near Swatara Creek Headwaters (C1 - C3) Limestone Diversion Wells (LDZ) implemented ture .28 .67 .92 .28 .95 .61 .75 .24 .00 (>) .92 .60 .33 .00 (>) .40 .64 .88 .22 upstream monitoring site over a range of hydrologic conditions, June 1996 - July 2000 June 1996 hydrologic conditions, over a range of monitoring site upstream Limestone Diversion Wells (LDL) implemented November 19, 1998 near Lorberry Creek Headwaters (E2-244 - E2-0) Limestone Diversion Wells (LDL) implemented Tempera- 1.00 n.p. n.p. 0.82 1.00 1.00 1.00 Limestone-Compost Base Wetlands (WLR) implemented December 15, 1998 at mouth of Lower Rausch Creek (E3-1 - E3-2) Limestone-Compost Base Wetlands (WLR) Rate .50 .01 (<) .01 (<) .04 (<) .09 (>) .01 (>) .01 (>) 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 0.00 (>) 1.00 1.00 1.00 1.00 Streamflow 2 Hydrologic condition for date of sampling at reference site on Swatara Creek at Ravine, Pa. Low-, normal-, and high-flow conditi Hydrologic condition for date of sampling at See figure 1 for location of sites and text for description and discussion of treatment. The p-value is the two-tailed probability that constituent values for the matched pairs of upstream and downstream samples are t The p-value is the two-tailed probability that 2 3 1 Condition Hydrologic Normal, n=6 Normal, n=6 High, n=4 Low, n=11 Normal, n=18 High, n=11 Low, n=2 High, n=5 Low, n=2 High, n=5 Normal, n=14 Low, n=2 Low, n=8 Normal, n=12 High, n=10 Low, n=9 High, n=10 Normal, n=10 [For p-values less than 0.1, “>” or “<” indicates mean rank at downflow site was greater than or less than that for the upstream or less than that was greater than at downflow site mean rank “>” or “<” indicates less than 0.1, [For p-values All samples, n=13 All samples, n=13 All samples, n=40 All samples, n=33 All samples, n=30 All samples, n=16 the p-value is shown as 1.0 if ties > 2/3 total count and all differences are <0 or >0 or it is shown as n.p. with trend if ties the p-value is shown as 1.0 if ties > 2/3 total matched-pair signed-rank test (P-STAT, Inc. 1989; Helsel and Hirsch 1992). For each pair, the constituent value for the upstream matched-pair signed-rank test (P-STAT, Inc. of the difference ranks is zero. The smaller the p-value, the stronger the evi site. The null hypothesis, Ho, is that the median for the period June 1996 - July 200 greater than 75th percentile values, respectively, of mean-daily flow frequency distribution Table 1. Probability (p-value Table 1.

16 10,000 Anoxic Limestone Drain (A2 - A3).) Before 5,000 construction of the anoxic limestone drain 2,000 (ALD) at site A2 near the headwaters of 1,000 500 Swatara Creek, acidic mine drainage flowed 200 PINE GROVE from site A2 for about 600 yd (550 m) to site A3 100 at Swatara Creek (fig. 1). The untreated AMD s

/ 50

3 was anoxic (DO <0.5 mg/L) and had elevated t f

, 20 RAVINE concentrations of dissolved Fe and Al (fig. 9B).

W 10

O 5 As this water became aerated flowing L

F downstream from sites A2 to A3, the pH and

M 2

A 1 Fe concentrations declined due to oxidation E NEWTOWN R 0.5 and hydrolysis of Fe (reactions 2 and 3); T

S 0.2 however, dissolved Al remained elevated (>0.3 0.1 mg/L) because of low pH (<4). After 0.05 construction of the ALD, the concentration of 0.02 0.01 dissolved Fe was unchanged at the ALD Low Normal High 0.005 outflow compared to the untreated mine 0.002 discharge, whereas the pH and net alkalinity 0.001 0 0.25 0.50 0.75 1.00 increased and dissolved Al concentration PROBABILITY decreased (<0.1 mg/L) (fig. 9B; mine discharge, pre- vs. post-implementation). As Figure 10. Probability of streamflow for the treated water flowed downstream from Swatara Creek at Newtown, Ravine, and sites A2 to A3, the pH increased while the net Pine Grove, Pa., June 1996 - July 2000. alkalinity and Fe concentrations decreased Vertical lines at 25th and 75th percentiles due to aeration and the resulting exsolution of divide “low”, “normal”, and “high” flow carbon dioxide and oxidation and hydrolysis of intervals. Fe. The result was significantly improved, near-neutral water quality at site A3. Even though high flows diluted the effects at the downstream site(s), the ALD treatment was effective for neutralization of AMD and attenuation of metals over the range of flow conditions (table 1).

Open Limestone Channel (B1 - B3). Before construction of the open limestone channel (OLC) near the headwaters of Swatara Creek, the stream water at sites B1 and B3 was acidic, with pH <4.5 but relatively low concentrations of Al and Fe (0.5-3 mg/L) for a mining-affected stream (fig. 9C). These conditions and preliminary field experiments indicating initially rapid increases in the pH of stream water in contact with limestone warranted the construction of the OLC. The preliminary experiments indicated the rate of limestone dissolution decreased with increased pH, which is consistent with other reports (e.g. Cravotta and Trahan 1999). However, the conditions of preliminary testing at the OLC were reflected for only one set of post-implementation samples before the ALD was constructed on an upstream . The ALD, which was constructed 2 months after the OLC was constructed, produced near-neutral stream water at site B1 above the OLC (fig. 9C; upstream, pre- vs. post-implementation). This near-neutral water was not aggressive toward limestone in the stream channel compared to acidic stream water. Hence, differences in pH and concentrations of net alkalinity, and dissolved Fe and Ca between the upstream site (B1) and downstream site (B3) at the OLC were not significant for most flow conditions (table 1, fig. 9C). Nevertheless, decreases in concentrations of dissolved SO4, Al, and Mn were significant for some conditions (table 1).

17 Limestone Diversion Wells (C1 - C3). Before installation of the limestone diversion wells near the headwaters of Swatara Creek (LDZ), the stream water at sites C1 and C3 had pH <4.5 and dissolved Al >1.5 mg/L (fig. 9B). After installation, the pH of stream water at site C3 150 yd (137 m) below the diversion wells usually was 1-2 units higher than that at site C1 150 yd (137 m) above the wells (figs. 7B, 10B). Considering intermittent data, only, the limestone diversion wells increased pH, net alkalinity, and calcium for all flow conditions and decreased iron and aluminum at low- to normal-flow conditions. (fig. 9B). However, during extreme stormflow conditions associated with tropical storms in September 1999 and spring storms in March-May 2000, the pH of downstream water (continuously monitored) was not effectively increased (fig. 7B). During such stormflow conditions, the effectiveness of the limestone diversion wells was diminished because a smaller proportion of total streamflow was treated and metallic deposits from the streambed could be resuspended and (or redissolved). Additional or larger diversion wells would be needed to treat these large stormflows.

Limestone Diversion Wells (E2-244 - E2-0). Before installation of the limestone diversion wells below Rowe Drainage Tunnel (LDL), Lorberry Creek at sites E2-244 and E2-0 had variable pH (4-7), but typically was net acidic (median net alkalinity = -13 mg/L) and had high concentrations of Fe (median Fe = 5.8 mg/L) and Al (median Al = 0.2 mg/L) (fig. 9B). The diversion wells below Rowe Drainage Tunnel (LDL) increased pH and net alkalinity and decreased dissolved Fe and Al concentrations over variable flow conditions (fig. 9B, table 1). Nevertheless, the continuous monitoring data indicated that pH was not effectively changed during extremely high flow associated with tropical-storm conditions, probably because only a small proportion of total streamflow was treated. For the period March 1999 - July 2000, the pH of stream water at site E2-0 about 300 yd (274 m) below the diversion wells usually was higher than that for site E2-244 about 300 yd (274 m) above the wells (fig. 8B). During baseflow conditions through August 1999, the diversion wells effectively increased pH in the reach below the wells by about 1-2 units. However, during stormflow conditions associated with tropical storms in September 1999 and spring storms in March-May 2000, the pH of downstream water was not effectively increased. Additional or larger diversion wells would be needed to treat these large stormflows.

Limestone-Sand Dosing (C4 - C6). The limestone-sand dosing at Coal Run (LSC) was aptly called dumping, whereby several truck loads of sand were spilled over the bank into the channel. As the mound was eroded at the base, sand spilled into the channel. Although slight increases in net alkalinity and dissolved Ca concentrations were indicated by the matched-pair tests (table 1), the pH and metals concentrations were not different between the upstream site (C4) and downstream sites (C6) for most conditions (fig. 9C, table 1).

Limestone-Compost-Based Wetlands (E3-1 - E3-2). The wetlands at Lower Rausch Creek (WLR) did not affect pH, net alkalinity, Ca, or SO4 concentrations, but were effective decreasing dissolved and total Fe and Mn and total Al over most flow conditions (fig. 9C, table 1). The implication is that the wetlands acted more as settling basins than as neutralization treatment systems. Although most pairs of samples indicated declines in metal concentrations and transport from the upstream to downstream monitoring sites, three pairs collected during different stormflow conditions indicated concentrations of total metals and suspended solids were greater at the downstream site than the upstream site.

18 Summary and Conclusions Monitoring on the mainstem of Swatara Creek indicates watershed-scale effects owing primarily to changes in mining practices and secondarily to watershed-wide implementation of treatment systems. For Swatara Creek at Ravine, Pa., long-term data collected since 1959 indicate SO4 concentration declined from about 150 mg/L in 1959 to 75 mg/L in 1999; pH increased sharply from 3.5-4.4 (median ~4) to 4.6-7.0 (median ~6) after 1975. These trends resulted from a decline in pyrite oxidation because of mine flooding and the onset of carbonate buffering attributed to natural processes (carbonate dissolution, SO4 reduction) and recently installed limestone treatment systems.

A variety of limestone treatment systems were implemented during 1995-98 for the neutralization of acidity and the removal of dissolved metals in Swatara Creek and its tributaries. Continuous measurements of pH and SC and periodic sampling for chemical constituents upstream and downstream of each treatment system show that (1) the open limestone channel and limestone- sand dosing generally had negligible effects on water quality and (2) limestone diversion wells, limestone drains, and limestone-compost based wetlands generally were effective at attenuating dissolved and suspended metals during baseflow conditions but were less effective during stormflow conditions. Generally, stormflow tends to be acidic, and, as streamflow volume increases, a smaller fraction of total flow is treated and (or) residence time in the treatment system is reduced. Furthermore, metal-rich sediments commonly can be scoured and resuspended from the streambed during stormflow conditions.

Despite the implementation of a variety of limestone treatment systems during 1996-98 for the neutralization of acidity and the removal of dissolved metals, acidification of streamflow in the Swatara Creek Basin persisted locally and periodically, particularly during extreme hydrologic conditions associated with stormflow or drought conditions. Concentration of SO4, SC, and pH for Swatara Creek were correlated and generally declined with increasing streamflow at Newtown and Ravine, indicating episodic acidification and dilution from acidic rainfall and runoff during stormflow. In contrast, SC and pH were inversely correlated at Lorberry Creek, achieving lowest and highest values, respectively, during the drought of July-August 1999 and tropical- storm period of September 1999. During the drought and subsequent high-flow conditions in 1999, the Rowe Drainage Tunnel discharge changed from near-neutral (pH 5-6; SC <400 mS/cm) to acidic (pH 3-4; SC 500-950 mS/cm) causing acidification of Lorberry Creek to its mouth on Swatara Creek. Acidification of the mine discharge probably resulted from a decline in the water table, oxidation of previously inundated pyritic material, and subsequent dissolution of oxidation products by infiltrating water or as the water table rebounded. Nevertheless, the pH of Swatara Creek at Ravine remained near-neutral during this particular event.

The recent monitoring on the tributaries and mainstem of Swatara Creek indicates the limestone diversion wells and limestone drains mitigated extreme fluctuations in pH during most storm events; however, additional buffering capacity could be needed to maintain near-neutral pH of Swatara Creek for the extreme range of conditions. Generally, to maintain stream pH during storms, additional or larger limestone diversion wells could be constructed to begin or increase alkalinity production as the stream stage rises and/or additional or larger limestone drains could be constructed to produce greater amounts of alkalinity and enhance the buffering capacity of baseflow. Increasing the buffering capacity of baseflow also will mitigate acidification effects during drought conditions. Neutralization and pH buffering alone will not remedy the problem of

19 metals transport. Alkalinity-producing systems such as limestone diversion wells or limestone drains combined with wetlands could attenuate metals transport.

Acknowledgments Kovaldas “KB” Balciauskas formerly at USGS, Jeffrey J. Chaplin, John J. Rote, Michael D. Bilger, and Robin A. Brightbill at USGS are acknowledged for critical assistance with field work and data analysis. This research was funded by the PaDEP, Schuylkill County Conservation District, and USGS under the Federal-State Cooperative Water-Resources Program.

References Cited Arnold, D. E. 1991. Diversion wells--a low-cost approach to treatment of . In Proceedings of the 12th Annual West Virginia Surface Mine Drainage Task Force Symposium: Charleston, W.Va., West Virginia Mining and Reclamation Association, p. 39-50. Baker, J.P., and C.L. Schofield. 1982. Aluminum toxicity to fish in acidic waters: Water Air Soil Pollution 18:289-309. Berg, T.M., J.H. Barnes, W.D. Sevon, and others. 1989. Physiographic provinces of Pennsylvania: Pennsylvania Geol. Surv., 4th ser. Map 13, scale 1:2,000,000. Bigham, J.M., U. Schwertmann, S.J. Traina, R.L. Winland, and M. Wolf. 1996. Schwertmannite and the chemical modeling of iron in acid sulfate waters: Geochim. Cosmochim. Acta 60:2111-2121. Bilger, M.D., R.A. Brightbill, and H.L. Campbell. 1999. Occurrence of organochlorine compounds in whole fish tissue from streams of the Lower Susquehanna River Basin, Pennsylvania and Maryland, 1992: U.S. Geol. Surv. Water-Res. Inv. Rep. 99-4065, 17 p. Brady, K.B.C., R.J. Hornberger, and G. Fleeger. 1998. Influence of geology on post-mining water quality-Northern Appalachian Basin. In Coal Mine Drainage Prediction and Pollution Prevention in Pennsylvania, eds. Brady, K.B.C., M.W. Smith, and J.H. Schueck: Harrisburg, Pennsylvania Department of Environmental Protection, 5600-BK-DEP2256, p. 8.1-8.92. Cravotta, C.A. III. 1994. Secondary iron-sulfate minerals as sources of sulfate and acidity--The geochemical evolution of acidic ground water at a reclaimed surface coal mine in Pennsylvania. In Environmental geochemistry of sulfide oxidation, eds. Alpers, C.N., and D.W. Blowes: American Chemical Society Symposium Series 550, p. 345-364. ———. 2000. Relations among sulfate, metals, sediment, and streamflow data for a stream draining a coal-mined watershed in east-central Pennsylvania. In ICARD 2000 Proceedings Fifth International Conference on Acid Rock Drainage: Littleton, Co., Society for Mining, Metallurgy, and Exploration, Inc., 1:401-410. ———, and M.D. Bilger. 2001, Water-quality trends for a stream draining the Southern Anthracite Field, Pennsylvania: Geochemistry-Exploration, Environment, Analysis 1:33-50. ———, K.B.C. Brady, A.W. Rose, and J.B. Douds. 1999. Frequency distribution of the pH of coal-mine drainage in Pennsylvania. In U.S. Geological Survey Toxic Substances Hydrology Program--Proceedings of the Technical Meeting, Charleston, South Carolina, March 8-12, 1999, eds. Morganwalp, D.W., and H. Buxton: U.S. Geol. Surv. Water-Resour. Inv. Rep. 99- 4018A, p. 313-324. ———, and M.K. Trahan. 1999. Limestone drains to increase pH and remove dissolved metals from acidic mine drainage: Appl. Geochem. 14:581-606. Durlin, R.R., and Schaffstall, W.P.1998, Water resources data, Pennsylvania, water year 1997, vol. 2, Susquehanna and Potomac River Basins: U.S. Geol. Surv. Water-Data Report PA-97- 2, p. 214-269, 342-392.

20 ———. 1999. Water resources data, Pennsylvania, water year 1998, vol. 2, Susquehanna and Potomac River Basins: U.S. Geol. Surv. Water-Data Report PA-98-2, p. 212-257, 346-400. ———. 2000. Water resources data, Pennsylvania, water year 1999, vol. 2, Susquehanna and Potomac River Basins: U.S. Geol. Surv. Water-Data Report PA-99-2, p. 212-269, 342-386. Earle, J., and T. Callaghan. 1998. Effects of mine drainage on aquatic life, water uses, and man- made structures. In Coal Mine Drainage Prediction and Pollution Prevention in Pennsylvania, eds. Brady, K.B.C., M.W. Smith, and J.H. Schueck: Harrisburg, Pennsylvania Department of Environmental Protection, 5600-BK-DEP2256, p. 4.1-4.10. Elder, J.F. 1988. Metal biogeochemistry in surface-water systems--a review of principles and concepts: U.S. Geol. Surv. Circular 1013, 43 p. Fishel, D.K. 1988. Preimpoundment hydrologic conditions in the Swatara Creek (1981-84) and estimated postimpoundment water quality in and downstream from the planned Swatara State Park reservoir, Lebanon and Schuylkill Counties, Pa: U.S. Geol. Surv. Water-Resour. Inv. Rep. 88-4087, 108 p. Fishman, M.J., and L.C. Friedman, eds. 1989. Methods for determination of inorganic substances in water and fluvial sediments: U.S. Geol. Surv. Tech. Water-Resour. Inv., Book 5, Chapter A1, 545 p. Francis, A.J., D.J. Dodge, A.W. Rose, and A.J. Ramirez. 1989. Aerobic and anaerobic microbial dissolution of toxic metals from coal wastes--mechanism of action: Env. Sci. Tech. 23:435- 441. Greenberg, A.E., L.S. Clesceri, A.D. Eaton, and M.A.H. Franson. eds. 1992. Standard methods for the examination of water and wastewater (18th): Washington, D.C., American Public Health Association. Growitz, D.J., L.A. Reed, and M.M Beard. 1985. Reconnaissance of mine drainage in the coal fields of eastern Pennsylvania: U.S. Geol. Surv. Water-Resour. Inv. Rep. 83-4274, 54 p. Hedin, R.S., R.W. Nairn, and R.L.P. Kleinmann. 1994a. Passive treatment of coal mine drainage: U.S. Bureau of Mines Information Circular IC 9389, 35 p. Hedin, R.S., G.R. Watzlaf, and R.W. Nairn. 1994b. Passive treatment of acid mine drainage with limestone: J. Environ. Qual. 23:1338-1345. Helsel, D.R., and R.M. Hirsch. 1992. Statistical methods in water resources: New York, Elsevier Publishers, Inc., Studies in Environmental Science no. 49, 522 p. Hoffman, G.L., M.J. Fishman, and J.R. Garbarino. 1996. Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory--In-bottle acid digestion of whole-water samples: U.S. Geol. Surv. Open-File Rep. 96-225, 28 p. Ladwig, K.J., P.M. Erickson, R.L.P. Kleinmann, and E.T. Posluszny. 1984. Stratification in water quality in inundated anthracite mines, eastern Pennsylvania: U.S. Bureau Mines Rep. Inv. RI 8837, 35 p. McCarren, E.F., J.W. Wark, and J.R. George. 1964. Water quality of the Swatara Creek Basin, Pa.: U.S. Geol. Surv. Open-File Rep., 88 p. Olyphant, G.A., E.R. Bayless, and D. Harper. 1991. Seasonal and weather-related controls on solute concentrations and acid drainage from a pyritic coal-refuse deposit in southwestern Indiana, U.S.A.: J. Contam. Hydrol. 7:219-236. Pennsylvania Department of Environmental Protection. 1998. Commonwealth of Pennsylvania 1998 Water Quality Assessment 305(b) Report: Pennsylvania Department of Environmental Protection, Bureau Watershed Management (www.dep.state.pa.us/dep/deputate/watermgt/wc/Subjects/WQStandards/305_wq98). P-STAT, Inc. 1989. P-STAT user’s manual: Princeton, N.J., P-STAT, Inc., vols. I and II.

21 Rantz, S.E., and others. 1982a. Measurement and computation of streamflow--1. Measurement of stage and discharge: U.S. Geol. Surv. Water-Supply Paper 2175, vol. 1, 284 p. ———. 1982b. Measurement and computation of streamflow--2. Computation of discharge: U.S. Geol. Surv. Water-Supply Paper 2175, vol. 2, 631 p. Rose, A.W. and C.A. Cravotta, III. 1998. Geochemistry of coal-mine drainage. In Coal Mine Drainage Prediction and Pollution Prevention in Pennsylvania, eds. Brady, K.B.C., M.W. Smith, and J.H. Schueck: Pennsylvania Department of Environmental Protection, 5600-BK- DEP2256, p. 1.1-1.22 Skelly & Loy, Inc. 1987. A watershed pollution study of the Swatara Creek: Pennsylvania Department of Environmental Resources Bureau of Abandoned Mine Reclamation, Final Report, n.p. Skousen, J.G., A.W. Rose, G. Geidel, J. Foreman, R. Evans, W. Hellier, and others. 1998. Handbook of technologies for avoidance and remediation of acid mine drainage: Morgantown, W.Va., National Mine Land Reclamation Center, 131 p. Smith, K.S., J.R. Ranville, G.S. Plumlee, and D.L. Macalady. 1998. Predictive double-layer modeling of metal sorption in mine-drainage systems. In Metal adsorption by geomedia, ed. E.A. Jenne: San Diego, Academic Press. Stuart, W.T., W.J. Schneider, and J.W. Crooks. 1967. Swatara Creek Basin of southeastern Pennsylvania--An evaluation of its hydrologic system: U.S. Geol. Surv. Water-Supply Paper 1829, 79 p., 3 plates. Stumm, W., and J.J. Morgan. 1996. Aquatic chemistry--chemical equilibria and rates in natural waters (3rd): New York, Wiley-Interscience, 1022 p. U.S. Geological Survey. 2000. Land cover: U.S. Geological Survey Map 66, scale 1:500,000 (http://edcwww.cr.usgs.gov/programs/lccp/nationallandcover.html). Way, J.H. 2000. Appalachian Mountain section of the Ridge and Valley province. In The Geology of Pennsylvania, ed. C.H. Schultz: Pennsylvania Geol. Surv. Special Pub. 1, p. 352-361. Webster, J.G., P.J. Swedlund, and K.S. Webster. 1998. Trace metal adsorption onto an acid mine drainage iron(III) oxy hydroxy sulfate: Env. Sci. Tech. 32:1361-1368. Wilde, F.D., D.B. Radtke, J. Gibs, and R.T. Iwatsubo. 1998. National field manual for the collection of water-quality data: U.S. Geol. Surv. Tech. Water-Resour. Inv., Book 9, Handbooks for Water-Quality Investigations, variously paged. Wood, C.R. 1996. Water quality of large discharges from mines in the anthracite region of eastern Pennsylvania: U.S. Geol. Surv. Water-Resour. Inv. Rep. 95-4243, 69 p. Wood, G.H., Jr., T.M. Kehn, and J.R. Eggleston. 1986. Deposition and structural history of the Pennsylvania Anthracite region. In Paleoenvironmental and tectonic controls in coal-forming basins of the United States, eds. P.C. Lyons and C.L. Rice: Geological Society of America Special Paper 210, p. 31-47. Wood, G.H., Jr., J.P. Trexler, and T.M. Kehn. 1968. Geologic maps of anthracite-bearing rocks in the west-central part of the southern Anthracite Field Pennsylvania, eastern area: U.S. Geol. Surv. Misc. Geologic Inv. Map I-528. Ziemkiewicz, P.F., J.G. Skousen, D.L. Brant, P.L. Sterner, and R.J. Lovett. 1997. Acid mine drainage treatment with armored limestone in open limestone channels: J. Environ. Qual. 26:1017-1024.

22