Tree-ring carbon isotope discrimination (Δ13C) and oxygen isotopes (δ18O) collected from white pine (Pinus strobus) trees adjacent to Lake Superior show potential to produce the first winter-specific paleoclimate reconstruction with inter-annual resolution for this region. Isotopic signatures from 1976 to 2015 were strongly linked to antecedent winter minimum temperatures (Tmin), Lake Superior peak ice cover, and regional to continental-scale atmospheric winter pressure variability including the North American Dipole. The immense thermal inertia of Lake Superior underlies the unique connection between winter conditions and tree-ring Δ13C and δ18O signals from the following growing season in trees located near the lake. By combining these signals, we demonstrate feasibility to reconstruct variability in Tmin, ice cover, and continental-scale atmospheric circulation patterns (r ≥ 0.65, P < 0.001).
Trees growing in cold environments do not directly record winter conditions in their tree-rings because they are dormant during this period. Isotopic signals imprinted upon tree-ring cellulose can serve as robust climate proxies that are conventionally known to record growing season conditions1. Indeed, the stable isotope composition of tree-ring cellulose has been shown to reflect canopy-integrated leaf responses to environmental drivers that are further modified by downstream ecophysiological processes2,3,4,5,6,7. More specifically, in temperate, near-coastal locations, the primary environmental signals recorded by tree-ring carbon and oxygen isotopes include warm season temperature, vapor pressure deficit, irradiance, and cloud cover8,9,10,11,12,13,14. Therefore, to our knowledge, no tree-ring signals of any type have been linked to winter temperatures or ice cover of lakes or oceans via well-understood meteorological and ecophysiological phenomena. To overcome this limitation, herein we first describe aspects of coastal “lake-effect” climate, and show how climate conditions adjacent to Lake Superior during spring and early summer are seasonally-lagged and driven by antecedent winter conditions. Thereafter we demonstrate how isotopic signals fixed in the central portion of each tree-ring appear to record winter temperature anomalies for the Great Lakes region of North America.
Lake Superior is the largest freshwater lake in the world by surface area and the third largest by volume. The resulting thermal mass produces “lake-effect” climate conditions that strongly regulate minimum and maximum air temperatures (Tmin and Tmax, respectively) and other meteorological variables near the lake-land-atmosphere interface15,16,17,18,19. Ice forms around the shallow perimeter of Lake Superior every winter, but the lake surface rarely freezes over completely due to the large heat storage capacity of the lake20,21,22. Winter areal ice cover generally peaks near the first week of March at about 40%, but thin ice forms and melts earlier compared to lags of approximately two and four weeks for medium and thicker ice, respectively (Supplementary Fig. 1a). Inter-annual variation in peak ice cover across thickness classes is controlled primarily by winter minimum temperatures (Supplementary Fig. 1b).
Data from NOAA buoys deployed since 1979 reveal that Lake Superior water surface temperatures and over-lake air temperatures during summer are inversely related to ice coverage during the previous winter23 (Supplementary Fig. 1c). Lake-cooled air temperatures occur near the shoreline and at some distance inland, as demonstrated by spatial patterns of significant correlations between antecedent winter Tmin and warm-season Tmax during June, July and August (Supplementary Fig. 1d). The inland-advection of cool lake breezes is driven in part by greater overland warming and convection that can modify air temperature regimes up to 40 km from one of the Great Lakes17,24,25,26. Moreover, it is notable that inter-annual variability in surface water temperatures is amplified compared to regional air temperatures, corresponding to decreases in ice cover of 0.49% yr−1 over the same period (Supplementary Fig. 2). Amplified water temperature variability has been thought to represent ice albedo-effects27 interacting with the timing of stratification of lake water temperatures23.
Climate responses of tree-ring carbon and oxygen isotopes are widely acknowledged to record growing season conditions1, and this includes climate conditions along the US West Coast10,28,29. The climate of the US West Coast is influenced by ocean temperatures, but compared to Lake Superior, ocean temperature variability is muted and influenced primarily by coastal upwelling and large-scale atmospheric and ocean interactions. Due to the strong variability in Lake Superior water temperatures, we hypothesized that by sampling stable isotope signals fixed in cellulose during the spring and early summer, trees will have recorded winter season conditions rather than growing season conditions as have been demonstrated in other coastal locations. Here we test this hypothesis using tree-ring carbon and oxygen isotopes collected from white pines (Pinus strobus L.) growing 2.3 km distance from the Lake Superior shoreline, at a location where strong summer lake-effect air temperature gradients have been documented using a high-density network of temperature sensors26.
We sampled five dominant or codominant white pine trees 2.3 km from the Lake Superior shoreline. Each tree was approximately 130 years in age and was located on the southern shoreline of Rush Lake, on the property of the Huron Mountain Club (Fig. 1). Three 12 mm-diameter cores were collected from each tree. Tree cores were mounted, surfaced and visually cross-dated. Thereafter rings were measured with a linear encoder (Velmex Inc., Bloomfield, NY) and MeasureJ2X software (Voortech Consulting) and statistically cross-dated against local and regional white pine chronologies (R. Fahey and S. Voelker, unpublished data) using COFECHA software30,31. After cross-dating, the middlewood28 (MW) and latewood components of each ring were excised with a Dremel tool and saved for analyses whereas the earliest cells formed in each ring were discarded (approximating 20% of the anatomically distinct earlywood). This practice minimized the potential for physiological carryover effects from previous growing season conditions28. Here, we only report data from the MW component of each ring – which should represent growth and isotope signals during late spring through early summer. After stable isotope analyses (described below), tree-ring isotope series from each tree were detrended with ARSTAN software32 using a 100-year spline. Detrending had the two-fold goals of (1) removing low frequency trends in the individual isotope time series that we expected were not due to climate, but which can result from changes in tree height, competition and rooting depth, and (2) retaining inter-annual to decadal-scale variation in the isotope time series that we expected to be related to inter-annual climate variability.
Isotopic composition is expressed using “delta” notation as δ13C or δ18O = (Rsample/Rstandard − 1) × 1000, where δ13C or δ18O is the molar ratio of heavy to light isotopes and Rstandard is Vienna Pee Dee Belemnite (VPDB) or Vienna Standard Mean Ocean Water (VSMOW), respectively. Wood was ground to a powder, heat-sealed within a polyester filter bag (ANKOM Technology, Macedon, NY), extracted to α-cellulose33,34, and weighed out on a microbalance. Carbon isotope ratios were obtained using standard high temperature combustion in a vario-Pyrocube elemental analyzer interfaced with IsoPrime/Elementar IsoPrime 100 gas phase isotope ratio mass spectrometer (IsoPrime Ltd., Manchester, UK) at the Roden laboratory at Southern Oregon University. Oxygen isotope ratios were determined by pyrolizing α-cellulose in an elemental analyzer (TC/EA, IsoPrime/Elementar vario-Pyrocube) and analyzing the resulting gas with an isotope ratio mass spectrometer (IsoPrime100) at the Dawson laboratory at UC-Berkeley. At these laboratories, the long-term precision is less than 0.1‰ for δ13C and 0.2‰ for δ18O. δ13C data were converted to Δ13C following Farquhar (1989) as: Δ13C = (δ13Cair − δ13Cplant)/(1 + δ13Cplant/1000) where δ13Cair was estimated from McCarroll and Loader (2004) through 2003 and merged seamlessly with a similarly smoothed record of annual δ13Cair from Mauna Loa, Hawaii (http://scrippsco2.ucsd.edu/data/atmospheric_co2/mlo) for 2004 to 2015.
Daily resolution temperature data, from buoys and meteorological stations near the lake were from the United States National Ocean and Atmospheric Association National Buoy Data Center (http://www.ndbc.noaa.gov/maps/WestGL.shtml) and the Canadian Government (http://climate.weather.gc.ca/). Daily resolution cloud cover data were from the NOAA Great Lakes Environmental Research Laboratory (https://www.glerl.noaa.gov/). The Big Bay and Marquette, MI meteorological stations are 15 and 55 km to the southwest of the tree-ring collection site, respectively. Both stations are located <0.1 km from Lake Superior. After summarizing temperature data by month at each site, missing and/or incomplete monthly values for either record were gap-filled using linear regressions developed between the two sites. Hereafter we define this site-averaged data set to be the “local” temperature signal. We also summarized meteorological data from all stations within 7 km of the Lake Superior, which we here define as the “regional near-lake” mean. Likewise, regional near-lake cloud cover data were summarized from airports within 7 km of Lake Superior.
For regional and local winter Tmin, we averaged monthly values across the previous November to the current March to provide a winter seasonal Tmin whereas summer seasonal data were averaged across the current June to August. One potential concern of using temperature data from lake-effect locations is that summer Tmax is significantly (P < 0.01) influenced by winter Tmin, whereas upwind locations that are outside of the influence of Lake Superior do not show these patterns (Supplementary Fig. 1d). To remove covariation among these signals, we obtained residuals, from a linear regression relating summer Tmax to winter Tmin using data summarized from Marquette and Big Bay (data not shown), to produce a warm-season temperature record independent of the effect of winter Tmin, which are hereafter referred to as Tmax*.
Weekly ice cover estimates were obtained from the Canadian Ice Service (http://iceweb1.cis.ec.gc.ca/IceGraph/page1.xhtml). Thin, medium, and thick ice cover were averaged into a single value per year across weeks 7–12, 9–14 and 11–16, respectively, with the week of January 1st defined as week zero. This procedure combines ice cover signals while accounting for lags in the formation of medium and thick ice (Supplementary Fig. 1a).
Winter cold outbreaks across the Midwest, southern Canada and Eastern United States have been associated with low atmospheric pressure or geopotential height (HGT) anomalies, generally centered between Lake Superior and the eastern edge of Hudson Bay, Canada35,36,37. The fluctuation of such winter “stationary waves” in the upper atmosphere modulates surface air temperature across the Great Lakes Region. Indeed, the “Polar Vortex” event of the 2013–14 winter easily set a new record of areal ice cover on Lake Superior since records have been kept. This included a peak coverage of 100%, with some ice remaining into early June visible from MODIS images (Fig. 1a). The winter climate over North America and particularly the Great Lakes is affected by mean-state stationary atmospheric waves, characterized by a pair of circulation features in the upper troposphere: A high-pressure ridge over western North America and a low-pressure trough over the Hudson Bay area38, referred to as the dipole36,37. To characterize this continental-scale circulation pattern, we obtained HGT anomaly data at 250 hPa from the two centers of this dipole: over the Northeastern Pacific Ocean and the eastern edge of Hudson Bay, Canada, respectively. The difference between these anomalies has been called the “dipole index”35,36,39 and the HGT anomaly values forming this dipole index were derived from data with a 2.5° latitude × longitude resolution from the National Centers for Environmental Prediction/National Center for Atmospheric Research Reanalysis40. Linear regression and multiple linear regression models were constructed and compared using the “lm” or “glm” packages in the R statistical computing environment version 3.4.341.
Winter Tmin, the strongest driver of ice formation on Lake Superior, differs greatly from year to year, as demonstrated by MODIS images from early March that contrast ice conditions during the relatively warm winter of 2011–2012 (Fig. 1a) versus the cold and so-called “Polar Vortex” winter of 2013–2014 (Fig. 1b). In spring and summer, the thermal inertia of Lake Superior causes the lake surface to be much cooler than adjacent land surfaces, resulting in subsidence of air passing over the lake and advection of cool lake breezes inland while stabilizing the surrounding atmosphere. These effects can be visualized collectively through a snapshot of cloud cover surrounding the perimeter of Lake Superior, but essentially no cloud cover over the lake (Fig. 1c).
White pine tree-ring δ18O and Δ13C chronologies were characterized by series inter-correlation of 0.76 and 0.54, respectively, while expressed population signals42 were 0.94, and 0.85, respectively. These characteristics are a hallmark of chronologies that have a robust common signal driven by climatic variability. To determine the strongest climate drivers of Δ13C and δ18O, multiple regression analyses were conducted that included local (winter Tmin and summer Tmax* from the two closest meteorological stations), regional (summer cloud cover and lake-wide peak ice cover %,), and large-scale variables (winter HGT and winter dipole index) for the current year (Table 1).
In contrast to most other studies of tree-ring temperature signals, only winter season variables were found to be significant and thereafter retained in the final multiple regression models for white pine MW Δ13C and δ18O (Table 1). Ice cover and local winter Tmin values had the strongest influence on the MW Δ13C and δ18O, respectively, whereas HGT and dipole index also described significant variation in the same models (Table 1). The modeling process was then inverted, with the goal of identifying targets for climate reconstruction utilizing MW Δ13C and δ18O.
All winter-season variables were strongly predicted by individual isotope chronologies and their dual-isotope combinations (Fig. 2). For the variables included in climate reconstructions (Fig. 2), variance inflation factors were low, <2, indicating that multicollinearity did not contribute to over-fitting the model. Across all variables there were no significant relationships with summer Tmax* or summer near-lake cloud cover (Table 2). It is notable for interpreting these results that summer Tmax* is the summer temperature record after having removed the impact of winter Tmin due to the lagging response of lake-effect climate (See Methods). This means that stable isotope variation did not reflect variability in regional summer climate conditions that were independent of from lake-effect climate variability. The strength of the relationships noted in Table 2 and that between ice cover and the dipole index increased only weakly through time (Supplementary Fig. 3). These results clearly indicated that, for this species, spring and early summer conditions near the lakeshore were dominated by the influence of winter season climate signals.
Compared to models using current and antecedent year isotopic signals (Table 2), when only current year signals were considered, climate variables that operate at local to regional scales, such as winter Tmin and ice cover, were predicted somewhat better compared to larger-scale influences on winter atmospheric circulation, such as winter HGT and the winter dipole index (Supplementary Table 1). Current year Δ13C predicted variability in all of the winter climate signals reasonably well, with little additional explanatory power provided by current year δ18O for either of the winter HGT or dipole index (Supplementary Table 1). However, predictions of Tmin and peak ice cover were substantially improved by including Δ13C of the subsequent year, and predictions of HGT and the winter dipole index improved substantially by including δ18O of the subsequent year (cf. Table 2 and Supplementary Table 1). Atmospheric pressure patterns and attendant “waviness” of the jet stream produces winter weather conditions at a given location but these large scale phenomena do not directly affect isotopic signatures. Likewise, ice cover on Lake Superior is largely a product of winter Tmin (Supplementary Fig. 1b). Therefore, ultimately, winter temperatures drove tree-ring isotopic variation over time, with spatial signatures of this relationship being strongest over south-central Canada, the upper peninsula of Michigan and Northern Wisconsin (Supplementary Fig. 4).
A regression map of the peak ice cover regressed with the 250 hPa geopotential height eddies (HGT, i.e. with the zonal mean removed for the depiction of stationary waves in the upper atmosphere) elucidates large-scale circulation patterns producing cold anomalies and associated inter-annual variability in ice cover (Fig. 3a). As expected, an east-west oriented wave-train pattern is revealed over North America that spans the North Pacific to North Atlantic. For comparative purposes, the climatological stationary waves are overlaid as contours, outlining the “normal” position of the stationary eddies during winter (Fig. 3a). Over North America, the regression pattern of the circulation associated with the inter-annual variation in ice cover corresponds closely to the winter-mean stationary wave patterns, signaling that amplified stationary waves produce cold anomalies across the northeast U.S. and southern Canada36,39,43 and associated increases in ice cover for Lake Superior. Similarly, in Fig. 3b, the regression map of HGT with the dipole index predicted from tree-ring Δ13C and δ18O (i.e., from Fig. 2d) also shows a dipole pattern, despite a slight difference in the ridge intensity over Alaska. In contrast to the similar dipole patterns in Fig. 3a,b, the regression map of HGT with the well-known Pacific-North America pattern (PNA)44 shows a distinctly different and spatially shifted pattern compared to the dipole position. As a result, tree-ring Δ13C and δ18O chronologies and peak ice cover were not significantly correlated with the PNA (data not shown).
Our results provide evidence that tree-ring Δ13C and δ18O of white pines growing adjacent to Lake Superior show potential to provide novel proxy records of winter temperature anomalies with inter-annual resolution. Although we sampled only five trees at one location, the resulting tree-ring Δ13C and δ18O chronologies were able to predict previous winter temperature and ice conditions locally and regionally, as well as depicting the dominant large-scale atmospheric circulation pattern associated with these winter temperature and ice conditions (Figs 2, 3; Table 2). There are limitations to using a data set from only one site, since climatic controls on tree-ring isotopic variation can differ depending on microsite conditions, and specifically with distance to water bodies45. It is possible that the sampling location was more responsive to lake effect climate than others at similar distances from Lake Superior or that adding additional sites may further strengthen winter climate signals affecting this region. Additional research would logically address whether comparable winter climate-isotope signals can be obtained from other locations around Lake Superior, to what extent distance from Lake Superior or variation in local geomorphology can modify these signals. Toward this end, tree-ring isotope carbon chronologies from other species from this region have demonstrated that upland sites with relatively thin soils overlying bedrock tend to more closely reflect the degree of summer drought stress (S. Voelker, unpublished data). Therefore, to obtain coherent winter climate signals from tree-ring isotopes it will likely be critical to sample trees from mesic but well-drained sites, which often occur next to lakes. Sampling trees from such sites has the added benefit of the potential to extend chronologies back in time by locating sub-fossil wood preserved in lakes that grew under similar climate and hydrologic conditions46. Finally, at any given site, tree ring-width and stable isotope records may not display long-term stability with climate47, so it will be important to extend isotope chronologies further back in time to ascertain this possibility over the full instrumental climate record. Despite the importance of sampling considerations for future studies aimed at exploiting the unique lake-effect climate signals demonstrated here, it is important to first address which physiological and climatic processes best explain how tree-ring carbon and oxygen isotopic signals record antecedent winter temperatures in our current data set.
According to carbon isotope theory, Δ13C is primarily controlled by the ratio of photosynthetic assimilation rate to stomatal conductance (A/gs)48, which results in two physiological processes which can be modified by climate variability. Tree-ring Δ13C is often strongly influenced by gs, as indicated by tree-ring Δ13C relationships with soil moisture deficits and/or atmospheric drought stress14,48,49,50,51. However, gs is unlikely to constrain Δ13C near Lake Superior during spring and early summer because the air is cool and humid, the soil should still be wet from snowmelt, and precipitation >5 mm occurs on about 50% of days. Indeed, monthly Palmer Drought Severity Index showed no significant relationships with MW Δ13C (Supplementary Table 2). Theory relating Δ13C to A/gs45 suggests the primary remaining variable that could potentially explain inter-annual variation in MW Δ13C, is how late spring and early summer air temperatures affect A, as influenced by antecedent winter conditions. Near Marquette, average air temperatures during June and July include a Tmin of 11.8 °C and Tmax of 22.4 °C. A responds positively across this range of air temperatures, peaking near 20 °C for most temperate species52,53. Therefore, close to the lake, Δ13C signals from early in the growing season are fixed in the MW component of tree-rings in association with frigid previous winters. This is caused by near-lake temperatures during this period being inversely proportional to ice cover such that A is constrained by cooler Tmax adjacent to the lake26 and more strongly so after winters with colder Tmin.
Tree-ring δ18O is primarily influenced by source water isotopic signatures and the degree to which leaf water is either enriched by evaporation or influenced by atmospheric water vapor2,3,4,14,51. Source water isotopic signals near Lake Superior surely contributed to some of the inter-annual variation in the tree-ring δ18O patterns observed due to how temperatures influence precipitation δ18O54,55,56. Monthly temperatures can explain about 67% of the variation in monthly-resolution precipitation δ18O data from 1995 to 2015 located 180 km North across Lake Superior from our collection site at Sibley, Ontario (F. Longstaffe, unpublished data). White pine MW δ18O also had a strong relationship (R2 = 0.47, P < 0.001) with the same Sibley precipitation δ18O data centered on the late winter to early spring period over the previous two years (Supplementary Fig. 5). Finally, Sibley precipitation δ18O averaged across the previous November to current March was significantly correlated with HGT and the Dipole index (r > 0.49, P < 0.027, data not shown), thereby linking winter atmospheric circulation, precipitation δ18O and tree-ring δ18O. The apparent integration of late winter precipitation δ18O signals across two previous years rather than just one likely reflects the time it takes for winter snow to melt and move down the slope of the north face of Huron Mountain, to the base of the mountain where seeps and springs commonly occur along the shoreline of Rush Lake near where the white pine trees were sampled for this study. Hence, this link to conditions across two previous winters for tree-ring δ18O may not be as strong at sites with less varied local topography. On the other hand, at other sites that lack a strong influence of lateral soil water flow, we would expect the current year precipitation δ18O signals to have a relatively stronger influence on tree-ring δ18O.
Lake-effect precipitation is known to affect large spatial gradients in water isotopes, across distances of 50 kilometers or more in lower Michigan57. However, there is only a very weak relationship between winter Tmin and winter precipitation amounts from data compiled across all near-lake stations including the years 1900–2015 (R2 = 0.02, P = 0.1574) or across only the Marquette and Big Bay locations including the years 1966–2015 (R2 = 0.01, P = 0.4114). Tree-ring δ18O also had very weak correlations with winter precipitation records from all near-lake stations or from Marquette and Big Bay (r = −0.11 and 0.06, respectively, data not shown) and summer precipitation (r = −0.16 and −0.20, respectively, data not shown). Monthly Palmer Drought Severity Index also showed no significant relationships with MW δ18O (Supplementary Table 2). The PNA pattern contributes to continent-wide patterns in precipitation δ18O from the Northwest to Southeast across central North America, but interpolations of isotopic variation suggest there is no difference in precipitation δ18O between PNA phases for Lake Superior and the upper Great Lakes Region58.
Another influence on inter-annual variation in tree-ring δ18O surely derives from leaf water δ18O dynamics during the spring and early summer. Indeed, the dominant control over tree-ring δ18O was air temperatures during the early growing season, which were in turn strongly influenced by antecedent winter Tmin and other unmeasured winter conditions producing variation in winter ice cover (Table 2). This pattern can largely be attributed to a reduction in leaf temperature and leaf water evaporative enrichment by cold and humid air adjacent to the lake, particularly under conditions where lake breezes advect cool air inland. Overall, tree-ring δ18O formed during the spring and early summer is influenced by winter temperatures via source water (late winter precipitation δ18O) and indirectly via leaf water enrichment (lake-effect lags in lake temperature due to previous winter conditions). In contrast, the amount of lake-effect precipitation or phases of the PNA appear to have had little influence on tree-ring δ18O.
Agreement between winter-mean stationary wave eddies and the wave-train patterns produced by regression maps of HGT with Lake Superior ice cover and the dipole index predicted from tree-ring Δ13C and δ18O (Fig. 3a,b) contrasts sharply with long-term mean atmospheric circulation patterns associated with the PNA pattern (Fig. 3c). Physically, this translates to the dipole pattern reflecting amplification or attenuation of the mean stationary atmospheric wave pattern59, whereas the PNA describes atmospheric circulation patterns that are comparatively shifted in space.
The potential for using tree-ring isotopes to better understand past variation in winter conditions is strongly evident from the data and analyses presented here. Both carbon and oxygen isotopic signals fixed in the MW of white pines reflected winter conditions through the lagging lake-effect climate near the shore of Lake Superior interacting with well-known ecophsyiological drivers of how trees and leaves record environmental influences. Hence, there appears to be strong potential for reconstructing winter temperatures or atmospheric circulation patterns influencing this region using an intra-annual tree-ring sampling and a dual isotope approach. These same climate signals are also highly correlated with winter ice cover for Lake Superior, which raises the possibility of better understanding past ice dynamics that are a major determinant of annual evaporation and associated lake levels for Lake Superior60 and can also affect the limnological metabolism of this important ecosystem61.
McCarroll, D. & Loader, N. J. Stable isotopes in tree rings. Quat. Sci. Rev. 23, 771–801 (2004).
Roden, J. S., Lin, G. & Ehleringer, J. R. A mechanistic model for interpretation of hydrogen and oxygen isotope ratios in tree-ring cellulose. Geochim. Cosmochim. Acta 64, 21–35 (2000).
Barbour, M. M., Roden, J. S., Farquhar, G. D. & Ehleringer, J. R. Expressing leaf water and cellulose oxygen isotope ratios as enrichment above source water reveals evidence of a Péclet effect. Oecologia 138, 426–435 (2004).
Kahmen, A. et al. Cellulose δ18O is an index of leaf-to-air vapor pressure difference (VPD) in tropical plants. Proc. Natl. Acad. Sci. 108, 1981–1986 (2011).
Cernusak, L. A. et al. Environmental and physiological determinants of carbon isotope discrimination in terrestrial plants. New Phytol. 200, 950–965 (2013).
Cernusak, L. A. et al. Stable isotopes in leaf water of terrestrial plants. Plant, Cell Environ. 39, 1087–1102 (2016).
Wieloch, T. et al. Intramolecular 13C analysis of tree rings provides multiple plant ecophysiology signals covering decades. Sci. Rep. 8, 5048 (2018).
Szymczak, S., Joachimski, M. M., Brauning, A., Hetzer, T. & Kuhlemann, J. A 560 yr summer temperature reconstruction for the Western Mediterranean basin based on stable carbon isotopes from Pinus nigra ssp. Laricio (Corsica/France). Clim. Past 8, 1737–1749 (2012).
Young, G. H. F. et al. Changes in atmospheric circulation and Arctic Oscillation preserved within millennial length reconstruction of summer cloud cover from northern Fennoscandia. Clim. Dyn. 39, 495–507 (2012).
Johnstone, J. A., Roden, J. S. & Dawson, T. E. Oxygen and carbon stable isotopes in coast redwood tree rings respond to spring and summer climate signals. J. Geophys. Res. Biogeosci. 118, 1–13 (2013).
Konter, O. et al. Climatic sensitivity and parameter coherency in annually resolved δ13C and δ18O from Pinus uncinata tree-ring data in the Spanish Pyrenees. Chem. Geol. 377, 12–19 (2014).
Urrutia-Jalabert, R. et al. Increased water use efficiency but contrasting tree growth patterns in Fitzroya cupressoides forests of southern Chile during recent decades. J. Geophys. Res. Biogeosci. 120, 2505–2524 (2015).
Sakashita, W. et al. Hydroclimate reconstruction in central Japan over the past four centuries from tree-ring cellulose δ18O. Quat. Int. 455, 1–7 (2017).
Voelker, S. L., Roden, J. S. & Dawson. T. E. Millennial-scale tree-ring isotope chronologies from coast redwoods provide insights on controls over California hydroclimate variability. Oecologia. 187, 897-909 (2018).
Changnon, S. A. & Jones, D. M. A. Review of the influences of the Great Lakes on weather. Water Resour. Res. 8, 360–371 (1972).
Phillips, D. W. Environmental climatology of Lake Superior. J. Great Lakes Res. 4, 288–309 (1978).
Scott, R. W. & F. A. Huff. Lake effects on climatic conditions in the Great Lakes Basin. Illinois State Water Survey, Champaign, IL, 73 p. (1997).
Notaro, M. et al. Influence of the Laurentian Great Lakes on regional climate. J. Clim. 26, 789–804 (2013).
Kopec, R. J. Continentality around the Great Lakes. Bull. Am. Meteorol. Soc. 46, 54–57 (1965).
Assel, R., Cronk, K. & Norton, D. Recent trends in Laurentian Great Lakes ice cover. Clim. Change 57, 185–204 (2003).
Wang, J. et al. Temporal and spatial variability of Great Lakes ice cover, 1973–2010. J. Clim. 25, 1318–1329 (2012).
Clites, A. H. et al. Cold water and high ice cover on Great Lakes in spring 2014. EOS 95, 305–306 (2014).
Austin, J. A. & Colman, S. M. Lake Superior summer water temperatures are increasing more rapidly than regional air temperatures: A positive ice-albedo feedback. Geophys. Res. Lett. 34, L06604 (2007).
Cole, H. S. & W. A. Lyons. The impact of the Great Lakes on the air quality of urban shoreline areas; Some practical applications with regard to air pollution control policy and environmental decision making. In: Proceedings, Fifteenth Conference on Great Lakes Research. Ann Arbor: International Association for Great Lakes Research, pp. 436–463 (1972).
Ryznar, E. & Touma, J. S. Characteristics of true lake breezes along the eastern shore of Lake Michigan. Atmos. Environ. 15, 1201–1205 (1981).
Hinkel, K. M. & Nelson, F. E. Spatial and temporal aspects of the lake effect on the southern shore of Lake Superior. Theor. and Appl. Climatol. 109, 415–428 (2012).
Curry, J. A., Schramm, J. L. & Ebert, E. E. Sea ice-albedo climate feedback mechanism. J. Clim. 8, 240–247 (1995).
Roden, J. S., Johnstone, J. A. & Dawson, T. E. Intra-annual variation in the stable oxygen and carbon isotope ratios of cellulose in tree rings of coast redwood (Sequoia sempervirens). Holocene 19, 189–197 (2009).
Saffell, B. J. et al. Tree-ring stable isotopes record the impact of a foliar fungal pathogen on CO2 assimilation and growth in Douglas-fir. Plant, Cell Environ. 37, 1536–1547 (2014).
Holmes, R. Computer assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 44, 69–75 (1983).
Grissino-Mayer, H. D. Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA. Tree-Ring Res. 57, 205–221 (2001).
Cook, E. R. & Krusic, P. J. ARSTAN version 44h3: A tree-ring standardization program based on detrending and autoregressive time series modeling, with interactive graphics. Tree-Ring Laboratory, Lamont-Doherty Earth Observatory of Columbia University, Palisades, New York, USA (2014).
Leavitt, S. W. & Danzer, S. R. Methods for batch processing small wood samples to holocellulose for stable-carbon isotope analysis. Anal. Chem. 65, 87–89 (1993).
Sternberg, L. S. L. Oxygen and hydrogen isotope measurements in plant cellulose analysis., In: Modern methods of plant analysis, Volume 10, Plant fibers (eds Linskens H. F. & Jackson J. F.) pp 89–99, Heidelberg, Springer-Verlag (1989).
Wang, S.-Y., Hipps, L., Gillies, R. R. & Yoon, J.-H. Probable causes of the abnormal ridge accompanying the 2013-14 California drought: ENSO precursor and anthropogenic warming. Geophys. Res. Lett. 41, 3220–3226 (2014).
Wang, S.-Y., Huang, W.-R. & Yoon, J.-H. The North American winter ‘dipole’ and extremes activity: A CMIP5 assessment. Atmos. Sci. Lett. 16, 338–345 (2015).
Wang, S.-Y., Yoon, J.-H., Becker, E. & Gillies, R. California from drought to deluge. Nat. Clim. Change 7, 465–468 (2017).
Held, I. M., Ting, M. & Wang, H. Northern winter stationary waves: Theory and modeling. J. Clim. 15, 2125–2144 (2002).
Singh, D. et al. Recent amplification of the North American winter temperature dipole. J. Geophys. Res. Atmos. 121, 9911–9928 (2016).
Kalnay, E. et al. The NCEP/NCAR 40-Year Reanalysis Project. Bull. Am. Meteorol. Soc. 77, 437–471 (1996).
R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, https://www.R-project.org/ (2017).
Wigley, T. M. L., Briffa, K. R. & Jones, P. D. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. J. Appl. Meteorol. Climatol. 23, 201–213 (1984).
Schulte, J. A. & Lee, S. Strengthening North Pacific influences on United States temperature variability. Sci. Rep. 7, 124 (2017).
Horel, J. D. & Wallace, J. M. Planetary-scale atmospheric phenomena associated with the Southern Oscillation. Mon. Weather Rev. 109, 813–829 (1981).
Esper, J. et al. Site-specific climatic signals in stable isotope records from Swedish pine forests. Trees 32, 855–869 (2018).
Naulier, M. et al. Carbon and oxygen isotopes of lakeshore black spruce trees in northeastern Canada as proxies for climatic reconstruction. Chem. Geol. 374, 37–43 (2014).
Naulier M. et al Temporal instability of isotopes-climate statistical relationships – A study of black spruce trees in northeastern Canada. Dendrochronologia 34, 33–42.
Farquhar, G. D., Ehleringer, J. R. & Hubick, K. T. Carbon isotope discrimination and photosynthesis. Ann. Rev. Plant Biol. 40, 503–537 (1989).
Roden, J. S. & Ehleringer, J. R. Summer precipitation influences the stable oxygen and carbon isotopic composition of tree-ring cellulose in Pinus ponderosa. Tree Physiol. 27, 491–501 (2007).
Voelker, S. L., Meinzer, F. C., Lachenbruch, B. J., Brooks, R. & Guyette, R. P. Drivers of radial growth and carbon isotope discrimination of bur oak (Quercus macrocarpa Michx.) across continental gradients in precipitation, vapor pressure deficit and irradiance. Plant, Cell Environ. 37, 766–779 (2014).
Ratcliff, C. J., Voelker, S. L. & Nolin, A. W. Tree-ring carbon isotope records from the western Oregon Cascade Mountains primarily record summer maximum temperatures. Tree-Ring Res. 74, 185-195 (2018).
Berry, J. & Björkman, O. Photoysnthetic response and adaptation to temperature in higher plants. Ann. Rev. Plant Physiol. 31, 491–543 (1980).
Helliker, B. & Richter, S. L. Subtropical to boreal convergence in tree-leaf temperatures. Nature 454, 511–514 (2008).
Voelker, S. L. et al. Isolating relative humidity: dual isotopes δ18O and δD as deuterium deviations from the global meteoric water line. Ecol. Appl. 24, 960–975 (2014).
Dansgaard, W. Stable isotopes in precipitation. Tellus 16, 436–468 (1964).
Kohn, M. J. & Welker, J. M. On the temperature correlation of δ18O in modern precipitation. Earth Planet. Sci. Lett. 231, 87–96 (2005).
Bowen, G. J. Spatial analysis of intra-annual variation of precipitation isotope ratios and its climatological corollaries. J Geophys. Res. 113, D05113 (2008).
Bowen, G. J., Kennedy, C. D., Henne, P. D. & Zhang, T. Footprint of recycled water subsidies downwind of Lake Michigan. Ecosphere 3, 1–16 (2012).
Liu, Z., Yoshimura, K., Bowen, G. J. & Welker, J. M. Pacific-North American teleconnection controls on precipitation isotopes (δ18O) across the contiguous United States and adjacent Regions: a GCM-based analysis. J. Clim. 27, 1046–1061 (2014).
Gronewold, A. D. et al. Coasts, water levels, and climate change: A Great Lakes perspective. Clim. Change 120, 697–711 (2013).
Brothers, S. & Sibley, P. Light may have triggered a period of net heterotrophy in Lake Superior. Limnol. Oceanogr. 63, 1785–1798 (2018).
We gratefully acknowledge the support of the Huron Mountain Wildlife Foundation for providing travel funds, lodging, and access to the Huron Mountain Club property. We also thank two anonymous reviewers who provided comments that have greatly improved this manuscript. This is Western’s Laboratory for Stable Isotope Science Contribution #371. S.L.V. was supported in part by the National Science Foundation Paleo Perspectives on Climate Change awards AGS-1003050 and AGS-1003601 and the Utah Agriculture Experiment Station Project 1304.
The authors declare no competing interests.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Voelker, S.L., Wang, S.Y.S., Dawson, T.E. et al. Tree-ring isotopes adjacent to Lake Superior reveal cold winter anomalies for the Great Lakes region of North America. Sci Rep 9, 4412 (2019). https://doi.org/10.1038/s41598-019-40907-w
Linking canopy‐scale mesophyll conductance and phloem sugar δ 13 C using empirical and modelling approaches
New Phytologist (2020)