r% rpA United States Cl#\ Environmental Protection Agency Region 4 Water Division EPA 904-R21 -002 June 2021 Everglades Regional Environmental Monitoring and Assessment Program - REMAP Water Quality and Management, Nutrients, Mercury, Soils and Habitat Monitoring for Adaptive Management: A Phase IV Update U. S. Environmental Protection Agency Region 4 Water Division Atlanta, Georgia June 2021 ------- Everglades Regional Environmental Monitoring and Assessment Program - REMAP Water Quality and Management, Nutrients, Mercury, Soils and Habitat Monitoring for Adaptive Management: A Phase IV Update U.S. Environmental Protection Agency Region 4 Water Division Atlanta, Georgia June 2021 This document is available at: https://www.epa.aov/everalades/environmental-monitorina-everalades ------- Please reference this report as follows: Scheldt, D. J., P. I. Kalla and D. D. Surratt. 2021. Everglades Regional Environmental Monitoring and Assessment Program - REMAP. Water quality and management, nutrients, mercury, soils and habitat. Monitoring for adaptive management: a phase IV update. EPA 904-R21-002. USEPA Region 4, Atlanta, Georgia. ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) TABLE OF CONTENTS SUMMARY i ABBREVIATIONS iv ACKNOWLEDGMENTS vi INTRODUCTION AND PURPOSE... 1 BACKGROUND 4 The Everglades 4 An Altered System 5 THE COMPREHENSIVE EVERGLADES RESTORATION PLAN (CERP) 8 EVERGLADES REGIONAL ENVIRONMENTAL MONITORING AND ASSESSMENT PROGRAM (REMAP) 11 Program Design 11 Probability Samples 12 Everglades REMAP History 15 Phase IV Sampling 16 Data Quality Assurance 17 Data Uses 18 Data Analysis and Presentation 22 WATER CONDITIONS 25 WATER QUALITY. 28 Nutrients 28 Phosphorus 29 Soil Phosphorus 35 Nitrogen 39 Specific Conductance 42 Chloride 47 Sulfur 50 Carbon 58 Dissolved Oxygen 61 pH 63 SOILS AND SOIL SUBSIDENCE 65 MACROPHYTES AND PERIPHYTON 75 Plant Communities 75 Periphyton 77 Floe 78 MERCURY 79 LITERATURE CITED 88 EVERGLADES REMAP PROGRAM REPORTS AND PUBLICATIONS 107 APPENDICES 110 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) SUMMARY The United States Environmental Protection Agency's Everglades Regional Environmental Monitoring and Assessment Program (REMAP) is a comprehensive, long-term monitoring and assessment effort. Its goal is to provide critical scientific information needed for management decisions on the Everglades ecosystem and its restoration. Since 1993, four phases of marsh sampling and one phase of canal sampling have been conducted throughout the Everglades landscape at 1250 different locations. REMAP is unique to the Everglades in consistently combining several aspects of scientific study: a probability-based sampling design which results in quantitative statements across space about the condition of the Everglades; multi-media sampling; and extensive spatial coverage that includes all of the freshwater Everglades. Everglades REMAP: contributes to understanding the effectiveness of phosphorus control efforts; contributes to the joint federal-state Comprehensive Everglades Restoration Plan (CERP) by quantifying conditions in three physiographic regions: Everglades ridge and slough; Everglades marl prairie/rocky glades; and Big Cypress Swamp; provides information on four groups of Everglades restoration success indicators: surface water, soil and sediment, vegetation, and fish; provides a baseline against which future conditions and the effectiveness of restoration efforts can be measured; assesses the effects and potential risks of multiple environmental stressors on the Everglades ecosystem, such as water management, soil loss, water quality degradation, habitat loss, and mercury contamination; and provides data with many applications: - document water conditions, soil thickness and subsidence; - document landscape patterns of water quality conditions and soil nutrients and contaminants; - understand water management impacts on water quality; - determine which portions of the Everglades are phosphorus-impacted; - understand sulfur sources and distribution; - document mercury conditions, landscape patterns, mass balances, biomagnification factors and identify environmental conditions associated with high mercury in prey fish; - map vegetation and determine vegetation classes, biomass, standing stocks and presence of exotic species; - determine landscape patterns of periphyton; - and understand spatial variation in aquatic food webs. This report summarizes results for REMAP's September 2014 wet season Phase IV biogeochemical sampling, which documented conditions during a two-week sampling snapshot throughout the 2,098-square-mile freshwater portion of the Everglades. The 2014 conditions are also compared to the conditions observed by REMAP during the 1990s and 2005. ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Key findings: Less phosphorus enrichment and improved phosphorus conditions: Surface water phosphorus was lower in 2014 as compared to 2005. During 2014, 14.8 ± 4.8% of the Everglades marsh had a surface water phosphorus concentration greater than 10 parts per billion, as compared to 27.2 ± 7.5% in 2005. This improvement is due to the combination of agricultural best management practices and Florida's Stormwater Treatment Areas. REMAP data also indicate no change in soil phosphorus conditions in 2014 compared to 2005. During 2014, 21.1 ± 5.6% of the Everglades had soil phosphorus exceeding 500 milligrams per kilogram (mg/kg), Florida's definition of "impacted", as compared to 24.5 ± 6.4% in 2005. REMAP previously found that the area with soil phosphorus >500 mg/kg expanded from the 1990s to 2005. Mercury was lower in mosquitofish and water, but was still elevated: The mercury concentration in mosquitofish was lower during the 2014 REMAP wet season sampling event, as compared to 2005 and 1995. Mosquitofish are a key food item for other fish, which then are consumed by Everglades gamefish as well as wading birds. Therefore, mosquitofish are relevant to human health and ecological health. During the 2014 wet season, 13.0 ± 5.7% of the marsh had mosquitofish that exceeded 77 parts per billion, the maximum concentration USEPA recommends in trophic level 3 fish as being protective of top predators such as birds and mammals, in contrast to 64.7 ± 7.3% in 2005 and 70.5 ± 7.1% in 1995. Concentrations of methylmercury and total mercury in surface water were also lower in 2014 than in 2005. The biomagnification of methylmercury from surface water to mosquitofish continues to be high. Florida's fish consumption advisories for gamefish to protect human health are still in effect. No further indication of soil loss in the public Everglades: There was no indication of further soil loss overall in the Everglades in 2014 as compared to 2005 and 1995-96. REMAP previously found that from 1946 to 1996 about one-half of the peat soil was lost from approximately 200,000 acres of the public Everglades that were subjected to drier conditions. An inch of peat soil that took centuries to form can be lost within a few years, or within hours if dry soils are subjected to fire. In 2014, about 25% of the greater Everglades had 1.0 feet or less of soil, as did 52% of Everglades National Park. The median soil thickness in the Everglades was 2.3 feet. The volume of soil remaining in the Everglades was 4.7. x 109 cubic meters. Water management structural and operational changes would help to maintain the remaining marsh soils in drier areas, along with the plant communities and wildlife habitats of these wetlands. The northern portion of WCA3 must dry out less frequently if further soil loss is to be prevented. This is a focus of the Central Everglades Planning Project. Pronounced water quality gradients: There are pronounced spatial gradients in surface water sulfate, chloride, specific conductance, carbon, nitrogen and phosphorus in the Everglades marsh. These gradients are due to the relative contribution of rainwater, stormwater and groundwater and the proximity to canals. The highest concentrations typically occurduring the wet season in WCA2, due to its proximityto the Everglades Agricultural Area and stormwater discharges. Concentrations generally decrease to natural background levels as water moves downstream through the Everglades. Years with higher discharge tend to have higher concentrations. Location, time of year, rainfall and water management practices are important factors that affect water quality. ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Canals are a conduit for transport of degraded water: The canal system, constructed primarily for flood control and water supply, is an inadvertent conduit for the transport of degraded water with elevated sulfate, chloride, specific conductance, carbon, phosphorus and nitrogen into the Everglades marsh. Water management practices affect water quality. Downstream water quality could be improved if canals were eliminated, or if they were operated to minimize the influence of canal flow and maximize the sheet flow of marsh surface water, with the diluting influence of rainfall and cleaner marsh water, as is proposed by various restoration efforts. Regardless, pollutant control is usually most effective at the source. Sulfate was lower in surface water: Surface water sulfate was lower during the 2014 REMAP wet season sampling event, as compared to 2005 and 1995. In 2014, about 37.1 ± 6.0% of the Everglades marsh had a surface water sulfate concentration exceeding 1.0 parts per million (ppm), the Everglades restoration goal, as compared to 57.3 ± 6.0% in 2005. Interior portions of the Everglades, distant from stormwater discharges from the Everglades Agricultural Area, had concentrations as low as 0.02 milligrams per liter, although some elevated concentrations were still found as far south as Shark River Slough within the Park. Water quality and soil conditions vary by location andtime: Waterquality conditions in the Everglades vary greatly with location. Rainfall-driven areas that are distant from the influence of canal water, such as the interior of the Refuge and the southwest portion of WCA3, have good water quality and low soil phosphorus. The interior of the Refuge tends to have good waterquality and the lowest phosphorus concentrations observed in peat soils. In contrast, northern WCA3 has thinner soil because of drainage, elevated soil phosphorus, and extensive cattail encroachment. WCA2 has phosphorus enrichment and cattail encroachment, along with higher surface water sulfate, chloride, conductivity, nitrogen and organic carbon. Water depth at any given location varies with season and year. Environmental threats are interrelated: Environmental stresses such as water management, soil loss, waterquality degradation, cattail expansion, and mercury contamination are often interrelated. Efforts to manage water quantity and restore the Everglades should be integrated with efforts to manage or control phosphorus, mercury and sulfur. Management and restoration options should be assessed from a holistic perspective. For three decades beginning in the 1990s, Everglades REMAP has provided monitoring and assessment data for measuring ecosystem health and the effectiveness of Everglades restoration. As CERP restoration efforts and Everglades phosphorus and mercury control efforts proceed, this probability-based sampling can be repeated to quantitatively document the condition of the Everglades and the effectiveness of these efforts. This report describes the condition of the Everglades during the intensive 2014 marsh sampling effort. All REMAP data and reports are available at: https://www.epa.aov/everalades/environmental-monitorina-everalades. iii ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) ABBREVIATIONS cc = cubic centimeters, cm3 cm = centimeters ft = feet g = grams g/cc = grams per cubic centimeter in = inches km = kilometers Ibs/ac = pounds per acre n = number or count ppb = parts per billion (|jg/L) ppm = parts per million (mg/L) or (mg/kg) mg/kg = milligrams per kilogram (parts per million, ppm) mg/L = milligrams per liter (parts per million, ppm) ng/g = nanograms per gram (parts per billion, ppb) ng/L = nanograms per liter (parts per trillion, ppt) jjg/cc = micrograms per cubic centimeter |jg/g = micrograms per gram (parts per million, ppm) jjg/kg = micrograms per kilogram (parts per billion, ppb) |jM = micromolar = 1 micromole per liter (jjmol/L) jjmhos/cm = micromhos per centimeter AA = Alligator Alley (Interstate 75) ASR = Aquifer Storage and Recovery BCFm = Bioconcentration Factor BCNP = Big Cypress National Preserve BMPs = Best Management Practices CD = Consent Decree CDF = Cumulative Distribution Function CERP = Comprehensive Everglades Restoration Plan CEPP = Central Everglades Planning Project CI = Confidence Interval DO = Dissolved Oxygen DOC = Dissolved Organic Carbon DOM = Dissolved Organic Matter EAA = Everglades Agricultural Area ECP = Everglades Construction Project EDEN = Everglades Depth Estimation Network EFA= Everglades Forever Act ENP = Everglades National Park EMAP = Environmental Monitoring and Assessment Program EMM = Eastern Marl Marsh; E = East; W = West EPA = Everglades Protection Area FDEP = Florida Department of Environmental Protection FDHRS = Florida Department of Health and Rehabilitative Services FDOH = Florida Department of Health FIU = Florida International University FWM = Flow-Weighted Mean GIS = Geographic Information System I-75 = Interstate 75 (Alligator Alley) iv ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) LNWR = Arthur R. Marshall Loxahatchee National Wildlife Refuge LSASD = USEPA Region 4 Laboratory Services and Applied Science Division MDL = Method Detection Limit MeHg = Methylmercury N = Nitrogen NADP = National Atmospheric Deposition Program NRC = National Research Council NESS = Northeast Shark Slough OFW = Outstanding Florida Water OMM = Ochopee Marl Marsh P = Phosphorus Park = Everglades National Park QA= Quality Assurance QAPP = Quality Assurance Project Plan RECOVER = Restoration Coordination and Verification Refuge = Arthur R. Marshall Loxahatchee National Wildlife Refuge REMAP = Regional Environmental Monitoring and Assessment Program S = Sulfur SESD = USEPA Region 4 Science and Ecosystem Support Division SRS = Shark River Slough SFWMD = South Florida Water Management District SSAC = Site-Specific Alternative Criterion STA = Stormwater Treatment Area THg = Total Mercury TOC = Total Organic Carbon TMDL = Total Maximum Daily Load TN = Total Nitrogen TP = Total Phosphorus TS = Taylor Slough USACE = United States Army Corps of Engineers USDOI = United States Department of Interior USEPA = United States Environmental Protection Agency USFWS = United States Fish and Wldlife Service USGS = United States Geological Survey WCA = Everglades Water Conservation Area 2A, 2B, 3A, 3B; 3AN = WCA 3A north of 1-75; 3AS = WCA 3A south of 1-75 WY = Water Year, May to April (e.g., WY 16 is May 1, 2015 to April 30, 2016) v ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) ACKNOWLEDGMENTS The United States Environmental Protection Agency's Regional Environmental Monitoring and Assessment Program (REMAP) is a comprehensive, long-term monitoring and assessment effort. Its goal is to provide critical scientific information needed for management decisions on the Everglades ecosystem and its restoration. This large undertaking would not have been possible without the efforts and collaborative support of many people. The United States Department of Interior and Everglades National Park provided helicopter operations support. Sampling permits or access were received from the Miccosukee Tribe of Indians of Florida, the Arthur R. Marshall Loxahatchee National Wildlife Refuge, Big Cypress National Preserve and Everglades National Park. This report was reviewed by USEPA scientists and program experts (Chris Decker, Susan Dye, Morris Flexner, Cecelia Harper, Michele Wetherington). This report benefited from peer reviews provided by scientists at the South Florida Water Management District (Dr. Binhe Gu, Dr. Nenad Iricanin, Dr. Susan Newman, Dr. Fred Sklar), Florida Department of Environmental Protection (Dr. Paul Julian), Arthur R. Marshall Loxahatchee National Wildlife Refuge (Dr. Rebekah Gibble, Rolf Olson), United States Army Corps of Engineers (Dr. Gretchen Ehlinger), and the United States Geological Survey (Dr. Nicholas Aumen, Dr. David Krabbenhoft, Dr. Wlliam Orem). Funding: Extramural funding for REMAP IV was provided by the South Florida Ecosystem Office of the National Park Service under Interagency Agreement # P13PG7. USEPA funding was provided by the Office of Water and the Region 4 Water Division's South Florida Geographic Initiative. The following individuals helped make Everglades REMAP IV possible: Alion, Inc. Don Fortson - field support team leader Sue Jones - data reviewer Michael Keller - data reviewer Yvette Lane-Walcott - laboratory analysis Nathan Mangle - instrumentation manager Ray Popovic - statistician Louie Pounds - field chemist Florida International University Yong Cai - principal investigator, mercury Daniel Gann - principal investigator, remote sensing Diana Johnson - analytical laboratory Ruth Justiniano - quality assurance Mark Kershaw - analytical laboratory Guangliang Liu - principal investigator, mercury Pedro Lorenzo - analytical laboratory Jenny Richards - principal investigator, vegetation Len Scinto - principal investigator, nutrients and soil properties HMC Helicopters, Inc. Gary Freeman - pilot Jorge Gomes - pilot Carlos Luque - pilot Miccosukee Tribe of Indians of Florida James Erskine - acting director, water resources US Department of the Interior US Fish and Wildlife Service - Arthur R. Marshall Loxahatchee National Wildlife Refuge Rebekah Gibble - helicopter operations Marcie Kapsch - permitting vi ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) National Park Service - Big Cypress National Preserve Steve Schulze - permitting National Park Service - Everglades National Park Rick Anderson - helicopter operations NickAumen - administration Clayton Camblin - helicopter operations and safety training GaiV Carnall - helicopter operations Joffre Castro - project officer Andrew Gill - helicopter operations Kim Gomez - administration Bob Johnson - administration Carol Mitchell - administration Donatto Surratt - project officer Heather Walker - administration P. J. Walker - permitting United States Geological Survey Bryan McCloskey - Everglades Depth Estimation Network (EDEN) Office of Aircraft Services Margaret Gallagher - helicopter safety training USEPA Office of Research and Development National Health and Environmental Effects Research Laboratory Tom Kincaid - program sampling design and statistical consulting Tony Olsen - program sampling design and statistical consulting Region 4 Office of Regional Counsel Michele Wetherington - report review Region 4 Water Management Division Jon Becker - GIS John Davis - data statistics Cecelia Harper - report review Tara Houda - field base support and sampling Robin Polinsky - GIS Daniel Scheidt - associate program leader Jennifer Shadle - financial contracting Eric Somerville - field base support and sampling Region 4 Science & Ecosystem Support Division Jerry Ackerman - field sampling, equipment design and fabrication, instrumentation, logistics Danny Adams - laboratory analysis Sandra Aker - data quality assurance (QA) Nathan Barlet - field data and statistics review Pam Betts - laboratory analysis Tony Carroll - laboratory analysis Chris Decker - field sampling, equipment preparation, project QAPP, QA sample booking, report review Lonnie Dorn - preparation of field supplies, field base support Sue Dye - project QAPP, field data form and stepwise procedure, field sampling, report review Morris Flexner- procurement and preparation of field supplies, field sampling, report review Cornell Gayle - field base support Denise Goddard - data QA Jeff Hendel - project laboratory QA officer vii ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Hunter Johnson - field sampling Peter Kalla - program leader Bobby Lewis - field data quality Derek Little - field equipment procurement and preparation, field sampling, data management Art Masters - field QA officer Jon McMahan - procurement, field logbook and equipment preparation, mapping, field sampling Landon Pruitt - field QA officer Michael Roberts - field sampling John Ruiz - field equipment preparation and fabrication, training, sampling, QA Tim Simpson - procurement, field sampling Kevin Simmons - instrumentation, sample manager Brain Striggow - instrumentation, field sampling Ray Terhune - project laboratory QA officer Linda Watson - administration Greg White - field logbook and equipment, aviation safety plan, communications, field sampling John Thomason - laboratory and data QA Jeff Wlmoth - data QA Photograph Acknowledgments: Cover upper right and Figure 15 lower right: Derek Little and Tim Simpson; Cover lower right and Figures 1 and 2 bottom: Susan Dye and Hunter Johnson; Table of Contents person in cattail and Figure 66: Phyllis Meyer; Figure 2 top: Lonnie Dorn and Morris Flexner; Figures 3 and 23: Everglades National Park; Figure 9 satellite base map and Figure 24: South Florida Water Management District; Figure 16 upper right: Mel Parsons; Figure 58: Jerry Ackerman, Susan Dye, Morris Flexner, Derek Little, Hunter Johnson, Michael Roberts, John Ruiz and Brian Striggow; Figure 64 right: Morris Flexner and Derek Little; all other photographs Daniel Scheidt. viii ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) INTRODUCTION AND PURPOSE The United States Environmental Protection Agency (USEPA) Everglades Regional Environmental Monitoring and Assessment Program (REMAP; Program) is a unique, long-term, monitoring, assessment and research effort. Its goal is to provide scientific information that is needed for decisions about the protection, restoration and management of the Everglades ecosystem. REMAP data have been used by over 30 entities including federal and state agencies, Indian Tribes, agriculture, the public, non-governmental organizations, and the National Academies of Sciences to document conditions and determine whether they are improving, are unchanged, or are worsening. Data also help to assess restoration progress and understand the effectiveness of control programs for phosphorus and mercury. Since 1993, one phase of canal sampling and four phases of marsh sampling have been conducted throughout the Everglades landscape at 1250 sampling locations (Figures 1 and 2). The purpose of this report is to document conditions in the Everglades during 2014, the fourth phase of marsh sampling, and make statements about whether conditions in the Everglades were better, unchanged, or worse compared to 1995 and 2005. REMAP is unique to the Everglades in that it consistently combines several aspects of scientific study - a probability-based sampling design which results in quantitative statements across space about the condition of the Everglades, multi-media sampling (water, soil, biota) and an extensive spatial coverage that includes all of the freshwater Everglades. Figure 1. Everglades wet prairie and sawgrass marsh mosaic. Tree islands are visible on the horizon. Numerous environmental issues threaten the Everglades "River of Grass," such as water management, soil loss, water quality degradation, and habitat alteration. 1 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Goal: Provide environmental information that contributes to environmental management decisions on Everglades protection and restoration. Everglades REMAP: documented pre-restoration conditions in the Everglades during the 1990s, as well as conditions during more recent decades after initiation of restoration efforts; contributes to understanding the effectiveness of phosphorus and mercury control efforts; contributes to the joint federal-state Comprehensive Everglades Restoration Plan (CERP) by quantifying conditions in three physiographic regions: Everglades ridge and slough (Figures 1 and 2); Everglades marl prairie/rocky glades; and Big Cypress Swamp; provides information on four groups of Everglades restoration success indicators (media): surface water, soil and sediment, vegetation, and fish; simultaneously assesses the effects and potential risks of multiple environmental stressors on the Everglades ecosystem, such as water management, soil loss, water quality degradation, nutrient enrichment, habitat loss, and mercury contamination; permits spatial analyses and identifying associations that provide insight into relationships among environmental stressors and observed ecological responses; and provides data with many applications: - document water conditions, soil thickness and subsidence; - document landscape patterns of water quality conditions and soil nutrients and contaminants; - understand water management impacts on water quality; - determine which portions of the Everglades are phosphorus-impacted; - understand sulfur sources and distribution; - document mercury conditions, landscape patterns, mass balances, biomagnification factors and identify environmental conditions associated with high mercury in prey fish; - map vegetation and determine vegetation classes, biomass, standing stocks and presence of exotic species; - determine landscape patterns of periphyton; - understand spatial variation in aquatic food webs. REMAP has been carried out in cooperation and coordination with the United States Fish and Wildlife Service, National Park Service, Florida Department of Environmental Protection, United States Army Corps of Engineers, Miccosukee Tribe of Indians of Florida, Seminole Tribe of Indians, Florida International University, United States Geological Survey, and the South Florida Water Management District. 2 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 2. The Everglades marsh sawgrass and open water mosaic: ground view (top) and aerial view (bottom). Forthree decades beginning in the 1990s, Everglades REMAP has provided monitoring and assessment data for measuring ecosystem health and the effectiveness of Everglades restoration. As CERP restoration efforts and Everglades phosphorus and mercury control efforts proceed, this probability-based sampling can be repeated to quantitatively document the condition of the Everglades and the effectiveness of these efforts. This report describes the condition of the Everglades during the intensive 2014 marsh sampling effort. All REMAP data and reports are available at https://www.epa.aov/everalades/environmental-monitorina- everalades. 3 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) BACKGROUND The Everglades "Here are no lofty peaks seeking the sky, no mighty glaciers or rushing streams wearing away the uplifted land. Here is land, tranquil in its quiet beauty, serving not as a source of water but as a last receiver of it." "The Everglades were not really set aside for any kind of geological wonders or scenic features. It's the first national park set aside simply for its wildlife and the plants and trees - for its biological diversity." President Harry Truman, Everglades National Park dedication, 1947. The Florida Everglades is one of the largest freshwater marshes in the world (Ramsar2006), consisting of a unique mosaic of sawgrass, wet prairies, sloughs, and tree islands (Figures 1 and 2). Just over 100 years ago, this vast wilderness encompassed over 4,000 square miles, extending 100 miles from the shore of Lake Okeechobee south to Florida Bay. The intermingling of temperate and Caribbean flora created habitat for a variety of fauna, including alligators, Florida panthers and hundreds of thousands of wading birds. The Everglades were defined by several characteristics: How the water flowed. Water connected the ecosystem, from north to south. Surface water flowed freely and slowly across the flat landscape. Rainfall during the wet season took months to move through the marsh. The large water storage capacity and the slow flow made wetlands and coastal waters less vulnerable to South Florida's variable and often intense rainfall (USACE and SFWMD 1999). Vastness. The large area provided a variety of wildlife habitats. Millions of acres of wetlands provided large feeding ranges and diverse habitat for wildlife. The vastness produced abundant aquatic life while facilitating recovery from hurricanes, fires, and other natural disturbances (USACE and SFWMD 1999). There was no development, so there was no need to provide flood control for agricultural or urban areas. Diverse mosaic of landscapes. The Everglades was a complex system of plants and animals dictated in part by varied water regime - minimum, average, and maximum water depths, the duration of surface water inundation, and the slow, important flow of water that determined the ridge and slough landscape. This complex water regime resulted in diverse, expansive areas of sloughs, wet prairies, sawgrass marshes, cypress swamps, mangrove swamps, coastal lagoons and bays (USACE and SFWMD 1999). Natural water quality conditions. Centuries ago, there were no external sources of pollutants, either from surface water or the atmosphere. There was no urban development or agriculture. Nutrients such as phosphorus and nitrogen, ions such as sulfate and chloride, and metals such as mercury all occurred at natural background levels. Clean rainfall recharged groundwater during the dry season and generated surface water, which defined soils and the natural mosaic of plant communities. There were no canals to connect groundwater and surface water, or to serve as a conduit for pollutants. The slow, shallow flow of 4 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) surface water across the landscape provided ample opportunity for cleansing by extensive wetlands. The wet prairies and sawgrass marshes of the Everglades developed under and were defined by extremely low phosphorus concentrations. The mosaic of habitats, their vastness and the variety of water patterns supported the long-term survival of wildlife under a range of seasonal and annual water conditions, and across dry and wet years and occasional extreme events such as tropical storms. An Altered System Today, 50% of historic Everglades wetlands have been irreversibly drained. The Everglades ecosystem has been altered by extensive agricultural and urban development (Figures 5 to 7). South Florida's human population, which was 8.7 million in 2020 (SFWMD https://www.sfwmd. gov/who-we-are/chairmans-messaae accessed 10/30/2020), continues to increase, requiring more water and an expanding area needing flood control. This population is projected to increase to 15 million within a few decades (USACE and SFWMD 1999) (Figure 8). The Everglades landscape changed during the twentieth century as drainage canals were dug to facilitate development. Most of the remaining Everglades are in the Everglades Protection Area (EPA): Arthur R. Marshall Loxahatchee National Wildlife Refuge (LNWR or the Refuge), the Water Figure 4. Everglades wading bird populations significantly declined during the 1900s. A focus of Everglades restoration is getting the water quantity, timing, distribution and quality right so wading birds can thrive. One century ago, the greatest threat to wading bird populations was hunting (Figures 3 and 4). During the 1900s, the Everglades became an altered system. In response to periods of drought in the 1930s and 1940s, and severe flooding with loss of human life in the 1920s and 1940s, the Central and Southern Florida Project for Flood Control and other Purposes (the Project) was authorized in 1948 by federal legislation. The Project's often conflicting purposes include flood control, water supply, water conservation, prevention of salt water intrusion, and preservation offish and wildlife. The Project is one of the world's most extensive public water management systems, consisting of over 1,800 miles of levees and canals, 25 major pumping stations, and over 200 large and 2,000 smaller water control structures. When the Project began its design in the 1950s, about 500,000 people lived in the region and it was estimated that there might be two million people by 2000 (USACE and SFWMD 1999). The Project has effectively provided the flood control and water supply that has facilitated urban and agricultural growth and the resulting economic benefits. Figure 3. Decorating women's hats with wading bird plumage led to the decimation of Everglades wading bird populations around 1900. 5 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) During the last century, the Everglades became subject to multiple, often interrelated, environmental threats. Effective ecosystem protection and restoration requires addressing these threats holisticaliy. Figure 7. Residential development on former Everglades wetlands in western Dade County, 2005. Figure 8. South Florida population in millions from 1900-2050 (USACE and SFVVMD 1999 and US Census Bureau). Conservation Areas (WCAs), and Everglades National Park (ENP or the Park) (Figure 9). The Park, which was established in 1947, includes only one-fifth of the original "River of Grass" that once encompassed over 4,000 square miles (2 million acres) (Davis and Ogden 1994). One-fourth of the historic Everglades is in agricultural production within the 1000-square-mile Everglades Agricultural Area (EAA), where sugarcane and vegetables are grown on the rich, productive peat soils of drained sawgrass marshes. Another one-fourth of the historic Everglades has been drained and converted into urban areas along Florida's lower east coast. The Everglades watershed begins near Orlando 100 miles north of Lake Okeechobee. Although one- third of the 16,000-square-mile Everglades watershed is in public ownership, there are many environmental issues, often interrelated, that must be resolved to protect and restore the Everglades ecosystem. These include: water management complexities; water supply and timing conflicts; loss of water storage capacity; soil loss; water quality degradation and eutrophication; unnatural changes in plant communities; habitat alteration and loss (Sklar et al. 2019); mercury contamination of game fish and wildlife such as wading birds and Florida panthers; protection of endangered species; and introduction and spread of nuisance exotic species of plants and animals. Figure 5. Urban expansion into drained Everglades wetlands within western Broward County, 1995.. Note the black peat soil. Figure 6. About 400,000 acres of sugarcane are grown on the peat soils of former Everglades wetlands within the EAA. 6 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 9. Satellite image of South Florida, circa 1995, with the areas sampled by REMAP outlined in yellow: Everglades Agricultural Area (EAA); Arthur R. Marshall Loxahatchee National Wildlife Refuge (Refuge); Everglades Water Conservation Area 2 (WCA2A and WCA2B); Everglades Water Conservation Area 3 north of I-75 (WCA3A N); Everglades Water Conservation Area 3 south of I-75 (WCA3A S and WCA3B); the eastern portion of Big Cypress National Preserve, and the freshwater portion of Everglades National Park (Park). Light areas on the east are urban development. The black line approximates the extent of the pre-1900 Everglades marsh. The Everglades watershed begins north of Lake Okeechobee. 7 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Historic Flow Figure 10. Historic Everglades water flow patterns (left) and present flow patterns (right) (adapted from USAGE and SFWMD 1999)- Figure 11. An extensive system of canals, levees, and water control structures has modified Everglades water conditions and provides an inadvertent conduit for pollutant transport. The S9 pump station (foreground) in Broward County discharges stormwater westward from an urban basin into the Everglades (background). THE COMPREHENSIVE EVERGLADES RESTORATION PLAN (CERP) The Central and Southern Florida Project has provided flood protection and water supply for urban and agricultural lands, as intended, and fostered growth. Simultaneously, the Project has altered the Everglades ecosystem. Some of the Everglades no longer receives the right quality or quantity of water at the right place or the right time. The remnant Everglades no longer exhibits the water regimes and patterns of flow, vast area, and mosaic of habitats that defined the pre-drainage, natural ecosystem. Wildlife habitat has been lost or changed, and the number of wading birds (wood stork, great egret, snowy egret, tricolored heron, and white ibis) nesting in the Everglades decreased markedly during the twentieth century (Ogden 1994). Historically, most water slowly flowed across or soaked into the region's vast wetlands. Today, over one-half of Everglades wetlands have been irreversibly drained, with the loss of their functions including water storage and water quality filtration. Flows into the Everglades marsh are at times too much or too little, and at the wrong time (Figures 9 and 10). Some areas are too wet while other areas are too dry. Overland sheet flow of water is interrupted by levees and canals that crisscross the Everglades and can provide a conduit for pollutant transport and release into the marsh (Figure 11). The canal system can quickly drain water from developed areas. Releases of nutrient-rich water from Lake Okeechobee present water quality challenges whether they flow east to the St. Lucie Estuary, west to the Caloosahatchee Estuary, or south to the Everglades. As the human population continues to increase, conflicts over water among natural resources, agriculture, industry, and a growing population may intensify. Present Flow 8 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) GUQ 1 Figure 12. The right quality, quantity, timing and distribution of water are all critical to South Florida ecosystem protection and restoration. ® V r-r-^l U \ £vT I v\T \ nA \ n 1 \ V N \ - niAp I V 1 pi V i \ r \Jr Figure 13. The anticipated effect of the original Comprehensive Everglades Restoration Plan (CERP). Without the Plan (left) restoration targets are not expected to be met (red). With the Plan fully implemented (right) restoration targets are more likely to be met (green). Yellow indicates uncertainty in meeting restoration targets (USACE and SFWMD 1999). Protection and Restoration During the 1970s and 1980s there was growing recognition of threats to the Everglades. Many of the problems with declining ecosystem health revolve around water: water quantity, quality, timing, and distribution (Figure 12). Consequently, a major goal of restoration is to deliver the right amount of water, that is clean enough, to the right places and at the right time. Since water largely defined the natural system, it is expected that the natural system will respond to improvements in water management (Figure 13). The federal Water Resources Development Acts of 1992 and 1996 directed the U.S. Army Corps of Engineers (USACE) to review the Project and develop a comprehensive plan to restore and preserve south Florida's natural ecosystem, while providing for other water-related needs of the region, including urban and agricultural water supply and flood protection. The result is the Comprehensive Everglades Restoration Plan (CERP, or the Plan, http://www.everaladesrestoration.gov/). which was authorized by the United States Congress in the Water Resources Development Act of 2000. The development of the Plan was led by the Army Corps and the South Florida Water Management District (SFWMD) and over 100 ecologists, hydrologists, planners, engineers and other professionals from over 30 federal, state, tribal, local agencies and the public. The Plan includes: about 180,000 acres of surface water storage areas; about 36,000 acres of man-made wetlands to treat urban or agricultural runoff from basins other than the EAA; wastewater reuse; extensive aquifer storage and recovery; water management operational changes; and structural changes to improve how and when water is delivered to 9 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Evaluating restoration success requires reliable pre-restoration and post-restoration information on ecosystem condition. the Everglades, including removal of some of the canals or levees that prevent natural overland sheet flow. If nothing is done, the health of the Everglades will continue to decline, water quality will degrade further, some plant and animal populations will be stressed further, water shortages for urban and agricultural users will become more frequent with economic consequences, and the ability to protect people and their property from flooding will be compromised (USACE and SFWMD 1999). CERP includes over 50 projects. The initial focus is on regaining some of the water storage capacity that was lost with wetland drainage by building water storage reservoirs, restoring water quality with treatment wetlands, restoring surface water sheet flow, and enhancing water management options. The entire Plan will take many decades to implement and cost over $23 billion as of 2020 (USACE and USDOI 2020; $16 billion as of 2014 and $8 billion as of 1999, Congressional Research Service 2017). The cost is generally split equally by the taxpayers of Florida and the United States. A key to Everglades restoration is the Central Everglades Planning Project (CEPP) which includes: increasing storage, treatment and conveyance of water south of Lake Okeechobee; removing canals and levees within the central Everglades; and retaining water within Everglades National Park and protecting urban and agricultural areas to the east from flooding (USACE 2014). Example CERP Everglades Ecosystem Restoration Performance Measures (RECOVER 2007) Water Management Reinstate system-wide natural hydropatterns and sheet flow Habitat Alteration Increase spatial extent of habitat and wildlife corridors Eutrophication Water total phosphorus is not to exceed 10 micrograms per liter and should not exceed OFW concentration levels. Decrease the areal extent of areas with soil phosphorus exceeding 500 milligrams per kilogram (mg/kg) and maintain or reduce to 400 mg/kg or less. Surface water total nitrogen should be < 1994-2004 baseline. Mercury Contamination No statistically significant increase in levels of mercury in fish tissue; state water quality standards will be met Sulfate Contamination Maintain or reduce surface water sulfate to 1 milligram per liter or less Surface Water Specific Conductance No more than 25% increase above background, maintain low specific conductance in the Refuge Increase aerial coverage of habitats that reflect Natural Systems Model Periphyton Soil Loss Restore natural soil formation processes and rates 10 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) EVERGLADES REGIONAL ENVIRONMENTAL MONITORING AND ASSESSMENT PROGRAM (REMAP) Given the multi-billion dollartaxpayer investment in CERP and efforts to control phosphorus and mercury, monitoring and assessment are important. Monitoring data are needed to determine ecosystem conditions, identify threats, detect changes, and evaluate environmental restoration efforts in a holistic manner. Monitoring objectives include: Documenting status and trends; Determining baseline variability; Detecting responses to management actions; Improving the understanding of cause and effect relationships. Accordingly, CERP has adopted an integrated monitoring and assessment plan that includes key performance measures as indicators of ecosystem health. Performance measures are tools to allow the public and managers to predict system-wide performance of alternative plans and to assess actual performance following implementation. Achieving targets for performance measures is expected to result in system-wide sustainable restoration (RECOVER 2007). Everglades REMAP data have been used by CERP for performance measures, the monitoring and assessment plan, the system status report, and the 5-year Report to Congress (USACE and USDOI 2015). Program Design The attention and funding devoted to Everglades ecosystem protection and restoration are unprecedented. It is important to regularly and comprehensively assess ecosystem health in a cost-effective, quantitative manner. Such an assessment identifies resource conditions and restoration needs, and allows one to determine the effectiveness of restoration and pollution control efforts. A defining feature of the Everglades is its large spatial area. Therefore, to document conditions and track restoration it is important to accurately determine the proportion of the current Everglades that is subject to various impacts or stressors. The Everglades Regional Environmental Monitoring and Assessment Program (REMAP) permits a holistic view of indicators of ecological condition throughout the freshwater Everglades landscape. An indicator is a measurable characteristic of the environment, abiotic or biotic, that can provide information on the condition of ecological resources. REMAP's large-scale perspective is critical to understanding the impacts of different factors (such as phosphorus, mercury, sulfur, habitat alteration, or hydropattern modification) on the entire Everglades landscape, rather than at individual locations or in small areas. Looking only at isolated sites in a specific area and extrapolating to the larger ecosystem can be misleading. 11 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Probability Samples REMAP employs a statistical, probability-based sampling strategy, similar to polls or surveys. This approach was initiated throughout the United States in the early 1990s by USEPA (USEPA1993, Thornton et al. 1994, Diaz-Ramos et al. 1996, Stevens 1997). Indicators of pollutant exposure and habitat condition can be used to identify associations between human-induced stresses and ecological condition. This design has been reviewed by the National Academies of Sciences, and USEPA has applied it to lakes, rivers, streams, wetlands, estuaries, forests, arid ecosystems and agro-ecosystems throughout the United States (Olsen et al. 1999, USEPA 1993, 1995). In a probability sample, or survey, every member of the population has a known chance of being selected as the samples are drawn at random. The probability design USEPA uses to sample the Everglades was developed from the Environmental Monitoring and Assessment Program (EMAP) base grid, a generalized random-tessellation stratified approach (Stevens and Olsen 2004). The design includes locating stations separately within each of the four major Everglades subareas in order to assure a sufficient number of stations in the smaller subareas (the Refuge and WCA2, as compared to WCA3 or the Park). Every location within each subarea had an equal chance of being sampled. Estimates forthe entire study area are made possible by accounting for unequal sample size among subareas. Estimates also can be made for individual subareas. The sampling design is not biased to favor one marsh type over another (e.g., selecting a specific location because it looks good or bad, sampling next to a road or airboat trail because it is easier, or avoiding tall, dense sawgrass or cattail (Figures 15 and 16) because it is difficult to sample). Probabilistic designs have two strengths: a) the results represent the spatial distributions of the environmental parameters that were measured; and, b) the results can be used to estimate, with known confidence, the proportion of the study area that was in any given condition, and how much the proportion changed. For example, the proportion of the Everglades marsh having a surface water sulfate concentration greater than 1.0 milligram per liter (mg/L, the CERP goal) was 37.1 ± 5.0 % during the 2014 sampling, and this proportion was statistically significantly smaller than the 57.3 ± 6.0 % measured during the 2005 sampling (Wald F test). This change indicates an improved condition in 2014. During 2014, 118 stations were sampled in the 2098 square-mile study area, and each station represents an average of about 18 square miles of marsh area. REMAP design disadvantages at this sampling density include the assumption of minimal heterogeneity throughout the 18-square mile area, and stations may not fall in localized areas with gradients near and downstream of water control structures, such as the phosphorus gradient in WCA2A. A transect sampling design along gradients, or a greater density of stations, are preferable in these areas. In the Phase IV 2014 sampling, USEPA utilized a design approximating a 50-50 mix of new random points and points from the previous phase (phase III, 2005). USEPA's Office of Research and Development, Western Ecology Division, National Health and Environmental Effects Research Laboratory, provided the statistical design and sample draw. The 2014 statistical design is a probability survey design that consists of two parts: a) 50% of the sites are a probability subsample of the prior survey design (a resample of 58 stations sampled in 2005); and b) 50% of the sites (60 stations) are a new probability sample. Since the 12 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) REMAP data provide a snap-shot of conditions throughout the Everglades during a two-week sampling window. REMAP detects change among survey events. Statements or inferences should not be made about conditions between surveys, especially for quickly-changing media such as surface water. two designs are completed independently, the combined survey design is also a probability survey design. The combined design has two objectives: estimate the current status of the Everglades across space; and detect change, or lack thereof, from prior time periods (2014 versus 1995 or 2005). An advantage of this combined approach is that the power of detecting a change is increased by visiting some sites in both time periods (Breidt and Fuller 1999, USEPA 2015). Throughout REMAP Program planning, a focus has been to assure that data meet key information needs of managers and scientists involved with Everglades protection and restoration. REMAP Program leaders met with Florida and Federal managers and scientists involved with CERP and Everglades phosphorus and mercury control efforts, including Everglades National Park, Arthur R. Marshall Loxahatchee National Wildlife Refuge, the Army Corps, the SFWMD and the Florida Department of Environmental Protection (FDEP). REMAP study plans have been subjected to external scientific peer review by these agencies and by the USEPA REMAP national program office. Efforts have been coordinated across agencies to avoid redundancy and maximize the information gained. Funding for REMAP phases has been provided by the Park, the Refuge, the Army Corps, FDEP, and USEPA. Elements such as carbon, nitrogen, phosphorus, oxygen, hydrogen and sulfur are critical components of living organisms. These elements cycle through ecosystems, plants, animals, air, water, and soil. These cycles are called biogeochemical cycles, because they include a variety of biological, geological, and chemical processes. A major focus of REMAP is sampling for these elements, and these data are referred to as biogeochemical data. The biogeochemical data are the focus of this report. REMAP provides a snapshot of conditions across the Everglades landscape during a two-week sampling window. Conclusions about conditions are based on the freshwater Everglades study area as a whole. REMAP detects change, or lack of change, by comparing survey events. Because the frequency of sampling is low (marsh sampling 5 years out of 20), statements about change should not be construed as traditional trend analysis. This is especially true for surface water constituents such as sulfate, specific conductance and nutrients, which can change quickly due to local rain events and water management operations. The concentrations observed during a two-week REMAP sampling snapshot may not represent conditions observed during other weeks in the same wet season. Consequently, statements or inferences should not be made about conditions during time periods that were not sampled. Long-term monitoring networks at fixed stations sampled more frequently, such as annually (the USGS Greater Everglades Priority Ecosystems Program) or weekly or monthly (the SFWMD's data that are reported by the SFWMD and FDEP in the annual South Florida Environmental Report), are better suited for making statements about surface water conditions at a location throughout the year, and identifying trends across consecutive years. 13 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) N A 0 2.5 5 10 Miles ill REMAP Stations ~ CANAL 1993-1995 O MARSH 1994 O MARSH 1995 O MARSH 1996 ® MARSH 1999 MARSH 2005 MARSH 2014 Canals Stormwater Treatment Areas Seminole Tribe - Big Cypress Federal Reservation Miccosukee Tribe - Everglades Federal Reservation Figure 14. All 1250 stations sampled by REMAP from 1993 to 2014. Canals -199 randomly located stations (4 sampling events 1993-1995), and marsh - 45 transects stations (1994) and 1006 randomly located stations (10 sampling events across 4 phases 1995-2014). See Table 1. 14 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Everglades REMAP History Everglades REMAP began sampling in the freshwater portion of the Everglades in 1993. REMAP was the first effort in the Everglades to sample canals at randomly located probability-based locations away from water control structures. Canals were sampled during the wet season in September 1993 and 1994, and during the dry season in May 1994 and 1995 (about 50 sites per sampling cycle) (Table 1, Figure 14, Stoberet al. 1995, 1998; Scheidt et al. 2000). Four marsh transects (45 stations) along phosphorus gradients downstream of water discharge structures were sampled during April 1994. Marshes were sampled at random locations in Phase I during the dry season (April 1995 and May 1996) and wet season (September 1995 and 1996), at about 120 sites per sampling cycle (Stober et al. 1998). The eastern portion of Big Cypress Swamp within Big Cypress National Preserve, the Seminole Tribe of Indians' Big Cypress Federal Reservation, and the Miccosukee Tribe of Indians of Florida's Federal Reservation in the Everglades were also sampled during Phase I. During Phase II, the freshwater Everglades marsh was sampled during May 1999 and September 1999 at another 119 sites per cycle (Stober et al. 2001a, 2001 b). Phase III was conducted in May 2005 and November 2005 at another 228 Everglades marsh sites (Scheidt and Kalla 2007). Phase IV was conducted Figure 15. Biogeochemical sampling included surface water (top left), floe and soil (bottom left), and mosquitofish (top right). The surface water sampling apparatus (top left) and soil coring device (bottom left) were designed and constructed for Everglades REMAP. Crews sample in all habitats, including dense sawgrass and cattail mix (bottom right). 15 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Phase I II III IV Year(s) 1993 - 1996 1999 2005 2013 - 2014 Distinguishing features Baseline data. Multiple stressors. Big Cypress & canals included. Change detection. Added periphyton assessment & plant studies. Omitted Big Cypress & canals. Change detection. Added food web studies & invasive plant survey. Change detection. Wet season only. Omitted aquatic community ecology & some macrophytic plant studies. Canal Stations 199 0 0 0 Marsh Stations 240 dry season 240 wet season 480 total 119 dry season 119 wet season 238 total 109 dry season 119 wet season 228 total 51 wet season 2013 118 wet season 2014 169 total Biogeochemical Media Surface water X X X X Floe X X X Porewater X X Bottom water X Soil X X X X Periphyton X X X X Mosquitofish X X X X Macrophytic Vegetation X X Macrophytic Plants Qualitative habitat categorization X X X X Species frequency X X Classified vegetation mapping X X X Invasive plant survey X Aquatic Community Ecology Periphyton assemblage X X Mosquitofish food habits X Macroin vertebrate assemblage X Isotope studies X Table 1. Everglades REMAP history showing phases, media and indicators. during September 2013 at 51 marsh sites (USEPA 2014a; Richards et al. 2017) and September of 2014 at 118 marsh sites (Kalla and Scheidt 2017). As of 2014, REMAP has sampled 1051 different marsh locations and 199 canal locations throughout the freshwater Everglades or Big Cypress, representing the ecological condition in 3,000 square miles of freshwater marsh and over 750 miles of canals. Phase IV Sampling In late September 2013, the USEPA Region 4 Science and Ecosystem Support Division (SESD) initiated a Phase IV wet season sampling event at 125 stations, and collected biogeochemical data at 51 stations 16 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) within ENP and WCA3 (USEPA 2014a, Richards et al. 2017). Due to a federal government shutdown during the sampling period, this event was not completed. Sampling forthe entire study area was reinitiated and completed in September 2014, and those results are presented here. The September 2014 Phase IV sampling included two efforts: WorldView-2 satellite imagery of a 1 km2 area centered on 65 REMAP sampling stations was used to produce classified vegetation community maps. Vegetation mapping provides a landscape context for REMAP biogeochemical and biotic information. Standing stocks of carbon, nitrogen and phosphorus in sawgrass, water, soil, floe, and periphyton were also estimated (Richards et al. 2017); multi-media biogeochemical sampling was conducted to understand water quality and soil conditions. This report focuses on the biogeochemical sampling. The biogeochemical effort included seven media (Table 1, Appendix I) that were sampled concurrently and consistently throughout the freshwater Everglades marsh. USEPA field crews sampled 118 stations during a 16-day period in September 2014 (Figure 15). Not all media were present at each station: surface water (n=116 stations); bottom water (the water at the bottom of the water column, immediately above the soil-water interface) (116); marsh soil or sediment (0 to 10 cm profile) (117); floe (flocculent detrital material found at the surface water-soil interface (96)); prey fish (mosquitofish) (104); composite periphyton (visible and floating within the water column plus epiphytic or 'sweaters' attached to plants) (71); benthic periphyton (discrete layer at the soil surface) (42). In addition, sawgrass, the most common plant in the Everglades, was sampled at a subset of stations: sawgrass leaf clippings (the middle 20 centimeters of a representative leaf from three sawgrass plants) (60); and whole sawgrass plants (27). These media are important for elucidating the cycling of nutrients and mercury. Unlike some long-term sampling programs, REMAP does not have a minimum surface watersampling depth due to the importance of shallow conditions in understanding cycling processes for mercury and nutrients. Digital photographs were taken to document conditions at each sampling location- ground view of the area sampled to the left of the helicopter, nine panoramic photos at 45-degree increments, each of the three soil cores, and an aerial view at 100 to 200 feet. For each station, photodocumentation, classified vegetation maps, and biogeochemical data are available at: http://diair.fiu.edu/amaps/EverREMAP.php. Data Quality Assurance Everglades REMAP has defined data quality objectives to ensure that data are of known and documented quality that satisfy predefined uses and requirements. An independently reviewed Quality Assurance Program Plan (QAPP) was developed in accordance with USEPA protocol (USEPA 2002, 2014b). Data quality is an essential component of the work, throughout planning, field sampling and laboratory analyses, and final data review and validation. USEPA Region 4 SESD Field Branch Standard Operating Procedures were followed as applicable (https://www.epa.aov/aualitv/aualitv-svstem-and-technical-procedures-lsasd-field-branches'). 17 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 16. The probability-based sampling desian ensures that Fi9ure 17- Surface water and Pore water samples in the chain of all habitats, such as dense cattail, are sampled. custody lab during 2005 at the end of a day's sampling. During 2014 samples were distributed to four analytical laboratories for determination of nutrient, ion and mercury content. During 2014 SESD was accredited for both field and laboratory operations by the American National Standards Institute - American Society for Quality National Accreditation Board. Some field sampling equipment and procedures have been developed specifically for Everglades REMAP (USEPA 2014a). During the September 2014 sampling, USEPA field personnel collected about 2140 samples (Figure 17). Four analytical laboratories performed the analyses for these samples. Data that potentially could be used for regulatory purposes, such as phosphorus, sulfur, and mercury, were obtained from analytical laboratories that are accredited by the National Environmental Laboratory Accreditation Program. Laboratory analytical methods and minimum detection limits are identified in the Project QAPP (USEPA 2014a). There were about 60 laboratory test methods performed on the samples, including forms of phosphorus, nitrogen, sulfur, carbon, mercury and physical parameters (Appendix I). Approximately 4,930 laboratory analytical sample results were produced, 100% of which were subjected to an independent quality assurance review. All analytical results met data quality objectives and none of the results were rejected. About 10% of the budget was invested in data quality assurance. Including field measurements, approximately 7,000 data values were generated. All biogeochemical and physical data are available at: https://www.epa.aov/everalades/environmental- monitorina-everalades. Data Uses REMAP data have been used for many purposes by environmental decision makers and scientists from over 30 Florida or federal agencies, Indian Tribes, non-governmental organizations, and agricultural and environmental interests. Program data have been used in over 40 peer-reviewed scientific journal 18 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) The REMAP Program consistently documents conditions throughout the Everglades in a quantitative manner. This documentation helps to assess the effectiveness of Everglades protection and restoration activities. publications or agency reports (see Program reports after Literature Cited section). REMAP Program reports or publications with REMAP data have been cited thousands of times. Key users include those that have provided funding, support, or sampling access: Everglades National Park, Arthur R. Marshall Loxahatchee National Wildlife Refuge, the Miccosukee Tribe of Indians of Florida, the Seminole Tribe of Indians, the U. S. Army Corps of Engineers, Florida Department of Environmental Protection and USEPA. The National Academies of Sciences also have relied on REMAP data in their scientific reviews of CERP and Everglades restoration progress (National Research Council (NRC) 2008, 2010, 2012, 2014). REMAP data can be used to help answer questions about multiple issues: Water management (e.g., water depth) Water quality and eutrophication (e.g., phosphorus concentrations in water, soil, vegetation and periphyton; cattail distribution) Habitat alteration (e.g., wet prairie and sawgrass marsh distribution, presence of exotic plant species) Mercury contamination (water, soil, algae, and prey fish) Specific Everglades restoration questions that REMAP helps to answer include: How much of the marsh has surface water sulfate concentrations that exceed 1 part per million (ppm), the CERP performance measure for Everglades marsh restoration? How much of the marsh has a soil total phosphorus concentration greater than 500 milligrams per kilogram (mg/kg), Florida's definition of phosphorus-impacted for Everglades soils, or 400 mg/kg, the CERP restoration target? How much of the marsh is dominated by sawgrass? Wet prairie? In what percent of the Everglades is cattail present? How much of the marsh has a natural oligotrophic periphyton community? How much of the marsh area is dry, and where? How much of the marsh soil has been lost due to subsidence? Has this loss rate changed overtime? How much of the marsh has prey fish with mercury levels that exceed 77 micrograms per kilogram, a level that presents an increased risk to top predators such as wading birds? What water quality conditions are associated with high mercury in fish? 19 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Everglades REMAP Data Uses 1- Support Everglades restoration efforts (used by many Florida and federal agencies, and in scientific publications): South Florida Water Management District (SFWMD); U. S. Army Corps of Engineers (USACE) 1996, 2014; USACE and U. S. Department of Interior (USDOI) 2015; RECOVER 2007, 2011; National Academies of Sciences (NRC) 2008, 2010, 2012, 2014; Wetzel et al. 2017. 2- Assess water management practices and drought-related ecological risk in the Everglades (SFWMD, Smith et al. 2003). 3- Document soil thickness and subsidence in the Everglades and identify areas that need to be restored with more water (USACE, SFWMD, National Park Service - Everglades National Park, Scheidt et al. 2000, Scheidt and Kalla 2007; NRC 2010, 2012, 2014; USACE 2014; Dreschel et al. 2018). 4- Model water quality implications of Everglades water restoration alternatives (Naja et al. 2017). 5- Determine which portions of the Everglades are phosphorus-impacted according to Florida's Everglades phosphorus criterion rule and assessing the effectiveness of Florida's multi-billion dollar phosphorus control effort (Florida Department of Environmental Protection (FDEP), SFWMD, United States Fish and Wldlife Service (USFWS) - Arthur R. Marshall Loxahatchee National Wldlife Refuge, Everglades National Park, USEPA; (Payne et al. 2001, 2002) 6- Determine landscape patterns of soil nutrients and contaminants (USACE, USEPA, SFWMD, Everglades National Park, Loxahatchee National Wldlife Refuge; Osborne et al. 2011a; NRC 2008,2012). 7- Understand morphological response of plant species to phosphorus (Richards and Ivey 2004). 8- Determine the effects of phosphorus on periphyton, invertebrates and fish (Sargeant et al. 2010, 2011). 9- Understand spatial variations in aquatic food webs and food web structure and drivers (Abbey-Lee et al. 2013). 10- Map vegetation and determine vegetation classes, associations, biomass, nutrient standing stocks, and presence of exotic plants (Richards et al. 2008, 2017; Zweig and Newman, 2015). 20 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 11- Understand sulfur cycling, the sources and distribution of surface water sulfate, the penetration of water with high sulfate content into the Everglades marsh, and whether there should be water quality standards for sulfate (FDEP, SFWMD, USEPA; SFWMD 2007, Axelrad et al. 2007, 2008; NRC 2010, 2012; Corrales et al. 2011, Gabriel et al. 2014a, Orem et al. 2011, Julian et al. 2016b, Orem et al. 2019). 12- Understand the penetration of water with high ionic content from agricultural areas into the low ionic content marsh in the Refuge and its potential impacts on natural periphyton communities (USFWS - Loxahatchee National Wildlife Refuge; McCormick and Harvey 2007; NRC, 2010, 2012; McCormick et al. 2011). 13- Document water quality background conditions and develop performance measures for soil phosphorus and surface water conductivity and sulfate for the multi-billion dollar Florida-Federal Comprehensive Everglades Restoration Plan (USACE, SFWMD, Everglades National Park, Loxahatchee National Wldlife Refuge; RECOVER 2007). 14- Determine landscape patterns of periphyton, and the role of periphyton mats in aquatic consumer community structure (USACE, SFWMD, Everglades National Park; Loxahatchee National Wldlife Refuge; Gaiser et al. 2011; Trexler et al. 2014). 15- Document mercury conditions, landscape patterns, and mass balance in soil, water, plants and biota (Liu et al. 2008a, 2008b, 2015; Richards et al. 2017). 16- Document mercury biomagnification from water or sediment to fish and understand the ecological risk to Everglades birds (FDEP, SFWMD, USEPA, USACE, Loxahatchee National Wldlife Refuge, Everglades National Park; PTI 1995; USACE 1996; USEPA 1997; Axelrad et al. 2007; Julian 2013a, 2013b). 17- Elucidate the environmental conditions, such as sulfur, carbon and phosphorus, that may influence mercury in prey fish (FDEP, SFWMD, USACE, USEPA, Loxahatchee National Wldlife Refuge, Everglades National Park; Liu et al. 2009, 2011; Julian et al. 2013a, 2013b, 2016b, 2017a; Dierberg et al. 2014; Pollman 2012, 2014, 2019; Pollman and Axelrad 2014; Jiang et al. 2018, Rumbold 2019a; Kalla et al. 2019) 18- Evaluate the distribution, potential sources and controlling factors of toxic metals (Li et al. 2015). 19- Determine carbon and organic matter characteristics and drivers across the Everglades landscape (Yamashita et al. 2010). REMAP data have been presented by Program personnel in over 40 peer-reviewed scientific journal publications or agency reports. 21 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Data Analysis and Presentation Kriging A strength of Everglades REMAP is its spatial coverage. Kriging is a geostatistical method of generating contour maps to understand and display data across space (Figure 21). The contours are isopleths, lines of equal estimated value for the measurement. Kriging algorithms interpolate between actual data points, producing a grid of estimated values from which the contours are drawn. The krigs in this report are true to the data - i.e., the data value at each sampling station matches the color of the contour interval at that point. Sampling stations are shown on krigs with black dots. Krigs were made by estimating a value for each node (intersection of lines) of the grid using the linear variogram model (no nugget effect). A variogram is an expression of how quickly the actual values change over space, on average, while taking into account the overall variability of the data set. An underlying assumption of variograms is that values from points closer together are more similar than those from points farther apart. Fifty-eight of the stations sampled in 2014 were also sampled in 2005. For the soil parameter krigs, data from both years were used. Krigs are included in this report to provide visual information across space. Conclusions about the spatial extent of conditions, or that conditions were different during sampling events, are made not from the krigs, but rather from rigorous statistics (the cumulative distribution function and the Wald F-test, see below). In the Everglades, there are physical barriers to water sheet flow, such as levees and roads. However, during the peak of the wet season (generally the time of wet-season REMAP samplings reported here), the subareas of the Everglades are hydrologically connected by surface water flowing through numerous water control structures. The 2014 sampling took place from September 4-20. During this time, and the month before, all of the Stormwater Treatment Areas were discharging into the Everglades, and there was water moving across the subareas: from the Refuge to WCA2A, WCA2A to WCA3A, and from WCA3A to the Park at Shark River Slough (Abtew and Ciuca 2016). The same was true of the November 2005 sampling, except that there was no flow from the Refuge into WCA2A. While wet season connectivity is not complete, isolation of the subareas is not complete either. Reality lies somewhere in between, and is dependent upon location, ground relief, water management operations, and proximity to canals, levees and water control structures. Some water quality constituents, such as surface water total mercury, are driven more by atmospheric deposition than by water flow, minimizing the influence of physical barriers such as levees. For these parameters, kriging was performed on the entire study area as one unit. Other constituents, such as surface water chloride, conductivity, sulfate and pH, are influenced by surface water inflows and proximity to canals and levees, especially in the vicinity of the Refuge and WCA2. For these constituents krigs were done by subarea. Spearman rank order correlations Correlation is a statistical tool for determining the strength of association between two variables. Parametric statistical tests assume that the data values are random, and independent. In the Everglades, values for some constituents, such as sulfate in surface water, occur mostly nearthe low end of their range, while some large values also occur. These data are not normally distributed, but instead are skewed. In addition, the 22 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) sulfate concentration at one location tends to be similar to the concentrations observed at nearby locations, so these values are not independent (see Figures 40-41). Nonparametric statistics, such as the Spearman rank order correlation, are preferable for such data. If two variables are perfectly correlated, the change in one variable is accompanied by equivalent change in the other variable. The correlation coefficient can be any number between -1 and +1. Positive values indicate direct correlation, where one variable increases as the other increases, whereas negative values indicate inverse correlation (one variable decreases as the other variable increases). Coefficients near 0 indicate lack of correlation. REMAP measurements vary over space (station to station). Data from the same sampling event can be analyzed for correlation because the measurements at a station were obtained at the same time. For the 2014 sampling, correlation coefficients and their statistical significance have been reported (Kalla and Scheidt 2017). The significance (p) of a statistic is a measure of the reliability of the sample data set as representative of the entire population of possible data points. The value of p is the chance that the true coefficient in the population is 0, or in other words, that instead of a strong or even a weak correlation, there is none at all. A correlation coefficient of p<0.001 indicates a 1 in 1000 chance that the two variables are not associated. For example, the correlation coefficient between surface water sulfate and conductivity was 0.71, which was significant with a probability < 0.001 (Kalla and Scheidt 2017). Correlation or association between two variables does not necessarily indicate that there is an underlying explanation or cause. Results are highlighted here only if p<0.001 and if there is a plausible mechanism or process to explain why the two variables would be correlated. CDFs and area estimates One way to visually portray survey statistics is to plot the cumulative distribution function (CDF) of the data. A CDF curve can be used to estimate the area (proportion) of the Everglades where a given analyte was found at a concentration above or below any value of interest. This is a key strength of REMAP's probability-based sample design. The estimates of area by CDF curve in this report were generated from the original data, using algorithms and scripts for the R statistical package developed by EMAP statisticians at the USEPA Office of Research and Development. Krigs were not used to estimate CDF curves. In this report, the CDF curve is shown in bold. By reading up to the CDF curve from any concentration of interest on the x-axis, and then across from the curve horizontally to the y-axis, one can read the corresponding proportion of the Everglades on that axis. Bounding the CDF are two lines representing the upper and lower 95% confidence limits, respectively, calculated using the Horvitz-Thompson estimator. These limits show the confidence interval (CI) around the area estimate. This interval, expressed as percent values above and below the estimate, indicates the precision of the estimate: narrower intervals represent more precise estimates. Estimates tend to become more precise as the number of samples increases. At the 95% confidence level, there is a 1 in 20 chance that the true value for the Everglades study area was outside the range defined by the confidence interval. For example, looking at the CDF for mosquitofish mercury in 2014 (Figure 71), in 2014 87 ± 6% of the 2,098-square-mile Everglades study area had mercury concentrations in mosquitofish below 77 |jg/kg, a threshold to protect birds and mammals that feed on fish (USEPA 1997). The 95% CI around the 87% estimate 23 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Probability-based sampling design is an assessment approach that provides unbiased estimates of ecosystem condition with known confidence. is ± 6%, which is within the ± 10% REMAP data quality objective tolerance limits for 95% CIs. Previous experience in the national EMAP Program and in the earlier phases of Everglades REMAP (Scheidt and Kalla 2007) indicated that approximately 125 stations was a sufficient sample size to meet this objective. Statistical testing for differences across sampling year (F test) Any pair of REMAP data sets for a given variable, represented by their respective CDF curves, can be tested statistically to indicate a difference, or lack thereof, between them. The Wald F test was used to test 2014, 2005, and 1995 CDF curves against each other to determine whether there was a significant change in conditions among sampling events. This test allows for statistical inference to the sampled population, or study area. If curves from different time periods differ, the underlying condition of the resource can be said to have changed, although nothing can be inferred about the intervening years. Statements about change are made with a specified degree of confidence, typically no more than a 1 in 20 chance of being wrong (probability < 0.05). The source of such an error is that the supposed difference is due merely to samples that inaccurately represented the population, instead of being a real change in the resource. The random REMAP sample spread out over the entire Everglades study area can be used to draw conclusions about changed conditions across the whole area. In this report, the largest p-value from the F test used to support a conclusion that there was a changed condition is 0.03. In addition, a weight of evidence approach was used. While comparing krigs may suggest changed conditions, this assertion was corroborated by not only an F test with p<0.05 (Wasserstein et al. 2019), but also by lack of overlap of the CDF curve 95% confidence intervals across years, lack of overlap of the boxes in the box and whisker plots, and a plausible explanation for a change in condition. Box and whisker plots A box and whisker plot (for example, Figures 31 or 36) is a graphical method of displaying the variation, or general shape, of a data set. The large box contains the middle 50% (the interquartile range between the 25th and 75th percentiles) of the data values. The median is shown by the square symbol or horizontal line within the box. Half of the values are greaterthan the median, or 50th percentile, and half are less. The whiskers or vertical lines include values outside the interquartile range that are not considered outliers or extremes. Outliers are defined as values that are smaller than the 25th percentile or larger than the 75th percentile, respectively, by at least 1.5 times the interquartile range. Extremes are values that are smallerthan the 25th percentile or largerthan the 75th percentile, respectively, by at least 2 times the interquartile range. 24 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) WATER CONDITIONS I I Figure 18. Rainfall at an Everglades wet prairie-slough in WCA3 during 2005. Surface water depth is about 2 feet. Water is the lifeblood of the Everglades. The historic Everglades was defined by: highly seasonal rainfall; shallow, unimpeded, sheet-like surface flow; and a large storage capacity that prolonged wetland flooding. These characteristics, along with subtle changes in ground surface elevation of only a few feet over tens of miles, produced a variety of water depths and hydroperiods (durations of surface water inundation). Slowly- flowing water helped create the ridge and slough landscape that is a defining feature of the Everglades. Water depth and hydroperiod affect fish populations, as well as wading bird feeding habitat, prey availability, and nesting success. During the 1900s, anthropogenic changes in surface water depth, distribution and hydroperiod caused many harmful impacts to the Everglades. Water is key to ecosystem preservation and restoration. South Florida water management operations are complex and involve balancing the water supply and flood control needs of urban areas, agriculture, and the environment. Canals can move Lake Okeechobee water east to the St. Lucie Estuary, west to the Caloosahatchee Estuary, or south through the Everglades Agricultural Area (EAA) into the Everglades. These southern canals then pass through the Everglades marsh on their way to outlets along the east coast (Figures 10 and 14). Water in these canals can also flow through water control structures into the Everglades marsh. In the EAA, the canals are used for irrigation, drainage and flood control, depending on the season and local rainfall. In drier years, less water flows from Lake Okeechobee and the EAA downstream into the Everglades. Each of the Everglades WCAs and Lake Okeechobee have regulation schedules that are used by the Army Corps and SFWMD to manage water levels and releases. The regulation schedules have water level thresholds that vary with time of year, and trigger water releases that are made to protect the integrity of 25 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Everglades Region Annual Rainfall (inches) Average WY 1991-2020 WY 1996 WY 2006 WY 2015 51.6 60.8 48.1 47.7 3 Month Inflow water volume (1000 acre-feet) September 1995 November 2005 September 2014 LNWR 432.6 90.9 124.1 WCA2 542.6 265.2 422.4 WCA3 1,068.3 487.3 628.9 Park 659.4 802.4 498.1 Table 2. Annual rainfall for REMAP sampling years and period of record, and three-month discharge into the Refuge, WCA2, WCA3 and the Park. .E -3- <0 a: WY1991-2020 m WY1996 WY2006 WY2015 Figure 19. Everglades region annual rainfall (inches) for water years 1996, 2006 and 2015 as compared to the WY1991 - WY 2020 annual average. 7000 WY1996 Figure 20. Annual flow (thousand acre-feet) into the Refuge, WCA2, WCA3 and the Park during water years 1996, 2006 and 2015. levees and developed areas, and to lower water levels in preparation for wet season rainfall and inflows. The regulation schedules are developed through public processes, and are designed to balance competing objectives including water supply, flood control, and environmental enhancement. Rainfall in the Everglades is highly seasonal, with about 80% falling during the May to October wet season (Figure 18). Climatic events such as tropical storms, droughts and El Nino all have a profound effect on rainfall. In South Florida, a wateryear (e.g., WY15 is May 1, 2014 to April 30, 2015) is used to summarize rainfall and hydrologic data. For the years with wet season REMAP sampling events, annual rainfall throughout the Everglades was comparable in WY06 (48.1 inches, 93% of WY91-WY20 annual average) and WY15 (47.7 inches, 92%), and highest in WY96 (60.8 inches, 118%) (Figure 19, Table 2). Water flow into and within the Everglades is highly seasonal and varies annually and by location. During the month of REMAP sampling and the prior two months, flows at all structures into the Refuge, WCA2 and WCA3 were highest in 1995, lowest in 2005 and intermediate in 2014. In contrast, for the Park inflow was highest in 2005 (Table 2, Abtew et al. 2007, Abtew and Ciuca 2016). Annual inflows into the Refuge, WCA2, WCA3, and the Park were highest in WY96, lower in WY06 and lowest in WY15 (Figure 20). Inflow does not dictate marsh water depth since outflow may balance inflow. The Stormwater Treatment Areas (STAs) discharged 1.32 million acre-feet of treated water into the Everglades WCAs from the north in WY15 (Pietro et al. 2016), as compared to 1.48 million acre-feet in WY06 (Pietro et al. 2007). During the September 2014 sampling and the month before, all of the STAs were discharging, and there was water moving from the Refuge to WCA2A, 26 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 21. Krigs of measured surface water depth encountered in the Everglades marsh during the 1995, 2005 and 2014 REMAP wet season sampling events. Figure 22, Wet season average surface water depth, September to November, 1991-2014. From Everglades Depth Estimation Network (EDEN), courtesy of USGS. Figure 23. Dry marsh in Everglades National Park. The Eastern Marl Marsh within the Park is dry on average about 40% of the year. 27 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) WCA2Ato WCA3A, and from WCA3Ato the Park at Shark River Slough (Abtew and Ciuca 2016). The same water movement was true for the November 2005 sampling, except that there was no flow from the Refuge to WCA2A (Abtew et al. 2007). Marsh water depths vary greatly with season, year and location in response to short-term and long-term rainfall and flow from water control structures (Figures 18, 21-23). During all REMAP wet season sampling events the entire marsh was inundated. Water depths are deepest immediately upstream of levees that impede the natural flow of water, such as in the southern portions of the Refuge, WCA2 and WCA3A (Figures 21-22). Unnaturally deep water (over five feet) was observed in 1995 within eastern WCA3A where the L67 levee prevents water sheet flow to the south. Short hydroperiod portions of the marsh dry out in most years, such as the marl marshes within the Park (Eastern Marl Marsh and Ochopee Marl, see soil section) and northern WCA3A (yellow areas in Figure 21). Figure 22 shows Everglades Depth Estimation Network (EDEN; Haider et al. 2020) model projections by USGS of average wet season water depth for September to November, 1991-2014. The overall result of rainfall and water management operations was that for most of the Everglades water depths during the 1995 sampling were deeper than a typical wet season. This was the year with the highest rainfall, and the highest water flow into the WCAs. For the Park and WCA3A the shallowest conditions were observed in 2014, the year with the lowest flow into the Park at Shark River Slough, which may in part reflect new Everglades Restoration Transition Plan water management operations that were initiated in 2012. WATER QUALITY Nutrients Nutrients are elements that are essential for plant growth. Phosphorus (P), nitrogen (N), and potassium are considered the primary nutrients, and they are the most common nutrients in fertilizer. The Everglades marsh is phosphorus-limited (McCormick et al. 1999; McCormick et al. 2002), which means that very low concentrations of P naturally limit undesirable biological growth. The consequence of P limitation is that elevated P concentrations cause undesirable ecological responses, such as changes in algal and plant communities, and loss of dissolved oxygen in water. Much of the attention about Everglades water quality has been focused on controlling P in order stop its deleterious ecological effects. Although N enrichment has not been identified as a concern in the freshwater Everglades, it is a concern in coastal waterbodies including Florida Bay and the Caloosahatchee and St. Lucie Estuaries. Florida has adopted protective numeric nutrient criteria (total phosphorus, total nitrogen and chlorophyll a, Chapter 62-302.532(d) and (z) F. A. C.) for the St. Lucie and Caloosahatchee Estuaries, and has identified the loads of pollutants such as nutrients that are required to meet water quality standards in Lake Okeechobee, and the St. Lucie and Caloosahatchee Estuaries. 28 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Phosphorus The natural Everglades marsh has very low total phosphorus (TP) concentrations of less than 10 |jg/L (micrograms per liter, or parts per billion (ppb)) in surface water (McCormick et al. 1999, Noe et al. 2001). Concentrations can be as low as the laboratory analytical method detection limit of 2 ppb (Julian et al. 2016d). Historically, rainfall was the dominant source of external phosphorus, and the hydrology of the marsh was rainfall-driven, with slow, overland, shallow sheet flow supplying water to downstream wetlands. There were no canals in the Everglades priorto the early 1900s. This naturally nutrient-poor condition resulted in a unique mosaic of wetland habitats, such as sawgrass marshes, wet prairies and sloughs, and algal communities (periphyton) that developed under low TP conditions. Well-developed periphyton communities are a defining characteristic of the Everglades, and they play a critical role in food webs, habitat, soil formation and water quality (McCormick 2011). The first canals that dissected the Everglades were completed by 1917. From 1954-59 levees forming the boundary of the EAAwere constructed, along with pump stations to provide flood control. From 1960-63 levees were completed that compartmentalized the Everglades into WCAs, along with structures allowing flow into and out of the WCAs (Light and Dineen 1994). The canal system became an inadvertent conduit for transport of high levels of nutrients in water pumped out of the EAA into the Everglades. Water quality became a concern (Cornwell et al. 1970, Marshall 1971, Storch 1971, US Department of Interior 1971), and by the 1970s scientists began to observe phosphorus impacts in Everglades marshes, especially in the Refuge and WCA2A (Dineen 1973, Gleason 1975a, 1975b, SFWMD 1977). From 1979-88 the flow-weighted TP concentrations discharged into the Refuge from the north and west were 198 ppb and 144 ppb at pump stations S5Aand S6 respectively, and they were 129 ppb at the S10 structures that discharge from the Refuge into WCA2A, as compared to marsh interior concentrations of about 10 ppb or less (SFWMD 1992a). At this time there were no numeric phosphorus requirements for water within or flowing into the Everglades. In the Everglades, annual or long-term means for phosphorus at water control structures are typically expressed as a flow-weighted mean. When calculating a flow-weighted mean, the TP concentration data are weighted by the flow during the time period that the concentration represents in order to provide an accurate estimate of average conditions. During the 1980s, TP impacts in the Everglades became well-documented (Swift 1981, Dineen 1986, Frydenborg and Ross 1987, Swift and Nicholas 1987, LOTAC II 1988, Belangeret al. 1989, Scheidt et al. 1989, SFWMD 1992a, Davis 1994). Phosphorus impacts are most evident downstream of canal and water structure inflows to the Refuge, WCA2A, the Miccosukee Tribe of Indians' Everglades Federal Reservation in WCA3A, and the Park (Doren et al.1996, Childers et al. 2003, Wright et al. 2009, Osborne et al. 2013). Among the undesirable, cascading eutrophic impacts to the Everglades are altered periphyton communities, loss of water column dissolved oxygen, increased soil phosphorus content, conversion of the wet prairie- sawgrass mosaic to tall, dense single-species stands of cattail with minimal open water, and subsequent loss 29 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) of wading bird foraging habitat. Vegetation impacts can be observed from satellite imagery. These collective changes impact the structure and function of the aquatic ecosystem (McCormick et al. 1999, Payne et al. 2001, McCormick et al. 2002, Gaiser et al. 2005, McCormick et al. 2009). By about 1990, over 40,000 acres of the Everglades were impacted by cattail expansion (Davis 1994). How clean the water should be, how to control phosphorus, and who should pay for these efforts became the subject of federal litigation in 1988 (Godfrey and Catton 2006, Rizzardi 2001a, 2001b). This resulted in the first requirements for TP control programs (Consent Decree (CD), U. S. District Court 1992), which became incorporated into Florida's 1994 Everglades Forever Act (EFA). Today, all of the Everglades has a numeric water quality criterion of 10 ppb for phosphorus. In 1999 the Miccosukee Tribe of Indians of Florida adopted, and under the Clean Water Act USEPA approved, a 10 ppb TP criterion forthe Everglades portion of the Tribe's Federal Reservation in WCA3A (Figure 14). In 2005 Florida adopted and USEPA approved a 10 ppb water quality criterion for TP in the Everglades Protection Area (EPA), which includes the Refuge, Park and WCAs (FAC 62-302.540). The objective of both of these water quality criteria is to prevent nutrient-induced imbalances in natural populations of aquatic flora orfauna. Florida's TP criterion is applied in the EPA as a long-term average, with achievement of the criterion within the Everglades waterbody determined by data collected monthly at a network of fixed long-term marsh sampling locations. In the Everglades, geometric means are routinely used to average marsh phosphorus data to minimize the influence of extreme high values on the average. Compliance is determined by a 4-parttest specifying that: 1) the five-year geometric mean averaged across all stations is < 10 ppb; 2) the annual geometric mean averaged across all stations is < 10 ppb for three of five years; 3) the annual geometric mean averaged across all stations is < 11 ppb; and 4) the annual geometric mean at each individual station is < 15 ppb. Each of the four parts must be met to achieve the criterion (FAC 62-302.540(4)). The test is intended to simultaneously: allow forthe natural temporal and spatial variability that is observed at marsh reference sites; be sensitive enough to detect long-term increases in TP above 10 ppb; and place an upper limit on phosphorus at individual marsh locations. The same test is applied separately at 'impacted' and 'unimpacted' stations, with 'impacted' areas defined as those where the total phosphorus concentration in the upper 10 centimeters of the soil isgreaterthan 500 rng/kg. Forthe Park, achievement of the criterion is determined not by the 4-part test at marsh stations within the Park, but rather by P Figure 24. Location of phosphorus control program treatment wetlands (Stormwater Treatment Areas, or STAs) operational in 2017. In combination with agricultural best management practices, they are to discharge phosphorus into the Everglades at a long- term concentration equivalent to the 10 part per billion water quality standard. A-1 is a shallow Flow Equalization Basin that equalizes water flow into the STAs (adapted from SFWMD). 30 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) concentration requirements at Park inflow structures ('at the door')(U. S. District Court 1992, SFWMD 1992b). Because REMAP data are not collected monthly throughout the year at fixed sampling locations, it is not possible or appropriate to apply Florida's 4-part test for assessing the 10 ppb criterion to REMAP marsh data. However, because of REMAP's probability-based design, it is appropriate to make statements about the area of the marsh that exceeds a concentration of interest, such as 10 ppb, during individual REMAP sampling events, although these statements have no regulatory bearing. The Park and Refuge have an additional level of water quality protection because they have been designated by Florida as Outstanding Florida Waters (OFWs). This antidegradation designation requires that the quality of water that existed the year prior to March 1,1979 must be maintained. To protect these OFWs, the CD requires a long-term average geometric mean TP concentration of 7 ppb at a network of 14 interior marsh stations in the Refuge, and a long-term flow-weighted mean (FWM) of 8 ppb at inflows to the Park at Shark River Slough, and 6 ppb at inflows to the Park at Taylor Slough, respectively (SFWMD 1992b, U. S. District Court 1992, Walker 2000). These Park requirements were also adopted in the Florida Administrative Code. In addition, CERP has a performance measure for surface water TP: TP concentration is not to exceed 10 ppb for both the annual geometric mean at marsh stations and the flow-weighted annual geometric mean at water control structures, and should not exceed OFW concentration levels (RECOVER 2007). Phosphorus control programs, required by the CD and Florida's EFA, were initiated in the 1990s in order to prevent further loss of Everglades plant communities and wildlife habitat. The initial phase of this program required that TP in discharges from the EAA into the Everglades be at 50 ppb or less. Control was achieved by a combination of agricultural Best Management Practices (BMPs) and about 47,000 initial acres of constructed treatment wetlands within the EAA (the Everglades Construction Project, ECP) referred to as Stormwater Treatment Areas (STAs) (Figure 24). Agricultural BMPs were required to be in place by 1995. The first prototype STA (the Everglades Nutrient Removal Project, 3700 acres) began treating water in 1994 (Chimney and Goforth, 2006). The REMAP Phase I sampling period (1993-96) corresponds to the phase-in period for EAA BMPs, as during these years the percentage of EAA farms with phosphorus control BMPs in place went from 0 to 100. Full BMP implementation began in 1996 with a 25% TP load reduction required at the EAA basin level. # 2.0 to g.g # 10.0 to 29.9 30.0 to 49.9 # 50.0 to 99.9 # 100.0 to 467.0 Figure 25. Total phosphorus (pg/L) in surface water at 99 random locations in canals sampled by REMAP during the wet season, September 1993 and 1994 (Scheidt et al. 2000). At this time there were no phosphorus controls. Phosphorus conditions have improved since the 1990s due to the STA and EAA agricultural BMP programs. 31 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) During the 1990s, REMAP documented that canal TP concentrations had strong north to south gradients due to stormwater pumping, with the highest TP concentrations in canals in the EAA during the wet season (median of 149 ppb). During the 1993-94 wet season, about 80% of the canal miles in the EAA had TP concentrations greater than the initial TP control target of 50 ppb, and overall, 44% of Everglades canal miles had water TP concentrations greater than 50 ppb (Figure 25, Scheidt et al. 2000). These data represent conditions before TP controls were initiated. TP concentrations in surface water during September 2014 are shown in Figure 26. About 70% of the Everglades marsh stations sampled by SFWMD had TP < 10 ppb, and 35% had TP < 5 ppb. The highest marsh TP of 42 ppb was in WCA2A near canal inflows. For the 28 days prior to REMAP sampling, STA inflows were 61-142 ppb FWM (provisional data) and discharges were 13-21 ppb, except for STA5/6 which had an inflow of 240 ppb and a discharge of 28 ppb (SFWMD 2014). For WY15, annual geometric mean TP concentrations at marsh stations within the Park, away from canals and anthropogenic nutrient sources, were as low as the laboratory analytical method detection limit of 2 ppb (Julian et al. 2016d). In contrast, TP concentrations in water flowing into the Miccosukee Tribe of Indians of Florida's Federal Reservation in WCA3A were 76 and 78 ppb, greater than the Tribe's 10 ppb water quality criterion (Figure 26). SFWMD's extensive long-term data, as frequent as bi-weekly at fixed locations, show improving TP conditions in the Everglades (Julian et al. 2016a). SFWMD average annual TP data at the inflows to the Canals Total Phosphorus ~~ Stormwater Treatment Areas (M9'L) O S10 Miccosukee Tribe Federal Reservation O >10-50 10 20 Miles ~ USEPA REMAP Study Area q >5o.-|oO O >100 J i i i I Figure 26. Total phosphorus (|jg/L) in surface water during September 2014 at 116 REMAP Everglades marsh stations. Also shown are SFWMD data for 12 locations in Lake Okeechobee and about 100 locations in the Everglades marsh (geometric mean), and 33 water control structures in canals (geometric mean when flowing) including selected STA inflows and outflows. 32 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) foursubareas ranged from 12-21 ppb for WY15 (2014), as compared to 10- 67 ppb during WY06 (2005, Table 3). Interior marsh annual geometric mean TP concentrations also improved from 6-18 ppb in 2005 to 4-9 ppb in 2014. REMAP data also show progressively lower surface water phosphorus concentrations from 1995 to 2005 to 2014 (depicted in Figure 27 and evaluated in Figure 28). For September 2014,14.8 ± 4.8% of the marsh had a TP concentration higher than 10 ppb, an improvement from the 27.2 ± 7.5% observed during November 2005 (Figure 28). Conditions were worse during September 1995 when 57.8 ± 7.8% of the Everglades marsh had TP greater than 10 ppb. The proportion of the marsh above 10 ppb was cut in half twice for these three years, and the improvements indicated by the CDF curves are significant (Wald F, p=0.04). INFLOW OUTFLOW INTERIOR Flow-weighted mean Geometric mean 2005 2014 2005 2014 2005 2014 EAA BMP Basins 117-221 62-144 STAs 144 (104-213) 99 (72-230) 44 (21-146) 17 (15-32) Refuge 67 20 15 9 WCA2 27 15 18 9 WCA3 24 21 10 5 Park 10 12 6 4 Table 3. Annual surface watertotal phosphorus concentration (ppb) from selected Everglades locations for WY06 (2005) and WY15 (2014) (data from SFWMD). TOTAL PHOSPHORUS SURFACE WATER NOVEMBER 2005 Figure 27. REMAP total phosphorus concentration in surface water (|jg/L) during November 2005 (left) and September 2014 (right). 33 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 28. Surface water total phosphorus concentration estimates of marsh area for 1995, 2005 and 2014 wet season REMAP samples. The areas greater than 10 ppb were progressively smaller, reflecting improved conditions (Wald F test). The improvement in Everglades water phosphorus concentrations since 1995 is due to the EAA BMP program (Izuno et al. 1999, Daroub et al. 2009, 2011, Lang et al. 2010, Yoder et al. 2020) and Florida's STA program. From 1996 to 2015, it is estimated that the BMP program removed 3001 metric tons of TP (a 56% TP load reduction) from the EAA basin as compared to the load that would have been expected without BMPs. Post-BMP EAA TP concentrations for WY15 were 93 ppb basin-wide, with a 79% load reduction (SFWMD 2016a, ludicello et al. 2016). In addition, the area of STAs treating post-BMP water increased by 73% from 32,980 acres during the REMAP 2005 sampling to 57,045 acres during the 2014 sampling. This is the largest wetland treatment system in the world. For WY15, flow-weighted annual mean TP inflow to the STAs ranged from 72 ppb for STA2 to 230 ppb for STA5/6, with an average inflow for all STAs of 99 ppb. STA discharge concentrations into the Everglades ranged from 15 ppb for STA3/4 to 32 ppb for STA5/6, with an average discharge across all STAs of 17 ppb, as compared to 44 ppb for WY06 and the previous REMAP sampling (Table 3). The overall WY15 load reduction for the STAs was 83%. The cumulative amount of phosphorus retained by the STAs from 1994 to 2015 was about 2000 metric tons, which is 75% of the inflow load (Pietro et al. 2016). This is about two-thirds of the amount removed by agricultural BMPs. It is easier to remove P near the source where concentrations tend to be higher. In WY06, the flow-weighted TP concentration discharged into the Refuge, WCA2 and WCA3 was 64 ppb (Payne et al. 2007). In WY15 it improved to 19 ppb (Xue 2016). Some of the cells orflowways within the better performing STAs, STA2 and STA3/4, have discharged at 13 ppb for several years. Complex factors affect STA TP removal and outflow TP, including calcium (Juston and DeBusk2011), annual hydraulic and phosphorus loading rates, hydraulic residence time, inflow TP concentration, water depth, vegetation coverage and condition, soil characteristics, antecedent land use, and extreme weather (Chen et al. 2015, Zhao and Piccone 2020). While there is considerable 34 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Water phosphorus concentrations were much lower in 2014 than in 2005 and 1995-96 due to agricultural BMPs and the Stormwater Treatment Areas effort underway to reduce future phosphorus loads into the Everglades, areas of the Everglades that are impacted by P often are dominated by cattails, which are resilient (Bansal et al. 2019). Active management in these P-impacted areas may speed their recovery (Hagerthey et al. 2014). During 2012, USEPAand Florida reached consensus on additional Everglades water quality Restoration Strategies (SFWMD 2012). STA discharge permits and consent orders issued by FDEP require: a) TP at the discharge from each STAshall not exceed 13 ppb as an annual flow-weighted mean (equivalent to 10 ppb as a geometric mean at the STA discharge) in more than 3 out of 5 water years on a rolling basis, and shall not exceed 19 ppb as an annual flow-weighted mean in any water year; b) expansions to STA1W totaling 6,500 acres; c) about 110,000 acre-feet of water storage areas (flow equalization basins) that will slowly release water to the STAs in order to improve their performance; d) an enforceable compliance schedule for $880 million of projects with completion dates of 2018 to 2025; and e) a monitoring and research plan to confirm that the performance of each STAis optimized and restoration is moving forward. The Regional Administrator for USEPA Region 4 and the Secretary of the FDEP meet bi-annually to ensure that progress is being made until each STA meets the 13 ppb long-term phosphorus discharge limit. Technical representatives from state and federal agencies also meet at least twice annually to review research, evaluate operation of the STAs, and assess water quality and progress in achieving the deadlines. Florida has also developed a comprehensive long-term plan for achieving water quality goals forthe other basins (non-ECP) that discharge into the Everglades (Burns and McDonnell 2003). This plan to treat large volumes of stormwater down to 10 ppb TP is unprecedented. The plan recognizes that additional control measures are necessary to ensure that discharges to the Everglades from all basins meet water quality standards (ludicello et al. 2016), including the basins flowing into western WCA3A. Soil Phosphorus Excess phosphorus in marsh soils is an aggregate indicator of TP enrichment over a longer time scale than TP in the water column. In portions of the Everglades with peat soils there is an association between increasing soil TP and cattail encroachment. Various elevated soil TP concentrations have been used as indicators of enrichment or where cattails occur: 700 milligrams TP per kilogram of soil (mg/kg) (Childers et al. 2003); 610 mg/kg (Walker and Kadlec 1996); and 600 mg/kg (Craft and Richardson 1993, Payne et al. 2001). Florida's Everglades TP criterion rule defines P-impacted areas as being where soil TP exceeds 500 mg/kg (FAC 62-302.540(3)(d)). CERP has a restoration goal of decreasing the areal extent of the Everglades with soil TP > 500 mg/kg, along with maintaining or reducing long-term average concentrations to 400 mg/ kg or less (RECOVER 2007). 35 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) REMAP data (average of three cores, 0-10 centimeter soil depth) show spatial gradients in soil phosphorus (Figures 29 and 30). These results are similar to landscape data obtained by others for the Everglades in 2003 (Corstanje et al. 2006, Sklar et al. 2006, Bruland et al. 2007, Osborne et al. 2011a) and in 2005 (Osborne et al. 2015). REMAP data indicate no change in TP concentration in soil system-wide in 2014 as compared to 2005 (Wald F, p=0.82; Figure 31). There also was no change for the Park, Refuge, WCA2 and WCA3 subareas. This is in contrast to the 2005 sampling, when REMAP documented higher soil TP in 2005 than in 1995-96 (Wald F, p < 0.05) (Scheidt and Kalla 2007). Other efforts also documented increases in Everglades soil TP from 1990-2003. Spatial expansion of elevated soil TP within WCA2A occurred from 1990-98, such that the WCA2A median changed from 516 mg/'kg to 860 mg/kg (Grunwald et al. 2004). Soii TP data within WCA3A collected from 1992 and 2003 indicated that the area with soil TP > 500 mg/kg increased from 21% to 30% (Bruland et al. 2007). Transect sampling downstream of inflow structures along TP gradients in the Refuge and WCA2A in 1989 and 1999 also indicated expansion of the area with soil TP > 700 mg/kg (Childers et al. 2003). Figure 29. REMAP data from 1995, 1996, 2005, and 2014 for total phosphorus in soil expressed as milligrams phosphorus per kilogram of soil (left) and as micrograms of phosphorus per cubic centimeter of soil (right). 36 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Š Median ~ 25%-75% * X Non-Outlier Range o Outliers * Extremes i LNWR WCA3AN WCA3B ENP OM WCA2 WCA3AS ENP SRS ENP EMM ? II ' b r Š Median ~ 25%-75% I Non-Outlier Range o Outliers * Extremes 1 T LNWR WCA3AN WCA3B ENPOM WCA2 WCA3AS ENP SRS ENP EMM Figure 30. REMAP data by subarea for total phosphorus in soil expressed as micrograms of phosphorus per cubic centimeter of soil (left) and milligrams phosphorus per kilogram of soil (right). Data are from 1995, 1996, 2005, and 2014. Soil Total Phosphorus (mg/kg) Figure 31. Soil total phosphorus (mg/kg) estimates of marsh area for 2005 and 2014 wet season REMAP samples. The percent of marsh area greater than 500 mg/kg did not change significantly from 2005 to 2014 (Wald F test). Landscape-wide, the median REMAP soil TP concentration was 390 mg/kg in both 2005 and 2014, as compared to 343 mg/kg in the 1995-96 wet season, and 325 mg/kg in the 1995 wet season. For 2014, the area of the Everglades with soil TP concentrations exceeding the 500 mg/kg impacted threshold was 21.1 ± 5.6%, which is not statistically different than the 24.5 ±6.4% observed in November 2005 (Wald Ftest, Figure 31). This contrasts with 16.3 ± 4.1% exceeding 500 mg/kg in 1995-96. In 2014, 45.3 ± 7.1% of the marsh 37 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) exceeded the 400 mg/kg CERP restoration goal, as compared to 49.3 ± 7.1% in 2005 and 33.7 ± 5.4% in 1995-96 (Scheidt and Kalla 2007). Small-scale heterogeneity exists, particularly in Park soils. During 2014 each sampling station represented 17.8 square miles of marsh. Because of the sampling density and the random design, gaps in the sampling stations can result in no information for key areas of interest with high soil phosphorus, such as in portions of WCA2A near inflow structures (Figures 27 and 29). Documenting change at locations with the highest vulnerability to enrichment, such as near inflow structures, and documenting restoration success in these areas, necessitates a spatially-focused intensive monitoring design (Cohen et al. 2009a, Osborne et al. 2011b). These designs include transects along phosphorus gradients, or an increased density of sampling stations. Soil percent organic matter and bulk density vary greatly throughout the Everglades due to differences between organic peat soils and inorganic marl soils. Soil bulk density is low in peat soils (typically < 0.12 g/cc), and high in marl soils (median of 0.36 g/cc for Park marl soils, Figures 55-56). Soil phosphorus can be expressed on a mass basis as milligrams of phosphorus per kilogram of soil, or on a volume basis as micrograms of phosphorus per cubic centimeter of soil. When soil TP is adjusted for bulk density (Figure 29 right), TP is more uniform throughout the areas with peat soils, and the locations in WCA3A south of I-75 that were above 500 mg/kg ('impacted') no longer have high TP. Peat soils with higher TP are generally limited to WCA2A and the edges of the Refuge. The places with the lowest bulk density have the lowest soil phosphorus. These observations are consistent with monthly surface water TP data from fixed marsh stations at these locations, which have annual geometric mean TP concentrations <10 |jg/L (McCormick et al. 1999; Julian et al. 2016a, 2016d). Higher soil TP concentrations at some locations may reflect peat versus marl soil types, and may be influenced by soil loss due to oxidation and subsidence and higher bulk density. Soil TP concentrations in the 300-600 mg/kg range may not be an appropriate indicator of enrichment for mineral (marl) soils within the Everglades (USEPA2000). 38 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Nitrogen Nitrogen (N) is the most abundant nutrient in the Everglades and it plays an important role in nutrient cycling. In phosphorus-impacted areas, N can become the nutrient that limits undesirable biological growth (Inglett et al. 2011). Although the concentrations of N that are discharged from the EAA and the STAs are higher than those found in the downstream marsh, overall, N enrichment has not been identified as a concern in the Everglades. Most of the nitrogen flowing into the Everglades is in the organic form, which is not as reactive as inorganic forms such as nitrite, nitrate and ammonia. Nitrogen cycles within water bodies in organic and inorganic forms. Nitrification is the oxidation of ammonium to nitrate, the form of N that is most easily assimilated by many plant species. Denitrification is the reduction of nitrate to N gas, which can leave the water body and enter the atmosphere. Although denitrification is a potential pathway for N loss from Everglades surface waters, this pathway is not major. Ten percent of nitrogen was lost from peat soils to denitrification, while 34% was lost from marl soils (Gordon et al. 1986), and the rate of removal increased as soils became more phosphorus-enriched (Gordon et al. 1986; White and Reddy 2003). The water quality criterion for nitrogen that applies to the Everglades is a narrative: nutrient concentrations shall not be altered so as to cause an imbalance in natural populations of aquatic flora or fauna (FAC 62- 302.530(48)(b)). CERP adopted an Everglades restoration goal for N of less than or equal to the baseline mean observed during 1994-2004; however, the N concentrations for this baseline have not been defined (RECOVER 2007). Most coastal marine systems are N-limited, and therefore N is important as the nutrient that limits undesirable biological growth in the marine portions of the Everglades, such as Florida Bay. It is important to holistically understand restoration efforts and the potential effects that increased freshwater flows may have on nitrogen processes (Inglett et al. 2011). Surface water total nitrogen (TN) during the September 2014 sampling had a spatial gradient, with the highest concentrations above 0.9 mg/L observed in WCA2A and the northern portion of the Refuge, along with locations in the Park closest to the L67 extended canal (Figure 32, right). The median for all 116 sampling sites was 0.77 mg/L. Dissolved nitrite, nitrate and ammonia accounted for a very small amount of the TN (median = 2%). During WY15, TN concentration varied greatly at water control structures, depending upon proximity to the EAA. The WY15 flow-weighted mean annual TN entering the Park was 1.15 mg/L, as compared to 5.46 mg/L for the Refuge. The annual median TN interior marsh concentrations for the Refuge, WCA2, WCA3 and Park were 1.1, 1.6, 1.3 and 1.1 mg/L, respectively (Julian et al. 2016a). Surface water nitrogen gradients throughout the Everglades have been reported since the 1970s. During 1978-82, the five-year mean nitrate concentration entering the Park in western Shark River Slough was 0.012 mg/L, as compared to 0.938 mg/L at the S8 structure in the Miami Canal which was discharging untreated stormwater from the EAA into WCA3A (Mattraw et al. 1987). For 1978-87, the mean TN concentration at pumps discharging stormwater from the EAA were 3.4-6.0 mg/L, with inflows to the Park at 2.0 mg/L (Scheidt 39 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 32. REMAP soil total nitrogen (percent) during 2005 and 2014 (left), and surface water total nitrogen (mg/L) during September 2014 (right). et ai. 1989). A decreasing trend in TN throughout the Everglades during WY79-WY15 has been reported, especially at inflows to the Refuge (-0.14 mg/L/WY) (Julian et al. 2016a). During WY15 the inflows to STAs 1E, 1W, 2, 3/4 and 5/6 had flow-weighted mean TN concentrations of 2.4, 2.8, 2.5, 2.9, and 1.6 mg/L, respectively, while the outflow concentrations were 1.6,1.8,1.6,1.5 and 1.4 rng/L. The STAs removed a minimal to moderate amount of TN, with WY15 treatment efficiencies of 13% to 49%. The STAs removed most of the nitrite plus nitrate, which was a small portion of TN (SFWMD 2016b). The five-year TN treatment efficiency for STA-1W was 26%, as compared to 79% for TP (Gu et al. 2006). The median TN content of soils sampled during REMAP 2005 and 2014 was 2.9% (n=344, Figures 32-33). Content greater than 3.5% tended to occur in the highly organic peat soils found in the interior of the Refuge and the longest hydroperiod portion of southern WCA3A. In contrast, the marl soils of the Park generally have a TN content of <1.5%. These REMAP data have the same spatial gradients documented throughout the Everglades by SFWMD in 2003 (Orem et al. 2014b, Osborne et al. 2015), and the Park 40 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) a> o> o z 3 2 ,o o iD i Š T Š I Š Median ~ 25%-75% I Non-Outlier Range o Outliers * Extremes T T 7 LNWR WCA3AN WCA3B NESS EMM WCA2 WCA3AS SRS OM Figure 33. and 2014. Soil total nitrogen concentration (percent) by subarea, 2005 results are consistent with the SFWMD results at 342 locations (Osborne et al. 2011b). The peat soils of the EAA, which originated from decomposing Everglades sawgrass prior to drainage in the 1900s, also have a TN content of about 3% (Inglett et al. 2011). Soils within STA2 and STA3/4 have mean TN content of 2.8% and 2.4% (Pietro and Ivanoff2015). The major source of agricultural nitrogen in the EAA is the soil itself. No fertilizer additions of nitrogen are necessary for sugarcane and minimal additions are necessary for vegetables (de Camargo Santos 2020, Morgan etal. 2009). Drainage water from the EAA has been reported to export TN at rates ranging from 30-46 kg N/hectare/year (Porter and Sanchez 1994) and 12-40 kg N/hectare/year (Gilbert and Rice 2006). 41 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Specific Conductance Specific conductance is an indirect measure of the overall concentration of dissolved ions or minerals in water, and is defined as the electrical conductance of 1 cubic centimeter (cm3) of a solution at 25 degrees Celsius. The terms "specific conductance," "electrical conductivity," or simply "conductivity" have been used interchangeably in recent decades to report the same measurement. The units of electrical conductivity are mhos per centimeter (mhos/cm), which are equivalent to Siemens per centimeter. Pure water has a very low electrical conductivity of a few hundredths of a micromhos per centimeter (ie., less than 0.05 jjmhos/cm) (Hem 1985, USGS 2019). Specific conductance provides information about the total dissolved ion content of water, but not the ion composition. The terms 'soft water' and 'hard water' are used to describe waters containing relatively low or high concentrations of dissolved ions. As the concentration of dissolved ions increases, so does the conductivity. In the Everglades, conductivity is useful for understanding the source of water and its flow path. The canals that provide flood control for the EAAcut into the shallow aquifer, which is highly mineralized and begins at a depth below the ground surface of only six to ten feet. Conductivities in this aquifer at a depth of 20 feet vary from about 500 to several thousand jjmhos/cm. From the 1940s to 1980s, there was an increase in the mineral content of the shallow aquifer due to the upward migration of groundwater, a response to removal of surface water by pumping for flood control (Miller 1988). During 1997-2003, the median conductivity at 10 farm canals within the EAA ranged from 770-1670 jjmhos/cm, as compared to 600 jjmhos/cm for Lake Okeechobee. The highest values within the EAA occur in the eastern portion within the S5A and S6 basins. During the wet season, farm canals upstream of the Refuge have conductivities above 1200 jjmhos/cm due to the seasonal input of groundwater during drainage pumping for flood control operations (Bhada et al. 2014). During 1974, when water from the EAA was pumped north into Lake Okeechobee, surface water conductivity was about 1000-1400 jjmhos/cm in canals within the EAA, whereas lake conductivity decreased at a gradient to about 500-800 jjmhos/cm toward the interior (Brezonik and Federico 1975). Conductivities in canals surrounding the Refuge and in the EAA are higher due in part to the upward migration of residual (connate) seawater (Miller 1988). The major ions in EAA farm canal water contributing to conductivity are chloride, sulfate, carbonate, sodium and calcium (Chen et al. 2006). Pronounced spatial gradients in surface water conductivity, sulfate and chloride throughout the canal and marsh system vividly demonstrate that the canal system can be a conduit for transport of degraded water. Conductivity varies with the number and types of ions in solution. The water in the interior marsh of the Refuge is soft, slightly acidic (median pH = 6.1, REMAP data from September 2014), and strongly influenced by rainfall (precipitation pH=5.1, NADP 2014a). The limestone (calcium carbonate) substrate underlying the Refuge is overlain by several feet of peat soil, so surface water is not in contact with the limestone. In contrast, the rest of the Everglades marsh has hard water with a neutral pH (median = 7.4). In the shorter hydroperiod portions of the Park there is little soil, so surface water has higher ionic content 42 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) because it is influenced by the limestone substrate (Figures 53 and 58). Conductivity of water is closely related to its hardness, because calcium, the major contributor to hardness in the Everglades, also aids in conductance. Conductivity is of ecological interest in the Everglades because it is one of the determinants of the composition of periphyton communities, which play a critical role in food webs, habitat, soil formation and water quality (McCormick et al. 2011). REMAP data document pronounced conductivity gradients and transport of high conductivity water by canals well into the Everglades marsh and the Park (Figures 34-35, Stober et al. 1998, Scheidt et al. 2000, Scheidt and Kalla 2007). These gradients are due to the relative influence of rainwater, groundwater, and stormwater inflows, and they indicate pathways of water flow throughout the canal-marsh system and the extent to which the water management infrastructure and operations influence water quality. Precipitation in the Everglades has very low ionic content, with 2014 average annual specific conductivity of 10.0 |jrnhos/cm within the Park (NADP 2014a). In contrast, the conductivity of water discharged from the EAA into STAs 1E, 1W and 2 during the September 2014 REMAP sampling was 1000-1200 pmhos/cm, with STA5/6 having the lowest inflow conductivity of about 540 p mhos/cm. The highest marsh conductivities of about 900 pmhos/cm were within the Refuge and WCA2 near the canals due to proximity to the EAA and the influence of canal water and groundwater. The WY15 annual flow-weighted mean conductivity in water discharged from the Stormwater Treatment Areas (STAs) into the Refuge was 930 pmhos/cm from STA1E and 863 pmhos/cm from STA1W. The average conductivity in annual discharges from the other STAs were: STA2 887 pmhos/ 10 I 20 Miles Canals I I Stormwater Treatment Areas I I Miccosukee Tribe Federal Reservation D USEPA REMAP Study Area Conductivity |jmhos/cm O £200 o >200 - 400 O >400 - 600 O >600-800 O >800 Figure 34. Surface water conductivity (specific conductance, pmhos/cm) during September 2014 at 116 REMAP marsh stations. Also shown are SFWMD data for 11 iocations in Lake Okeechobee and 70 locations in the Evergiades marsh (geometric mean), and 32 water control structures in canals (geometric mean when flowing) including selected STA inflows and outflows. 43 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) cm, STA3/4 735 (jmhos/cm, and 456 |jmhos/cm from STA5/6. As water passes through the STA treatment wetlands, designed and operated for phosphorus removal, there is little change in conductivity. Conductivity in the STA discharges was within 10% of the inflows, except STA1E where the discharge conductivity was 21% higher than in the inflow (SFWMD 2016b). Lower wet season conductivity in the interior of the Refuge (about 70 jjmhos/cm) and the western portion of WCA3A (about 250-300 pmhos/cm) reflect the greater influence of rainwater as compared to canal water (Figures 35 and 36). Marsh conductivity increases in the dry season due to lessening dilution by rainwater, evapo-concentration as the marsh dries out, and greater influence of canal flows (Scheidt et al. 2000). From 1959 to 1974, as inflow to the Park at Shark River Slough changed from being dominated by marsh sheetflow to canal discharge at the new S12 structures, wet season mean marsh conductivity rose from 270 to over 500 pmhos/cm (Flora and Rosendahl 1982a, 1982b). During 1978 to 1982, conductivity varied spatially such that at structure S12A, a gated spillway that discharges water into the Park at western Shark Slough, conductivity averaged 303 pmhos/cm, as compared to 1184 pmhos/cm entering WCA3A in the Miami Canal at the S8 structure 39 miles to the north (Mattraw et al. 1987). Conductivities observed in WCA3 and WCA2 during the 2014 REMAP sampling were lower than those observed during 2005. The maximum in WCA2A during 2014 was 780 pmhos/cm, as compared to 1423 ^mhos/cm in 2005 (Figure 35). WCA2 had a median of 1041 |jmhos/cm in 2005 versus 676 pmhos/cm in 2014, and the Refuge had a median of 151 |jmhos/cm in 2005 versus 77 |jmhos/cm in 2014. Figure 35. Surface water conductivity (|jmhos/cm) in the marsh during REMAP September 1995, November 2005 and September 2014 wet season sampling events. 44 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 1500 E S 1000 500 LNWR WCA2 WCA3N WCA3S ENP Wet/Dry season: E|3 Dry ^ Wet Figure 36. Box and whisker plot of surface water conductivity (|jmhos/cm) during the wet and dry seasons for all REMAP marsh sampling events. Surface water conductivity, chloride, sulfate and dissolved organic carbon were all lower in 2014 than in 2005. These differences could be explained by less discharge from the EAA into the Everglades preceding the 2014 sampling. September 2014 REMAP data indicate that as expected marsh conductivity is correlated with surface water chloride, sulfate and dissolved organic carbon [Spearman rank- order correlation coefficients of 0.95, 0.71, and 0.65 respectively, p<0.001 (Kalla and Scheidt 2017)]. A long-term decrease in the conductivity at inflows to the Refuge and WCA2A since 1979 has been reported (Julian etal. 2016a), as has a decrease in conductivity from 1977 to 2005 at some of the structures that release water into the Park at Shark River Slough (Fan et al. 2010). The discharge of high conductivity water into the Refuge marsh has been a concern since the 1970s, when the penetration into the Refuge marsh of highly mineralized water at 10 to 20 times background conditions was documented (Gleason and Spackman 1974, McPherson et al. 1976). Concern was raised about the impacts of this mineralized water on Refuge biota, such as periphyton, that are adapted to low conductivity, soft water conditions (Gleason and Spackman 1974, Gleason et al. 1975a). Florida's Class III water quality criterion for conductivity that applies to the Everglades is: "shall not be increased more than 50% above background or to 1275 micromhos/cm, whichever is greater" (Chapter 62-302.530(22) F.A.C). In the annual South Florida Environmental Report, Florida reported conductivity excursions for 2014 that exceed the Class III criterion and are of concern for inflows to the Refuge and the WCA2A marsh (Julian et al. 2016a). The Park and Refuge are also Outstanding Florida Waters, which requires that the water quality condition that existed in these waterbodies during the year prior to March 1, 1979 must be maintained. Background conductivity within the interior of the Refuge is less than 100 jjmhos/ cm. Highly mineralized water penetrates several kilometers into the Refuge marsh (Harwell et al. 2008; Surratt et al. 2008; McCormick et al. 2011). Data from transects in the Refuge marsh downstream of the STAs indicate that during WY15 water entered the Refuge downstream of STA1E at about 1000 jjmhos/ cm (annual geometric mean), and dropped to 200 jjmhos/cm at about 3.5 kilometers (km) into the marsh. 45 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Water entered the Refuge from STA1W at about 800 jjmhos/cm, and dropped to 200 jjmhos/cm by about 1.0-2.5 km into the marsh, depending upon the transect. In contrast, water entered WCA2A at about 900 jjmhos/cm and by 5 km into the marsh there was no decrease along the three transects (SFWMD 2016b). The value of periphyton communities as a food source is affected by conductivity, in that increases in water ionic content can shift periphyton community structure (Browder et al. 1993, Sklar et al. 2005b, McCormick and Harvey 2007, McCormick 2011, McCormick et al. 2011, Gaiser et al. 2011, Gottlieb et al. 2015). Periphyton communities in the Refuge tend to be attached to plants (epiphytic) and are dominated by desmid and diatom species, while the extensive periphyton mats (Figure 63 and 64) in hard water portions of the Everglades are dominated by calcium-precipitating cyanobacteria with a high calcium carbonate content (McCormick and Harvey 2007; McCormick et al. 2011). CERP developed an Everglades protection and restoration performance measure for conductivity of no more than a 25% increase above background while taking into consideration natural seasonal and annual variation (RECOVER 2007). A restoration goal is that the natural soft, low-conductivity surface water will be maintained in the Refuge, while the hard, higher conductivity water consistent with natural background levels in much of the rest of the Everglades will also be maintained. Everglades restoration should take into consideration ways of minimizing the intrusion of mineral-rich water into the Refuge interior (Newman and Hagerthey 2011). However, given the inevitable groundwater-surface water interaction due to the presence of canals, elevated surface water conductivity is unavoidable to some extent. Water management operations have been developed to help reduce canal water intrusion in the Refuge interior (Surratt et al. 2008). 46 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Chloride Chloride is an ion that is found in surface water and groundwater. Chloride ions can come from sodium chloride or from other chloride salts. Saltwater has very high natural concentrations of chloride. In coastal areas, estuarine waters contain more chloride than freshwater due to the influence of saltwater. Chloride is a conservative ion, meaning its concentration does not readily change. Therefore, chloride can be a useful indicator of a water's source. The concentration of chloride varies greatly throughout the Everglades depending upon the relative influence of rainwater, groundwater and stormwater. During 2014, precipitation in the Everglades had a precipitation-weighted annual chloride concentration of 1.0 milligram per liter (mg/L) (NADP 2014a). (The chloride concentration was measured in weekly samples of precipitation, and when calculating the average annual chloride concentration, the chloride results are weighted by the volume of precipitation during that week in order to get an accurate annual average.) During the September 2014 wet season sampling, the lowest surface water chloride concentrations of <15 mg/L were observed in the interior of the Refuge and the western portion of the Park and southwestern WCA3 away from canal inflows. The highest concentrations, exceeding 100 mg/L, were observed at some STA inflows and discharges and in the Refuge and WCA2A marshes near inflows (Figures 37-38). The median concentration in the Refuge was 12 mg/L, as compared to 83 mg/L in WCA2. While these spatial patterns are similar to those observed in 2005, the 2014 medians Canals Chloride Stormwater Treatment Areas mg/L O S40 Miccosukee Tribe Federal Reservation 0 5 10 20 Miles Q USEPA REMAP Study Area O >70 100 1 i i i l i i i l w O >40 - 70 >70 - >100 Figure 37. Surface water chloride concentration (mg/L) during September 2014 at 116 REMAP marsh stations. Also shown are SFWMD data for 6 locations in Lake Okeechobee, 75 locations in the Everglades marsh and 32 water control structures in canals (geometric mean when flowing) including selected STA inflows and outflows. 47 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Pronounced spatial gradients in surface water conductivity, chloride, and sulfate throughout the canal and marsh system vividly demonstrate that the canal system is a conduit for transport of degraded water. This transport is an unintended consequence of the flood control project. Figure 38. Surface water chloride (mg/L) in the marsh during the REMAP November 2005 and September 2014 wet season sampling events. were lower (2005 Refuge median 26 mg/L and WCA2 median 140 mg/L, Scheidt and Kalia 2007). Chloride concentrations in the northern portion of WCA3A near the Miami Canal were also lower in 2014 than 2005, and may reflect somewhat greater discharge from the north and influence of canal water in 2005 (121,788 ac-ft at inflows into WCASAat S8 during Octoberto November 2005 versus 106,358 ac-ft in August to September 2014, Abtew et al. 2007; Abtew and Ciuca 2016). Over the duration of REMAP sampling events, wet season chloride concentrations within each subarea were lowerthan dry season concentrations, presumably because of the diluting effect of rainfall (Figure 39). The groundwater chloride concentration in the shallow aquifer within the EAA at a depth of 20 feet is typically reported between 100-200 mg/L. The chloride concentration within this shallow aquifer increased 48 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) from the 1940s to the 1980s due to the upward migration of groundwater in response to pumping for flood control (Miller 1988). During 1999-2003 the median chloride concentration at 10 farm canals within the EAA was 72- 174 mg/L (Chen et al. 2006). Chloride concentrations reported for 1972-73 in the Hillsboro Canal and the downstream WCA2 marsh both averaged about 170 mg/L, with some concentrations as high as 500 mg/L (Gleason 1974). From 1959 to 1974, as inflow to the Park at Shark River Slough changed over time from being dominated by marsh sheetflow to canal discharge at the new S12 structures, canal chloride concentration rose from about 20 mg/L to 60 mg/L (Waller 1982). Pronounced gradients in chloride concentration within the Refuge and WCA2 marshes have existed since the early 1970s (Gleason 1974, Gleason et al. 1975b). From 2000-06, marsh stations within 1 km of the canal had median chloride concentrations of 100-130 mg/L, as compared to about 20 mg/L for marsh stations father from the canal. A hydrodynamic and constituent transport modelling effort found that interior marsh chloride originated from pumped inflows ratherthan rainfall (Chen et al. 2012). Chloride concentrations in the Refuge perimeter canal during 1995-2007 ranged from about 50-170 mg/L, while concentrations in the interior marsh typically ranged from 10-50 mg/L, but exceeded 100 mg/L during particular events. Chloride was used as a conservative surface water constituent to develop a water quality and hydrology model for the Refuge. Pumped inflow was found to be the dominant factor responsible for the substantially increased chloride and sulfate concentrations in the marsh (Wang et al. 2012). Chloride is one of the constituents of highly mineralized water that is of ecological concern to the naturally soft-water Refuge. The STA wetland treatment systems that are designed and managed to remove phosphorus generally do not remove a large proportion of the chloride. The chloride concentration in the STA discharges into the Refuge from STA1 Wgenerally varied between 100 to 200 mg/L from 1994-1999 (Gu etal. 2006). ForWY06, during the 2005 REMAP sampling, STA 1W, STA 2, STA 3/4, and STA 6 discharged dissolved chloride at concentrations of 142, 157, 73 and 22 mg/L respectively, and outflow concentrations generally were only slightly lower or higher than inflow concentrations (Pietro et al. 2007). 200 Wet/Dry season: Dry ^ Wet Figure 39. Wet and dry season surface water chloride (mg/L) by subarea in the marsh during all REMAP sampling events. 49 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Sulfur Sulfur is an element that exists in several forms in water bodies. Sulfur generally occurs in surface water in the oxidized state as sulfate, an ion that is common in nature. Sulfate is a natural ingredient of rainfall, surface water and groundwater. The form ofsulfurthat exists in the environment where there is no oxygen is the reduced form, sulfide, which is associated with sulfate reduction by anaerobic bacteria. Sulfur is applied to condition soils in the EAAfor crops (Morgan et al. 2009). Sulfur is of interest in the Everglades for three reasons: sulfate and sulfide are associated with mercury methylation and subsequent biomagnification in gamefish and wildlife (Fink and Rawlik2000, Jeremiason et al. 2006, Orem et al. 2019, Pollman 2014,2019; Rumbold 2019a); elevated sulfate can mobilize phosphorus in some water bodies (Smolders et al. 2006, Lamers et al. 1998, Lamers et al. 2002, Beltman et al. 2000); and sulfide at elevated concentrations can be toxic to or adversely affect plants (Smolders et al. 2006, Lamers et al. 1998, Lamers et al. 2002; Li et al. 2009) and animals (USEPA 1986). Because of these ecological concerns CERP adopted a restoration performance measure for surface water sulfate: maintain or reduce sulfate concentration to 1 milligram per liter (mg/L) or less throughout the Everglades marsh (RECOVER 2007). Florida does not have specific water quality criteria for sulfate or sulfide in the Everglades, and USEPA does not have a recommended surface water criterion for sulfate. For sulfide in surface water, USEPA recommends a continuous concentration criterion of less than 0.002 mg/L for protection of aquatic life (USEPA 1986). Florida water quality standards require that "Substances in concentrations which injure, are chronically toxic to, or produce adverse physiological or behavioral response in humans, plants or animals - none shall be present." (Chapter 62-302(62) F. A. C.). This is referred to as the 'free from' requirement. In addition, Florida also designated the Park and Refuge as Outstanding Florida Waters, requiring that the water quality that existed as of March 1, 1979 be maintained. Sulfate concentration varies throughout the Everglades marsh depending upon proximity to the EAAand canals, and the relative influence of rainwater, stormwater and groundwater. Pronounced marsh and canal surface water sulfate gradients have been documented independently by USEPA REMAP during 1993 to 1996, 1999 and 2005 (Scheidt et al. 2000, Stober et al. 2001b, Scheidt and Kalla 2007), USGS (Orem et al. 2011) and SFWMD (Julian et al. 2016b). Earlier studies documenting sulfate gradients during the 1970s are summarized in Scheidt and Kalla (2007). Surface water sulfate is highest during the wet season with the movement ofwaterthroughoutthe Everglades forflood control. Elevated sulfate levels in 2014 followed this same landscape pattern (Figures 40-42). Concentrations varied from <0.022 mg/L (the analytical laboratory method detection limit (MDL) in 2014) at marsh locations away from canals, to 75 mg/L and 80 mg/L at the inflows to STAs 1W and 1E near the Refuge (Figure 40). Concentrations in southern Lake Okeechobee averaged 40 mg/L. Figure 41 shows marsh stations sampled by REMAP during the 1995, 2005 and 2014 wet seasons. The highest marsh concentrations are at locations that are proximate to canals or stormwater flows from the EAA. Concentrations in the Everglades progressively decrease to the south and west. These landscape patterns indicate that the canal system delivers sulfate from the north into Everglades marshes, 50 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Canals Sulfate Stormwater Treatment Areas m9"- I I Miccosukee Tribe Federal Reservation O O >1-10 10 20 Miles ~ USEPA REMAP Study Area O >10-35 _J i i i I O >35 Canals Sulfate (mg/L) I I Stormwater Treatment Areas 0*1 . , O >1 - 1° | | Miccosukee Tribe Federal Reservation >< ^ ? I 4'J5 9;5 , , 1|9Mlles USEPA REMAP Study Area Q >35 Figure 40. Sulfate (mg/L) in surface water during September 2014 at 116 REMAP Everglades marsh stations. Also shown are SFWMD data for 8 locations in Lake Okeechobee, 60 locations in the Everglades marsh (geometric mean) and 16 water control structures (geometric mean when flowing) in canals, including selected STA inflows and outflows. The top view is from the east, the bottom view is from the west. 51 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) with penetration into the Shark River Slough marsh within the Park (Figure 40). The influence of canal water flowing into the marsh in the wet season was more apparent in 1995 and 2005 than in 2014. Sulfate concentrations in rainfall and interior marsh locations away from canals are lower. The annual precipitation-weighted sulfate concentration in rainfall within the Park for 2014 was 0.69 mg/L (NADP 2014a). It was less (0.47 mg/L) during August and September of 2014 around the time of the REMAP sampling (NADP 2014b). During the September 2014 REMAP sampling, interior portions of the Park, Refuge and WCA3 that are most influenced by rainfall had 21 stations with sulfate concentrations in surface water reported at the analytical laboratory MDL of about 0.02 mg/L (Figures 40-41), indicating that the true concentration may be even lower. There are 72 stations sampled in 2014 in Figure 41 within the <2 mg/L isopleth that were <1 mg/l, with a median concentration of 0.08 mg/L. Sulfate is one of the constituents of highly mineralized water that is of ecological concern to the naturally soft-water Refuge. There are pronounced surface water sulfate gradients from the canals to the Refuge interior marsh (Harwell et al. 2008). Water and soil chemistry data collected along an east-west transect show an extensive zone of sulfur enrichment associated with episodic canal-water intrusion. Natural background concentrations are <1 mg/L for surface water sulfate and <1% for soil sulfur. Concentrations are elevated to as high as 9 mg/L and 2% at sites within about 5 km of the western and eastern rim canals (McCormick et al. 2011). Soil sulfur is lower in the Refuge interior, with much higher concentrations along the marsh periphery (Osborne et al. 2011a). Pumped inflow is the dominant factor responsible for the substantially Figure 41. REMAP surface water sulfate concentration (mg/L) during September 1995 (left), November 2005 (center) and September 2014 (right). 52 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) increased sulfate and chloride concentrations in the interior marsh (Wang et al. 2012), ratherthan atmospheric deposition (Wang et al. 2009). Residual seawaterthat contains high sulfate also influences canal sulfate concentrations in the vicinity of the Refuge and WCA2A (Pollman 2012). The wetland STAs designed and managed to remove phosphorus remove little sulfate (Orem et al. 2011). During WY15, forthe STAs the flow-weighted sulfate inflow and outflow concentrations, and percent removal were as follows: STA1W 68, 66 mg/L, 3%; STA1E 56, 57 mg/L-3%; STA2 51, 45 mg/L, 12%; STA3/4 48, 46 mg/L, 5%; and STA5/6 10, 2 mg/L, 77% (SFWMD 2016b). Sulfate removal by the STAs was about 10% or less, with the exception of STA5/6 which had a much lower inflow sulfate concentration. STA1W exhibited moderate removal of sulfate from 1994-99 (Gu et al. 2006). The eastern EAA (the S5A, S2 and S6 basins) has consistently had the highest concentrations of sulfate in groundwater and surface water. If the high sulfate within the STAs mobilizes phosphorus, this may hinder phosphorus removal by STAs, especially for STAs 1W, 2 and 3/4. The concentration of sulfate in Everglades groundwater is higher than surface water natural background levels. Sampling of the surficial aquifer underlying the EAA at about 20 locations in 1983-84 found sulfate concentrations of 25-580 mg/L at a groundwater depth of 45 feet, about 20 feet below the depth that the major canals penetrate. These high concentrations may be due in partto residual seawater (Miller 1988). The highest concentrations were in the eastern EAA in the area of the S-2 and S-6 basins. In 1976-77, sulfate was 20- 490 mg/L in shallow groundwater in the EAA, with mean concentrations of 153 mg/L beneath sugarcane and 199 mg/L beneath vegetables. The mean surface water concentrations ranged from 40-459 mg/L (CH2MHILL 1978). In contrast, the median groundwater concentration in 189 wells tapping the Biscayne Aquifer was 17 mg/L (Radell and Katz 1991). The Biscayne Aquifer is the shallow, unconfined, highly-permeable aquifer underlaying the Everglades and southeast Florida. Agricultural sulfur (S) has been applied to EAA soils for various purposes. EAA soils have been prone to copperdeficiency, which has been addressed by treatment with copper sulfate. Magnesium has been supplemented by the use of fertilizer blends containing potassium-magnesium sulfate (Anderson 1990). The sulfur content of EAA peat soils is considered WCA2 WCA3N WCA3S Wet/Dry season: Ej3 Dry E$3 Wet Figure 42. REMAP surface water sulfate concentration (mg/L) by subarea during the wet and dry seasons for all REMAP marsh sampling events. 53 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) adequate to supply some S requirements. However, surface application of S has been recommended when soil pH is > 6.6 in order to increase plant nutrient availability, with a recommended application rate of 500 pounds S per acre (Coale 1994, Schueneman and Sanchez 1994). It has been reported that grower S application rates are lower than recommended rates (Schueneman, 2001). Application rates have been estimated at 10 pounds per acre (Ibs/ac) to sugarcane and 78 Ibs/ac to vegetables (CH2MHILL 1978). Recent agricultural studies suggest that S application may have minimal benefits for increasing nutrient availability (Ye et al. 2009, 2011a). As EAA soils subside, the pH of the remaining soil tends to increase, which may result in a need for higher S application rates (Ye et al. 2011 b). This could increase S export to the Everglades (Ye et al. 2010). Investigators analyzing sulfur concentrations and isotopic ratios for rainwater, EAA groundwater, and EAA fertilizer concluded that excess sulfate in the Everglades originates from canals draining the EAA (Bates et al. 2002). The sulfate concentration and isotopic data appear to exclude rainwater and some ground water as major contributors. Isotopic evidence implicates agricultural applications as a major contributor to the sulfate load. This fertilizer could be recent additions, legacy deposits, or both. However, EAA groundwater and oxidation of agricultural soil may also contribute sulfate (Bates et al. 2002). It has been reported that, based on isotopic composition, groundwater is not a major source of sulfate to surface water in WCA2A (Gilmour et al. 2007). Sulfate in Lake Okeechobee was about 40 mg/L from 1975-2001, and declined over the last three decades (James and McCormick 2012). About 97% of the sulfate in Lake Okeechobee originates from the watersheds that discharge to the lake, including backpumping from the EAA, with rainfall and atmospheric deposition contributing only 3% of the sulfate load. Lake water provides about 20% of the sulfate to the EAA, and smaller loads to the Everglades. There are several mass balance estimates of S inputs to the EAA and outputs to the Everglades (Corrales et al. 2011, Gabriel et al. 2014a, Landing 2015). Sulfur outputs to the Everglades either equal or exceed inputs to the EAA, depending upon assumptions (see summary in Landing 2015). The importance of controlling sulfate to the lake increases if more lake water is moved south to the EAA for water supply and environmental enhancement because increased flow from the lake could result in higher loads to the EAA and Everglades (James and McCormick 2012), especially since STAs do little to remove sulfate. Elevated surface water sulfate concentrations were found during all REMAP sampling events, with varying extent. Sulfate was highest during the 1995 REMAP wet season sampling event, lower in 2005 and lowest in 2014 (Figures 41,43-44). The laboratory analytical MDL for sulfate improved greatly between 1995 and 2005, from 2.0 mg/L down to 0.02 mg/L, so the apparent differences between years are exaggerated. In order to account for this, the minimum value for all years was changed to the highest MDL of 2 mg/L (Figure 43). Despite these data being left-censored at 2 mg/L due to laboratory MDLs, and the overlap of confidence intervals, all three CDF curves are different (Wald F, p<0.05). The Park, the Refuge, and WCA3 had less sulfate in surface water in 2014 than in 2005 (Wald F, p<0.02). Based on the CDF of REMAP data, during September 2014 the proportion of the Everglades marsh where sulfate exceeded the 1.0 mg/LCERP restoration goal was 37.1 ± 6.0%, lower than 57.3 ± 6.0% in 2005. 54 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 43. Cumulative Distribution Function (CDF) showing surface water sulfate concentration (mg/L) estimates of marsh area for 1995, 2005 and 2014 wet season REMAP samples. The areas reflect different conditions (Wald F test). A significant decreasing trend forwet (rainfall) and dry atmospheric deposition of sulfate has been reported from WY95-WY15, the duration of REMAP sampling. Air emissions controls may be a factor in this decrease (Julian et al. 2017a). No overall trends are apparent in annual sulfate at surface water inflows to the Refuge, WCA2, WCA3 or the Park from WY95-WY15 (Julian et al. 2016b). At a marsh location in central WCA3A sampled by USGS from 1995 to 2007, sulfate declined significantly beginning in 1998 and remained below 1.0 mg/L through 2007. The cause for the decline was not known (Orem et al. 2020). The lower surface water sulfate concentration observed by REMAP in 2005 compared to 1995 is not likely due to dilution because the lower concentrations observed during the 2005 wet season occurred in shallower water than in 1995 (Figure 21). Less loading from stormwater is a possible explanation, as stormwater inflow to the EPA in the 60 days prior to the 1995 wet season sampling was double the inflow during the 60 days prior to the 2005 wet season sampling (Scheidt and Kalla 2007). The lower concentrations in WCA2A during 2014 as compared to 2005 are suggestive of less stormwater discharge. Concentrations in water may vary within the wet season depending upon particular rainfall and discharge events, so the sulfate concentrations observed during a two-week REMAP sampling snapshot may not represent conditions observed during other weeks in the same wet season. This transient surface water quality characteristic is also applicable to other water quality constituents such as nutrients and ions. USGS found that from 2008-2013 sulfate concentrations within the Park at SRS were influenced by canal water, and there was high interannual variation due to flow (Maglio et al. 2015). 55 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 26 24 22 20 o» 18 J 16 m $ 14 o « 12 tr W 10 Š Median ~ 25%-75% X Non-Outlier Range Š Š Figure 44. Box-and-whisker plots of sulfate in surface water (mg/L) iri the 1995, 2005 and 2014 wet season REMAP samples by year. In 1995 the lowest possible value was the MDL of 2 mg/L. are no water quality criteria for sulfide in pore water. Under anaerobic conditions, sulfide is formed from sulfate. Sulfur speciation and isotopic composition of Everglades plant materials suggests that sulfate reduction occurs in the periphyton mat (Bates et al. 1998). Higher sulfide can have the benefit of inhibiting mercury methylation (Jeremiason et al. 2006, Benoit et al. 1999), but it can also have the detriment of being toxic to macrophytes (Lamers et al. 1998). Sulfide levels in Everglades soils within WCA2 are consistent with concentrations that negatively affect sawgrass (Li et al. 2009). Others suggest that P limitation could override effects of sulfide toxicity, although the effect of sulfide accumulation on P-enriched areas is unknown (DeBusk et al. 2015). There In 2014, REMAP measured the sulfide concentration in water at the bottom of the surface water profile immediately above the soil (bottom water, Kalla and Scheidt 2017). In 2005, sulfide was measured in the water within the top portion of the soil (pore water). In the Everglades, sulfide in bottom water can be predicted by sulfide in pore water (r2 = 0.82, p < 0.001, Kalla et al. 2017). During 2005, sulfide concentration in pore water had pronounced spatial gradients, with concentrations > 5 mg/L in WCA2A. A lab study found that small Everglades sawgrass plants were adversely affected by sulfide at concentrations above 7 mg/L (Li et al. 2009). During 2014 the highest bottom water sulfide concentration in WCA2A was only 0.56 mg/L, and 90% of the values were < 0.1 mg/L. These lower sulfide concentrations are consistent with the lower surface water sulfate concentrations observed in 2014. Bottom water sulfide was correlated with surface water sulfate and dissolved organic carbon (DOC) (Kalla and Scheidt 2017). In 2005, pore water sulfide also was positively correlated with pore water DOC and methylmercury in surface water, periphyton and soil. Fewer correlations observed in 2014 may be due to the much lower range in sulfide concentrations. Others also found that surface water sulfate was positively correlated with porewater sulfide due to reducing conditions favoring conversion of sulfate to sulfide, and porewater DOC significantly influenced porewater sulfide concentrations, suggesting that organic substrate supply could be a factor that affects sulfate reduction (Julian et al. 2016e). Sulfur has been raised as a concern in Everglades restoration regarding Aquifer Storage and Recovery (ASR), a technology proposed near Lake Okeechobee to help meet the ecological and water supply needs of south Florida and the Everglades. During the wet season surface water is pumped into the aquifer and recovered later during the dry season to supplement water supply. The water recovered from ASR is anticipated 56 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) to have high sulfate (200-550 mg/L), which could stimulate microbial sulfate reduction and methylmercury production. Modeling indicates that ASR release would temporarily elevate sulfate concentrations in Lake Okeechobee by 60%, but that this would have little impact on methylmercury production in the lake. The Lake's sandy soils are not conducive to methylation. ASR water released to surface water would have minimal impacts on sulfate loading to the Everglades, primarily due to the much higher sulfate loading from other sources within the EAA. Overall impacts due to increased methylmercury production risk were predicted to be low, although locations in the Everglades near canals or STA discharges could experience significantly higher sulfate loading from ASR and some change in mercury risk (Orem et al. 2014a). The role of sulfate in the biogeochemical cycling of mercury in the Everglades is complex and has been under investigation since the 1990s. REMAP data from 1995-2005 were used to explore the relationship between mercury in prey fish (mosquitofish) and many biogeochemical analytes by using structural equation modeling. This approach provides a framework for selecting or rejecting hypotheses or assumed relationships using data acquired by observation (empirical data). The role of sulfate was found to be complex in terms of its total effect on mosquitofish mercury and was one of several factors that were found to be important (Pollman 2014). Modeling that used REMAP data found that sulfate was one of several factors that influenced mercury concentration in mosquitofish (Kalla et al. 2019, 2021). In an Everglades STA that is managed to remove phosphorus, inflow sulfate along with chloride had a significant correlation with mosquitofish mercury, suggesting that sulfate input may be a factor influencing mercury (Feng et al. 2014). Inflow sulfate, chloride, DOC, and dissolved oxygen were factors related to outflow methylmercury and total mercury (Zheng et al. 2013). Sulfide reacts with dissolved organic matterto form organic sulfur, which impacts mercury methylation and subsequent bioconcentration (Poulin et al. 2017). Modeling has been conducted to predict changes in mosquitofish mercury that would result from potential changes in sulfate export from the EAA. Reductions in excess sulfate were projected to result in, depending on location, either increases or decreases in mosquitofish mercury, with the overall shifts in mosquitofish mercury expected to be small, regardless of the magnitude of reduction in sulfate (Pollman 2012). Whether to manage, mitigate or control sulfate because of concern about impacts to the Everglades has been debated since the 1990s. Many scientists state that it is clear sulfate contributes to mercury methylation and although high and low fish mercury are found across the spectrum of surface water sulfate concentrations, peak fish mercury occurs from about 1 to 12 or 20 mg/L sulfate, and decreasing sulfate loading to background sulfate, or <1 mg/L, would reduce methylmercury risk (Gabriel et al. 2014a, 2014b; Orem et al. 2011, 2019, 2020; Corrales et al. 2011; Pollman and Axelrad 2014, Rumbold 2019a). Other scientists argue that: a) sulfate is one of many factors that can influence mercury accumulation in gamefish such as largemouth bass; b) multiple water quality factors such as pH, specific conductivity and alkalinity also may be factors; c) there is not enough quantitative information to justify sulfur management strategies (Julian and Gu 2014, Julian et al. 2015); d) the mercury-related end products of these complexities must be predictable and quantified before an effective control or management strategy can be considered; and e) it is uncertain that reduction of sulfur inputs can reduce mercury methylation or shift methylation hot spots in the Everglades landscape or on a regional scale (Julian et al. 2016b). 57 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Carbon Carbon is present in all plants and animals. Organic matter refers to the carbon-based compounds found in terrestrial and aquatic environments, and it comes from the remains of plants and animals. Dissolved Organic Matter (DOM) in Everglades water plays a role in many biogeochemical and ecological processes, such as light adsorption, precipitation of minerals, and transport and reactivity of metals such as mercury. Factors that control DOM include peat drainage and oxidation, water movement and management, marsh drydown and re-wetting, and microbial and geochemical processes. DOM sources include peat, vegetation, periphyton, and detritus. A portion of DOM is dissolved organic carbon (DOC), the form that is most commonly measured (Aiken et al. 2011, Graham 2019). One of the reasons that carbon is abundant in the Everglades is because of the extensive peat soils that are as much as 90% organic matter (Figures 52, 56). During 1993 to 1996, REMAP documented spatial gradients in surface water organic carbon in canals and in the marsh, with the highest concentrations (40-69 mg/L) observed in canals within the EAA (Figure 45, Scheidt et al. 2000). An examination of the DOM characteristics of 2005 wet season REMAP samples indicates that the origin of this carbon is the peat soils of the EAA, with export in stormwater in canals due to flood control pumping and subsequent transport into the Everglades marsh (Yamashita et al. 2010). During 1974, when water from the EAA was backpumped into Lake Okeechobee, surface water Total Organic Carbon (TOC) was about 90-106 mg/L in canais within the EAA, with a decreasing gradient with distance into the lake such that TOC decreased to about 20-50 mg/L toward the interior (Brezonik and Federico 1975). During 2014, TOC and DOC (DOC, 0.45 pM filter) were measured in REMAP surface water samples. Essentially all of the carbon was in the dissolved form, with the average ratio = 1.00 for DOC/TOC (n=116). During 2014, as in 2005, REMAP data show that DOC in the Everglades exhibited a spatial gradient with the highest concentrations in WCA2 and the Refuge at locations near the EAA (Figures 46 and 47). The lowest DOC concentrations of 11 mg/L or less were all found in the Park in areas with marl soils of low organic content (about 20%, Figures 56 and 57). These spatial gradients throughout the Everglades are consistent with those reported by other investigators (Aiken et al. 2011, Julian et al. 2017b). DOC in WCA2 was lower during the 2014 wet season sampling (median = 25.0 mg/L) compared to 2005 (median = 30.5 mg/L). DOC and TOC concentrations tend to be higher during the dry season (Figure 47) (Scheidt and Kalla 2007; Ding et al. 2014). near the EAA, and had higher carbon content. Figure 45. Surface water samples collected during canal sampling in the 1990s. Samples with more color were collected at locations within or 58 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 46. Surface water dissolved organic carbon (mg/L) during November 2005 (left) and September 2014 (right). Carbon also is of interest because it plays a complicated role in mercury cycling (Graham 2019). DOM binds mercury, affects mercury solubility, and can influence the availability of mercury to microbes that methylate mercury. Areas strongly influenced by EAA stormwater have higher DOM concentrations and these areas have been reported to be more reactive with mercury than more pristine areas of the Everglades (Aiken et al. 2006). For STA2, from 2000-2011 an association was found between several water quality parameters including DOC at the STA inflow, and total mercury and methyl mercury at the STA outflow (Zheng et al. 2013). Using 2005 REMAP data, DOC was found to be associated with distribution of total and methyl mercury throughout the Everglades (Liu et al. 2009). A significant negative association was found between DOC and the mercury biomagnification factor between water and mosquitofish (Liu et al. 2008b). in the Everglades, DOM has been shown to enhance mercury methylation under low sulfide conditions (Graham et al. 2013). REMAP data from 1995-2005 were used to explore the relationship between mercury in prey fish (mosquitofish) and many biogeochemical analytes by using structural equation modeling. The modelling 59 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 1. indicated that DOC had a significant negative direct effect on mosquitofish mercury, yet DOC also had a direct positive significant effect on methylmercury in surface water and in periphyton, each of which had a direct positive significant effect on mosquitofish mercury (Pollman 2014). In contrast, during the November2005 REMAP sampling, DOC had a significant negative association with mercury biomagnification factor (Scheidt and Kalla 2007). From 1994 to 1999, STA1W exhibited no net removal of carbon, and about 93% of the surface water TOC was in the dissolved fraction (Gu et al. 2006). Water management and primary productivity processes in the marsh are important drivers controlling DOM dynamics (Yamashita et al. 2010), and water management is a key driver of DOC flux from Shark River Slough in the Park to the coastal estuary (Regier et al. 2016). A focus of restoration is increasing the flow of water into areas of the Everglades and coasts without adversely impacting these areas with water of a lesser quality (Aiken et al. 2011, Graham 2019). Wet/Dry season: ^ Dry ^ wet Figure 47. Wet and dry season surface water total organic carbon (mg/L) by subarea in the marsh during all REMAP sampling events. 60 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Dissolved Oxygen The concentration of dissolved oxygen (DO) within the water column is an important indicator of water quality. Aquatic life such as fish depend on oxygen for survival. A common default minimum DO water quality criterion for surface water is 5.0 mg/L. During photosynthesis, plants use the sun's energy to convert carbon dioxide and water into oxygen and cellular material (growth). During respiration, animals and algae remove oxygen from the water and use it to produce energy, releasing carbon dioxide and water as by-products. In natural wet prairie and open water slough habitats of the Everglades, plants such as bladderwort and associated algal communities produce oxygen throughout the day during photosynthesis. During daylight, there is more photosynthesis than respiration, and DO levels increase. During night, respiration exceeds photosynthesis and DO decreases. Natural background areas of the Everglades exhibit a strong daily fluctuation in DO, and at a single location DO can range from 1 mg/L. around sunrise to 9 mg/L around sunset. DO may not reach 5.0 mg/L until the afternoon. Dense emergent aquatic vegetation such as sawgrass and cattail contribute little oxygen to the water. In contrast, open water sloughs produce a surplus of water DO, which can flow into adjacent sawgrass habitats. DO concentrations in sawgrass marsh are lowerthan in open water, but oxygen remains throughout the night. In contrast, phosphorus-enriched cattail areas become void of DO during the night (Belanger et al. 1989). At phosphorus-enriched locations DO may never exceed 1 mg/L (McCormick et al. 1997, McCormick and Laing 2003). WCA2A, which has the largest area of phosphorus enrichment and impact, also has the highest sulfate concentrations and most extensive area of sulfate enrichment (Figures 40-41). DO concentrations also can be affected by elevated sulfate, in that sulfate loading increases microbial sulfate reduction in soils. This leads to reducing conditions, increased sulfide, *| 8.. MI--: and lower water column DO concentrations (Orem et al. 2011). Figure 48. Daily water column dissolved oxygen (DO) at background marsh locations in the Everglades. DO peaks late in the afternoon at concentrations higher than 100% saturation, with the lowest concentrations occurring shortly after sunrise (from Weaver 2004). Because of the natural daily variation in DO in the Everglades due to photosynthesis and respiration, in 2014 FDEP adopted a Site-Specific Alternative Criterion (SSAC) for DO to replace the 5.0 mg/L minimum requirement. Data from marsh reference sites with low TP and away from canals and groundwater influence were used to define the natural background levels and variation in DO (Figure 48, Weaver 2004). DO concentration varies naturally not only with time, but also with water temperature; as water temperature increases, the water is able to 61 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) WET PRAIRIE SAWGRASS CATTAIL SLOUGH WET PRAIRIE SAWGRASS CATTAIL SLOUGH Figure 49. REMAP water column dissolved oxygen (DO) percent saturation by habitat (left), and DO concentration minus the SSAC by habitat (right). hold less oxygen (lower saturation capacity). The DO concentration expected for a sampling event at a station is calculated from the water temperature and sampling time. The annual SSAC for the station is assessed based on comparing the annual average measured DO concentration and the average of the corresponding DO SSAC limits. In 2014 the annual minimum DO concentration required by the SSAC varied from about 1.0 mg/Lto 4.0 mg/L depending upon location (Julian 2016c). As expected, phosphorus-impacted marsh sites had lower DO concentrations and they usually did not meet the SSAC (Julian et al. 2016a). REMAP DO data are shown in Figure 49 by habitat at 833 marsh locations sampled from 1995 to 2014. DO concentration ranged from 0.3-13.6 mg/L, and from 3.5-189% saturation. Cattail had the lowest DO percent saturation, and was the only habitat that never exceeded 100% (Figure 49, left). Because REMAP data have only single measurements at a location, it is not appropriate to apply the DO SSAC to REMAP data for regulatory purposes. However, for informational purposes only, the expected DO concentration (SSAC) for the particular temperature and time that DO was measured was subtracted from the observed DO concentration (Figure 49, right). Once again, cattail was the habitat with the lowest DO concentration, and this difference was very highly significant (F = 15.52, p < .000001). DO in cattail was lower than the SSAC for most of the data. 62 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) pH The pH of water or soil can influence many water quality and chemical processes. The pH is defined as the logarithm (base 10) of the reciprocal of hydrogen ion activity or concentration. The pH scale ranges from 0 to 14, with pure water having a pH of 7.00, or neutral. Increased hydrogen ion activity lowers the pH toward acidity, while decreased activity increases the pH toward becoming basic. The pH of natural waters is usually between 6.0 and 8.5 (Hem 1985). Surface water pH and soil pH varied spatially during September 2014, with the lowest values found in the Refuge (Figures 50-51). Rainwater in the Everglades had a precipitation-weighted mean pH of 5.1 for 2014 (NADP 2014a). The soft (low ionic content) water of the Refuge has low capacity to buffer against acidity (annual median alkalinities at interior locations as low as 8 mg as calcium carbonate per liter), while the hard (high ionic content) waters of the Park have high buffering capacity (annual median alkalinities of about 200 mg as calcium carbonate per liter) (Weaver et al. 2007). Its low pH and soft water conditions make the Refuge vulnerable to mineral enrichment from groundwater or surface water (McCormick et al 2011). For 2014, the median surface water pH by subarea was 6.11 for the Refuge, 7.23 for WCA3, 7.46 for WCA2 and 7.52 for the Park. Photosynthesis by aquatic organisms removes carbon dioxide from the water column during daylight hours, resulting in an increase in surface water pH (Hem 1985, Gleason and Spackman 1974). For example, in a natural wet prairie marsh in the Park, dominated by spikerush and bladderwort 63 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) plants and an extensive calcareous periphyton mat, the pH at a single location fluctuated over 24 hours from 7.1 at midnight to 8.5 late in the afternoon while the corresponding dissolved oxygen (DO) went from about 1.0 mg/Lto 10.0 mg/L (Scheidt et al. 1985). The REMAP pH values >8 all occurred in the Park after 1400 hours. Given that the September 2014 REMAP measurements of in-situ water pH occurred between 0830 and 1730 hours, and that REMAP sampling took place from south to north over a 17-day period, the observed spatial pattern in pH cannot simply be explained by diurnal fluctuations. Florida's water quality criterion for pH generally requires that pH shall be between 6.0 and 8.5 and shall not vary more than one unit above or below natural background (Chapter 62-302.530(52)(c) F. A. C.). From 1995 to 2014, REMAP had only 17 of 951 (2%) pH measurements that were less than 6.0, all of which were in the interior of the Refuge. Florida reports pH values <6.0 within the Refuge at a similar frequency. These excursions below the criterion are viewed as a consequence of the Refuge's naturally low alkalinity and are not of ecological concern (Julian et al. 2016a). The lowest soil pH (median = 6.72) was within the Refuge interior, which has highly organic peat soils. The highest soil pH (median = 7.74) generally occurred within the mineral enriched soils in WCA2A, and in the Park's marl soils, which contain more calcium than peat soils. The marl soil found throughout much of the Park (Figures 56-57) contributes to buffering capacity and results in these higher pH values. Soil calcium, which contributes to higher soil pH, is five times higher in the soils within the Park's eastern marl prairies (227 g/kg) than in the peat soils of Shark River Slough (54 g/kg) (Osborne et al. 2011 b). WCA3N WCA3S Wet/Dry season: Dry S Wet Figure 51. Wet and dry season marsh surface water pH (Standard Units) by subarea during all REMAP sampling events. 64 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) SOILS and SOIL SUBSIDENCE Soil is a defining characteristic of an ecosystem, and soil preservation is important for ecosystem protection and restoration (Scheidt and Kalla 2007, Nungesseret al. 2014, Orem et al. 2014b). The Comprehensive Everglades Restoration Plan has adopted objectives, performance measures, and performance targets in order to define restoration goals, track ecosystem status, and measure restoration effectiveness. Among these is restoring the natural rates of organic soil and marl soil accretion in the Everglades and stopping soil subsidence (RECOVER 2007). There are two major soil types in the Everglades (Figure 52). The wetland soils of the central Everglades are primarily peat formed by slowly decaying plant matter in areas that are flooded much of the year. The other major soil type is calcitic mud or marl commonly found in the shallower peripheral marshes of the Park that have shorter periods of surface water inundation. Marl is found in association with thick, calcitic algal mats (periphyton) (Figure 58, Eastern Marl Marsh), which precipitate calcium carbonate from the water column in these hard water areas (Gleason and Stone 1993). Figure 52. Everglades peat (top, organic matter 78%, WCA2A) arid marl (bottom, organic matter 14%, the Park). Bottom photo also shows a benthlc periphyton mat overlaying the soil core. Historically, the Everglades contained the largest body of organic peat soils in the world, covering over 3,000 square miles and accumulating to a thickness of up to 17 feet in what would become the EAA (Stephens and Johnson 1951). The origin and perpetuation of peat and marl soils are dependent upon water depth, the duration of surface water inundation, and the resulting wetland vegetative communities. Shortened surface water inundation can cause soil loss or changes in soil composition, which may in turn result in altered vegetative communities. These altered plant communities in turn may cause further changes in soil type and thickness as this different plant community eventually decomposes and forms altered soil. Some soil cores collected by REMAP have alternating peat and marl layers within the 0-10 cm profile. Soil loss within the Everglades is largely due to water management practices during the last 100 years. The major canals draining the EAA were completed by 1917 and extend southeast through the Everglades to the Atlantic Ocean. However, unimpeded surface water flow from the EAA southward through the Everglades to the Park, Florida Bay, and the Gulf of Mexico still occurred until the late 1950s, when levees were constructed forming the southern boundary of the EAA. During the early 1960s, additional levees were completed that 65 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) compartmentalized the Everglades into the Water Conservation Areas. By the 1960s, Everglades surface water depths, flow, and inundation periods had been greatly altered (Light and Dineen 1994). Peat soils are subject to subsidence and loss of surface elevation when drained. Oxidation, burning and compaction are considered the dominant subsidence forces. An inch of Everglades peat that takes a century to form can be lost within a few years, or within hours if dry soils are subjected to fire. In the early 1900s, the deep peat soils (mostly formed by decaying sawgrass) of the 700,000-acre EAA were drained to facilitate agricultural production. The process of soil formation was reversed in 1906 when the first drainage canals were cut from Lake Okeechobee through the EAA to the coast (Stephens 1956). Subsequent subsidence within the EAA and efforts to control it on agricultural lands are well documented (Shih et al. 1998, Ingebritsen et al. 1999, Snyder 2005). In 1912, much of the EAA had soils thicker than 10 feet, or 120 inches (Stephens and Johnson 1951, Stephens 1956). By 1988, only 18% of the EAA had soil thicker than 51 inches, while 53% ofthe area had soils less than 36 inches thick, and 11% had soils less than 20 inches thick (Coxet al. 1988). By 2050, under current agricultural practices, up to 93% ofthe EAA is projected to have soils less than 36 inches thick, with 82% less than 20 inches thick and 53% less than 8 inches thick. The decrease in soil volume from 1988 to 2050 is projected to be 57% or 1.2 x 109 m3 (Snyder 2005). Geospatial techniques have been used to estimate that during the 1900s about two-thirds ofthe peat volume within the EAA (4.5- 4.9 x 109 m3) was lost due to subsidence (Aich et al. 2013). Within the EAA, production of agricultural crops such as vegetables and the more prevalent varieties of sugarcane require that the water table be maintained below the ground surface. The ground surface of the EAA basin, which historically was sawgrass marsh that flooded most ofthe year, is now an average of 2 meters lowerthan it was in 1910 due to subsidence (Aich etal. 2013). Frequent wet season rain events necessitate recurring pumping in orderto maintain the watertable below the ground surface, which continues to subside further. Each flood control pumping event has the potential to leach and export soil constituents, such as phosphorus, nitrogen, sulfur and carbon, in the water pumped southward to the Everglades. Agricultural Best Management Practices (BMPs) are directed at phosphorus removal (ludicello et al. 2016). The STAs are more effective at removing phosphorus than nitrogen, sulfur or carbon. Projections that one-half of the EAA may have less than 8 inches of soil by 2050 bring into question the viability of agriculture with current practices (Snyder 2005). Conversion from agriculture to residential land use could result in the need to export greater volumes of stormwater and its constituents to the Everglades, if residential land use requires that the water table be maintained at even lower levels to provide flood control. The rate of soil subsidence slowed from 1.12 inches per year between 1924 and 1967 to 0.55 inches per year between 1967 and 2009. Implementation of BMPs since the 1990s has led to more water storage on EAA fields, which has slowed subsidence (Wright and Snyder 2016). Raising water levels and increasing biomass also would slow subsidence (Rodriguez et al. 2020), as would using agricultural lands for water storage (Ouellette et al. 2018). During the 1995-96, 1999, 2005, and 2014 REMAP sampling events, a metal probe was inserted to the point of refusal to measure soil thickness (0.00 to 12.00 ± 0.05 feet) at 976 locations (Figure 53). REMAP is the only source of directly measured soil thickness data throughout the Everglades post-1940s. Comparisons 66 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 53. Soil thickness at 9/6 locations measured by REMAP for all sampling events from 1995 to 2014. The iriset shows soil thickness as reported in 1946 (Davis). 67 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) From the 1940s to the 1990s, over one-half of the soil was lost from drier portions of the Everglades. Longer periods of surface water inundation would help to maintain marsh soils and the plant communities and wildlife habitat of these wetlands. of soil thicknesses measured by REMAP in 1995-96 to those reported by Davis in 1946 indicate that the drier short hydroperiod portions of the Everglades, such as WCA3 north of I-75 (Figure 21), lost 39% to 65% (2.0 to 6.0 x 108 m3) of their soil. Soil thicknesses of 3 to 5 feet in the 1940s had diminished to only 1 to 3 feet by 1995-96, with less than 1 foot remaining in some areas. WCA3B and the Northeast Shark Slough portion of the Park were found to have lost up to 3 feet of soil, representing a 42% and 53% loss of volume, respectively. These three portions of the Everglades, about 200,000 acres, have been subjected to decreased surface water inundation since completion of the levees that surround the WCAs in the 1950s and 1960s. From the 1940s to 1990s, the entire Everglades Protection Area lost up to 28% of its soil volume due to oxidation and subsidence (Scheidt et al. 2000). Comparison of REMAP data from 2014 to 2005 do not indicate any further subsidence of soils landscape- wide (Figure 54, Wald Ftest p=0.68). Soil thickness data from 1995-2014 were combined and the study area was divided into 10 subareas (Figure 55). The WCAs were divided into five subareas: the Refuge, WCA2, WCA3A north of I-75, WCA3A south of I-75, and WCA3B. The Park also was divided into five subareas: Shark RiverSlough, Northeast Shark Slough, Eastern Marl Marsh, Taylor Slough, and Ochopee Marl Marsh. The deepest soils are the peat deposits within the Refuge, with a median soil thickness of 8.6 feet (Table 4, Figures 53 and 57). The maximum soil thickness measurable with the probes was 12.00 feet, so deeper peat in the Refuge is possible. The deepest median soil thicknesses for the other portions of the Everglades were 4.2 feet in WCA2 and about 3 feet in WCA3A south of I-75 and WCA3B, areas which typically stay inundated year-round. Median soil thickness in Shark River Slough was 1.5 feet, and the marl marshes in the Park that dry out more frequently than other parts of the Everglades had median soil thicknesses of less than 1.0 feet. The overall median soil thickness for the Everglades was 2.3 feet. Figure 54. Soil thickness (feet) estimates of marsh area for 2005 and 2014. Soil thickness was not different (Wald F test, p=0.68). 68 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 0 2.5 5 10 Miles 1 i i i I i i i I Canals rj Stormwater Treatment ' Areas _ Miccosukee Tribe - I I Everglades Federal Reservation Subareas o LNWR o WCA2 o WCA3N Q WCA3B O WCA3S O ENP SHARK SLOUGH o ENP NORTHEAST SHARK SLOUGH © ENP OCHOPEE MARL o ENP EASTERN MARL o ENP TAYLOR SLOUGH Figure 55. Ten subareas used for soil thickness and volume calculations. Soil volume in the freshwater Everglades was calculated by subarea as the area times the median soil thickness. Summing the subareas results in an overall soil volume of 4.7 x 10" m3 (Table 4). Other investigators used historical data sets to calculate soil volume by subtracting the current ground surface elevation from the bedrock surface elevation. They also obtained a soil volume estimate of 4.7 x 109 m3 although some of 69 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) 2 Area (km ) Count Median Soil Thickness (feet) Median Soil Thickness (m) 3 Volume (m ) LNWR 567.0 104 8.57 2.61 1.49 x 109 WCA2 539.0 104 4.15 1.26 0.69 x 109 WCA3AN 715.8 129 1.60 0.49 0.35 x 109 WCA3AS 1288.0 224 2.90 0.88 1.14 x 109 WCA3B 401.4 78 3.25 0.99 0.40 x 109 ENP Ochopee Marl Marsh (OMM) 437.6 79 0.83 0.25 0.11 x 109 ENP Shark River Slough (SRS) 357.6 67 1.50 0.46 0.16 x 109 ENP Northeast Shark Slough (NESS) 251.1 44 1.64 0.50 0.12 x 109 ENP Eastern Marl Marsh (EMM) 830.5 134 0.70 0.21 0.18 xlO9 ENP Taylor Slough (TS) 59.1 13 3.00 0.91 0.05 x 109 TOTAL 5447.1 976 2.30 0.70 4.69 x lo' Table 4. Soil thickness and volume by 10 subareas, REMAP data 1995 to 2014. the subarea volumes differed. Today the Everglades covers about one-half of the area than it did historically, but less than one-quarter of the peat remains. One-half of the carbon in the peat soils has been lost (Hohner and Dreschel 2015, Dreschel et al. 2017). This peat loss is a concern for Everglades restoration (Nungesser et al. 2014). Managing water is essential to preventing overdrainage and further loss of peat soils. The soil organic matter content observed by REMAP from 1995 to 2014 ranged from 1 % to 100% (Figures 56 to 57), with a median of 80%. Peat soils are highly organic, while marl soils are primarily mineral. The highest organic matter content was found in the thick peat soils of the Refuge, having a median of 94%. WCA2A and WCA3 south of I-75 also had soils exceeding 75% organic matter. These highly organic zones coincide with the longer hydroperiod portions of the Everglades. The area of maximum soil loss within WCA3A north of I-75 had a median soil organic matter content of 75%, the lowest in the Water Conservation Areas. The peat soils in the Shark River Slough trough of the Park had a median organic matter of 78%, in contrast to the marl soils of the Park which had a median of only 27%. Soil bulk density, the mass of dry soil per unit of bulk volume, ranged from 0.04 to 0.90 g/cc (Figures 56 to 57). The highly organic peat soils of the Refuge had the lowest bulk density, with a median of 0.06 g/ cc, in contrast to the marl soils of the Park which had a median of 0.36 g/cc. The median soil bulk density for WCA3A north of I-75 was 0.17 g/cc, the highest in the WCAs. Within the WCAs, this portion of northern WCA3A had the lowest organic matter content, the highest bulk density, and the greatest soil loss. All of these observations are suggestive of formerly deeper peat soils of the 1940s being subjected to drier conditions due to water management changes. Surface water inundation has been decreased, and consequently soils have subsided and become less organic due to increased biochemical oxidation and more frequent wildfires. A focus of Everglades restoration and the Central Everglades Planning Project (CEPP) is restoring more natural water flow, depth, and duration into and within WCA3 and the Park by: increasing storage, 70 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 56. Soil percent organic matter (left) and bulk density (g/cc) (right), 0-10 cm profile. Soils in the Park with organic matter < 30% typically are marl, and organic matter > 50% indicates peat soils. treatment and conveyance of water south of Lake Okeechobee; removing canals and levees within the central Everglades; and retaining waterwithin the Park while protecting urban and agricultural areas to the east from flooding (USACE 2014). It is expected that returning the water flows of the central Everglades to a more natural state will decrease peat subsidence, increase soil accretion and return the central Everglades to a net carbon sink (Richardson et al. 2014). Hydroperiod and nutrient availability drive plant litter decomposition and are important to peat formation and accumulation (Pisani et al. 2018). It is also essential to reduce nutrient inputs if natural peat accumulation is to be restored (Sklar et al. 2005a, Sklar et al. 2010, Osborne et al. 2017). Further decreases in water depth and inundation periods would alter the Everglades ecosystem through drought, peat loss and carbon emissions, wildfires, loss of the unique ridge and slough patterns, shifts in plant and animal communities, and spread of exotic species. Adaptive restoration planning incorporates climatic and environmental uncertainties into long-term ecosystem restoration plans, structural design, and management (Nungesser et al. 2014, Flower et al. 2019). 71 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Š Median ~ 25%-75% X Non-Outlier Range o Outliers * Extremes 1 O 40 LNWR WCA 3AN WCA3B NESS EMM WCA2 WCA3AS SRS OM TS Š Median ~ 25%-75% X Non-Outlier Range o Outliers * Extremes 11 1 ? WCA3AN WCA3B NESS EMM WCA 2 WCA3AS SRS OM TS Š Median ~ 25%-75% I Non-Outlier Range o Outliers * Extremes $ * I | o o i S :» 7 T ' 1 LNWR WCA 3AN WCA 3B NESS EMM WCA 2 WCA3AS ENP SRS OM TS Figure 57. Soil thickness (feet, top), percent organic matter (middle) and bulk density (g/cc,bottom) by 10 subareas. 72 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Loxahatchee National Wildlife Refuge Peat soil Soil Thickness > 8.6 ft (median, n=104) Bulk density = 0.06 g/cc Organic matter = 94% Water depth (1994-2014)* = 1.0 ft* Days Dry = 2%* * water level gage 1-7, data from EDEN WCA 3A South Peat soil Soil Thickness = 2.9 ft (median, n=224) Bulk density = 0.11 g/cc Organic matter = 88% Water depth (1994-2014)* = 2.2 ft Days Dry = 0%* * water level gage 3-65, data from EDEN Figure 58. Contrasting soil and vegetation characteristics of selected locations in the Refuge (top), WCA3A south (bottom) and the Eastern Marl Marsh in the Park (next page). 73 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Eastern Marl Marsh - Rocky Glades Marl soil Soil Thickness = 0.7 ft (median, n=134) Bulk density = 0,45 g/cc Organic matter = 21% Water depth (1994-2014)* = 0.14 ft Days dry* = 42% * water level gage NP206, data from EDEN Figure 59. Interstate 75 at the eastern edge of the Everglades looking west. Northern WCA3A is to the right. Shortened periods of surface water inundation in northern WCA3A have led to soil loss since the 1940s. 74 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) MACROPHYTES and PERIPHYTON Figure 60. Mosaic of tree islands, sawgrass marsh arid wet prairies within Shark River Slough, Everglades National Park. The brownish color is the periphyton mat at the water surface in wet prairies. This photo was taken during the wet season when water depths were about 3 feet. 5lant Communities The Everglades are defined by a unique mosaic of vegetation community types such as tree islands, sloughs, wet prairies and sawgrass marshes (Figures 1, 2, 60). Factors driving vegetation community composition include hydroperiod, water depth, water velocity, nutrients, invasive plant species that are not native, disease, and disturbances such as fire, frosts, and hurricanes. Wet prairies and open water sloughs without dense plants growing out of the water (emergent macrophytes) serve as preferred habitats for foraging wading birds (Bancroft et al. 1992). These areas are the marsh habitats with the greatest diversity of native flora and fauna (Gunderson and Loftus 1994). Classified vegetation maps were produced for 1 km2 areas centered on each REMAP station during 1999 (Stober et al. 2001b, Welch and Madden 2000) and 2005 (Richards and Philippi 2005). In addition, during 1999 and 2005 plant species frequencies were determined along transects (Stober et al. 2001b, Richards and Philippi 2005). In 2005 the presence and distribution of non-native or invasive plant species were determined using several methods (Richards and Philippi 2005). During the September 2014 REMAP event there were several efforts to document vegetation. World View-2 satellite imagery of a 1 km2 area centered on 65 REMAP sampling stations was used to create classified vegetation community maps. Vegetation mapping provides a landscape context for REMAP biogeochemical and biotic information. Standing stocks of carbon, nitrogen and phosphorus in sawgrass, water, soil, floe, and periphyton were estimated (Richards et al. 2017). Digital photographs document plant communities during September 2014 at all 118 sampling locations: a ground view of the area sampled to the left of the helicopter, 75 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) TOTAL PHOSPHORUS SOIL q Cattail present Figure 61. Cattail presence during September 2014. nine panoramic photos at 45-degree increments, each of the three soil cores, and an aerial view at 100-200 feet. For each station, photodocumentation, classified vegetation maps, and biogeochemical data from 2014 and 2005 are available at: http://diair.fiu.edu/ amaps/EverREMAP.Dhp. Field crews also recorded the dominant plant community at each sample point based on visual observation. In addition, vegetation was sampled at a subset of stations for chemical analyses: sawgrass leaf clippings at 60 stations for carbon, nitrogen, and phosphorus; and whole sawgrass plants at 27 stations for total mercury and methylmercury. Figure 62. Soil total phosphorus (mg/kg) and cattail presence based on REMAP data 1995 to 2014. mg/kg 1000 800 700 600 500 400 300 Cattail is a native species that can respond to phosphorus enrichment and replace sawgrass, water lily in sloughs, and wet prairie plants such as spikerush and bladderwort. Cattail expansion was one of the first visual consequences of P enrichment observed during the 1970s (Davis 1994). The rate of cattail expansion in WCA2 from 2003 to 2011 was only 20% of the rate of expansion from 1996 to 2003 (Zweig and Newman 2015). Conversion of wet prairies with open waterto dense cattail constitutes a loss of the preferred foraging habitat for wading birds (Turner etal. 1999). There is an association between cattail presence and elevated soil phosphorus or proximity to canals (Figures 61-62, Scheidt and Kalla 2007). 76 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Low phosphorus conditions are required for natural Everglades periphyton and plant communities to be maintained. Periphyton Periphyton is a complex mixture of algae, bacteria, microbes, and detritus that is attached to submerged surfaces. Periphyton can be found in abundance throughout the natural Everglades ecosystem, and includes: loose, flocculent aggregations in the soft-water Refuge; thick, calcareous mats in the central Everglades sloughs; and benthic mats in the Park's marl prairies (Gaiser et al. 2011). Well-developed, attached or floating calcareous periphyton mats are a conspicuous and defining characteristic of the hard-water Everglades, particularly in wet prairies and deeper slough areas (Figures 63-64). The natural periphyton mats across mineral-rich portions of the Everglades (see section on specific conductance) are dominated by calcium-precipitating (calcareous) cyanobacteria and have a high calcium carbonate content, while those in the soft-water Refuge are largely organic (non-calcareous) (McCormick et al. 2011). These assemblages of microscopic plants serve multiple functions such as providing oxygen to the water column for fish, removing calcium carbonate from the water and depositing it as soil, removing phosphorus from the water to very low concentrations, and serving as the base of the food web (McCormick et al. 1999). Figure 63. Presence of floating, benthic, and attached (epiphytic) caicitic periphyton communities during September 2014. Hydroperiod and water depth, water ions, and phosphorus concentration all affect periphyton extent and community structure (Browder et al. 1993). Periphyton communities are sensitive to very slight increases in nutrient concentrations, with increases in phosphorus causing changes to the periphyton assemblage, including species composition and biomass, or even the disappearance of the entire mat. Consequently, periphyton is a sensitive and important indicator of Everglades marsh ecosystem status (Gaiser 2009, Gaiser et al. 2004, McCormick et al. 2002, Surratt et al. 2012). In the Refuge, shifts in periphyton community Figure 64. Types of tverglades caicitic periphyton communities sampled: epiphytic - attached to plants such as bladderwort (left) or 'sweaters' attached to spike rush shown on a gloved hand (middle); and benthic at the soil surface within a 0.5 meter wide quadrat (right). 77 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) composition and function provide an especially sensitive indicator of mineral enrichment (McCormick et al. 2011). Monitoring and species-based approaches provide early warning signs of environmental change that are only later identified using traditional water quality approaches, and that have the potential to be propagated into higher trophic levels. Changes in periphyton biomass, nutrient content and composition can have cascading effects upon the Everglades food web (Gaiser et al. 2005, 2015). During 2014 REMAP documented three types of easily-seen calcareous periphyton growth forms (Figures 63-64): periphyton floating on the surface or within the water column; epiphytic periphyton or 'sweaters' attached to plants; and benthic periphyton as a discrete layer at the soil surface. The most common form of periphyton was epiphytic, which was observed at 53% of the stations, followed by benthic (31%) and floating (23%). Benthic mats were most common in the marl, short-hydroperiod portions of the Park. Percent periphyton cover, which was documented within a randomly located 0.25 m2 frame, ranged from 0-100%. There were no periphyton mats or sweaters sampled in the soft-water Refuge due to the loose, fiocculent nature of the periphyton there and time constraints on REMAP field sampling crews. Generally, there was no periphyton observed or sampled in most of northern WCA3A, an area with shortened hydroperiod and higher risk of fire. Floe FLOC THICKNESS A conspicuous layer of fiocculent matter (floe) is present at the sediment-water interface in much of the Everglades (Figures 65-66). Floe in the Everglades is believed to consist of an assembly of plant detritus, periphyton, carbonates, and other remains of aquatic organisms (Neto et al. 2006). Floe plays a role in nutrient and carbon cycling and resuspension. Detrital remains of plants such as bladderwort (Utricularia species) comprise the primary components of floe in deep sloughs in WCA3 (Troxlerand Richards 2009). In 2014, REMAP floe thickness ranged from 0-19 cm, with the thicker layers generally occurring at the longer hydroperiod locations. Mass balance estimates using REMAP 2014 data show that fioc contains a significant pool of carbon, nitrogen and phosphorus (Richards et al. 2017). The TP content of fioc is a sensitive indicator of enrichment at marsh transects downstream of structures discharging into the Refuge, WCA2, WCA3A within the Miccosukee Tribe's Federal Reservation, and Taylor Slough in the Park (Wright et al. 2009). 0 to 1 A 1 to 3 ~ 3 to 6 ~ 6 to 19 Figure 65. Floe thickness (cm) during September 2014. Figure 86. Peat soil and an extensive floe layer at the surface water - soil interface. 78 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) MERCURY Mercury contamination of gamefish and wildlife in the Everglades has been a management concern and focus of research and monitoring since the 1990s. The primary sources are atmospheric. Mercury falls on the Everglades as gaseous divalent mercury (Hgll). Biogeochemical processes within the Everglades marsh convert this mercury to methylmercury, a toxic form. Methylmercury accumulates (bioaccumulates) in the tissue of aquatic organisms in the Everglades. Mercury also increases in concentration (biomagnifies) the higher an animal is in the food chain, and at higher concentrations it presents a risk to wildlife and humans through fish consumption. During the early 1970s, mercury was documented in birds and gamefish in Everglades National Park, and mercury was identified as a potential ecological concern (Ogden et al. 1974). In the late 1980s, unexpectedly high levels of mercury were found in gamefish and two dead Florida panthers in the Park (Roelke et al. 1991). The mercury sources were unknown. In order to protect human health, in 1989 Florida issued a consumption advisory either restricting or recommending no consumption of gamefish such as largemouth bass from the Everglades (FDHRS 1989) (Figure 67). Consumption advisories have remained in place throughout the Everglades since, with the current advisory recommending that women of child-bearing age and young children should not eat largemouth bass, and should limit consumption of 17 other fish species (FDOH 2020). The existence of these advisories means that the Everglades waterbody does not meet the "fishable" portion of its designated use underthe Clean Water Act. In addition, ecological risk assessments and mercury dosing studies have indicated that populations of top predators in the Everglades could be adversely affected by mercury contamination in that mercury accumulation through the food web has the potential to reduce the health or breeding success of wading birds (Rumbold 2005, 2019b; Rumbold et al. 2008, Spalding et al. 2000, Duvall and Barron 2000, Zabala et al. 2020) and the Florida panther (Barron et al. 2004, Rumbold 2019b). In the early 1990s, the South Florida Mercury Science Program was established. This program was a collaborative effort by FDEP, USEPA, SFWMD, Florida Game and Freshwater Fish Commission and the USGS. The program was designed to determine: the Figure 67. Young fisherman and fish consumption advisory to protect human health at the main Refuge boat ramp. Depending upon the location in the Everglades, Florida's present advisory recommends that women of child-bearing age and young children should not eat largemouth bass, and should limit consumption of 17 other fish species (FDOH 2020). The Florida Department of Health and Rehabilitative Services has issued a health advisory urging limited consumption of largemouth bass and warmouth caught in certain portions of the Everglades due to excessive accumulation at the element mercury. Fish caught in Arthur R. Marshall Loxahatchee National Wildlife Refuge (Water Conservation Area 1} should not be eaten more than once per week by adults and not more than once per month by children under 15 and pregnant women. Fish caught in Water Conservation Areas 2a and 3 should not be eaten at all. For additional information, contact the Florida Department of Health and Rehabilitative Services at (405) 355-3018. 79 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) sources of mercury; potential risks to humans and wildlife; how mercury enters the aquatic food chain and bioaccumulates; how inorganic mercury is transformed into methylmercury; how mercury cycles through air, water and soil; and what actions could be taken to decrease mercury levels in wildlife and gamefish (FDEP1997). Potential mercury sources were: stormwater pumped into the Everglades from the Everglades Agricultural Area; release from Everglades soils into waterduring wet season rewetting; burning of sugarcane during harvesting; medical and municipal waste incinerators; and other sources of emissions into the atmosphere. It was determined that atmospheric deposition was by far the predominant mercury source to the Everglades, originating from both within and outside of Florida (FDEP 2013). Florida's class III surface water criterion for total mercury is 12.0 nanograms per liter (ng/L or parts per trillion). Since 1995, including all sampling events, REMAP has sampled 851 locations within the Everglades marsh for total mercury in surface water, and only 6 samples exceeded 12.0 ng/L. All 6 samples were collected during the dry season (1990s and 2005) at shallow marsh sites (water depths from 0.1 to 0.7 feet). Biomagnification and bioaccumulation of mercury to unacceptable levels in gamefish occurs even though surface water concentrations are below the 12 ng/L criterion. Divalent mercury in rainfall that is deposited into surface water can be converted to methylmercury (MeHg) by bacteria in the presence of sulfate and organic carbon (Orem et al. 2011, Aiken et al. 2011). Methylmercury is the form of mercury that bioaccumulates and biomagnifies in the aquatic food chain. Concentrations in largemouth bass are 10 million times higher than concentrations in surface water. There is no numeric water quality criterion for MeHg in surface water in the Everglades. Many factors affect the bioaccumulation of mercury in aquatic life, such as the length of the aquatic food chain, soil type, pH, and dissolved organic material (USEPA2001). During the last two decades, about 30 factors have been suggested by scientists as affecting mercury bioaccumulation in the Everglades. Interrelationships among the factors are complex and may be waterbody-specific. Because of these complexities, USEPA concluded that in orderto protect human health it is more appropriate to have a water quality criterion for MeHg based on fish tissue concentrations, rather than on water concentrations. The MeHg water quality criterion USEPA recommended in 2001 is a fish tissue residue criterion of less than 300 micrograms per kilogram (|jg/kg, or parts per billion), or 0.3 mg/ kg (USEPA 2001). Accordingly, Florida uses fish consumption advisories based on mercury in gamefish exceeding 0.3 mg/kg to determine that a waterbody is impaired (is not meeting its designated use, i.e., is not fishable) (FDEP 2013). In 2013 Florida adopted a statewide mercury Total Maximum Daily Load (TMDL) to establish the allowable loadings and needed reductions of mercury into Florida's fresh and marine waters that would restore these waterbodies so that the human health concern associated with the elevated mercury in fish tissue impairment will be addressed. The TMDL calls for an 86% reduction in mercury sources, which arrive in Florida waters predominantly by atmospheric deposition, and are from both within and outside Florida. Over 95% of wet and dry deposition of mercury to the Everglades originates outside of Florida (FDEP 2013), and it is estimated that 85% to 95% of the mercury deposited on the Everglades is from the long- range transport of mercury from sources outside of the United States (Vijayaraghavan and Pollman 2019). 80 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Previously, Florida's anthropogenic atmospheric mercury emissions were significantly reduced, from about 160,000 pounds per year in 1988 to 3,000 pounds in 2009, due to air pollution emission reductions required by the federal Clean Air Act and Florida's implementing rules (FDEP 2013). Mosquitofish (Figure 68) are a small preyfish (maximum length 1.5 inches, n=1841, REMAP 2005 and 2014) that have been sampled by REMAP since the 1990s because they are an ideal indicator of mercury contamination: a) they are the most abundant fish in the Everglades and are found throughout the canals and in all marsh habitats; b) they are easily sampled; c) they are a prey fish in the food web for gamefish and wading birds, so they provide insights for both human health and ecological health; and, d) because of their lifespan of only several months and a small home range, they integrate mercury exposure over a short time frame in a discrete area. During the five REMAP wet season sampling events, mosquitofish were collected at 94% of the 532 Everglades marsh sites, including wet prairie, sawgrass and cattail habitats. Everglades mosquitofish are a secondary consumer reported to be at trophic level 2.0 to 3.0 (Loftus et al. 1998) and 4.0 to 4.5 (Williams and Trexler 2006). They consume animal prey (crustaceans, insects, arachnids), algae, detritus and plant matter (Loftus et al. 1998). USEPA has recommended a mercury concentration of 77 |jg/kg at trophic level 3 for protection of birds and mammals (USEPA 1997), while the United States Fish and Wildlife Service has recommended a level of 100 jjg/kg in prey fish in order to protect top predators such as wading birds from mercury contamination (Eisler 1987). During 1995-96 and 2005, REMAP documented a pronounced spatial gradient in mosquitofish mercury, with the highest concentrations in remote portions of WCA3A and extending into Shark River Slough in the Park (Figure 69) (Stober et al. 1998, Stober et al. 2001b, Scheidt and Kalla 2007). These results are consistent with those for other biota indicating that the highest mercury in the Everglades occurs in the Park or WCA3 for largemouth bass and great egrets (Axelrad et al. 2007) and alligators (Rumbold et al. 2002). A risk assessment on the effects of MeHg on great egrets concluded that birds foraging in the Park have a high probability of exceeding both the daily and cumulative acceptable dose levels necessary to protect nestlings and pre-nesting females (Rumbold et al. 2008). A regional-scale ecological risk assessment and synthesis concluded that mercury remains a risk to a variety of birds and mammals in certain South Florida sub-regions or hotspots including portions of the Everglades (Rumbold 2019b). Mercury in mosquitofish was lower in 2014 than in 2005, which was lower than in 1995. This result is suggested by krigs, and confirmed by box-and-whisker plots and the Wald F test for differences between means in the CDF curves (Figures 69-71). In 2014, the proportion of the Everglades marsh that was above the 77 jjg/kg predator protection level was 13.0 ± 5.7% (Figure 71), as compared to 64.7 ± 7.3% in 2005 and 70.5 ± 7.1% in 1995. In 2014, 6.5 ± 4.2% of the marsh had mosquitofish mercury concentrations that Figure 68. Mosquitofish, such as most of the fish shown here, are the most abundant fish in the Everglades and are an ideal indicator of mercury. 81 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 69. Total mercury concentration in mosquitofish (|jg/kg) during Everglades REMAP wet season sampling events: September 1995 (left), November 2005 (middle) and September 2014 (right). Š Median ~ 25%-75% X Non-Outlier Range o Outliers * Extremes I Predator Protection Level is 77 ng/g o X Figure 70. Box and whisker plot of total mercury concentration in mosquitofish (ng/g) during Everglades REMAP wet season sampling events: September 1995 (left), November 2005 (middle) and September 2014 (right). 82 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Mercury concentrations in preyfish in 2014 were lower than in 2005, which was lower than 1995. Concentrations in some locations were still too high to protect top carnivores such as birds. Concentrations in largemouth bass remain high. 0 50 100 150 200 250 300 Total Mercury in Mosquitofish (ug/kg) Figure 71. Estimates of marsh area for wet season total mercury concentration in mosquitofish (|jg/kg) during Everglades REMAP wet season sampling events: September 1995, November 2005 and September 2014. The estimates indicate improved conditions for mosquitofish in 2005 and 2014 (Wald F test). exceeded the 100 |jg/kg protection level for predators, as compared to 40.1 ±6.7% in 2005 and 59.9 ± 7.3% in 1995. The differences among the curves are statistically significant (Wald F, p<0.05). The changes for the entire Everglades also apply to each subarea (ENP and the three WCAs), as the CDFs for each subarea in 2014 are all different than in 2005 and 1995 (Wald F, p<0.04) (Kalla and Scheidt 2017). Analysis of variance indicated that the lower concentrations observed in 2014 compared to 2005 cannot be explained by fish length or weight. WY15 mercury concentrations reported by FDEP and SFWMD in mosquitofish from 13 locations in the Everglades marsh had an overall median concentration of 44 |jg/kg (Julian et al. 2016b), which is consistent with the median of 33.5 |jg/kg for 2014 REMAP data. In spite of the drop in mosquitofish mercury in REMAP data for the three years sampled, overall annual concentrations of mercury in largemouth bass in the Everglades were unchanged from WY99 to WY14 (Julian et al. 2016b), indicating that biomagnification in the food chain remained high. Bass at 6 of 9 locations sampled in 2014 exceeded the 0.3 mg/kg fish tissue criterion to protect human health, with the highest median concentration of 1.66 mg/kg in bass from Shark River Slough in the Park. From 2000 to 2017, 79% of the largemouth bass across the WCAs and nearly 99% of the bass in the Park in Shark River Slough exceeded the 0.3 mg/kg human health criterion (Lange 2019). 83 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 72. Total mercury (top) arid rnethylmercury (bottom) iri surface water (ng/L) during Everglades REMAP wet season sampling events: November 2005 (left) and September 2014 (right). Total mercury in surface water was also lower during 2014 (Figures 72-73). As compared to 1995, there was a slight increase in 2005 and a larger decrease in 2014, with both differences significant (Wald F, p<0.05). The drop over the whole study area comparing 2014 to 2005 also applies to all four subareas (Wald F, p<0.01) 84 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Methylmercury in Surface Water (ng/L) Total Mercury in Surface Water (ng/L) Figure 73. Estimates of marsh area for wet season methylmercury (top) and total mercury (bottom) concentrations in surface water (ng/L) during Everglades REMAP wet season sampling events: September 1995, November 2005 and September 2014. The estimates indicate improved conditions in 2014 (Wald F test). (Kalla and Scheidt 2017). Methylated mercury is present in concentrations that are an order of magnitude less than total mercury (2014 median of 0.1 ng/L vs. 1.6 ng/L). There was also less methylmercury in surface water over the course of the 1995, 2005 and 2014 REMAP surveys (Figures 72-73), and the differences among the CDF curves are statistically significant (Wald F, p<0.05). The changes for the whole study area also apply to the three WCAs for 2014 compared to 2005 (Wald F, p<0.02), and for 2014 compared to 1995 (Wald F, p<0.01) (Kalla and Scheidt 2017). USGS data for the Park from 2008 to 2013 show interannual variation with no apparent trend in surface water total mercury or methylmercury (Maglio et al. 2015). 85 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Figure 74. Total mercury in soil for REMAP data from all sampling events, 1995 to 2014, expressed as micrograms of mercury per kilogram of soil (left), and micrograms of mercury per cubic centimeter of soil (right). Lower mercury in surface water during 2014 cannot be easily attributed to less atmospheric deposition. Most total mercury in surface water consists of inorganic mercury deposited from the atmosphere (reviewed in Liu et al. 2008a). Atmospheric deposition of mercury is influenced by precipitation, local sources, global sources and air circulation patterns. Although a decline in global atmospheric mercury emissions has been reported in recent years (Zhang et al. 2016), annual wet deposition of mercury in the Everglades remained relatively constant from WY94 to WY14 (Julian et al. 2016b). Likewise, data from the Mercury Deposition Network of the National Atmospheric Deposition Program monitoring station within Everglades National Park indicate no difference in mercury loading for 2005 versus 2014 (Kalla and Scheidt 2017). Mass budget estimates indicate that by farthe soil contains the largest pool of mercury in the Everglades, as compared to water, periphyton, floe, or macrophytes (Stoberet al. 2001b: Liu et al. 2008a, 2015). Total 86 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Š Median ~ 25%-75% I Non-Outlier Range o Outliers * Extremes I Š I V LNWR WCA 3AN WCA3B NESS EMM WCA2 WCA3AS SRS OM TS Š Median ~ 25%-75% I Non-Outlier Range o Outliers * Extremes JL 4 0 o y 61 LNWR WCA 3AN WCA3B NESS EMM WCA 2 WCA3AS SRS OM TS Figure 75. Box and whisker plots of soil total mercury (|jg/kg, left) and methylmercury (pg/kg, right) by subarea, all sampling events from 1995 to 2014. mercury in soil (average of three cores, 0-10 centimeter soil depth) can be expressed on a mass basis (micrograms of mercury per kilogram of soil, median = 130 |jg/kg), or on a volume basis (micrograms of mercury per cubic centimeter of soil, median = 17.9 jjg/cc) (Figures 74-75). When expressed by mass, the higher mercury concentrations occur in the peat soils with higher organic matter content (Figures 56 and 74). These results have a similar landscape pattern as data obtained during 2003-04 by others (Cohen et al. 2009b). The mass of total mercury in the Everglades 0-10 cm soil profile has been reported at 13,482 kg using REMAP data (Liu et al. 2015), and 13,000 kg by others (Cohen et al. 2009b). REMAP data continue to indicate high mercury biomagnification. September 2014 median concentrations fortotal mercury were: surface water, 1.6 parts per trillion; water column periphyton, 19 ppb; benthicperiphyton, 24 ppb; mosquitofish, 33.5 ppb; floe, 120 ppb; and soil 150 ppb (Appendix II). Median MeHg concentrations were: surface water 0.1 parts per trillion, water column periphyton 1.8 ppb; benthic periphyton 0.51 ppb; floe 2.6 ppb; and soil 0.77 ppb. The bioconcentration factor expresses the degree to which mercury accumulates in fish compared to its concentration in surface water. Bioconcentration factors were calculated as the concentration of mercury in mosquitofish divided by the concentration of MeHg in surface water (BCFm). In 2014 the median BCFm was 340,000 as compared to 410,000 in 2005. Nonparametric Spearman rank correlations showed that in 2014 no single REMAP variable had a statistically robust association with mercury in mosquitofish (rho coefficient > 0.7, p < .001) Kalla and Scheidt 2017), although several associations were found in 2005 (Scheidt and Kalla 2007). This result may be due to mosquitofish being exposed to mercury by somewhat different pathways in 2005 versus 2014. This result also could simply reflect the decrease and tighter range in mosquitofish mercury concentrations in 2014. The potential influences of water quality parameters, such as DOC and sulfur, on mercury methylation and mercury bioconcentration to fish and wildlife are discussed in the carbon and sulfur sections of this report. 87 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) LITERATURE CITED Abbey-Lee, R. N., E. N. Gaiserand J. C. Trexler. 2013. Relative roles of dispersal dynamics and competition in determining the isotopic niche breadth of a wetland fish. Freshwater Biology 58:780-792 doi:10.1111/fwb. 12084. Abtew, W., R. S. Huebner, C. Pathak and V. Ciuca. 2007. Appendix 2-8: Water Year 2006 Monthly Inflows and Outflows. In 2016 South Florida Environmental Report. South Florida Water Management District and Florida Department of Environmental Protection. Abtew, W. and V. Ciuca. 2016. Appendix 2-5: water Year 2015 Monthly Inflows and Outflows. In 2016 South Florida Environmental Report. South Florida Water Management District and Florida Department of Environmental Protection. Aich, S., C. W. McVoy, T. W. D. and F. Santamaria. 2013. Estimating soil subsidence and carbon loss in the Everglades Agricultural Area, Florida using geospatial techniques. Agriculture, Ecosystems & Environment 171:124- 133. Aiken, G., D. P. Krabbenhoft, W. H. Orem and C. C. Gilmour. 2006. Dissolved organic matter and mercury in the Everglades: implications for ecosystem restoration, p. 5 in "2006 Greater Everglades Ecosystem Restoration Conference." http://sofia.usas.aov/aeer/2006/index.html Aiken, G. R., C. C. Gilmour, D. P. Krabbenhoft and W. Orem. 2011. Dissolved Organic Matter in the Florida Everglades: Implications for Ecosystem Restoration. Reviews in Environmental Science and Technology 41(S1):217- 248. Anderson, D. L. 1990. A review: soils, nutrition and fertility practices of the Florida sugarcane industry. Soil and Crop Science Society of Florida 49:78-87. Axelrad, D. M., T. D. Atkeson, T. Lange, C. D. Pollman, C. C. Gilmour, W. H. Orem, I. A. Mendelssohn, P. C. Frederick, D. P. Krabbenhoft, G. R.Aiken, D. G. Rumbold, D. J. Scheidt, and P. I. Kalla. 2007. Chapter 3B: Mercury Monitoring, Research and Environmental Assessment in South Florida. In 2007 South Florida Environmental Report. South Florida Water Management District and Florida Department of Environmental Protection. Axelrad, D.M., T. Lange, M. Gabriel, T.D. Atkeson, C.D. Pollman, W.H. Orem, D.J. Scheidt, P.I. Kalla, PC. Frederick and C.C. Gilmour. 2008. Chapter 3B: Mercury and Sulfur Monitoring, Research and Environmental Assessment in South Florida. In: 2008 South Florida Environmental Report. South Florida Water Management District, West Palm Beach, FL. Bancroft, G. T., W. Hoffman, R. Sawicki and J. C. Ogden. 1992. The Importance of the Water Conservation Areas in the Everglades to the Endangered Wood Stork (Mycteria americana). Conservation Biology 6(3):392-398. Bansal, S., S. Lishawa, S. Newman, B. Tangen, D. Wilcox, D. Albert, M. Anteau, M. Chimney, R. Cressey, E. DeKeyser, K. Elgersma, S. Finkelstein, J. Freeland, R. Grosshans, P. Klug, D. Larkin, B. Lawrence, G. Linz, J. Marburger, G. Noe, C. Otto, N. Reo, J. Richards, C. Richardson, L. Rodgers, A. Schrank, D. Svedarsky, S. Travis, N. Tuchman, L. Windham-Myers. 2019. Typha (Cattail) Invasion in North American Wetlands: Biology, Regional Problems, Impacts, Ecosystem Services, and Management. Wetlands (2019) 39:645-684. Barron, M., S. Duvall and K. Barron. 2004. Retrospective and current risks of mercury to panthers in the Florida Everglades. Ecotoxicology 13(3):223-229. Bates, A. L, E. C. Spikerand C. W. Holmes. 1998. Speciation and isotopic composition of sedimentary sulfur in the Everglades, Florida, USA. Chemical Geology 146:155-170. Bates, A. L, W. H. Orem, J. W. Harvey and E. C. Spiker. 2002. Tracing sources of sulfur in the Florida Everglades. Journal of Environmental Quality 31:287-299. 88 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Belanger, T. V., D. J. Scheidt, and J.R. Platko II. 1989. Effects of Nutrient Enrichment on the Florida Everglades. Lake and Reservoir Management 5:101-111. Beltman, B., T. G. Rouwenhorst, M. B. Van Kerkhoven, T. Van Der Krift and J. T. A. Verhoeven. 2000. Internal eutrophication in peat soils through competition between chloride and sulphate with phosphate for binding sites. Biogeochemistry 50(2): 183-194. Benoit, J., R. Mason and C. Gilmour. 1999. Estimation of mercury-sulfide speciation in sediment pore waters using the octanol-water partitioning and implications for availability to methylating bacteria. Environmental Toxicology and Chemistry 18(10):2138-2141. Bhadha, J.H., T. A. Lang and S. M. Daroub. 2014. Seasonal delivery of organic matter and metals to farm canals: effect on sediment phosphorus storage capacity. Journal of Soils and Sediments 14:991-1003. Breidt, F.J., and W.A. Fuller. 1999. Design of supplemented panel surveys with application to the national resources inventory. Journal of Agricultural, Biological, and Environmental Statistics 4(4):391-403. Brezonik, P. L. and A. Federico. 1975. Effects of backpumping from agricultural drainage canals on water quality in Lake Okeechobee. Florida Department of Environmental Regulation Technical Series Volume 1, Number 1. Tallahassee, Florida. 64 pages. Browder, J. A., P. J. Gleason and D. R. Swift. 1993. Periphyton in the Everglades: spatial variation, environmental correlates and ecological implications. Pages 379-418 in Everglades: The Ecosystem and Its Restoration. Davis, Steven M. and John C. Ogden (editors). St. Lucie Press. Delray Beach, Florida. 826 pages. Bruland, G. L, T. Z. Osborne, K. R. Reddy, S. Grunwald, S. Newman and W. F. DeBusk. 2007. Recent changes in soil total phosphorus in the Everglades: Water Conservation Area 3. Environmental Monitoring and Assessment 129:379- 395. Burns and McDonnell. 2003. Everglades Protection Area Tributary Basins Long-Term Plan for Achieving Water Quality Goals. Submitted to the SFWMD, West Palm Beach, FL. October 2003. 402 pages. CH2MHILL. 1978. Water quality studies in the Everglades Agricultural Area. Submitted to the Florida Sugarcane League. Gainesville, Florida. 136 pages. Chen, M., S. Daroub, T. Lang and O. Diaz. 2006. Specific conductance and ionic characteristics of farm canals in the Everglades Agricultural Area. Journal of Environmental Quality 35:141-150. Chen, C., E. Meselhe and M. Waldon. 2012. Assessment of mineral concentration impacts from pumped stormwater on an Everglades Wetland, Florida, USA - Using a spatially-explicit model. Journal of Hydrology 452-453:25-39. Chen, H., Ivanoff, D., Pietro, K., 2015. Long-term phosphorus removal in the Everglades stormwater treatment areas of South Florida in the United States. Ecological Engineering 79:158-168. Childers, D. L., R. F. Doren, R. Jones, G. B. Noe, M. Rugge and L. J. Scinto. 2003. Decadal change in vegetation and soil phosphorus pattern across the Everglades landscape. Journal of Environmental Quality 32:344-362. Chimney, M. J. and G. Goforth. 2006. History and description of the Everglades Nutrient Removal Project, a subtropical constructed wetland in south Florida (USA). Ecological Engineering 27(4):268-278. Coale, F. J. 1994. Sugarcane production in the EAA. Chapter 9 (pages 224-237) In Everglades Agricultural Area (EAA): Water, Soil, Crop and Environmental Management. A. B. Bottcher and F. T. Izuno, editors. University Press of Florida, Gainesville. 322 pages. 89 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Cohen, M. J., T. Z. Osborne, S. Lamsal and M. W. Clark. 2009a. Regional Distribution of Soil Nutrients - Hierarchical Soil Nutrient Mapping for Improved Ecosystem Change Detection. Final Report. Central Everglades Restoration Project Monitoring and Assessment Plan Greater Everglades (GE) 3.1.3.2. University of Florida. Gainesville, FL. Cohen M.J., S. Lamsal, T. Z. Osborne, J. C. Bonzongo, K. R. Reddy and S. Newman. 2009b. Soil total mercury concentrations across the greater Everglades. Journal of the Soil Science Society of America 73:675-685. Congressional Research Service. 2017. Everglades Restoration: Federal Funding and Implementation Progress. 16 pages. Cornwell, G.W., R.L. Downing, A.R. Marshall, J.N. Layne, and C.M. Loveless. 1970. Everglades water and its ecological implications. Report of the Special Study Team of the Florida Everglades. 42 pages. Corrales, J., G.M. Naja, C. Dziuba, R.G. Rivero and W. Orem. 2011. Sulfate threshold target to control methylmercury levels in wetland ecosystems. Science of The Total Environment 409(11):2156-2162. Corstanje, R., S. Grunwald, K. R. Reddy, T. Z. Osborne and S. Newman. 2006. Assessment of the spatial distribution of soil properties in a northern Everglades marsh. Journal of Environmental Quality 35:938-949. Cox, S., L. Douglas, S. McCollum, M. Bledsoe, and R. Marrotte. 1988. Subsidence Study of the Everglades Agricultural Area. USDASoil Conservation Service Greenacres Field Office. Greenacres, Florida. 27 pages. Craft, C. B. and C. J. Richardson. 1993. Peat accretion and phosphorus accumulation along a eutrophication gradient in the northern Everglades. Biogeochemistry 22:133-156. Daroub, S. H., T. A. Lang, O. A. Diaz and S. Grunwald. 2009. Long-term Water Quality Trends after Implementing Best Management Practices in South Florida. Journal of Environmental Quality 38:1683-1693. Daroub, S. H., S. Van Horn, T. A. Lang and O. A. Diaz, 2011. Best Management Practices and Long-Term Water Quality Trends in the Everglades Agricultural Area. Critical Reviews in Environmental Science and Technology, 41 :(S1) 608-632, DPI: 10.1080/10643389.2010.530905. Davis, J. H., Jr. 1946. The Peat Deposits of Florida: Their Occurrence, Development and Uses. Geological Bulletin No. 30. Florida Geological Survey. Tallahassee, Florida. 247 pages. Davis, S. M. 1994. Phosphorus inputs and Vegetation Sensitivity in the Everglades. Pages 357-378 in Everglades: The Ecosystem and its Restoration. S. M. Davis and J. C. Ogden, editors. St. Lucie Press, Delray Beach, Florida. Davis, S. M. and J. C. Ogden. 1994. Everglades: The Ecosystem and Its Restoration. St. Lucie Press. Delray Beach, Florida. 826 pages. de Camargo Santos, A., J. M. McCray, S. H. Daroub, D. L. Rowland, S. Ji and H. Sandhu. 2020. Nitrogen Assessment of Shallow Florida Histosols, Communications in Soil Science and Plant Analysis, 51:14, 1916-1929, DPI: 10.1080/00103624.2020.1798990. DeBusk, T. A., F. E. Dierberg, W. F. DeBusk, S. D. Jackson, J. A. Potts, S. C. Galloway, D. Sierer-Finn, B. Gu. 2015. Sulfide concentration effects on Typha domingensis Pers. (cattail) and Cladium jamaicense Crantz (sawgrass) growth in Everglades marshes. Aquatic Botany 124:78-84. Diaz-Ramos, S., D.L. Stevens, Jr., and A.R. Plsen. 1996. EMAP statistical methods manual. USEPA, Corvallis, PR. EPA/620/R-96/002. Dierberg, F. E., T. A. DeBusk, M. Jerauld and B. Gu. 2014. Appendix 3B-1: Evaluation of Factors Influencing Methylmercury Accumulation in South Florida Marshes. In: 2014 South Florida Environmental Report. South Florida Water Management District, West Palm Beach, FL. 90 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Dineen, J. W. 1973. Memorandum from Walter J. Dineen, Environmental Sciences Division, to Bill Storch. North Springs Improvement District. April 10, 1973. South Florida Flood Control District. Dineen, J. W. 1986. Memorandum from J. W. Dineen, Director, Environmental Sciences Division, to Pat Bidol, Executive Program Director, Executive Office. Water Conservation Areas. August 28, 1986. South Florida Flood Control Management District. Ding, Y., K. M. Cawley, C. Nunes da Cunha and R. Jaffe. 2014. Environmental dynamics of dissolved black carbon in wetlands. Bioaeochemistrv DPI 10.1007/s10533-014-9964-3. Doren, R. F., T. V. Armentano, L. D. Whiteaker, and R. D. Jones. 1996. Marsh Vegetation Patterns and Soil Phosphorus Gradients in the Everglades Ecosystem. Aquatic Botany 56:145-163. Dreschel, T., S. Hohner, S. Aich and C. McVoy. 2018. Peat Soils of the Everglades of Florida, USA. In Peat, B. Topcuoglu and M. Turan, editors. DPI: 10.5772/intechopen.72925 Duvall, S. and M. Barron. 2000. A screening level probabilistic risk assessment of mercury in Florida Everglades food webs. Ecotoxicology and Environmental Safety 47:298-305. Eisler, R. 1987. Mercury hazards to fish, wildlife, and invertebrates: A synoptic review. U.S. Fish and Wildlife Service Biological Report 85 (1.10). 90 pages. Florida Administrative Code Chapter 62-302. Florida Surface Water Quality Standards. April 30, 2018 version. Fan, X., B. Gu, E. A. Hanlon, Y. Li, K. Migliaccio and T. W. Dreschel. 2010. Investigation of long-term trends in selected physical and chemical parameters of inflows to Everglades National Park, 1977-2005. Environmental Monitoring and Assessment DPI 10.1007/s10661-010-1710-2. Feng, S., Z. Ai, S. Zheng, B. Gu and Y. Li. 2014. Effects of Dryout and Inflow Water Quality on Mercury Methylation in a Constructed Wetland. Water Air & Soil Pollution 225:1929. DPI 10.1007/s11270-014-1929-6. Fink, L. and P. Rawlik. 2000. The Everglades Mercury Problem. Chapter 7 in 2000 Everglades Consolidated Report. SFWMD. West Palm Beach, FL. Flora, M. D. and P. C. Rosendahl. 1982a. The Response of Specific Conductance to Environmental Conditions in the Everglades National Park, Florida. Water, Air and Soil Pollution 17:51-59. Flora, M. D. and P. C. Rosendahl. 1982b. Historical changes in the conductivity and ionic characteristics of the source water for the Shark River Slough, Everglades National Park, Florida, U. S. A. Hydrobiologia 97:249-254. FDEP 1997. The South Florida Mercury Science Program. Tallahassee, FL. 4 pages. FDEP 2013. Mercury Total Maximum Daily Load (TMDL) for the State of Florida. Pctober 24, 2013. Florida Department of Environmental Protection, Tallahassee, FL. 120 pages. Florida Department of Health and Rehabilitative Services. 1989. Health advisory for Florida fisherman. FDPH. 2020. Your guide to eating fish caught in Florida, May 2020 v2. 61 pages. http://www.floridahealth.aov/proarams-and-services/prevention/healthv-weiaht/nutrition/seafood-consumption/ documents/fish-advisorv-bia-book2020.pdf Flower H., M. Rains, H. C. Fitz, W. Prem, S. Newman, T. Z. Psborne, K. R. Reddy, J. Pbeysekera and J. Barnes. 2019. Shifting ground: landscape-scale modeling of soil biogeochemistry under climate change in the Florida Everglades. Environmental Management 64:416-435. 91 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Frydenborg, R. and L. Ross. 1987. Memorandum through Jerry Brooks to Pam McVety. Florida Department of Environmental Regulation. Projects- general; Everglades Water Conservation Areas. 5 pages. Gabriel, M.C., N. Howard and T.Z. Osborne. 2014a. Fish mercury and surface water sulfate relationships in the Everglades Protection Area. Environmental Management 53(3):583-593. Gabriel, M. C., D. Axelrad, W. Orem and T. Z. Osborne. 2014b. Response to Julian et al. (2014) "Comment on and Reinterpretation of Gabriel et al. (2014) 'Fish Mercury and Surface Water Sulfate Relationships in the Everglades Protection Area'. Environmental Management 55(6): 1227-1231. Gaiser, E. F., L. J. Scinto, J. H. Richards, K. Jayachandran, D. L. Childers, J. C. Trexlerand R. D. Jones. 2004. Phosphorus in periphyton mats provides the best metric for detecting low-level P enrichment in an oligotrophic wetland. Water Research 38:507-516. Gaiser, E. E., J. C. Trexler, J. H. Richards, D. L. Childers, D. Lee, A. L. Edwards, L. J. Scinto, K. Jayachandran, G. B. Noe and R. D. Jones. 2005. Cascading ecological effects of low-level phosphorus enrichment in the Florida Everglades. Journal of Environmental Quality 34:1-7. Gaiser, E.E. 2009. Periphyton as an indicator of restoration in the Florida Everglades. Ecological Indicators 9S:S37-S45. Gaiser, E. E. P. V. McCormick, S. E. Hagerthey and A. D. Gottlieb. 2011. Landscape Patterns of Periphyton in the Florida Everglades. Reviews in Environmental Science and Technology 41 (S1):92-120. Gaiser, E., A. D. Gottlieb, S. S. Lee and J. C. Trexler. 2015. The Importance of Species-Based Microbial Assessment of Water Quality in Freshwater Everglades Wetlands. Chapter 6 in Microbiology of the Everglades Ecosystem. James A. Entry, Andrew D. Gottlieb, Krish Jayachandran, Andrew Ogram, editors. CRC Press. Gilbert, R. A. and R. W. Rice. 2006. Nutrient requirements for sugarcane production on Florida muck soils. University of Florida Institute of Food and Agricultural Sciences Extension publication SS-AGR-226. 4 pages. Gilmour, C. C., W. Orem, D. Krabbenhoft, and I. Mendelssohn. 2007. Appendix 3B-3: preliminary assessment of sulfur sources, trends and effects in the Everglades. In 2007 South Florida Environment Report. South Florida Water Management District and Florida Department of Environmental Protection. West Palm Beach, Florida. 83 pages. Gleason, P. J. 1974. Chemical Quality of Water in Conservation Area 2A and Associated Canals. South Florida Water Management District Technical Publication 74-1. West Palm Beach, Florida. 72 pages plus appendices. Gleason, P. J. and W. Spackman, Jr. 1974. Calcareous periphyton and water chemistry in the Everglades. Pages 146-181 in Environments of south Florida: present and past edited by Patrick J. Gleason. Miami Geological Society. Miami, Florida. Gleason, P., P. Stone, D. Hallett and M. Rosen. 1975a. Preliminary Report on the Effect of Agricultural Runoff on the Periphytic Algae of Conservation Area 1. SFWMD. West Palm Beach, Florida. 69 pages. Gleason, P., P. A. Stone, P., S. M. Davis, M. Zaffke and L. Harris. 1975b. The impact of agricultural runoff on the Everglades marsh located in the Conservation Areas of the Central and Southern Florida Flood Control District. CSFFCD Resource Planning Department, West Palm Beach, Florida. 119 pages plus appendices. Gleason, P. J. and P. Stone. 1993. Age, Origin and Landscape Evolution of Everglades Peatland. Pages 149-197 in Everglades: The Ecosystem and Its Restoration. Davis, Steven M. and John C. Ogden (editors). St. Lucie Press. Delray Beach, Florida. 92 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Godfrey. M. C. and T. Catton. 2006. River of Interests: Water Management in South Florida and the Everglades, 1948-2000. Historical Research Associates, Inc. 360 pages. Gordon, A. S., W. J. Cooper and D. J. Scheidt. 1986. Denitrification in marl and peat sediments in the Florida Everglades. Applied and Environmental Microbiology 52(5):987-991. Gottlieb, A., E. Gaiser and S. Lee. 2015. Changes in Hydrology, Nutrient Loading and Conductivity in the Florida Everglades, and Concurrent Effects on Periphyton Community Structure and Function. Chapter 7 in Microbiology of the Everglades Ecosystem. J. A. Entry, A. D. Gottlieb, K. Jayachandran, A. Ogram (editors). CRC Press. Boca Raton, Florida. Graham, A. M., G. R. Aiken and C. C. Gilmour. 2013. Effect of Dissolved Organic Matter Source and Character on Microbial Hg Methylation in Hg-S-DOM Solutions. Environmental Science and Technology 47:5746-5754. doi. ora/10.1021 /es400414a. Graham, A. M. 2019. Dissolved Organic Matter Interactions with Mercury in the Florida Everglades. Chapter 4 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration. https://doi. ora/10.1007/978-3-030-32057-7 4. Grunwald, S., K. R. Reddy, S. Newman and W. F. DeBusk. 2004. Spatial variability, distribution and uncertainty assessment of soil phosphorus in a south Florida wetland. Environmetrics 15:811-825. Gu, B., M. J. Chimney, J. Newman and M. K. Nungesser. 2006. Limnological characteristics of a subtropical constructed wetland in south Florida (USA). Ecological Engineering 27:345-360. Gunderson, L. H., and W. F. Loftus. 1994. The Everglades. Pages 199-255 In Biodiversity of the southeastern United States. W. H. Martin, S. C. Boyce, and A. C. Echternacht (editors). John Wiley, New York, New York, USA. Hagerthey, S. E., M. I. Cook, R. M. Kobza, S. Newman and B. J. Bellinger. 2014. Aquatic faunal responses to an induced regime shift in the phosphorusDimpacted Everglades. Freshwater Biology 59(7):1389-1405. Haider, S., Swain, E., Beerens, J., Petkewich, M., McCloskey, B., and Henkel, H. 2020. The Everglades Depth Estimation Network (EDEN) surface-water interpolation model, version 3: USGS Scientific Investigations Report 2020-5083. 31 pages. https://doi.ora/10.3133/sir20205083. Harwell, M. C., D. D. Surratt, D. M. Barone and N. G. Aumen. 2008. Conductivity as a tracer of agricultural and urban runoff to delineate water quality impacts in the northern Everglades. Environmental Monitoring and Assessment 147: 445-462. Hem, J.D., 1985, Study and interpretation of the chemical characteristics of natural water (3d ed.): USGS Water- Supply Paper 2254, 264 pages, https://pubs.usas.aov/wsp/wsp2254/. Hohner, S. M., and T. W. Dreschel. 2015. Everglades peats: using historical and recent data to estimate pre-drainage and current volumes, masses, carbon contents and water retention capacities. Mires and Peat 16:1-15. Ingebritsen, S. E., C. McVoy, B. Glaz, and W. Park. 1999. Florida Everglades: Subsidence threatens agriculture and complicates ecosystem restoration. Pages 95-106 in "Land Subsidence in the United States". USGS Circular 1182. Reston, Virginia. Inglett, P. W., V. H. Rivera-Monroy and J. R. Wozniak. 2011. Biogeochemistry of Nitrogen Across the Everglades Landscape. Critical Reviews in Environmental Science and Technology 41(S1):187-216. ludicello, J., T. Davison and P. Wade. 2016. Chapter 4: Nutrient Source Control Programs. In 2016 South Florida Environmental Report. SFWMD and FDEP 93 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Izuno, F.T., Rice, R.W. and L. T. Capone. 1999. Best Management Practices to Enable the Coexistence of Agriculture and the Everglades Environment. HortScience 34:27-33. James, R. T. and R V. McCormick. 2012. The sulfate budget of a shallow subtropical lake. Fundamentals of Applied Limnology 181(4):253-269. Jeremiason, J., D. Engstrom, E. Swain, E. Nater, B. Johnson, J. Almendinger, B. Monson and R. Kolka. 2006. Sulfate addition increases methylmercury production in an experimental wetland. Environmental Science and Technology 40:3800-3806. Jiang, R, G. Liu, W. Cui and Y. Cai. 2018. Geochemical modeling of mercury speciation in surface water and implications on mercury cycling in the Everglades wetland. Science of the Total Environment 640-641:454-465. Julian, P. 2013a. Mercury bio-concentration factor in mosquito fish (Gambusia spp.) in the Florida Everglades. Bulletin of Environmental Contamination and Toxicology 90:329-332. Julian, P. 2013b. Mercury hotspot identification in Water Conservation Area 3, Florida USA. Annals of GIS 19:79-88. Julian, P. and B. Gu. 2014. Mercury accumulation in largemouth bass (Micropterus salmoides Lacepede) within marsh ecosystems of the Florida Everglades, USA. Ecotoxicoloav DPI 10.1007/s10646-014-1373-9. Julian P., B. Gu and G. Redfield. 2015. Comment on and reinterpretation of Gabriel et al. (2014) fish mercury and surface water sulfate relationships in the Everglades Protection Area. Environmental Management 55:1-5. Julian, P., G. G. Payne and S. K. Liu. 2016a. Chapter 3A. Water quality in the Everglades Protection Area. In 2016 South Florida Environmental Report. South Florida Water Management District and Florida Department of Environmental Protection. 50 pages. Julian, P., B. Gu, G. Redfield and K. Weaver, editors. 2016b. Chapter 3B: Mercury and sulfur environmental assessment for the Everglades. In 2016 South Florida environmental report. South Florida Water Management District and Florida Department of Environmental Protection. 47 pages. Julian, P. 2016c. Water Year 2015 Attainment of the Everglades Dissolved Oxygen Site-Specific Alternative Criteria at Individual Stations. Appendix 3A-3 In 2016 South Florida Environmental Report. South Florida Water Management District and Florida Department of Environmental Protection. 9 pages. Julian, P. 2016d. Appendix 3A-6: Water Year 2011-2015 Annual Total Phosphorus Criteria Compliance Assessment. In 2016 South Florida Environmental Report. Florida Department of Environmental Protection. Julian, P., A. L. Wright and T. Z. Osborne. 2016e. Iron and Sulfur Porewater and Surface Water Biogeochemical Interactions in Subtropical Peatlands. Soil Science Society of America Journal 80:794-802. doi: 10.2136/ sssai2015.11.0418. Julian, P., B. Gu and K. Weaver, editors. 2017a. Chapter 3B: Mercury and Sulfur Environmental Assessment for the Everglades. In 2017 South Florida Environmental Report. South Florida Water Management District. West Palm Beach, FL. 58 pages. Julian, P., S. Gerber, A. L. Wright, B. Gu and T. Z. Osborne. 2017b. Carbon pool trends and dynamics within a subtropical peatland during long-term restoration. Ecological Processes 6:43. https://doi.ora/10.1186/s13717-017- 0110-8. Juston, J. M. and T. A. DeBusk. 2011. Evidence and implications of the background phosphorus concentration of submerged aquatic vegetation wetlands in Stormwater Treatment Areas for Everglades restoration. Water Resources Research 47(W01511), doi: 10.1029/2010WR009294. 94 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Kalla, P.I., and D.J. Scheidt. 2017. Everglades ecosystem assessment - Phase IV, 2014: Data reduction and initial synthesis. USEPA, Science and Ecosystem Support Division, Athens, GA. SESD 14-0380. 58 pages. Kalla, P., M. Parsons, and J. Ackerman. 2017. Operating procedure for bottom water sampling for sulfide. SESDPROC-515-RO. USEPA, Region 4 Science and Ecosystem Support Division, Athens, GA. 22 pages. Kalla, P., M. Cyterski, J. Minucci and D. Scheidt. 2019. Modelling of mercury bioaccumulation in mosquitofish (Gambusia sp.) for the Everglades Regional Monitoring and Assessment Program. LSASD 20-0085. USEPA, Laboratory Services and Applied Science Division, Athens, GA. 35 pages. Kalla, P., M. Cyterski, D. Scheidt and J. Minucci. 2021. Spatiotemporal Effects of Interacting Water Quality Constituents on Mercury in a Common Prey Fish in a Large, Perturbed, Subtropical Wetland. Science of the Total Environment 792 (2021) 148321. https://doi.Org/10.1016/j.scitotenv.2021.148321. Lamers, L. P. M., H. B. M. Tomassen, and J. G. M. Roelofs. 1998. Sulfate-induced eutrophication and phytotoxicity in freshwater wetlands. Environmental Science and Technology 32:199-205. Lamers, L. P. M., S. Falla, E. M. Samborska, I. A. R. van Dulken, G. van Hengstum and J. G. M. Roelofs. 2002. Factors controlling the extent of eutrophication and toxicity in sulfate-polluted freshwater wetlands. Limnology and Oceanography 47(2):585-593. Landing, W. 2015. Appendix 3B-2: Peer-Review Report on the Everglades Agricultural Area Regional Sulfur Mass Balance: Technical Webinar. In 2015 South Florida Environmental Report. SFWMD and FDEP West Palm Beach, Florida. 45 pages. Lang, T.A., O. Oladeji, M. S. Josan, and S. H. Daroub. 2010. Environmental and management factors that influence drainage water P loads from Everglades Agricultural Area farms of South Florida. Agriculture, Ecosystems and Environment 138:170-180. Lange, T. 2019. Comparison of Everglades Fish Tissue Mercury Concentrations to Those for Other Fresh Waters. Chapter 9 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration, D. G. Rumbold et al. (eds.). https://doi.orq/10.1007/978-3-030-32057-7 9. Li, S., I.A. Mendelssohn, H. Chen and W.H. Orem. 2009. Does sulphate enrichment promote the expansion of Typha domingensis (cattail) in the Florida Everglades. Freshwater Biology doi: 10.1111/i. 1365-2427.2009.02242.x. Li, Y., Z. Duan, G. Liu, P. Kalla, D. Scheidt and Y. Cai. 2015. Evaluation of the possible sources and controlling factors of toxic metals/metalloids in the Florida Everglades and their potential risk of exposure. Environmental Science and Technology DQI:10.1021/acs.est.5b01638. Light, S. S., and J.W. Dineen. 1994. Water control in the Everglades: A historical perspective. Pages 47-84 in Everglades: the ecosystem and its restoration. S. M. Davis and J. C. Ogden, editors. St. Lucie Press, Delray Beach, Florida, USA. Liu, G., Y. Cai, P. Kalla, D. Scheidt, J. Richards, L. J. Scinto, E. Gaiser, and C.Appleby. 2008a. Mercury Mass Budget Estimates and Cycling Seasonality in the Florida Everglades. Environmental Science and Technology 42:1954-1960. Liu, G., Y. Cai, T. Philippi, P. Kalla, D. Scheidt, L. Scinto, J. Richards and J. Trexler. 2008b. Distribution of total and methylmercury in different ecosystem compartments in the Everglades: Implications for mercury bioaccumulation. Environmental Pollution 153:257-265. Liu, G., Y. Cai, Y. Mao, D. Scheidt, P. Kalla, J. Richards, L. Scinto, G. Tachiev, D. Roelant and C. Appleby. 2009. Spatial Variability in Mercury Cycling and Relevant Biogeochemical Controls in the Florida Everglades. Environmental Science and Technology 43(12):4361-4366. DOI: 10.1021/es803665c. 95 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Liu, G., G. M. Naja, P. Kalla, D. Scheidt, E. Gaiser, and Y. Cai. 2011. Legacy and fate of mercury and methylmercury in the Florida Everglades. Environmental Science and Technology, 45(2):496-501. Supporting information S1-S41. Liu, G., Y. Cai, G. Tachiev and L. Lagos. 2015. Mercury Mass Budget Estimates and Cycling in the Florida Everglades. Chapter 4 in Microbiology of the Everglades Ecosystem. J. A. Entry, A. D. Gottlieb, K. Jayachandran, A. Ogram (editors). CRC Press. Loftus, W. F., J. C. Trexler and R. D. Jones. 1998. Mercury transfer through an Everglades aquatic food web. Final Report. Contract SP-329 to the Florida Department of Environmental Protection. LOTAC II. 1988. Lake Okeechobee Technical Advisory Council (LOTAC) II. Volume II. Interim Report to the Governor, State of Florida Legislature. South Florida Water Management District, West Palm Beach, Florida. 74 pages. Maglio, M., D. Krabbenhoft, M. Tate, J. DeWild, J. Ogorek, C. Thompson, G. Aiken, W. Orem, J. Kline, J. Castro and C. Gilmour. 2015. Drivers of Geospatial and Temporal Variability in the Distribution of Mercury and Methylmercury in Everglades National Park. Conference on Greater Everglades Ecosystem Restoration (GEER). Coral Springs, Florida. Marshall, A. R. 1971. Can the Everglades damage be repaired? Florida Naturalist. October 1971 pages 104-106, 118. Mattraw, H. C. Jr., D. J. Scheidt and A. C. Federico. 1987. Analysis of trends in water quality data for Water Conservation Area 3A, the Everglades, Florida. USGS Water Resources Investigation 87-4142. Tallahassee, Florida. 52 pages. McCormick, P.V., M.J. Chimney, and D.R. Swift. 1997. Diel oxygen profiles and water column community metabolism in the Florida Everglades, USA. Archives Hydrobiologie 140:117-129. McCormick, P. C., S. Newman, S. Miao, R. Reddy, D. Gawlik, C. Fitz, T. Fontaine and D. Marley. 1999. Ecological Needs of the Everglades. Chapter 3 in: 1999 Everglades Consolidated Report. South Florida Water Management District. McCormick., P. V., S. Newman, S. Miao, D. E. Gawlik, D. Marley, K. R. Reddy, T. D. Fontaine. 2002. Effects of anthropogenic phosphorus inputs on the Everglades. Pages 83-126 in The Everglades, Florida Bay, and coral reefs of the Florida Keys: An ecosystem sourcebook. J. W. Porter and K. G. Porter (editors). CRC Press LLC, Boca Raton, FL. McCormick, P. V. and J. A. Laing. 2003. Effects on increased phosphorus loading on dissolved oxygen in a subtropical wetland, the Florida Everglades. Wetlands Ecology and Management 11:199-216. McCormick, P. V. and J. W. Harvey. 2007. Influence of changing water sources and mineral chemistry on the Everglades ecosystem. U. S. Geological Survey Administrative Report. Reston, Virginia. 54 pages. McCormick, P. V., S. Newman, and L. Vilchek. 2009. Landscape responses to wetland eutrophication: loss of slough habitat in the Florida Everglades, USA. Hydrobiologia 621:105-114. McCormick, P. V. 2011. Soil and periphyton indicators of anthropogenic water-quality changes in a rainfall-driven wetland. Wetlands Ecology and Management 19:19-34. DPI 10.1007/s11273-010-9196-9. McCormick, P. V., J. W. Harvey and E. S. Crawford. 2011. Influence of Changing Water Sources and Mineral Chemistry on the Everglades Ecosystem. Critical Reviews in Environmental Science and Technology 41(S1):28-63. McPherson, B. F., B. G Waller and H. C. Mattraw. 1976. Nitrogen and phosphorus uptake in the Everglades Conservation Areas, Florida, with special reference to the effects of backpumping runoff. USGS Water Resources Investigation 76-29. Tallahassee, Florida. 120 pages. 96 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Miller, W.L. 1988. Description and Evaluation of the Effects of Urban and Agricultural Development on the Surficial Aquifer System, Palm Beach County, Florida. Report 88-4056. USGS, Tallahassee, FL. 66 pages. Morgan, K.T, J. M. McCray, R. W. Rice, R. A. Gilbert, and L. E. Baucum. 2009. Review of Current Sugarcane Fertilizer Recommendations: A Report from the UF/IFAS Sugarcane Fertilizer Standard Task Force. Soil and Water Department, University of Florida IFAS Extension Publication SL295. 5 pages. Naja, M., D. L. Childers and E. E. Gaiser. 2017. Water quality implications of hydrologic restoration alternatives in the Florida Everglades, United States. Restoration Ecology doi: 10.1111/rec.12513. National Atmospheric Deposition Program (NADP). 2014a. National Atmospheric Deposition Program, 2014 annual data summary, site FL11. http://nadp.sws.uiuc.edu/default.html Accessed March 30, 2017. National Atmospheric Deposition Program (NADP), 2014b. National Atmospheric Deposition Program, 2014 monthly data, site FL11. http://nadp.sws.uiuc.edu/default.html Accessed March 30, 2017. National Research Council. 2008. Progress Toward Restoring the Everglades. The second biennial review-2008. National Academies Press. Washington, D. C. 340 pages. National Research Council. 2010. Progress Toward Restoring the Everglades. The third biennial review-2010. National Academies Press. Washington, D. C. 277 pages. National Research Council. 2012. Progress Toward Restoring the Everglades. The fourth biennial review-2012. National Academies Press. Washington, D. C. 227 pages. National Research Council. 2014. Progress Toward Restoring the Everglades. The fifth biennial review-2014. National Academies Press. Washington, D. C. 254 pages. Neto, R. R., R. N. Mead, J. W. Louda and R. Jaffe. 2006. Organic biogeochemistry of detrital flocculent material (floe) in a subtropical, coastal wetland. Biogeochemistry 77:283-304. Newman S. and Hagerthey S.E. 2011. Water Conservation Area 1- a case study of hydrology, nutrient, and mineral influences on biogeochemical processes. Critical Reviews in Environmental Science and Technology 41:702-722. Noe, G. B., D. L. Childers, and R. D. Jones. 2001. Phosphorus biogeochemistry and the impact of phosphorus enrichment: Why is the Everglades so unique? Ecosystems 4:603-624. Nungesser, M. K., C. Saunders, C. A. Coronado-Molina, J. T. B. Obeysekera, J. R. Johnson, C. W. McVoy, B. Benscoter. 2014. Potential effects of climate change on Florida's Everglades. Environmental Management DOI: 10.1007/S00267-014-0417-5. Ogden, J. C., W. B. Robertson, G. E. Davis, and T. Schmidt. 1974. Pesticides, Polychlorinated Biphenols and Heavy Metals in Upper Food Chain Levels, Everglades National Park and Vicinity. Prepared for the Department of Interior, South Florida Environmental Project: Ecological Report No. DI-SFEP-74-16. Distributed by National Technical Information Service, U.S. Department of Commerce, Springfield, Virginia. 27 pages. Ogden, J. C. 1994. AComparison of Wading Bird Nesting Colony Dynamics (1931-1946 and 1971-1989) as an Indication of Ecosystem Conditions in the Southern Everglades. Pages 533-570 in Everglades: The Ecosystem and Its Restoration. Davis, Steven M. and John C. Ogden (editors). St. Lucie Press. Delray Beach, Florida. Olsen, A. R., Sedransk, J., Edwards, D., Gotway, C. A., Liggett, W. 1999. Statistical Issues for Monitoring Ecological and Natural Resources in the United States. Environmental Monitoring and Assessment 54:1-45. Orem, W., C. Gilmour, D. Axelrad, D. Krabbenhoft, D. Scheidt, P. Kalla, P. McCormick, M. Gabriel and G. Aiken. 2011. Sulfur in the South Florida ecosystem: Distribution, sources, biogeochemistry, impacts, and management for restoration. Critical Reviews in Environmental Science and Technology 41 (Supplement 1):249-288. 97 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) OremW., H. C. Fitz, D. Krabbenhoft, M. Tate, C. Gilmourand M. Shafer. 2014a. Modeling sulfate transport and distribution and methylmercury production associated with Aquifer Storage and Recovery implementation in the Everglades Protection Area. Sustainability of Water Quality and Ecology 3-4:33-46. Orem, W., S. Newman, T. Z. Osborne and K. R. Reddy. 2014b. Projecting Changes in Everglades Soil Biogeochemistry for Carbon and Other Key Elements, to Possible 2060 Climate and Hydrologic Scenarios. Environmental Management 55: 776-798. Orem, W. H., D. P. Krabbenhoft, B. A. Poulin, and G. R. Aiken. 2019. Sulfur Contamination in the Everglades, a Major Control on Mercury Methylation. Chapter 2 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration. D. G. Rumbold et al. (eds.). https://doi.ora/10.1007/978-3-030-32057-7 2. Orem, W. H., C. Fitz, D. P. Krabbenhoft, B. A. Poulin, M. S. Varonka and G. R. Aiken. 2020. Ecosystem-Scale Modeling and Field Observations of Sulfate and Methylmercury Distributions in the Florida Everglades: Responses to Reductions in Sulfate Loading. Aquatic Geochemistry https://doi.ora/10.1007/s10498-020-09368-w. Osborne, T. Z., S. Newman, D. J. Scheidt, P. I. Kalla, G. L. Bruland, M. J. Cohen, L. J. Scinto, and L. R. Ellis. 2011a. Landscape patterns of significant soil nutrients and contaminants in the greater Everglades ecosystem: past, present, and future. Critical Reviews in Environmental Science and Technology 41:121-148. Osborne T. Z., G. L. Bruland, S. Newman, K. R. Reddy and S. Grunwald. 2011b. Spatial distributions and eco- partitioning of soil biogeochemical properties in the Everglades National Park. Environmental Monitoring and Assessment 183(1-4): 395-408. Osborne, T. Z, K. R. Reddy, L. R. Ellis, N. G. Aumen, D. D. Surratt. 2013. Evidence of Recent Phosphorus Enrichment in Surface Soils of Taylor Slough and Northeast Everglades National Park. Wetlands DPI 10.1007/s13157-013-0381-5. Osborne, T., S. Newman, K. R. Reddy, L. R. Ellis and M. Ross. 2015. Spatial Distribution of Soil Nutrients in the Everglades Protection Area. Chaper 3 in Microbiology of the Everglades Ecosystem. James A. Entry, Andrew D. Gottlieb, Krish Jayachandran, Andrew Ogram (editors). CRC Press. Boca Raton, Florida. 67 pages. Osborne, T., H. C. Fitz and S. Davis. 2017. Restoring the foundation of the Everglades ecosystem: assessment of edaphic responses to hydrologic restoration scenarios. Restoration Ecology 25(S1): S59-S70. Ouellette, K, K. Alsharif, J. Capece and H. Torres. 2018. The concept of water storage on agriculture lands: Exploring the notion in South Florida. Water Science 32:138-150. Payne, G., T. Bennett and K. Weaver. 2001. Chapter 3: Ecological effects of phosphorus enrichment in the Everglades. In 2001 Everglades Consolidated Report. SFWMD, West Palm Beach, FL. Payne, G., T. Bennett and K. Weaver. 2002. Chapter 5: development of a numeric phosphorus criterion for the Everglades Protection Area. In 2002 Everglades Consolidated Report. SFWMD, West Palm Beach, FL. Payne, G., K. Weaver and S. K. Xue. 2007. Chapter 3C: status of phosphorus and nitrogen in the Everglades Protection Area. In 2007 Everglades Consolidated Report. SFWMD, West Palm Beach, FL. 28 pages. Pietro, K, R. Bearzotti, M. Chimney, G. Germain, N. Iricanin and T. Piccone. 2007. Chapter 5: STA Performance, Compliance and Optimization. In 2007 South Florida Environmental Report. SFWMD and FDEP West Palm Beach, Florida. Pietro, K. C. and D. Ivanoff. 2015. Comparison of long-term phosphorus removal performance of two large-scale constructed wetlands in South Florida, U.S.A. Ecological Engineering 79:143-157. Pietro, K. C., editor. 2016. Chapter 5B: Performance and Operation of the Everglades Stormwater Treatment Areas. In 2016 South Florida Environmental Report. SFWMD and FDEP. West Palm Beach, Florida. 98 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Pisani, O., M. Gao, N. Maie, T. Miyoshi, D. L. Childers and R. Jaffe. 2018. Compositional aspects of herbaceous litter decomposition in the freshwater marshes of the Florida Everglades. Plant Soil 423:87-98. Pollman, C.D. 2012. Modeling sulfate and gambusia mercury relationships in the Everglades (Technical No. SP696). Florida Department of Environmental Protection, Tallahassee, FL. 109 pages. Pollman, C. D. 2014. Mercury cycling in aquatic ecosystems and trophic-state related variables - implications from structural equation modeling. Science of the Total Environment 499:62-73. Pollman, C. D. and D. M. Axelrad. 2014. Mercury Bioaccumulation and Bioaccumulation Factors for Everglades Mosquitofish as Related to Sulfate: A Re-analysis of Julian II (2013). Bulletin of Environmental Contamination and Toxicology 93:509-516. DPI 10.1007/s00128-014-1384-5. Pollman, C. D. 2019. Major Drivers of Mercury Methylation and Cycling in the Everglades: A Synthesis. Chapter 6 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration. D. G. Rumbold et al. (eds.). https://doi.ora/10.1007/978-3-030-32057-7 6. Porter, P. S. and C. A. Sanchez. 1994. Nitrogen in the organic soils of the EAA. Pages 42-61 in Everglades Agricultural Area (EAA): water, soil, crop and environmental management. A. B. Bottcher and F. T. Izuno (editors). University press of Florida. Gainesville, Florida. 319 pages. PTI Environmental Services. 1995. Ecological risks of wading birds of the Everglades in relation to phosphorus reductions in water and mercury bioaccumulation in fishes. Prepared for Sugar Cane Growers Cooperative. PTI Contract C663-01-01. 41 pages. Poulin, B. A., J. N. Ryan, K. L. Nagy, A. Stubbins, T. Dittmar, W. Orem, D. P. Krabbenhoft and G. R. Aiken. 2017. Spatial Dependence of Reduced Sulfur in Everglades Dissolved Organic Matter Controlled by Sulfate Enrichment. Environmental Science and Technology 51(7):3630-3639. https://doi.ora/10.1021/acs.est.6b04142 Radell, M. and B. G. Katz. 1991. Major-ion and selected trace metal chemistry of the Biscayne Aquifer, Southeast Florida. USGS Water Resources Investigations Report 91-4009. Tallahassee, Florida. 18 pages. Ramsar. 2006. Ramsar Convention. The list of wetlands of international importance, version of 7 December 2006. http://www.ramsar.ora/sitelist.pdf. Accessed December 13, 2006. RECOVER. 2007. Comprehensive Everglades Restoration Plan System-wide Performance Measures. Restoration Coordination and Verification, U.S. Army Corps of Engineers, Jacksonville, FL, and SFWMD, West Palm Beach, FL. RECOVER. 2011. Scientific and Technical Knowledge Gained in Everglades Restoration (1999-2009). Restoration Coordination and Verification, U.S. Army Corps of Engineers, Jacksonville, FL, and SFWMD, West Palm Beach, FL. Regier, P., H. Briceno and R. Jaffe. 2016. Long-term environmental drivers of DOC fluxes: Linkages between management, hydrology and climate in a subtropical coastal estuary. Estuarine, Coastal and Shelf Science 182 (A5): 112-122. Richards, J.H. and C.T. Ivey. 2004. Morphological plasticity of Sagittaria lancifolia L. in response to phosphorus. Aquatic Botany 80:53-67. Richards, J. and T. Philippi. 2005. Macrophyte sampling, South Florida R-EMAP project. 2005 dry season. Final Report. Florida International University Department of Biological Sciences. Miami, Florida. 54 pages. Richards, J. H., T. Philippi, P. Kalla and D. Scheidt. 2008. Everglades marsh vegetation: species associations and spatial distributions from R-EMAP Phase III. Florida International University Department of Biological Sciences. Miami, Florida. 37 pages. 99 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Richards, J., L. J. Scinto, Y. Cai, D. Gann and G. Liu. 2017. Regional Environmental Monitoring and Assessment Program (REMAP) IV: Greater Everglades Whole-Ecosystem Monitoring and Assessment. Florida International University. September 28, 2017. Miami, Florida. 208 pages. Richardson, L., K. Keefe, C. Huber, L. Racevskis, G. Reynolds, S. Thourot and I. Miller. 2014. Assessing the value of the Central Everglades Planning Project (CEPP) in Everglades restoration: an ecosystem services approach. Ecological Economics 107:366-377. Rizzardi, K. 2001a. Translating science into law: phosphorus standards in the Everglades. Journal of Land Use and Environmental Litigation 17(1): 149-168. Rizzardi. K. W. 2001b. Alligators and litigators: a recent history of Everglades regulation and litigation. Florida Bar Journal 75(3):18. Rodriguez, A. F., S. Gerber and S. H. Daroub. 2020. Modeling soil subsidence in a subtropical drained peatland. The case ofthe Everglades Agricultural Area. Ecological Modelling 108859. https://doi.ora/10.1016/i. ecolmodel.2019.108859 Roelke, M.E., D.P Schultz, C.F. Facemire, S.F. Sundlof and H.E. Royals. 1991. Mercury Contamination in Florida Panthers (A Report of the Florida Panther Technical Subcommittee to the Florida Panther Interagency Committee). Florida Game and Fresh Water Fish Commission, Eustis, FL. Rumbold, D.G., L.E. Fink, K.A. Laine, S.L. Niemczyk, T. Chandrasekhar, S.D. Wankel and C. Kendall. 2002. Levels of mercury in alligators (Alligator mississippiensis) collected along a transect through the Florida Everglades. Science of The Total Environment 297(1-3):239-252. Rumbold, D. G. 2005. A probabilistic risk assessment of the effects of methylmercury on great egrets and bald eagles foraging at a constructed wetland in south Florida relative to the Everglades. Human and Ecological Risk Assessment 11:365-388. Rumbold D. G., T. R. Lange, D. M. Axelrad and T. D. Atkeson. 2008. Ecological risk of methylmercury in Everglades National Park, Florida, USA. Ecotoxicology 17:632-641. Rumbold, D. G. 2019a. A Causal Analysis for the Dominant Factor in the Extreme Geographic and Temporal Variability in Mercury Biomagnification in the Everglades. Chapter 3 In Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration, https://doi.ora/10.1007/978-3-030-32057-7 3. Rumbold, D. G. 2019b. Regional-Scale Ecological Risk Assessment of Mercury in the Everglades and South Florida. Chapter 10 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration. https://doi. ora/10.1007/978-3-030-32057-7 10. Sargeant, B., J.C. Trexler, E.E. Gaiser. 2010. Biotic and Abiotic Determinants of Intermediate-Consumer Trophic Diversity in the Florida Everglades. Marine & Freshwater Research 61:11-22. Sargeant, B. L., E. E. Gaiser, and J. C. Trexler. 2011. Indirect and direct controls of macroinvertebrates and small fish by abiotic factors and trophic interactions in the Florida Everglades. Journal of Freshwater Biology DPI: 10.1111/i. 1365-2427.2011,02663.x Scheidt, D. J., D. R. Walker, R. G. Rice and M. D. Flora. 1985. Diel dissolved oxygen levels under experimental nutrient loading conditions in Shark Slough, Everglades National Park, Florida. Florida Scientist 48: supplement 1:36. Scheidt, D. J., M. D. Flora and D. R. Walker. 1989. Water quality management for Everglades National Park. Pages 377-390 in Wetlands: Concerns and Successes. American Water Resources Association. Scheidt, D., J. Stober, R. Jones and K. Thornton. 2000. South Florida ecosystem assessment: water management, soil loss, eutrophication and habitat. EPA 904-R-00-003. USEPA Region 4, Athens, GA. 38 pages. 100 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Scheidt, D.J., and P.I. Kalla. 2007. Everglades ecosystem assessment: water management and quality, eutrophication, mercury contamination, soils, and habitat: monitoring for adaptive management: a REMAP status report. USEPA Region 4, Athens, GA. EPA904-R-07-001. 104 pages. Schueneman, T. J. and C. A. Sanchez. 1994. Vegetable Production in the EAA. Pages 238-277 in Everglades Agricultural Area (EAA): Water, Soil, Crop, and Environmental Management. University Press of Florida. Gainesville, Florida. Schueneman, T. J. 2001. Characterization of sulfur sources in the EAA. Soil and Crop Sciences Society of Florida Proceedings 60:49-52. Shih, S. F., B. Glaz, and R. E. Barnes Jr. 1998. Subsidence of organic soils in the Everglades agricultural area during the past 19 years. Soil and Crop Science Society of Florida Proceedings 57:20-29. Sklar, F. H., M. Chimney, S. Newman, P. McCormick, D. Gawlik, S Miao, C. McVoy, W. Said, J. Newman, C. Coronado, G. Crozier, M. Korvela and K. Ken Rutchey. 2005a. The ecological-societal underpinnings of Everglades restoration. Frontiers in Ecology and the Environment 3:161-169. Sklar, F., K. Rutchey, S. Hagerthey, M. Cook, S. Newman, S. Miao, C. Coronado-Molina, J. Leeds, L. Bauman, J. Newman, M. Korvela, R. Wanvestraut and A. Gottlieb. 2005b. Chapter 6: Ecology of the Everglades Protection Area. In 2005 South Florida Environmental Report. SFWMD and FDEP West Palm Beach, Florida. Sklar, F., M. Cook, E. Call, R. Shuford, M. Kobza, R. Johnson, S. Miao, M. Korvela, C. Coronado, L. Bauman, J. Leeds, B. Garrett, J. Newman, E. Cline, S. Newman, K. Rutchey and C. McVoy. 2006. Chapter 6: Ecology of the Everglades Protection Area. In 2006 South Florida Environmental Report. SFWMD and FDEP. West Palm Beach, Florida. 65 pages. Sklar, F., T. Dreschel and K. Warren. 2010. Chapter 6: Ecology of the Everglades Protection Area. In 2010 South Florida Environmental Report. SFWMD and FDEP. West Palm Beach, Florida. 91 pages. Sklar, F. H., J. F. Meeder, T. G. Troxler, T. Dreschel, S. E. Davis and P. L. Ruiz. 2019. The Everglades: At the Forefront of Transition. Coasts and Estuaries 277-292. https://doi.ora/l0.1016/8978-0.12-814003-1.00016-2. Smith, S. M., D. E. Gawlik, K. Rutchey, G. E. Crozier, and S. Gray. 2003. Assessing drought-related ecological risk in the Florida Everglades. Environmental Management 68:355-366. Smolders, A. J. P., L. P. Lamers, C., H. Lucassen, G. van der Velde and J. G. Roelofs. 2006. Internal eutrophication: how it works and what to do about it - a review. Chemistry and Ecology 22(2):93-111. Snyder G. H. 2005. Everglades agricultural area soil subsidence and land use projections. Soil Crop Sci. Soc. Fla. Proc. 64:44-51. SFWMD. 1977. Water Use and Supply Development Plan. Volume IIIA. Lower East Coast Planning Area Technical Exhibits A - H. South Florida Water Management District. West Palm Beach, Florida. SFWMD. 1992a. Surface Water Improvement and Management Plan for the Everglades. Supporting Information Document. West Palm Beach, Florida. 470 pages. SFWMD. 1992b. Surface Water Improvement and Management Plan for the Everglades. Appendix E. West Palm Beach, Florida. 60 pages. SFWMD. 2007. South Florida Environmental Report. March 1, 2007. Executive Summary. West Palm Beach, Florida. SFWMD. 2012. Restoration Strategies Regional Water Quality Plan. April 27, 2012. Final Version. West Palm Beach, Florida. 101 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) SFWMD. 2014. Weekly STA Performance Summary. September 09, 2014. SFWMD. 2016a. South Florida Environmental Report Highlights, https://www.sfwmd.aov/sites/default/files/ documents/2016 sfer hiahliahts.pdf SFWMD. 2016b. South Florida environmental report: Appendix 3-1: annual permit report for the Everglades Stormwater Treatment Areas, permit report (May 1, 2014 - April 30, 2015, permit 0311207). Spalding, M.G., P.C. Frederick, H.C. McGill, S.N. Bouton and L.R. McDowell. 2000. Methylmercury accumulation in tissues and its effects on growth and appetite in captive great egrets. Journal of Wildlife Diseases 36(3):411-422. Stephens, J. C. 1956. Subsidence of Organic Soils in the Florida Everglades. Soil Science Society Proceeding, pages 77-80. Stephens, J. C. and L. Johnson. 1951. Subsidence of Peat Soils in the Everglades Region of Florida. United States Department of Agriculture Soil Conservation Service. 47 pages. Stevens, D. L., Jr. 1997. Variable Density Grid-based Sampling Designs for Continuous Spatial Populations. Environmetrics 8:167-195. Stevens, D.L., Jr. and A.R. Olsen. 2004. Spatially-balanced sampling of natural resources. Journal of the American Statistical Association 99(465):262-278. DPI: 10.1198/016214504000000250. Stober, Q. J., R. D. Jones and D. J. Scheidt. 1995. Ultra-trace level mercury in the Everglades Ecosystem: A Multi- media Pilot Study. Water, Air and Soil Pollution 80:991-1001. Stober, J., D. Scheidt, R. Jones, K. Thornton, R. Ambrose, and D. France. 1998. South Florida Ecosystem Assessment. Monitoring for Adaptive Management: Implications for Ecosystem Restoration. Final Technical Report - Phase I. EPA-904-R-98-002. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 478 pages. Stober, Q. J., K. Thornton, R. Jones, J. Richards, C. Ivey, R. Welch, M. Madden, J. Trexler, E. Gaiser, D. Scheidt and S. Rathbun. 2001a. South Florida Ecosystem Assessment: Phase l/ll - Everglades Stressor Interactions: Hydropatterns, Eutrophication, Habitat Alteration, and Mercury Contamination (Summary). EPA904-R-01-002. September 2001. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 63 pages. Stober, Q.J., K. Thornton, R. Jones, J. Richards, C. Ivey, R. Welch, M. Madden, J. Trexler, E. Gaiser, D. Scheidt, and S. Rathbun. 2001b. South Florida ecosystem assessment: Phase l/ll (technical report) - Everglades stressor interactions: Hydropatterns, eutrophication, habitat alteration, and mercury contamination. EPA904-R-01-003. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 1625 pages. Storch, W. V. 1971. Pollution control aspects of water management. Central and Southern Florida Flood Control District. SFWMD Report DRE-09. West Palm Beach, Florida. 19 pages. Surratt, D., M. Waldon, M. Harwell and N. Aumen. 2008. Time-series and spatial tracking of polluted canal water intrusion into wetlands of a national wildlife refuge in Florida, USA. Wetlands 28:176-183. Surratt, D., D. Shinde and N. Aumen. 2012. Recent cattail expansion and possible relationships to water management: changes in upper Taylor Slough (Everglades National Park, USA). Environmental Management 49:720- 733. Swift, D. R. 1981. Preliminary Investigations of Periphyton and Water Quality Relationships in the Everglades Water Conservation Areas. South Florida Water Management District Technical Publication 81-5. West Palm Beach, Florida. 83 pages plus appendices. 102 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Swift, D. and R. B. Nicholas. 1987. Periphyton and Water Quality Relationships in the Everglades Water Conservation Areas 1978-1982. South Florida Water Management District Technical Publication 87-2. West Palm Beach, Florida. 44 pages plus appendices. Thornton, K. W., Saul, G. E. and Hyatt, D. E. 1994. Environmental Monitoring and Assessment Program Assessment Framework. USEPA Report EPA/620/R-94/016. Research Triangle Park, North Carolina. 47 pages. Trexler, J. C., E. E. Gaiser, J. S. Kominoski, and J. Sanchez. 2014. The role of periphyton mats in consumer community structure and function in calcareous wetlands: Lessons from the Everglades. Pages 155 -179 In Microbiology of the Everglades Ecosystem, J. A. Entry, A. D. Gottlieb, K. Jayachandrahan and A. Ogram (editors). Science Publishers, CRC Press. Troxler, T. G. and J. H. Richards. 2009. d13C, d15N, carbon, nitrogen and phosphorus as indicators of plant ecophysiology and organic matter pathways in Everglades deep slough, Florida, USA. Aquatic Botany 91:157-165. Turner, A. M., J. C. Trexler, F. Jordan, S. J. Slack, P. Geddes, J. H. Chick and W. F. Loftus. 1999. Targeting Ecosystem Features for Conservation: Standing crops in the Florida Everglades. Conservation Biology 13(4):898- 911. USACE. 1996. Florida's Everglades Program. Everglades Construction Project. Final programmatic environmental impact statement. Jacksonville, Florida. USACE and SFWMD. 1999. Rescuing an Endangered Ecosystem: The Plan to Restore America's Everglades. The Central and Southern Florida Project Restudy. 28 pages. USACE. 2014. Comprehensive Everglades Restoration Plan. Central Everglades Planning Project. Final integrated project implementation report and programmatic environmental impact statement. Jacksonville, Florida. 366 pages. USACE and USDOI. 2015. Comprehensive Everglades Restoration Plan (CERP), Central and Southern Florida Project, 2015 Report to Congress. US Army Corps of Engineers, Jacksonville, Florida, and US Department of the Interior, Washington, D. C. 94 pages. USACE and USDOI. 2020. Draft Comprehensive Everglades Restoration Plan (CERP), Central and Southern Florida Project, 2020 Report to Congress. US Army Corps of Engineers, Jacksonville, Florida, and US Department of the Interior, Washington, D. C. 133 pages. USDOI. 1971. Appraisal of Water Quality Needs and Criteria for Everglades National Park. USDI National Park Service, Washington, D. C. June 1971. 50 pages. U.S. District Court. 1992. Consent Decree. USA vs. SFWMD and FDER. Case No. 88-1886-CIV-HOEVELOR. Southern District of Florida. USEPA. 1986. Quality Criteria for Water ("Gold Book"). Office of Water, Regulations and Standards. Washington, DC. EPA 440/5-86-001. USEPA. 1993. R-EMAP: Regional Environmental Monitoring and Assessment Program. EPA/625/R-93/012. September 1993. USEPA Office of Research and Development. Washington, D. C. USEPA. 1995. Environmental Monitoring and Assessment Program (EMAP) Cumulative Bibliography. USEPA, Office of Research and Development. EPA/620/R-95/006. Research Triangle Park, North Carolina. 44 pages. USEPA. 1997. Mercury study report to Congress. Volume VI: an ecological assessment for anthropogenic mercury emissions in the United States. USEPA Office of Air Quality Planning & Standards and Office of Research and Development. EPA-452/R-97-008. 103 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) USEPA. 2000. Ambient water quality criteria recommendations: information supporting the development of state and tribal nutrient criteria. Wetlands in nutrient ecoregion XIII. Publication EPA 822-B-00-023. Office of Water, Office of Science and Technology, Health and Ecological Criteria Division. Washington, D. C. 78 pages. USEPA. 2001. Water Quality Criterion for the Protection of Human Health: Methylmercury. EPA-823-R-1005 01-001. USEPA, Washington, DC. USEPA. 2002. Guidance for Quality Assurance Project Plans (EPA QA/G-5), EPA/240/R-02/009, Office of Environmental Information. Washington, D. C. USEPA. 2014a. Everglades ecosystem assessment phase IV. Miami, Florida. September 23-29, 2013. SESD Project Identification Number: 13-0513. USEPA Region 4 SESD, Athens, GA. 22 pages. USEPA. 2014b. Everglades ecosystem assessment phase IV quality assurance project plan. USEPA Region 4 SESD, Athens, GA. 135 pages. USEPA. 2015. National rivers and streams assessment 2013-14: Quality assurance project plan. EPA-841-B-12-007. USEPA, Office of Water, Washington, DC. USGS, 2019, Specific conductance: U.S. Geological Survey Techniques and Methods, book 9, chapter A6.3. 15 pages. https://doi.ora/10.3133/tm9A6.3. Vijayaraghavan, K. and C. D. Pollman. 2019. Mercury Emission Sources and Contributions of Atmospheric Deposition to the Everglades. Chapter 5 in Mercury and the Everglades. A Synthesis and Model for Complex Ecosystem Restoration. C. D. Pollman et al. (eds.). https://doi.ora/10.1007/978-3-030-20070-1 5. Walker, W. W. Jr. 2000. Interim phosphorus standards for the Everglades. Pages B-34 to B-46 in A Nutrient Criteria Technical Guidance Document. Lakes and Reservoirs. First Edition. EPA822-B00-001. USEPA Office of Water, Office of Science and Technology. Washington, D. C. Walker, W. W. Jr. and R. H. Kadlec. 1996. A model for simulating phosphorus concentrations in waters and soils downstream of Everglades Stormwater Treatment Areas. 108 pages. Waller, B. 1982. Water Quality Characteristics of Everglades National Park, 1959-1977, with reference to the effects of Water Management. USGS Water Resources Investigation 82-34. Tallahassee, Florida. 51 pages. Wang, H., Waldon, M. G., Meselhe, E.A., Arceneaux, J. C., Chen, C., Harwell, M. C. 2009. Surface water sulfate dynamics in the northern Florida Everglades, USA. Journal of Environment Quality 38:1-8. Wang, H., E. A. Meselhe, M. G. Waldon, M. C. Harwell and C. Chen. 2012. Compartment-based hydrodynamics and water quality modeling of a Northern Everglades Wetland, Florida, USA. Ecological Modelling 247:273-285. Wasserstein, R. L., A. L. Schirm and N. A. Lazar. 2019. Moving to a World Beyond "p < 0.05" The American Statistician, 73:sup1, 1-19, DOI:1Q. 1080/00031305.2019.1583913 Weaver, K. 2004. Everglades Marsh Dissolved Oxygen Site Specific Alternative Criterion Technical Support Document. FDEP Tallahassee, Florida. 61 pages. Weaver, K., G. Payne and S. Xue. 2007. Chapter 3A: Status of Water Quality in the Everglades Protection Area. In 2007 South Florida Environmental Report. SFWMD and FDEP. West Palm Beach, Florida. Welch, R. and M. Madden. 2000. Aerial photo vegetation assessment in the Everglades ecosystem. Final report. Submitted to U. S. Department of Interior, October 19, 2000. Center for Remote Sensing and Mapping Sciences, Department of Geography, University of Georgia. Athens, Georgia. 285 pages. 104 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Wetzel, P. R., S. E. Davis, T. van Lent, S. M. Davis, H. Henriquez. 2017. Science synthesis for management as a way to advance ecosystem restoration: evaluation of restoration scenarios for the Florida Everglades. Restoration Ecology 25(S1): S4-S17. White, J. R. and K. R. Reddy. 2003. Nitrification and denitrification rates of Everglades wetland soils along a phosphorus-impacted gradient. Journal of Environmental Quality 32:2436-2443. Williams, A. J., and J. C. Trexler. 2006. A preliminary analysis of the correlation of food-web characteristics with hydrology and nutrient gradients in the southern Everglades. Hydrobiologia 569:493-504. Wright, A. L., K. R. Reddy and S. Newman. 2009. Microbial Indicators of Eutrophication in Everglades. Soil Sci. Soc. Am. J. 73:1597-1603. doi: 10.2136/sssai2009.0083. Wright, A. L. and G. H. Snyder. 2016. Soil Subsidence in the Everglades Agricultural Area. Everglades Research and Education Center, Soil and Water Science Department. Issue SL551. 4 pages. Xue, S. K. 2016. Appendix 3A-5: Water Year 2015 and Five-Year (Water Years 2011-2015) Annual Flows and Total Phosphorus Loads and Concentrations by Structure and Area. In 2016 South Florida Environmental Report. SFWMD and FDEP West Palm Beach, Florida. Yamashita, Y., L. J. Scinto, N. Maie and R. Jaffe. 2010. Dissolved Organic Matter Characteristics Across a Subtropical Wetland's Landscape: Application of Optical Properties in the Assessment of Environmental Dynamics. Ecosystems 13:1006-1019. Ye, R., A L. Wright, J. M. McCray, K. R. Reddy and L. Young. 2009. Sulfur-induced changes in phosphorus distribution in Everglades Agricultural Area soils. Nutrient Cycling in Agrosystems. DOI 10.1007/s10705-009-9319-y. Ye, R., A L. Wright, W. H. Orem and J. M. McCray. 2010. Sulfur distribution and transformations in Everglades Agricultural Area soil influenced by sulfur amendment. Soil Science 175(6):263-269. Ye, R., A.L. Wright, and J.M. McCray. 2011a. Seasonal changes in nutrient availability for sulfur-amended Everglades soil under sugarcane. Journal of Plant Nutrition 43:2095-2113. Ye, R., A.L. Wright, and J.M. McCray. 2011b. Microbial response of a calcareous Histosolto sulfur amendment. Soil Science 176:479-486. Yoder, L., R. R. Chowdhury and C. Hauck. 2020. Watershed restoration in the Florida Everglades: Agricultural water management and long-term trends in nutrient outcomes in the Everglades Agricultural Area 302:107070. https://doi. ora/10.1016/i.aaee.2020.107070 Zabala, J., J. C. Trexler, N. Jayasena, and P. Frederick. 2020. Early Breeding Failure in Birds Due to Environmental Toxins: A Potentially Powerful but Hidden Effect of Contamination. Environmental Science and Technology, https:// dx. doi. ora/10.1021/acs. est. 0c04098 Zhang, Y, D.J. Jacob, H.M. Horowitz, L. Chen, H.M.Amos, D.P Krabbenhoft, F. Slemr, V.L. St. Louis, and E. Sunderland. 2016. Observed decrease in atmospheric mercury explained by global decline in anthropogenic emissions. Proc. National Academy of Sciences 113f3V526-531. doi: 10.1073/pnas.1516312113. Zhao, H. and T. Piccone. 2020. Large scale constructed wetlands for phosphorus removal, an effective nonpoint source pollution treatment technology. Ecological Engineering 145:105711. Zheng, S., B. Gu, Q. Zhou and Y. Li. 2013. Variations of mercury in the inflow and outflow of a constructed treatment wetland in south Florida, USA. Ecological Engineering 61:419-425. 105 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Zweig, C. L. and Newman, S. 2015. Using landscape context to map invasive species with medium-resolution satellite imagery. Restoration Ecology 23: 524-530. doi: 10.1111/rec. 12214. ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) EVERGLADES REMAP PROGRAM REPORTS AND PUBLICATIONS USEPA Everglades REMAP: https://www.epa.gov/everglades/environmental-monitoring-everglades Abbey-Lee, Robin N., Evelyn N. Gaiserand J. C. Trexler. 2013. Relative roles of dispersal dynamics and competition in determining the isotopic niche breadth of a wetland fish. Freshwater Biology 58:780-792. doi:10.1111/fwb. 12084. Axelrad, Donald M., Thomas D. Atkeson, Ted Lange, Curtis D. Pollman, Cynthia C. Gilmour, William H. Orem, Irving A. Mendelssohn, Peter C. Frederick, David P. Krabbenhoft, George R. Aiken, Darren G. Rumbold, Daniel J. Scheidt, and Peter I. Kalla. 2007. Chapter 3B: Mercury Monitoring, Research and Environmental Assessment in South Florida. In 2007 South Florida Environmental Report. SFWMD and FDEP 57 pages. Axelrad, Donald M., Ted Lange, Mark Gabriel, Thomas D. Atkeson, Curtis D. Pollman, William H. Orem, Daniel J. Scheidt, P. I. Kalla, P. C. Frederick and C. C. Gilmour. 2008. Chapter 3B: Mercury and Sulfur Monitoring, Research and Environmental Assessment in South Florida. In 2008 South Florida Environmental Report. SFWMD and FDEP. Cai, Yong, Rudolf Jaffe, and Ronald Jones. 1997. Ethylmercury in the Soils and Sediments of the Florida Everglades. Environmental Science and Technology 31(1):302-305. Doren, Robert, John C. Volin and Jennifer H. Richards. 2009. Invasive exotic plant indicators for ecosystem restoration: An example from the Everglades restoration program. Ecological Indicators 9(6) Supplement: S29-S36. Gaiser, Evelyn E. Paul V. McCormick, Scot E. Hagerthey and Andrew D. Gottlieb. 2011. Landscape Patterns of Periphyton in the Florida Everglades. Critical Reviews in Environmental Science and Technology 41(S1):92-120. Ivey, C. and J.H. Richards. 2001. Genotypic diversity and clonal structure of Everglades sawgrass, Cladium jamaicense (Cyperaceae). International Journal of Plant Science 162:1327-1335. Ivey, C. and J.H. Richards. 2001. Genetic diversity of Everglades sawgrass (Cladium jamaicense, Cyperaceae). International Journal of Plant Science 162: 817-825. Jiang, P., G. Liu, W. Cui and Y. Cai. 2018. Geochemical modeling of mercury speciation in surface water and implications on mercury cycling in the Everglades wetland. Science of the Total Environment 640-641:454-465. John, A., H. R. Fuentes and F. George. 2021. Characterization of the water retention curves of Everglades wetland soils. Geoderma 381:114724. https://www.sciencedirect.eom/science/iournal/00167061/381/supp/C Kalla, P. I. and D. J. Scheidt. 2017. Everglades ecosystem assessment - Phase IV, 2014: Data reduction and initial synthesis. United States Environmental Protection Agency, Science and Ecosystem Support Division. SESD Project 14-0380. Athens, Georgia. 58 pages. Kalla, P., M. Cyterski, J. Minucci and D. Scheidt. 2019. Modelling of mercury bioaccumulation in mosquitofish (Gambusia sp.) for the Everglades Regional Monitoring and Assessment Program. USEPA Laboratory Services and Applied Science Division, Athens, GA. LSASD 20-0085. 35 pages. Kalla, P., M. Cyterski, D. Scheidt and J. Minucci. 2021. Spatiotemporal Effects of Interacting Water Quality Constituents on Mercury in a Common Prey Fish in a Large, Perturbed, Subtropical Wetland. Science of the Total Environment 792 (2021) 148321. https://doi.Org/10.1016/j.scitotenv.2021.148321. Li, Yanbin, Yuxiang Mao, Guangliang Liu, Georgio Tachiev, David Roelant, Xinbin Feng and Yong Cai. 2011. Degradation of Methylmercury and Its Effects on Mercury Distribution and Cycling in the Florida Everglades. Environmental Science and Technology 44:6661-6666. 107 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Li, Yanbin, Zhiwei Duan, Guangliang Liu, Peter Kalla, Daniel Scheidt and Yong Cai. 2015. Evaluation of the possible sources and controlling factors of toxic metals/metalloids in the Florida Everglades and their potential risk of exposure. Environmental Science & Technology. DPI: 10.1021/acs.est.5b01638. Liu, G., Y. Cai, P. Kalla, D. Scheidt, J. Richards, L. Scinto, E. Gaiserand C.Appleby. 2008a. Mercury mass budget estimates and cycling seasonality in the Florida Everglades. Environmental Science and Technology 42:1954-1960. Liu, G., Y. Cai, T. Philippi, P. Kalla, D. Scheidt, L. Scinto, J. Richards and J. Trexler. 2008b. Distribution of total and methylmercury in different ecosystem compartments in the Everglades: Implications for mercury accumulation. Environmental Pollution 153:257-265. Liu, G., Yong Cai, Yuxiang Mao, Daniel Scheidt, Peter Kalla, Jennifer Richards, Leonard Scinto, Georgio Tachiev, David Roelant and Charlie Appleby. 2009. Spatial Variability in Mercury Cycling and Relevant Biogeochemical Controls in the Florida Everglades. Environmental Science and Technology 43(12):4361-4366. DPI: 10.1021/ es803665c Liu, G., G. M. Naja, P. Kalla, D. Scheidt, E. Gaiser, and Y. Cai. 2011. Legacy and fate of mercury and methylmercury in the Florida Everglades. Environmental Science and Technology, 45(2):496-501. Supporting information S1-S41. Liu, Guangliang, Yong Cai, Georgio Tachiev and Leone Lagos. 2015. Mercury Mass Budget Estimates and Cycling in the Florida Everglades. Chapter 4 in Microbiology of the Everglades Ecosystem. James A. Entry, Andrew D. Gottlieb, Krish Jayachandran, Andrew Pgram (editors). CRC Press. Prem, William, Cynthia Gilmour, Donald Axelrad, David Krabbenhoft, Daniel Scheidt, Peter Kalla, P, McCormick, Mark Gabriel and George Aiken. 2011. Sulfur in the South Florida Ecosystem: Distribution, Sources, Biogeochemistry, Impacts, and Management for Restoration. Critical Reviews in Environmental Science and Technology 41 (S1 ):249- 288. Psborne, T. Z., S. Newman; D. J. Scheidt, P. I. Kalla, G. L. Bruland, M. J. Cohen, L. J. Scinto and L. R. Ellis. 2011. Landscape Patterns of Significant Soil Nutrients and Contaminants in the Greater Everglades Ecosystem: Past, Present, and Future. Critical Reviews in Environmental Science and Technology 41 (S1): 121 -148. Richards, J.H. and C.T. Ivey. 2004. Morphological plasticity of Sagittaria lancifolia L. in response to phosphorus. Aquatic Botany 80:53-67. Richards, Jennifer and Tom Philippi. 2005. Macrophyte sampling, South Florida R-EMAP project. 2005 dry season. Final Report. Florida International University Department of Biological Sciences. Miami, Florida. 54 pages. Richards, J. H., T. Philippi, P. Kalla and D. Scheidt. 2008. Everglades marsh vegetation: species associations and spatial distributions from R-EMAP Phase III. Florida International University Department of Biological Sciences. Richards, Jennifer, Leonard J. Scinto, Yong Cai, Daniel Gann and Guangliang Liu. 2017. Regional Environmental Monitoring and Assessment Program (REMAP) IV: Greater Everglades Whole-Ecosystem Monitoring and Assessment. September 28, 2017. Florida International University. Miami, Florida. 208 pages. Sargeant, B., J.C. Trexler, E.E. Gaiser. 2010. Biotic and Abiotic Determinants of Intermediate-Consumer Trophic Diversity in the Florida Everglades. Marine & Freshwater Research 61:11-22. Sargeant, Brooke L., Evelyn E. Gaiser, and Joel C. Trexler. 2011. Indirect and direct controls of macroinvertebrates and small fish by abiotic factors and trophic interactions in the Florida Everglades. Journal of Freshwater Biology DPI: 10.1111 /j. 1365-2427.2011.02663.x Scheidt, Daniel, J. Stober, R. Jones and K. Thornton. 2000. South Florida ecosystem assessment: water management, soil loss, eutrophication and habitat. EPA 904-R-00-003. USEPA Region 4. Atlanta, Georgia. 38 pages. 108 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Scheidt, D. J. and P. I. Kalla. 2007. Everglades ecosystem assessment: water management, water quality, eutrophication, mercury contamination, soils and habitat. Monitoring for adaptive management: A R-EMAP status report. EPA904-R-07-001. USEPA Region 4. Athens, Georgia. 104 pages. Stober, Q. J., R. D. Jones and D. J. Scheidt. 1995. Ultra-trace level mercury in the Everglades Ecosystem: A Multi- media Pilot Study. Water, Air and Soil Pollution 80:991-1001. Stober, Jerry, Daniel Scheidt, Ron Jones, Kent Thornton, Robert Ambrose, and Danny France. 1996. South Florida Ecosystem Assessment. Monitoring for Adaptive Management: Implications for Ecosystem Restoration. Interim Report. EPA-904-R-96-008. USEPA Region 4 SESD and ORD. Athens, Georgia. 26 pages. Stober, Jerry, Daniel Scheidt, Ron Jones, Kent Thornton, Robert Ambrose, and Danny France. 1998. South Florida Ecosystem Assessment. Monitoring for Adaptive Management: Implications for Ecosystem Restoration. Final Technical Report - Phase I. EPA-904-R-98-002. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 478 pages. Stober, Q. J., K. Thornton, R. Jones, J. Richards, C. Ivey, R. Welch, M. Madden, J. Trexler, E. Gaiser, D. Scheidt and S. Rathbun. 2001. South Florida Ecosystem Assessment: Phase l/ll- Everglades Stressor Interactions: Hydropatterns, Eutrophication, Habitat Alteration, and Mercury Contamination (Summary). EPA904-R-01-002. September 2001. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 63 pages. Stober, Q. J., K. Thornton, R. Jones, J. Richards, C. Ivey, R. Welch, M. Madden, J. Trexler, E. Gaiser, D. Scheidt and S. Rathbun. 2001. South Florida Ecosystem Assessment: Phase l/ll (Technical Report)- Everglades Stressor Interactions: Hydropatterns, Eutrophication, Habitat Alteration, and Mercury Contamination. EPA904-R-01-003. September 2001. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 1625 pages. Tarpey, Thaddeus and Christopher T. Ivey. 2006. Allometric Extension for Multivariate Regression. Journal of Data Science 4:479-495. Trexler, J. C., E. E. Gaiser, J. S. Kominoski, and J. Sanchez. 2015. The role of periphyton mats in consumer community structure and function in calcareous wetlands: Lessons from the Everglades. Pages 155-179 In: Microbiology of the Everglades Ecosystem. J. A. Entry, A. D. Gottlieb, K. Jayachandrahan and A. Ogram (editors). Science Publishers, CRC Press. USEPA. 1993. R-EMAP: Regional Environmental Monitoring and Assessment Program. EPA/625/R-93/012. September 1993. USEPA Office of Research and Development. Washington, D. C. 96 pages. USEPA. 1997. Mercury study report to Congress. Volume VI: an ecological assessment for anthropogenic mercury emissions in the United States. USEPA Office of Air Quality Planning & Standards and Office of Research and Development. EPA-452/R-97-008. USEPA. 2014a. Everglades Ecosystem Assessment Phase IV. Miami, Florida. September 23-29, 2013. SESD Project Identification Number: 13-0513. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 22 pages. USEPA. 2014b. Everglades ecosystem assessment phase IV quality assurance project plan. USEPA Region 4 Science and Ecosystem Support Division. Athens, Georgia. 135 pages. Welch, R. and M. Madden. 2000. Aerial photo vegetation assessment in the Everglades ecosystem. Final report. Submitted to U. S. Department of Interior, October 19, 2000. Center for Remote Sensing and Mapping Sciences, Department of Geography, University of Georgia. Athens, Georgia. 285 pages. Yamashita, Youhei, Leonard J. Scinto, Nagamitsu Maie and Rudolf Jaffe. 2010. Dissolved Organic Matter Characteristics Across a Subtropical Wetland's Landscape: Application of Optical Properties in the Assessment of Environmental Dynamics. Ecosystems 13:1006-1019. 109 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix I. Everglades REMAP 2014 Measurements and Analytes by Media. SURFACE WATER (SW) FLOC (FC) Depth Thickness Temperature Ash Free Dry Weight Dissolved Oxygen Bulk Density pH Water Content Conductivity Total Carbon Turbidity Total Nitrogen Total Phosphorus Total Phosphorus Soluble Reactive Phosphorus Total Mercury Total Nitrogen Methyl mercury Total Inorganic Nitrogen Chlorophyll a Total Organic Nitrogen Filtered Ammonia PERIPHYTON [PB (benthic); Filtered Nitrate PC = PF (floating) + PE Filtered Nitrite Ash Free Dry Weight Filtered Nitrate-Nitrite Bulk Density Dissolved Organic Carbon Total Carbon Total Organic Carbon Total Nitrogen Chlorophyll a Total Phosphorus Sulfate Methyl mercury Chloride Total Mercury Methyl mercury Chlorophyll a Total Mercury VEGETATION BOTTOM WATER (BW) Sawgrass Total Mercury Sulfide Sawgrass Methylmercury Sawgrass Total Phosphorus SOIL (SD) Sawgrass Total Carbon Type Sawgrass Total Nitrogen Thickness Vegetation Type pH Dominant Macrophyte Ash Free Dry Weight Cattail Presence Bulk Density Water Content MOSQUITOFISH Total Phosphorus Total Mercury Total Carbon Total Nitrogen Total Mercury Methyl mercury 110 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix II. Everglades REMAP Median Values of Selected Wet Season Parameters. Analyte Medium 1995 Year 2005 2014 Units* Total Phosphorus surface water 11.6 a 7.50 b 6.6 c ug/L Total Phosphorus soil 325 a 390b 390 b mg/kg Total Mercury mosquitofish 141a 87 b 33.5 c ng/g Total Mercury surface water 1.9 a 2.2 b 1.6° ng/L Methylmercury surface water 0.28 a 0.21 b 0.1 c ng/L Total Mercury water column periphyton 18 a 19a ng/g Methylmercury water column periphyton 1.66 a 1.8a ng/g Total Mercury floe 130 a 120b ng/g Methylmercury floe 2.95 a 2.55 a ng/g Total Mercury benthic periphyton 9.70 a 24 a ng/g Methylmercury benthic periphyton 0.466 a 0.505b ng/g Total Mercury soil 120 a 140a 150a ug/kg Methylmercury soil 0.44 a 0.49 a 0.77b ng/g ab 0 Distributions with medians having different letters are different (P < 0.05). Units are nanograms per gram [(ng/g), parts per billion], micrograms per liter [(ug/L), parts per billion], nanograms per liter [(ng/L), parts per trillion], milligrams per kilogram [(mg/kg), parts per million], micrograms per kilogram [(ug/kg), parts per billion], milligrams per liter [(mg/L), parts per million]. 111 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix III. Water quality characteristics by subarea for REMAP dry season sampling events 1995, 1996, 1999 and 2005. Subarea Analyte Number of Samples Minimum Maximum Median Mean ± Standard Deviation Standard Error 95% Confidence Interval DRY SEASON LNWR Chloride (mg/L) 11 21 75 39.0 37.7 ± 16.4 4.9 11.0 WCA2 11 73 200 140.0 132.6 ±40.7 12.3 27.3 WCA3N 9 38 72 58.0 55.9 ± 11.2 3.7 8.6 WCA3S 48 15 84 50.0 50.4 ± 18.4 2.7 5.3 ENP 16 35 110 60.0 60.8 ± 20.0 5.0 10.7 LNWR Conductivity (umhos/cm) 30 78.1 717 183.5 226.5 ± 148.9 27.2 55.6 WCA2 27 539 1610 894.0 941.6 ±243.2 46.8 96.2 WCA3N 23 148 1281 520.0 607.7 ± 273.1 56.9 118.1 WCA3S 111 231 807 527.0 520.4 ± 132.4 12.6 24.9 ENP 58 307 823 587.5 572.7 ± 119.4 15.7 31.4 LNWR PH 30 5.33 7.33 6.3 6.4 ±0.5 0.1 0.2 WCA2 27 6.85 8.06 7.5 7.53 ±0.31 0.1 0.1 WCA3N 23 6.91 8.05 7.4 7.4 ±0.3 0.1 0.1 WCA3S 111 6.46 7.88 7.3 7.3 ±0.3 0.0 0.0 ENP 58 6.79 8.25 7.4 7.4 ±0.3 0.0 0.1 LNWR Sulfate (mg/L) 30 0.039 14 2.0 2.49 ± 3.0 0.6 1.1 WCA2 27 4.3 73 38.0 37.6 ± 19.3 3.7 7.7 WCA3N 23 2 42 7.0 12.5 ± 12.8 2.7 5.5 WCA3S 111 0.029 36 2.0 4.0 ±6.5 0.6 1.2 ENP 59 0.16 83 2.0 4.46 ± 10.9 1.4 2.8 LNWR Total organic carbon (mg/L) 30 15.12 57.14 25.0 28.5 ±9.8 1.8 3.7 WCA2 27 25.81 50.38 38.6 37.5 ± 7.9 1.5 3.1 WCA3N 23 18 62.98 25.4 28.4 ± 11.0 2.3 4.8 WCA3S 105 11 54.68 22.9 23.8 ± 7.1 0.7 1.4 ENP 59 11.53 49 25.0 25.1 ± 7.7 1.0 2.0 112 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix III (continued). Water quality characteristics by subarea for REMAP wet season sampling events 1995, 1996, 1999, 2005 and 2014. Subarea Analyte Number of Samples Minimum Maximum Median Mean ± Standard Deviation Standard Error 95% Confidence Interval WET SEASON LNWR Chloride (mg/L) 36 5.9 140 19.0 29.5 ± 31.4 5.2 10.6 WCA2 34 46 260 89.5 101.9 ±48.7 8.3 17.0 WCA3N 46 12 97 38.5 45.7 ± 20.2 3.0 6.0 WCA3S 102 4.6 80 28.0 33.7 ± 18.2 1.8 3.6 ENP 106 5.4 89 32.5 36.5 ± 20.5 2.0 3.9 LNWR Conductivity (umhos/cm) 57 37.1 1205 102.0 223.7 ± 281.9 37.3 74.8 WCA2 56 302 1423 706.0 764.1 ±291.5 39.0 78.1 WCA3N 69 165 801 456.0 490.7 ± 165.8 20.0 39.8 WCA3S 161 125 704 383.0 406.2 ± 136.5 10.8 21.3 ENP 175 72 655 372.0 378.9 ± 119.1 9.0 17.8 LNWR PH 57 5.66 7.39 6.4 6.5 ±0.5 0.1 0.1 WCA2 56 6.77 7.88 7.3 7.3 ±0.3 0.0 0.1 WCA3N 70 6.69 7.88 7.2 7.2 ±0.2 0.0 0.1 WCA3S 161 6.77 8.36 7.3 7.3 ±0.3 0.0 0.0 ENP 175 6.9 8.39 7.7 7.7 ±0.3 0.0 0.0 LNWR Sulfate (mg/L) 57 0.022 66 1.2 6.9 ± 15.1 2.0 4.0 WCA2 56 2.9 110 26.5 34.3 ± 26.3 3.5 7.1 WCA3N 70 0.12 52 5.7 13.1 ± 14.8 1.8 3.5 WCA3S 161 0.012 39 1.8 4.0 ±6.2 0.5 1.0 ENP 175 0.009 16 1.4 2.2 ±2.7 0.2 0.4 LNWR Total organic carbon (mg/L) 57 10.11 44.53 18.0 20.4 ± 7.9 1.0 2.1 WCA2 56 17.19 51.6 29.0 30.3 ±9.2 1.2 2.5 WCA3N 70 12 35.39 20.6 21.9 ± 6.0 0.7 1.4 WCA3S 161 7.9 29.63 16.7 16.7 ±4.3 0.3 0.7 ENP 175 4.6 31 12.0 13.2 ±5.1 0.4 0.8 113 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix IV. Soil characteristics by subarea for REMAP 1995, 1996, 1999, 2005 and 2014. Subarea Analyte Number of Samples Minimum Maximum Median Mean ± Standard Deviation Standard Error LNWR Thickness (feet) 105 1.30 14.00 8.57 8.7512.59 0.25 WCA 2 103 0.93 8.60 4.10 4.25 ± 1.60 0.16 WCA 3AN 129 0.20 5.20 1.60 1.74 ± 1.04 0.09 WCA 3AS 223 0.05 7.40 2.90 2.85 ± 1.33 0.09 WCA3B 78 1.00 6.53 3.25 3.51 ± 1.33 0.15 SRS 67 0.10 5.70 1.50 1.73 ± 1.12 0.14 NESS 44 0.10 4.30 1.64 1.67 ± 1.09 0.17 OM 79 0.00 4.77 0.83 1.1010.88 0.10 EMM 134 0.00 7.33 0.70 1.01 ± 1.11 0.10 TS 13 0.90 7.83 3.00 3.26 ± 2.07 0.58 LNWR Organic Matter (percent) 101 25.00 98.83 93.77 91.51 ± 10.51 1.05 WCA 2 103 20.58 100.00 87.00 82.92 ± 13.61 1.34 WCA 3AN 128 9.79 98.55 75.30 64.96 ± 24.57 2.17 WCA 3AS 224 1.19 98.70 87.85 80.52 ± 17.41 1.16 WCA3B 78 30.48 93.65 82.34 74.98 ± 16.88 1.91 SRS 67 12.00 91.24 78.00 66.58 ± 22.78 2.78 NESS 44 13.51 90.55 44.86 51.58 ± 25.87 3.90 OM 78 10.07 79.64 31.89 33.82 ± 15.66 1.77 EMM 133 1.57 85.14 20.57 24.80 ± 15.62 1.35 TS 13 10.54 93.73 21.84 39.24 ± 31.45 8.72 LNWR Bulk Density (g/cc) 79 0.04 0.16 0.06 0.0710.02 0.00 WCA 2 79 0.05 0.43 0.10 0.1110.06 0.01 WCA 3AN 99 0.08 0.90 0.17 0.2210.16 0.02 WCA 3AS 171 0.06 0.50 0.11 0.1210.06 0.00 WCA3B 60 0.05 0.32 0.14 0.1410.04 0.01 SRS 53 0.07 0.44 0.13 0.1510.08 0.01 NESS 35 0.08 0.88 0.21 0.2510.16 0.03 OM 59 0.11 0.73 0.30 0.3210.14 0.02 EMM 104 0.06 0.79 0.45 0.4410.16 0.02 TS 8 0.11 0.51 0.27 0.3010.13 0.05 114 ------- Everglades Regional Environmental Monitoring and Assessment Program (REMAP) Appendix IV (continued). Soil characteristics by subarea for REMAP 1995, 1996, 1999, 2005 and 2014. Number Mean ± 95% of Standard Standard Confidence Subarea Analyte Samples Minimum Maximum Median Deviation Error Interval LNWR 87 0.02 9.47 1.70 2.54 ±2.41 0.26 0.51 WCA 2 89 0.02 8.40 0.46 1.00 ± 1.53 0.16 0.32 WCA 3AN 126 0.02 7.90 1.17 1.67 ± 1.62 0.14 0.29 WCA 3AS Methyl Mercury (ug/kg) 214 0.02 6.54 0.57 0.97 ± 1.20 0.08 0.16 WCA3B 74 0.02 5.40 0.44 0.81 ±0.94 0.11 0.22 SRS 63 0.02 6.74 0.72 1.43 ± 1.57 0.20 0.40 NESS 44 0.02 4.69 0.29 0.62 ±0.89 0.13 0.27 OM 73 0.02 4.50 0.24 0.51 ±0.76 0.09 0.18 EMM 121 0.02 3.71 0.24 0.37 ±0.50 0.05 0.09 TS 9 0.11 2.00 0.31 0.51 ±0.59 0.20 0.45 LNWR 104 32.00 250.0 160.0 161.7 ±40.0 3.92 7.78 WCA 2 102 64.00 350.0 160.0 165.6 ±49.0 4.85 9.63 WCA 3AN 128 13.00 290.0 110.0 108.0 ±54.2 4.79 9.49 WCA 3AS Total Mercury (ug/kg) 222 28.00 290.0 160.0 159.9 ±52.4 3.52 6.93 WCA3B 77 51.00 290.0 140.0 138.9 ±46.4 5.28 10.52 SRS 67 42.00 290.0 140.2 150.3 ±59.2 7.23 14.44 NESS 44 47.00 210.0 125.0 120.3 ±50.0 7.54 15.22 OM 79 9.30 190.0 65.0 76.7 ±42.6 4.79 9.53 EMM 133 17.22 220.0 48.0 63.9 ±43.8 3.80 7.52 TS 13 22.00 160.0 46.0 65.3 ±45.5 12.62 27.50 LNWR 38 2.40 4.60 3.10 3.28 ±0.58 0.09 0.19 WCA 2 36 0.58 3.70 3.10 2.84 ±0.75 0.13 0.25 WCA 3AN 46 0.61 3.90 3.00 2.81 ±0.76 0.11 0.23 WCA 3AS Total 77 1.20 4.72 3.30 3.30 ±0.66 0.07 0.15 WCA3B Nitrogen 31 1.21 3.90 2.90 2.65 ±0.77 0.14 0.28 SRS (percent) 27 0.90 3.80 2.80 2.61 ±0.89 0.17 0.35 NESS 15 0.93 3.30 1.80 1.98 ±0.91 0.24 0.51 OM 25 0.64 3.10 1.40 1.49 ±0.68 0.14 0.28 EMM 43 0.21 1.80 0.78 0.86 ±0.35 0.05 0.11 LNWR 79 183.6 1400.0 327.3 387.4 ± 187.9 21.14 42.10 WCA 2 77 110.0 1591.6 410.0 484.2 ±274.3 31.25 62.25 WCA 3AN 101 97.2 1701.8 480.0 523.7 ±286.6 28.52 56.58 WCA 3AS Total Phosphorus (mg/kg) 172 115.7 880.0 390.4 413.9 ± 135.7 10.35 20.43 WCA3B 60 125.6 1700.0 353.5 399.7 ±239.0 30.86 61.75 SRS 53 170.7 771.8 381.9 401.8 ± 130.2 17.88 35.88 NESS 35 128.4 750.0 315.5 333.7 ± 140.0 23.66 48.08 OM 59 152.1 956.9 339.1 391.0 ± 196.9 25.64 51.32 EMM 104 55.7 1332.2 240.0 303.3 ±250.5 24.56 48.71 TS 8 83.7 610.0 192.4 220.2 ± 167.7 59.28 140.17 115 ------- |