Abstracts
Abstract
In the Great Lakes region, the vertical motion of crustal rebound since the last glaciation has decelerated with time, and is described by exponential decay constrained by observed warping of strandlines of former lakes. A composite isostatic response surface relative to an area southwest of Lake Michigan beyond the limit of the last glacial maximum was prepared for the complete Great Lakes watershed at 10.6 ka BP (12.6 cal ka BP). Uplift of sites computed using values from the response surface facilitated the transformation of a digital elevation model of the present Great Lakes basins to represent the paleogeography of the watershed at selected times. Similarly, the original elevations of radiocarbon-dated geomorphic and stratigraphic indicators of former lake levels were reconstructed and plotted against age to define lake level history. A comparison with the independently computed basin outlet paleo-elevations reveals a phase of severely reduced water levels and hydrologically-closed lakes below overflow outlets between 7.9 and 7.0 ka BP (8.7 and 7.8 cal ka BP) in the Huron-Michigan basin. Severe evaporative draw-down is postulated to result from the early Holocene dry climate when inflows of meltwater from the upstream Agassiz basin began to bypass the upper Great Lakes basin.
Résumé
Dans la région des Grands Lacs, le soulèvement isostatique lié à la dernière glaciation a ralenti selon une courbe de décroissance exponentielle établie à partir du gauchissement observé dans les anciennes lignes de rivage lacustres. Une surface de référence composite de la réponse isostatique datant de 10,6 ka BP (12,6 cal ka BP) a été préparée pour l’ensemble du bassin-versant des Grands Lacs, par rapport à une région au sud-ouest du lac Michigan située au-delà de la limite du dernier maximum glaciaire. Le calcul du soulèvement des sites en fonction des altitudes de la surface de référence a facilité la conversion d’un modèle altimétrique de terrain du bassin actuel des Grands Lacs en des cartes paléo-géographiques pour différents âges choisis. De plus, afin de définir l’évolution des niveaux lacustres, les altitudes initiales des indicateurs géomorphologiques et stratigraphiques des paléo-niveaux lacustres, datés au 14C, ont été reconstituées puis reportées en fonction de l’âge. La comparaison de cette courbe à celles des paléo-altitudes des exutoires lacustres, calculée indépendamment, révèle, entre 7,9 et 7 ka BP (8,7 et 7,8 cal ka BP), une phase d’abaissement majeur des niveaux lacustres au-dessous de ceux des exutoires et la fermeture hydrologique des lacs dans le bassin des lacs Huron-Michigan. La forte évaporation nécessaire à l’abaissement du niveau d’eau est attribuée à un climat sec peu après le début de l’Holocène, dans un contexte de détournement progressif hors des Grands Lacs des eaux de fonte du lac Agassiz.
Article body
Introduction
The Great Lakes of North America consist of five major basins whose water surfaces comprise 32% of a total watershed area of 766 000 km2 (The Great Lakes Environmental Atlas, 1995) (Fig. 1). The watershed forms the headwaters of the St. Lawrence River that drains to the Gulf of St. Lawrence and Atlantic Ocean. At the last glacial maximum, the watershed was completely covered by the southern flank of the Laurentide Ice Sheet about 18-21 ka BP (21.3-25.4 cal ka BP). Deglaciation of the basins occurred as the ice margin retreated generally in a northerly direction in a series of oscillations, first exposing the Erie basin about 15.5 ka BP (18.8 cal ka BP), and finally receding to the northern Superior basin about 9.5 ka BP (10.7 cal ka BP) (Table I; Dyke et al., 2003; Dyke, 2004). During this retreat, a series of proglacial lakes formed shorelines of different ages that are upwarped today towards the north-northeast in the direction of thicker and longer-lasting ice (Fig. 2). This deformation is the cumulative isostatic adjustment of Earth’s crust since formation of the shorelines as a result of removal of the former ice sheet load. This differential rebound proceeded throughout the period of ice retreat and postglacial time at decelerating rates of uplift. It is continuing today, as evidenced by tilting of the Great Lakes basins measured in long-term records of lake level gauges between southern and northern shores of the Great Lakes (CCGLBHHD, 1977; Tushingham, 1992; Mainville and Craymer, 2005).
Figure 1
Shaded-relief map of the present-day Laurentian Great Lakes with their drainage areas outlined in white from the National Atlas of Canada (1985). The map file is of 30 arc-second resolution, and is from the Canadian Hydrographic Service (1996) with bathymetry added for Georgian Bay and the area west of LaCloche Lowland from United States Lake Survey (1965). Numbers refer to lake-level indicator sites listed in Table II. Lake outlets shown by black arrows.
Carte du relief de la région des Grands Lacs et des limites de ses sous-bassins de drainage (en blanc) établie à partir de l’Atlas national du Canada (1985). La base de données possède une résolution de 30 secondes d’arc et provient du Service hydrographique du Canada (1996), où la bathymétrie ajoutée de la baie Géorgienne et du secteur à l’ouest des basses terres de LaCloche provient de la Commission des Lacs des États-Unis (1965). Les chiffres renvoient aux sites d’indicateurs du niveau des lacs présentés dans le tableau II. Les flèches noires localisent l’exutoire des lacs.
Table I
Sources of isobase and deglacial information
The differential and continuing nature of rebound has raised indicators of former lake levels to differing elevations, and this makes their correlation difficult and uncertain, particularly for those that are at scattered locations within the basin. This difficulty is especially evident for lake-level indicators that are now below the surface of the present Great Lakes. To facilitate the reconstruction of original elevations of former lake-level indicators, a reference isostatic response surface with an age of 10.6 ka BP is constructed for the entire basin. This surface is used in conjunction with an exponential function to describe the vertical motion in the time domain of any location in the entire watershed of the Great Lakes. This motion is constrained by, and consistent with, the observed empirical data provided by well-known upwarped strandlines in the Great Lakes basins. With this reference response surface and the exponential expression for uplift, the original elevations of specific lake-level indicator sites and potential overflow outlets are reconstructed and compared. The paleogeography of individual basins and of the entire Great Lakes watershed may likewise be reconstructed for any desired age. Applications of the isostatic response surface and the exponential model of uplift are illustrated for: (1) digital reconstruction of Great Lakes paleogeography in a geographic information system (GIS) environment, and (2) the assessment of former lake-level indicators in relation to possible overflow outlets, leading to the conclusion that the Michigan, Huron and Georgian Bay basins, if not all the Great Lakes, were once hydrologically closed.
Figure 2
A
B
Maps of isobases of selected reference paleoshorelines in basins of the Great Lakes, showing the radiocarbon ages (ka BP) and trends of isobases with their present lowest and highest elevations. Sources of isobases given in Table I. C and D mark end points of a section to which isobases of the Iroquois, Algonquin and Nipissing lakes were projected for illustration in Figure 3. (A) Map of Laurentide maximum ice margin and isobases of lakes Wisconsin, Washburn, Minong, Algonquin, Glenwood II, Whittlesey and Iroquois. (B) Map of isobases for the Nipissing Great Lakes, and for an older phase of glacial Lake Wisconsin.
Cartes des isobases des anciennes lignes de rivage de référence de la région des Grands Lacs, où les âges au 14C (ka BP) et le soulèvement différentiel par rapport aux altitudes actuelles sont illustrées. Les sources bibliographiques sont présentées au tableau I. Les points C et D sont les limites du profil vertical projeté à la figure 3 des lacs Iroquois, Algonquin et Nippising. (A) Carte de l’extension glaciaire maximale de l’inlandsis laurentidien et des isobases des lacs Wisconsin, Washburn, Minong, Algonquin, Glenwood II, Whittlesey et Iroquois. (B) Carte des isobases de la phase Nipissing des Grands Lacs et pour une phase plus ancienne du lac Wisconsin.
Figure 3
A plot of relative uplift vs. age for the Nipissing, Algonquin and Iroquois strandlines showing a relatively good fit to the exponential uplift equation for relaxation times (τ) between 3000 and 5000 years. The relative uplifts are the differences in projected strandline elevations at C and D shown on Figure 2.
Graphique temporel du soulèvement relatif des lignes de rivage des lacs Nipissing, Algonquin et Iroquois montrant un très bon ajustement à l’équation exponentielle de soulèvement calculée pour des temps de relaxation (τ) de 3000 à 5000 ans. Les soulèvements relatifs correspondent aux différences entre les altitudes projetées entre les points C et D de la figure 2.
Uplift in the Time Domain
Isobases of Glacio-Isotatic Rebound
In the continental interior, isobases of glacio-isostatic rebound are usually defined by the elevations of a differentially uplifted shoreline of a former lake (Goldthwait, 1907, 1910; Leverett and Taylor, 1915; Hough, 1958; Walcott, 1972; Lewis and Anderson, 1989; Schaetzl et al., 2002). Sets of isobases selected for this study for basins within the Great Lakes watershed are illustrated in Figure 2A and 2B in which the lowest and highest elevations, trends of isobases, the name and uncalibrated radiocarbon age are shown for each lake. Although the Algonquin and Nipissing highstands were confluent in three of the basins about 10.6 ka BP (12.6 cal ka BP) and 5.0 ka BP (5.7 cal ka BP), respectively, most sets of isobases are confined to a single basin. These former, once-level lake surfaces are all warped upward in a north to northeasterly direction (Fig. 2A-B) as a result of differential glacio-isostatic recovery of Earth’s crust following deglaciation of the Laurentide Ice Sheet. The sources of isobase and related deglacial information for basins within and adjacent to the Great Lakes are listed in Table I.
Uplift as a Function of Time
Whereas isobases portray the present configuration of a rebounding surface, it is often desirable to describe the uplift of a given location through time, or to construct surfaces at intermediate times, especially when reference shorelines are widely spaced in time. Following Andrews (1970) and others, the exponential function is adopted to describe relative uplift (Ut) with time where time is expressed as age t (cal years BP) for a landscape previously loaded by an ice sheet. From Peltier (1994, 1998) we use:
where A and τ (tau) are parameters of the equation. A is a site-specific amplitude factor, and is evaluated as:
for a known relative uplift, age and relaxation time. Tau (τ) is the relaxation time or period in years for which decelerating uplift is reduced by 1/exp (1/2.7183 or 36.8%) in successive periods. It is evaluated by solving (Equation 1) at sites where Ut is known for at least two sets of isobases of different ages as shown below. As a first-order approximation, τ is assumed to be time invariant and similar in value throughout the Great Lakes region. Figure 3 shows the reasonably good fit of the exponential uplift curve to Great Lakes relative rebound data for the Nipissing, Algonquin and Iroquois shorelines between sites C and D for trial values of relaxation time between 3000 and 5000 years (see Fig. 2 for location and isobase values).
Figure 4
Map of the Great Lakes basin showing locations of transects where isobases of two shorelines of known age (shown as 14C ka BP) provide a basis for computing the relaxation time of the uplifting process. The mean value of these estimates, 3700 ± 700 years, is used as the relaxation time throughout the basin.
Carte du bassin des Grands Lacs montrant l’emplacement des transects, où les isobases de deux lignes de rivage d’âge connu (montrés en 14C ka BP) ont servi comme base de calcul du temps de relaxation du processus de soulèvement isostatique. La valeur moyenne de ces estimations, 3700 ± 700 ans, est utilisée comme temps de relaxation sur l’ensemble du bassin.
Evaluation of parameters (τ and A) in the relative uplift equation (1)
Relaxation time (τ) was determined on 20 specific short transects (Fig. 4) where relative uplifts, U1, U2, and their ages t1, t2, are known for two shorelines, for example the Algonquin and Nipissing isobases whose domains largely overlap one another. Then, on each transect:
As the amplitude factor is identical on a specific short transect, these equations were rearranged in terms of A and subtracted to yield a single equation which was solved for τ:
From Figure 4, a mean value of τ, rounded to 3700 ± 700 years, is used for computations of relative rebound in the Great Lakes basin. Similar values of 3400 years and 3500 ± 400 years, respectively, were obtained for the relaxation time of an exponential fit to relative sea-level changes in rapidly-uplifting James Bay south of Hudson Bay (Fig. 3A and 3C; Peltier, 1998), and for relaxation of glacial rebound in the Lake Winnipeg area, Manitoba (Lewis et al., 2000; Brooks et al., 2005).
The amplitude factor is then determined by evaluating Equation 2 using the known decay time τ:
Relative uplift (Ux) or shoreline slope (Sx) at other times
With τ and A known, relative uplift (Ux) or shoreline slope (Sx) for any given age tx (cal years BP) in the same transect can be computed:
Ux can be either larger or smaller than U1 or U2, and is the basis for removing the effects of rebound from a modern DEM for construction of the watershed paleogeography at a given time tx, as described in a later section. Average shoreline slope (Sx) = relative uplift (Ux) / transect length (d) or:
Construction of a Multi-Basin Response Surface of Algonquin-Age Rebound About 10.6 Ka Bp (12.6 Cal Ka Bp)
Profiles perpendicular to isobases were selected on which an uplifted surface of Algonquin age could be constructed throughout the Great Lakes watershed (see locations of profiles on inset map of Fig. 5). Starting with area B on the Michigan-eastern Superior profile, the slope of the Elderdon phase shoreline of Lake Wisconsin (about 14 ka BP or 16.7 cal BP) was adjusted to the selected reference age, 10.6 ka BP (about 12.6 cal ka BP) using Equation 5 (Fig. 5). Then the reference profile was extended to the northeast using the Glenwood II shoreline. The uplift of the second isobase of the Glenwood II shore relative to its lowest isobase was adjusted for uplift using Equation 5. This procedure was repeated for each Glenwood II isobase until a complete shore profile, adjusted to the reference age, could be plotted. The adjusted Glenwood II profile was then shifted vertically to connect with the uplifted end of the Wisconsin profile. Next, the Algonquin profile, already at the selected reference age, was shifted vertically to connect with the end of the adjusted Glenwood II profile. The Minong isobases were then adjusted from an age of 9.5 ka BP (10.7 cal ka BP) to 10.6 ka BP (12.6 cal ka BP), and the adjusted shore profile was shifted vertically to connect with the end of the Algonquin profile. The continuous curve rising from zero to about 200 m at 800 km distance and beyond represents uplift on the Michigan-eastern Superior profile since 10.6 ka BP (12.6 cal ka BP). This rebound is relative to area B just beyond the Laurentide ice limit in the southwestern corner of the Great Lakes map area.
Figure 5
Michigan-east Superior profile showing construction of the gradients of shorelines representing uplift since the reference age of 10.6 ka BP (12.6 cal ka BP) relative to area B beyond the limit of the last glaciation. Dots are observed shoreline isobases. Dotted, dashed and solid lines are shore gradients adjusted to an age of 10.6 ka BP (12.6 cal ka BP). Inset map shows locations of reference profiles for contouring the reference response surface.
Profil du Michigan-est Supérieur qui montre le soulèvement des lignes de rivage depuis l’âge de référence de 10,6 ka BP (12,6 cal ka BP) par rapport à la région B située au-delà de la limite de la dernière glaciation. Les points sont les isobases des lignes de rivages observées, les lignes pointillées, les lignes tiretées et les lignes continues représentent les gradients de rivage ajustés à un âge de 10,6 ka BP (12,6 cal ka BP). L’insertion montre l’emplacement des profils utilisés pour délimiter la surface de référence de la réponse isostatique.
The adjusted Algonquin-age uplift values were transferred along connecting isobases to the other profiles (inset map on Fig. 5). After populating these profiles with Algonquin-age uplift values, contours were drawn throughout the region honouring the profile data and isobase trends. These contours constitute the reference response surface (Fig. 6) for the Great Lakes basin, representing isostatic rebound of the region from 10 600 BP (12 600 cal BP) to the present relative to area B. For a response surface at any other age, tx, the contour lines remain the same, but each uplift contour value becomes Utx using Equation 5, i.e.
where Acontour is the amplitude factor of the contour, such that:
and U10.6 ka is the value of the contour from the reference response surface.
Paleogeographic Reconstruction
The use of an isostatic response surface in paleogeographic reconstructions is described by Leverington et al. (2002) for Arctic Canada, and is similar to the use of relative-sea-level isobase maps as employed for reconstruction of Atlantic Canada paleogeography (Shaw et al., 2002). For paleogeographic reconstruction of the Great Lakes basin at a specific age, values for contours of an isostatic response surface at the specific age were first computed as outlined above in Equations 7 and 8. The new reference uplift contours, expressed as a vector data set of isobases, were transformed to a gridded data set or surface using the Triangulated Irregular Network (TIN) interpolation method. The interpolated values from this surface were subtracted from each corresponding pixel value of the modern Great Lakes digital elevation model (DEM) (Fig. 1) to generate a paleo-DEM for the desired age. Paleo-lake shorelines were determined by contouring the paleo-DEM within individual basins according to the elevation of their outlet sill or constriction, which is known from geological data. Paleo-DEM pixel values for areas within the shoreline contours were subtracted from the shoreline elevations to express paleo-lake water depths, which in turn were used to calculate lake area and volume. When present in the map area, an ice cover was superposed using information from the glacial geological literature (Table I) and from syntheses of deglaciation such as Barnett (1992) and Dyke et al. (2003). Twelve paleogeographic reconstructions, ranging from the Kirkfield Algonquin phase at 11.4 ka BP (13.3 cal ka BP) to the Nipissing Great Lakes at 5 ka BP (5.7 cal ka BP), were compiled using a GIS; images of two of these reconstructions are shown in Figure 7.
Figure 6
The reference isostatic response surface portraying glacial rebound since 10.6 ka BP (12.6 cal ka BP) in the Great Lakes region relative to an area southwest of Lake Michigan beyond the limit of the last glaciation.
Surface de référence de la réponse isostatique montrant le relèvement depuis 10,6 ka BP (12,6 cal ka BP) pour la région des Grands Lacs par rapport à un secteur au sud-ouest du lac Michigan se situant au-delà de la limite maximale de la dernière glaciation.
Application to Paleohydrological Modeling
The 12 GIS paleogeographic reconstructions based on the empirical model of isostatic adjustment were used to measure land and water areas, and lake volumes for studies of the paleohydrology and meltwater flow through the Great Lakes system (Moore et al., 2000). Variations in the basin attributes in the 11.4 to 5.0 ka BP (13.3 to 5.7 cal ka BP) period were substantial (Fig. 8). Maximum variation in individual basins under overflow conditions ranges from +72% to ‑95% for lake area, and from +200% to ‑97% for lake volume, compared to the present Great Lakes.
Reconstruction of Former Lake Levels and Discovery of Closed Lowstands, 8900-7800 CAL BP (7.9-7.0 KA BP)
Seventy-seven radiocarbon-dated, and two other, indicators of former lake levels in the Georgian Bay, Huron and Michigan basins (Table II and examples shown in Fig. 9) have been restored to their original elevations using a site uplift equation:
Here, the elevation at time t cal BP is equated to the present elevation (Ep) minus the uplift of the site since time t. This equation, based on Equations 1 or 5, is used to compare the relative altitudes of two or more lake-level indicators or outlets at various times during their isostatic adjustment.
Original elevations for the reported error limits of each dated lake-level indicator were computed using Equation 9 and are shown in Table II. Most indicators are plotted in Figure 10 for the Huron and Georgian Bay basins, and in Figure 11 for the Michigan basin for the interval 11.7 to 6.2 ka BP (13.5 to 7.1 cal ka BP). On these figures, pairs of symbols joined by tie lines indicate the age and reported error range (X axis) in the date of each water surface indicator (Table II). The original altitude of an indicator is the Y-axis value of the plotted symbol; the tilt of the tie line represents uplift during the reported error range of its age. See Figure 1 for locations of the numbered indicators. The lake level histories were interpreted from the restored elevations of the radiocarbon-dated geologic indicators of former water levels in Huron and Michigan basins. For the times when lakes in the Michigan and Huron basins were confluent (relatively high levels), water levels in the Michigan basin were transferred from the Huron basin diagram.
Figure 7
A
B
Paleogeographic maps showing reconstructions of the bathymetry and topography of the Great Lakes basin. (A) Main Lake Algonquin phase about 10.55 ± 0.1 ka BP (12.4-12.7 cal ka BP). At a late stage of Main Lake Algonquin, shown here, the lake had expanded to a possible maximum area in Superior (dashed line) and Huron basins by calving of icebergs from glacier margins in deep water. This lake overflowed via Port Huron to the Erie basin. Early Lake Erie, controlled by the Lyell-Johnson sill near the present Niagara Falls, overflowed via Niagara River to Early Lake Ontario which discharged via the emerging St. Lawrence River to Champlain Sea. (B) Mattawa highstand phase about 8.7-8.8 ka BP (9.6-9.9 cal ka BP), possibly a phase of high discharge when lake level in the upper Great Lakes basins was hydraulically dammed by constrictions downstream from North Bay at either the Rutherglen Moraine or the Rankin Constriction or both (Fig. 1). Discharge continued via Ottawa River to St. Lawrence River. There was no drainage from the upper to lower Great Lakes at this time; Lake Erie, Niagara River and Lake Ontario drained a separate watershed to the upper St. Lawrence River. Legend as for Figure 7A.
Cartes paléo-géographiques de la reconstitution de la bathymétrie et de la topographie du bassin des Grands Lacs. (A) Phase du lac Algonquin principal en 10,55 ± 0,1 ka BP (12,4-12,7 cal ka BP). Au cours du stade tardif de cette phase, le lac Algonquin a atteint son extension maximale dans les bassins du lac Supérieur (tireté) et Huron par le vêlage des icebergs issus des marges glaciaires en eau profonde. Ce lac se déversait dans le bassin du lac Érié par le Port Huron. Le lac Érié initial, dont le niveau était contrôlé par le seuil de Lyell-Johnson situé près des chutes actuelles du Niagara, se déversait dans le lac Ontario initial par la rivière Niagara, se déversant à son tour (le lac Ontario) dans le Saint-Laurent naissant par la Mer de Champlain. (B) Phase de haut niveau lacustre de Mattawa en 8,7-8,8 ka BP (9,6-9,9 cal ka BP), une phase de fort débit probablement liée à la fermeture hydraulique des Grands Lacs par un étranglement topographique situé en aval de North Bay provoqué par la moraine de Rutherglen, l’étranglement de Rankin, ou les deux. Le déversement dans le Saint-Laurent s’effectuait par la rivière des Outaouais. Il n’y avait pas de drainage entre l’amont et l’aval des Grands Lacs à cette époque : les lacs Érié et Ontario et la rivière Niagara drainaient un bassin-versant distinct vers le haut Saint-Laurent. La légende est la même qu’en A.
Figure 8
Bar graphs showing variation in the land drainage areas, lake surface areas, and lake volumes of selected Great Lakes overflowing phases measured from reconstructions of the paleo-Great Lakes in a GIS computing environment. Triangles indicate adjacent lakes that were and are confluent. Names at right refer to lake phases in Huron and other upper Great Lakes basins.
Diagrammes de la variation des superficies de drainage, des surfaces et des volumes des Grands Lacs au cours des différentes phases de déversement obtenues à partir des reconstitutions des paléo-Grands Lacs par un système d’information géographique (SIG). Les triangles indiquent s’il y a des lacs adjacents ayant ou encore connectés. Les noms à la droite font référence aux phases lacustres du bassin du lac Huron et des autres Grands Lacs.
The accuracy for estimating lake levels from empirical indicators is generally considered here to be at best about ±1 m based on survey error and variations in elevations of strandline bluffs, beaches, bars and spits (Schaetzl et al., 2002). The computed initial elevations for lake-level indicators are considered accurate within about 2 m in areas of good isobase control, and somewhat >2 m in areas where constraining isobases have been extrapolated.
A second set of data (Table III) comprising the uplift history of lowest-possible potential overflow outlets for the Huron and Michigan basins, was also computed in the same manner and plotted on the diagrams of Figures 10 and 11 as dotted bands and dashed lines. These constraints include sills at North Bay (dotted band), Dalles Rapids, Deane-Tovell Saddle, and the head of Mackinac River (Lewis and Anderson, 1989). The North Bay sill controlled the overflow level of the water in the combined Georgian Bay, Huron and Michigan basins during the Nipissing transgression (Karrow 1980; Monaghan et al., 1986; Colman et al., 1994a, 1994b). The Dalles Rapids sill controlled overflow levels of the Georgian Bay basin during earlier phases while the North Bay sill was isostatically depressed at lower elevations. The Lucas and Fitzwilliam channel sills were selected as representative of constraints on overflow water levels in the northern Huron basin by the Deane-Tovell Saddle between Bruce Peninsula and Manitoulin Island. The Michigan basin overflow water levels were controlled by a sill at the head of the now-submerged Mackinac River channel between Michigan and Huron basins beneath northern Lake Michigan and the Straits of Mackinac (Stanley, 1938b).
Interpretation of Huron and Georgian Bay Water-Level History
13 500-11 000 cal BP
In this synthesis, the interpreted water-level history is shown as a thick line in Figures 10A-B and 11, and begins with the transgression of the Kirkfield Algonquin lake level to the Main Algonquin level (sites 4, 2, 1c, 1b, 3) under influence of the isostatically-rising Kirkfield outlet (Fig. 1). After 12 600 cal BP (10 600 BP), lake levels fell, indicated by the initiation of gyttja sedimentation in isolated basins (sites 12, 13, 10), through the Post-Algonquin phases and below, evidenced by the onset of organic accumulation on Manitoulin Island (23) and in Georgian Bay (7b). Land emergence by this fall of lake levels through the Post-Algonquin phases is recorded further by additional, but somewhat delayed, organic sediment accumulation in small basins south of Georgian Bay (5, 6), near North Bay (24, 25) and on Manitoulin Island (8).
Figure 9
A
B
C
D
E
F
G
H
Types of evidence of former water levels in the Huron, Michigan, and Erie basins. Abandoned geomorphic shore features: (A) Former lake surfaces are inferred from erosional coastal features at the change in gradient between the gently-sloping shoreface and the steeply-sloping shorebluff, illustrated here for the Algonquin shore near Kirkfield, Ontario. (B) Storm beach ridges above Lake Huron on Manitoulin Island (air photo courtesy of Natural Resources Canada). (C) and (D) Beach ridges 53 m below Georgian Bay imaged by sidescan sonar and multibeam sonar, respectively (Blasco 2001). (E) Submerged beach and shoreface 21-28 m below Lake Erie (from Coakley and Lewis, 1985). Isolation basins, tree stumps, and unconformities: (F) Air photo (courtesy of Natural Resources Canada) view of a small basin isolated during the regression of a large lake, Manitoulin Island. A radiocarbon date of the contact zone between gray clastic large lake sediment and dark-coloured organic small-lake-sediment (photo courtesy of T.W. Anderson) provides the age at which the Great Lake surface passed below the sill elevation of the small isolated basin (Lewis, 1971). (G) In situ tree stump on lakefloor at entrance to Georgian Bay. (H) Subsurface reflections in a seismic profile indicate offshore unconformities and sequence boundaries caused by episodes of reduced lake level (Moore et al., 1994).
Indicateurs des anciens niveaux lacustres dans les bassins des lacs Huron, Michigan et Érié. Formes littorales abandonnées : (A) les surfaces des anciens lacs sont estimées à partir de formes d’érosion côtière qui sont situées à la rupture de pente entre la face de la côte en pente douce et le talus littoral en pente abrupte, comme le montre cet exemple du rivage du lac Algonquin près de Kirkfield, Ontario. (B) Crêtes de plage de tempête en bordure du lac Huron sur l’Île Manitoulin (courtoisie de Ressources Naturelles Canada). (C) et (D) Crêtes de plage à 53 m de profondeur dans la baie Géorgienne identifiées respectivement à l’aide du sonar latéral et d’un sonar multi-faisceaux (Blasco, 2001). (E) Plage et avant-côte submergées entre 21 et 28 m sous la surface du lac Érié (d’après Coakley et Lewis, 1985). Bassins isolés, souches d’arbres et discontinuités : (F) Photo aérienne d’un petit bassin isolé durant la régression d’un grand lac, Île Manitoulin (Source : Ressources Naturelles Canada). Une datation au radiocarbone de la zone de contact entre les sédiments clastiques gris du grand lac et les sédiments organiques foncés du petit lac (photo fournie par T.W. Anderson) fournit l’âge auquel le plan d’eau du grand lac est passé sous le niveau du seuil du petit bassin isolé (Lewis, 1971). (G) Souche d’arbre in situ au fond du lac à l’entrée de la baie Géorgienne. (H) Profil sismique montrant les discontinuités et les limites causées par des épisodes de bas niveau lacustre (Moore et al., 1994).
Table II
Indicators of former lake levels in the Michigan, Huron and Georgian Bay basins, and Georgian Bay-North Bay lowland
Surface elevation (m asl) of sampled lake or bog.
Present elevation (m asl) of feature representing former level of a Great Lake.
Uplift (m) since 10.6 ka BP (12.6 cal ka BP) from reference isostatic surface (Fig. 6).
Elevations (m asl) computed using equation (9) for the age limits of the calibrated age range of the lake-level indicator as listed in column 7.
Ages at limits of the calibrated range (1σ) of the lake-level indicator date and its uncertainty reported in column 3. Calibration by the Calib 5.0.1 program
(Stuiver and Reimer, 1993) using the INTCAL98 calibration dataset (Stuiver et al., 1998).
Shell age uncorrected for reservoir effect.
n.r. = not reported.
PSV = Age by comparing paleomagnetic secular variation with a 14C-dated standard for northeastern U.S.A. (King and Peck, 2001).
Table III
Present and early Holocene elevations1 of North Bay, Georgian Bay, Huron and Michigan sills (m above present sea level)
Computed using Equation 9.
Sill elevation; full discharge water level is 9 m. higher. Estimated from Lewis (1969).
Drom Lewis and Anderson (1989).
This study
11 000-10 000 cal BP
Relative lake levels continued to fall toward the level of driftwood deposited in Georgian Bay (31b). During the lowstand, trees grew in the Straits of Mackinac (28b), and organic sediments accumulated at relatively low-elevation at northern and southern Huron basin sites (29c, 30f). A subsequent abrupt lake-level rise and inundation is recorded by a shift to more aquatic species in Georgian Bay peat at site 7 (Lewis and Anderson, 1989). The lake level rise at about 10 800 cal BP (9500 BP) was possibly in response to the onset of outburst floods from upstream glacial Lake Agassiz through its eastern outlets to Superior basin (Clayton, 1983; Teller and Thorleifson, 1983; Farrand and Drexler, 1985; Lewis and Anderson, 1989; Breckenridge et al., 2004; Breckenridge and Johnson, 2005). The rise of lake level during high-discharge events may be related to constrictions in the Mattawa River valley at the Rutherglen Moraine (Harrison, 1972; Chapman, 1975) and in Ottawa River valley at the Rankin Constriction (Lewis and Anderson, 1989) (Fig. 1). The maximum of the lake-level rise possibly exceeded the altitude of the Sheguiandah site on Manitoulin Island (9c, 9d) where swamp organic sediment began accumulating just after the lake-level rise. Some sites with thick clay beneath gyttja (41, 33) suggest the waters of this inundation carried abundant fine-grained sediment in suspension. This pulse of higher water (Lake Mattawa) was apparently short-lived, and its decline may account for the emergence of small basins and onsets of gyttja sedimentation or peat accumulation at several sites (9c, 9d, 29c, 30b, 30c, 30d, 30f, 33, 34b, 38, 40, 41).The lake level drop was low enough to expose the eastern Huron basin lakebed for tree growth at site 66. Within a century or so, rising water level allowed for deposition of driftwood (64 then 32c), and terminated peat accumulation at slightly higher sites in southern Huron (39) and Georgian Bay (7a) basins.
10 000-9000 cal BP
The rising lake level that terminated peat accumulation (7a) continued to rise to a second high-level Lake Mattawa phase (about 9800 cal BP, 8700 BP) that inundated Pure Lake (50), and possibly culminated at a high enough level to induce sediment aggradation in Amable du Fond River (49) and to construct a barrier lagoon north of Lake Huron (43) (Fig. 10B). Thick clay deposits beneath gyttja in Pure Lake suggest that this water-level rise was also accompanied by a high influx of fine-grained suspended sediment. An extended period of moderately low water level following the Mattawa highstand is indicated by tree growth on the saddle between Huron and Georgian Bay basins (67, 68a, 68b), and by peaty marsh sediments in southwestern Huron basin (30a) and South Bay, Manitoulin Island (42a, 42b, 42c). Oddly, an erosional event and sequence boundary in the deep water sediments of Huron basin is not recorded at this time of low water level (Fig. 10; Moore et al., 1994; Rea et al., 1994a).
At about 9050 cal BP (8100 BP) water levels again rose, transforming sedimentation in South Bay from plant detritus to clastic silty clay (42a, 42b), and possibly overflowing the North Bay outlet where delta construction occurred at nearby Trout Mills (48). Water levels then rose abruptly to the last Mattawa highstand, likely as a result of overflow from large volumes of subglacial meltwater discharged during the Nakina ice advance over the Great Lakes-Hudson Bay divide, and recorded throughout the Superior basin as a set of 36-40 thick varves dated 9035 ± 170 cal BP or about 7950-8250 BP (Breckenridge et al., 2004). Although the maximum lake rise in Huron and Michigan basins is not known, it apparently fell short of the Sheguiandah archeological site on Manitoulin Island as it is not recorded in the peat sequence there between dated samples 9b, 9c and 9d (Anderson, 2002). Land emergence following this short-lived high Mattawa phase lake is documented by onsets of organic sedimentation at Wye Marsh, southern Georgian Bay (72), gyttja deposition in small basins northeast of Georgian Bay (44, 45, 46, 47), and tree growth at low elevations now under water in Huron basin (75a, 75c) and possibly in the Straits of Mackinac (28a).
Figure 10
A
Plot of original elevations of numbered water-level indicators (Table II) in Huron and Georgian Bay basins relative to area B (Fig. 5), inferred lake surface, and potential overflow outlet elevations (Table III) vs. age in calibrated and radiocarbon years BP. Downward-pointing triangles indicate data implying a lake level below the plotted site, as for in situ tree stumps. Upward-pointing triangles are for data implying a lake level above the plotted site, as for shallow water mollusks or driftwood. Right-pointing triangles indicate data implying isolation of a small basin by a falling large lake level. Left-pointing triangles indicate data implying a transgression by a large lake. See Figure 10B for a plot of all data in the 10 000-7500 cal BP interval. The inferred lake level for Huron basin is shown by the thick continuous line (dotted for Georgian Bay basin where different). The elevation history of potential overflow outlets (Table III) is portrayed by a dotted band (North Bay–NB) or by dashed curves (Huron basin sills–H, Georgian Bay basin sill–G, and Michigan basin–M). The Algonquin lake phases are Kirkfield Algonquin (K), Main Algonquin (M), and Post-Algonquin (P). MS and ES are middle and early Stanley lowstands in Huron basin, respectively. In Georgian Bay basin, the equivalent lowstands are named Hough. The horizontal lines above the X-axis, named by colours, indicate age and duration of lowstand sequence boundaries inferred from seismic reflection and core data from Moore et al. (1994) and Rea et al. (1994a).
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres (tabl. II) pour les bassins du lac Huron et de la baie Géorgienne relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires (tabl. III), l’âge en années étalonnées et conventionnelles. Les triangles pointant vers le bas font référence à un niveau lacustre inférieur au site reporté, comme dans le cas des souches d’arbres in situ. Les triangles pointant vers le haut font référence à un niveau lacustre supérieur au site reporté, comme dans le cas des mollusques en eau peu profonde ou de bois flottant. Les triangles pointant vers la droite font référence à l’isolement d’un petit bassin suite à baisse du niveau d’un grand lac. Les triangles pointant vers la gauche font référence à la transgression d’un grand lac. Voir la figure 10B pour le diagramme des données comprises dans l’intervalle 10 000-7500 cal BP. Le niveau lacustre reconstitué du lac Huron est représenté par un trait gras continu (en discontinu pour la baie Géorgienne lorsqu’il diffère). L’historique de l’altitude des exutoires potentiels (tabl. III) est représenté par une bande en pointillé (North Bay–NB) ou par les courbes en tireté (seuils du bassin du lac Huron–H, seuil de la baie Géorgienne–G et bassin du lac Michigan–M). Les phases du Lac Algonquin sont les suivantes : Kirkfield Algonquin (K), Algonquin principal (M) et Post-Algonquin (P). MS et ES identifient respectivement les niveaux moyen et ancien de Stanley du lac Huron. Dans le bassin de la baie Géorgienne, les niveaux équivalents sont nommés Hough. Les lignes horizontales au-dessus de l’abscisse, nommées par une couleur, indiquent l’âge et la durée des limites déduites à partir de données de réflexion sismique et de forage issues des travaux de Moore et al. (1994) et de Rea et al. (1994a).
B
Expanded plot of original elevations of water-level indicators for the 10 000-7500 cal BP period in Huron and Georgian Bay basins relative to area B (Fig. 5), inferred lake surface, and potential overflow outlets vs. age in calibrated and radiocarbon years BP. Symbols as described in Figure 10A.
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres pour la période 10 000-7 500 cal BP pour les bassins du lac Huron et de la baie Géorgienne relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires, l’âge en années étalonnées et conventionnelles. Les symboles sont les mêmes qu’en A.
9000-7000 cal BP
Lake-level decline from the last high Mattawa phase continued down to levels defined by the Stanley unconformity (77) in Huron basin (Lewis et al., in press) and the Flowerpot beach (FB) in Georgian Bay basin (Fig. 10). This lowstand phase, which continued for several centuries, formed the Light Blue reflector and sequence boundary in deep water sediments (Moore et al., 1994; Rea et al., 1994a), and allowed tree growth on the Deane-Tovell saddle between Huron and Georgian Bay basins (65, 69, 70). Lake-level rise from this low phase by 7.5 ka BP (8.3 cal ka BP) is suggested by the chronology of the Light Blue reflector in northern Huron and Georgian Bay basins (Moore et al., 1994).
The final increase in lake levels from the low phases of late Lake Stanley and late Lake Hough in Huron and Georgian Bay basins, respectively, probably occurred about 8000 cal BP, late enough to allow tree growth in eastern (71) and western (52) Huron basin, but early enough to transgress the Rains Lake site (53b) in northwestern Huron basin. At this elevation, the water surface had risen to the level of the Nipissing beach at North Bay, indicating that the lake in Michigan-Huron-Georgian Bay basins was overflowing the North Bay outlet at full discharge. As this outlet uplifted faster than other parts of the lake basin, relative lake levels rose throughout most of the upper Great Lakes basins as the well-known Nipissing transgression, manifested in the present data by transgression of the Rains Lake (53b) and Smoky Hollow Lake (34a) sites (Fig. 10A).
Interpretation of Water-Level History of the Michigan Basin
Because of their connection via the Indian River lowland and Straits of Mackinac (Fig. 1), high water levels were always at common elevations in the Huron and Michigan basins, as at present, following retreat of ice from their northern regions (Eschman and Karrow, 1985; Hansel et al., 1985). Consequently, the history of lake-level variation above the Michigan basin sill defined in the Huron basin is applicable to the Michigan basin as shown in Figure 11. Only evidence of lowstands and other unique indicators of Michigan basin lake elevations are shown and discussed here.
Submerged in situ tree stumps in southern Lake Michigan (the Olson Forest of Chrzastowski et al., 1991) have been traditionally interpreted as having been drowned in the Nipissing transgression (Chrzastowski and Thompson, 1992, 1994; Colman et al., 1994b). However, in this analysis, the Olson tree stumps (sites 74a, 74b at 8.2-8.4 ka BP, 9.2-9.4 cal ka BP, in Figs. 1 and 11) appear tens of metres above the North Bay outlet, and could not be affected by the Nipissing transgression. A rise of water level to the final highstand of Lake Mattawa follows the Olson Forest after 8.2 ka BP (Fig. 11). This flooding event seems to have been the cause of forest drowning.
Figure 11
Michigan basin plot of original elevations of water-level indicators relative to area B (Fig. 5), inferred lake surface, and potential overflow outlets vs. age in calibrated years BP and radiocarbon age in ka BP. Water-level indicators for the Huron basin above the Michigan sill are applied to the Michigan basin as these highstand lakes were confluent in both basins. Black symbols indicate data specific to the Michigan basin. MC and EC represent middle and early Chippewa lowstands, respectively. Otherwise symbols and lines as for Figure 10A.
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres pour le bassin du lac Michigan relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires (tabl. III), l’âge en années étalonnées et conventionnelles. Les indicateurs du niveau lacustre du lac Huron situés au-dessus du seuil Michigan sont appliqués au lac Michigan puisque ces plans d’eau élevés étaient connectés aux deux bassins. Les symboles noirs font référence aux données spécifiques du lac Michigan. MC et EC représentent les niveaux moyen et ancien de Chippewa. Les autres symboles et les lignes sont identiques à ceux de la figure 10A.
As in the Huron and Georgian Bay basins, an extreme decline in lake level following the last Mattawa highstand is suggested in the Michigan basin by the occurrence of shallow-water sediments and molluscan fauna in the deepwater environment of central Lake Michigan (51a, 51b, 51c) (Lewis et al., in press). A final recovery of lake level to an overflowing condition at the North Bay outlet probably occurred about 8.1-7.8 cal ka BP, allowing for tree growth in northern Michigan basin (54). From this time forward, relative lake level rose under control of the North Bay outlet as the Nipissing transgression.
Closed Lowstand Conditions
During the period 8.95-8.3 cal ka BP (8.05-7.4 ka BP) lake levels (late Stanley in Huron basin, late Hough in Georgian Bay basin, and late Chippewa in Michigan basin) were below the sill of the North Bay outlet, the lowest possible overflow outlet at the time (Figs. 10B‑11). Thus the late Stanley, late Chippewa, and late Hough phases were hydrologically closed lakes at their lowest level, and as such may have resulted from the impact of a severe dry climate in which evaporative water losses exceeded water inflows by precipitation and runoff. The inference of dry climate is supported by the presence of thecamoebians that indicate a more saline lake environment, and by climate transfer function analysis of pollen assemblages which suggests less precipitation and warmer temperatures at 7.7 ka BP than today in the Georgian Bay area (Blasco, 2001). Other lowstands around 9.3 and 9.8 ka BP (10.5 and 11.2 cal ka BP) may have also been closed briefly.
Figure 12
Paleogeographic lowstand map of Huron, Erie and Ontario basins at about 7.8 ka BP (8.6 cal ka BP) showing the reduced area of inferred closed (terminal) lakes whose shorelines are well offshore from the present lake shores. E = Erie basin, GB = Georgian Bay basin, H = Huron basin, and O = Ontario basin.
Carte paléogéographique du bas niveau lacustre à environ 7,8 ka BP (8,6 cal ka BP) pour les bassins des lacs Huron, Érié et Ontario. La carte montre la superficie réduite des lacs fermés dont les anciens rivages étaient bien au large des rivages actuels. E = bassin du lac Érié, GB = bassin de la baie Géorgienne, H = bassin du lac Huron et O = bassin du lac Ontario.
Possible shorelines for the dry climate-induced closed lowstands at 7.8 ka BP (8.6 cal ka BP) are illustrated in a paleogeographic reconstruction of the Huron, Erie and Ontario basins on Figure 12. The northern Huron lowstand shore is tied to the Stanley unconformity beneath northwestern Lake Huron (Hough 1962; Lewis et al., in press). The Georgian Bay lowstand shore is tied to the low-level Flowerpot beach in the entrance to Georgian Bay (Blasco, 2001). Water bodies in these basins and those beneath southern Lake Huron were isolated closed lowstands, well offshore from the present lake boundaries. The Erie lowstand shore is tied to a submerged beach in eastern Lake Erie (Fig. 9E; Coakley and Lewis, 1985; Lewis et al., 2004) at a lake level that would have extended into the central Lake Erie area before that sub-basin became mostly infilled with sediment (Sly and Lewis, 1972). The Ontario basin is assumed to have been impacted by the same dry climate that affected the other basins. Under these conditions, its lowstand shoreline was inferred at a level that made its basin area-to-lake area ratio equal to that of the Huron basin (Bengsston and Malm, 1997).
Discussion
Geophysical Models of Isostatic Adjustment
With knowledge or estimates of the elastic and viscous properties of Earth’s lithosphere and mantle, respectively, geophysical models compute the crustal isostatic response following a known or inferred history of ice sheet loading and pattern of deglaciation. Clark et al. (1994) used this approach to evaluate isostatic movements during deglaciation of the Great Lakes basin. They successfully illustrated the increasing amplitude of postglacial isostatic rebound towards the north and northeast in the direction of past thicker ice and ice retreat in accordance with the evidence of deformed paleo-lake shorelines. They applied an innovative approach to constrain and calibrate their model by using geological knowledge of large-lake drainage transfers from northern to southern outlets as basins were differentially tilted during the Algonquin and Nipissing phases. Lake-level history was derived by tracking the upward movement of overflow sills through time. Lake-level indicators were not tracked separately from the sills with the result that possible episodes of hydrologic closure could not be detected.
Gravitationally self-consistent models of glacio-isostatic adjustment on a global scale have been progressively developed and improved over the past few decades (Peltier, 1998). These models are constrained and calibrated to relative sea-level histories at ocean-continent boundaries. Model estimates of isostatic adjustment are produced for continental interiors, and one of these, the ICE‑3G model (Tushingham and Peltier, 1991), has been compared favourably with short-term evidence of tilting of the Great Lakes basins derived from trend analysis of lake-level gauge records (Tushingham, 1992). However, over the longer term, the geophysical model results are not in complete agreement with past events in the watershed. The relative performances of two geophysical models and the empirical model were assessed by comparing the elevations of the northern and southern outlets from the upper Great Lakes basins during the Algonquin and Nipissing drainage transfers which are known to occur at about 10.5 ka BP and 5 ka BP, respectively, from independent geological evidence (Karrow et al., 1975; Karrow, 1980; Monaghan et al., 1986). Outflows were transferred from the Kirkfield outlet to the Port Huron and possibly Chicago outlets during the Algonquin transfer (Fig. 7A), and from the North Bay to the same southern outlets during the Nipissing transfer. Estimates of the elevations of the overflow sills and timing for these drainage transfers when northern outlets rose isostatically above southern outlets were provided for the ICE‑4G glacio-isostatic model (Peltier, 1995, 1996) courtesy of W.R. Peltier. Differences between Port Huron and the other outlets which should have been near zero at the times of the drainage transfers are shown in Figure 13. The empirical model performs best with northern and southern sills being within two metres of each other during transfers. The Clark et al. (1994) model is similar, although the best agreement for the Algonquin transfer was obtained at 10 ka BP, rather than the expected 10.5 ka BP, and the Nipissing-age transfer at Chicago differed by 3 m. Variances in sill elevations during transfers for the ICE‑4G model are larger, and ranged from 3.5 to 10 m.
Figure 13
Deviations of sill altitudes when they should have been of equal elevation during drainage transfers from northern to southern outlets of the Algonquin and Nipissing lake phases for the Clark et al. (1994) and ICE‑4G geophysical models and the empirical model used in this paper. Sill elevations, reduced relative to Port Huron at the southern end of the Huron basin, show that the empirical uplift model used in this study best describes (lowest deviations) the relative elevations of outlets during drainage transfers.
Écarts entre les altitudes des seuils lorsqu’elles devraient être égales aux autres durant les transferts du drainage des exutoires septentrionaux aux exutoires méridionaux des phases lacustres Algonquin et Nipissing pour les modèles géophysiques de Clark et al. (1994) et de ICE‑4G et le modèle empirique utilisé dans cette étude. Les altitudes des seuils, ajustées à celles du Port Huron à l’extrémité sud du bassin du lac Huron, montrent que le modèle empirique de soulèvement glacio-isostatique utilisé dans cette étude rend mieux compte (écarts plus faibles) des altitudes relatives des seuils durant les transferts de drainage.
Near and beyond the maximum margins of the ice sheets, the geophysical models predict crustal uplift as a forebulge of modest relief compared with the amplitude of depression beneath the centre of the ice load. This has been demonstrated in ice-marginal Atlantic and Arctic coastal regions (Barnhardt et al., 1995; Dyke, 1998). The forebulge effect has not been recognized from empirical evidence in the continental region adjoining or south of the Great Lakes basin. However, ice marginal areas in this region could have undergone uplift and subsidence associated with the growth, migration and decay of a glacio-isostatic forebulge (Colman et al., 1994a), and some evidence for subsidence exists, for example, anomalies in modern tilting of the Erie and southern Michigan basins (Mainville and Craymer, 2005).
Benefits of geophysical models are their ability to compute both vertical and horizontal earth movements associated with the isostatic process, and to relate these estimates to an absolute datum such as present sea level. Confidence in predictions of past crustal movements in the Great Lakes region by geophysical models, especially deformation over several millennia, will be increased markedly when these models are calibrated and constrained by the empirical observations of glacio-isostatic adjustments within the same region.
Previous Recognition of Low Lake Levels in Huron Basin
An extreme lowering of lake level in the Huron and Michigan basins after the Algonquin phases was predicted by Stanley (1936) when the ice sheet margin receded from the isostatically-depressed lowland drainage route past North Bay, Ontario, to the Mattawa and Ottawa river valleys. Sediment unconformities, discovered by Hough (1955, 1962), confirmed lake-level lowstands in each of the Michigan and Huron basins; these lowstands were named Chippewa and Stanley, respectively. Hough (1962) inferred a delay in isostatic adjustment so that, in his model, the North Bay sill remained at its depressed Algonquin level, and Lake Stanley could be interpreted as an open overflowing water body (Lewis et al., in press).
Although a regional advance (Marquette readvance) of the Laurentide Ice Sheet about 500 km wide reached the southern coast of the Superior basin at 10 ka BP (11.5 cal ka BP) (Lowell et al., 1999), and ice retreat was slowed during the Younger Dryas in eastern Ontario and western Québec (Simard et al., 2003), there is no evidence of ice advance or of delayed isostatic recovery in the Georgian Bay-North Bay outlet region, about 300-500 km east of the Superior basin (Dyke et al., 2003). As a result, it is highly unlikely that rebound was delayed significantly and, at the time of late Lake Stanley (Fig. 10), the recovered outlet was above water level, as was also interpreted by Lewis et al. (in press). Closed lake conditions for late Lake Stanley are also consistent with the preliminary paleoecological findings based on thecamoebian analysis of somewhat higher salinity in late Lake Hough, the equivalent closed water body in the Georgian Bay basin, as reported by Blasco (2001).
Working Hypothesis for Climate Change and the Hydrologically Closed Lowstands
Major shifts in Great Lakes water levels have long been understood in terms of overflowing lakes. Lake elevations change as a result of shifting ice dams during glacial retreat or advance, outlet erosion, or by differential glacio-isostatic adjustment of the outlets (Eschman and Karrow, 1985; Hansel et al., 1985; Larsen, 1987; Barnett, 1992; Larson and Schaetzl, 2001). These mechanisms apply to the Algonquin highstand and the transgression to the Nipissing Great Lakes. The intervening Mattawa high phases (Lewis and Anderson, 1989) are also considered to be overflowing lakes, but with variable levels controlled by resistance at hydraulic constrictions downstream of North Bay to variable high-discharge flows from upstream Lake Agassiz or subglacial drainage.
The discovery of low lake levels below the lowest possible overflow outlet in their basins implies a controlling factor of excess evaporation (water loss) over precipitation and runoff (water supply) during a period of dry climate. These lowstands, the late Chippewa, late Stanley, and late Hough lakes in the Michigan, northern Huron, and Georgian Bay basins, respectively, are thought to reflect a phase of climatically-driven hydrologic closure that has not been recognized previously. The earlier pre‑9.5 ka BP (pre‑10.9 cal ka BP) lowstands, middle and early lakes Chippewa, Stanley and Hough, are inferred to be close in level to, or slightly below, their basin outlets, and thus were possibly also subject to enhanced evaporative losses of water. During these early, low lake phases, sustaining upstream inflows were not available, as Lake Agassiz discharge was diverted away from the Great Lakes by advances of ice in the Superior and Nipigon basins (Thorleifson and Kristjansson, 1993; Lewis et al., 1994). Dry air from the glacial atmospheric circulation over the nearby Laurentide Ice Sheet (Bryson and Wendlund, 1967; David, 1988; Anderson and Lewis, 2002; Wolfe et al., 2004) probably exerted continuous evaporative stress on nearby water surfaces, causing drawdown of these lakes during phases of reduced inflow.
By the time of the latest and longest-duration lowstands (late Chippewa, Stanley and Hough) the effects of glacial atmospheric circulation would have diminished greatly owing to the small size of the remaining Laurentide Ice Sheet, and to its more distant location relative to the Great Lakes basin. At this time (8 ka BP, 8.9 cal ka BP), meltwater drainage from the merged glacial lakes Agassiz-Ojibway began bypassing the upper Great Lakes basins directly into the Ottawa Valley (Veillette, 1994; Teller et al., 2002; Teller and Leverington, 2004). Thus, the Great Lakes watershed suddenly lost a significant source of water supply and was susceptible to the early Holocene dry climate (Edwards et al., 1996).
Atmospheric water supply today can be thought of as a function of the relative time spent over the Great Lakes basin by three major air masses over North America. These masses are the Arctic air from the north (dry and cold), the Pacific air from the west (dry and warm), and the Maritime Tropical air bringing moist, warm air north from the Gulf of Mexico (Bryson and Hare, 1974; Bradbury and Dean, 1993). Once glacial lake drainage began bypassing the Great Lakes, and as the Laurentide Ice Sheet downwasted over Hudson Bay, southward incursions of dry Arctic air, previously blocked by the high Laurentide Ice Sheet, were likely becoming increasingly frequent (Yu and Wright, 2001), and this dry air may have initiated or intensified draw down of lake levels to the late Chippewa, Stanley, and Hough lowstands. The closed lowstands may have been maintained by enhanced evaporation into increasingly strong flows of dry, warm Pacific air from the west. These flows are indicated by abundant evidence of vegetation shifts to drought-resistant plant taxa west of the Michigan basin (Baker et al., 1992; Wright et al., 2004). By 7 ka BP (7.8 cal ka BP), increasing incursions of the Maritime Tropical air mass were delivering sufficient precipitation, indicated by the appearance of mesic forest species in pollen diagrams of the Great Lakes region (Webb III et al., 1998) to convert the Michigan, Huron and Georgian Bay water bodies to open, overflowing lakes, as at present. These changes are consistent with paleovegetation maps which show a rapid northward migration through the Great Lakes region of the mixed-boreal forest biome boundary between 8.0 and 7.0 ka BP (8.9 and 7.8 cal ka BP) (Dyke et al., 2004).
The foregoing hypothesized climatic history is consistent with changes inferred for small Elk Lake in Minnesota, based on a comprehensive study of proxy climatic and limnological indicators (Bradbury and Dean, 1993). Similarly, this climatic history is supported by a coeval depletion in the 18O composition of precipitated carbonate in Deep Lake, Minnesota, attributed by Yu and Wright (2001) to the blocking of southern air masses by more frequent presence of dry Arctic air.
The phase of reduced and closed lakes at the onset of the present hydrologic regime offers a unique opportunity to evaluate the sensitivity of the Great Lakes system to high-amplitude climatic change. Understanding the sensitivity of the lakes to high-amplitude and long-duration change would be a distinct benefit in the light of the need to project and adapt to future changes under global warming which may drive the lakes below instrumentally-observed variability (Mortsch et al., 2000).
Isotopic Composition of Lake Water
As the evaporation process favours concentration of water molecules containing the heavier 18O isotope, previous findings of high concentrations of the lighter 16O isotope, similar to that of glacial meltwater, in fossil valves of Huron basin benthic ostracodes at about 7600 ka BP (Rea et al., 1994a, 1994b; Dettman et al., 1995) appear to contradict the coeval presence of evaporatively-driven lowstands as postulated in this paper and in Lewis et al. (in press). Although surface lake water isotopic composition undoubtedly became concentrated in 18O during seasonal periods of rapid evaporation, it apparently did not greatly influence bottom water composition, similar to the isotopic stratification found for glacial Lake Agassiz (Buhay, 1998; Birks et al., in press). Alternatively, adjustments in the chronology of isotopic events in the Huron and Michigan basins suggest that the low 18O inflow occurred prior to the lowstands (Breckenridge and Johnson, 2005; Breckenridge, in press), and may have remained as bottom water during the evaporative phase. Additional research is needed to resolve the origin of bottom water in the Huron and Michigan basins.
Summary
Postglacial isostatic adjustment in the Great Lakes region is described here in the time domain using an exponential decay expression constrained by the observed cumulative differential deformation of dominant paleo-lake strandlines in individual or groups of basins. Vertical earth movement throughout and following the last deglaciation into the middle Holocene is characterized as progressive, differential uplift relative to an area southwest of southern Lake Michigan, beyond the limit of the last glacial maximum. Rates and amplitudes of uplift increase towards the north-northeast in the direction of deglacial retreat and thicker ice. A mean relaxation time for the uplifting process of 3700 ± 700 years was obtained by averaging solutions for the decay time parameter at 20 transects throughout the basin where isobases (gradients) of two strandlines of different age were known. This relaxation time was used in the exponential expression to adjust the isobase gradients of dominant strandlines in the Great Lakes basins to 10.6 ka BP (12.6 cal ka BP), the approximate age of the well-known Algonquin phase.
Collectively, the Algonquin and adjusted isobases constitute a reference response surface for isostatic adjustment throughout the Great Lakes region. Interpolated values from this surface and the mean relaxation time were used in the exponential uplift expression to determine an ‘amplitude’ factor for uplift at any desired location. With these values and the exponential expression, uplift since any desired age could be computed. The original elevation of a site at any desired age could also be determined by subtracting the computed uplift from the present elevation of the site. This approach was used to transform values of pixels in a DEM for the present Great Lakes region to new DEMs at previous ages. A total of 12 paleogeographic reconstructions for the topography and bathymetry of the Great Lakes basins were prepared for ages between 11.4 and 5.0 ka BP (13.3 and 5.7 cal ka BP). These reconstructions showed that water surface areas ranged from +72% to ‑95%, and lake volumes from +200% to ‑97%, relative to the present lakes. Improvements in the estimation of glacio-isostatic effects can be expected from geophysical models when these are calibrated to the available observations of differential rebound in the Great Lakes basin.
The same empirical approach was applied to reconstruct the original elevations of dated indicators of former lake levels in the Michigan, Huron and Georgian Bay basins for the interval between 11.7 and 6.2 ka BP (13.5 and 7.1 cal ka BP). Original elevations for 79 dated indicators, comprising fossils from beach lagoon sediments, basal organic sediment from isolation basins, shallow-water fossils in unconformable zones within deepwater sediment, and submerged tree stumps in growth position, and others, were reconstructed to form the basis for interpreting lake-level history in these basins. In parallel fashion, the history of potential overflow sills was also reconstructed. Comparison of sill elevations with lake levels revealed a period (about 8.05 to 7.4 ka BP, or 8.95 to 8.3 cal ka BP) in which water surfaces were up to several tens of metres below lowest possible overflow outlets. The low lake levels are postulated to reflect the increased impact of early Holocene dry climate when upstream Agassiz overflow and/or subglacial floods were diverted around the Great Lakes basin, and water supply to the upper Great Lakes was reduced. Lake evaporation and draw down were likely enhanced by frequent incursions of dry cold Arctic air and, later, warm dry Pacific air, possibly related to atmospheric reorganization associated with the demise of the Laurentide Ice Sheet. This period of closed lakes, early in the present hydrological regime of the Great Lakes, offers an opportunity to probe and understand the sensitivity of the Great Lakes system to high-amplitude, long-duration climate change. Such information could improve confidence in projections of, and adaptations to, future levels of the Great Lakes under global warming.
Appendices
Acknowledgements
We thank the Canadian Hydrographic Service for the Great Lakes DEM. W.R. Peltier and R. Drummond, University of Toronto, kindly provided output from the ICE‑4G geophysical glacio-isostatic model. We are grateful to reviewers whose suggestions helped us improve the paper, colleagues T.W. Anderson, J. Shaw, B.J. Todd and S. Occhietti, and S.M. Colman and T. James, reviewers for the journal. S. Occhietti and M. Parent kindly provided the French translations. The research for this paper has been aided in many ways by support from the Climate Change Program of Natural Resources Canada, Earth Sciences Sector, and from the National Science Foundation, Grant ATM‑0354762. Figures were prepared by K. Hale and P. O’Regan, Electronic Publishing Unit, GSC Atlantic.
Note
-
[*]
Geological Survey of Canada contribution number 2005546
References
- Anderson, T.W., 1979. Stratigraphy, age, and environment of a Lake Algonquin embayment site at Kincardine, Ontario p. 147-152. Geological Survey of Canada, Ottawa, Current Research 79-1B.
- Anderson, T.W., 2002. Pollen stratigraphy and vegetation history, Sheguiandah archaeological site, p. 179-194. In P.J. Julig, ed., The Sheguiandah Site-Archaeological, Geological and Paleobotanical Studies at a Paleoindian Site on Manitoulin Island, Ontario. Canadian Museum of Civilization, Hull-Québec, Archaeological Survey of Canada Paper 161, 314 p.
- Anderson, T.W. and Lewis, C.F.M., 2002. Upper Great Lakes climate and water-level changes 11 to 7 ka: effect on the Sheguiandah archeological site, p. 195-234. In P.J. Julig, ed., The Sheguiandah Site-Archaeological, Geological and Paleobotanical Studies at a Paleoindian Site on Manitoulin Island, Ontario. Canadian Museum of Civilization, Hull-Québec, Archaeological Survey of Canada Paper 161, 314 p.
- Andrews, J.T., 1970. A geomorphological study of postglacial uplift with particular reference to Arctic Canada. Institute of British Geographers, London, Special Publication 2, 156 p.
- Attig, J.W., Clayton, L. and Mickelson, D.M., 1985. Correlation of late Wisconsin glacial phases in the western Great Lakes area. Geological Society of America Bulletin, 96: 1585-1593.
- Baker, R.G., Maher, L.J., Chumbley, C.A. and Van Zant, K.L., 1992. Patterns of change in the midwestern United States. Quaternary Research, 37: 379-389.
- Barnett, P.J., 1979. Glacial Lake Whittlesey: The probable ice frontal position in the eastern end of the Erie basin. Canadian Journal of Earth Sciences, 16: 568-574.
- Barnett, P.J., 1988. History of the northwestern arm of the Champlain Sea, p. 25-36. In N.R. Gadd, ed., The Quaternary Development of the Champlain Sea Basin. Geological Association of Canada, St. John’s, Special Paper 35, 312 p.
- Barnett, P.J., 1992. Quaternary geology of Ontario, p. 1011-1088. In P.C. Thurston, H.R. Williams, R.H. Sutcliffe, and G.M. Scott, ed., Geology of Ontario. Ontario Geological Survey, Sudbury, Special Volume 4, part 2, 1525 p.
- Barnhardt, W.A., Gehrels, W.R., Belknap, D.F. and Kelley, J.T., 1995. Late Quaternary relative sea-level change in the western Gulf of Maine: Evidence for a migrating glacial forebulge. Geology, 23: 317-320.
- Bengtsson, L. and Malm, J., 1997. Using rainfall-runoff modeling to interpret lake level. Journal of Paleolimnology, 18: 235-248.
- Birks, S.J., Edwards, T.W.D. and Remenda, V.H., in press. Isotopic evolution of glacial Lake Agassiz: new insights from cellulose and porewater archives. Palaeogeography, Palaeoclimatology and Palaeoecology.
- Blasco, S.M., 2001. Geological history of Fathom Five National Marine Park over the past 15,000 years, p. 45-62. In S. Parker, and M. Munawar, ed., Ecology, Culture and Conservation of a Protected Area: Fathom Five National Marine Park, Canada. Backhuys Ecovision World Monograph Series, Leiden, 306 p.
- Blasco, S.M., Janusa, S.E., McClellan, S. and Amos, A., 1997. Prehistoric drainage across the submerged Niagara Escarpment north of Tobermory, p. 218-225. In S. Carty, R. Murzin, S. Powell and D. Ramsay, ed., Leading Edge ’97, The Edge and The Point, Niagara Escarpment and Long Point. Conference Proceedings, October 16-18, 1997, Burlington, Ontario, 391 p.
- Bradbury, J.P. and Dean, W.E., ed., 1993. Elk Lake, Minnesota: Evidence for rapid climate change in the north-central United States. Geological Society of America, Boulder, Special Paper 276, 336 p.
- Breckenridge, A., in press. The Lake Superior varve stratigraphy and implications for eastern Lake Agassiz outflow from 10,700 to 8,900 cal BP (9.5-8.0 14C ka). Palaeogeography, Palaeoclimatology and Palaeoecology.
- Breckenridge, A. and Johnson, T.C., 2005. Paleohydrology of Lake Agassiz and the upper Great Lakes from 10,700 to 8,800 cal BP (9,500-7,900 14C BP). Program and abstracts, CANQUA 2005 conference, Canadian Quaternary Association, June 5-8, 2005, Winnipeg, Manitoba, p. A7.
- Breckenridge, A., Johnson, T.C., Beske-Diehl, S., and Mothersill, J.S. 2004. The timing of regional lateglacial events and post-glacial sedimentation rates from Lake Superior. Quaternary Science Reviews, 23: 2355-2367.
- Brooks, G.R., Thorleifson, L.H. and Lewis, C.F.M., 2005. Influence of loss of gradient from postglacial uplift on the Red River flood hazard, Manitoba, Canada. The Holocene, 15: 347-352.
- Bryson, R.A. and Wendlund, W.M., 1967. Tentative climatic patterns for some late glacial and post-glacial episodes in central North America, p. 271-297. In W.J. Mayer-Oakes, ed., Life, Land and Water. University of Manitoba Press, Winnipeg, 414 p.
- Bryson, R.A. and Hare, F.K., 1974. The climates of North America, p. 1-47. In R.A. Bryson and F.K. Hare, ed., Climates of North America, World Survey of Climatology, Volume 11. Elsevier, Amsterdam, 420 p.
- Buhay, W.M. and Betcher, R.N., 1998. Paleohydrologic implications of 18O enriched Lake Agassiz water. Journal of Paleolimnology, 19: 285-296.
- Coordinating Committee on the Great Lakes Basic Hydraulic and Hydrologic Data (CCGLBHHD), 1977. Apparent vertical movement over the Great Lakes. United States Army Corps of Engineers, Chicago, and Environment Canada, Cornwall, 70 p.
- Calkin, P.E., 1970. Strand lines and chronology of the glacial Great Lakes in northwestern New York. Ohio Journal of Science, 70: 78-96.
- Calkin, P.E. and Feenstra, B.H., 1985. Evolution of the Erie-Ontario Great Lakes, p. 149-170. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Canadian Hydrographic Service, 1996. Great Lakes Basin Regional Bathymetry Map. Fisheries and Oceans Canada, Canadian Hydrographic Service Map 880-A, scale 1:1 500 000.
- Chapman, L.J., 1975. The physiography of the Georgian Bay-Ottawa Valley area of southern Ontario. Ontario Geological Survey, Sudbury, Ontario Division of Mines, Geoscience Report 128, 33 p.
- Chapman, L.J. and Putnam, D.F., 1984. The physiography of southern Ontario. Ontario Geological Survey, Sudbury, Special Volume 2, 270 p.
- Chittenden, R., 1990. The origin of Wye Marsh. Wye Marsh Wildlife Centre Publication, Midland, 50 p.
- Chrzastowski, J. and Thompson, T.A., 1992. Late Wisconsinan and Holocene coastal evolution of the southern shore of Lake Michigan, p. 397-413. In J.F. Wehmiller and C.H. Fletcher, ed., Quaternary Coasts of the United States: Marine and Lacustrine Systems. Society for Sedimentary Geology, Special Publication 48, 450 p.
- Chrzastowski, J. and Thompson, T.A., 1994. Late Wisconsinan and Holocene geologic history of the Illinois-Indiana coast of Lake Michigan. Journal of Great Lakes Research, 20: 9-26.
- Chrzastowski, M.J., Pranschke, F.A. and Shabica, C.W., 1991. Discovery and preliminary investigation of the remains of an early Holocene forest on the floor of southern Lake Michigan. Journal of Great Lakes Research, 17: 543-552.
- Clark, J.A., Hendriks, M., Timmermans, T.J., Struck, C. and Hilverda, K.J., 1994. Glacial isostatic deformation of the Great Lakes region. Geological Society of America Bulletin, 106: 19-30.
- Clayton, L., 1983. Chronology of Lake Agassiz drainage to Lake Superior, p. 291-301. In J.T. Teller and L. Clayton, ed., Glacial Lake Agassiz. Geological Association of Canada, St. John’s, Special Paper 26, 451 p.
- Clayton, L. and Attig, J.W., 1989. Glacial Lake Wisconsin. Geological Society of America, Boulder, Memoir 173, 80 p.
- Coakley, J.P. and Lewis, C.F.M., 1985. Postglacial lake levels in the Erie basin, p. 195-212. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Coleman, A.P. 1937. Lake Iroquois, p. 1-36 and map. Ontario Geological Survey, Sudbury, Ontario Department of Mines Annual Report for 1936, vol. 45, part 7, 125 p.
- Colman, S.M., Clark, J.A., Clayton, L., Hansel, A.K. and Larsen, C.E., 1994a. Deglaciation, lake levels and meltwater discharge in the Lake Michigan basin. Quaternary Science Reviews, 13: 879-890.
- Colman, S.M., Forester, R.M., Reynolds, R.L., Sweetkind, D.S., King, J.W., Gangemi, P., Jones, G.A., Keigwin, L.D. and Foster, D.S., 1994b. Lake-level history of Lake Michigan for the past 12 000 years: the record from deep lacustrine sediments. Journal of Great Lakes Research, 20: 73-92.
- Crane, H.R. and Griffin, J.B., 1961. University of Michigan radiocarbon dates VI. Radiocarbon, 3: 105-125.
- Crane, H.R. and Griffin, J.B., 1965. University of Michigan radiocarbon dates X. Radiocarbon, 7: 123-152.
- Crane, H.R. and Griffin, J.B., 1968. University of Michigan radiocarbon dates XII. Radiocarbon, 10: 61-114.
- Crane, H.R. and Griffin, J.B., 1970. University of Michigan radiocarbon dates XIII. Radiocarbon, 12: 161-180.
- Crane, H.R. and Griffin, J.B., 1972. University of Michigan radiocarbon dates XV. Radiocarbon, 14: 195-222.
- David, P.P., 1988. The coeval eolian environment of the Champlain Sea episode, p. 291-305. In N.R. Gadd, ed, The Late Quaternary Development of the Champlain Sea Basin. Geological Association of Canada, St. John’s, Special Paper 35, 312 p.
- Deane, R.E., 1950. Pleistocene geology of the Lake Simcoe district, Ontario. Geological Survey of Canada, Ottawa, Memoir 256, 108 p.
- Dettman, D.L., Smith, A.J., Rea, D.K., Moore, Jr., T.C. and Lohman, K.C., 1995. Glacial meltwater in Lake Huron during early postglacial time as inferred from single-valve analysis of oxygen isotopes in ostracodes. Quaternary Research, 43: 297-310.
- Drexler, C.W., Farrand, W.R. and Hughes, J.D., 1983. Correlation of glacial lakes in the Superior basin with eastward discharge events from Lake Agassiz, p. 309-329. In J.T. Teller and L. Clayton, ed., Glacial Lake Agassiz. Geological Association of Canada, St. John’s, Special Paper 26, 451 p.
- Dyke, A.S., 1996. Preliminary paleogeographic maps of glaciated North America. Geological Survey of Canada, Ottawa, Open File 3296.
- Dyke, A.S., 1998. Holocene delevelling of Devon Island, Arctic Canada: implications for ice sheet geometry and crustal response. Canadian Journal of Earth Sciences, 35: 885-904.
- Dyke, A.S., 2004. An outline of North American deglaciation with emphasis on central and northern Canada, p. 373-424. In J. Ehlers and P.L. Gibbard, ed., Quaternary glaciations–Extent and Chronology, Part II: North America, Developments in Quaternary Science, vol. 2b, Elsevier, Amsterdam, 514 p.
- Dyke, A.S., Moore, A. and Robertson, L., 2003. Deglaciation of North America, Geological Survey of Canada, Ottawa, Open File 1574, 2 map sheets, 1 CD-ROM.
- Dyke, A.S., Giroux, D. and Robertson, L., 2004. Paleovegetation maps of Northern North America, 18 000 to 1 000 BP. Geological Survey of Canada, Ottwawa, Open File 4682, 1 map sheet, 1 CD-ROM.
- Edwards T.W.D., Wolfe B.B., MacDonald G.M., 1996. Influence of changing atmospheric circulation on precipitation 18O–temperature relations in Canada during the Holocene. Quaternary Research, 46: 211-218.
- Eschman, D.F. and Karrow, P.F., 1985. Huron basin glacial lakes: a review, p. 79-93. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Evenson, E.B., 1973. Late Pleistocene shorelines and stratigraphic relations in the Lake Michigan basin. Geological Society of America Bulletin, 84: 2282-2298.
- Farrand, W.R. and Drexler, C.W., 1985. Late Wisconsinan and Holocene history of the Lake Superior basin, p. 17-32. In P.F. Karrow, and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Fitzgerald, W.D., 1985. Postglacial history of the Minesing basin, Ontario, p. 133-146. In P.F. Karrow, and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Fraser, G.S., Larsen, C.E. and Hester, N.C. 1990. Climatic control of lake levels in the Lake Michigan and Lake Huron basins, p. 75-89. In A.F. Schneider and G.S. Fraser, ed., Late Quaternary history of the Lake Michigan basin. Geological Society of America, Boulder, Special Paper 251, 123 p.
- Futyma, R.P., 1981. The northern limits of Lake Algonquin in Upper Michigan. Quaternary Research, 15: 291-310.
- Futyma, R.P. and Miller, N.G., 1986. Stratigraphy and genesis of the Lake Sixteen peatland, northern Michigan. Canadian Journal of Botany, 64: 3008-3019.
- Goldthwait, J.W., 1907. Abandoned shorelines of eastern Wisconsin. Wisconsin Geological and Natural History Survey, Bulletin 17, 134 p.
- Goldthwait, J.W., 1910. An instrumental survey of the shorelines of extinct lakes in southwestern Ontario. Geological Survey of Canada, Ottawa, Memoir 10, 57 p.
- Hansel, A.K., Mickelson, D.M., Schneider, A.F. and Larsen, C.E., 1985. Late Wisconsinan and Holocene history of the Lake Michigan basin, p. 39-53. In P.F. Karrow, and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Harrison, J.E., 1972. Quaternary geology of the North Bay-Mattawa region, Ontario. Geological Survey of Canada, Ottawa, Paper 71-26, 37 p.
- Hough, J.L., 1955. Lake Chippewa, a low stage of Lake Michigan indicated by bottom sediments. Geological Society of America Bulletin, 66: 957-968.
- Hough, J.L., 1958. Geology of the Great Lakes. University of Illinois Press, Urbana, 313 p.
- Hough, J.L., 1962. Lake Stanley, a low stage of Lake Huron indicated by bottom sediments. Geological Society of America Bulletin, 98: 113-120.
- Julig, P.J. and Mahaney, W.C., 2002. Geoarchaeological studies of the Sheguiandah Site and analysis of museum collections, p. 101-137. In P.J. Julig, ed., The Sheguiandah Site-Archaeological, Geological and Paleobotanical Studies at a Paleoindian Site on Manitoulin Island, Ontario. Canadian Museum of Civilization, Hull-Québec, Archaeological Survey of Canada Paper 161, 314 p.
- Karrow, P.F., 1980. The Nipissing transgression around southern Lake Huron. Canadian Journal of Earth Sciences, 17: 1271-1279.
- Karrow, P.F., 1986. Valley terraces and Huron basin water levels, southwestern Ontario. Geological Society of America Bulletin, 97: 1089-1097.
- Karrow, P.F., 1987. Glacial and glaciolacustrine events in northwestern Lake Huron, Michigan and Ontario. Geological Society of America Bulletin, 98: 113-120.
- Karrow, P.F., Anderson, T.W., Clarke, A.H., Delorme, L.D. and Sreenivasa, M.R., 1975. Stratigraphy, paleontology, and age of Lake Algonquin sediments in southwestern Ontario, Canada. Quaternary Research, 5: 49-87.
- Kaszycki, C.A., 1985. History of glacial Lake Algonquin in the Haliburton region, south central Ontario, p. 109-123. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Kehew, A.E., 1993. Glacial-lake outburst erosion of the Grand Valley, Michigan, and impacts on glacial lakes in the Lake Michigan basin. Quaternary Research, 39: 36-44.
- King, J. and Peck, J., 2001. Use of paleomagnetism in studies of lake sediments, p. 371-389. In W.M. Last and J.P. Smol, ed., Tracking Environmental Change Using Lake Sediments. Volume 1: Basin Analysis, Coring, and Chronological Techniques. Kluwer Academic Publishers, Dordrecht, 548 p.
- Larsen, C.E., 1987. Geological history of Glacial Lake Algonquin and the Upper Great Lakes. United States Geological Survey, Reston, Bulletin 1801, 36 p.
- Larson, D.W. and Kelly, P.E., 1994. Holocene remnants of Niagara escarpment Thuja Occidentalis forests: Ecological, hydrological and geological implications. Geological Association of Canada, Program with Abstracts, 19: A61.
- Larson, G. and Schaetzl, R., 2001. Origin and evolution of the Great Lakes. Journal of Great Lakes Research, 27: 518-546.
- Leverett, F. and Taylor, F.B., 1915. The Pleistocene of Indiana and Michigan and the history of the Great Lakes. United States Geological Survey, Reston, Monograph 53, 529 p.
- Leverington, D.W., Teller, J.T. and Mann, J.D., 2002. A GIS method for reconstruction of late Quaternary landscapes from isobase data and modern topography. Computers and Geosciences, 28: 631-639.
- Lewis, C.F.M., 1969. Late Quaternary history of lake levels in the Huron and Erie basins, p. 250-270. Proceedings of the 12th Conference Great Lakes Research, May 5-7, 1969, University of Michigan, Ann Arbor, Michigan. International Association of Great Lakes Research, Ann Arbor, 923 p.
- Lewis, C.F.M., 1971. Recent uplift of Manitoulin Island, Ontario. Canadian Journal of Earth Sciences, 7: 665-675.
- Lewis, C.F.M. and Anderson, T.W., 1989. Oscillations of levels and cool phases of the Laurentian Great Lakes caused by inflows from glacial Lakes Agassiz and Barlow-Ojibway. Journal of Paleolimnology, 2: 99-146.
- Lewis, C.F.M., Moore, Jr., T.C., Rea, D.K., Dettman, D.L., Smith, A.J. and Mayer, L.A., 1994. Lakes of the Huron basin: their record of runoff from the Laurentide Ice Sheet. Quaternary Science Reviews, 13: 891-922.
- Lewis, C.F.M., Todd, B.J., Forbes, D.L., Nielsen, E. and Thorleifson, L.H., 2000. Architecture, age and lithology of sediments in Lake Winnipeg: seismostratigraphy, long-core lithostratigraphy, and basin evolution, p. 127-188. In B.J. Todd, C.F.M. Lewis, D.L. Forbes, L.H. Thorleifson and E. Nielsen, ed., 1996 Lake Winnipeg project: cruise report and scientific results. Geological Survey of Canada, Ottawa, Open File 3070, 842 p.
- Lewis, C.F.M. Coakley, J.P., King, J.W., Heil, C.W. Jr., Hubeny, J.B., Cameron, G.D.M., and McCarthy, F.M.G., 2004. A submerged beach and shore face below eastern Lake Erie implies early Holocene climate-driven, closed-lake conditions. Joint Annual Meeting of Geological Association of Canada and Mineralogical Association of Canada, May 12-14, 2004, St. Catharines, Ontario, Abstracts, 29: 191.
- Lewis, C.F.M., Heil, Jr., C.W., Hubeny, J.B., King, J.W., Moore, Jr., T.C. and Rea, D.K., in press. The Stanley unconformity in Lake Huron basin, evidence for a climate-driven closed lowstand about 7900 14C BP, with similar implications for the Chippewa lowstand in Lake Michigan basin. Journal of Paleolimnology.
- Lowdon, J.A. and Blake, Jr., W., 1968. Geological Survey of Canada radiocarbon dates VII. Radiocarbon, 10: 207-245.
- Lowdon, J.A. and Blake, Jr., W., 1975. Geological Survey of Canada radiocarbon dates XV. Geological Survey of Canada, Ottawa, Paper 75-7, 32 p.
- Lowdon, J.A., Fyles, J.G. and Blake, Jr., W., 1967. Geological Survey of Canada radiocarbon dates VI. Geological Survey of Canada, Ottawa, Paper 67-2, Part B, 42 p.
- Lowdon, J.A., Robertson, I.M. and Blake, Jr., W., 1971. Geological Survey of Canada radiocarbon dates XI, reprinted from Radiocarbon, 13: 255-324. Geological Survey of Canada, Ottawa, Paper 71-7, 70 p.
- Lowell, T.V., Larson, G.J., Hughes, J.D. and Denton, G.H., 1999. Age verification of the Lake Gribben forest bed and the Younger Dryas Advance of the Laurentide Ice Sheet. Canadian Journal of Earth Sciences, 36: 383-393.
- Mainville, A. and Craymer, M., 2005. Present-day tilting of the Great Lakes region based on water level gauges. Geological Society of America Bulletin, 117: 1070-1080.
- Monaghan, G.W., Lovis, W.A. and Fay, L., 1986. The Lake Nipissing transgression in the Saginaw Bay region, Michigan. Canadian Journal of Earth Sciences, 23: 1851-1854.
- Moore, Jr., T.C, Rea, D.K., Mayer, L.A., Lewis, C.F.M. and Dobson, D.M., 1994. Seismic stratigraphy of Lake Huron-Georgian Bay and post-glacial lake level history. Canadian Journal of Earth Sciences, 31: 1606-1617.
- Moore, Jr., T.C., Walker, J.C.G., Rea, D.K., Lewis, C.F.M., Shane, L.C.K. and Smith, A.J., 2000. The Younger Dryas interval and outflow from the Laurentide Ice Sheet. Paleoceanography, 15: 4-18.
- Mortsch, L., Hengeveld, H., Lister, M., Lofgren, B., Quinn, F.H., Slivitsky, M. and Wenger, L., 2000. Climate change impacts on the hydrology of the Great Lakes-St. Lawrence system. Canadian Water Resources Journal, 25: 153-179.
- Muller, E.H. and Prest, V.K., 1985. Glacial lakes in the Ontario basin, p. 213-229. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- National Atlas of Canada, 1985. Canada Drainage Basins, Map 5.1. The National Atlas of Canada, Energy Mines and Resources Canada, Ottawa, Map MCR 4055, scale 1:7 500 000.
- Pair, D.L. and Rodrigues, C.G., 1993. Late Quaternary deglaciation of the southwestern St. Lawrence Lowland, New York and Ontario. Geological Society of America Bulletin, 105: 1151-1164.
- Pair, D., Karrow, P.F. and Clark, P.U. 1988. History of the Champlain Sea in the central St. Lawrence lowland, New York, and its relationship to waterlevels in the Lake Ontario basin, p. 107-123. In N.R. Gadd, ed., The Late Quaternary Development of the Champlain Sea Basin. Geological Association of Canada, St. John’s, Special Paper 35, 312 p.
- Peltier, W.R., 1994. Ice age paleotopography. Science, 265: 195-201.
- Peltier, W.R., 1995. VLBI baseline variations for the ICE‑4G model of postglacial rebound. Geophysical Research Letters, 22: 465-468.
- Peltier, W.R., 1996. Mantle viscosity and ice-age ice sheet topography. Science, 273: 1359-1364.
- Peltier, W.R., 1998. Postglacial variations in the level of the sea: Implications for climate dynamics and solid-earth geophysics. Reviews of Geophysics, 36: 603-689.
- Peterson, R.L., 1965. A well-preserved grizzly bear skull recovered from a late glacial deposit near Lake Simcoe, Ontario. Nature, 208: 1233-1234.
- Rea, D.K., Moore, Jr., T.C., Lewis, C.F.M., Mayer, L.A., Dettman, D.J., Smith, A.J. and Dobson, D.M., 1994a. Stratigraphy and paleolimnologic record of lower Holocene sediments in northern Lake Huron and Georgian Bay. Canadian Journal of Earth Sciences, 31: 1586-1605.
- Rea, D.K., Moore, Jr., T.C., Anderson, T.W., Lewis, C.F.M., Dobson, D.M., Dettman, D.L., Smith, A.J. and Mayer, L.A., 1994b. Great Lakes paleohydrology: Complex interplay of glacial lake levels and sill depths. Geology, 22: 1059-1602.
- Saarnisto, M., 1974. The deglaciation history of the Lake Superior region and its climatic implications. Quaternary Research, 4: 316-339.
- Schaetzl, R.J., Drzyzga, S.A., Weisenborn, B.N., Kincare, K.A., Lepczyk, X.C., Shein, K., Dowd, C.M. and Linker, J., 2002. Measurement, correlation, and mapping of glacial Lake Algonquin shorelines in northern Michigan. Annals of the Association of American Geographers, 92: 399-415.
- Shaw, J., Gareau, P. and Courtney, R.C., 2002. Palaeogeography of Atlantic Canada 13-0 kyr. Quaternary Science Reviews, 21: 1861-1878.
- Simard, J., Occhietti, S. and Robert, F. 2003. Retrait de l’inlandsis sur les Laurentides au début de l’Holocene: Transect de 60 km entre le Saint-Maurice et le Témiscamingue (Québec). Géographie physique et Quaternaire, 57: 189-204.
- Sly, P.G. and Lewis, C.F.M., 1972. The Great Lakes of Canada–Quaternary Geology and Limnology. 24th International Geological Congress, Montréal, Guidebook for Field Excursion A43, 92 p.
- Sly, P.G. and Sandilands, R.G., 1988. Geology and environmental significance of sediment distributions in an area of the submerged Niagara escarpment, Georgian Bay. Hydrobiologica, 163: 47-76.
- Stanley, G.M., 1936. Lower Algonquin beaches of Penetanguishene Peninsula. Geological Society of America Bulletin, 47: 1933-1959.
- Stanley, G.M., 1937. Lower Algonquin beaches of Cape Rich, Georgian Bay. Geological Society of America Bulletin, 48: 1665-1686.
- Stanley, G.M., 1938a. Impounded early Algonquin beaches at Sucker Creek, Grey County, Ontario. Michigan Academy of Science, Arts, and Letters Papers, 23: 477-495.
- Stanley, G.M., 1938b. The submerged valley through Mackinac Straits. Journal of Geology, 46: 966-974.
- Stuiver, M. and Reimer, P.J., 1993., Extended 14C database and revised CALIB radiocarbon calibration program. Radiocarbon, 35: 315-230.
- Stuiver, M., Reimer, P.J., Bard, E., Beck, J.W., Burr, G.S., Hughen, K.A., Kromer, B., McCormac, F.G., Plicht, J. and Spurk, M. 1998. INTCAL98 radiocarbon age calibration 24,000–0 cal BP. Radiocarbon, 40: 1041-1083.
- Taylor, L.D., 1990. Evidence of high glacial-lake levels in the northeastern Lake Michigan basin and their relation to the Glenwood and Calumet phases of glacial Lake Chicago, p. 91-109. In A.F. Schneider and G.S. Fraser, ed., Late Quaternary history of the Lake Michigan basin. Geological Society of America, Boulder, Special Paper 251, 123 p.
- Teller, J.T., 1985. Glacial Lake Agassiz and its influence on the Great Lakes, p. 1-16. In P.F. Karrow and P.E. Calkin, ed., Quaternary Evolution of the Great Lakes. Geological Association of Canada, St. John’s, Special Paper 30, 258 p.
- Teller, J.T. and Thorleifson, L.H., 1983. The Lake Agassiz-Lake Superior connection, p. 261-290. In J.T. Teller and L. Clayton, ed., Glacial Lake Agassiz. Geological Association of Canada, St. John’s, Special Paper 26, 451 p.
- Teller, J.T, Leverington, D.W. and Mann, J.D., 2002. Freshwater outbursts to the oceans from glacial Lake Agassiz and their role in climate change during the last deglaciation. Quaternary Science Reviews, 21: 879-887.
- Teller, J.T. and Leverington, D.W. 2004. Glacial Lake Agassiz: A 5000 yr history of change and its relationship to the 18O record of Greenland. Geological Society of America Bulletin, 116: 729-742.
- The Great Lakes Environmental Atlas, 1995. The Great Lakes, An Environmental Atlas and Resource Book, Third edition. Government of Canada, Toronto and United States Environmental Protection Agency, Chicago, 46 p.
- Thorleifson, L.H., 1996. Review of Lake Agassiz history, p. 55-84. In J.T. Teller, L.H. Thorleifson, G.Matile, and W.C. Brisbin, Compilers, Sedimentology, geomorphology and history of the central Lake Agassiz basin (Field Trip B2), Geological Association of Canada/Mineralogical Association of Canada Annual Meeting, Winnipeg, Manitoba, May 27-29, 1996, 101 p.
- Thorleifson, L.H. and Kristjansson, F.J., 1993. Quaternary geology and drift prospecting, Beardmore-Geraldton area, Ontario. Geological Survey of Canada, Ottawa, Memoir 435, 146 p.
- Tovell, W.M and Deane, R.E., 1966. Grizzly bear skull: site of a find near Lake Simcoe, Ontario. Science, 154: 158.
- Tushingham, A.M., 1992. Postglacial uplift predictions and historical water levels of the Great Lakes. Journal of Great Lakes, 18: 440-455.
- Tushingham, A.M. and Peltier, W.R., 1991. ICE‑3G: A new global model of late Pleistocene deglaciation based upon geophysical predictions of postglacial relative sea level. Journal of Geophysical Research, 96: 4497-4523.
- United States Lake Survey, 1965. General Chart of the Great Lakes. United States Corps of Engineers, United States Army, United States Lake Survey Chart No 0. Scale 1:1 200 000.
- Veillette, J.J., 1994. Evolution and paleohydrology of glacial Lakes Barlow and Ojibway. Quaternary Science Reviews, 13: 944-971.
- Vincent, J.-S. and Hardy, L., 1979. The evolution of glacial Lakes Barlow and Ojibway, Québec and Ontario. Geological Survey of Canada, Ottawa, Bulletin 316, 18 p.
- Walcott, R.I., 1972. Late Quaternary vertical movements in eastern North America: quantitative evidence of glacio-isostatic rebound. Reviews of Geophysics and Space Physics, 10: 849-884.
- Warner, B.G., Hebda, R.J. and Hann, B.J., 1984. Postglacial paleoecological history of a cedar swamp, Manitoulin Island, Ontario, Canada. Palaeogeography, Palaeoclimatology and Palaeoecology, 45: 301-345.
- Webb III, T., Anderson, K.H., Bartlein, P.J. and Webb, R.S., 1998. Late Quaternary climate change in eastern North America: a comparison of pollen-derived estimates with climate model results. Quaternary Science Reviews, 17: 587-606.
- Wolfe, S.A., Huntley, D.S. and Ollerhead, J., 2004. Relict late Wisconsinan dune fields of the northern Great Plains, Canada. Géographie physique et Quaternaire, 58, 323-336.
- Woodend, S.L., 1983. Glacial and post-glacial history of Lake Huron as defined by lithological and palynological analysis of a core in southern Lake Huron with special regard to a low-level phase. B.Sc. thesis, Carleton University, 74 p.
- Wright, Jr., H.E., Stefanova, I., Tian, J., Brown, T.A. and Hu, F.S., 2004. A chronological framework for the Holocene vegetational history of central Minnesota: the Steel Lake pollen record. Quaternary Science Reviews, 23: 611-626.
- Yu, Z. and Wright, Jr., H.E., 2001. Response of interior North America to abrupt climate oscillations in the North Atlantic regions during the last deglaciation. Earth Science Reviews, 52: 333-369.
List of figures
Figure 1
Shaded-relief map of the present-day Laurentian Great Lakes with their drainage areas outlined in white from the National Atlas of Canada (1985). The map file is of 30 arc-second resolution, and is from the Canadian Hydrographic Service (1996) with bathymetry added for Georgian Bay and the area west of LaCloche Lowland from United States Lake Survey (1965). Numbers refer to lake-level indicator sites listed in Table II. Lake outlets shown by black arrows.
Carte du relief de la région des Grands Lacs et des limites de ses sous-bassins de drainage (en blanc) établie à partir de l’Atlas national du Canada (1985). La base de données possède une résolution de 30 secondes d’arc et provient du Service hydrographique du Canada (1996), où la bathymétrie ajoutée de la baie Géorgienne et du secteur à l’ouest des basses terres de LaCloche provient de la Commission des Lacs des États-Unis (1965). Les chiffres renvoient aux sites d’indicateurs du niveau des lacs présentés dans le tableau II. Les flèches noires localisent l’exutoire des lacs.
A
B
Maps of isobases of selected reference paleoshorelines in basins of the Great Lakes, showing the radiocarbon ages (ka BP) and trends of isobases with their present lowest and highest elevations. Sources of isobases given in Table I. C and D mark end points of a section to which isobases of the Iroquois, Algonquin and Nipissing lakes were projected for illustration in Figure 3. (A) Map of Laurentide maximum ice margin and isobases of lakes Wisconsin, Washburn, Minong, Algonquin, Glenwood II, Whittlesey and Iroquois. (B) Map of isobases for the Nipissing Great Lakes, and for an older phase of glacial Lake Wisconsin.
Cartes des isobases des anciennes lignes de rivage de référence de la région des Grands Lacs, où les âges au 14C (ka BP) et le soulèvement différentiel par rapport aux altitudes actuelles sont illustrées. Les sources bibliographiques sont présentées au tableau I. Les points C et D sont les limites du profil vertical projeté à la figure 3 des lacs Iroquois, Algonquin et Nippising. (A) Carte de l’extension glaciaire maximale de l’inlandsis laurentidien et des isobases des lacs Wisconsin, Washburn, Minong, Algonquin, Glenwood II, Whittlesey et Iroquois. (B) Carte des isobases de la phase Nipissing des Grands Lacs et pour une phase plus ancienne du lac Wisconsin.
Figure 3
A plot of relative uplift vs. age for the Nipissing, Algonquin and Iroquois strandlines showing a relatively good fit to the exponential uplift equation for relaxation times (τ) between 3000 and 5000 years. The relative uplifts are the differences in projected strandline elevations at C and D shown on Figure 2.
Graphique temporel du soulèvement relatif des lignes de rivage des lacs Nipissing, Algonquin et Iroquois montrant un très bon ajustement à l’équation exponentielle de soulèvement calculée pour des temps de relaxation (τ) de 3000 à 5000 ans. Les soulèvements relatifs correspondent aux différences entre les altitudes projetées entre les points C et D de la figure 2.
Figure 4
Map of the Great Lakes basin showing locations of transects where isobases of two shorelines of known age (shown as 14C ka BP) provide a basis for computing the relaxation time of the uplifting process. The mean value of these estimates, 3700 ± 700 years, is used as the relaxation time throughout the basin.
Carte du bassin des Grands Lacs montrant l’emplacement des transects, où les isobases de deux lignes de rivage d’âge connu (montrés en 14C ka BP) ont servi comme base de calcul du temps de relaxation du processus de soulèvement isostatique. La valeur moyenne de ces estimations, 3700 ± 700 ans, est utilisée comme temps de relaxation sur l’ensemble du bassin.
Figure 5
Michigan-east Superior profile showing construction of the gradients of shorelines representing uplift since the reference age of 10.6 ka BP (12.6 cal ka BP) relative to area B beyond the limit of the last glaciation. Dots are observed shoreline isobases. Dotted, dashed and solid lines are shore gradients adjusted to an age of 10.6 ka BP (12.6 cal ka BP). Inset map shows locations of reference profiles for contouring the reference response surface.
Profil du Michigan-est Supérieur qui montre le soulèvement des lignes de rivage depuis l’âge de référence de 10,6 ka BP (12,6 cal ka BP) par rapport à la région B située au-delà de la limite de la dernière glaciation. Les points sont les isobases des lignes de rivages observées, les lignes pointillées, les lignes tiretées et les lignes continues représentent les gradients de rivage ajustés à un âge de 10,6 ka BP (12,6 cal ka BP). L’insertion montre l’emplacement des profils utilisés pour délimiter la surface de référence de la réponse isostatique.
Figure 6
The reference isostatic response surface portraying glacial rebound since 10.6 ka BP (12.6 cal ka BP) in the Great Lakes region relative to an area southwest of Lake Michigan beyond the limit of the last glaciation.
Surface de référence de la réponse isostatique montrant le relèvement depuis 10,6 ka BP (12,6 cal ka BP) pour la région des Grands Lacs par rapport à un secteur au sud-ouest du lac Michigan se situant au-delà de la limite maximale de la dernière glaciation.
A
B
Paleogeographic maps showing reconstructions of the bathymetry and topography of the Great Lakes basin. (A) Main Lake Algonquin phase about 10.55 ± 0.1 ka BP (12.4-12.7 cal ka BP). At a late stage of Main Lake Algonquin, shown here, the lake had expanded to a possible maximum area in Superior (dashed line) and Huron basins by calving of icebergs from glacier margins in deep water. This lake overflowed via Port Huron to the Erie basin. Early Lake Erie, controlled by the Lyell-Johnson sill near the present Niagara Falls, overflowed via Niagara River to Early Lake Ontario which discharged via the emerging St. Lawrence River to Champlain Sea. (B) Mattawa highstand phase about 8.7-8.8 ka BP (9.6-9.9 cal ka BP), possibly a phase of high discharge when lake level in the upper Great Lakes basins was hydraulically dammed by constrictions downstream from North Bay at either the Rutherglen Moraine or the Rankin Constriction or both (Fig. 1). Discharge continued via Ottawa River to St. Lawrence River. There was no drainage from the upper to lower Great Lakes at this time; Lake Erie, Niagara River and Lake Ontario drained a separate watershed to the upper St. Lawrence River. Legend as for Figure 7A.
Cartes paléo-géographiques de la reconstitution de la bathymétrie et de la topographie du bassin des Grands Lacs. (A) Phase du lac Algonquin principal en 10,55 ± 0,1 ka BP (12,4-12,7 cal ka BP). Au cours du stade tardif de cette phase, le lac Algonquin a atteint son extension maximale dans les bassins du lac Supérieur (tireté) et Huron par le vêlage des icebergs issus des marges glaciaires en eau profonde. Ce lac se déversait dans le bassin du lac Érié par le Port Huron. Le lac Érié initial, dont le niveau était contrôlé par le seuil de Lyell-Johnson situé près des chutes actuelles du Niagara, se déversait dans le lac Ontario initial par la rivière Niagara, se déversant à son tour (le lac Ontario) dans le Saint-Laurent naissant par la Mer de Champlain. (B) Phase de haut niveau lacustre de Mattawa en 8,7-8,8 ka BP (9,6-9,9 cal ka BP), une phase de fort débit probablement liée à la fermeture hydraulique des Grands Lacs par un étranglement topographique situé en aval de North Bay provoqué par la moraine de Rutherglen, l’étranglement de Rankin, ou les deux. Le déversement dans le Saint-Laurent s’effectuait par la rivière des Outaouais. Il n’y avait pas de drainage entre l’amont et l’aval des Grands Lacs à cette époque : les lacs Érié et Ontario et la rivière Niagara drainaient un bassin-versant distinct vers le haut Saint-Laurent. La légende est la même qu’en A.
Figure 8
Bar graphs showing variation in the land drainage areas, lake surface areas, and lake volumes of selected Great Lakes overflowing phases measured from reconstructions of the paleo-Great Lakes in a GIS computing environment. Triangles indicate adjacent lakes that were and are confluent. Names at right refer to lake phases in Huron and other upper Great Lakes basins.
Diagrammes de la variation des superficies de drainage, des surfaces et des volumes des Grands Lacs au cours des différentes phases de déversement obtenues à partir des reconstitutions des paléo-Grands Lacs par un système d’information géographique (SIG). Les triangles indiquent s’il y a des lacs adjacents ayant ou encore connectés. Les noms à la droite font référence aux phases lacustres du bassin du lac Huron et des autres Grands Lacs.
A
B
C
D
E
F
G
H
Types of evidence of former water levels in the Huron, Michigan, and Erie basins. Abandoned geomorphic shore features: (A) Former lake surfaces are inferred from erosional coastal features at the change in gradient between the gently-sloping shoreface and the steeply-sloping shorebluff, illustrated here for the Algonquin shore near Kirkfield, Ontario. (B) Storm beach ridges above Lake Huron on Manitoulin Island (air photo courtesy of Natural Resources Canada). (C) and (D) Beach ridges 53 m below Georgian Bay imaged by sidescan sonar and multibeam sonar, respectively (Blasco 2001). (E) Submerged beach and shoreface 21-28 m below Lake Erie (from Coakley and Lewis, 1985). Isolation basins, tree stumps, and unconformities: (F) Air photo (courtesy of Natural Resources Canada) view of a small basin isolated during the regression of a large lake, Manitoulin Island. A radiocarbon date of the contact zone between gray clastic large lake sediment and dark-coloured organic small-lake-sediment (photo courtesy of T.W. Anderson) provides the age at which the Great Lake surface passed below the sill elevation of the small isolated basin (Lewis, 1971). (G) In situ tree stump on lakefloor at entrance to Georgian Bay. (H) Subsurface reflections in a seismic profile indicate offshore unconformities and sequence boundaries caused by episodes of reduced lake level (Moore et al., 1994).
Indicateurs des anciens niveaux lacustres dans les bassins des lacs Huron, Michigan et Érié. Formes littorales abandonnées : (A) les surfaces des anciens lacs sont estimées à partir de formes d’érosion côtière qui sont situées à la rupture de pente entre la face de la côte en pente douce et le talus littoral en pente abrupte, comme le montre cet exemple du rivage du lac Algonquin près de Kirkfield, Ontario. (B) Crêtes de plage de tempête en bordure du lac Huron sur l’Île Manitoulin (courtoisie de Ressources Naturelles Canada). (C) et (D) Crêtes de plage à 53 m de profondeur dans la baie Géorgienne identifiées respectivement à l’aide du sonar latéral et d’un sonar multi-faisceaux (Blasco, 2001). (E) Plage et avant-côte submergées entre 21 et 28 m sous la surface du lac Érié (d’après Coakley et Lewis, 1985). Bassins isolés, souches d’arbres et discontinuités : (F) Photo aérienne d’un petit bassin isolé durant la régression d’un grand lac, Île Manitoulin (Source : Ressources Naturelles Canada). Une datation au radiocarbone de la zone de contact entre les sédiments clastiques gris du grand lac et les sédiments organiques foncés du petit lac (photo fournie par T.W. Anderson) fournit l’âge auquel le plan d’eau du grand lac est passé sous le niveau du seuil du petit bassin isolé (Lewis, 1971). (G) Souche d’arbre in situ au fond du lac à l’entrée de la baie Géorgienne. (H) Profil sismique montrant les discontinuités et les limites causées par des épisodes de bas niveau lacustre (Moore et al., 1994).
A
Plot of original elevations of numbered water-level indicators (Table II) in Huron and Georgian Bay basins relative to area B (Fig. 5), inferred lake surface, and potential overflow outlet elevations (Table III) vs. age in calibrated and radiocarbon years BP. Downward-pointing triangles indicate data implying a lake level below the plotted site, as for in situ tree stumps. Upward-pointing triangles are for data implying a lake level above the plotted site, as for shallow water mollusks or driftwood. Right-pointing triangles indicate data implying isolation of a small basin by a falling large lake level. Left-pointing triangles indicate data implying a transgression by a large lake. See Figure 10B for a plot of all data in the 10 000-7500 cal BP interval. The inferred lake level for Huron basin is shown by the thick continuous line (dotted for Georgian Bay basin where different). The elevation history of potential overflow outlets (Table III) is portrayed by a dotted band (North Bay–NB) or by dashed curves (Huron basin sills–H, Georgian Bay basin sill–G, and Michigan basin–M). The Algonquin lake phases are Kirkfield Algonquin (K), Main Algonquin (M), and Post-Algonquin (P). MS and ES are middle and early Stanley lowstands in Huron basin, respectively. In Georgian Bay basin, the equivalent lowstands are named Hough. The horizontal lines above the X-axis, named by colours, indicate age and duration of lowstand sequence boundaries inferred from seismic reflection and core data from Moore et al. (1994) and Rea et al. (1994a).
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres (tabl. II) pour les bassins du lac Huron et de la baie Géorgienne relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires (tabl. III), l’âge en années étalonnées et conventionnelles. Les triangles pointant vers le bas font référence à un niveau lacustre inférieur au site reporté, comme dans le cas des souches d’arbres in situ. Les triangles pointant vers le haut font référence à un niveau lacustre supérieur au site reporté, comme dans le cas des mollusques en eau peu profonde ou de bois flottant. Les triangles pointant vers la droite font référence à l’isolement d’un petit bassin suite à baisse du niveau d’un grand lac. Les triangles pointant vers la gauche font référence à la transgression d’un grand lac. Voir la figure 10B pour le diagramme des données comprises dans l’intervalle 10 000-7500 cal BP. Le niveau lacustre reconstitué du lac Huron est représenté par un trait gras continu (en discontinu pour la baie Géorgienne lorsqu’il diffère). L’historique de l’altitude des exutoires potentiels (tabl. III) est représenté par une bande en pointillé (North Bay–NB) ou par les courbes en tireté (seuils du bassin du lac Huron–H, seuil de la baie Géorgienne–G et bassin du lac Michigan–M). Les phases du Lac Algonquin sont les suivantes : Kirkfield Algonquin (K), Algonquin principal (M) et Post-Algonquin (P). MS et ES identifient respectivement les niveaux moyen et ancien de Stanley du lac Huron. Dans le bassin de la baie Géorgienne, les niveaux équivalents sont nommés Hough. Les lignes horizontales au-dessus de l’abscisse, nommées par une couleur, indiquent l’âge et la durée des limites déduites à partir de données de réflexion sismique et de forage issues des travaux de Moore et al. (1994) et de Rea et al. (1994a).
B
Expanded plot of original elevations of water-level indicators for the 10 000-7500 cal BP period in Huron and Georgian Bay basins relative to area B (Fig. 5), inferred lake surface, and potential overflow outlets vs. age in calibrated and radiocarbon years BP. Symbols as described in Figure 10A.
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres pour la période 10 000-7 500 cal BP pour les bassins du lac Huron et de la baie Géorgienne relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires, l’âge en années étalonnées et conventionnelles. Les symboles sont les mêmes qu’en A.
Figure 11
Michigan basin plot of original elevations of water-level indicators relative to area B (Fig. 5), inferred lake surface, and potential overflow outlets vs. age in calibrated years BP and radiocarbon age in ka BP. Water-level indicators for the Huron basin above the Michigan sill are applied to the Michigan basin as these highstand lakes were confluent in both basins. Black symbols indicate data specific to the Michigan basin. MC and EC represent middle and early Chippewa lowstands, respectively. Otherwise symbols and lines as for Figure 10A.
Diagramme de l’altitude initiale des indicateurs des anciens niveaux lacustres pour le bassin du lac Michigan relativement à la région B (Fig. 5), où les points donnent l’altitude des exutoires (tabl. III), l’âge en années étalonnées et conventionnelles. Les indicateurs du niveau lacustre du lac Huron situés au-dessus du seuil Michigan sont appliqués au lac Michigan puisque ces plans d’eau élevés étaient connectés aux deux bassins. Les symboles noirs font référence aux données spécifiques du lac Michigan. MC et EC représentent les niveaux moyen et ancien de Chippewa. Les autres symboles et les lignes sont identiques à ceux de la figure 10A.
Figure 12
Paleogeographic lowstand map of Huron, Erie and Ontario basins at about 7.8 ka BP (8.6 cal ka BP) showing the reduced area of inferred closed (terminal) lakes whose shorelines are well offshore from the present lake shores. E = Erie basin, GB = Georgian Bay basin, H = Huron basin, and O = Ontario basin.
Carte paléogéographique du bas niveau lacustre à environ 7,8 ka BP (8,6 cal ka BP) pour les bassins des lacs Huron, Érié et Ontario. La carte montre la superficie réduite des lacs fermés dont les anciens rivages étaient bien au large des rivages actuels. E = bassin du lac Érié, GB = bassin de la baie Géorgienne, H = bassin du lac Huron et O = bassin du lac Ontario.
Figure 13
Deviations of sill altitudes when they should have been of equal elevation during drainage transfers from northern to southern outlets of the Algonquin and Nipissing lake phases for the Clark et al. (1994) and ICE‑4G geophysical models and the empirical model used in this paper. Sill elevations, reduced relative to Port Huron at the southern end of the Huron basin, show that the empirical uplift model used in this study best describes (lowest deviations) the relative elevations of outlets during drainage transfers.
Écarts entre les altitudes des seuils lorsqu’elles devraient être égales aux autres durant les transferts du drainage des exutoires septentrionaux aux exutoires méridionaux des phases lacustres Algonquin et Nipissing pour les modèles géophysiques de Clark et al. (1994) et de ICE‑4G et le modèle empirique utilisé dans cette étude. Les altitudes des seuils, ajustées à celles du Port Huron à l’extrémité sud du bassin du lac Huron, montrent que le modèle empirique de soulèvement glacio-isostatique utilisé dans cette étude rend mieux compte (écarts plus faibles) des altitudes relatives des seuils durant les transferts de drainage.
List of tables
Table I
Sources of isobase and deglacial information
Table II
Indicators of former lake levels in the Michigan, Huron and Georgian Bay basins, and Georgian Bay-North Bay lowland
Surface elevation (m asl) of sampled lake or bog.
Present elevation (m asl) of feature representing former level of a Great Lake.
Uplift (m) since 10.6 ka BP (12.6 cal ka BP) from reference isostatic surface (Fig. 6).
Elevations (m asl) computed using equation (9) for the age limits of the calibrated age range of the lake-level indicator as listed in column 7.
Ages at limits of the calibrated range (1σ) of the lake-level indicator date and its uncertainty reported in column 3. Calibration by the Calib 5.0.1 program
(Stuiver and Reimer, 1993) using the INTCAL98 calibration dataset (Stuiver et al., 1998).
Shell age uncorrected for reservoir effect.
n.r. = not reported.
PSV = Age by comparing paleomagnetic secular variation with a 14C-dated standard for northeastern U.S.A. (King and Peck, 2001).
Table III
Present and early Holocene elevations1 of North Bay, Georgian Bay, Huron and Michigan sills (m above present sea level)