Edward S. Rutherford

Joint analysis of stressors and ecosystem services to enhance restoration effectiveness

With increasing pressure placed on natural systems by growing human populations, both scientists and resource managers need a better understanding of the relationships between cumulative stress from human activities and valued ecosystem services. Societies often seek to mitigate threats to these services through large-scale, costly restoration projects, such as the over one billion dollar Great Lakes Restoration Initiative currently underway. To help inform these efforts, we merged high-resolution spatial analyses of environmental stressors with mapping of ecosystem services for all five Great Lakes. Cumulative ecosystem stress is highest in near-shore habitats, but also extends offshore in Lakes Erie, Ontario, and Michigan. Variation in cumulative stress is driven largely by spatial concordance among multiple stressors, indicating the importance of considering all stressors when planning restoration activities. In addition, highly stressed areas reflect numerous different combinations of stressors rather than a single suite of problems, suggesting that a detailed understanding of the stressors needing alleviation could improve restoration planning. We also find that many important areas for fisheries and recreation are subject to high stress, indicating that ecosystem degradation could be threatening key services. Current restoration efforts have targeted high-stress sites almost exclusively, but generally without knowledge of the full range of stressors affecting these locations or differences among sites in service provisioning. Our results demonstrate that joint spatial analysis of stressors and ecosystem services can provide a critical foundation for maximizing social and ecological benefits from restoration investments.

Using cultural ecosystem services to inform restoration priorities in the Laurentian Great Lakes.

Ecological restoration programs often attempt to maintain or enhance ecosystem services (ES), but fine-scale maps of multiple ES are rarely available to support prioritization among potential projects. Here we use agency reports, citizen science, and social media as data sources to quantify the spatial distribution of five recreational elements of cultural ES (CES) – sport fishing, recreational boating, birding, beach use, and park visitation – across North Americas Laurentian Great Lakes, where current restoration investments exceed US

Forecasting the Impacts of Silver and Bighead Carp on the Lake Erie Food Web

1.5 billion. These recreational CES are widely yet unevenly distributed, and spatial correlations among all except park visitation indicate that many locations support multiple CES benefits. Collectively, these five service metrics correlate with tourism gross domestic product, indicating that local economies benefit from ecosystem conditions that support CES. However, locations of high recreational CES delivery are often severely affected by environmental stressors, suggesting that either ecosystem condition or human enjoyment of these recreational CES is resilient even to substantial levels of stress. Our analyses show that spatial assessments of recreational CES are an informative complement to ecosystem stress assessments for guiding large-scale restoration efforts.

Fine-scale spatial variation in ice cover and surface temperature trends across the surface of the Laurentian Great Lakes

Nonindigenous bigheaded carps (Bighead Carp Hypophthalmichthys nobilis and Silver Carp H. molitrix; hereafter, “Asian carps” [AC]) threaten to invade and disrupt food webs and fisheries in the Laurentian Great Lakes through their high consumption of plankton. To quantify the potential effects of AC on the food web in Lake Erie, we developed an Ecopath with Ecosim (EwE) food web model and simulated four AC diet composition scenarios (high, low, and no detritus and low detritus with Walleye Sander vitreus and Yellow Perch Perca flavescens larvae) and two nutrient load scenarios (the 1999 baseline load and 2£ the baseline [HP]). We quantified the uncertainty of the potential AC effects by coupling the EwE model with estimates of parameter uncertainty in AC production, consumption, and predator diets obtained using structured expert judgment. Our model projected mean § SD AC equilibrium biomass ranging from 52 § 34 to 104 § 75 kg/ha under the different scenarios. Relative to baseline simulations without AC, AC invasion under all detrital diet scenarios decreased the biomass of most fish and zooplankton groups. The effects of AC in the HP scenario were similar to those in the detrital diet scenarios except that the biomasses of most Walleye and Yellow Perch groups were greater under HP because these fishes were buffered from competition with AC by increased productivity at lower trophic levels. Asian carp predation on Walleye and Yellow Perch larvae caused biomass declines among all Walleye and Yellow Perch groups. Large food web impacts of AC occurred in only 2% of the simulations, where AC biomass exceeded 200 kg/ha, resulting in biomass declines of zooplankton and planktivorous fish near the levels observed in the Illinois River. Our findings suggest that AC would affect Lake Erie’s food web by competing with other planktivorous fishes and by providing additional prey for piscivores. Our methods provide a novel approach for including uncertainty into forecasts of invasive species’ impacts on aquatic food webs.

Annual variation in habitat-specific recruitment success: implications from an individual-based model of Lake Michigan alewife (Alosa pseudoharengus)

The effects of climate change on north temperate freshwater ecosystems include increasing water temperatures and decreasing ice cover. Here we compare those trends in the Laurentian Great Lakes at three spatial scales to evaluate how warming varies across the surface of these massive inland water bodies. We compiled seasonal ice cover duration (1973–2013) and lake summer surface water temperatures (LSSWT; 1994–2013), and analyzed spatial patterns and trends at lake-wide, lake sub-basin, and fine spatial scales and compared those to reported lake- and basin-wide trends. At the lake-wide scale we found declining ice duration and warming LSSWT patterns consistent with previous studies. At the lake sub-basin scale, our statistical models identified distinct warming trends within each lake that included significant breakpoints in ice duration for 13 sub-basins, consistent linear declines in 11 sub-basins, and no trends in 4 sub-basins. At the finest scale, we found that the northern- and eastern-most portions of each Great Lake, especially in nearshore areas, have experienced faster rates of LSSWT warming and shortening ice duration than those previously reported from trends at the lake scale. We conclude that lake-level analyses mask significant spatial and temporal variation in warming patterns within the Laurentian Great Lakes. Recognizing spatial variability in rates of change can inform both mechanistic modeling of ecosystem responses and planning for long-term management of these large freshwater ecosystems.

Investigation of interbasin exchange and interannual variability in Lake Erie using an unstructured‐grid hydrodynamic model

The identification of important spawning and nursery habitats for fish stocks can aid fisheries management, but is complicated by various factors, including annual variation in recruitment success. The alewife (Alosa pseudoharengus) is an ecologically important species in Lake Michigan that utilizes a variety of habitats for spawning and early life growth. While productive, warm tributary mouths (connected to Lake Michigan) may contribute disproportionately more recruits (relative to their habitat volume) to the adult alewife population than cooler, less productive nearshore habitats, the extent of interannual variation in the relative contributions of recruits from these two habitat types remains unknown. We used an individual-based bioenergetics simulation model and input data on daily temperatures to estimate alewife recruitment to the adult population by these different habitat types. Simulations suggest that nearshore lake habitats typically produce the vast majority of young alewife recruits. Howeve...

Rating impacts in a multi‐stressor world: a quantitative assessment of 50 stressors affecting the Great Lakes

Interbasin exchange and interannual variability in Lake Eries three basins are investigated with the help of a three-dimensional unstructured-grid-based Finite Volume Coastal Ocean Model (FVCOM). Experiments were carried out to investigate the influence of grid resolutions and different sources of wind forcing on the lake dynamics. Based on the calibrated model, we investigated the sensitivity of lake dynamics to major external forcing, and seasonal climatological circulation patterns are presented and compared with the observational data and existing model results. It was found that water exchange between the western basin (WB) and the central basin (CB) was mainly driven by hydraulic and density-driven flows, while density-driven flows dominate the interaction between the CB and the eastern basin (EB). River-induced hydraulic flows magnify the eastward water exchange and impede the westward one. Surface wind forcing shifts the pathway of hydraulic flows in the WB, determines the gyre pattern in the CB, contributes to thermal mixing, and magnifies interbasin water exchange during winter. Interannual variability is mainly driven by the differences in atmospheric forcing, and is most prominent in the CB.

Use of structured expert judgment to forecast invasions by bighead and silver carp in Lake Erie

Ecosystems often experience multiple environmental stressors simultaneously that can differ widely in their pathways and strengths of impact. Differences in the relative impact of environmental stressors can guide restoration and management prioritization, but few studies have empirically assessed a comprehensive suite of stressors acting on a given ecosystem. To fill this gap in the Laurentian Great Lakes, where considerable restoration investments are currently underway, we used expert elicitation via a detailed online survey to develop ratings of the relative impacts of 50 potential stressors. Highlighting the multiplicity of stressors in this system, experts assessed all 50 stressors as having some impact on ecosystem condition, but ratings differed greatly among stressors. Individual stressors related to invasive and nuisance species (e.g., dreissenid mussels and ballast invasion risk) and climate change were assessed as having the greatest potential impacts. These results mark a shift away from the longstanding emphasis on nonpoint phosphorus and persistent bioaccumulative toxic substances in the Great Lakes. Differences in impact ratings among lakes and ecosystem zones were weak, and experts exhibited surprisingly high levels of agreement on the relative impacts of most stressors. Our results provide a basin-wide, quantitative summary of expert opinion on the present-day influence of all major Great Lakes stressors. The resulting ratings can facilitate prioritizing stressors to achieve management objectives in a given location, as well as providing a baseline for future stressor impact assessments in the Great Lakes and elsewhere.

Modeling the Transport of Larval Yellow Perch in Lake Michigan

Identifying which nonindigenous species will become invasive and forecasting the damage they will cause is difficult and presents a significant problem for natural resource management. Often, the data or resources necessary for ecological risk assessment are incomplete or absent, leaving environmental decision makers ill equipped to effectively manage valuable natural resources. Structured expert judgment (SEJ) is a mathematical and performance-based method of eliciting, weighting, and aggregating expert judgments. In contrast to other methods of eliciting and aggregating expert judgments (where, for example, equal weights may be assigned to experts), SEJ weights each expert on the basis of his or her statistical accuracy and informativeness through performance measurement on a set of calibration variables. We used SEJ to forecast impacts of nonindigenous Asian carp (Hypophthalmichthys spp.) in Lake Erie, where it is believed not to be established. Experts quantified Asian carp biomass, production, and consumption and their impact on 4 fish species if Asian carp were to become established. According to experts, in Lake Erie Asian carp have the potential to achieve biomass levels that are similar to the sum of biomasses for several fishes that are harvested commercially or recreationally. However, the impact of Asian carp on the biomass of these fishes was estimated by experts to be small, relative to long term average biomasses, with little uncertainty. Impacts of Asian carp in tributaries and on recreational activities, water quality, or other species were not addressed. SEJ can be used to quantify key uncertainties of invasion biology and also provide a decision-support tool when the necessary information for natural resource management and policy is not available.

Assessment of Top-Down and Bottom-Up Controls on the Collapse of Alewives (Alosa pseudoharengus) in Lake Huron

The transport of larval yellow perch (Perca flavescens) in Lake Michigan is studied with a 3D particle trajectory model. The model uses 3D currents generated by the Great Lakes version of the Princeton Ocean Model driven by observed momentum and heat fluxes in June-August 1998, 1999 and 2000. Virtual larvae were released in the nearshore region with the most abundant preferred substrate for yellow perch spawning, rocks. This paper also investigated the potential for physical transport mechanisms to affect recruitment of Lake Michigan yellow perch by coupling hydrodynamic models with individual-based particle models of fish larvae to study variation in larval distributions, growth rates, and potential recruitment. Larval growth rates were simulated using a bioenergetics growth model with fixed consumption rates. Results indicate that lake circulation patterns are critical for understanding interannual variability in Great Lake fish recruitment.