Stable isotope mass balance of the Laurentian Great Lakes

    We investigate the physical limnology of the Laurentian Great Lakes of North America using a new dataset of ¹⁸O/¹⁶O and ²H/¹H ratios from over 500 water samples collected at multiple depths from 75 stations during spring and summer of 2007. δ¹⁸O and δ²H values of each lake plot in distinct clusters along a trend parallel to, but offset from, the Global Meteoric Water Line, reflecting the combined effects of evaporative enrichment and the addition of precipitation and runoff along the chain lake system. We apply our new dataset to a stable-isotope-based evaporation model that explicitly incorporates downwind lake effects, including humidity build-up and changes to the isotope composition of atmospheric vapor. Our evaporation estimates are consistent with previous mass transfer results for Michigan, Huron, Ontario and Erie, but not for Superior, which has a much longer residence time. Calculated evaporation from Superior is ~300 mm per year, less than previous estimates of ~500 mm per year, likely arising from integration of the ‘isotopic memory' of lower evaporation rates under cooler climatic conditions with greater ice-cover than the present. Uncertainties in the estimates from the stable-isotope-based model are comparable to mass transfer results, offering an independent technique for evaluating evaporation fluxes.