2017 Study Plan for the Lower Tuscarawas River & Select Tributaries

Tuscarawas and Coshocton Counties

Division of Surface Water Ecological Assessment Section & District Water Quality July 11, 2017

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

2017 Study Plan for the Lower Tuscarawas River & Select Tributaries

Tuscarawas and Coshocton Counties

State of Environmental Protection Agency

Division of Surface Water

Modeling, Assessment & TMDL/NPDES/Standards and Technical Assistance Sections Lazarus Government Center 122 South Front Street, Columbus, OH 43215 P.O. Box 1049, Columbus, OH 43216-1049

Ecological Assessment Section 4675 Homer Ohio Lane, Groveport, OH 43125

Southeast District Office 2195 Front Street, Logan, OH 43138

July 11, 2017

John R. Kasich, Governor Craig W. Butler, Director

State of Ohio Ohio Environmental Protection Agency

1

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Signature Page

7/12/17 Tiffani Brian Kavalec Hall or Brian Hall, Chief/Assistant Chief, Division of Surface Date Water 7/12/2017 Cathy Alexander, Manager, Modeling and TMDL Section Date 7/12/2017

Jeff DeShon, Manager, Ecological Assessment Section Date

7/12/2017 Jeff Reynolds, Quality Assurance Coordinator, Division of Surface Water Date 7/18/2017 Marco DeShaies, Water Quality Manager/Study Team Leader, Division of Date Surface Water, SEDO District

2

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Introduction As part of Ohio’s statewide monitoring strategy, a biological and water quality assessment will be conducted in the Tuscarawas River basin during the 2017 field season. This study plan will focus specifically on the lower portion of the Tuscarawas River basin while a separate effort will detail the sampling strategy for the Upper Tuscarawas River. The lower Tuscarawas River study area is entirely contained within the 8-digit hydrologic unit code (HUC) 05040001 and is composed of thirteen 12-digit watershed assessment units (WAUs) and two large river assessment units (LRAUs). The HUCs and their description are listed in Table 1. Ambient biology, macrohabitat and water quality data will be collected from the sites listed in Table 2. Details regarding matrixes, parameters, and sample frequency are also included in the Table 2. NPDES permitted facilities within the study area are listed in Tables 3 and 4. An analytical parameter list with associated methods and reporting limits is presented in Table 5.

Sampling Objectives 1 Systematically survey and assess the mainstem and principal tributaries of the lower Tuscarawas River in support of Clean Water Act Sections 305(b) and 303(d); 2 Gather ambient environmental information (biological, chemical, and physical) from designated water bodies to assess current beneficial uses (aquatic life, recreation and water supply); 3 Verify existing uses and recommend appropriate uses for unverified designations; 4 Document any changes in biological, chemical, and physical conditions of the study areas where historical information exists, thus expanding the Ohio EPA data base for statewide trends analysis (e.g., 305[b]); 5 Collect fish tissue samples at selected stations in support of the Ohio Sport Fish Health and Consumption Advisory program; 6 Evaluate major and significant minor NPDES permitted facilities within the study area; and 7 Support mine drainage impact assessment, AMDAT development and audit the effectiveness of mine drainage abatement activities sponsored or implemented by Ohio DNR, DMRM.

Total Maximum Daily Load (TMDL) The Tuscarawas River basin has had one TMDL approved and this survey represents the second sampling effort in support of a TMDL. TMDLs were calculated for sediment, nutrients (total phosphorus), habitat factors, and pathogens and was approved by US EPA on September 15, 2009. A copy of the report can found on the Ohio EPA web page at: http://www.epa.state.oh.us/portals/35/tmdl/TuscarawasRiverTMDL_final_jul09_wo_app.pdf. In addition to the items listed above, information collected as part of this survey will support the refinement or modification of the existing TMDL and nutrient mass balance modeling at or greater than the HUC 10 scale.

Beneficial Use Assessment The status of all beneficial uses within the study area will be evaluated. Current beneficial use designation for waters that together form the lower Tuscarawas basin are listed in Ohio Administrative Code (OAC) Chapter 3745-1-24, Table 24-1.

Including the mainstem, 31 lower Tuscarawas waterbodies are listed the Ohio WQS. Regarding aquatic life uses, all or portions of these are designated warmwater habitat (WWH). Only the lower 47 miles of the mainstem, from Stillwater Creek down to the mouth, are designated EWH. Of these, the aquatic life use designation of 18 streams or stream segments are field verified. The remaining 13 waterbodies were designated prior to the adoption of biocriteria and thus represent desktop designations, per the 1978 WQS.

3

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

All of the streams within the study area are presently listed as primary contact for the recreation use. As with the aquatic life use, 18 of the 31 streams or stream segments are field verified, and the remaining 13 are unverified. The agricultural water supply (AWS) and industrial water supply (IWS) uses apply to all of the streams within the lower Tuscarawas basin.

The 2017 field work will also support the Ohio Sport Fish Health and Consumption Advisory programs. Presently, all waters have a statewide/nationwide advisory in place of 1 meal per week for sensitive populations (young children and women of childbearing age) due to mercury.

SAMPLING ACTIVITIES

Chemical/Physical Water and Sediment Chemical sampling locations within the study area are listed in Table 2. Conventional chemical/physical water quality samples will be collected 5 times at each designated location. Analysis at these stations will include organics. Sediment samples will be collected at 8 locations. Chemical parameters to be tested are listed in Table 5. Surface water sampling will occur across a variety of flow conditions, from lower flows to moderate and higher flows.

Bacteriological Sampling Water samples for bacteriological analyses will be collected to determine PCR use status. Each WAU will have at least 1 site sampled, provided access is secured. Most of this effort will focus on waterbodies with public access and documented recreational use. The Tuscarawas River mainstem is a larger river and is subject to regular use by the general public, this facilited by multiple canoe liveries and numerous points of public access. Ideally, 10 E. coli samples should be collected per station from the mainstem if resources allow. Each tributary bacteria site will have a minimum of 5 sets of E. coli samples collected within a 90- day period, between the Memorial Day and Labor Day holidays.

Stream Nutrient Assessment Specialized sampling will be done at locations designated as nutrient sites in Table 2. The objectives of these sites are to determine the influence of nutrients on the aquatic life beneficial use and further support nutrient criteria development. There are two designations for nutrient assessment, one for large rivers and another for wadeable streams. The sampling requirements for these protocols were detailed in a memo to Ohio EPA water quality staff (Dripps and Alexander, personal communication, February 22, 2017). A combination of benthic and sestonic chlorophyll-a samples will be collected. Benthic chlorophyll samples will be collected at least once and should be timed to coincide with deployment of multiparameter water quality sondes during stable, baseflow conditions. The sondes capture the diel dynamics of dissolved oxygen in the stream, reflecting the balance of primary production and respiration in the stream.

Biological Assemblage and Habitat Fish and macroinvertebrate assemblage samples will be collected as listed in Table 2. Qualitative Habitat Evaluation Index (QHEI) scores will be calculated by evaluating habitat at all fish sampling locations.

Fish Tissue Ten stations on the lower Tuscarawas River mainstem are scheduled for fish tissue collection. Based on the availability of target species, access, or other factors that may emerge over the course of the summer sampling season, field staff responsible for this work may expand, contact, or otherwise modify the fish

4

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017 tissue work plan as field conditions necessitate. Fillet samples of edible size sport or target species will be tested for organochlorine pesticides, PCBs, mercury, lead, cadmium, arsenic, and selenium. Results will be used in in support of the Ohio Sport Fish Health and Consumption Advisory Programs. Final sampling locations may vary from those identified in Table 2, and parameters to be analyzed are listed in Table 5.

Property Access Of the 54 sites scheduled for evaluation as part of the 2017 lower Tuscarawas River survey, full access has yet to be secured for nine locations. Sites so identified are highlighted in Table 2. Specifics regarding these sites (type of access still needed, landowner and contact information, summary of efforts thus far, potential alternatives, etc.) are detailed in the master access table dated June 27, circulated to all study team members via e-mail and uploaded to the agency’s lower Tuscarawas share point site. Outstanding sites include three mainstem locations and six tributaries. Of these, seven are identified as being high priority, in that every effort should be made by field staff to secure permission while actively working in the study area. To that end, all field staff are advised to have permission forms and prepaid envelopes at hand while in the study area so they may attempt to obtain written consent for these sites when and if time permits. The remaining three low priority sites may be dropped from the study area entirely or converted to chem. only, per input from study team as field work unfolds. This is not to say that no effort should be made to gain lawful entry at these locations, thus allowing for a full assessment; rather, it is in recognition that in many instances, three or more attempts have already been made at these and other sites, thus we may be approaching the point of diminishing returns. If access is secured following the distribution of the final study plan, please forward a scanned copy of the consent form to Charles Boucher (acting study plan coordinator) so the master access document may be updated. For other changes or developments related to the outstanding sites, promptly inform Charles Boucher as well.

QUALITY ASSURANCE/SAMPLING METHODS Ohio EPA Manuals All biological, chemical, data processing, and data analysis methods and procedures adhere to those specified in the Surface Water Field Sampling Manual for water quality parameters and flows (Ohio EPA 2015c), Biological Criteria for the Protection of Aquatic Life, Volumes II - III (Ohio EPA 1987, 1989a, 1989b, 2015a, 2015b), The Qualitative Habitat Evaluation Index (QHEI); Rationale, Methods, and Application (Ohio EPA 1989b, 2006) for habitat assessment, Ohio EPA Sediment Sampling Guide and Methodologies (Ohio EPA 2015e) and Ohio EPA Fish Tissue Collection Guidance Manual (Ohio EPA 2012b).

Use Attainment Attainment/non-attainment of aquatic life uses will be determined by using biological criteria codified in OAC 3745-1-07, Table 7-1. Numerical biological criteria are based on multi-metric biological indices including the Index of Biotic Integrity (IBI) and modified Index of Well-Being (MIwb), indices measuring the response of the fish community, and the Invertebrate Community Index (ICI), which indicates the response of the macroinvertebrate community. Performance expectations for the basic aquatic life uses (Warmwater Habitat [WWH], Exceptional Warmwater Habitat [EWH], and Modified Warmwater Habitat [MWH]) were developed using the regional reference site approach (Hughes et al. 1986, Omernik 1987). This fits the practical definition of biological integrity as the biological performance of the natural habitats within a region (Karr and Dudley 1981). Attainment of an aquatic life use is FULL if all three indices (or those available) meet the applicable criteria, PARTIAL if at least one of the indices did not attain and performance did not fall below the fair category, and NON if all indices either fail to attain or any index indicates poor or very poor performance. The results will be compared to the appropriate tiered WWH biocriteria based on ecoregion. 5 Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Attainment/non-attainment of recreational uses will be determined using E. coli criteria codified in OAC 3745-1-37, Table 37-2. Water quality must meet a 90-day geometric mean and a statistical threshold not to be exceeded more than 10% of the time.

Stream Habitat Evaluation Physical habitat is evaluated using the QHEI developed by the Ohio EPA for streams and rivers in Ohio (Ohio EPA 1989b, 2006). Various attributes of the available habitat are scored based on their overall importance to the establishment of viable, diverse aquatic faunas. Evaluations of type and quality of substrate, amount of instream cover, channel morphology, extent of riparian canopy, pool and riffle development and quality and gradient are among the metrics used to evaluate the characteristics of a stream segment, not just the characteristics of a single sampling site. As such, individual sites may have much poorer physical habitat due to a localized disturbance yet still support aquatic communities closely resembling those sampled at adjacent sites with better habitat, provided water quality conditions are similar. QHEI scores from hundreds of segments around the state have indicated that values higher than 60 were generally conducive to the establishment of warmwater faunas while those which scored in excess of 75-80 often typify habitat conditions which have the ability to support exceptional faunas.

Biological Community Assessment A combination of quantitative and qualitative methods will be employed to monitor benthic macroinvertebrate communities. Quantitative collections are made through the use of modified Hester- Dendy multiple-plate artificial substrate samplers, deployed at all biomonitoring sites draining more than 20 mi2, or at reference sites regardless of size. Once deployed, artificial substrates are left to colonize, in- stream, for a six-week period. Qualitative sampling will be conducted at all biomonitoring stations. This sampling method consists of a basic inventory of macroinvertebrate taxa from natural substrates, noting dominant taxa among major habitat types (e.g., riffle, run, pool, and margin. Macroinvertebrate sampling protocols are detailed in Ohio EPA (2015b).

Fish will be sampled at each sampling location using pulsed DC headwater, wading, or boat electrofishing methods depending on watershed size at each sampling zone. Sites with drainage areas greater than 20 mi2 or at reference site locations will be sampled twice during the sampling index period. The number of passes may be adjusted as necessary based on best professional judgment of the fish crew leader. Reasons for a single pass monitoring at sites otherwise identified as needing multiple passes may include, extremely difficult and time consuming access, work delays related to weather, or the emergence of alterations (natural or otherwise) to the point of access or sampling reach, rendering replication of the initial effort hazardous or costly, or both. Fish are processed in the field which includes identifying each individual to species, counting individuals at all sites, weighing individuals at wading and boat sites, and recording external abnormalities. Fish sampling protocols are detailed in Ohio EPA (2015b).

Surface Water Surface water grab samples will be collected and preserved using appropriate methods as outlined in the Surface Water Field Sampling Manual for water quality parameters and flow (Ohio EPA 2015c) and delivered to the Ohio EPA Division of Environmental Services (DES) for analyses. Field measurements of dissolved oxygen, pH, temperature and conductivity will be made using YSI Professional Plus meters along with all grab samples for surface water chemistry. Datasonde® continuous recorders will be placed at select locations to evaluate diel measurements of dissolved oxygen, pH, temperature, and conductivity.

6

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Bacteria Water samples will be collected into appropriate containers, cooled to 4°C, and transported to DES lab in Reynoldsburg, Ohio within 6 hours of sample collection. All samples will be analyzed for E. coli bacteria using U.S.EPA approved methods.

Chlorophyll Benthic and sestonic chlorophyll a samples will be collected and preserved using appropriate methods, as outlined in Appendix II of the Surface Water Field Sampling Manual (Ohio EPA 2015d) and delivered to the Ohio EPA DES for analyses. Alkalinity must be requested as a routine water quality parameter at all study sites along with the routine field parameters, especially temperature and pH.

Sediment Fine grained multi-incremental sediment samples will be collected in the upper 4 inches of bottom material using either decontaminated stainless steel scoops or dredges. Collected sediment will be placed into appropriate containers, placed on ice (to maintain 4oC) and shipped to the Ohio EPA DES for analysis. Sampling and decontamination protocols will follow those listed in Appendix III of the Surface Water Field Sampling Manual (Ohio EPA 2015e).

Fish Tissue Tissue fillet samples will be collected from fish of edible size. Preferred species for analysis include spotted bass, largemouth bass, smallmouth bass, flathead catfish, walleye, saugeye, white bass, common carp, freshwater drum, buffalo, and channel catfish. Fillet composite samples should be comprised of three or more similarly sized fish of the same species that together yield at least 150 grams of tissue. At each monitoring location, the collection of tissue samples from five species is desirable and may include multiple samples (size classes and sample type or both) for each taxa.

Results from previous investigations on the lower Tuscarawas River found high levels of PCB contamination in channel catfish, flathead catfish and common carp. Therefore, every effort must be made to collect multiple samples (whole-body, composite, or a combination thereof) of each of these special target species. Three samples of each species are needed to allow for delisting/listing determinations per the fish consumption advisory program. This is not to suggest that field investigators should forgo the collection of other preferred species, rather, it is intended to convey that priority should be given to these special target species should sampling, processing or storage constraints arise.

As fish tissue investigations are limited to the lower Tuscarawas mainstem, samples will be collected with standard boat electrofishing methods (Ohio EPA 1987b). Scheduled sampling locations are listed in Table 2 and parameters to be analyzed are listed in Table 5. Fish captured for tissue analysis will be filleted in the field using decontaminated or otherwise clean materials (stainless steel fillet knives, cutting boards, and nitrile, PCV, or powder-free latex glove). Filleted samples will be wrapped in aluminum foil, placed in a sealed plastic bag with data label, and placed on dry ice if samples are to be stored in the field for a period greater than 48 hours or wet ice, if field storage will be less than 48 hours. Sampling, decontamination and storage protocols will follow those listed in the Ohio EPA Fish Collection Guidance Manual (2012b). Fish tissue samples will be stored in chest freezers at the Ohio EPA Groveport Field Facility prior to delivery to DES.

7

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Field Quality Control Samples Ten percent of the total water samples will be submitted to the lab as field quality control samples. About 5% will be duplicates, including replicates if necessary, and about 5% will be blanks, including equipment blanks. Matrix spike duplicates will be collected for organic water samples at a minimum of 5%. Field instruments will be calibrated daily using manufacturer guidelines. One sonde recorder site will have two instruments placed in the river as a duplicate. All field quality control requirements and data validation methods are detailed in the Surface Water Field Sampling Manual (Ohio EPA 2015c).

8

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Table 1. Description of Lower Tuscarawas River & Selected Tributaries watershed assessment units (12- digit HUCs)

Watershed Description Assessment Unit 05040001-12-04 Wolf Run-Tuscarawas River 05040001-17-01 Stone Creek 05040001-17-02 Oldtown Creek 05040001-17-03 Beaverdam Creek 05040001-17-04 Pone Run-Tuscarawas River 05040001-18-01 Dunlap Creek 05040001-18-02 Mud Run-Tuscarawas River 05040001-18-03 Buckhorn Creek 05040001-18-04 Blue Ridge Run-Tuscarawas River 05040001-19-01 Evans Creek 05040001-19-02 West Fork White Eyes Creek 05040001-19-03 White Eyes Creek 05040001-19-04 Morgan Run-Tuscarawas River 05040001-90-02 Tuscarawas River ( to Stillwater Creek) 05040001-90-03 Tuscarawas River (Stillwater Creek to )

9

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Table 2. Lower Tuscarawas River study area sampling stations.

Area STATION LANDMARK RM LRAU-HUC 12 USGS 7.5’ Quad LATITUDE LONGITUDE SAMPLING (miles2) Tuscarawas River (17-500-000)

R06K03 NEAR BOLIVAR, DST. SANDY CREEK 72.6 1090 05040001 90 02 Bolivar 40.655000 -81.445300 NLR, BLR, C, F2, MQ, FT ZOAR DAM - RM 68.8 (Breached) R06K01 @ TWP. RD. 387 DST. ZOAR DAM 68.72 1103 05040001 90 02 Dover 40.609200 -81.428600 BLR,C, F2, MQ DOVER DAM - RM 63.8 R10K18 NE OF DOVER @ POWER LINES DST. DOVER DAM 63.2 1404 05040001 90 02 Dover 40.551900 -81.420800 NLR, F2, MQ, FT

R10W17 UPST. DOVER @ DST. END OF ISLAND 62.3 1405 05040001 90 02 Dover 40.540300 -81.429400 NLR, BLR, C R10K15 JUST DST. EMPIRE DETROIT STEEL 59.2 1413 05040001 90 02 Dover 40.521700 -81.466400 ES,

R10S22 AT DOVER @ ST. RT. 800 58.8 1413 05040001 90 02 Dover 40.520000 -81.473300 ES, NLR, BLR, C, FT DOVER-NEW PHILIDELPHIA LOWHEAD DAM - RM 58.3

R10S03 AT NEW PHILADELPHIA @ BROADWAY AVE 54.9 1814 05040001 90 02 New Philadelphia 40.481100 -81.448100 NLR, BLR, C, FT NEW PHILADELPHIA LOWHEAD DAM - ~RM 54.55

R10G01 UPST. NEW PHILADELPHIA WWTP @ GAGE 54.2 1814 05040001 90 02 New Philadelphia 40.483300 -81.438900 NLR, F2, MQ 1.0 MI dst. US 250/DST. NEW PHILIDELPHIA R10K12 52.3 1816 05040001 90 02 Urichsville 40.4642 -81.4308 NLR, BLR, C, F2, MQ, FT WWTP STATE DAM?? (per local tackle shop, Impoundment, relic dam, or other navigation hazard between RMs 52.3 and 47.4)

R10S17 UPST. STILLWATER CREEK @ MOOSE ISLAND HP 47.4 1870 05040001 90 02 New Philadelphia 40.429400 -81.374400 NLR, BLR, C, F2, MQ

R10K10 AT TUSCARAWAS @ CR 60 44.5 2364 05040001 90 03 New Philadelphia 40.394400 -81.390300 NLR, BLR, C, F2, MQ, FT

R10G02 UPST. GNADENHUTTEN WWTP 38.7 2381 05040001 90 03 Gnadenhutten 40.359000 -81.442000 NLR, BLR, C, F2, MQ

R10S13 AT PORT WASHINGTON @ CR 14 30.9 2399 05040001 90 03 Newcomerstown 40.326100 -81.511700 NLR, BLR, C, F2, MQ, FT R 611790 AT NEWCOMERSTOWN @ RIVER ST. 21.2 2443 05040001 90 03 Newcomerstown 40.261100 -81.609700 NLR, BLR, C, Sd, O, F2, MQ

DST. NEWCOMERSTOWN, 0.5 MI. DST BLUE RIDGE NLR, BLR, C, Sd, O, F2, R10G04R 17.6 2474 05040001 90 03 Fresno 40.276700 -81.650600 RUN MQ, FT R 601840 AT WEST LAFAYETTE @ SR 93 HP 10.7 2566 05040001 90 03 Fresno 40.291900 -81.750000 NLR, BLR, C, Sd, O, F2, MQ R R10K04 W OF WEST LAFAYETTE, ADJ. US 36 HP 6.9 2574 05040001 90 03 Coshocton 40.290800 -81.801400 NLR, BLR C, Sd, O, FT

10

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Area STATION LANDMARK RM LRAU-HUC 12 USGS 7.5’ Quad LATITUDE LONGITUDE SAMPLING (miles2) Tuscarawas River (17-500-000)

R10S34 JUST E OF CANAL LEWISVILLE, BETWEEN TRIBS. 3.8 2585 05040001 90 03 Coshocton 40.298300 -81.830000 NLR, BLR, C, F2, MQ

611730 AT COSHOCTON @ KIA BRIDGE 0.3 2596 05040001 90 03 Coshocton 40.278900 -81.870600 NLR, BLR, C, F2, MQ, FT Wolf Run (17-529-000)P1 R06G72 W OF ZOAR @ LANE OFF TR 411 0.25 3.49 05040001 12 04 Dover 40.619800 -81.445400 C, F, Mq Middle Run (17-528-000) P1 R06G71 JUST S OF ZOAR @ TR 411 0.23 2.1 05040001 12 04 Dover 40.620000 -81.445000 C, F, Mq Trib to Tuscarawas River (70.31) (17-500-030) 303823 At Towpath Rd 0.23 3.63 05040001 12 04 Dover 40.6242 -81.4527 B, C, F, Mq Crooked Run (17-526-000) R10G42 Kneuss Rd. Note: Moved site ust. for access 2.42 - 05040001 17 01 Stone Creek 40.48732 -81.52420 C, F, Mq R10K32 @ STONE CREEK RD. 0.28 7.6 05040001 17 01 New Philadelphia 40.483100 -81.485200 C, F, Mq Stone Creek (17-525-000)P1 R10G39 N OF STONE CREEK @ RICE RD. (TR 243) 6.37 9.1 05040001 17 01 Stone Creek 40.421000 -81.545000 N, C, F, Mq R10G40 NE OF STONE CREEK @ MATHIAS RACEWAY RD. 4.27 23 05040001 17 01 Stone Creek 40.443000 -81.518000 S, C, F2, Mq R10G41 SR 250/I77 off-ramp or Steele Hill Rd. SEE FOOTNOTE 0.09/0.9 38.9 05040001 17 01 New Philadelphia 40.478000 -81.483000 N, B, Sd, C, O, F2, MQ UN Stone Creek Trib. (RM 7.27) (17-525-003)P1 303824 Stone Creek Rd 0.08 3.5 05040001 17 01 Stone Creek 40.4144 -81.5549 C only Oldtown Creek (17-524-000) P4 R10G36 @ TR263 (AUBIHL RD.) 5.39 8.9 05040001 17 02 New Philadelphia 40.423000 -81.48800 C, F, Mq R10K26 Goshen Hill Rd. 0.9 19.3 05040001 17 02 New Philadelphia 40.4606 -81.4394 N, B, C, F, Mq Beaverdam Creek (17-523-000)P3 R10G34 E OF NEW PHILADELPHIA @ CR 313 6.9 9.5 05040001 17 03 New Philadelphia 40.494000 -81.395000 C, F, Mq R10P06 AT SCHOENBRUNN @ BALTZLEY VALLEY RD. 0.94 19.5 05040001 17 03 New Philadelphia 40.459700 -81.401400 N, B, C, F, Mq Pone Run (17-522-000)P1 R10G33 AT MIDVALE @ SOUTH ST. 0.1 4 05040001 17 04 Uhrichsville 40.437200 -81.372900 C, F, Mq Pike Run (17-521-000) P4 301825 State St., ust. Pone Run 0.6 6.2 05040001 17 04 New Philadelphia 40.434000 -81.378000 B, C, F, Mq

11

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Area STATION LANDMARK RM LRAU-HUC 12 USGS 7.5’ Quad LATITUDE LONGITUDE SAMPLING (miles2) Mud Run (17-520-000)P3 MUD RUN NEAR TUSCARAWAS @ MORAVIAN R10G31 1.61 6.5 05040001 18 02 New Philadelphia 40.407275 -81.410358 B, C, F, Mq CHURCH RD. Frys Creek (17-518-000) P3 R10G29 FRYS VALLEY RD. 0.1 5.7 05040001 18 02 Gnadenhutten 40.374000 -81.468000 C, F, Mq Dunlap Creek (17-515-000)P4 R10G25 Dog Rd./TR 118 Note: consolidated all at RM 2.05 HP 2.05 18.0 05040001 18 01 Newcomerstown 40.272 -81.543 N, B, C, F, Mq Browning Run (17-516-000)P4 R10G27 E OF NEWCOMERSTOWN @ SR 258 0.3 7.2 05040001 18 01 Newcomerstown 40.271100 -81.527500 C, F, Mq Buckhorn Creek (17-511-000)P4 R10G19 NE OF NEWCOMERSTOWN @ TR 215 4.88 10.3 05040001 18 03 Newcomerstown 40.318000 -81.572000 S, C, F, Mq R10G20 AT NEWCOMERSTOWN @ CROSS ST. HP 1.43 22 05040001 18 03 Newcomerstown 40.279 -81.608 N, B, C, F2, MQ Lick Run (17-512-000) P4 N OF NEWCOMERSTOWN @ FIRST LANE NEAR R10G21 0.4 2.5 05040001 18 03 Newcomerstown 40.290000 -81.605000 C, F, Mq MOUTH Blue Ridge Run (17-509-000) R10G17 NW OF NEWCOMERSTOWN @ TR 105 0.4 2.9 05040001 18 04 Fresno 40.282000 -81.640000 B, C, F, Mq Evans Creek (17-505-000) R10G11 S OF BAKERSFIELD @ TR 246 HP 5.09 12.2 05040001 19 01 Fresno 40.335000 -81.644000 S, C, F, Mq R10G13 AT ORANGE @ SR 751 0.67 23.5 05040001 19 01 Fresno 40.290000 -81.683000 N, B, Sd, C, O, F2, MQ White Eyes Creek (17-502-000)P1 R10G07 N OF FRESNO @ 2ND BRIDGE OFF SR 93 5.07 11.3 05040001 19 03 Fresno 40.369000 -81.728000 N, C, F, Mq R10G09R Adj. SR 93, ust. US 36 0.2 52 05040001 19 03 Coshocton 40.305000 -81.750000 N, B, Sd, C, O, F2, MQ West Fork White Eyes Creek (17-504-000) P1 R10P11 AT CHILI @ CR 10 3.68 11.7 05040001 19 02 New Bedford 40.375800 -81.753300 S, C, F, Mq R10P10 AT FRESNO @ MOUTH 0.01 20.9 05040001 19 02 Fresno 40.331400 -81.743100 N, B, C, F2, MQ UN W. Fk. White Eyes Cr. Trib. at RM 4.1 (17-504-001) P1 303827 TR 231 0.93 5.6 05040001 19 02 New Bedford 40.3959 -81.7544 C, F, Mq

12

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Area STATION LANDMARK RM LRAU-HUC 12 USGS 7.5’ Quad LATITUDE LONGITUDE SAMPLING (miles2) East Fork White Eyes Creek (17-503-000) P1 R10P09 W OF BAKERSVILLE @ CR 2 5.93 5.6 05040001 19 03 Fresno 40.352500 -81.682500 C, F, Mq R10P08 TR 172 Note: site moved dst. due to access issues 0.1 12.6 05040001 19 03 Fresno 40.327951 -81.737107 N, C, F, Mq Morgan Run (17-501-000) P4 R10G06 E OF COSHOCTON @ CR 16 0.6 3.9 05040001 19 04 Coshocton 40.276000 -81.807000 C, F, Mq R - Ecoregional Reference Station - Stations where access has not yet been secured are highlighted in yellow. Ohio DNR MRM, AMDAT Prioritization HP - Stations so identified are of the highest high priority RE obtaining access over the course of the study. Non P1 - AMDAT a high priority, greatest recovery potential priority sites may be dropped or converted to chem. only, per input from study team as field work unfolds. P2 - AMDAT under consideration, moderate recovery - Changes in station or meta data, or both compared against study plan version 4.2 are highlighted in blue. potential P3 - AMDAT likely not feasible, low recovery potential * - Lower Stone Creek alternative: Steele Hill Rd. RM 0.9, 40.478196/-81.482985. As of 7/11/17, bridge access only. P4 - AMDAT not warranted due to multiple factors.

Sample Type Key # Sites # Passes Stream Water Chemistry Totals C 49 245  subset with larger river nutrient protocol NLR 17 -  subset with nutrient protocol N 11 -  subset with organics O 8 -  subset with sonde S 4 - Macroinvertebrate Totals 48 30  quantitative MQ 20 20  qualitative Mq 26 9 HD eqvl. Fish Totals 48 66  fish 2 pass F2 20 40  fish 1 pass F 26 26 Fish Tissue (potential) FT 10 E. coli tributaries B 12 60 E. coli large river BLR 16 160 Sediment Sd 8 Enhanced Sampling ES 2

13

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Table 3. Facilities with Individual NPDES permits located in the study area (majors highlighted in red).

Permit Facility Name No. Receiving Stream RM Description Wilkshire Hills WWTP 0PJ00008 UT to Tuscarawas River (RM 70.01) 0.4 Activated Sludge (0.75 MGD) Greer Steel Co 0ID00003 Tuscarawas River 59.64 Steel Mill Dover Municipal Light & Power Plant 0IB00016 Tuscarawas River* 59.1 Coal-Fired Power Plant Dover WWTP 0PD00005 Tuscarawas River 58.28 Oxidation Ditch BADCT (3.0 MGD) Arizona Chemical Co LLC 0IF00041 Tuscarawas River* 57.82 Organic Chemical Plant (1.48 MGD) Stone Creek WWTP 0PA00104 Stone Creek 8.1 Extended Aeration Plant (30,000 GPD) UTs to Crooked Creek (RM Kimble Sanitary Landfill 0IN00159 2.29/1.91) 0.64/0.29 Sanitary Landfill Storm Water Control New Philadelphia WTP 0IX00050 Tuscarawas River 55.19 Water Treatment - Ion Exchange System Activated Sludge- Extended Aeration (4.5 New Philadelphia WWTP 0PD00012 Tuscarawas River 53.63 MGD) Penn-Ohio Coal Co - Young Property 0IP00849 UT to Oldtown Creek (RM 1.81) 1.25 Acid Mine Drainage Treatment Wood's Tall Timber Lake Camp 0PR00155 UT to Beaverdam Creek (RM 7.78) 0.9 Seasonal Package Plant (17,000 GPD) Kerr Allotment WWTP 0PG00031 UT to Beaverdam Creek (RM 0.92) 0.1 Package Plant (8,000 GPD) Tuscarawas WWTP 0PB00083 Tuscarawas River 44.72 Lagoon Plant (0.17 MGD) Gnadenhutten WWTP 0PB00017 Tuscarawas River 38.56 Activated Sludge (0.25 MGD) Boltaron, Inc 0IQ00019 UT to Tuscarawas River (RM 22.84) 0.22 NCCW-Storm Water Newcomerstown WWTP 0PD00024 Tuscarawas River 19.8 Trickling Filter (1.25 MGD) Activated Sludge-Extended Aeration West Lafayette WWTP 0PB00048 Tuscarawas River 11.2 (0.32 MGD) Ohio Central Railroad - Morgan Run Shop 0IT00006 Tuscarawas River 6.02 Oil/Water Separator - Storm Water WestRock CP LLC 0IA00005 Tuscarawas River* 1.04 Pulp and Paper Mill - Closed * primary discharge outfall location

14

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Table 4: Facilities with General Surface Coal Mine NPDES permit location in study area

Permit Number HUC 12 Permittee and Mine Name Stream(s) and Discharge RM(s) 0504000 Penn-Ohio Coal Co. – Yorktown 0GM000364 11 701 UNT to Stone Creek (RM 2.05) Penn-Ohio Coal Co. – Shell 0GM00442 11 701 UNT to Stone Creek (RM 3.0) Harmon Hill IBR 0GM00478 11 701 UNT to Stone Creek (RM 4.1) Kimble Clay and Limestone – 0GM00470 11 701 Harmon Hill UNT to Stone Creek (RM 4.5) Harmon Hill 2 Mine Site 0GM00470 11 701 UNT to Stone Creek (RM 4.5) Penn-Ohio Coal Co. – Young 0IP00849*AD 11 702 UNT to Tuscarawas River (RM 55.6) Property 11 704 & UNT to Oldtown Creek (RM 1.81) Penn-Ohio Coal Co. – Young 0GM00469*BG 11 702 UNT to Tuscarawas River (RM 55.6) Property 11 704 & UNT to Oldtown Creek (RM 1.81) Oxford Mining Co. – Schupp 0GM00529 11 802 UNT to Oldtown Creek (RM 3.0) Oxford Mining Co. – Garrett 0GM00500 11 802 UNT to Tuscarawas River (RM 36.8) Quigley Surface Mine 0GM00491 11 801 UNTs to Dunlap Creek (RM 4.1) Oxford Mining Co. – Hunt 0GM00548 11 801 UNT to Browning Run (RM 1.80) Penn-Ohio Coal Co. – Port 0GM00435*AG 11 803 Washington UNT to Buckhorn Creek (RM 6.63)

15

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Figure 1: Schematic of the Lower Tuscarawas main stem showing sampling sites and significant wastewater dischargers.

16

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Table 5. List of parameters that will be analyzed by matrix with analytical method and reporting limits presented. Temperature, dissolved oxygen, pH and conductivity will be measured in the field.

Parameter Method Stream Lake Sediment Fish Tissue Oxygen Demand BOD, 5-dayA SM 5210B 2 mg/L cBOD, 20-dayB OEPA 310.2 2 mg/L Chemical Oxygen Demand SM 5220D 20 mg/L Physical Properties Alkalinity USEPA 310.1 5 mg/L 5 mg/L Acidity SM 2320B 5 mg/L 5 mg/L Hardness USEPA 200.7 10 mg/L 10 mg/L Specific Conductance SM 2510B 1 µS/cm 1 µS/cm Total Dissolved Solids SM 2540C 10 mg/L 10 mg/L Total Suspended Solids SM 2540D 5 mg/L 5 mg/L Turbidity OEPA 180.1 2 NTU 2 NTU % Solids SM 2540G 0 % % Lipids OEPA 581.5 0 % Nutrients Ammonia USEPA 350.1 0.05 mg/L 0.05 mg/L Bromide USEPA 300.1 20 µg/L Nitrate-Nitrite USEPA 350.1 0.5 mg/L 0.5 mg/L Nitrite USEPA 353.2 0.02 mg/L 0.02 mg/L Total Kjeldahl Nitrogen USEPA 351.2 0.2 mg/L 0.2 mg/L Total Phosphorus USEPA 365.4 0.01 mg/L 0.01 mg/L 50 mg/kg OrthophosphateB USEPA 365.4 0.01 mg/L 0.01 mg/L Total Organic CarbonA SM 5310B 2 mg/L 2 mg/L 0.1 % Anions Carbonate/Bicarbonate SM 2320B 5 mg/L 5 mg/L Chloride USEPA 325.1 5 mg/L 5 mg/L Sulfate USEPA 375.2 10 mg/L 10 mg/L Cations Aluminum USEPA 200.7 200 µg/L 200 µg/L Barium USEPA 200.7 15 µg/L 15 µg/L Calcium USEPA 200.7 2 mg/L 2 mg/L Iron USEPA 200.7 50 µg/L 50 µg/L Magnesium USEPA 200.7 1 mg/L 1 mg/L Manganese USEPA 200.7 10 µg/L 10 µg/L Potassium USEPA 200.7 2 mg/L 2 mg/L Sodium USEPA 200.7 5 mg/L 5 mg/L Strontium USEPA 200.7 30 µg/L 30 µg/L Toxic Metals Zinc USEPA 200.7 10 µg/L 10 µg/L 8 mg/kg 0.05 Arsenic USEPA 200.8/SM 3113B 2 µg/L 2 µg/L 0.8 mg/kg mg/kg .004 Cadmium USEPA 200.8/SM 3113B 0.2 µg/L 0.2 µg/L 0.08 mg/kg mg/kg Chromium USEPA 200.8 2 µg/L 2 µg/L 0.8 mg/kg Copper USEPA 200.8 2 µg/L 2 µg/L 0.8 mg/kg

17

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Parameter Method Stream Lake Sediment Fish Tissue 0.04 Lead USEPA 200.8/SM 3113B 2 µg/L 2 µg/L 0.8 mg/kg mg/kg Nickel USEPA 200.8 2 µg/L 2 µg/L 0.8 mg/kg 0.05 Selenium USEPA 200.8/SM 3113B 2 µg/L 2 µg/L 0.8 mg/kg mg/kg 0.02 Mercury USEPA 245.1/SM 3113B 0.02 mg/kg mg/kg

Bacteria Escherichia coliform USEPA 1603 10 CFU 10 CFU Algal Biomass/Toxins Chlorophyll-aC USEPA 445.0 0.3 µg/L 0.3 µg/L Microcystins OEPA 701.0 0.3 µg/L Saxitoxin OEPA 702.0 0.022 µg/L Cylindrospermopsin OEPA 703.0 0.05 µg/L Organic Compounds Atrazine OEPA 704.0 0.2 µg/L 0.2 µg/L Glyphosate USEPA 547 5 µg/L Carbomates USEPA 531.1 0.5 µg/L Semi-volatile organics USEPA 625 2-20 µg/L Semi-volatile organics USEPA 8270C .4-2 mg/kg Aldrin USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg a-BHC USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg b-BHC USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg d-BHC USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg y-BHC USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg 4,4’-DDD USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg 4,4’-DDE USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg 4,4’-DDT USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Dieldrin USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Endosulfan I USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Endosulfan II USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Endosulfan Sulfate USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Endrin USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Endrin Aldehyde USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Heptachlor USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Heptachlor Epoxide USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Methoxychlor USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Mirex USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg Hexachlorobenzene USEPA 8082A/OEPA 590.1 4 µg/kg 10 µg/kg PCB-1016 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1221 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1232 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1242 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1248 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1254 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg PCB-1260 USEPA 8082A/OEPA 590.1 20 µg/kg 50 µg/kg A – Sample parameter for streams if greater than 500mi2 B - Sample parameter for streams greater than 500mi2 and all nutrient sites (MAS will sample CBOD20) C- Sample sestonic Chl-a for streams greater than 700mi2 and nutrient sites (MAS will sample nutrient sites < 700mi2) and sample benthic Chl-a for nutrient sites only (MAS will sample)

18

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Study Team Members and Contact Information

Chuck Boucher (EAS Fish) 614-836-8776

Chuck McKnight (EAS Macroinvertebrate) 614-836-8784

Rachel Taulbee (Team Leader) 740-380-5433

Kelly Capuzzi (District WQ) 740-380-5283

Randy Spencer (District WQ) 740-380-5240

Josh Griffin (Modeling) 614-644-2874

Emily Poling (District NPDES) 740-380-5425

Tim Fulks (District NPDES) 740-380-5418

Gary Stuhlfauth (CO NPDES) 614-644-2026

Jennifer Witte (District Enforcement) 740-380-5206

Gregg Sablak (TMDL Coordinator) 614-644-4132

Darla Peelle (PIC) 614-728-0035

Dina Pierce (PIC) 614-728-0029

Chris Skalski (Standards) 614-644-2144

Jeff Reynolds (Quality Management) 614-705-1011

Maggie Selbe (401/Wetlands) 740-380-5226

Gary Klase (Fish Tissue) 614-644-2865

Janet Barth (District Drinking Water) 740-380-5250

Heather Raymond (CO Drinking Water) 614-644-2911

19

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Law Enforcement Contact Information

Tuscarawas County Sheriff (330) 339-7743 Coshocton County Sheriff (740) 622-2411

Tuscarawas County Wildlife (330) 245-3047 Coshocton County Wildlife (740) 589-9982

Local Hospitals

Trinity Hospital Twin City, 819 N 1st Street, Dennison, Ohio 44621

20

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

Local Hospitals

Coshocton County Memorial Hospital, 1460 Orange Street, Coshocton, Ohio 43812

21

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

References

Hughes, R. M., D. P. Larsen, and J. M. Omernik. 1986. Regional reference sites: a method for assessing stream pollution. Env. Mgmt. 10(5): 629-635.

Karr, J.R. and D.R. Dudley. 1981. Ecological perspective on water quality goals. Env. Mgmt. 5(1): 55-68.

Ohio Environmental Protection Agency. 2015a. Updates to Biological criteria for the protection of aquatic life: Volume II and Volume II Addendum. User’s manual for biological field assessment of Ohio surface waters. Div. of Surface Water, Ecological Assessment Section, Columbus, Ohio. http://epa.ohio.gov/Portals/35/documents/BioCrit88_Vol2Updates_May2015.pdf

Ohio Environmental Protection Agency. 2015b. Biological criteria for the protection of aquatic life: Volume III. Standardized biological field sampling and laboratory methods for assessing fish and macroinvertebrate communities. Div. of Surface Water, Ecological Assessment Section, Columbus, Ohio. http://epa.ohio.gov/portals/35/documents/BioCrit15_Vol3.pdf

Ohio Environmental Protection Agency. 2015c. Surface water field sampling manual. Version 5.0. July 31, 2015. Div. of Surface Water, Columbus, Ohio. 45pp. http://epa.ohio.gov/Portals/35/documents/SW%20Sampling%20Manual%202015%20Update%20Final%2 0Main.pdf

Ohio Environmental Protection Agency. 2015d. Surface water field sampling manual- Appendix II. Version 5.0. July 31, 2015. Div. of Surface Water, Columbus, Ohio. 37pp. http://epa.ohio.gov/Portals/35/documents/SW%20Sampling%20Manual%202015%20Update%20Final%2 0Appendix%20II.pdf

Ohio Environmental Protection Agency. 2015e. Surface water field sampling manual- Appendix III. Version 5.0. July 31, 2015. Div. of Surface Water, Columbus, Ohio. 52pp. http://www.epa.ohio.gov/Portals/35/documents/SW%20Sampling%20Manual%202015%20Final%20App endix%20III.pdf

Ohio Environmental Protection Agency. 2015f. Surface water field sampling manual- Appendix I. Version 5.0, July 31, 2015. Div. of Surface Water, Columbus, Ohio. 59pp. http://epa.ohio.gov/Portals/35/documents/Inland_Lake_Sampling_Manual.pdf

Ohio Environmental Protection Agency. 2012. Fish Tissue Collection Manual, cooperative fish tissue monitoring program, State of Ohio. Div. of Surface Water. April 12, 2012. 21pp. www.epa.state.oh.us/portals/35/fishadvisory/FishCollectionGuidanceManual12.pdf

Ohio Environmental Protection Agency. 2006. Methods for assessing habitat in flowing waters: Using the Qualitative Habitat Evaluation Index (QHEI). Ohio EPA Tech. Bull. EAS/2006-06-1. Div. of Surface Water, Ecological Assessment Section, Columbus, Ohio.

Ohio Environmental Protection Agency. 1997. Draft. Biological Criteria for the Protection ofAquatic Life: Volume IV: Fish and Macroinvertebrate Indicies for Ohio’s Lake Erie Nearshore Waters, Harbors, and Lacustuaries. Div. of Surface Water, Ecological Assessment Section, Columbus, Ohio.

22

Version 4.4 2017 Lower Tuscarawas River & Select Tributaries Study Plan July 11, 2017

References

Ohio Environmental Protection Agency. 1987. Biological criteria for the protection of aquatic life: Volume II. User’s manual for biological field assessment of Ohio surface waters. Div. of Water Quality Monitoring and Assessment, Surface Water Section, Columbus, Ohio. http://epa.ohio.gov/Portals/35/documents/Vol2.pdf

Ohio Environmental Protection Agency. 1989a. Addendum to biological criteria for the protection of aquatic life: User’s manual for biological field assessment of Ohio surface waters. Div. of Water Quality Planning and Assessment, Surface Water Section, Columbus, Ohio. http://epa.ohio.gov/portals/35/documents/BioCrit88_Vol2Addendum.pdf

Ohio Environmental Protection Agency. 1989b. The qualitative habitat evaluation index (QHEI): rationale, methods, and application. Div. of Water Quality Planning and Assessment, Columbus, Ohio. http://epa.ohio.gov/Portals/35/documents/QHEI_1989.pdf

Omernik, J.M. 1987. Ecoregions of the conterminous . Ann. Assoc. Amer. Geogr. 77(1): 118-125.

23