Mid Holocene lake level and shoreline behavior during the Nipissing phase of the upper Great Lakes at Alpena, Michigan, USA

The Nipissing phase was the last pre-modern high-water stage of the upper Great Lakes. Represented as either a one- or two-peak highstand, the Nipissing occurred following a long-term lake-level rise. This transgression was primarily an erosional event with only the final stage of the transgression preserved as barriers, spits, and strandplains of beach ridges. South of Alpena, Michigan, mid to late Holocene coastal deposits occur as a strandplain between Devils Lake and Lake Huron. The landward part of this strandplain is a higher elevation platform that formed during the final stage of lake-level rise to the Nipissing peak. The pre-Nipissing shoreline transgressed over Devils Lake lagoonal deposits from 6.4 to 6.1 ka. The first beach ridge formed ~ 6 ka, and then the shoreline advanced toward Lake Huron, producing beach ridges about every 70 years. This depositional regression produced a slightly thickening wedge of sediment during a lake-level rise that formed 20 beach ridges. The rise ended at 4.5 ka at the Nipissing peak. This peak was short-lived, as lake level fell > 4 m during the following 500 years. During this lake-level rise and subsequent fall, the shoreline underwent several forms of shoreline behavior, including erosional transgression, aggradation, depositional transgression, depositional regression, and forced regression. Other upper Great Lakes Nipissing platforms indicate that the lake-level change observed at Alpena of a rapid pre-Nipissing lake-level rise followed by a slower rise to the Nipissing peak, and a post-Nipissing rapid lake-level fall is representative of mid Holocene lake level in the upper Great Lakes.     

Factors influencing microplastic abundances in nearshore,tributary and beach sediments along the Ontario shoreline of Lake Erie

Sediment samples were collected from nearshore, tributary and beach environments within and surrounding the northern part of Lake Erie, Ontario to determine the concentrations and distribution of microplastics. Following density separation and microscopic analysis of 29 samples, a total of 1178 microplastic particles were identified. Thirteen nearshore samples contained 0–391 microplastic particles per kg dry weight sediment (kg^?1), whereas 4 tributary samples contained 10–462?kg^?1 and 12 beach samples contained 50–146?kg^?1. The highest concentrations of nearshore microplastics were from near the mouths of the Detroit River in the western basin and the Grand River in the eastern basin, reflecting an urban influence. The highest microplastic concentrations in beach samples were determined from Rondeau Beach in the central basin where geomorphology affects plastics concentration. The Welland Canal sample in the eastern basin contained the greatest concentration of microplastics of the tributary samples, which is consistent with high population density and shipping traffic. The overall abundance of microplastic in northern Lake Erie nearshore, tributary and beach samples is 6 times lower than in sediment sampled from northern Lake Ontario. The nearshore and beach sample results potentially reflect the transport patterns of floating plastics modeled for Lake Erie, which predict that the majority of plastic particles entering the lake are transported to southern shoreline regions rather than northern areas.     

Beach sediment magnetism and sources: Lake Erie,Ontario, Canada

Measurements were made along the northwestern shore of Lake Erie, Canada to determine whether grain magnetic properties can be used to identify and distinguish sources of beach sediment. Although surface magnetic susceptibilities were highly variable, ranging from 56 to 9867 × 10⁻⁵ SI (Bartington MS2D), there was generally a gradual increase from the low beach (near the waterline) towards the high beach; there were also narrow, shore-parallel bands with high susceptibility at various points on the beach surface. Magnetic mineralogy on the beaches was dominated by low-Ti magnetite (570° < Tc < 580 °C), and the effective grain-size varied from pseudosingle domain in the low beach to multidomain on the high beach. Sandy bluff sediments in the eastern part of the study area had magnetic properties (e.g. S-ratios, hysteresis loops, thermomagnetic curves) that were similar to those on the beaches, whereas the magnetic properties of the extensive till bluffs and river basin sediments were quite different. The data suggested that, whereas the beaches in the western part of the study area are supplied with sediment from bluffs several tens of kilometres to the east, the source of the high magnetic concentrations on the eroding beaches of eastern Point Pelee remains to be determined.     

Freshwater mussel community response to warm water discharge in western Lake Erie

Native unionid mussels are endangered in the Laurentian Great Lakes due to habitat degradation and biofouling by invasive dreissenids. However, a robust community was discovered living within the thermal discharge of a power plant at Oregon, Ohio, on the south shore of Lake Erie. Our study compared this community to nearby communities outside the thermal plume, and examined habitat characteristics that may affect unionids. Unionids were sampled from the exposed lake bed at three sites during a seiche in 2011: (1) within the thermal plume, (2) at Bayshore Park (2.0 km east of the plant), and (3) at the University of Toledo's Lake Erie Center (4.0 km east). In 2010, sediment samples were collected along a 2 km transect extending east from the plant discharge roughly parallel to the south shore of Lake Erie. Results indicated that the community within the thermal plume had higher densities, higher diversity (H′), more small individuals but overall larger sizes than communities outside the plume. Both the rate and intensity of fouling by dreissenids were lower within the plume. Dry mass of coarse surface sediment and sediment organic matter content were negatively related to distance from the plant (R² = 0.497, and 0.479, respectively). An unexpected discovery was that the bulk of the coarse sediment was comprised of shell material from Asian clams and dreissenid mussels, suggesting a contribution of these exotic species to sediment accumulation. In total, our results suggest that several habitat characteristics close to the power plant are favorable to unionids.     

Ground-penetrating radar stratigraphy and depositional model for evolving Late Holocene aeolian dunes on the Lake Huron coast,Ontario

One of the largest (36 km²) Late Holocene systems anywhere in Ontario lies in Pinery Provincial Park on the southeast coast of Lake Huron and postdates a 12‐km‐long beach barrier system left by Lake Nipissing some 5000 years ago. Dunes lie parallel to the lake shoreline, oblique to the dominant west–northwest wind. Linear dunes in the north of the field regardless of age are stable and of moderate height (< 10 m) having grown in situ from foredunes on emergent beach ridges. Southwards however, the same dunes evolve into wider, higher (up to 20 m) parabolic forms with large blow outs. This change in form from north to south appears to have been a persistent evolutionary trend over the past 5000 years. Some 5 km of high resolution radar profiles identifies the changing internal stratigraphy of dunes as they evolve in shape. Linear dunes are composed of a ‘vertical aggradation sequence’ (VAS) that records initial formation of a foredune on top of shore-parallel beach ridges, and upwards growth of the dune by sand trapped by vegetation. Southwards over some 12 km, dunes grow in height and assume a parabolic form expressed on radar profiles by a lowermost VAS overlain by cross-bedded sand comprising a landward accretion sequence (LAS) recording landward migration. A depositional model relates the change in dune shape and stratigraphy to southward increasing sand supply within a large persistent littoral cell in Lake Huron.     

Fecal indicator organism density in beach sands: Impact of sediment grain size,uniformity, and hydrologic factors on surface water loading

Fecal loading to beach sands and subsequent transport to surface water may result in the degradation of surface water quality. To better understand the relationship between Escherichia coli in sands, beach hydrology, and recreational water quality, we collected surface water, groundwater, and sand samples from three Lake Michigan beaches with varying substrates (fine sand to pebbles, July 2005–June 2006). Sediment cores collected within transects perpendicular to and at fixed distances from the shoreline (0 m, 10 m, 20 m) and surface water samples collected at a depth of 1 m were analyzed for E. coli. Grain size analysis was performed on duplicate core samples to assess the relationship between E. coli density and mean grain size and uniformity. Groundwater samples, collected from shallow drive-point piezometers within the test area, were also analyzed for E. coli. E. coli density in beach sands differed significantly with distance from shore with the highest density occurring at the berm crest (0 m). Mean grain size and uniformity accounted for variation in E. coli density with fine sand of uniform distribution having the highest content. E. coli density in surface water was correlated to E. coli density in beach sand samples at the berm crest. E. coli in groundwater was < 10 to 579 MPN/100 ml (2005); none was detected in 2006. Management interventions, including altered beach grooming practices and slope assessments, may be effective in reducing E. coli content at beaches comprised of fine sands of uniform grain size, hence reducing water quality advisories.     

The effect of Lake Erie water level variations on sediment resuspension

Variability in Lake Erie water levels results in variations of the fluid forces applied to the lake bed by free-surface gravity wind-waves. An increase in the bed stress may re-suspend sediment deposited years earlier. This study identifies areas of possible non-cohesive sediment mobilization in response to the forcing conditions and water levels present in Lake Erie. Observations from NOAA buoy 45005 were used to identify wave events generated by a variety of atmospheric forcing conditions. For each event, numerical predictions of significant wave height, wave period, and water level from the Great Lakes Forecasting System (GLFS) were used to characterize the wave event variability over the lake. The Shields parameter was estimated at each 2 km × 2 km grid cell with the local wave forcing as predicted by GLFS assuming an estimate of the wave-induced friction factor. In the Cleveland harbor region of the central basin, the Shields parameter was also estimated by assuming uniform wave conditions as observed by NOAA buoy 45005. The “contour of incipient motion” for both variable and uniform wave events was defined as the offshore contour where the Shields parameter exceeds the critical limit for motion. Comparisons with a radiometrically corrected image from Landsat-7 showed that the spatially varying wave events from GLFS were in qualitative agreement with the satellite observations. A sensitivity analysis of wave height, wave period, and grain size showed the contour of incipient motion to be the most sensitive to wave period. Calculations performed for record high and low water levels showed that the incipient motion of non-cohesive sediments in the relatively flat central basin to be the most sensitive to the historic hydrologic variability present in Lake Erie.     

A modeling study of ice–water processes for Lake Erie applying coupled ice-circulation models

A hydrodynamic model that includes ice processes and is optimized for parallel processing was configured for Lake Erie in order to study the ice–water coupling processes in the lake. A hindcast from April 2003 to December 2004 with hourly atmospheric forcing was conducted. The model reproduced the seasonal variation of ice cover, but the development of ice extent in January and its decay in March somewhat preceded the observations. Modeled lake circulation in ice-free seasons is consistent with previous studies for Lake Erie. Thermal structure of the lake was reasonably comparable to both satellite-derived observations and in-situ measurements, with mean differences ranging from − 2 °C to 4 °C, depending on the season. The impacts of ice–water stress coupling and basal melting of ice were examined based on numerical experiments. The results show that: 1) ice–water stress coupling significantly dampens the subjacent lake circulation in winter due to packed ice cover that slows down the surface water, and 2) basal melting of ice contributes to widespread ice cover in the lake. The demonstrated model validity could lead to further studies of ice–water processes in the lake, including interannual variation and impacts on ecosystems.     

Urban pocket-beach morphodynamics along the wave-dominated southwest coast of Lake Michigan: An analysis of shoreline and sand volumetric changes

Many of the world’s beaches are embayed, but while a large body of work addresses the geomorphology of pocket beaches in oceanic settings, little is known about urban analogs, especially within the Great Lakes of North America. Groins and jetties shelter these systems from direct interaction with littoral processes, which elsewhere can influence how changes in lake level, winter-ice cover, and wave climate impact beach evolution. We address the direct controls of these forcing parameters on beach morphodynamics over a 33-yr period at North Point Beach, which is confined to an engineered ‘container’ along Lake Michigan’s wave-dominated SW margin. Analysis of near-annual beach change suggests lake-level change is the dominant geomorphic driver over inter-annual to decadal timeframes, with winter ice playing a secondary role. Pocket-beach shoreline positions were found to be unreliable indicators of sand volumetric changes. Lake-level rise facilitated shoreline retreat and overwash-induced beach accretion while high lake levels created the accommodation for additional sands to enter the embayment. This is important for coastal managers to consider when developing mitigation strategies for ongoing lake-level fluctuations and anticipated regional climate impacts. This foundational assessment has implications for embayed beaches of the greater Chicago coastal margin (n > 20), where many other site-specific variables (e.g., groin orientations and shoreline aspect) may factor into nearshore and onshore beach morphodynamics. Continued research into urban pocket beaches of the Great Lakes stands to offer useful information on the impacts of littoral fragmentation on coastal sediment routing during different lake-level phases and degrees of littoral interconnectivity.     

Sediment Oxygen Demand in the Central Basin of Lake Erie

Three separate procedures were used to estimate the sediment oxygen demand (SOD) in the central basin of Lake Erie and were compared with other estimates determined previously and with historical data. First, whole core incubations involved sealing sediment cores at 12°C to ensure no interaction between the overlying water and the atmosphere and monitoring continuously to define the linear disappearance of oxygen. Second, sediment plugs were placed inside flow-through reactors and the influent and effluent concentrations were monitored to obtain steady-state reaction rates. Third, an extensive data set for the central basin of Lake Erie was compiled for input into the diagenetic BRNS model, and the SOD was calculated assuming all primary redox reactions, but no secondary reactions. All three procedures produced estimates of SOD that were in reasonable agreement with each other. Whole core incubations yield an average SOD of 7.40 × 10⁻¹² moles/cm²/sec, the flow-through experiments had an average SOD of 4.04 × 10⁻¹² moles/cm²/sec, and the BRNS model predicts an SOD of 7.87 × 10⁻¹² moles/cm²/sec over the top 10 cm of sediment and appears to be calibrated reasonably well to the conditions of the central basin of Lake Erie. These values compare reasonably well with the 8.29 × 10⁻¹² moles/cm²/sec obtained from diffusion modeling of oxygen profiles (Matisoff and Neeson 2005). In contrast, values reported from the 1960s to 1980s ranged from 10.5–32.1 × 10⁻¹² moles/cm²/sec suggesting that the SOD of the central basin has decreased over the last 35 years, presumably, in response to the decrease in phosphorus loadings to Lake Erie. However, since hypoxia in the hypolimnion persists these results suggest that improvement in hypolimnetic oxygen concentrations may lag decreases in loadings or that the hypolimnion in the central basin of Lake Erie is simply too thin to avoid summer hypoxia during most years.