GLERL Publication Abstracts: FY 2014

Capitalized names represent GLERL authors.

 


Alderstein, S.A., T.F. NALEPA, H.A. VANDERPLOEG, and G.L. FAHNENSTIEL. Trends in phytoplankton, zooplankton, and macroinvertebrates in Saginaw Bay relative to zebra mussel (Dreissena polymorpha) colonization: A generalized linear model approach. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 525-543 pp. (2013).

We quantify temporal and spatial trends in densities of main taxonomic groups of phytoplankton, zooplankton, and benthic macroinvertebrates relative to the colonization of Saginaw Bay, Lake Huron, by zebra mussels (Dreissena polymorpha) in 1990-1996. We used data from bay-wide surveys and generalized linear models (GLMs). Mussels were first found in 1991, peaked in 1992, and declined to stable levels in 1993 at stations with hard substrates where they were most abundant. Annual trends in phytoplankton and zooplankton were negatively correlated with zebra mussel trends at different time lags. Most phytoplankton taxa declined within three years of mussel colonization, with disappearance of photosensitive cyanophytes in the first year of colonization. All phytoplankton taxa tended to recover within 5 years except for cyanophytes and chlorophytes. Diatom annual densities were least affected, but species composition changed. Major zooplankton groups declined after 1 year of mussel colonization. Cyclopoids and cladocerans exhibited lowest densities in 1993, and calanoids and rotifers continued to decline through 1996. Spatial distributions of both plankton groups were fairly homogeneous and remained stable despite changes in densities and patchy mussel distributions. Responses of macroinvertebrate taxa depended on their life history, mobility, and spatial proximity to mussel colonies. Most taxa except for Gammarus and Diporeia were most abundant in silt and silty-sand substrate in the inner bay where mussel densities were lowest. Oligochaetes declined between 1991 and 1993 but then increased after 1994, chironomids showed no strong patterns but tended to increase after 1992, sphaeriids decreased in the inner bay after 1991, Gammarus increased between 1991 and 1993, and Diporeia declined after 1991. Our results highlight the importance of prompt analysis of monitoring data collected during the initial-colonization period. Such data can then provide information to optimize survey design during ongoing assessments of zebra mussel impacts.

ANDERSON, E.J., and D.J. SCHWAB. Predicting the oscillating bi-directional exchange flow in the Straits of Mackinac. Journal of Great Lakes Research 39(4):663-671 (DOI:10.1016/j.jglr.2013.09.001) (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130038.pdf

The Straits of Mackinac are a unique feature that connects Lake Michigan and Lake Huron into a single hydraulically linked system. With currents of up to 1 m/s and oscillating volumetric transport up to 80,000 m3/s, they play an important role in water quality, contaminant transport, navigation, and ecological processes. We present the first three-dimensional hydrodynamic model of the combined Lake Michigan–Huron, including the Straits of Mackinac at high-resolution, that is able to simulate the three dimensional structure of the oscillating flows at the Straits. In comparison with individual lake models for Michigan and Huron (no connection at the Straits), we are able to isolate the effects of the bi-lake oscillation and have found that although the oscillation (Helmholtz mode) is the dominant forcing mechanism, the flow can be modulated when atmospheric systems are in-phase with water level fluctuations. Furthermore, the area of influence of the Straits is found to extend up to 70 km into each lake, underscoring the need for realistic predictions within the Straits. For the first time, this combined lake hydrodynamic model provides the capability to investigate and accurately predict flow at the Straits of Mackinac and its effect on Lake Michigan and Huron. This model forms the basis for the next generation of real-time hydrodynamic models being developed for the Great Lakes Coastal Forecasting System, a suite of models designed by the National Oceanic and Atmospheric Administration Great Lakes Environmental Research Laboratory (NOAA/GLERL) that predict hydrodynamic conditions such as currents, temperatures, and water levels in three dimensions.

ANDERSON, E.J., and D.J. SCHWAB. Spill reference tables for the St. Clair River. NOAA Technical Memorandum GLERL-162. NOAA, Great Lakes Environmental Research Laboratory, Ann Arbor, MI, 16 pp. (2014). https://www.glerl.noaa.gov/pubs/tech_reports/glerl-162/tm-162.pdf

Serving as the international border between the United States and Canada, the St. Clair River provides a source for public drinking water as well as a hub for petrochemical refineries and commercial shipping (Fig. 1). This dichotomy creates the potential threat to public health via drinking water contamination in the event of a toxic spill from commercial freighters or refineries within the waterway. Due to the high flow rates experienced in the river (average discharge of 5,200 m3/s; Holtschlag and Koschik 2002), once released, contaminants can travel to downstream public water intakes within minutes (Anderson and Schwab, 2012; Tsanis et al., 1996; Derecki 1983; Sun et al., 2013). As a result, water intake managers may need to react immediately in order to mitigate contaminant uptake. However, the common approach for spill response is to use real-time operational hydrodynamic models to predict currents in the system and then simulate spill transport forecasts via Lagrangian particle transport simulations, a process that may not provide decision makers with the necessary spill information in time.

BAI, X., J. WANG, J. Austin, D.J. SCHWAB, R.A. ASSEL, A.H. CLITES, J.F. BRATTON, M.C. COLTON, J. Lenters, B.M. LOFGREN, T. Wohlleben, S. Helfrich, H.A. VANDERPLOEG, L. LUO, and G.A. LESHKEVICH. A record-breaking low ice cover over the Great Lakes during winter 2011/2012: Combined effects of a strong positive NAO and La Nina. Climate Dynamics:27 pp. (DOI:10.1007/s00382-014-2225-2) (2014).

A record-breaking low ice cover occurred in the North American Great Lakes during winter 2011/2012, in conjunction with a strong positive Arctic Oscillation/ North Atlantic Oscillation (AO/NAO) and a La Nina event. Large-scale atmosphere circulation in the Pacific/ North America (PNA) region reflected a combined signal of La Nina and NAO. Surface heat flux analysis shows that sensible heat flux contributed most to the net surface heat flux anomaly. Surface air temperature is the dominant factor governing the interannual variability of Great Lakes ice cover. Neither La Nina nor NAO alone can be responsible for the extreme warmth; the typical midlatitude response to La Nina events is a negative PNA pattern, which does not have a significant impact on Great Lakes winter climate; the positive phase of NAO is usually associated with moderate warming. When the two occurred simultaneously, the combined effects of La Nina and NAO resulted in a negative East Pacific pattern with a negative center over Alaska/Western Canada, a positive center in the eastern North Pacific (north of Hawaii), and an enhanced positive center over the eastern and southern United States. The overall pattern prohibited the movement of the Arctic air mass into mid-latitudes and enhanced southerly flow and warm advection from the Gulf of Mexico over the eastern United States and Great Lakes region, leading to the record-breaking low ice cover. It is another climatic pattern that can induce extreme warming in the Great Lakes region in addition to strong El Nino events. A very similar event occurred in the winter of 1999/2000. This extreme warm winter and spring in 2012 had significant impacts on the physical environment, as well as counterintuitive effects on phytoplankton abundance.

BELETSKY, D., N. HAWLEY, and Y. Rao. Modeling summer circulation and thermal structure of Lake Erie. Journal of Geophysical Research 118:1-15 (DOI:10.1002/2013JC008854) (2013).

A three-dimensional primitive equation numerical model was applied to Lake Erie on a 2 km grid to study its summer circulation and thermal structure. Model results were compared to long-term observations of currents and temperature made in 2005 at several locations, mostly in its central basin. In the shallow and mostly unstratified western basin circulation is driven by Detroit River inflow (modified to some extent by wind) and is from west to east. In the central basin (which is of intermediate depth and has a relatively flat bottom), the modeled circulation is anticyclonic (clockwise), driven by anticyclonic vorticity in the surface wind, and the thermocline is bowl-shaped, in line with observations. In the deep part of the eastern basin, the thermocline is dome-shaped and circulation is cyclonic (counterclockwise), due to density gradients (a configuration typical for other large deep lakes), while shallower areas are occupied by anticyclonic circulation driven by anticyclonic wind vorticity. In the central basin, modeled temperature and circulation patterns are quite sensitive to the specification of the wind field. Anticyclonic wind vorticity leads to thinning of the hypolimnion in the central basin and earlier destratification in the fall.

BELETSKY, D., H. HU, J. WANG, and N. HAWLEY. Modeling thermal structure, circulation, and ice in Lake Erie. Proceedings, 22nd IAHR International Symposium on Ice 2014, Singapore, August 11-15, 2014. International Association of Hydro-Environment Engineering and Research, 449-451 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140053.pdf

Winter circulation and thermal structure in Lake Erie is studied with a three-dimensional coupled Great Lakes Ice-circulation Model (GLIM) during 1979-1980 and 2010-2011. The hydrodynamic model has 20 vertical levels and a uniform horizontal grid size of 2 km. The model uses time-dependent wind stress and heat flux forcing at the surface which are calculated from the hourly meteorological observations obtained from National Weather Service land stations and NOAA buoys. The model reproduces several month long ice periods in both winters with maximum ice thicknesses up to 45 cm in 1979-1980. The model is validated with hourly temperature and current measurements at several moorings that were deployed in central basin of Lake Erie during both winters and additionally with ice concentration and thickness measurements in 2010-2011.

BELETSKY, R., D. BELETSKY, N. HAWLEY, and J. WANG. Interannual variability of winter circulation and ice in Lake Erie. Proceedings, 22nd IAHR International Symposium on Ice 2014, Singapore, August 11-15, 2014. International Association of Hydro-Environment Engineering and Research, 905-907 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140055.pdf

Winter circulation in the Laurentian Great Lakes is less known than summer circulation and impacts of ice on lake circulation are poorly understood. Lake Erie is at least partially covered with ice from December until April, and its normal peak ice cover is about 90%. To study thermal structure, circulation, and ice thickness, 7 current profilers, 4 ice profilers and 37 temperature sensors were deployed in Lake Erie in October 2010 - May 2011, and again in September 2011 - May 2012. Wind speed observations were obtained from National Weather Service and Environment Canada meteorological stations. Significant interannual variability in the extent of ice cover was observed during the two field years, and that allowed us to study the impact of ice on thermal structure and circulation in Lake Erie.

CHA, Y.K., S.S. Park, K. Kim, M. Byeon, and C.A. STOW. Probabilistic prediction of cyanobacteria abundance in a Korean reservoir using a Bayesian Poisson model. Water Resources Research 50:2518-2532 (DOI:10.1002/2013WR014372) (2014).

There have been increasing reports of harmful algal blooms (HABs) worldwide. However, the factors that influence cyanobacteria dominance and HAB formation can be site-specific and idiosyncratic, making prediction challenging. The drivers of cyanobacteria blooms in Lake Paldang, South Korea, the summer climate of which is strongly affected by the East Asian monsoon, may differ from those in well studied North American lakes. Using the observational data sampled during the growing season in 2007– 2011, a Bayesian hurdle Poisson model was developed to predict cyanobacteria abundance in the lake. The model allowed cyanobacteria absence (zero count) and nonzero cyanobacteria counts to be modeled as functions of different environmental factors. The model predictions demonstrated that the principal factor that determines the success of cyanobacteria was temperature. Combined with high temperature, increased residence time indicated by low outflow rates appeared to increase the probability of cyanobacteria occurrence. A stable water column, represented by low suspended solids, and high temperature were the requirements for high abundance of cyanobacteria. Our model results had management implications; the model can be used to forecast cyanobacteria watch or alert levels probabilistically and develop mitigation strategies of cyanobacteria blooms.

CHA, Y.K., and C.A. STOW. A Bayesian network incorporating observation error to predict phosphorus and chlorophyll a in Saginaw Bay. Environmental Modelling and Software:11 pp. (DOI:10.1016/j.envsoft.2014.02.010) (2014).

Empirical relationships between lake chlorophyll a and total phosphorus concentrations are widely used to develop predictive models. These models are often estimated using sample averages as implicit surrogates for unknown lake-wide means, a practice than can result in biased parameter estimation and inaccurate predictive uncertainty. We develop a Bayesian network model based on empirical chlorophyll-phosphorus relationships for Saginaw Bay, an embayment on Lake Huron. The model treats the means as unknown parameters, and includes structure to accommodate the observation error associated with estimating those means. Compared with results from an analogous simple model using sample averages, the observation error model has a lower predictive uncertainty and predicts lower chlorophyll and phosphorus concentrations under contemporary lake conditions. These models will be useful to guide pending decision-making pursuant to the 2012 Great Lakes Water Quality Agreement.

CLITES, A.H., J.P. SMITH, T.S. HUNTER, and A.D. GRONEWOLD. Visualizing relationships between hydrology, climate, and water level fluctuations on Earth's largest system of lakes. Journal of Great Lakes Research 40:807-811 (DOI:10.1016/j.jglr.2014.05.014) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140026.pdf

Understanding drivers behind monthly, annual, and decadal water level fluctuations on the North American Great Lakes is a high priority for regional research and water resource management planning. The need for improved understanding of these relationships is underscored by a series of recent unprecedented extreme water level patterns, including (but not limited to) record low water levels on Lakes Michigan and Huron in December 2012 and January 2013. To address this need, we developed the Great Lakes Hydro-Climate Dashboard (GLHCD), a dynamic flash-based web interface that builds upon the previously-released Great Lakes Water Level Dashboard (GLWLD). In addition to including water level data and projections from the GLWLD, the GLHCD presents a range of hydrological and climatological data through an improved graphical user interface specifically designed to manage, and display simultaneously, a variety of data time series from different sources. By serving as a common portal to critical regional hydro-climate and water level data, the GLHCD helps visualize and explain lake level phenomena including water level declines across all of the Great Lakes in the early 1960s and their relationship to changes in regional precipitation, as well as the abrupt water level declines in the late 1990s and their relationship to remarkable changes in over-lake evaporation. By providing insight into these, and other important regional hydro-climate events, the GLHCD helps practitioners, researchers, and the general public improve their understanding of the drivers behind Great Lakes water levels, and to employ that understanding in prudent water resource management planning.

CLITES, A.H., J. WANG, K.B. CAMPBELL, A.D. GRONEWOLD, R.A. ASSEL, X. BAI, and G.A. LESHKEVICH. Cold water and high ice cover on Great Lakes in spring 2014. EOS 95(34):305-306 (DOI:10.1002/2014EO340001) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140042.pdf

Very cold temperatures across much of North America caused by the recent anomalous meridional upper air flow—commonly referred to in the public media as a polar vortex (for details, see Blackmon et al. [1977] and National Climatic Data Center, State of the climate: Synoptic discussion for January 2014, https:// www.ncdc.noaa.gov/sotc/synoptic/2014/1)—have contributed to extreme hydrologic conditions on the Great Lakes. The Great Lakes are the largest system of lakes and the largest surface of freshwater on Earth—Lake Superior alone is the single largest lake by surface area.

Cooke, R.M., M.E. Wittmann, D.M. Lodge, J.D. Rothlisberger, E.S. RUTHERFORD, H. ZHANG, and D.M. MASON. Out-of-sample validation for structured expert judgment of Asian Carp establishment in Lake Erie. Integrated Environmental Assessment and Management 10(4):522-528 (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140047.pdf

Structured expert judgment (SEJ) is used to quantify the uncertainty of nonindigenous fish (bighead carp [Hypophthalmichthys nobilis] and silver carp [H. molitrix]) establishment in Lake Erie. The classical model for structured expert judgment model is applied. Forming a weighted combination (called a decision maker) of experts' distributions, with weights derived from performance on a set of calibration variables from the experts' field, exhibits greater statistical accuracy and greater informativeness than simple averaging with equal weights. New methods of cross validation are applied and suggest that performance characteristics relative to equal weighting could be predicted with a small number (1–2) of calibration variables. The performance-based decision maker is somewhat degraded on out-of-sample prediction, but remained superior to the equal weight decision maker in terms of statistical accuracy and informativeness.

Cosens, B., and C.A. STOW. Chapter 5: Resilience and water governance: Addressing fragmentation and uncertainty in water allocation and water quality law. In Social-Ecological Resilience and Law. A.S. Garmestani and C.R. Allen (Eds.). Columbia University Press, New York, NY, 142-175 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140003.pdf

The U.S. EPA reports that almost half of the nation's rivers and two thirds of its lakes are use-impaired due to poor water quality (U.S. EPA 1998, 2002, 2010; Houck 2002). The Western Water Policy Review Advisory Commission identified both poor water quality and unhealthy aquatic systems among the water challenges facing the West (Western Water Policy Review Advisory Commission 1998). The water quality impairment is caused both by chemical pollution and physical alteration of streams. Nutrients and excess sediment impair water quality in 30 percent of the nation's streams (U.S. EPA 2011). In the Great Basin nearly two-thirds of the native fish are either listed under the Endangered Species Act (ESA) or considered of concern by the U.S. Fish and Wildlife Service (USFWS). Water development is considered second only to the introduction of nonnative fish in causing these problems (Doremus 2001).

Cross, V.A., J.F. BRATTON, H.A. Michael, K.D. Kroeger, A. Green, and E. Bergeron. Continuous resistivity profiling and seismic-reflection data collected in April 2010 from Indian River Bay, Delaware. Open File Report 2011-1039. U.S. Geological Survey, Reston, VA, 32 pp. (2014). https://pubs.usgs.gov/of/2011/1039/ofr2011-1039-title_page.html

A geophysical survey to delineate the fresh-saline groundwater interface and associated sub-bottom sedimentary structures beneath Indian River Bay, Delaware, was carried out in April 2010. This included surveying at higher spatial resolution in the vicinity of a study site at Holts Landing, where intensive onshore and offshore studies were subsequently completed. The total length of continuous resistivity profiling (CRP) survey lines was 145 kilometers (km), with 36 km of chirp seismic lines surveyed around the perimeter of the bay. Medium-resolution CRP surveying was performed using a 50-meter streamer in a baywide grid. Results of the surveying and data inversion showed the presence of many buried paleochannels beneath Indian River Bay that generally extended perpendicular from the shoreline in areas of modern tributaries, tidal creeks, and marshes. An especially wide and deep paleochannel system was imaged in the southeastern part of the bay near White Creek. Many paleochannels also had high-resistivity anomalies corresponding to low-salinity groundwater plumes associated with them, likely due to the presence of fine-grained estuarine mud and peats in the channel fills that act as submarine confining units. Where present, these units allow plumes of low-salinity groundwater that was recharged onshore to move beyond the shoreline, creating a complex fresh-saline groundwater interface in the subsurface. The properties of this interface are important considerations in construction of accurate coastal groundwater flow models. These models are required to help predict how nutrient-rich groundwater, recharged in agricultural watersheds such as this one, makes its way into coastal bays and impacts surface-water quality and estuarine ecosystems.

DAVIS, T.W., S.B. Watson, M.J. Rozmarynowycz, J.J.H. Ciborowski, R.M. McKay, and G.S. Bullerjahn. Phylogenies of microcystin-producing cyanobacteria in the lower Laurentian Great Lakes suggest extensive genetic connectivity. PLOS One 9(9):9 pp. (DOI:10.1371/journal.pone.0106093) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140048.pdf

Lake St. Clair is the smallest lake in the Laurentian Great Lakes system. MODIS satellite imagery suggests that high algal biomass events have occurred annually along the southern shore during late summer. In this study, we evaluated these events and tested the hypothesis that summer bloom material derived from Lake St. Clair may enter Lake Erie via the Detroit River and represent an overlooked source of potentially toxic Microcystis biomass to the western basin of Lake Erie. We conducted a seasonally and spatially resolved study carried out in the summer of 2013. Our goals were to: 1) track the development of the 2013 summer south-east shore bloom 2) conduct a spatial survey to characterize the extent of toxicity, taxonomic diversity of the total phytoplankton population and the phylogenetic diversity of potential MC-producing cyanobacteria (Microcystis, Planktothrix and Anabaena) during a high biomass event, and 3) compare the strains of potential MC-producers in Lake St. Clair with strains from Lake Erie and Lake Ontario. Our results demonstrated a clear predominance of cyanobacteria during a late August bloom event, primarily dominated by Microcystis, which we traced along the Lake St. Clair coastline downstream to the Detroit River’s outflow at Lake Erie. Microcystin levels exceeded the Province of Ontario Drinking Water Quality Standard (1.5 mg L21) for safe drinking water at most sites, reaching up to five times this level in some areas. Microcystis was the predominant microcystin producer, and all toxic Microcystis strains found in Lake St. Clair were genetically similar to toxic Microcystis strains found in lakes Erie and Ontario. These findings suggest extensive genetic connectivity among the three systems.

Deal, C.J., N. Steiner, J. Christian, J. Clement Kinney, K. Denman, S. Elliott, G. Gibson, M. Jin, D. Lavoie, S. Lee, W. Lee, W. Maslowski, J. WANG, and E. Watanabe. Chapter 12. Progress and challenges in biogeochemical modeling of the Pacific Arctic region. In The Pacific Arctic Region: Ecosystem Status and Trends in a Rapidly Changing Environment. J.M. Grebmeier and W. Maslowski (Eds.). Springer Science+Business Media, Dordrecht, 393-445 pp. (DOI:10.1007/978-94-017-8863-2_12) (2014).

At this early stage of modeling marine ecosystems and biogeochemical cycles in the Pacific Arctic Region (PAR), numerous challenges lie ahead. Observational data used for model development and validation remain sparse, especially across seasons and under a variety of environmental conditions. Field data are becoming more available, but at the same time PAR is rapidly changing. Biogeochemical models can provide the means to capture some of these changes. This study introduces and synthesizes ecosystem modelling in PAR by discussing differences in complexity and application of one-dimensional, regional, and global earth system models. Topics include the general structure of ecosystem models and specifics of the combined benthic, pelagic, and ice PAR ecosystems, the importance of model validation, model responses to climate influences (e.g. diminishing sea ice, ocean acidification), and the impacts of circulation and stratification changes on PAR ecosystems and biogeochemical cycling. Examples of modeling studies that help place the region within the context of the Pan-Arctic System are also discussed. We synthesize past and ongoing PAR biogeochemical modelling efforts and briefly touch on decision makers’ use of ecosystem models and on necessary future developments.

Dufour, M.R., J.J. Pritt, C.M. Mayer, C.A. STOW, and S.S. QIAN. Bayesian hierarchical modeling of larval walleye (Sander vitreus) abundance and mortality: Accounting for spatial and temporal variability on a large river. Ecological Applications:12 pp. (DOI:10.1016/j.jglr.2014.08.001) (2014).

Larval fish are extremely variable in space and time while sampling of populations is generally restricted and incomplete. However, estimates of abundance and mortality are important for understanding population dynamics, habitat quality, and anthropogenic impacts. Acknowledging and addressing variability during sampling and data analysis are imperative to producing informative estimates. A combination of spatially and temporally distributed ichthyoplankton sampling and Bayesian hierarchical and state-space modeling was used to partition variance and estimate abundance and mortality of larval walleye (Sander vitreus) in the Maumee River during 2010 and 2011. System variability and degree of sampling coverage have a direct impact on the quality of abundance estimates. Small scale factors (i.e., within site and day-to-day) accounted for the most variation in larval walleye densities, therefore sampling should concentrate on capturing these sources. Bayesian state-space modeling can improve estimates by sharing information through time, properly accounting for uncertainty, and producing probability distribution based estimates. Larval fish are highly variable and difficult to sample; however, the application of Bayesian methods during the data analysis process can lead to improved estimates of abundance and informed management actions.

Filstrup, C.T., T. Wagner, P.A. Soranno, E.H. Stanley, C.A. STOW, K.E. Webster, and E.J. Downing. Regional variability among nonlinear chlorophyll-phosphorus relationships in lakes. Limnology and Oceanography 59(5):1691-1703 (DOI:10.4319/lo.2014.59.5.1691) (2014).

The relationship between chlorophyll a (Chl a) and total phosphorus (TP) is a fundamental relationship in lakes that reflects multiple aspects of ecosystem function and is also used in the regulation and management of inland waters. The exact form of this relationship has substantial implications on its meaning and its use. We assembled a spatially extensive data set to examine whether nonlinear models are a better fit for Chl a–TP relationships than traditional log-linear models, whether there were regional differences in the form of the relationships, and, if so, which regional factors were related to these differences. We analyzed a data set from 2105 temperate lakes across 35 ecoregions by fitting and comparing two different nonlinear models and one log-linear model. The two nonlinear models fit the data better than the log-linear model. In addition, the parameters for the best-fitting model varied among regions: the maximum and lower Chl a asymptotes were positively and negatively related to percent regional pasture land use, respectively, and the rate at which chlorophyll increased with TP was negatively related to percent regional wetland cover. Lakes in regions with more pasture fields had higher maximum chlorophyll concentrations at high TP concentrations but lower minimum chlorophyll concentrations at low TP concentrations. Lakes in regions with less wetland cover showed a steeper Chl a–TP relationship than wetland-rich regions. Interpretation of Chl a–TP relationships depends on regional differences, and theory and management based on a monolithic relationship may be inaccurate.

Foley, C.J., G.J. Bowen, T.F. NALEPA, M.S. Sepulveda, and T.O. Hook. Stable isotope patterns of benthic organisms from the Great Lakes region indicate variable dietary overlap of Diporeia spp. and dreissenid mussels. Canadian Journal of Fisheries and Aquatic Sciences 71:12 pp. (DOI:10.1139/cjfas-2013-0620) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140050.pdf

Competition between native and invasive species may bring about a suite of ecological and evolutionary outcomes, including local extirpations. In the Laurentian Great Lakes, competition for food may explain the dramatic decline of Diporeia spp. amphipods following the introduction of dreissenid mussels. This hypothesis has not been confirmed, in part because dreissenids and Diporeia appear to co-exist and flourish in other systems, including the Finger Lakes of New York. We used carbon, nitrogen, hydrogen, and oxygen stable isotope ratios to examine resource use by Diporeia from three spatially distinct populations (Lake Michigan, Lake Superior, and Cayuga Lake), dreissenids from areas where they co-occur with Diporeia (Lake Michigan and Cayuga Lake), and Diporeia from Lake Michigan collected before and after dreissenid invasion (1986–2009). Our results suggest that dreissenids may affect resource use by Diporeia in areas of co-occurrence, but the extent to which those effects are positive or negative is unclear. Terrestrial inputs may provide an important subsidy for Diporeia populations in small systems but may not be substantial enough in the Great Lakes to ensure that both taxa thrive.

Francoeur, S.N., K.A. Peters-Winslow, D. Miller, C.A. STOW, Y.K. CHA, and S.D. Peacor. Spatial and temporal patterns in macroscopic benthic primary producers in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 40(Supplement 1):53-63 (DOI:10.1016/j.jglr.2014.01.011) (2014).

We investigated spatial and temporal patterns in macroscopic benthic primary producer biomass, production, and composition in inner Saginaw Bay in 2009 and 2010. Charophytes and filamentous algae (FA) were relatively abundant, and vascular macrophytes were less common. The probability of benthic primary producer presence increased with the proportion of benthic substrate composed of rock. Most benthic primary producer biomass occurred at depths of 2–4 m, with very little biomass observed beyond 4 m deep. Charophyte and vascular macrophyte abundances displayed consistent patterns related to distance from the mouth of the Saginaw River. FA abundance also displayed such patterns, but they reversed between 2009 and 2010. Macrophytic benthic primary producer communities were generally dominated by charophytes. Three genera of vascular macrophytes, including Myriophyllum, were also observed. Filamentous algal communities were composed of a mixture of FA taxa. Ten FA genera were observed, including the red alga Compsopogon. Dominance of Compsopogon was related to low water clarity and low TP. Biomass-based benthic production estimates indicated that charophytes and FA strongly dominated macroscopic benthic production; production of vascular macrophytes was relatively low. The observed relationships of abundance and environmental conditions suggested regulation of benthic producer biomass by a shifting mosaic of substratum, nutrient, and light availabilities. The diverse nature of the benthic producer community could complicate understanding and management of excess benthic biomass and beach fouling in Saginaw Bay.

FRY, L.M., A.D. GRONEWOLD, V. Fortin, S. Buan, A.H. CLITES, C. Luukkonen, D. Hotsschlag, D. L., T.S. HUNTER, F. Seglenieks, D. Durnford, M. Dimitrijevic, C. Subich, E. Klyszejko, K. KEA, and P. Restrepo. The Great Lakes Runoff Intercomparison Project, Phase 1: Lake Michigan (GRIP-M). Journal of Hydrology:18 pp. (DOI:10.1016/j.jhydrol.2014.07.021) (2014).

We assembled and applied five models (one of which included three different configurations) to the Lake Michigan basin to improve our understanding of how differences in model skill at simulating total runoff to Lake Michigan relate to model structure, calibration protocol, model complexity, and assimilation (i.e. replacement of simulated discharge with discharge observations into historical simulations), and evaluate historical changes in runoff to Lake Michigan. We found that the performance among these models when simulating total runoff to the lake varied relatively little, despite variability in model structure, spatial representation, input data, and calibration protocol. Relatively simple empirical, assimilative models, including the National Oceanic and Atmospheric Administration (NOAA) Great Lakes Environmental Research Laboratory (GLERL) area ratio-based model (ARM) and the United States Geological Survey (USGS) Analysis of Flows in Networks of CHannels (AFINCH) model, represent efficient and effective approaches to propagating discharge observations into basin-wide (including gaged and ungaged areas) runoff estimates, and may offer an opportunity to improve predictive models for simulating runoff to the Great Lakes. Additionally, the intercomparison revealed that the median of the simulations from non-assimilative models agrees well with assimilative models, suggesting that using a combination of different methodologies may be an appropriate approach for estimating runoff into the Great Lakes. We then applied one assimilative model (ARM) to the Lake Michigan basin and found that there was persistent reduction in the amount of precipitation that becomes runoff following 1998, corresponding to a period of persistent low Lake Michigan water levels. The study was conducted as a first phase of the Great Lakes Runoff Intercomparison Project, a regional binational collaboration that aims to systematically and rigorously assess a variety of models currently used (or that could readily be adapted) to simulate basin-scale runoff to the North American Laurentian Great Lakes.

FUJISAKI, A., H. Mitsudera, J. WANG, and M. Wakatsuchi. How does the Amur River discharge flow over the northwestern continental shelf in the Sea of Okhotsk? Progress in Oceanography 126:8-20 (DOI:10.1016/j.pocean.2014.04.028) (2014).

The paths of the Amur River discharge on the continental shelf in the Sea of Okhotsk are still unknown despite their significance in transporting dissolved and particulate iron. In this study, we conduct a coupled ice-ocean simulation for the northern Sea of Okhotsk from June 1998 to September 2000 to answer the question: Does the Amur River discharge deposit materials to the pathway of the dense shelf water? In a series of numerical experiments, we identified two routes (the western and eastern routes) that could transport the river water more than 100 km offshore over the northwestern continental shelf. The two routes share the clockwise gyre in the Sakhalin Gulf and the northeastward flow on the northwestern continental shelf. These features are connected through the westward jet along the slope from the Sakhalin Gulf (the western route) and the northward transport over the shelf break canyon (the eastern route). The river water, the dense shelf water, and the easterly wind are in a fine geophysical balance for those features, and all are required for the formation of the two routes. The model results show that these unique joint effects in the Sea of Okhotsk allow the Amur River discharge to be effectively transported over the northwestern continental shelf, unlike a general river discharge that flows along the coast, and deposit materials into the pathway of the dense shelf water.

GRONEWOLD, A.D. Great Lakes water levels. In Michigan State of the Great Lakes. Michigan Department of Environmental Quality, Office of the Great Lakes, Lansing, MI, 6 pp. (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130061.pdf

Changes in the water levels of the Great Lakes impact humans and environmental systems across a variety of time and space scales. Storm events, for example, can lead to damaging and life-threatening water level surges along the Great Lakes coastline that are not only greater than the tidal fluctuations of marine coastlines, but are also more difficult to predict.

GRONEWOLD, A.D., and C.A. STOW. Unprecedented seasonal water level dynamics on one of the Earth's largest lakes. Bulletin of the American Meteorological Society 95:15-17 (DOI:10.1175/BAMS-D-12-00194.1) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140011.pdf

The North American Great Lakes (Fig. 1) contain roughly 20% of the Earth’s unfrozen fresh surface water and cover a massive area (Lake Superior alone is the largest unfrozen freshwater surface on the planet). Water levels on the Great Lakes have been recorded continuously for more than 150 years, representing one of the longest sets of direct hydroclimate measurements. This dataset, synthesized by Quinn (1981) and Lenters (2001), among many others, indicates that water levels on each of the Great Lakes follow a strong seasonal pattern closely linked with the timing and magnitude of the major components of the regional water budget, with relatively low water levels in the winter months, rising water levels in the spring, and decreasing water levels in the late summer and early fall. Water-level measurements on Lake Erie during the 2011 and 2012 water years (October 2010 through September 2011, and October 2011 through September 2012, respectively), however, reflect dramatic and unexpected changes in the seasonal water level cycle and in the Great Lakes regional water budget.

GRONEWOLD, A.D., and C.A. STOW. Water loss from the Great Lakes. Science 343:1084-1085 (DOI:10.1126/science.1249978) (2014). https://www.sciencemag.org/content/343/6175/1084.full?ijkey=zuj0Taz7dtJI6&keytype=ref&siteid=sci

As marine coastal populations experience and plan for rising ocean levels, residents along the coasts of Earth’s largest lake system are encountering the opposite problem: persistent low water levels and a receding shoreline. In January 2013, federal agencies from the United States and Canada documented the lowest water levels ever recorded on lakes Michigan and Huron. Only 6 years earlier, historically low water levels were recorded on Lake Superior, which feeds into the Lake Michigan-Huron system. These low water levels are symptoms of an imbalance in the water budget of the Great Lakes. Adapting to, and potentially mitigating, low water level conditions requires improved quantification of the factors that drive the imbalance.

HAWLEY, N., T. Redder, R. BELETSKY, E. Verhamme, D. BELETSKY, and J.V. DePinto. Sediment resuspension in Saginaw Bay. Journal of Great Lakes Research 40(Supplement 1):18-27 (DOI:10.1016/j.jglr.2013.11.010) (2014).

An integrated hydrodynamic and sediment transport model was applied to Saginaw Bay for the ice-free portions of 2009 and 2010. Observations of surface waves and suspended sediment concentration made during the spring of both years were used to constrain the model and to validate the model output. The results show that sediment resuspension in both the inner and outer bay is due almost entirely to surface wave action, and that the bulk of the resuspension events occur during the fall of each year. Although the model accurately predicted the occurrence of resuspension events, it did not always accurately simulate the amount of material resuspended. Because resuspension mixes bottom sediment into the water column and makes it and associated nutrients available to the biota, the effects of sediment resuspension need to be accounted for in any water quality model of the bay. Better specification of both the surface waves and the initial specification of the bottom sediment would probably improve the performance of the model.

He, C., L. Zhang, C. DeMarchi, and T.E. CROLEY II. Estimating point and nonpoint sources pollution loadings in the Saginaw Bay watersheds. Journal of Great Lakes Research 40(Supplement 1):11-17 (DOI:10.1016/j.jglr.2014.01.013) (2014). https://dx.doi.org/10.1016/j.jglr.2014.01.013

Databases of point sources including combined sewer overflows (CSOs) were acquired from the governmental agencies to map the occurrences and magnitude of the CSOs. Multiple databases of land use, topography, hydrography, soils, and agricultural statistics were used to estimate nonpoint source loading potential in the Saginaw Bay Basin, Michigan. Animal manure production was computed from tabulations of animals by 5-digit zip code area for the census years of 1987, 1992, 1997, and 2002. Fertilizer applications for both urban and agricultural land uses were calculated from county fertilizer estimates for the same periods. Results indicate that point sources from municipalities, industrial sectors and business entities contribute approximately 25% of the total phosphorus load to Saginaw Bay, with the remainder being accounted for by nonpoint source contributions. While the total amount of nutrients (N and P) from animal manure and fertilizer applications and atmospheric deposition declined in the Saginaw Bay Basin, fertilizer applications in non-farmland increased significantly. Estimation of nutrient loading potential at 5-digit zip code level reveals more detailed spatial variation and critical areas of nutrient loading than county level data for implementation of targeted water quality programs.

Hefferman, J.B., P.A. Soranno, M.J. Angilletta, L.B. Buckley, D.S. Gruner, T.H. Keitt, J.R. Kellner, J.S. Kominoski, A.V. Rocha, J. Xiao, T. Harms, S.J. Goring, L.E. Koenig, W.H. McDowell, H. Powell, A.D. Richardson, C.A. STOW, R. Vargas, and K.C. Weathers. Macrosystems ecology: Understanding ecological patterns and processes at continental scales. Frontiers in Ecology and the Environment 12(1):5-14 (DOI:10.1890/130017) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140005.pdf

Macrosystems ecology is the study of diverse ecological phenomena at the scale of regions to continents and their interactions with phenomena at other scales. This emerging subdiscipline addresses ecological questions and environmental problems at these broad scales. Here, we describe this new field, show how it relates to modern ecological study, and highlight opportunities that stem from taking a macrosystems perspective. We present a hierarchical framework for investigating macrosystems at any level of ecological organization and in relation to broader and finer scales. Building on well-established theory and concepts from other subdisciplines of ecology, we identify feedbacks, linkages among distant regions, and interactions that cross scales of space and time as the most likely sources of unexpected and novel behaviors in macrosystems. We present three examples that highlight the importance of this multiscaled systems perspective for understanding the ecology of regions to continents.

JOHENGEN, T.H., H.A. VANDERPLOEG, and J.R. LIEBIG. Effects of algal composition, seston stoichiometry, and feeding rate on zebra mussel (Dreissena polymorpha) nutrient excretion in two Laurentian Great Lakes. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 445-459 pp. (2013).

Nutrient excretion rates of zebra mussels from Saginaw Bay, Lake Huron, and Hatchery Bay, western Lake Erie, were compared over a 3-year period. Rates were estimated from nutrient accumulation in bottles of filtered lake water to which mussels were added for short (3-4 h) incubations. Excretion rates of soluble reactive phosphorus and ammonium were related to gross and net feeding rates measured in companion feeding experiments. Both feeding and excretion rates were related to seston quantity, seston nutrient concentrations, algal composition, and temperature. Excretion rates were significantly different between lakes and were significantly correlated to net feeding (assimilation) rate and algal composition as measured by the biomass percentage of cryptophytes, flagellates, and greens. Excretion rates were also significantly correlated to soluble reactive phosphorus concentrations and the assimilation rate of phosphorus, based on chlorophyll assimilation rate and phosphorus content of the seston. Nitrogen and phosphorus excretion rates showed distinct seasonal patterns and distinct differences between lakes. Notably, nitrogen excretion was most strongly correlated to temperature, while phosphorus excretion was most strongly related to assimilation of chlorophyll and phosphorus.

Kao, Y.C., S.A. Alderstein, and E.S. RUTHERFORD. The relative impacts of nutrient loads and invasive species on a Great Lakes food web: An Ecopath with Ecosim analysis. Journal of Great Lakes Research 40(Supplement 1):35-52 (DOI:10.1016/j.jglr.2013.11.012) (2014).

Excessive nutrient loads and species invasions pose significant threats to productivity and function of Great Lakes aquatic ecosystems. We used an Ecopath and Ecosim model to analyze impacts of changes in phosphorus loads, and dreissenid mussels and alewife biomass on the Saginaw Bay food web, Lake Huron. We configured the food web model in Ecopath with pre-dreissenids (1990) data on organism biomass, production, consumption, and diet from federal and state agency surveys and other sources. We conducted 70-year simulations in Ecosim of single factors (nutrients, alewives, and dreissenids) and their combinations. Phosphorus load scenarios were run under high (1970s), average (current), and low (target) levels; alewife scenarios were run under 2-fold the 1990�2003 average biomass, 1990�2003 average biomass, and alewife absence; dreissenid scenarios were run under the 1990�1996 average biomass, current (2009�2010) biomass, and dreissenid absence. Results indicated phosphorus loads were positively correlated with simulated biomass of most food web groups, alewife biomass were negatively correlated with biomass of most fish groups and macrozooplankton, while dreissenids impacts were most severe on lower trophic levels, but were relatively minor for fish groups compared with nutrient and alewife scenarios. Dreissenids had little effect on fish because Chironomids, the main component of fish diets that fed mostly on detritus, were not affected by dreissenids. Our results suggest that under current conditions of absence of alewives and reduced dreissenid biomass, obtaining target nutrient loads established in 1978 would not sustain current fishery harvests in Saginaw Bay given food web changes caused by invasive species.

Lavrentyev, P.J., H.A. VANDERPLOEG, G. Franze, D.H. Chacin, J.R. LIEBIG, and T.H. JOHENGEN. Microzooplankton distribution, dynamics, and trophic interactions relative to phytoplankton and quagga mussels in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 40(Supplement 1):95-105 (DOI:10.1016/j.jglr.2013.11.012) (2014).

Invasive quagga mussels have recently replaced zebra mussels as the dominant filter-feeding bivalves in the Great Lakes. This study examined microzooplankton (i.e., grazers b200 μm) and their trophic interactions with phytoplankton, bacteria, and bivalve mussels in Saginaw Bay, Lake Huron, following the zebra to quagga mussel shift. Microzooplankton distribution displayed strong spatial and temporal variability (1.73–28.5 μg C/L) relative to phytoplankton distribution. Ciliates were the dominant component, especially in the spring and early summer. Rotifers and dinoflagellates increased toward late summer/fall in the inner and outer parts of the bay, respectively. Microzooplankton grazing matched bacterial growth rates and removed ca. 30% of the phytoplankton standing stock in the b100 μm size fraction per day. The greatest herbivory occurred at the site dominated by colonial cyanobacteria. Microzooplankton, which comprised b4% of the quagga mussels prey field (i.e. available prey), contributed 77% and 34% to the quagga carbon-based diet during Microcystis and diatom blooms, respectively. Feeding on microzooplankton could buffer mussels during lean periods, or supplement other consumed resources, particularly during noxious cyanobacterial blooms. The results of this study demonstrate that microzooplankton are a resilient and critical component of the Saginaw Bay ecosystem.

Lesht, B.M., and N. HAWLEY. Video Clip 5: Zebra mussel movements on the bottom of Lake Michigan. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 755 pp. (2013).

This movie on the companion CD was created by concatenating a sequence of images made with a modified first-generation, time-lapse digital camera (Kodak DC260) controlled with a timing system of our own design. The camera was located on a bottom-resting tripod in Lake Michigan. Our original purpose for including a time-lapse camera among instruments mounted on the tripod was to document occurrences of sediment transport under stratified conditions in Lake Michigan. Although very little sediment movement occurred during the observation period (July 31 - August 26, 2002), the tripod happened to be positioned over preexisting clusters of zebra mussels, and our photographs captured what we believe to be the first long-term in situ observations of zebra mussel movements in the Great Lakes. As such, these observations complement laboratory studies described by Toomey et al. (2002) and descriptions of dreissenid movements given in Chapter 23 by Kobak (2013). Details of the deployment and brief comments on observations follow.

Li, J., Y. Qi, Z. Jing, and J. WANG. Enhancement of eddy-Ekman pumping inside anticyclonic eddies with wind-parallel extension: Satellite observations and numerical studies in the South China Sea. Journal of Marine Systems 132:150-161 (DOI:10.1016/j.jmarsys.2014.02.002) (2014).

Eddy-Ekman pumping was recently proposed as a potential mechanism to explain the high chlorophyll (Chl) concentrations inside anticyclonic eddies (AEs). In this paper, the influence of eddy-shape changes on the strength of eddy-Ekman pumping is investigated for AEs in the South China Sea (SCS) through multi-satellite data and numerical experiments. The results mainly indicate that wind-parallel (-perpendicular) extensions may enhance (reduce) eddy-Ekman pumping inside AEs. First, an AE detected in northern SCS winter shows significant Chl-wind relation and uplifted near-surface isopycnals when its shape extension becomes parallel with the wind direction. In addition, satellite observations indicate significant eddy-Ekman-pumping signals appear in the central-northern SCS during the winter monsoon and in the central-eastern SCS during the winter and summer monsoons, whenever and wherever AEs' shapes show wind-parallel extensions. This regional–seasonal characteristic is further confirmed by the numerical experiments, which display a positive linear relationship between wind-parallel projective length and the strength of eddy-Ekman pumping. This study may contribute to extend the theory of eddy-Ekman pumping with emphasizing the importance of eddy shapes and wind directions, especially in the situations that high spatial-temporal variability appears in eddies and/or the wind field.

Li, Y., S. Peng, J. WANG, and J. Yan. Impacts of non-breaking wave-stirring-induced mixing on the upper ocean thermal structure and typhoon intensity in the South China Sea. Journal of Geophysical Research: Oceans 119(8):5052-5070 (DOI:10.1002/2014JC009956) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140046.pdf

To investigate the effect of nonbreaking wave-induced mixing caused by surface wave stirring on the upper ocean thermal structure (UOTS) and the typhoon intensity, a simple nonbreaking wave stirring- induced mixing parameterization (WMP) scheme is incorporated into a regional coupled atmosphere-ocean model for the South China Sea (SCS), which couples the Princeton Ocean Model (POM) to the Weather Research and Forecasting (WRF) model using the OASIS3 coupler. The results of simulating two selected typhoon cases indicate that the nonbreaking wave-stirring-induced mixing has significant impacts on UOTS and the typhoon intensity, and the incorporation of the simple WMP scheme in the coupled model helps to improve the simulation of UOTS and thus the typhoon intensity. In the case that the typhoon intensity is underestimated by the atmosphere model alone, the improvement of initial UOTS by the ocean model with the WMP included can deepen the initial thermocline depth, reduce the effect of SST cooling, and prevent the typhoon intensity from undesired weakening. In the case that the typhoon intensity is overestimated (with strong winds), including the WMP in the ocean model significantly enhances the total vertical mixing rate in the upper ocean, which in turn enhances the SST cooling and thus reduces the typhoon intensity as desired. The results obtained in this study make a contribution to the ongoing efforts of improving the typhoon intensity forecast using a regional atmosphere-ocean coupled model by worldwide researchers and forecasters, especially for the typhoons in the SCS regions.

LIU, P.C. In search of realism on the study of freaque ocean waves. Proceedings, 35th Ocean Engineering Conference, Taiwan, China, November 21-22, 2013. Taiwan Society of Ocean Engineering, 115-118 pp. (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130054.pdf

In this presentation I expect to discuss the lesser known or underexplored aspects of the phenomenon of freaque ocean waves, what realities that might exist and what challenges we are facing. The current study of freaque waves has been an active research field over the last two decades or more. There have been significant advancements especially in connection with the study of nonlinear physics. Upon exploring the tangibility of what we do or do not know, there are still unrelenting challenges remain to be delved into.

LIU, P.C., J.C. Nieto Borge, G. Rodriquez, K.R. MacHutchon, and H.S. Chen. From single point gauge to spatio-temporal measurement of ocean waves - prospects and perspectives. Proceedings of the ASME 2014 33rd International Conference on Ocean, Offshore and Arctic Engineering, San Francisco, California, June 8-13, 2014. OMAE2014, 6 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140027.pdf

With the recent advancement of spatial measurements of ocean waves, we are clearly facing new challenges regarding how to handle an expanded new data system when it becomes widely available. In this paper we wish to present a preliminary attempt at confronting these prospects. Because the data is still very limited at present and also conceptually new, it’s a new, unfamiliar, and unrelenting world to pursue. We need a paradigm shift away from our familiar single-point conceptualization in order to effective approach the new world of truly spatial ocean waves.

LOFGREN, B.M., and A.D. GRONEWOLD. Reconciling alternative approaches to projecting hydrologic impacts of climate change. Bulletin of the American Meteorological Society 94(10):ES133-ES135 (DOI:10.1175/BAMS-D-13-00037.1) (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130064.pdf

Progress in the projection of hydrologic impacts of anthropogenic climate change (ACC) call for reassessing and documenting the state of the art. At this workshop, we placed a particular emphasis on understanding how consistency in the surface energy budget can be maintained, or lost, depending on how general circulation models’ internal calculations of hydrologic variables mesh with off line hydrologic models driven by their climatic variables (Lofgren et al. 2011; Sheffield et al. 2012). The science of projecting hydrologic impacts of climate change links the disciplines of meteorology (applied to long-term climate) and hydrology, yet these disciplines typically employ different perspectives on surface hydrological processes, leading to different modeling methods and means of linking models (Gronewold and Fortin 2012). To begin to resolve these inconsistencies, the workshop was structured around the following organizing questions: 1) How do we bridge the gap between climate projection and hydrologic projection? 2) How do we serve the data needs of the hydrologic and meteorological communities in a mutually consistent way? 3) What is the role of empirical and process-based models in a nonstationary regime? 4) How do we educate researchers and the general public about relevant caveats in simulations of hydrologic impacts of climate change?

LOFGREN, B.M., and A.D. GRONEWOLD. Chapter 12: Water resources. In Climate Change in the Midwest: A Synthesis Report for the National Climate Assessment. J.A. Winkler, J.A. Andresen, J.L. Hatfield, D. Bidwell, and D. Brown (Eds.). Island Press, 224-237 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140033.pdf

The water resources of the Midwest, and how they are managed under a future climate, have a significant collective impact on multiple economic sectors in the U.S., North America, and the world. The North American Laurentian Great Lakes, for example, hold nearly 20% of the earth’s accessible surface fresh water supply and have a coastline, and a coastal population, on the same order of magnitude as frequently-studied ocean coasts around the world (Fuller et al. 1995). In light of growing demands for clean water, access to coastal resources, and an improved understanding of climate dynamics in the Midwest region, a significant amount of research has recently been focused on understanding climate impacts on the lakes (both large and small), rivers, and streams in this region.

LOFGREN, B.M., A.D. GRONEWOLD, A. ACCIAIOLI, J. Cherry, A.L. Steiner, and D. Watkins. Methodological approaches to projecting the hydrologic impacts of climate change. Earth Interactions 17:19 pp. (DOI:10.1175/2013EI000532.1) (2013).

Climate change due to anthropogenic greenhouse gases (GHG) is expected to have important impacts on water resources, with a variety of societal impacts. Recent research has shown that applying different methodologies to assess hydrologic impacts can lead to widely diverging projections of water resources. The authors classify methods of projecting hydrologic impacts of climate change into those that estimate potential evapotranspiration (PET) based on air temperature and those that estimate PET based on components of the surface energy budget. In general, air temperature–based methods more frequently show reductions in measures of water resources (e.g., water yield or soil moisture) and greater sensitivity than those using energy budget–based methods. There are significant trade-offs between these two methods in terms of ease of use, input data required, applicability to specific locales, and adherence to fundamental physical constraints: namely, conservation of energy at the surface. Issues of uncertainty in climate projections, stemming from imperfectly known future atmospheric GHG concentrations and disagreement in projections of the resultant climate, are compounded by questions of methodology and input data availability for models that connect climate change to accompanying changes in hydrology. In the joint atmospheric–hydrologic research community investigating climate change, methods need to be developed in which the energy and moisture budgets remain consistent when considering their interaction with both the atmosphere and water resources. This approach should yield better results for both atmospheric and hydrologic processes

Lottig, N.R., T. Wagner, E.N. Henry, K.S. Cheruvelil, K.E. Webster, E.J. Downing, and C.A. STOW. Long-term citizen-collected data reveal geographical patterns and temporal trends in lake water clarity. PLOS One 9(4):1-8 (DOI:10.1371/journal.pone.0095769) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140020.pdf

We compiled a lake-water clarity database using publically available, citizen volunteer observations made between 1938 and 2012 across eight states in the Upper Midwest, USA. Our objectives were to determine (1) whether temporal trends in lake-water clarity existed across this large geographic area and (2) whether trends were related to the lake-specific characteristics of latitude, lake size, or time period the lake was monitored. Our database consisted of .140,000 individual Secchi observations from 3,251 lakes that we summarized per lake-year, resulting in 21,020 summer averages. Using Bayesian hierarchical modeling, we found approximately a 1% per year increase in water clarity (quantified as Secchi depth) for the entire population of lakes. On an individual lake basis, 7% of lakes showed increased water clarity and 4% showed decreased clarity. Trend direction and strength were related to latitude and median sample date. Lakes in the southern part of our study-region had lower average annual summer water clarity, more negative long-term trends, and greater interannual variability in water clarity compared to northern lakes. Increasing trends were strongest for lakes with median sample dates earlier in the period of record (1938–2012). Our ability to identify specific mechanisms for these trends is currently hampered by the lack of a large, multi-thematic database of variables that drive water clarity (e.g., climate, land use/cover). Our results demonstrate, however, that citizen science can provide the critical monitoring data needed to address environmental questions at large spatial and long temporal scales. Collaborations among citizens, research scientists, and government agencies may be important for developing the data sources and analytical tools necessary to move toward an understanding of the factors influencing macro-scale patterns such as those shown here for lake water clarity.

Lu, Y., P.A. Meyers, J.A. ROBBINS, B.J. EADIE, N. HAWLEY, and K.H. Ji. Sensitivity of sediment geochemical proxies to coring location and corer type in a large lake: Implications for paleolimnological reconstruction. Geochemistry, Geophysics, Geosystems 15:17 pp. (DOI:10.1002/2013GC004989) (2014).

We compared a suite of geochemical proxies in sediment cores collected in 1982, 1988, 1991, and 2003 from sites near the depocenter of Lake Erie to evaluate the reliability of paleoenvironmental reconstructions derived from lacustrine sediments. Our proxies included the concentrations and carbon isotopic compositions of organic and inorganic carbon (TOC, CaCO3, d13Corg, and d13CCaCO3), augmented by organic C to total N ratios (Corg:Ntot), d15N, and carbonate d18O values (d18OCaCO3). The three coring sites were clustered within 12 km; two types of corers—a Box corer and a Benthos gravity corer—were used for the 1991 sampling campaign. The variance of most proxies was accounted for not only by temporal environmental changes but also by coring locations and corer type, indicating that sediment spatial heterogeneity and differences in sediment recovery due to the use of different corers also played a part in determining the geochemical compositions of these cores. The TOC, d13Corg, and d13CCaCO3 values showed decadal temporal patterns that were consistent between the multiple sampling campaigns. In contrast, the d15N, Corg:Ntot, CaCO3, and d18OCaCO3 exhibited across-core differences in their temporal variations, making it difficult to extract consistent environment information from different cores. Our findings suggest that in addition to temporal environmental changes, high-resolution paleolimnological reconstruction is sensitive to many factors that could include spatial sediment heterogeneity, discontinuous sedimentation processes, bioturbation, sediment dating uncertainty, and artifacts associated with analytical and coring procedures. Therefore, multiple-core sampling and analysis are important in reliably reconstructing environmental changes, particularly for large, heterogeneous lacustrine basins.

Ma, J., Y. Du, H. Zhan, H. Liu, and J. WANG. Influence of oceanic Rossby waves on phytoplankton production in the southern tropical Indian Ocean. Journal of Marine Science 134:12-19 (DOI:10.1016/j.jmarsys.2014.02.003) (2014).

Using Sea-viewing Wide Field-of-view Sensor (SeaWiFS) ocean color data, we investigated the biological responses to oceanic Rossby waves in the southern tropical Indian Ocean (TIO) during 1997–2007. The findings indicate that, during the developing phase of El Niño/La Niña events, usually in boreal fall, and triggered by anomalous anticyclonic/cyclonic wind circulations in the southeast TIO, downwelling/upwelling Rossby waves form and then propagate westward. After a few months, downwelling/upwelling Rossby waves interface with the thermocline dome in the southern TIO, and suppress/enhance the upwelling. Correspondingly, less/more nutrient-rich waters enter the mixed layer and result in lower/higher chlorophyll concentrations. Due to the asymmetric effects on the thermocline dome between downwelling and upwelling Rossby waves, higher chlorophyll concentrations appear earlier and in the eastern part of the dome, whereas lower chlorophyll concentrations appear later and in the central part of the dome. Moreover, when El Niño/La Niña–Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) events coincide, the biological responses are stronger.

Madenjian, C.P., E.S. RUTHERFORD, C.A. STOW, E.F. Roseman, and J.X. He. Trophic shift, not collapse. Environmental Science & Technology 47:11915-11916 (DOI:10.1021/es404089y) (2013). https://pubs.acs.org/doi/full/10.1021/es404089y

Jerald Schnoor’s editorial describes the recent changes in Lake Huron’s aquatic ecosystem as a trophic collapse and attributes this collapse to invasive species dominating energy and nutrient flows in the food web. As state and federal scientists who are closely monitoring Lake Huron’s food web, we believe that the ongoing changes are more accurately characterized as a trophic shift in which benthic pathways have become more prominent. While decreases in abundance have occurred for some species, others are experiencing improved reproduction resulting in the restoration of several important native species.

Marko, K.M., E.S. RUTHERFORD, B.J. EADIE, T.H. JOHENGEN, and M.B. LANSING. Delivery of nutrients and seston from the Muskegon River watershed to nearshore Lake Michigan. Journal of Great Lakes Research 39(4):672-681 (DOI:10.1016/j.jglr.2013.08.002) (2013).

Drowned river mouth lakes are major features of coastal Great Lakes habitats and may influence nutrient and organic matter contributions from watersheds to near shore coastal zones. In May through October 2003, we measured loads of nutrients, surficial sediment, and seston to track the delivery of riverine-derived materials from the lower Muskegon River Watershed (MRW) into the near shore area of southeast Lake Michigan. Nutrient flux data indicated that seasonal loads of 1800 metric tons (MT) of particulate organic carbon, 3400 MT of dissolved organic carbon, and 24 MT of total phosphorus were discharged from the lower Muskegon River, with approximately 33% of TP load and 53% of the POC load intercepted within the drowned river mouth terminus, Muskegon Lake. Carbon: phosphorus molar ratios of seston in Muskegon River (C:P = 187) and Muskegon Lake (C:P = 176) were lower than in Lake Michigan (C:P = 334), indicating phosphorus limitation of phytoplankton in near shore Lake Michigan. Isotopic signatures of seston collected in Muskegon Lake were depleted in δ13C (−30.8 ± 1.6‰) relative to the isotope signatures of seston from Lake Michigan (−26.2 ± 1.3‰) or the mouth of the Muskegon River (−28.1 ± 0.5‰), likely due to the presence of biogenic methane in Muskegon Lake. Seston δ15N increased on a strong east-to-west gradient within Muskegon Lake, indicating significant microbial processing of nutrients. The extent of nutrient uptake in Muskegon Lake altered the chemical and isotopic characterization of seston flowing into Lake Michigan from Muskegon River.

MIEHLS, A.L.J., A.G. McAdam, P.E. Bourdeau, and S.D. Peacor. Plastic response to a proxy cue of predation risk when direct cues are unreliable. Ecology 94(10):2237-2248 (DOI:10.1890/12-2250.1) (2013).

Responses to proximate cues that directly affect fitness or cues directly released by selective agents are well-documented forms of phenotypic plasticity. For example, to reduce predation risk, prey change phenotype in response to light level (e.g., moon phase) when light affects predation risk from visual predators, and to chemical cues (kairomones) released by predators. Less well understood is the potential for organisms to perceive predation risk through ‘‘proxy cues’’: proximate cues that correlate with, but do not directly affect predation risk. Previous field studies indicate that body and spine length of an invasive cladoceran in Lake Michigan, Bythotrephes longimanus (the spiny water flea), increase during the growing season, coincident with a decrease in clutch size. Although the cause of seasonal trait changes is not known, changes are associated with warmer water temperature and increased predation risk from gape-limited fish (i.e., fish whose ability to consume Bythotrephes is limited by mouth size). Using a laboratory experiment, we found no effect of fish (Perca flavescens) kairomones on Bythotrephes morphology or life history. In contrast, higher water temperature led to longer absolute spine and body length, increased investment in morphological defense of offspring (measured as the ratio of spine-to-body length), and decreased clutch size and age at reproduction. These plastic responses are unlikely to be adaptive to temperature per se, but rather our findings indicate that temperature serves as a proxy cue of fish predation risk. Temperature correlates with risk of gape-limited fish predation due to growth of fish from larval stages incapable of consuming Bythotrephes early in the season, to larger sizes by midseason increasingly capable of consuming Bythotrephes, but limited by gape size to consuming smaller individuals. We argue that for Bythotrephes, temperature is a more reliable cue of predation risk than fish kairomones, because fish kairomones are present throughout the season due to continual presence of non-gape-limited adult fish, to which plastic response would have little effect. Organisms may, therefore, not only respond to changes in an environmental factor because the factor directly affects risk, but also when the environmental factor serves as a proxy signaling change in predation risk.

MIEHLS, A.L.J., S.D. Peacor, and A.G. McAdam. Gape-limited predators as agents of selection on the defensive morphology of an invasive invertebrate. Evolution 68(9):2633-2643 (DOI:10.1111/evo.12472) (2014).

Invasive species have widespread and pronounced effects on ecosystems and adaptive evolution of invaders is often considered responsible for their success. Despite the potential importance of adaptation to invasion, we still have limited knowledge of the agents of natural selection on invasive species. Bythotrephes longimanus, a cladoceran zooplankton, invaded multiple Canadian Shield lakes over the past several decades. Bythotrephes have a conspicuous caudal process (tail spine) that provides a morphological defense against fish predation. We measured viability selection on the longest component of the Bythotrephes spine, the distal spine segment, through a comparison of the lengths of first and second instar Bythotrephes collected from lakes differing in the dominance of gape-limited predation (GLP) and nongape-limited predation (NGLP) by fish. We found that natural selection varied by predator gape-limitation, with strong selection (selection intensity: 0.20–0.79) for increased distal spine length in lakes dominated by GLP, and no significant selection in lakes dominated by NGLP. Further, distal spine length was 17% longer in lakes dominated by GLP, suggesting the possibility of local adaptation. As all study lakes were invaded less than 20 years prior to our collections, our results suggest rapid divergence in defensive morphology in response to selection from fish predators.

MILLER, R., N. HAWLEY, and S.A. RUBERG. Video Clip 4: Visual documentation of quagga mussels (Dreissena rostriformis bugensis) at two depths in southeastern Lake Michigan. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 753 pp. (2013).

This video on the companion CD documents the extent of bottom coverage by quagga mussels (Dreissena rostriformis bugensis) at two sites in southeastern Lake Michigan near Muskegon, Michigan, in 2012. The sites were located at about 40 m (43 °11.259' N, 86°25.507' W) and 110 m (43°13.263' N, 86°34.9' W) along a depth-defined transect that is sampled regularly as part of a long-term monitoring program of the lower food web (Fahnenstiel et al. 2010). The video was obtained with cameras mounted on a Model 1000 Outland remotely operated vehicle (ROV) that collected high-definition 1080p video at 30 fps at the 40 m site, and standard-definition video at the 110m site. Videos at both sites were collected in the vicinity of research moorings that provide physical data (i.e., currents, temperature) as a component to the monitoring program.

NALEPA, T.F., J.S. Houghton, R. Paddock, and J. Janssen. Video Clip 3: Close-up view of inhalant siphons of quagga mussels (Dreissena rostriformis bugensis, deepwater morph) on the Midlake Reef Complex in Lake Michigan. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 751 pp. (2013).

This video clip on the companion CD provides a close-up view or the inhalant siphons of quagga mussels (Dreissena rostriformis bugensis, deepwater morph) found on the Midlake Reef Complex (MLRC) in Lake Michigan and is meant to supplement Nalepa et al. (2013). In Nalepa et al. (2013), relative siphon lengths (siphon length/shell length ratio) were determined on mussels observed in the laboratory and estimated from still-frame images of mussels observed in this video clip. The video was obtained with a remotely operated vehicle (ROV; Benthos MiniRover MK-11) in May 2010 at locations between 53 and 57 m on Northeast Reef (located about 43°15.0' N and 87°35.0' W), which is a part of the reef complex (see Houghton et al. 2013). The video begins with images of the underside of a rock ledge (0:0-0:23) in which attached Hydra can be seen extended downward and moving with the current. It then pans upward to show quagga mussels on the vertical faces of the same ledge (0:44-1:53). As evidenced by unoccupied spaces between individuals, mussel densities at this location were not as great as in other reef locations.

NALEPA, T.F., V. Pavlova, W.H. Wong, J. Janssen, J.S. Houghton, and K. MABREY. Variation in the quagga mussel (Dreissena rostriformis bugensis), with emphasis on the deepwater morphotype in Lake Michigan. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 315-329 pp. (2013).

Two morphological features (relative length of the inhalant siphon to shell length and relative shell height to shell length) of the deepwater morphotype of the quagga mussel (Dreissena rostriformis bugensis) were examined among individuals collected at three depths (25, 45, and 93 m) in Lake Michigan. Based on laboratory observations, mean siphon length to shell length ratios were 0.19, 0.44, und 0.20 at the three depths, respectively, and the mean ratio was significantly greater at 45 m than at the other two depths. Mean ratio of in situ deepwater morphs estimated from still-capture frames of underwater video taken at 53-57 m in Lake Michigan was 0.74 but ranged from 0.57 lo 0.78 depending on the location. In comparison, the ratio of quagga mussel shallow morphs from western Lake Erie, as observed in the laboratory, was 0.19. We speculate that longer siphons relative to shell lengths at 45 m are an adaptation to greater mussel densities at this depth compared to 25 and 93 m. Siphon lengths of in situ mussels also appeared to be longer relative to shell lengths at locations where mussels were more dense. A longer siphon provides more flexibility in placement and would thus offer a competitive advantage in obtaining food when population densities are high. The mean shell width to shell length ratio of the deepwater morph in Lake Michigan was 0.52, and there was no difference in ratios between the three depths. Water temperatures at the three depths in Lake Michigan were broadly different during seasonal stratification (June to October) and, hence, our results do not agree with previous studies that showed water temperature plays an important role in width to length ratios in shells of quagga mussels.

Nash, K.L., C.R. Allen, D.G. Angeler, C. Barichievy, T. Eason, A.S. Garmestani, N.A.J. Graham, D. Granholm, M. Knutson, R.J. Nelson, M. Nystrom, C.A. STOW, and S.M. Sundstrom. Discontinuities, cross-scale patterns, and the organization of ecosystems. Ecology 95(3):654-667 (DOI:10.1890/13-1315.1) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140016.pdf

Ecological structures and processes occur at specific spatiotemporal scales, and interactions that occur across multiple scales mediate scale-specific (e.g., individual, community, local, or regional) responses to disturbance. Despite the importance of scale, explicitly incorporating a multi-scale perspective into research and management actions remains a challenge. The discontinuity hypothesis provides a fertile avenue for addressing this problem by linking measureable proxies to inherent scales of structure within ecosystems. Here we outline the conceptual framework underlying discontinuities and review the evidence supporting the discontinuity hypothesis in ecological systems. Next we explore the utility of this approach for understanding cross-scale patterns and the organization of ecosystems by describing recent advances for examining nonlinear responses to disturbance and phenomena such as extinctions, invasions, and resilience. To stimulate new research, we present methods for performing discontinuity analysis, detail outstanding knowledge gaps, and discuss potential approaches for addressing these gaps.

Nguyen, T.D., P. Thupaki, E.J. ANDERSON, and M.S. Phanikumar. Summer circulation and exchange in the Saginaw Bay-Lake Huron system. Journal of Geophysical Research: Oceans 119:2713-2734 (DOI:10.1002/2014JC009828.) (2014). https://onlinelibrary.wiley.com/doi/10.1002/2014JC009828/abstract

We use a three-dimensional, unstructured grid hydrodynamic model to examine circulation and exchange in the Saginaw Bay-Lake Huron system during the summer months for three consecutive years (2009–2011). The model was tested against ADCP observations of currents, data from a Lagrangian drifter experiment in the Saginaw Bay, and temperature data from the National Data Buoy Center stations. Mean circulation was predominantly cyclonic in the main basin of Lake Huron with current speeds in the surface layer being highest in August. Circulation in the Saginaw Bay was characterized by the presence of an anticyclonic gyre at the mouth of the outer bay and two recirculating cells within the inner bay. New estimates are provided for the mean flushing times (computed as the volume of the bay divided by the rate of inflow) and residence times (computed as e-folding flushing times based on dye concentration modeling treating the bay as a continuously stirred tank reactor) for Saginaw Bay. The average flushing time (over the 3 months of summer and for all 3 years) was 23.0 days for the inner bay and 9.9 days for the entire bay. The mean e-folding flushing time was 62 days (2 months) for the inner bay and 115 days (3.7 months) for the entire bay for the summer conditions examined in this work. To characterize the behavior of river plumes in the inner Saginaw Bay, trajectory data from GPS-enabled Lagrangian drifters were used to compute the absolute diffusivity values in the alongshore and cross-shore directions.

Overland, J.E., J. WANG, R.S. Pickart, and M. Wang. Chapter 2. Recent and future changes in the meteorology of the Pacific Arctic. In The Pacific Arctic Region: Ecosystem Status and Trends in a Rapidly Changing Environment. J.M. Grebmeier and W. Maslowski (Eds.). Springer Science+Business Media, Dordrecht, 17-30 pp. (DOI:10.1007/978-94-017-8863-2_2) (2014).

The meteorology of the Pacific Arctic (the Bering Sea through the Chukchi and southern Beaufort Seas) represents the transition zone between the moist and relatively warm maritime air mass of the Pacific Ocean to the cold and relatively dry air mass of the Arctic. The annual cycle is the dominant feature shifting from near total darkness with extensive sea ice cover in winter to solar heating in summer that is equal to that of sub-tropic latitudes. Strong north-south gradients in air temperatures and sea level pressure are typical over the Pacific Arctic giving rise to climatological polar easterly winds (blowing from the east) throughout the year. Localized storms (regions of low sea level pressure) can propagate into the region from the south but high pressure regions are typical, connected to either northeastern Siberia or the southern Beaufort Sea. The northern portion of the Pacific Arctic has participated in the general Arctic-wide warming in all seasons over the last decade while the southern Bering Sea turned to near record cold temperatures after 2006. Future climate changes in the Pacific Arctic will come from shifts in the timing and extent of seasonal sea ice. Based on climate model projections, cold and dark conditions will still dominate over a climate warming scenario of +2° C by 2050. The northern Bering Sea will continue to have extensive sea ice through April, while the southern shelf will have on average less sea ice than currently observed but with large interannual variability. The largest change has the southern Chukchi Sea shifting from being sea ice free in September and October to becoming sea ice free from July through November within a decade or two, impacting shipping, oil exploration, and ecosystems.

POTHOVEN, S.A., and G.L. FAHNENSTIEL. Lake Michigan after dreissenid mussel invasion: Water quality and food web changes during the late winter/spring isothermal period. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 545-553 pp. (2013).

This chapter documents water quality and lower foodweb changes in the offshore region of southern Lake Michigan during the late winter/spring isothermal mixing period during three distinct decades: 1983-1989 (pre-dreissenids), 1995-1998 (post-zebra mussel/pre-quagga mussel), and 2007-2010 (post-quagga mussel). Total phosphorus, chlorophyll a, and phytoplankton primary productivity did not change or decreased minimally (<22% decrease) between 1983-1989 and 1995-1998 but decreased substantially between 1995-1998 and 2007-2009 (34%, 68%, and 71% decrease, respectively). Secchi depth transparency increased from 6 m in 1985-1989 and 7 m in 1995-1998 to 18 m in 2007-2010. These pronounced changes in water column properties during the isothermal period in 2007-2010 compared to earlier periods were primarily attributed to the filtering activities of the quagga mussel (Dreissena rovstriformis bugensis) and to a lesser extent to phosphorus load reductions. In contrast to substantial changes in phytoplankton, total zooplankton biomass did not differ between time periods (1986-1988, 1994-1998, and 2007- 2009). However, biomass of cyclopoid copepods declined to nearly negligible levels in 2007-2009, perhaps reflecting a more oligotrophic system that favored calanoid copepods. Trends in biomass of calanoid copepods were not readily apparent or possibly increased in 2007-2009. Densities of Mysis diluviana declined 60% between 1996-2000 and 2007-2010. These large changes in limnological parameters and the food web have contributed to significant changes in traditional nutrient stoichiometry: Increases in total phosphorus to chlorophyll ratios and decreases in phytoplankton carbon to chlorophyll ratios were noted in 2007-2009 as compared to 1983-1987 and 1995-1998. Our results show that the recent expansion of quagga mussels in Lake Michigan has affected almost all aspects of water column nutrient and foodweb dynamics during the late winter/spring isothermal mixing period. We believe that future actions to manage the Lake Michigan ecosystem must consider these unique impacts of dreissenids and such impacts will affect management actions elsewhere in the Great Lakes

POTHOVEN, S.A., T.O. Hook, and C.R. Roswell. Feeding ecology of age-0 lake whitefish in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 40(Supplement 1):148-155 (DOI:10.1016/j.jglr.2013.09.016) (2014).

Age-0 lake whitefish Coregonus clupeaformis (11–160 mm total length) were collected from Saginaw Bay, Lake Huron during April–November 2009 and 2010 for diet analysis and for the evaluation of ontogenetic changes in feeding ecology. Lake whitefish ≤ 50 mm ate mainly zooplankton, after which their diets switched mainly to benthic macroinvertebrates. Cyclopoida were the dominant prey consumed by very small lake whitefish (b17 mm) and the most frequently selected zooplankton type for individual small fish. Once lake whitefish reached 18–19 mm, Cyclopoida in the diet declined and cladocerans emerged as an important diet item. Daphnia were the most common cladoceran in the diets, but for fish 31–50 mm Bosminidae were also relatively important. Although the shift to Daphnia could represent an ontogenetic point when lake whitefish were large enough to effectively handle this prey, it also took place when the relative availability of Daphnia increased. Once lake  whitefish were N50 mm, Chironomidae larvae became a dominant prey item and this shift to benthivory coincided with a 55% increase in length-adjusted energy content between June and July. However, as fish grew (around 110–120 mm), Sphaeriidae and the benthic zooplankton Chydoridae became increasingly important in the diet. As these less energetically rich prey were incorporated into the diet, there were corresponding 21 and 15% decreases in length-adjusted energy content from July to August and September, respectively.

POTHOVEN, S.A., T.O. Hook, and C.R. Roswell. Energy content of young yellow perch and walleye in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 40(Supplement 1):133-138 (DOI:10.1016/j.jglr.2013.10.002) (2014).

We evaluated seasonal energy content of age-0 yellow perch Perca flavescens and walleye Sander vitreus in Saginaw Bay, Lake Huron in 2009 and 2010. We also determined the energy content of age-1 fish from the 2009 and 2010 cohorts the following spring (i.e., for fish that had survived one winter) to evaluate overwinter energy losses. As expected, larger fish within each species had disproportionately higher energy content (i.e., slope relating length and energy N3.0) than smaller conspecifics. By contrast to expectations, allometric slopes were N3.0 in nearly all months, not just the fall, and were higher for age-0 yellow perch than for walleye, even though increased allocation to growth would have seemingly been beneficial to even the largest yellow perch during summer. Seasonal energy allocation patterns differed between years. In 2009, length specific energy content increased from late summer to fall for both species. However, for the 2010 cohorts of fish, length specific energy content decreased between late summer and fall for yellow perch and did not change for walleye. There were 13–17% overwinter declines in length specific energy content between the fall (October or November) and the spring (May) with no major differences between cohorts within a species or between species for a given year. Because young yellow perch and walleye are similar physiologically but differ in size (i.e., yellow perch are smaller), it is possible that overwinter energy losses are more important for yellow perch than for walleye.

POTHOVEN, S.A., and C.P. Madenjian. Increased piscivory by Lake Whitefish in Lake Huron. North American Journal of Fisheries Management 33:1194-1202 (DOI:10.1080/02755947.2013.839973) (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130055.pdf

We evaluated the diet of Lake Whitefish Coregonus clupeaformis in Lake Huron during 2002–2011 to determine the importance of Round Goby Neogobius melanostomus and other fish as prey items. Lake Whitefish that had reached approximately 400 mm in length incorporated fish into their diets. The overall percentage of adult Lake Whitefish in Lake Huron that had eaten fish increased from 10% in 2002–2006 to 20% in 2007–2011, with a corresponding decrease in the frequency of Lake Whitefish that ate Dreissena spp. from 52% to 33%. During 2002–2006, Round Goby (wet mass, 38%), sculpins (Cottidae) (34%), and Ninespine Stickleback Pungitius pungitius (18%) were the primary fish eaten, whereas Round Goby accounted for 92% of the fish eaten in 2007–2011. Overall, Round Goby were found in the fewest Lake Whitefish stomachs in the north region of Lake Huron (6%) and in the most in the central (23%) and south (19%) regions of the lake. In the central region, Round Goby were eaten during all seasons that were sampled (spring through fall). In the south region, Round Goby were eaten only in the winter and spring but not in the summer when Dreissena spp. and spiny water flea Bythotrephes longimanus dominated the diet. Based on the 2007–2011 diet composition, an individual Lake Whitefish would need to have increased their consumption relative to that in 1983–1994 by 6% in the north region, 12% in the central region, and 41% in the southern region in order to achieve the same growth that was observed before dreissenid mussels arrived. However, Lake Whitefish weight adjusted for length only increased by 2% between 2002–2006 and 2007–2011 in the central region, decreased by 4% in the northern region, and remained constant in the southern region. This suggests that a shift toward more frequent piscivory does not necessarily improve the condition of a generalist feeder like Lake Whitefish.

REID, D.F., and D. Wilkinson. Catalyst for change: The little dreissenid that did (change national policy on aquatic invasive species). In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 177-185 pp. (2013).

The earliest and most extensive legislation that established laws for prevention and control of nonindigenous and invasive nuisance) species in the United States was for protection of crops and livestock. Federal legislation to protect aquatic ecosystems was absent until the late-twentieth century, and federal policies most often took a "benign neglect" approach, sometimes even encouraging introductions of new species in attempts to make natural ecosystems more responsive to human needs. Policies against species introductions were rare and primarily reactive to individual species that became nuisances after they were established. Each new species invasion was treated as a separate event. The discovery of the zebra mussel (Dreissena polymorpha) in 1988 led to recognition that aquatic invasions were also a problem. While discoveries of nonindigenous species raised concerns prior to the late-1980s, it was the discovery of the zebra mussel that catalyzed awareness and rapid policy changes, including establishment of new regulations at many levels of government across the Unites States and Canada. In the United States, the centerpiece of legislation was the Nonindigenous Aquatic Nuisance Prevention and Control Act, passed in 1990 and reauthorized and expanded as the National Invasive Species Act in 19%. It included legislative mandates for prevention (emphasizing ballast water), zebra mussel control, and other actions such as ecosystem surveys. It established the Aquatic Nuisance Species Task Force to encourage a coordinated federal approach and provided for establishment of regional panels; ultimately, it provided the impetus for Executive Order 13412, which resulted in the first comprehensive national plan for managing all invasive species in the United States. Although little additional progress was made after 1996 through 2011, the policy changes started in 1990 established aquatic invasive species as recognized national problems with regional and local implications. The legislation of the 1990s led to growth of sustained federal and state programs that continue to respond to new introductions while working toward preventing future introductions.

Roswell, C.R., S.A. POTHOVEN, and T.O. Hook. Patterns of age-0 yellow perch growth, diets, and mortality in Saginaw Bay, Lake Huron. Journal of Great Lakes Research 40(Supplement 1):123-132 (DOI:10.1016/j.jglr.2014.01.008) (2014).

Identifying mechanisms influencing early-life survival may elucidate recruitment variability of fish populations. Yellow perch (Perca flavescens), are economically and ecologically important in Saginaw Bay, Lake Huron, but have recently experienced low recruitment despite strong production of age-0 fish. Recent year classes have been characterized by slow first-year growth, as indexed by fall size of age-0 yellow perch; however, seasonal growth patterns and accompanying diet and survivorship patterns have not been documented for age-0 yellow perch in Saginaw Bay. To this end, we collected age-0 yellow perch weekly (larvae) and monthly (juveniles) throughout the first year of life during 2009 and 2010 to track changes in growth and diet composition. We also evaluated predation and over-winter energy-loss as potential mechanisms of size-selective mortality. Yellow perch growth, energy accumulation and size-specific condition decreased during late summer and fall. During larval and juvenile stages, predominant components of yellow perch diets transitioned from copepods to Daphnia and other zooplankton; however, we observed only weak ontogenetic shifts toward benthic prey. Smaller yellow perch a) were preferentially preyed upon by walleye (the bay's main piscivore) and b) displayed lower mass-specific energy content, potentially increasing overwinter starvation risk, suggesting that slow growth increases mortality risk. Our results are consistent with the hypothesis that recruitment dynamics are influenced by an interplay of size-selective mortality and diet-induced reductions in growth.

Rucinski, D.K., J.V. DePinto, D. Scavia, and D. BELETSKY. Modeling Lake Erie's hypoxia response to nutrient loads and physical variability. Journal of Great Lakes Research 40:11 pp. (DOI:10.1016/j.jglr.2014.02.003) (2014).

A 1-dimensional, linked hydrodynamic and eutrophication model was developed and calibrated with 19 years of observations (1987–2005) for the summer stratification period in the central basin of Lake Erie, corroborated by comparison with observed process rates and areal hypoxic extents, and confirmed with observations from the 1960s and 1970s. The model effectively captures observations of both vertical and temporal trends in dissolved oxygen, as well as temporal trends in chlorophyll-a, phosphorus, zooplankton biomass, and several key processes. The model was used to develop a relationship between external phosphorus load and hypolimnion oxygen conditions, and then to establish load–response envelopes that account for inter-annual variability in physical conditions driven by variation in meteorological drivers. The curves provide a valuable tool for reassessing phosphorus loading targets with respect to reducing hypoxia in Lake Erie.

Ryan, D.J., T.F. NALEPA, L.N. IVAN, M.S. Sepulveda, and T.O. Hook. A comparison of consumptive demand of Diporeia spp. and Dreissena in Lake Michigan based on bioenergetics models. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 713-723 pp. (2013).

Historically, Diporeia spp. were abundant and the dominant benthic consumer in offshore regions of the Great Lakes. Over the past several decades, however, Diporeia populations have declined, seemingly in response to Dreissena spp. introduction and expansion. As a result, consumption demand of food resources has likely undergone large changes. Bioenergetics models are commonly used to estimate consumption of aquatic organisms from observed growth patterns and environmental conditions. A Diporeia-specific bioenergetics model was parameterized from literature values of Diporeia growth, water temperature, energy density, and abundances in Lake Michigan. These data were used as input to derive model estimates of consumption. Specifically, the model was used to investigate Diporeia consumption of particulate matter in Lake Michigan in 1995, 2005, 2008, which covered a period of severe declines in Diporeia and increased abundances of dreissenids. Comparisons of Diporeia consumption before and after their declines indicated that consumption rate decreased by 91%. In contrast, conservatively-estimated consumption rates for dreissenids in Lake Michigan were up to 24 times greater than Diporeia consumption rates during 2005, and by 2008, consumption rates of dreissenids in Lake Michigan were more than 330 times greater than consumption by Diporeia before their decline in 1995. We hypothesize that higher abundance and consumption by dreissenid mussels, relative to Diporeia, altered the function of the Lake Michigan ecosystem.

Scavia, D., J.D. Allan, K.K. Arend, S. Bartell, D. BELETSKY, N.S. Bosch, S.B. Brandt, R.D. Briland, I. Daloglu, J.V. DePinto, D.M. Dolan, M.A. Evans, T.M. Farmer, D. Goto, H. Han, T.O. Hook, R. Knight, S.A. Ludsin, D.M. MASON, A.M. Michalak, R.P. Richards, J.J. Roberts, D.K. Rucinski, E.S. RUTHERFORD, D.J. SCHWAB, T. Sesterhenn, H. ZHANG, and Y. Zhou. Assessing and addressing the re-eutrophication of Lake Erie: Central basin hypoxia. Journal of Great Lakes Research 40:226-246 (DOI:10.1016/j.jglr.2014.02.004) (2014).

Relieving phosphorus loading is a key management tool for controlling Lake Erie eutrophication. During the 1960s and 1970s, increased phosphorus inputs degraded water quality and reduced central basin hypolimnetic oxygen levels which, in turn, eliminated thermal habitat vital to cold-water organisms and contributed to the extirpation of important benthic macroinvertebrate prey species for fishes. In response to load reductions initiated in 1972, Lake Erie responded quickly with reduced water-column phosphorus concentrations, phytoplankton biomass, and bottom-water hypoxia (dissolved oxygen b2 mg/l). Since the mid-1990s, cyanobacteria blooms increased and extensive hypoxia and benthic algae returned. We synthesize recent research leading to guidance for addressing this re-eutrophication, with particular emphasis on central basin hypoxia. We document recent trends in key eutrophication-related properties, assess their likely ecological impacts, and develop load response curves to guide revised hypoxia-based loading targets called for in the 2012 Great Lakes Water Quality Agreement. Reducing central basin hypoxic area to levels observed in the early 1990s (ca. 2000 km2) requires cutting total phosphorus loads by 46% from the 2003–2011 average or reducing dissolved reactive phosphorus loads by 78% from the 2005–2011 average. Reductions to these levels are also protective of fish habitat. We provide potential approaches for achieving those new loading targets, and suggest that recent load reduction recommendations focused on western basin cyanobacteria blooms may not be sufficient to reduce central basin hypoxia to 2000 km2.

Schloesser, D.W., J.A. ROBBINS, G.A. Matisoff, T.F. NALEPA, and N. MOREHEAD. A 200-year chronology of burrowing mayflies (Hexagenia spp.) in Saginaw Bay. Journal of Great Lakes Research 40(1):80-91 (DOI:10.1016/j.jglr.2013.12.016) (2014).

After an absence of 50 years, burrowing mayflies (Hexagenia spp.) colonized western Lake Erie which led to interest in whether this fauna can be used to measure recovery in nearshore waters throughout the Great Lakes. However, in many areas we do not know if mayflies were native/endemic and thus, whether recovery is a logical measure to assess progress of recovery. In the present study, we construct a chronologic record of relative abundance of burrowing mayflies in Saginaw Bay by the use of mayfly tusks and radionuclides in sediments (i.e., a paleoecologic record) and historic records of mayfly nymphs in the bay. These records reveal that mayflies: (1) were few before 1799, which indicates that nymphs were probably native/endemic in the bay, (2) increased between 1799 and 1807 and remained at relatively high levels between 1807 and 1965, probably in response to increased nutrient run-off from the watershed, (3) declined dramatically between 1965 and 1973, probably as a result of excessive eutrophication in the mid-1950s; and, (4) were few and highly variable between 1973 and 2001, probably as a result of low and unstable abundances of mayfly nymphs. Historic records verify that nymphs disappeared in the bay in the late-1950s to early-1960s which is in agreement with the paleoecologic record. Reoccurrence of low abundances of nymphs in the bay between 1991 and 2008 and comparison of chronologic records of nymphs in Saginaw Bay and western Lake Erie suggest that mayflies may return to Saginaw Bay in the early-21st century. Undoubtedly, watershed conservation and three decades of pollution abatement have set the stage for a recovery of burrowing mayflies in Saginaw Bay, and possibly in other areas of the Great Lakes.

Sesterhenn, T.M., C.R. Roswell, S.R. Stein, P. Klaver, E. Verhamme, S.A. POTHOVEN, and T.O. Hook. Modeling the implications of multiple hatching sites for larval dynamics in the resurgent Saginaw Bay walleye population. Journal of Great Lakes Research 40(Supplement 1):113-122 (DOI:10.1016/j.jglr.2013.09.022) (2014).

The early life environment experienced by most larval fish is largely dependent on a combination of hatch site and water currents. Until larvae are able to swim fast enough to overcome currents, they are largely passively transported and have limited control over ambient environmental conditions, including temperature and prey availability. These factors strongly influence growth and survival of larvae, with direct consequences for subsequent recruitment. Early life survival of Saginaw Bay walleye was formerly limited by alewife predation on larvae; but following the collapse of Lake Huron alewives, the walleye population has rebounded and recruitment success may now be influenced by other factors including spawning habitat. We sought to assess the implications of successful hatching at multiple locations in Saginaw Bay, using a hydrodynamics model, particle transport model, and an individual-based bioenergetics model in series. Model results were compared to locations of young larvae collected in Saginaw Bay during 2009–2010. Results suggest that larval growth is strongly influenced by hatch date, driven by seasonal variation in temperature between sites. Larvae hatched at any location could be transported extensively within inner Saginaw Bay before reaching a sufficient size to swim independently of currents, and retention within the productive inner bay varied among years and sites. Our results indicate multiple larval walleye origins in the field, augmenting the continued production from the Saginaw River system. Successful hatching at more locations would serve to buffer walleye recruitment variation through portfolio effects, supporting arguments for more emphasis on diverse spawning habitat management and restoration.

Soranno, P.A., K.S. Cheruvelil, E. Bissell, M.T. Bremigan, J.A. Downing, C.E. Fergus, C. Filstrup, E.N. Henry, N.R. Lottig, E.H. Stanley, C.A. STOW, P.-N. Tan, T. Wagner, and K.E. Webster. Cross-scale interactions: Quantifying multi-scaled cause-effect relationships in macrosystems. Frontiers in Ecology and the Environment 12(1):65-73 (DOI:10.1890/120366) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140006.pdf

Ecologists are increasingly discovering that ecological processes are made up of components that are multiscaled in space and time. Some of the most complex of these processes are cross-scale interactions (CSIs), which occur when components interact across scales. When undetected, such interactions may cause errors in extrapolation from one region to another. CSIs, particularly those that include a regional scaled component, have not been systematically investigated or even reported because of the challenges of acquiring data at sufficiently broad spatial extents. We present an approach for quantifying CSIs and apply it to a case study investigating one such interaction, between local and regional scaled land-use drivers of lake phosphorus. Ultimately, our approach for investigating CSIs can serve as a basis for efforts to understand a wide variety of multi-scaled problems such as climate change, land-use/land-cover change, and invasive species.

Staton, J.M., C.R. Roswell, D.G. Fielder, M.V. Thomas, S.A. POTHOVEN, and T.O. Hook. Condition and diet of yellow perch in Saginaw Bay, Lake Huron (1970-2011). Journal of Great Lakes Research 40(Supplement 1):139-147 (DOI:10.1016/j.jglr.2014.02.021) (2014).

In Saginaw Bay, Lake Huron, yellow perch (Perca flavescens) constitute an ecologically important component of the ecosystem and support both recreational and commercial fisheries. Over the past 40 years, Saginaw Bay has experienced multiple ecosystem-level changes (e.g., non-indigenous species introductions, reduced nutrient loading and variable temperatures). In turn, abundances and growth rates of yellow perch and their predators and prey have fluctuated. Recent changes to Saginaw Bay and Lake Huron foodwebs have potential to influence prey composition and subsequently, growth and condition for yellow perch; but a complete description of yellow perch diet composition across seasons has not been undertaken in recent years. We calculated mean relative weight (Wr), an index of condition, of age-1 and older yellow perch in Saginaw Bay annually for 1970–2011. We found high interannual variation in condition and documented low meanWr during 1978–1991.We developed regression models to explain this variation using phosphorus load, temperature, forage fish density, and yellow perch density as potential explanatory factors. Patterns of Wr were associated with changes in yellow perch densities, although interannual variation was not significantly associated with any of the available explanatory variables. Diet analysis of yellow perch collected in 2009 and 2010 demonstrated that age-1 and older yellow perch consumed a fundamentally different diet from a previous study (1986–1988), exhibiting a greater reliance on non-indigenous prey (e.g. Bythotrephes longimanus).

STOW, C.A. The news from Saginaw Bay: Where the mussels are strong, the walleye are good-looking, and all the phosphorus is above average. Journal of Great Lakes Research 40(Supplement 1):1-3 (DOI:10.1016/j.jglr.2014.04.006) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140023.pdf

Saginaw Bay has experienced some ups and downs since “Saginaw Michigan” was a hit song for Lefty Frizzell in 1964. By the 1960s the cumulative impact of years of anthropogenic stressors had taken its toll on the fishery and the ecosystem supporting it. The Bay was suffering from numerous water quality impairments including eutrophication, nuisance beach deposits (aka “muck”), and taste and odor problems at the drinking water intakes. A growing recognition of the loss of Saginaw Bay's ecosystem services, as they have since come to be called, helped catalyze a series of mitigation actions. In 1972 the US Congress passed the Clean Water Act which established defined water quality goals and began the process of regulating point and non-point discharges. The 1978 Great Lakes Water Quality Agreement set target total phosphorus loads for the Great Lakes, including a 440 MT/year target specific to Saginaw Bay, and a goal that the Bay should attain a mesotrophic state. Modeling efforts to support these developments established a template for subsequent efforts under the Clean Water Act, and the target loads became the prototype for what we now call TMDLs (Total Maximum Daily Loads). Initial assessments in the 1980s documented improvements in eutrophication-related water quality symptoms (Bierman et al., 1984); but then, attention was diverted to issues associated with toxic contaminants and nutrient monitoring lapsed.

STOW, C.A., Y.K. CHA, and S.S. Qian. A Bayesian hierarchical model to guide development and evaluation of substance objectives under the 2012 Great Lakes Water Quality Agreement. Journal of Great Lakes Research 40:7 pp. (DOI:10.1016/j.jglr.2014.07.005) (2014).

Under the 2012 Great Lakes Water Quality Agreement Canada and the United States are obliged to develop target concentrations for water quality constituents of particular concern. These “substance objectives” are closely analogous to numerical criteria under the US Clean Water Act. To develop effective substance objectives, it is important to consider how compliance with these objectives will be evaluated. Total phosphorus concentrations, for example, vary temporally and spatially, thus sample-based statistics will always be uncertain measures of the “true” underlying population characteristic. Using data from Saginaw Bay in Lake Huron, we develop a Bayesian hierarchical model that can be used to evaluate compliance with target concentrations on a temporally and spatially explicit basis. The “confidence of compliance” with targets can be assessed from the variance of the model parameter posterior distributions. This approach allows data to be grouped to represent spatial and temporal domains of particular interest, such as spring mean conditions in a certain area, and facilitates “partial pooling” of information so that regions with sparse data and high uncertainty can “borrow information” from more data-rich areas.

STOW, C.A., J. DYBLE, D.R. Kashian, T.H. JOHENGEN, K.P. Winslow, S.D. Peacor, S.N. Francoeur, A.M. BURTNER, D. PALLADINO, N. MOREHEAD, D.C. GOSSIAUX, Y.K. CHA, S.S. Qian, and D. Miller. Phosphorus targets and eutrophication objectives in Saginaw Bay: A 35-year assessment. Journal of Great Lakes Research 40(Supplement 1):4-10 (DOI:10.1016/j.jglr.2013.10.003) (2014).

An aggregated view of total phosphorus and chlorophyll a in Saginaw Bay indicates that concentrations of both constituents declined approximately in concert with declining total phosphorus (P) loads stabilizing by the late 1980s. A more spatially focused view reveals that total phosphorus declines outside of the Saginaw River plume, accompanied by more subtle chlorophyll a decreases. In contrast, soluble reactive phosphorus and ammonia have recently declined throughout the bay, while nitrate has remained relatively stable. Concentration data from nearshore transects do not exhibit large differences from open-water sample sites. The 440 tonne P/year target phosphorus load established in the 1978 amendments to the Great Lakes Water Quality Agreement has almost never been met, and total phosphorus concentrations regularly exceed the 15 μg/L concentration objective proposed in documentation supporting the 1978 amendments. Seasonal patterns in both total phosphorus and chlorophyll a are more pronounced in the most recent data, with peaks occurring in September–October. This apparently evolving seasonal pattern may result from seasonal changes in Saginaw River flow inputs, or seasonal variation in dreissenid mussel feeding and filtration rates. The adaptive management framework stipulated in the 2012 Great Lakes Water Quality Protocol should promote better monitoring of Saginaw Bay water quality into the future, with enhanced opportunities to better understand the factors that have maintained ongoing eutrophication symptoms.

STURTEVANT, R.A., J. LARSON, L. BERENT, M. McCARTHY, A. BOGDANOFF, A. Fusaro, and E.S. RUTHERFORD. An impact assessment of Great Lakes aquatic nonindigenous species. NOAA Technical Memorandum GLERL-161. NOAA, Great Lakes Environmental Research Laboratory, Ann Arbor, MI, 1022 pp. (2014). https://www.glerl.noaa.gov/pubs/tech_reports/glerl-161/tm-161.pdf

This purpose of this study is to provide a baseline assessment of realized, potential, and unknown impacts for established nonindigenous species in the Great Lakes. An organism impact assessment (OIA) tool was developed in order to analyze the extent of each species’ impact in a standardized manner. Following a thorough literature review, the OIA was used to rank the environmental impact, socio-economic impact, and beneficial effect of each species as high, moderate, low, or unknown. Importantly, this ranking system provides a method of identifying and comparing impacts across taxa and type of impact.

TANG, H., H.A. VANDERPLOEG, T.H. JOHENGEN, and J.R. LIEBIG. Quagga mussel (Dreissena rostriformis bugensis) selective feeding of phytoplankton in Saginaw Bay. Journal of Great Lakes Research 40(Supplement 1):83-94 (DOI:10.1016/j.jglr.2013.11.011) (2014).

Experiments from May to December measuring selective grazing and egestion of different phytoplankton taxa in natural Saginaw Bay (Lake Huron) seston by shallow-water morph quagga mussels (Dreissena bugensis rostriformis) showed that the mussels were highly selective filter feeders and that their net clearance rates on different species ranged widely, resulting in food consumption that was strongly driven by seasonal phytoplankton dynamics. Overall, net clearance rates by quagga mussels on the entire phytoplankton assemblage were similar to those observed for zebra mussels (Dreissena polymorpha) during the 1990s. Phytoplankton taxon, rather than size, was more important to food selection since quagga mussels cleared similar sized but different species of algae at very different rates. In contrast to many studies with zebra mussels, larger-sized algae such as Dinobryon divergens, Aulacoseira italica, Fragilaria crotonensis, and Anabaena were cleared at high rates and rejected at lower rates than many smaller species within the same broad taxonomic group. We suspect that these differences between dreissenid species do not stem from species differences but methodological factors and phytoplankton composition of systems studied. Small-sized diatoms, green algae with thick cell walls (Scenedesmus and Oocystis), and colonial cyanobacteria with gelatinous sheaths (Aphanocapsa, Chroococcus, and Microcystis) were cleared at low rates and rejected in high proportion in pseudofeces or feces during all seasons. We describe the likely mechanisms of pre- and post-ingestive behavior that explain these differences, which relate to phytoplankton size, morphology, cell wall characteristics, and chemical composition. Changes in the Great Lakes phytoplankton communities are consistent with mussel grazing preferences.

VANDERPLOEG, H.A., and J.R. Strickler. Video Clip 6: Behavior of zebra mussels exposed to Microcystis colonies from natural seston and laboratory cultures. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 757 pp. (2013).

Impacts of grazing by dreissenid mussels or zooplankton on the phytoplankton community are usually determined by traditional feeding experiments in which changes in phytoplankton concentration are measured over time in containers in which mussels or zooplankton are added (Vanderploeg et al. 2001, Dionisio Pires et al. 2005). These observations alone, however, may not provide much insight into feeding mechanisms and may have limited value for the development of general principles. Blooms of the cyanophyte Microcystis began to unexpectedly occur in Saginaw Bay, Lake Huron in 1994, which was 3 years after the bay was invaded by zebra mussels (Dreissena polymorpha). The blooms were unexpected because phosphorus (P), the usual driver of these blooms, had been decreasing in the bay for many years (Vanderploeg et al. 2001). Traditional feeding experiments with zebra mussels confined in small bottles with Saginaw Bay seston showed that Microcystis was removed at a lower rate than other phytoplankton (Lavrentyev et al. 1995); therefore, a link between the blooms and mussels was hypothesized. We had been using a sophisticated video system (Critter Cam; Bundy et al. 1998) to study feeding mechanisms in zooplankton, and when we conducted experiments to examine links between mussel feeding and the blooms of Microcystis, we felt further insights could be gained if we videotaped and observed mussel behavior at the same time we did the feeding experiments (Vanderploeg et al. 2001). Methods for making the videotapes were described in Vanderploeg et al. (2001, 2009). Like mussels in companion feeding experiments with Microcystis cultures, videotaped mussels received the same overnight acclimation to the suspensions made with these cultures (Vanderploeg et al. 2013). The videotapes were made in an aquarium that received no mixing from external sources; the water motion observed in the video sequences was caused by the filtering activity of the mussels. For point reference, total lengths of the mussels were 14-16 m.

VANDERPLOEG, H.A., A.E. Wilson, T.H. JOHENGEN, J. DYBLE BRESSIE, O. Sarnelle, J.R. LIEBIG, S.D. ROBINSON, and G.P. Horst. Role of selective grazing by dreissenid mussels in promoting toxic Microcystis blooms and other changes in phytoplankton composition in the Great Lakes. In Quagga and Zebra Mussels: Biology, Impacts, and Control, Second Edition. T.F. Nalepa, and D.W. Schlosser (Eds.). CRC Press, Boca Raton, FL, 509-523 pp. (2013).

We investigated the feeding response of zebra and quagga mussels to Microcystis aeruginosa strains from culture collection and from natural seston from Saginaw Bay (Lake Huron), western Lake Erie, and enclosures from Gull Lake, an inland lake in Michigan. These experiments were done lo evaluate the roles of strain identity, toxin concentration (microcystin), colony size, and environmental phosphorus concentrations as they affect ingestion or selective rejection of Microcystis in pseudofeces and potential Microcystis bloom promotion through the selective-rejection process.

WANG, J., X. BAI, J.F. BRATTON, B. Wu, and C.H. Greene. Accelerating Arctic summer sea ice decline driven by the dipole anomaly. Proceedings, 22nd IAHR International Symposium on Ice 2014, Singapore, August 11-15, 2014. International Association of Hydro-Environment Engineering and Research, 745-742 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20040056.pdf

Since the beginning of the satellite era (1979), record lows and highs of Arctic summer sea ice extent are found to be caused not only by dynamically-forced meridional wind advection, but also by the intensified ice/ocean albedo feedback, both of which are triggered by the Arctic atmospheric Dipole Anomaly (DA) pattern. This local, secondary mode of anomalous sea-level pressure (SLP) in the Arctic produces a strong anomalous meridional wind, which was responsible for driving more sea ice out of the Arctic Ocean from the western to the eastern Arctic into the northern Atlantic during the summers of 1995, 1999, 2002, 2005, and 2007-2012 during its positive phase. In September of 2012, a new record minimum in Arctic sea ice extent was caused by the combined dynamic wind forcing of +DA and thermodynamic forcing, accelerating melting, and advective loss of Arctic sea ice. Sea ice extent in August 2012 exceeded the previous record low of September 2007, occurring one month earlier in the summer season. In addition, high summer sea ice extent was associated with the –DA dynamic wind forcing and +AO’s divergence effect. A new feedback mechanism for the Pacific Arctic region’s atmosphere-ice-ocean climate system is hypothesized to be responsible for the accelerating decline of Arctic summer sea ice. In this hypothesis, the DA provides an intermittent, pulse-like forcing to the system, directly driving sea ice loss by wind forcing and enhancing northward transport of warm Pacific water, and indirectly driving its loss through thermal melting associated with positive ice/ocean albedo feedback and negative ice/cloud feedback. These processes contributed to a series of summer ice minima, particularly occurring since 1995, in which the +DA was associated with accelerating the Transpolar Drift Stream and sea ice export, and amplifying the melting of sea ice.

WANG, J., X. BAI, A. FUJISAKI-MANOME, H. HU, and D. BELETSKY. Great Lakes ice and climate research and forecast. Proceedings, 22nd IAHR International Symposium on Ice 2014, Singapore, August 11-15, 2014. International Association of Hydro-Environment Engineering and Research, 48-55 pp. (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140054.pdf

This paper describes recent progress made by a team from the NOAA Great Lakes Environmental Research Laboratory (GLERL) and the University of Michigan (UMich) Cooperative Institute for Limnology and Ecosystems Research (CILER). Over the past six years (since 2007), this team has studied Great Lakes ice and regional climate in response to global climate changes and how to transfer scientific research results into predictions of lake ice on the scales of several days to several months. The Great Lakes are located at the edge of the action centers of the North Atlantic Oscillation (NAO) and the nodal points of the action centers of the Pacific-North America (PNA) pattern. Aloft, there exists a strong highly-fluctuating annual south-north displacement of the westerly jet. Great Lakes ice responds linearly to the NAO and nonlinearly (quadratically) to the El Niño and South Oscillation (ENSO or Nino3.4). As a result, both NAO and ENSO have impacts on lake ice, but neither of them solely dominates the Great Lakes regional climate and lake ice cover. The combined effects of both NAO and ENSO on lake ice provide high predictability skills using statistical regression models. For the first time, fully-coupled Great Lakes Ice-circulation Models (GLIM) with both dynamics and thermodynamics have been developed at GLERL/CILER to simulate and investigate the lake ice variations on the synoptic, seasonal, interannual, and decadal time scales. The hindcast results were validated using in situ, airborne, and satellite measurements for various periods. The validated GLIM has been used since the 2010-2011 ice season to forecast Great Lakes ice cover concentration, thickness, velocity, and associated air-ice-sea variables for up to five days in advance.

WANG, J., H. Eicken, Y. Yu, X. BAI, J. Zhang, H. HU, D.R. Wang, M. Ikeda, K. Mizobata, and J. Overland. Chapter 4. Abrupt climate changes and emerging ice-ocean processes in the Pacific Arctic region and the Bering Sea. In The Pacific Arctic Region: Ecosystem Status and Trends in a Rapidly Changing Environment. J.M. Grebmeier and W. Maslowski (Eds.). Springer Science+Business Media, Dordrecht, 65-99 pp. (DOI:10.1007/978-94-017-8863-2_4) (2014).

The purpose of this chapter is to reveal several emerging physical ice-ocean processes associated with the unprecedented sea ice retreat in the Pacific Arctic region (PAR). These processes are closely interconnected under the scenario of diminishing sea ice, resulting in many detectable changes from physical environment to ecosystems. Some of these changes are unprecedented and have drawn the attention of both scientific and societal communities. More importantly, some mechanisms responsible for the diminishing sea ice cannot be explained by the leading Arctic Oscillation (AO), which has been used to interpret most of the changes in the Arctic for the last several decades. The new challenging questions are: (1) What is the major forcing? (2) Is the AO, the DA, or their combination, contributing to the sea ice minima in recent years? How do we use models to investigate the recent changes in the PAR. Is the heat transport through the Bering Strait associated with the DA? What processes accelerate sea ice melting in the PAR?

WANG, J., K. Mizobata, X. BAI, H. HU, M. Jin, Y. Yu, M. Ikeda, W. Johnson, W. Perie, and A. FUJISAKI. A modeling study of coastal circulation and landfast ice in the nearshore Beaufort and Chukchi Seas using CIOM. Journal of Geophysical Research: Oceans 19(6):3285-3312 (DOI:10.1002/2013JC009258) (2014).

This study investigates sea ice and ocean circulation using a 3-D, 3.8 km CIOM (Coupled Ice- Ocean Model) under daily atmospheric forcing for the period 1990–2008. The CIOM was validated using both in situ observations and satellite measurements. The CIOM successfully reproduces some observed dynamical processes in the region, including the Bering-inflow-originated coastal current that splits into three branches: Alaska Coastal Water (ACW), Central Channel branch, and Herald Valley branch. In addition, the Beaufort Slope Current (BSC), the Beaufort Gyre, the East Siberian Current (ESC), mesoscale eddies, and seasonal landfast ice are well simulated. The CIOM also reproduces reasonable interannual variability in sea ice, such as landfast ice, and anomalous open water (less sea ice) during the positive Dipole Anomaly (DA) years, vice versa during the negative DA years. Sensitivity experiments were conducted with regard to the impacts of the Bering Strait inflow (heat transport), onshore wind stress, and sea ice advection on sea ice change, in particular on the landfast ice. It is found that coastal landfast ice is controlled by the following processes: wind forcing, Bering Strait inflow, and sea ice dynamics.

Wittmann, M., R. Cooke, J. Rothlisberger, E.S. RUTHERFORD, H. ZHANG, D.M. MASON, and D. Lodge. Use of structured expert judgment to forecast invasions by Bighead and Silver carp in Lake Erie. Conservation Biology:10 pp. (DOI:10.1111/cobi.12369) (2014). https://www.glerl.noaa.gov/pubs/fulltext/2014/20140041.pdf

Identifying which nonindigenous species will become invasive, and forecasting their damages 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 expert combination where equal weights are assigned to experts, SEJ weights each expert on the basis of his or her statistical accuracy and informativeness, and can be used when the information required for decision-making is sparse or debated. We used SEJ to forecast impacts of nonindigenous Asian carp in Lake Erie, where it is believed not to be established. Experts quantified Asian carp biomass, production, consumption and the impact to four fish species if Asian carp were to become established. According to experts, Asian carp can establish in Lake Erie with 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 is estimated by experts to be small, with little uncertainty. Impacts in tributaries and to 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.

Wohlleben, T., A. Tivy, W. Stroeve, W. Meier, F. Fetterer, J. WANG, and R.A. ASSEL. Computing and representing sea ice trends: Toward a community consensus. EOS 94(40):352 (DOI:10.1002/2013EO400006) (2013). https://www.glerl.noaa.gov/pubs/fulltext/2013/20130034.pdf

Estimates of the recent decline in Arctic Ocean summer sea ice extent can vary due to differences in sea ice data sources, in the number of years used to compute the trend, and in the start and end years used in the trend computation. Compounding such differences, estimates of the relative decline in sea ice cover (given in percent change per decade) can further vary due to the choice of reference value (the initial point of the trend line, a climatological baseline, etc.). Further adding to the confusion, very often when relative trends are reported in research papers, the reference values used are not specified or made clear. This can lead to confusion when trend studies are cited in the press and public reports.

Wu, J., A.D. GRONEWOLD, R. Rodriquez, J.R. Stewart, and M.D. Sobsey. Integrating quantitative PCR and Bayesian statistics in quantifying human adenoviruses in small volumes of source water. Science of the Total Environment 470-471:255-262 (DOI:10.1016/j.scitotenv.2013.09.026) (2014).

Rapid quantification of viral pathogens in drinking and recreational water can help reduce waterborne disease risks. For this purpose, samples in small volume (e.g. 1 L) are favored because of the convenience of collection, transportation and processing. However, the results of viral analysis are often subject to uncertainty. To overcome this limitation, we propose an approach that integrates Bayesian statistics, efficient concentration methods, and quantitative PCR (qPCR) to quantify viral pathogens in water. Using this approach, we quantified human adenoviruses (HAdVs) in eighteen samples of source water collected from six drinking water treatment plants. HAdVs were found in seven samples. In the other eleven samples, HAdVs were not detected by qPCR, but might have existed based on Bayesian inference. Our integrated approach that quantifies uncertainty provides a better understanding than conventional assessments of potential risks to public health, particularly in cases when pathogens may present a threat but cannot be detected by traditional methods.

ZHANG, H., D.M. MASON, C.A. STOW, A.T. ADAMACK, S.B. Brandt, X. Zhang, D.G. Kimmel, M.R. Roman, W.C. Boicourt, and S.A. Ludsin. Effects of hypoxia on habitat quality of pelagic planktivorous fishes in the northern Gulf of Mexico. Marine Ecology Progress Series 505:209-226 (DOI:10.3354/meps10768) (2014).

To evaluate the impact of hypoxia (<2 mg O2 l−1) on habitat quality of pelagic prey fishes in the northern Gulf of Mexico, we used a spatially explicit, bioenergetics-based growth rate potential (GRP) model to develop indices of habitat quality. Our focus was on the pelagic bay anchovy Anchoa mitchilli and Gulf menhaden Brevoortia patronus. Positive GRP was considered high-quality habitat (HQH) and negative GRP was considered low-quality habitat (LQH). Models used water temperature, dissolved oxygen (DO), zooplankton biomass, and phytoplankton concentration collected during the peak periods of hypoxia in 2003, 2004, and 2006 to estimate fish GRP. Results showed that hypoxic areas were always LQH. However, with respect to the entire water column, hypoxia had only a minor impact on overall habitat quality, with habitat quality being driven primarily by prey availability followed by water temperature. These results are in contrast to other ecosystems, such as the Chesapeake Bay, where hypoxia affects a larger fraction of the water column than in the Gulf of Mexico and has a significant impact on overall habitat quality. Differences in the effect of hypoxia on habitat quality between these 2 ecosystems suggest that the vertical extent of hypoxia relative to water column depth (i.e. hypoxic volume) is a fundamental consideration when evaluating the impacts of hypoxia on pelagic fish production.

 

Ordering GLERL Publications:

To order a copy of GLERL publications not available for downloading at this site, please contact:

Nicole Rice Information Services

NOAA Great Lakes Environmental Research Laboratory
4840 S. State Rd.
Ann Arbor, MI 48108 US
734-741-2247
734-741-2055 (FAX