12 600 years of lake level changes, changing sills, ephemeral lakes and Niagara Gorge erosion in the Niagara Peninsula and Eastern Lake Erie basin

Over the last 12600 years, lake levels in the eastern Lake Erie basin have fluctuated dramatically, causing major changes in drainage patterns, flooding and draining ephemeral Lake Wainfleet several times and widening and narrowing the Niagara Gorge as the erosive effects of Niagara Falls waxed and waned. The control sill for Lake Erie levels was at first the Fort Erie/Buffalo sill, before the Lyell/Johnson sill in Niagara Falls took over due to isostatic rebound. This sill, in time, was eventually eroded by the recession of Niagara Falls and the Fort Erie/Buffalo sill regained control. The environmental picture is complicated by catastrophic outbursts from glacial Lake Agassiz and Lake Barlow-Ojibway, changes in outlet routes, isostatic rebound and climatic changes over the Great Lakes basins. Today, the flow of water into Lake Erie from the streams and rivers surrounding it only accounts for about 13% of the flow out of it, therefore, the importance of flow from the Upper Great Lakes, specifically the flow from Lake Huron, has a great effect on Lake Erie levels. While the changing control sills, Lyell/Johnson and Buffalo/Fort Erie would affect Lake Erie levels, overall they are mostly input driven by the amount of waters received from the Upper Great Lakes. Since Lake Erie's water level changes are so closely tied to Lake Huron's water level changes we have decided to use names assigned to Lake Huron such as the two Mattawa highstands and three Stanley lowstands rather than inflict a whole new set of names on the public. While the duration of each high and lowstand in Erie and Huron may not always be the same, they always happen within the same time frame. The datum elevations used for Lake Huron (175.8 m) and Lake Erie (173.3 m) are historically recorded averages. The Lake Erie levels proposed in this paper reflect Lake Hurons effects on Lake Erie and the levels occuring at the eastern end of the Erie Basin throughout the last 12600 years. All dates in this paper are uncorrected ¹⁴ C dates unless the date was obtained from shells, then the date has been corrected for hard-water effects. Also, all heights are given as modern day elevations and are not adjusted for isostatic rebound.     

Oscillations of levels and cool phases of the Laurentian Great Lakes caused by inflows from glacial Lakes Agassiz and Barlow-Ojibway

Two distinct episodes of increased water flux imposed on the Great Lakes system by discharge from upstream proglacial lakes during the period from about 11.5 to 8 ka resulted in expanded outflows, raised lake levels and associated climate changes. The interpretation of these major hydrological and climatic effects, previously unrecognized, is mainly based on the evidence of former shorelines, radiocarbon-dated shallow-water sediment sequences, paleohydraulic estimates of discharge, and pollen diagrams of vegetation change within the basins of the present Lakes Superior, Michigan, Huron, Erie and Nipissing. The concept of inflow from glacial Lake Agassiz adjacent to the retreating Laurentide Ice Sheet about 11–10 and 9.5–8.5 ka is generally supported, with inflow possibly augmented during the second period by backflooding of discharge from glacial Lake Barlow-Ojibway.Although greater dating control is needed, six distinct phases can be recognized which characterize the hydrological history of the Upper Great Lakes from about 12 to 5 ka; 1) an early ice-dammed Kirkfield phase until 11.0 ka which drained directly to Ontario basin; 2) an ice-dammed Main Algonquin phase (11.0–10.5 ka) of relatively colder surface temperature with an associated climate reversal caused by greater water flux from glacial Lake Agassiz; 3) a short Post Algonquin phase (about 10.5–10.1 ka) encompassing ice retreat and drawdown of Lake Algonquin; 4) an Ottawa-Marquette low phase (about 10.1–9.6 ka) characterized by drainage via the then isostatically depressed Mattawa-Ottawa Valley and by reduction in Agassiz inflow by the Marquette glacial advance in Superior basin; 5) a Mattawa phase of high and variable levels (about 9.6–8.3 ka) which induced a second climatic cooling in the Upper Great Lakes area. Lakes of the Mattawa phase were supported by large inflows from both Lakes Agassiz and Barlow-Ojibway and were controlled by hydraulic resistance at a common outlet — the Rankin Constriction in Ottawa Valley — with an estimated base-flow discharge in the order of 200000 m³s^–1. 6) Lakes of the Nipissing phase (about 8.3–4.7 ka) existed below the base elevation of the previous Lake Mattawa, were nourished by local precipitation and runoff only, and drained by the classic North Bay outlet to Ottawa Valley.Geological Survey of Canada Contribution 42488.This is the twelfth of a series of papers to be published by this journal that was presented in the paleolimnology sessions organized by R. B. Davis and H. Löffler for the XIIth Congress of the International Union for Quaternary Research (INQUA), which took place in Ottawa, Canada in August 1987. Dr. Davis is serving as guest editor of this series.     

A late Lake Minong transgression in the Lake Superior basin as documented by sediments from Fenton Lake, Ontario

The evolution of the early Great Lakes was driven by changing ice sheet geometry, meltwater influx, variable climate, and isostatic rebound. Unfortunately none of these factors are fully understood. Sediment cores from Fenton Lake and other sites in the Lake Superior basin have been used to document constantly falling water levels in glacial Lake Minong between 9,000 and 10,600 cal (8.1–9.5 ka) BP. Over three meters of previously unrecovered sediment from Fenton Lake detail a more complex lake level history than formerly realized, and consists of an early regression, transgression, and final regression. The initial regression is documented by a transition from gray, clayey silt to black sapropelic silt. The transgression is recorded by an abrupt return to gray sand and silt, and dates between 9,000 and 9,500 cal (8.1–8.6 ka) BP. The transgression could be the result of increased discharge from Lake Agassiz overflow or the Laurentide Ice Sheet, and hydraulic damming at the Lake Minong outlet. Alternatively ice advance in northern Ontario may have blocked an unrecognized low level northern outlet to glacial Lake Ojibway, which switched Lake Minong overflow back to the Lake Huron basin and raised lake levels. Multiple sites in the Lake Huron and Michigan basins suggest increased meltwater discharges occurred around the time of the transgression in Lake Minong, suggesting a possible linkage. The final regression in Fenton Lake is documented by a return to black sapropelic silt, which coincides with varve cessation in the Superior basin when Lake Agassiz overflow and glacial meltwater was diverted to glacial Lake Ojibway in northern Ontario.     

An early Holocene oxygen isotope record from the Olson buried forest bed, Southern Lake Michigan

Analysis of a 3.5 m vibracore from the Olson buried forest bed in the southern Lake Michigan basin provides new paleolimnological data for the early Holocene. The core records a rise in lake level from the Chippewa low water phase toward the Nipissing high water phase. Deepening of the water level at the core site is suggested by a trend toward decreasing organic carbon content up core that is interpreted as a response to increasing distance between terrestrial debris sources and the core site.Published data from deep water cores from the southern Lake Michigan basin suggest there had been an inflow of isotopically light water from glacial Lake Agassiz into the southern basin between 10.5-11 ka (A1 event). The data also indicate a second flood of isotopically light water between 8-9 ka (A2 event).Three new ¹⁴C dates from the Olson site core suggest that most of the sediment was deposited between 8.45 ka and 8.2 ka, an interval roughly coeval with the second pulse of ¹⁸O-depleted water (A2) from Lake Agassiz into the southern basin. Oxygen isotope ratio analysis of shell aragonite from the gastropods Probythinella lacustris and Marstonia deceptashows increasingly negative values up core. This trend in¹⁸O values suggests that ¹⁸O - depleted water entered the southern basin about 8.4 ka. The Olson site core thus provides a chronology of events in the southern Lake Michigan basin associated with the draining of glacial Lake Agassiz.     

Great Lakes response to catastrophic inflows from Lake Agassiz: some simulations of a hydraulic geometry for chained lake systems

Simulations (216) were undertaken to evaluate the impact of typical Lake Agassiz outbursts on the upper Great Lakes under plausible variations in lake surface areas and sill widths. Flows over sills out of lakes are modelled using the equation for a broad-crested weir, with the model time increment set to one day. The model was evaluated for Lake Agassiz outlet sill widths of 1, 4, and 10 km and with outbursts ranging from 100 000 m³ s^–1 to 600 000 m³ s^–1. The surface area of Lake Agassiz was evaluated for 182 000 km² ±20%. The surface area of the upper Great Lakes were modelled as either Lake Algonquin (Superior, Huron and Michigan basins =200 000 km²) or Lake Minong (Superior basin 87 000 km²) with sill widths of 0.5, 1.5, and 3 km.Downstream peak discharge modelled at the outlet sill of the upper Great Lakes, was normally between 20 and 60% of the initial outburst, with a lagtime to peak usually between 80 and 280 days. Upper Great Lakes water level rises of between 2 and 20 m are calculated with rises to 36 m for some configurations. Rise magnitude is inversely related to the width of the outlet sills at both lake systems and to the surface area of the receiving lake.The modeling implies that measuring outflow from the upper Great Lakes, or water level rises, does not in itself determine peak or total outflow from Lake Agassiz unless the dimensions of the Lake Agassiz and upper Great Lakes outflow sills are also known.Lake level rises probably coincided on the upper Great Lakes with meltout from the winter freeze-up. Lake levels re-attain equilibrium values with respect to through flow within three years of an outburst. Substantial episodic lake level rises in the upper Great Lakes may have had severe impacts on the lake biota, for example via the affect on spawning grounds.     

Lake levels in the Erie Basin of the Laurentian Great Lakes

Water levels in the Lake Erie basin are inferred from glacial lake times to present. An era of early to middle Holocene lowstands is defined below outlets by a submerged paleo-beach, and truncated reflectors in glaciolacustrine sediment beneath a mud-covered wave-cut terrace. Also, the glacial clay surface above the paleo-shore level has elevated shear strength because of porewater drainage during subaerial exposure. Below the paleo-shore where exposure did not occur, clay strength remained normal. Sedimentation rates were reduced during the lowstands. The distortion of once-level shore zone indicators by differential glacial rebound was removed by computing original elevations of the indicators using an empirical model of rebound based on observations of upwarped former lake shorelines. Erie water-level history was inferred from a plot of the original elevations of lake-level constraints and outlets versus age. The lake history was validated by reference to ~83 water-level indicators, not used as constraints. During the deglaciation, lake-crossing moraines were likely eroded by fluvial drainage into low-level Lake Ypsilanti and a subsequent unnamed low lake to produce the Lorain Valley and Pennsylvania Channel. Once inflow from the upper Great Lakes basins was directed to Ottawa Valley about 10,400 (12,270 cal BP), Erie water levels descended in a dry, evaporative climate to a closed lowstand during which ostracode δ¹⁸O increased ~2‰ above present values. Lake level began to rise 6,000 to 7,000 (6,830 to 7,860 cal) BP in response to increased atmospheric moisture and later, to northern inflow as the Nipissing Transgression returned upper Great Lakes drainage to Lake Erie by about 5,200 (6,000 cal) BP. At that time, the lake overflowed the uplifted Lyell–Johnson Sill north (downstream) of the present Niagara Falls at higher-than-present levels. After recession of the Falls breached this sill about ~3,500 (~3,770 cal) BP, Lake Erie fell 3–4 m to its present Fort Erie–Buffalo Sill. The extended low-water phase with its isolated sub-basins could have restricted migration of aquatic fauna. The early to middle Holocene closed-basin response highlights the sensitivity of Lake Erie to climatic reductions in its water budget.     

Holocene paleoclimatic evidence and sedimentation rates from a core in southwestern Lake Michigan

Preliminary results of a multidisciplinary study of cores in southwestern Lake Michigan suggest that the materials in these cores can be interpreted in terms of both isostatically and climatically induced changes in lake level. Ostracodes and mollusks are well preserved in the Holocene sediments, and they provide paleolimnologic and paleoclimatic data, as well as biogenic carbonate for stable-isotope studies and radiocarbon dating. Pollen and diatom preservation in the cores is poor, which prevents comparison with regional vegetation records. New accelerator-mass spectrometer ¹⁴C ages, from both carbon and carbonate fractions, provide basin-wide correlations and appear to resolve the longstanding problem of anomalously old ages that result from detrital organic matter in Great Lakes sediments.Several cores contain a distinct unconformity associated with the abrupt fall in lake level that occurred about 10.3 ka when the isostatically depressed North Bay outlet was uncovered by the retreating Laurentide Ice Sheet. Below the unconformity, ostracode assemblages imply deep, cold water with very low total dissolved solids (TDS), and bivalves have ¹⁸O (PDB) values as light as — 10 per mil. Samples from just above the unconformity contain littoral to sublittoral ostracode species that imply warmer, higher-TDS (though still dilute) water than that inferred below the unconformity. Above this zone, another interval with ¹⁸O values more negative than — 10 occurs. The isotopic data suggest that two influxes of cold, isotopically light meltwater from Laurentide ice entered the lake, one shortly before 10.3 ka and the other about 9 ka. These influxes were separated by a period during which the lake was warmer, shallower, but still very low in dissolved solids. One or both of the meltwater influxes may be related to discharge from Lake Agassiz into the Great Lakes.Sedimentation rates appear to have been constant from about 10 ka to 5 ka. Bivalve shells formed between about 8 and 5 ka have ¹⁸O values that range from-2.3 to-3.3 per mil and appear to decrease toward the end of the interval. The ostracode assemblages and the stable isotopes suggest changes that are climatically controlled, including fluctuating water levels and increasing dissolved solids, although the water remained relatively dilute (TDS < 300 mg/l).A dramatic decrease in sedimentation rates occurred at about 5 ka, about the time of the peak of the Nippissing high lake stage. This decrease in sedimentation rate may be associated with a large increase in effective wave base as the lake approached its present size and fetch. A dramatic reduction in ostracode and mollusk abundances during the late Holocene is probably due to this decrease in sedimentation rates, which would result in increased carbonate dissolution. Ostracode productivity may also have declined due to a reduction in bottom-water oxygen caused by increased epilimnion algal productivity.Woods Hole Oceanographic Institute Contribution No. 7492     

Paleohydrologic implications of 18O enriched Lake Agassiz water

Lake Agassiz water oxygen isotopic compositions inferred from sediment core organics and pore waters provide some additional insight into the paleohydrology of the Great Lakes and their drainage into the North Atlantic during the late glacial and early Holocene. Isotopically enriched Lake Agassiz water supports the hypothesis that high Huron Basin lake (Mattawa) phases, during the early Holocene (9600–9300 and 9100–8100 years BP) resulted from an influx of Lake Agassiz water and suggests that low lake (Stanley) phases (9800–9600, 9300–9100, 8100–7400 years BP) were influenced more by regional influxes of isotopically depleted glacial melt water. Eastward drainage of enriched early Lake Agassiz water supports an active Port Huron outlet between 11000 and 10500 years BP and also helps to explain the absence of an ¹⁸O depleted interval in North Atlantic foram records. This may be the result of a balance between the opposing isotopic effects of depleted Lake Agassiz water and lower sea surface temperatures on carbonate precipitation between 11000 and 10000 years BP.     

Modeling the history of Lake of the Woods since 11,000 cal yr B.P. using GIS

Because of differential isostatic rebound, many lakes in Canada have continued to change their extent and depth since retreat of the Laurentide Ice Sheet. Using GIS techniques, the changing configuration and bathymetry of Lake of the Woods in Ontario, Manitoba, and Minnesota were reconstructed for 12 points in time, beginning at 11,000 cal yr B.P. (9.6 ¹⁴C ka B.P.), and were also projected 500 years into the future, based on the assumption that Lake of the Woods continued to have a positive hydrological budget throughout the Holocene. This modeling was done by first compiling a bathymetric database and merging that with subaerial data from the Shuttle Radar Topography Mission (SRTM). This DEM file was then adjusted by: (1) isobase data derived from Lake Agassiz beaches prior to 9000 cal yr B.P. (8.1 ¹⁴C ka B.P.) and (2) modeled isostatic rebound trend analysis after 9000 cal yr B.P. Just after the end of the Lake Agassiz phase of Lake of the Woods, only the northernmost part of the basin contained water. Differential rebound has resulted in increasing water depth. In the first 3000 years of independence from Lake Agassiz, the lake transgressed >50 km to the south, expanding its area from 858 to 2857 km², and more than doubling in volume. Continued differential rebound after 6000 cal yr B.P. (5.2 ¹⁴C ka B.P.) has further expanded the lake, although today it is deepening by only a few cm per century at the southern end. In addition, climate change in the Holocene probably played a role in lake level fluctuations. Based on our calculation of a modern hydrological budget for Lake of the Woods, reducing runoff and precipitation by 65% and increasing evaporation from the lake by 40% would end overflow and cause the level of the lake to fall below the outlets at Kenora. Because this climate change is comparable to that recorded during the mid-Holocene warming across the region, it is likely that the area covered by the lake at this time would have been less than that determined from differential isostatic rebound alone.     

The Stanley unconformity in Lake Huron basin: evidence for a climate-driven closed lowstand about 7900 <Superscript>14</Superscript>C BP,with similar implications for the Chippewa lowstand in Lake Michigan basin

J.L. Hough in 1962 recognized an erosional unconformity in the upper section of early postglacial lake sediments in northwestern Lake Huron. Low-level Lake Stanley was defined at 70 m below present water surface on the basis of this observation, and was inferred to follow the Main Algonquin highstand and Post-Algonquin lake phases about 10 ¹⁴C ka, a seminal contribution to the understanding of Great Lakes history. Lake Stanley was thought to have overflowed from the Huron basin through the Georgian Bay basin and the glacio-isostatically depressed North Bay outlet to Ottawa and St. Lawrence rivers. For this overflow to have occurred, Hough assumed that post-Algonquin glacial rebound was delayed until after the Lake Stanley phase. A re-examination of sediment stratigraphy in northwestern Lake Huron using seismic reflection and new core data corroborates the sedimentological evidence of Hough’s Stanley unconformity, but not its inferred chronology or the level of the associated lowstand. Erosion of previously deposited sediment, causing the gap in the sediment sequence down to 70 m present depth, is attributed to wave erosion in the shoreface of the Lake Stanley lowstand. Allowing for non-deposition of muddy sediment in the upper 20 m approximately of water depth as occurs in the present Great Lakes, the inferred water level of the Stanley lowstand is repositioned at 50 m below present in northwestern Lake Huron. The age of this lowstand is about 7.9 ± 0.3¹⁴C ka, determined from the inferred ¹⁴C age of the unconformity by radiocarbon-dated geomagnetic secular variation in six new cores. This relatively young age shows that the lowstand defined by Hough’s Stanley unconformity is the late Lake Stanley phase of the northern Huron basin, youngest of three lowstands following the Algonquin lake phases. Reconstruction of uplift histories for lake level and outlets shows that late Lake Stanley was about 25–30 m below the North Bay outlet, and about 10 m below the sill of the Huron basin. The late Stanley lowstand was hydrologically closed, consistent with independent evidence for dry regional climate at this time. A similar analysis of the Chippewa unconformity shows that the Lake Michigan basin also hosted a hydrologically closed lowstand, late Lake Chippewa. This phase of closed lowstands is new to the geological history of the Great Lakes. This is the ninth in a series of ten papers published in this special issue of Journal of Paleolimnology. These papers were presented at the 47th Annual Meeting of the International Association for Great Lakes Research (2004), held at the University of Waterloo, Waterloo, Ontario, Canada. P.F. Karrow and C.F.M Lewis were guest editors of this special issue.     

The history and architecture of lacustrine depositional systems in the northern Lake Michigan basin

The post-glacial history of the Great Lakes has involved changes in lake levels that are equivalent in vertical extent to the Pleistocene changes in global sea level and changes in sediment accumulation by at least two orders of magnitude. In the sediments of the northern Lake Michigan basin, these radical changes in base level and sediment supply are preserved in detailed records of changing depositional environment and the impact of these changes on depositional architecture. The seismic sequences of the sediment fill previously described in Lake Huron have been carried into northern Lake Michigan and used to map the history and architecture of basinal deposition. As the Laurentide Ice Sheet retreated northward in the early Holocene, it opened progressively deeper channels to the east that allowed the larger lakes to drain through the North Channel, Huron, and Georgian Bay basins. At the end of the Main Algonquin highstand, about 10,200 (radiocarbon) yrs ago, the eastern drainage passage deepened in a series of steps that defined four seismic sequences and lowered lake levels by over 100 m. Near the same time a new source of sediment and meltwaters poured across the Upper Peninsula of Michigan and into the northern Lake Michigan basin from the Superior basin ice lobe. A marked increase in deposition is seen first in the northern part of the study area, and slightly later in the Whitefish Fan area at the southern end of the study area. Accumulation rates in the area gradually decreased even as lake levels continued to fall. Drainage directly from the Superior basin ended before the beginning of the main Mattawa phase about 9,200 (radiocarbon) yrs ago.Although individual lowstand systems tracts are at the most a few hundred yrs in duration, their geometries and seismic character are comparable to marine systems tracts associated with sea level falls of similar magnitudes. In some of the thicker lowstand deposits a second order cyclicity in sedimentation can be detected in the high resolution seismic records.     

A new water-level history for Lake Ontario basin: evidence for a climate-driven early Holocene lowstand

Piston cores from deep-water bottom deposits in Lake Ontario contain shallow-water sediments such as, shell-rich sand and silt, marl, gyttja, and formerly exposed shore deposits including woody detritus, peat, sand and gravel, that are indicative of past periods of significantly lower water levels. These and other water-level indicators such as changes in rates of sedimentation, mollusc shells, pollen, and plant macrofossils were integrated to derive a new water-level history for Lake Ontario basin using an empirical model of isostatic adjustment for the Great Lakes basin to restore dated remnants of former lake levels to their original elevations. The earliest dated low-level feature is the Grimsby-Oakville bar which was constructed in the western end of the lake during a near stillstand at 11–10.4 (12.9–12.3 cal) ka BP when Early Lake Ontario was confluent with the Champlain Sea. Rising Lake Ontario basin outlet sills, a consequence of differential isostatic rebound, severed the connection with Champlain Sea and, in combination with the switch of inflowing Lake Algonquin drainage northward to Ottawa River valley via outlets near North Bay and an early Holocene dry climate with enhanced evaporation, forced Lake Ontario into a basin-wide lowstand between 10.4 and 7.5 (12.3 and 8.3 cal) ka BP. During this time, Lake Ontario operated as a closed basin with no outlets, and sites such as Hamilton Harbour, Bay of Quinte, Henderson Harbor, and a site near Amherst Island existed as small isolated basins above the main lake characterized by shallow-water, lagoonal or marsh deposits and fossils indicative of littoral habitats and newly exposed mudflats. Rising lake levels resulting from increased atmospheric water supply brought Lake Ontario above the outlet sills into an open, overflowing state ending the closed phase of the lake by ~7.5 (8.3 cal) ka BP. Lake levels continued to rise steadily above the Thousand Islands sill through mid-to-late Holocene time culminating at the level of modern Lake Ontario. The early and middle Holocene lake-level changes are supported by temperature and precipitation trends derived from pollen-climate transfer functions applied to Roblin Lake on the north side of Lake Ontario.     

Evidence of early Holocene closed-basin conditions in the Huron-Georgian basins from within the North Bay outlet of the upper Great Lakes

Sub-bottom profiling and coring were undertaken at eight sub-basins along the lower French River and at five small lakes near North Bay, Ontario, to collect stratigraphical and chronological evidence to investigate whether lakes occupying the Huron–Georgian basins during the early- to mid-Holocene became hydrologically closed. All of the coring sites are located within the route of the North Bay outlet that carried outflow from the upper Great Lakes during this period. Sand beds containing organic detritus are present within five cores from Muskrat, Crombie and Deep bays that otherwise are composed of glaciolacustrine rhythmites or fine-grained lacustrine deposits. These sand beds are interpreted to represent intervals when water levels within the sub-basins were lower than present, based on chronology, sediment texture, and macrofossil assemblages. It is inferred that the water surface in the Huron–Georgian basins fell below the level of the Dalles Rapids sill isolating the lower French River sub-basins from the large lake. A core from Depensier Lake, North Bay, contains an organic-rich sand interval within a thicker sand unit barren of organic materials. Macrofossils within this organic-rich interval are interpreted to be evidence of substantially diminished flow through the North Bay outlet channel. Radiocarbon dates of terrestrial macrofossils provide correlation of the sand beds between the French River cores as well as with the organic-rich sand in the Depensier Lake core. The possibility that the sand beds in the French River cores represent flood deposits rather than evidence of hydrologically closed conditions is considered, but rejected, based on the occurrence of multiple peaty layers and the record of shallow water conditions inferred from macrofossils within the upper sand bed of core MUS1, Muskrat Bay, in combination with the evidence of quiescent depositional conditions from similarly aged macrofossils in the core from Depensier Lake. Eight radiocarbon dates from the French River cores are incorporated into an elevation-age plot of paleo-indicators of water levels in the Huron–Georgian basins, using additional data from the literature. This plot and stratigraphic evidence from the Muskrat Bay cores indicates that separate closed-basin intervals occurred between 9.0 and 8.4, and 9.5 and 9.3 ka cal BP (~ 8.1 and 7.6, and ~ 8.5 and 8.3 ka BP). The occurrence of these two closed-basin intervals between 9.6 and 8.4 ka cal BP (~ 8.7 and 7.6 ka BP) implies that run-off derived exclusively from precipitation within the non-glaciated portions of the upper Great Lakes drainage basins was likely insufficient at this time to support an open-basin lake hydrology during the contemporary climate, which was colder and drier than present, without being supplemented from glacial Lake Agassiz overflow and/or Laurentide Ice Sheet meltwater.     

Early postglacial history of the southeastern Assiniboine Delta,glacial Lake Agassiz basin

New stratigraphic evidence from the Rossendale area, Manitoba, Canada, provides insight into the early postglacial evolution of the southeastern Assiniboine Delta. In this region, much of the upper 13+ m of sediment accumulation is characterized by multiple cycles of sandy rhythmites interbedded with massive to laminated silt. These sediments were deposited rapidly by traction or turbidity currents and record the construction of the Assiniboine fan-delta during the deep-water Lockhart Phase of glacial Lake Agassiz (>10.8 ¹⁴C ka BP). Shortly before ∼10 ¹⁴C ka BP, fluvial incision into deltaic deposits occurred locally at the Rossendale Gully site in response to the regression of glacial Lake Agassiz during the Moorhead Phase. Plant macrofossils deposited in the gully by 10 ¹⁴C ka BP provide the first information on early postglacial plant colonization of the distal Assiniboine delta. These data suggest initial establishment of Scorpidium scorpioides, Potamogeton spp., Scirpus spp., and other wetland plants, followed by colonization of uplands by a Picea-Populus assemblage. Importantly, because the gully is located in a protected depression behind the Campbell beach, evidence of water table rise from aquatic macrophytes suggests that glacial Lake Agassiz transgressed to the Campbell level during the early Emerson Phase (∼10 ¹⁴C ka BP). Furthermore, no evidence exists for a post-Lockhart rise in Lake Agassiz above the Upper Campbell beach. If Agassiz stood at the Campbell level during the early Emerson Phase, then drainage through the southern outlet may have been possible at this time. This scenario, if true, may suggest that the northwestern outlet was temporarily closed by a glacial advance shortly before 10 ¹⁴C ka BP. This is the first in a series of ten papers published in this special issue of Journal of Paleolimnology. These papers were presented at the 47th Annual Meeting of the International Association for Great Lakes Research (2004), held at the University of Waterloo, Waterloo, Ontario, Canada. P.F. Karrow and C.F.M. Lewis were guest editors of this special issue     

Constraints on paleolake levels,spillways and glacial lake history,north-central Ontario,Canada

The recognition of ice-marginal deltas constructed during the formation of the Nakina II moraine and a previously unrecognized spillway, in the vicinity of Longlac, northern Ontario, indicates that existing concepts of ancestral lake level history and drainage systems in the Lake Superior–Lake Nipigon region is inadequate. Based on isostatically corrected digital elevation maps, ice-marginal deltas of the Nakina II moraine probably formed at the level of glacial Lake Minong, most likely Minong III, and not glacial Lake Nakina as has been commonly suggested. In addition, the presence of a spillway near Longlac indicates that lake water drained southward through the Mullet Outlet–Pic River system immediately following ice-marginal retreat from the Nakina II moraine and not eastward as previously proposed. Architectural-element analysis of exposures within the spillway indicates hyperconcentrated outbursts of meltwater produced thick channel-fill elements during flood conditions with peak-velocities exceeding 3 m/s. Subsequent retreat of ice from the Pic River valley to the east, may have allowed waters of Lake Agassiz, Lake Barlow–Ojibway, or both, to drain into post-Minong lake levels in the Lake Superior basin. These findings place major constraints on previously proposed concepts of northeastern or eastern outlets of Lake Agassiz.     

The volume and paleobathymetry of glacial Lake Agassiz

During the last retreat of the Laurentide Ice Sheet in North America, many proglacial lakes formed as continental drainage was impounded against the southern and western ice margin. Lake Agassiz was the largest of these lakes. The bathymetry of Lake Agassiz at the Herman and Upper Campbell beach levels – formed at about 11.5–11.0 ka and 9.9–9.5 ka, respectively – was computer modelled in this study by first collecting data for the isostatically-deformed paleowater planes of the two lake levels (derived from isobase lines constructed from beach elevations), and then subtracting these from the modern topography of the former lake floor. Pixels with dimensions of 1/30 × 1/30 of a degree were used in the model. Using these data, the area and volume of the lake were also calculated: at the Herman level these were 152 500 km2 and 13 100 km3 respectively; at the Upper Campbell level these were 350 400 km2 and 38 700 km3. Contour maps showing the paleobathymetry of both periods in the lake's history were also constructed. Determining the paleobathymetry and volume of Lake Agassiz is an important step in understanding the impact that the lake had on its surrounding environment and on the rivers, lakes, and oceans into which it flowed.     

Lake Winnipeg: geological setting and sediment seismostratigraphy

Lake Winnipeg, the seventh largest lake in North America, is located at the boundary between the Interior Plains and the Canadian Shield in Manitoba, Canada. Seismic profiles were obtained in Lake Winnipeg on two geoscientific cruises in 1994 and 1996. These data indicate the morphology of the bedrock surface. In most cases, a clear distinction between low relief Paleozoic carbonate rock and high relief Precambrian rock can be made. In northern Lake Winnipeg, the eastern limit of Paleozoic rock is clearly demarcated 30 km west of the previous estimate of its position. In southern Lake Winnipeg, all or most of the Paleozoic sequence terminates at a prominent buried escarpment in the centre of the lake. This indicates that Paleozoic rock on the eastern shore, known from drilling and outcrops, is an outlier. Major moraines are apparent as abrupt, large ridges having a chaotic internal reflection pattern. These include the Pearson Reef Moraine, the George Island Moraine and the offshore extension of The Pas Moraine. Little evidence for extensive or thick till was observed. Instead, fine-grained sediments deposited in glacial Lake Agassiz rest directly on bedrock over most of the lake basin. Hence an episode of erosion to bedrock was associated with glaciation and/or deglaciation. The Agassiz Sequence sediments are well-stratified, drape underlying relief and in some areas are over 100 m thick. In places, stratification in these sediments is disrupted, perhaps by dewatering. Evidence of erosion of Agassiz Sequence sediments by recent currents was observed. The contact between the Agassiz Sequence and the overlying Winnipeg Sequence sediments is a marked angular unconformity. The Agassiz Unconformity indicates up to 10 m of erosion in places. The low-relief character of this unconformity precludes subaerial erosion and the lack of till, moraines, or extensive deformation precludes glacial erosion. Waves appear to be the most likely erosional agent, either in waning Lake Agassiz or early Lake Winnipeg time. Winnipeg Sequence sediments, in places very thin, mantle most of the lakefloor. These sediments were deposited in the present Lake Winnipeg and are faintly stratified to massive and reach about 10 m in thickness in deep water. On the surface of the Winnipeg Sequence, vigorous, episodic currents are thought to contribute to the construction of flow-transverse sand waves as much as 6 m high in a deep, narrow constriction in the lake.     

Variation in the composition of Lake Bonneville marl: a potential key to lake-level fluctuations and paleoclimate

Lake Bonneville marl provides a stratigraphic record of lake history preserved in its carbonate minerals and stable isotopes. We have analyzed the marl in shallow cores taken at three localities in the Bonneville basin. Chronology for the cores is provided by dated volcanic ashes, ostracode biostratigraphy, and a distinctive lithologic unit believed to have been deposited during and immediately after the Bonneville Flood.A core taken at Monument Point at the north shore of Great Salt Lake encompasses virtually the entire Bonneville lake cycle, including the 26.5 ka Thiokol basaltic ash at the base and deposits representing the overflowing stage at the Provo shoreline at the top of the core. Two cores from the Old River Bed area near the threshold between the Sevier basin and the Great Salt Lake basin (the main body of Lake Bonneville) represent deposition from the end of the Stansbury oscillation ( 20 ka) to post-Provo time ( 13 ka), and one core from near Sunstone Knoll in the Sevier basin provides a nearly complete record of the period when Lake Bonneville flooded the Sevier basin (20–13 ka).In all cores, percent calcium carbonate, the aragonite to calcite ratio, and percent sand were measured at approximately 2-cm intervals, and ¹⁸O and ¹³C were determined in one core from the Old River Bed area. The transgressive period from about 20 ka to 15 ka is represented in all cores, but the general trends and the details of the records are different, probably as a result of water chemistry and water balance differences between the main body and the Sevier basin because they were fed by different rivers and had different hypsometries. The Old River Bed marl sections are intermediate in position and composition between the Monument Point and Sunstone Knoll sections. Variations in marl composition at the Old River Bed, which are correlated with lake-level changes, were probably caused by changes in the relative proportions of water from the two basins, which were caused by shifts in water balance in the lake.This is the second paper in a series of papers published in this issue on Climatic and Tectonic Rhythms in Lake Deposits.     

Late Pleistocene and Holocene lake fluctuations in the Sevier Lake basin,Utah, USA

Sevier Lake is the modern lake in the topographically closed Sevier Lake basin, and is fed primarily by the Sevier River. During the last 12 000 years, the Beaver River also was a major tributary to the lake. Lake Bonneville occupied the Sevier Desert until late in its regressive phase when it dropped to the Old River Bed threshold, which is the low point on the drainage divide between the Sevier Lake basin and the Great Salt Lake basin. Lake Gunnison, a shallow freshwater lake at 1390 m in the Sevier Desert, overflowed continuously from about 12 000 to 10 000 yr B.P., into the saline lake in the Great Salt Lake basin, which continued to contract. This contrast in hydrologic histories between the two basins may have been caused by a northward shift of monsoon circulation into the Sevier Lake basin, but not as far north as the Great Salt Lake basin. Increased summer precipitation and cloudiness could have kept the Sevier Lake basin relatively wet.By shortly after 10 000 yr B.P. Lake Gunnison had stopped overflowing and the Sevier and Beaver Rivers had begun depositing fine-grained alluvium across the lake bed. Sevier Lake remained at an altitude below 1381 m during the early and middle Holocene. Between 3000 and 2000 yr B.P. the lake expanded slightly to an altitude of about 1382.3 m. A second expansion, probably in the last 500 years, culminated at about 1379.8 m. In the mid 1800s the lake had a surface altitude of 1379.5 m. Sevier Lake was essentially dry (1376 m) from 1880 until 1982. In 1984–1985 the lake expanded to a 20th-century high of 1378.9 m in response to abnormally high snow-melt runoff in the Sevier River. The late Holocene high stands of Sevier Lake were most likely related to increased precipitation derived from westerly air masses.This is the first of a series of papers to be published by this journal that was presented in the paleolimnology sessions organized by R. B. Davis and H. Löffler for the XIIth Congress of the International Union for Quaternary Research (INQUA), which took place in Ottawa, Canada in August 1987. Drs. Davis and Löffler are serving as guest editors of this series.     

Paleogeography of Lake Agassiz and regional post-glacial uplift history of the Flin Flon region, central Manitoba and Saskatchewan

Detailed mapping and elevation measurements of glacial lake shorelines in the Flin Flon region has permitted the reconstruction of 6 well-defined levels of Lake Agassiz formed around 8.3 ka to 7.9 ka ¹⁴C BP. The Stonewall, The Pas, Gimli, Grand Rapids, Drunken Point and Ponton paleo-water planes have been tilted upward to the northeast in the Holocene, with gradients decreasing, from the highest to the lowest level, from about 0.34 m km^-1 to 0.22 m km^-1 in the study area. The Setting level, lower than the Ponton but less well defined, is also documented here for the first time. This mapping conclusively refutes the view, entrenched in the literature from the 1890's to the 1960's, that there has been negligible differential uplift in the region following final drainage of Lake Agassiz. The finding has major consequences regarding correlation of glacial lakes across the mid-continent, the post-glacial history of large lakes in the region, and for interpretation of earth rheology and its implications for ice sheet reconstruction.