Verification and Application of a Bio-optical Algorithm for Lake Michigan Using SeaWiFS: a 7-year Inter-annual Analysis

In this paper we utilize 7 years of SeaWiFS satellite data to obtain seasonal and interannual time histories of the major water color-producing agents (CPAs), phytoplankton chlorophyll (chl), dissolved organic carbon (doc), and suspended minerals (sm) for Lake Michigan. We first present validation of the Great Lakes specific algorithm followed by correlations of the CPAs with coincident environmental observations. Special attention is paid to the satellite observations of the extensive episodic event of sediment resuspension and calcium carbonate precipitation out of the water. We then compare the obtained time history of the CPA's spatial and temporal distributions throughout the lake to environmental observations such as air and water temperature, wind speed and direction, significant wave height, atmospheric precipitation, river runoff, and cloud and lake ice cover. Variability of the onset, duration, and spatial extent of both episodic events and seasonal phenomena are documented from the SeaWiFS time series data, and high correlations with relevant environmental driving factors are established. The relationships between the CPAs retrieved from satellite data and environmental observations are then used to speculate on the future of Lake Michigan under a set of climate change scenarios.     

An assessment of water quality in two Great Lakes connecting channels

Compared to the Great Lakes, their connecting channels are relatively understudied and infrequently assessed. To address this gap, we conducted a spatially-explicit water quality assessment of two connecting channels, the St. Marys River and the Lake Huron-Lake Erie Corridor (HEC) in 2014–2016. We compared the condition of the channels to each other and to the up- and downriver Great Lakes with data from an assessment of the Great Lakes nearshore. In the absence of channel-specific thresholds, we assessed the condition of the area of each channel as good, fair, or poor by applying the most protective water quality thresholds for the downriver lake. Condition of the St. Marys River was rated mostly fair for total phosphorus (TP, 56% of the area) and mostly good (61% of the area) for chlorophyll a. Area-weighted mean concentrations of these parameters were intermediate to Lake Superior and Lake Huron. Unlike Lake Superior and Lake Huron, a large proportion (97%) of the area of the St. Marys River was in poor condition for water clarity based on Secchi depth. Area-weighted mean concentrations of TP and chlorophyll a in the HEC were more like Lake Huron than Lake Erie. For these indicators, most of the area of the HEC was rated good (81% and 86%, respectively). Interpretation of assessment results is complicated by variation in thresholds among and within lakes. Appropriate thresholds should align with assessment objectives and in the case of connecting channels be at least as protective as thresholds for the downriver lake.     

A band-ratio algorithm for retrieving open-lake chlorophyll values from satellite observations of the Great Lakes

The U.S. Environmental Protection Agency's Great Lakes National Program Office (GLNPO) has collected water quality data from the five Great Lakes annually since 1993. We used the GLNPO observations made since 2002 along with coincident measurements made by the Sea-viewing Wide Field-of-View Sensor (SeaWiFS) and the Moderate-resolution Imaging Spectroradiometer (MODIS) to develop a new band-ratio algorithm for estimating chlorophyll concentrations in the Great Lakes from satellite observations. The new algorithm is based on a third-order polynomial model using the same maximum band ratios employed in the standard NASA algorithms (OC4 for SeaWiFS and OC3M for MODIS). The sensor-specific coefficients for the new algorithm were obtained by fitting the relationship to several hundred matched field and satellite observations. Although there are some seasonal variations in some lakes, the relationship between the observed chlorophyll values and those modeled using the new coefficients is fairly stable from lake to lake and across years. The accuracy of the satellite chlorophyll estimates derived from the new algorithm was improved substantially relative both to the standard NASA retrievals and to previously published algorithms tuned to individual lakes. Monte-Carlo fits to randomly selected subsets of the observations allowed us to estimate the uncertainty associated with the retrievals purely as a function of the satellite data. Our results provide, for the first time, a single simple band ratio method for retrieving chlorophyll concentrations in the offshore “open” waters of the Great Lakes from satellite observations.     

Laurentian Great Lakes Hydrology and Lake Levels under the Transposed 1993 Mississippi River Flood Climate

The Laurentian Great Lakes are North America's largest water resource, and include six large water bodies (Lakes Superior, Michigan, Huron, Erie, Ontario, and Georgian Bay), Lake St. Clair, and their connecting channels. Because of the relatively small historical variability in system lake levels, there is a need for realistic climate scenarios to develop and test sensitivity and resilience of the system to extreme high lake levels. This is particularly important during the present high lake level regime that has been in place since the late 1960s. In this analysis, we use the unique climate conditions which resulted in the 1993 Mississippi River flooding as an analog to test the sensitivity of Great Lakes hydrology and water levels to a rare but actual climate event. The climate over the Upper Mississippi River basin was computationally shifted, corresponding to a conceptual shift of the Great Lakes basin 10̊ west and 2̊ south. We applied a system of hydrological models to the daily meteorological time series and determined daily runoff, lake evaporation, and net basin water supplies. The accumulated net basin supplies from May through October 1993 for the 1993 Mississippi River flooding scenario ranged from a 1% decrease for Lake Superior to a large increase for Lake Erie. Water levels for each lake were determined from a hydro-logic routing model of the system. Lakes Michigan, Huron, and Erie were most affected. The simulated rise in Lakes Michigan and Huron water levels far exceeded the historically recorded rise with both lakes either approaching or setting record high levels. This scenario demonstrates that an independent anomalous event, beginning with normal lake levels, could result in record high water levels within a 6- to 9-month period. This has not been demonstrated in the historical record or by other simulation studies.     

Trace Element Concentrations in Near-Surface Waters of the Great Lakes and Methods of Collection,Storage, and Analysis

From 1980 through 1985, waters of the Great Lakes were sequentially sampled for dissolved, paniculate, and total trace elements. Major sampling occurred in 1980 for Lake Huron, in 1981 for Lakes Erie and Michigan, in 1983 for Lake Superior, and in 1985 for Lake Ontario. Great care was taken during collection, storage, and analysis to prevent sample contamination and to document any contamination occurring. Trace elements measured by atomic absorption techniques were silver, aluminum, arsenic, boron, barium, beryllium, bismuth, cadmium, cobalt, chromium, copper, iron, mercury, lithium, manganese, molybdenum, nickel, lead, antimony, selenium, tin, strontium, vanadium, and zinc. All results were field and laboratory blank corrected. Excluding aluminum, barium, iron, and strontium, concentrations of trace elements in most of the Great Lakes were a few ppb or less, with many elements being below one ppb. Element concentrations were highest in Lakes Erie and Michigan and lowest in Lakes Huron and Superior. All five Great Lakes had more than 50% of their total iron, aluminum, and manganese associated with paniculate matter.     

Results from the U.S. EPA's Biological Open Water Surveillance Program of the Laurentian Great Lakes: I. Introduction and Phytoplankton Results

The Great Lakes National Program Office of the U.S. EPA has been conducting biological monitoring of the Laurentian Great Lakes since 1983. This paper presents synoptic survey data of phytoplankton communities from all five lakes. These communities were highly diverse, each lake typi-cally supporting over 100 species during both the spring and summer surveys. Much of that diversity was contributed by diatoms, which dominated the plankton of all lakes except Lake Superior in the spring. Summer communities shifted away from diatoms, toward chrysophytes in the upper lakes and chloro-phytes in the lower lakes. Ordination analyses indicated the close similarity of communities in the upper lakes, in particular Lakes Huron and Michigan, and a diverse range of communities in Lake Erie. Floristically, Lake Ontario was fundamentally different from all other lakes.     

Concentrations of Polychlorinated Biphenyls in the Water Column of the Laurentian Great Lakes: Spring 1993

PCB concentrations were measured in the Laurentian Great Lakes of North America in the Spring of 1993. Quality control criteria were met for 97% of the dissolved phase samples and 90% of the particulate phase samples. Data are reported as total PCBs and as homolog groups for both of the operationally-defined phases at a 0.7 μM cutoff. Detection limit estimates were 21 pg/L for dissolved phase and 23 pg/L for particulate phase total PCBs, based on a sample volume of 190 L. Total PCB concentrations ranged from 100 pg/L in Lake Superior to 1.6 ng/L in the western basin of Lake Erie. Lakes Michigan and Ontario had the highest lakewide average total PCB concentrations, followed by Lake Huron, and then Lake Superior. Lakewide averaging across the three basins of Lake Erie was too variable to be useful.     

Are PCB Levels in Fish from the Canadian Great Lakes Still Declining?

Long- and short-term levels and trends of polychlorinated biphenyls (PCBs) in lake trout (Salvelinus namaycush) and walleye (Sander vitreus) from the Canadian waters of the Great Lakes are examined using the bootstrap resampling method in light of the Great Lakes Strategy 2002 (GLS-2002) objective of decrease in concentrations by 25% during 2000–2007. This objective has been set as an indicator of progress toward the long-term goal of all Great Lakes fish being safe to eat without restriction. Lake Superior lake trout and walleye PCB concentrations were almost unchanged between 1990-2006, and the bootstrap analysis suggests that the probability of achieving the GLS-2002 objective is negligible (< 2%). The PCB levels in Lake Huron lake trout and walleye are decreasing; the declines between 2000–2007 are estimated to be 25–35% and 5–30%, respectively. In contrast, Lake Erie walleye concentrations will likely increase by 25–50% between 2000–2007. For Lake Ontario lake trout, achieving the 25% reduction target seems highly probable with a likely decrease of 45–55%; for Lake Ontario walleye, the probability of achieving such a reduction is only 8% with an expected change of −13 to +15%. Although the targeted reduction may not be achieved for walleye from Lakes Superior, Huron, and Ontario, their best projected 2007 PCB levels are below the unlimited fish consumption guideline of 105 ng/g wet weight used by the Ontario Ministry of the Environment. In contrast, although there are high probabilities of achieving the goal for lake trout from Lakes Huron and Ontario, their best projected 2007 PCB levels (160 and 370 ng/g ww, respectively) will continue to result in consumption restrictions. Lake Superior lake trout concentrations may remain unchanged at the current elevated level of 160 ng/g ww. For Lake Erie fish, the projected 2007 concentrations and the increasing trends are both worrisome. Additional measurements beyond 2007 are necessary to confirm these estimates because of the observed periodic oscillations in the concentrations.     

Using satellite observations to assess the spatial representativeness of the GLNPO water quality monitoring program

The U.S. EPA's Great Lakes National Program Office (GLNPO) annual water quality survey (WQS) collects data at a relatively small number of stations in each lake. The survey was designed to measure conditions in the open-water regions of the lakes where an assumption of spatial homogeneity was thought likely to be met and the measured variables could be characterized by simple statistics. Here we use satellite observations to assess how well statistics based on samples collected in the GLNPO sampling network represent the lake-wide values of two variables, surface chlorophyll concentration and Secchi depth. We find strong linear relationships between the mean values calculated from the samples and the corresponding averages based on the subsets of the full satellite images. Although overall the means of the values from the sample locations agree well with means calculated from most of the non-coastal regions of the lakes, in terms of water depth, the GLNPO station averages best represent the regions of Lake Huron deeper than 30?m, of Lakes Michigan and Superior deeper than 90?m, and of Lake Ontario deeper than 60?m. When the lake regions are defined by distance offshore rather than by depth, the GLNPO station chlorophyll means in Lakes Huron, Ontario, and Superior are closest to the means for the area of the lakes >10?km offshore. In Lake Michigan the closest correspondence is with the >20?km offshore region. On a whole-lake basis in Lake Erie the GLNPO station chlorophyll averages are closest to the average calculated from the entire lake.     

The benthic community of the Laurentian Great Lakes: Analysis of spatial gradients and temporal trends from 1998 to 2014

We used the results of seventeen years of Great Lakes benthic monitoring conducted by the U.S. EPA's Great Lakes National Program Office to describe the spatial and temporal patterns of benthic communities, assess their status, trends, and main drivers, and to infer the potential impact of these community changes on ecosystem functioning. Benthic abundance and diversity were higher at shallow (<70?m in depth) stations with chlorophyll concentrations above 3?μg/L than at deeper sites (<1?μg/L). We infer that lake productivity, measured by chlorophyll was likely the major driver of benthic abundance and diversity across lakes. Consequently, benthic diversity and abundance were the highest in the most productive Lake Erie, followed by lakes Ontario, Michigan, Huron, and Superior. Multivariate analysis distinguished three major communities shared among lakes (littoral, sublittoral, and profundal) that differed in species composition and abundance, functional group diversity, and tolerance to organic pollution. Analysis of temporal trends revealed that the largest changes occurred in profundal communities, apparent in significant shifts in dominant taxa across all lakes except Lake Superior. In lakes Michigan, Huron, and Ontario, the former dominant Diporeia was replaced with Dreissena and Oligochaeta. Profundal species, with the exception of dreissenids, became less abundant, and their depth distribution has shifted. In contrast, density and diversity of native littoral and sublittoral communities increased. The invasion of dreissenids was among the most important drivers of changes in benthic communities. Continued monitoring is critical for tracking unprecedented changes occurring in the Great Lakes ecosystem.     

A cross-lake comparison of crustacean zooplankton communities in the Laurentian Great Lakes, 1997–2016

Spring and summer open-water crustacean zooplankton communities were examined across all five Laurentian Great Lakes from 1997 to 2016. Spring communities were dominated by calanoid (lakes Superior, Huron and Michigan) or cyclopoid (lakes Erie and Ontario) copepods. Volumetric biomass of summer communities increased along an assumed trophic gradient (Superior, Huron, Michigan, Ontario; eastern, central and western Erie), as did dominance by cyclopoids and cladocerans. Over the time series of the study, summer communities in lakes Michigan, Huron and Ontario shifted towards greater dominance by calanoids and greater similarity with Lake Superior. Trajectories of changes were different; however, reductions in cladocerans accounted for most of the change in lakes Michigan and Huron while reductions in cyclopoids and increases in Leptodiaptomus sicilis were behind the changes in Lake Ontario. Shifts in the predatory cladoceran community in Lake Ontario from Cercopagis pengoi to occasional dominance by Bythotrephes longimanus, a species much more vulnerable to planktivory, as well as the appearance of Daphnia mendotae in a daphnid community previously consisting almost exclusively of the smaller Daphnia retrocurva, suggest impacts of reduced vertebrate predation. In contrast, strong correlations between cladocerans and chlorophyll in lakes Michigan and Huron point to the possible importance of bottom-up forces in those lakes. Large interannual shifts in cladoceran community structure in the central and eastern basins of Lake Erie suggest intense but variable vertebrate predation pressure. The zooplankton communities of lakes Huron, Michigan and Ontario may be approaching a historic community structure represented by Lake Superior.     

Mean Circulation in the Great Lakes

In this paper new maps are presented of mean circulation in the Great Lakes, employing long-term current observations from about 100 Great Lakes moorings during the 1960s to 1980s. Knowledge of the mean circulation in the Great Lakes is important for ecological and management issues because it provides an indication of transport pathways of nutrients and contaminants on longer time scales. Based on the availability of data, summer circulation patterns in all of the Great Lakes, winter circulation patterns in all of the Great Lakes except Lake Superior, and annual circulation patterns in Lakes Erie, Michigan, and Ontario were derived. Winter currents are generally stronger than summer currents, and, therefore, annual circulation closely resembles winter circulation. Circulation patterns tend to be cyclonic (counterclockwise) in the larger lakes (Lake Huron, Lake Michigan, and Lake Superior) with increased cyclonic circulation in winter. In the smaller lakes (Lake Erie and Lake Ontario), winter circulation is characterized by a two-gyre circulation pattern. Summer circulation in the smaller lakes is different; predominantly cyclonic in Lake Ontario and anticyclonic in Lake Erie.     

Great Lakes total phosphorus revisited: 1. Loading analysis and update (1994–2008)

Phosphorus load estimates have been updated for all of the Great Lakes with an emphasis on lakes Superior, Michigan, Huron and Ontario for 1994–2008. Lake Erie phosphorus loads have been kept current with previous work and for completeness are reported here. A combination of modeling and data analysis is employed to evaluate whether target loads established by the Great Lakes Water Quality Agreement (GLWQA, 1978, Annex 3) have been and are currently being met. Data from federal, state, and provincial agencies were assembled and processed to yield annual estimates for all lakes and sources. A mass-balance model was used to check the consistency of loads and to estimate interlake transport. The analysis suggests that the GLWQA target loads have been consistently met for the main bodies of lakes Superior, Michigan and Huron. However, exceedances still persist for Saginaw Bay. For lakes Erie and Ontario, loadings are currently estimated to be at or just under the target (with some notable exceptions). Because interannual variability is high, the target loads have not been met consistently for the lower Great Lakes. The analysis also indicates that, because of decreasing TP concentrations in the lakes, interlake transport of TP has declined significantly since the mid-1970s. Thus, it is important that these changes be included in future assessments of compliance with TP load targets. Finally, detailed tables of the yearly (1994–2008) estimates are provided, as well as annual summaries by lake tributary basin (in Supplementary Information).     

Historical changes and current status of crayfish diversity and distribution in the Laurentian Great Lakes

Despite increasing recognition of the importance of invertebrates, and specifically crayfish, to nearshore food webs in the Laurentian Great Lakes, past and present ecological studies in the Great Lakes have predominantly focused on fishes. Using data from many sources, we provide a summary of crayfish diversity and distribution throughout the Great Lakes from 1882 to 2008 for 1456 locations where crayfish have been surveyed. Sampling effort was greatest in Lake Michigan, followed by lakes Huron, Erie, Superior, and Ontario. A total of 13 crayfish species occur in the lakes, with Lake Erie having the greatest diversity (n = 11) and Lake Superior having the least (n = 5). Five crayfish species are non-native to one or more lakes. Because Orconectes rusticus was the most widely distributed non-native species and is associated with known negative impacts, we assessed its spread throughout the Great Lakes. Although O. rusticus has been found for over 100 years in Lake Erie, its spread there has been relatively slow compared to that in lakes Michigan and Huron, where it has spread most rapidly since the 1990s and 2000, respectively. O. rusticus has been found in both lakes Superior and Ontario for 22 and 37 years, respectively, and has expanded little in either lake. Our broad spatial and temporal assessment of crayfish diversity and distribution provides a baseline for future nearshore ecological studies, and for future management efforts to restore native crayfish and limit non-native introductions and their impact on food web interactions.     

Warmer and Drier Climates that Make Terminal Great Lakes

A recent empirical model of glacial-isostatic uplift showed that the Huron and Michigan lake level fell tens of meters below the lowest possible outlet about 7,900 ¹⁴C years BP when the upper Great Lakes became dependent for water supply on precipitation alone, as at present. The upper Great Lakes thus appear to have been impacted by severe dry climate that may have also affected the lower Great Lakes. While continuing paleoclimate studies are corroborating and quantifying this impacting climate and other evidence of terminal lakes, the Great Lakes Environmental Research Laboratory applied their Advanced Hydrologic Prediction System, modified to use dynamic lake areas, to explore the deviations from present temperatures and precipitation that would force the Great Lakes to become terminal (closed), i.e., for water levels to fall below outlet sills. We modeled the present lakes with pre-development natural outlet and water flow conditions, but considered the upper and lower Great Lakes separately with no river connection, as in the early Holocene basin configuration. By using systematic shifts in precipitation, temperature, and humidity relative to the present base climate, we identified candidate climates that result in terminal lakes. The lakes would close in the order: Erie, Superior, Michigan-Huron, and Ontario for increasingly drier and warmer climates. For a temperature rise of T°C and a precipitation drop of P% relative to the present base climate, conditions for complete lake closure range from 4.7T + P > 51 for Erie to 3.5T + P > 71 for Ontario.     

Trends in Mysis diluviana abundance in the Great Lakes, 2006–2016

With the large Diporeia declines in lakes Michigan, Huron, and Ontario, there is concern that a similar decline of Mysis diluviana related to oligotrophication and increased fish predation may occur. Mysis density and biomass were assessed from 2006 to 2016 using samples collected by the Great Lakes National Program Office's biomonitoring program in April and August in all five Great Lakes. Summer densities and biomasses were generally greater than spring values and both increased with bottom depth. There were no significant time trends during these 10–11 years in lakes Ontario, Michigan, or Huron, but there was a significant increase in Lake Superior. Density and biomass were highest in lakes Ontario and Superior, somewhat lower in Lake Michigan, and substantially lower in Lake Huron. A few Mysis were collected in eastern Lake Erie, indicating a small population in the deep basin of that lake. On average, mysids contributed 12–18% (spring-summer, Michigan), 18–14% (spring-summer, Superior), 30–13% (spring-summer, Ontario), and 3% (Huron) of the total open-water crustacean biomass. Size distributions consisted of two peaks, indicating a 2-year life cycle in all four of the deep lakes. Mysis were larger in Lake Ontario than in lakes Michigan, Superior, and Huron. Comparisons with available historic data indicated that mysid densities were higher in the 1960s–1990s (5 times higher in Huron, 2 times higher in Ontario, and around 40% higher in Michigan and Superior) than in 2006–2016.     

A model for determining satellite-derived primary productivity estimates for Lake Michigan

A new MODIS based satellite algorithm to estimate primary production (PP) has been generated and evaluated for Lake Michigan. The Great Lakes Primary Productivity Model (GLPPM) is based on previous models that required extensive in situ data but it can utilize remotely sensed observations as input for some model variables and therefore allows greater spatial resolution for primary productivity estimates. The Color Producing Agent Algorithm (CPA-A) is utilized to obtain robust chlorophyll a values and the NASA KD2M approach is used to obtain the diffuse attenuation coefficient (K_d). Only incident PAR and carbon fixation rates are additionally needed to generate the primary productivity estimate. Comparisons of the satellite derived PP estimates from single monthly images to average monthly field measurements made by NOAA/GLERL found good agreement between estimates. Satellite derived PP estimates were used to estimate a preliminary Lake Michigan annual primary production of 8.5 Tg C/year. The new algorithm can be easily adapted to work on all the Great Lakes and therefore can be used to generate time series dating back to late 1997 (launch of SeaWiFs). These time series can contribute to improved assessment of Great Lakes primary productivity changes as a result of biological events, such as Dreissenid mussel invasions, climatic change and anthropogenic forcing.     

Toxaphene trends in the Great Lakes fish

As part of the U.S. Great Lakes Fish Monitoring and Surveillance Program (GLFMSP), more than 300 lake trout (Salvelinus namaycush) and walleye (Stizostedion vitreum vitreum) collected from the Laurentian Great Lakes each year from 2004 to 2009, have been analyzed for total toxaphene and eight selected congeners. The analytical results show fish toxaphene concentrations are quite different among lakes. Between 2004 and 2009, Lake Superior lake trout had the highest concentration (119 to 482 ng/g) and Lake Erie walleye had the lowest concentration (18 to 47 ng/g). Combining these results with the historical total toxaphene data (1977–2003), temporal changes were examined for each lake. Because of different analytical methods used in the previous studies, the historical data were adjusted using a factor of 0.56 based on a previous inter-method comparison in our laboratory. Trend analysis using an exponential decay regression showed that toxaphene in Great Lakes fish exhibited a significant decrease in all of the lakes with t_1/2 (confidence interval) of 0.9 (0.8–1.1) years for Lake Erie walleye, 3.8 (3.5–4.1) years for Lake Huron lake trout, 5.6 (5.1–6.1) years for Lake Michigan lake trout, 7.5 (6.7–8.4) years for Lake Ontario lake trout and 10.1 (8.2–13.2) years for Lake Superior lake trout. Parlars 26, 50 and 62 were the dominant toxaphene congeners accounting for 0.53% to 41.7% of the total toxaphene concentration. Concentrations of these congeners generally also decreased over time.     

Stability Effects on Great Lakes Evaporation

The effects of atmospheric stability on Great Lakes evaporation computed by the modified mass transfer method have been evaluated by analysis of stability effects on the variable mass transfer coefficient, land to lake data adjustments, and ice-cover reduction of evaporation during winter. The Great Lakes which produce extreme results, Lakes Superior and Erie, and a much smaller water body within the Great Lakes chain, Lake St. Clair, were studied. Comparison of these evaporation estimates with previous studies, which excluded variable stability effects, shows that the previous studies of Lake Superior produced generally similar total annual water loss from the lake, but significantly overestimated both the seasonal high evaporation and the condensation rates. These tended to balance each other. The atmospheric conditions over Lakes Erie and St. Clair do not become as strongly stable and they normally do not exhibit large condensation. Previous evaporation studies for these lakes indicate generally higher evaporation rates, with significant overestimation of the total annual water losses (25%).     

Alkalinity,pH, and pCO2 in the Laurentian Great Lakes: An initial view of seasonal and inter-annual trends

Ongoing human perturbations to the global inorganic carbon cycle can cause various changes in the pH and alkalinity of aquatic systems. Here seasonal and inter-annual trends in these variables were investigated in the five Laurentian Great Lakes using data from the U.S. EPA GLENDA database. These observations, along with temperature, were also used to predict the partial pressure of carbon dioxide in surface water (pCO₂). There are strong seasonal differences in pH in all five lakes, with higher pH levels in summer than in spring. All lakes show significantly higher pCO₂ values in spring than in summer. Michigan and Ontario show higher alkalinity values in spring; Huron shows lower spring values. Inter-annually, open-lake pH shows the highest values in all lakes around 2010, the time frame of lowest lake water levels, though only lakes Superior and Erie show statistically significant inflection points at that time. Inter-annual alkalinity trends differ considerably from those of pH. Superior’s alkalinity increases until ～2008 and then begins dropping; Ontario’s alkalinity decreases until ～2004 and then begins increasing, with the decrease coinciding with the introduction and establishment of Dreissenid mussels. The other lakes show much less clear inter-annual alkalinity trends. For pCO₂, inter-annual trends in each lake show either increases from 1992 to 2019 (for Superior, Michigan, and Huron) or shifts from slightly decreasing values to increasing values for the other lakes. The timing of this shift is from 2008 to 2010.