Methods

Data were collected from two lakes, Sparkling Lake (46.008, -89.701) and Peter Lake (46.253, -89.504), both located in the northern highlands Lake District of Wisconsin and the Upper Peninsula of Michigan over a 10 day period on each lake in July and August of 2007. Refer to Van de Bogert et al. 2011 for limnological characteristics of the study lakes.Measurements of dissolved oxygen and temperature were made every 10 minutes using multiple sondes dispersed horizontally throughout the mixed-layer in the two lakes (n=35 sondes for Sparkling Lake and n=27 sondes for Peter Lake). Dissolved oxygen measurements were corrected for drift.Conductivity, temperature compensated specific conductivity, pH, and oxidation reduction potential were also measured by a subset of sensors in each lake. Of the 35 sondes in Sparkling Lake, 31 were from YSI Incorporated: 15 of model 600XLM, 14 of model 6920, and 2 of model 6600). The remaining sondes placed in Sparkling Lake were 4 D-Opto sensors, Zebra-Tech, LTD. In Peter Lake, 14 YSI model 6920 and 13 YSI model 600XLM sondes were used.Sampling locations were stratified randomly so that a variety of water depths were represented, however, a higher density of sensors were placed in the littoral rather than pelagic zone. See Van de Bogert et al. 2012 for the thermal (stratification) profile of Sparkling Lake and Peter Lake during the period of observation, and for details on how locations were classified as littoral or pelagic. In Sparkling Lake, 11 sensors were placed within the shallowest zone, 12 in the off-shore littoral, and 6 in each of the remaining two zones, for a total of 23 littoral and 12 pelagic sensors. Similarly, 15 sensors were placed in the two littoral zones, and 12 sensors in the pelagic zone.Sensors were randomly assigned locations within each of the zones using rasterized bathymetric maps of the lakes and a random number generator in Matlab. Within each lake, one pelagic sensor was placed at the deep hole which is used for routine-long term sampling.Note that in Sparkling Lake this corresponds to the location of the long-term monitoring buoy. After locations were determined, sensors were randomly assigned to each location with the exception of the four D-Opto sensor is Sparkling Lake, which are a part of larger monitoring buoys used in the NTL-LTER program. One of these was located near the deep hole of the lake while the other three were assigned to random locations along the north shore, south shore and pelagic regions of the lake. Documentation: Van de Bogert, M.C., Bade, D.L., Carpenter, S.R., Cole, J.J., Pace, M.L., Hanson, P.C., Langman, O.C., 2012. Spatial heterogeneity strongly affects estimates of ecosystem metabolism in two north temperate lakes. Limnology and Oceanography 57, 1689-1700.

Methods

Raw data (in English units) were assembled by Douglas Clark - Wisconsin State Climatologist. Data were converted to metric units and adjusted for temporal biases by Dale M. Robertson. For adjustments applied to various parameters see Robertson, 1989 Ph.D. Thesis UW-Madison. Adjusted data represent the BEST estimated daily data and may be raw data. Data collected at Washburn observatory, 8-1-1883 to 9-30-1904. Data collected at North Hall, 10-1-1904 to 12-31-1947 Data collected at Truax Field (Admin BLDG), 1-1-1948 to 12-31-1959. Data collected at Truax Field (Center of Field), 1-1-1960 to Present. Much of the data after 1990 were obtained in digital form from Ed Hopkins, UW-Meteorology. Data starting in 2002-05 were obtained from Sullivan at <a href="http://www.weather.gov/climate/index.php?wfo=mkx%20">http://www.weather.gov/climate/index.php?wfo=mkx</a> ,then go to CF6 and download monthly data to Madison_sullivan_conversion. Since Robertson (1989) adjusted all historical data to that collected from 1884-1989; no adjustments were applied to the recent data except for (1) wind and (2) estimated vapor pressure:(1) Wind after January 1997, and only wind from the southwest after November 2007, was extended by Dale M. Robertson and Yvonne Hsieh.In 1996, a discontinuity in the wind record was caused by change in observational techniques and sensor locations (Mckee et al. 2000). To address the non-climatic changes in wind speed, data from MSN were carefully compared with those collected from the tower of the Atmospheric and Oceanic Science Building at the University of Wisconsin-Madison, see http://ginsea.aos.wisc.edu/labs/mendota/index.htm. Hourly data from both sites (UMSN,hourly and UAOS,hourly) during 2003&ndash;2010 were used to form a 4&times;12 (four components of wind direction &times; 12 months) matrix (K4,12) of wind correction factors, yielding UAOS,daily= Ki,j&times;UMSN,daily. The comparison results indicated that the MSN weather station reported a higher magnitude in winds out of the east by 5% and lower magnitude in winds out of the west and south by 30% and 10%. The adjusted wind data (=Ki,j&times;UMSN,daily) were therefore employed and used in the model simulation. After adjustments, there was a decrease in wind velocities starting shortly before 1996. Overall the adjusted wind data had a decline in wind velocities of 16% from 1988&ndash;93 to 1994&ndash;2009) compared to a 7% decline at a nearby weather station with no known observational changes (St. Charles, Illinois; 150 km southeast of Lake Mendota). (2) Estimated vapor pressure was updated (after April 2002) by using the equation from DYRESM for estimation of vapor pressure (a function of both air temperature and dew point temperature); where a=7.5, b=237.3, and c=.7858.