Methods

Ice and snow sampling was conducted weekly from 14 January to 30 March, 2010 on Lake Mendota when the ice was safe to walk on. A Kovacs Mark III core drill, manufactured by Ice Coring and Drilling Service (ICDS), Space Science and Engineering Center (SSEC) UW Madison, was used to collect ice cores. Snow depth was also measured at the locations where ice cores were sampled. All measurements were made in centimeters. Blue ice can be defined as the portion of the ice core that is strictly frozen lake water. White ice can be defined as “snow ice,” which occurs when water rushes through cracks in the ice and soaks the overlying snow, resulting in a mixture of ice and snow that subsequently freezes. Total ice is blue ice + snow ice. Finally, snow depth was calculated as the average of 10 snow depth samples at each sampling location.

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.