Lake ice records used to detect historical and future climatic changes

Historical ice records, such as freeze and breakup dates and the total duration of ice cover, can be used as a quantitative indicator of climatic change if long homogeneous records exist and if the records can be calibrated in terms of climatic changes. Lake Mendota, Wisconsin, has the longest uninterrupted ice records available for any lake in North America dating back to 1855. These records extend back prior to any reliable air temperature data in the midwestern region of the U.S. and demonstrate significant warming of approximately 1.5 °C in fall and early winter temperatures and 2.5 °C in winter and spring temperatures during the past 135 years. These changes are not completely monotonie, but rather appear as two shorter periods of climatic change in the longer record. The first change was between 1875 and 1890, when fall, winter, and spring air temperatures increased by approximately 1.5 °C. The second change, earlier ice breakup dates since 1979, was caused by a significant increase in winter and early spring air temperatures of approximately 1.3 °C. This change may be indicative of shifts in regional climatic patterns associated with global warming, possibly associated with the Greenhouse Effect.With the relationships between air temperature and freeze and break up dates, we can project how the ice cover of Lake Mendota should respond to future climatic changes. If warming occurs, the ice cover for Lake Mendota should decrease approximately 11 days per 1 °C increase. With a warming of 4 to 5 °C, years with no ice cover should occur in approximately 1 out of 15 to 30 years.     

A new method to estimate ice age temperatures

On glacial time scales, the waxing and waning of the Eurasian and North American ice sheets depend largely on variations in atmospheric temperature. As global sea level is primarily determined by the volume of these ice sheets, there is a direct (yet complex) relation between global sea level and the northern hemispheric (NH) temperature. This relation is essentially represented by a model of the NH ice sheets. We use a thermomechanical ice-sheet–ice shelf–bedrock model in conjunction with an inverse method to deduce a time series of NH temperature (from 120 kyr BP until present) that is consistent with the observed global sea level record. The advantage of this method is that it provides the annual mean surface air temperature averaged over the NH continents north of 40°N. The results reveal that ice age temperatures were 4–10°C lower than today, which agrees with other temperature reconstructions. However, reconstructed temperatures are comparitively low during the early stages of the glacial, a feature that is consistent with the rapid growth of the ice sheets. The sensitivity of the results for uncertainties in precipitation rate, in observed sea level and in some other model parameters is examined to quantify the error in reconstructed temperature. During the glacial period (120–15 kyr BP), surface air temperatures in the NH (north of 40°N) were 7.2±1.5°C lower than todays (interglacial) temperatures.     

Sensitivity of freezing drizzle formation in stably stratified clouds to ice processes

Summary This study examines the impact of ice formation and growth processes on freezing drizzle formation in stably stratified clouds. In particular we investigate the reason why freezing drizzle is rarely observed in clouds with top temperatures less than –15°C. We also investigate the sensitivity of freezing drizzle formation to the Hallett Mossop secondary ice process (Hallet and Mossop, 1974). The evaluation is performed by simulating cloud formation over a two-dimensional idealized mountain using a detailed microphysical scheme. The height and width of the two-dimensional mountain were designed to produce an updraft pattern with extent and magnitude similar to documented freezing drizzle cases. The simulations show that: (i) drizzle formation is very sensitive to the ice crystal concentration, with a significant reduction in the area over which drizzle forms and the maximum drizzle water content as the cloud top temperature decreases below –10°C, and (ii) secondary ice crystal formation has a significant effect on drizzle formation at cloud top temperatures below –10°C.The above two factors are likely the main cause for the lack of freezing drizzle at cloud top temperatures less than –15°C. We also found that neglecting the depletion of ice forming nuclei resulted in considerable overestimation of the ice crystal concentration and suppression of drizzle, even for the –10°C case.     

On the Response of the Greenland Ice Sheet to Greenhouse Climate Change

Numerical computations are performed with the three-dimensional polythermal ice-sheet model SICOPOLIS in order to investigate the possible impact of a greenhouse-gas-induced climate change on the Greenland ice sheet. The assumed increase of the mean annual air temperature above the ice covers a range from T = 1°C to 12°C, and several parameterizations for the snowfall and the surface melting are considered. The simulated shrinking of the ice sheet is a smooth function of the temperature rise, indications for the existence of critical thresholds of the climate input are not found. Within 1000 model years, the ice-volume decrease is limited to 10% of the present volume for T 3°C, whereas the most extreme scenario, T = 12°C, leads to an almost entire disintegration, which corresponds to a sea-level equivalent of 7 m. The different snowfall and melting parameterizations yield an uncertainty range of up to 20% of the present ice volume after 1000 model years.     

Comparative Study of Nucleation by Different Alcoholic Solutions of Benzoin, and Benzoin Dust

Nucleation properties of different alcoholic solutions of Benzoin and Benzoin dust have been studied by varying the seeding temperature in a laboratory Cold Chamber. In the present study, three different alcoholic solutions are being used and these are due to the standard Ethyl alcohol, one higher alcohol i.e. Isopropyle alcohol and one lower alcohol i.e. Methyl alcohol. In addition, Benzil dust and droplets of Acetone and Chloroform have been used to examine if these agents have any ice nucleation ability. It has been observed that all alcoholic solutions of Benzoin have sufficient ice nucleation ability at the negative side of 0^°C, with the maximum concentration occurring at −14.4^°C. In case of Benzoin dust, the nucleation ability is better on the positive side of 0^°C than the negative side. On the other hand, Benzil dust has insignificant nucleation ability and the droplets of two other agents have no nucleation ability at all.Considering the chemical as well as crystalline structure, one can conclude that the presence of OH group in the seeding molecule plays the central role in ice nucleation/hydrate crystal formation. In fact, a number of supplementary experiments like crystal growth in an incubator, XRD analysis and NMR study are performed to identify the nature of crystals formed by Benzoin dust above 0^°C. One can conclude that these are hydrate crystals of Benzoin with water caging around Benzoin crystals.     

Response of the Antarctic ice sheet to future greenhouse warming

Possible future changes in land ice volume are mentioned frequently as an important aspect of the greenhouse problem. This paper deals with the response of the Antarctic ice sheet and presents a tentative projection of changes in global sea level for the next few hundred years, due to changes in its surface mass balance. We imposed a temperature scenario, in which surface air temperature rises to 4.2° C in the year 2100 AD and is kept constant afterwards. As GCM studies seem to indicate a higher temperature increase in polar latitudes, the response to a more extreme scenario (warming doubled) has also been investigated. The mass balance model, driven by these temperature perturbations, consists of two parts: the accumulation rate is derived from present observed values and is consequently perturbed in proportion to the saturated vapour pressure at the temperature above the inversion layer. The ablation model is based on the degree-day method. It accounts for the daily temperature cycle, uses a different degree-day factor for snow and ice melting and treats refreezing of melt water in a simple way. According to this mass balance model, the amount of accumulation over the entire ice sheet is presently 24.06 × 10¹¹ m³ of ice, and no runoff takes place. A 1°C uniform warming is then calculated to increase the overall mass balance by an amount of 1.43 × 10¹¹ m³ of ice, corresponding to a lowering of global sea level with 0.36 mm/yr. A temperature increase of 5.3°C is needed for the increase in ablation to become more important than the increase in accumulation and the temperature would have to rise by as much as 11.4°C to produce a zero surface mass balance. Imposing the Bellagio-scenario and accumulating changes in mass balance forward in time (static response) would then lower global sea level by 9 cm by 2100 AD. In a subsequent run with a high-resolution 3-D thermomechanic model of the ice sheet, it turns out that the dynamic response of the ice sheet (as compared to the direct effect of the changes in surface mass balance) becomes significant after 100 years or so. Ice-discharge across the grounding-line increases, and eventually leads to grounding-line retreat. This is particularly evident in the extreme case scenario and is important along the Antarctic Peninsula and the overdeepened outlet glaciers along the East Antarctic coast. Grounding-line retreat in the Ross and Ronne-Filchner ice shelves, on the other hand, is small or absent.     

A thermomechanical model of ice flow in West Antarctica

 This study uses a three-dimensional thermo-mechanical model to investigate the internal flow dynamics of the West Antarctic Ice Sheet (WAIS). The model allows ice thickness, flow and temperature to interact freely. Its domain is prescribed as that of the present-day grounded WAIS. Realistic present-day climatic and topographical boundary conditions are employed. The analysis of a series of experiments pays particular attention to the location and dynamics of concentrations of ice flow (ice streams). Underlying topographic troughs are crucial in determining the strength and location of these concentrations of flow. The flow pattern generated by subglacial troughs is made more distinct by the inclusion of ice flow/temperature coupling. The inclusion of sliding leads to the generation of limit cycles in the ice flow. They are concentrated around the present-day ice streams B and C of the Siple Coast and have a period of 5 to 10 ky. There appears to be competition between several preferred ice flow pathways in this area. The two end members of the flow regime are a strong ice stream C with a weakened ice stream A/B complex, and strong ice streams A and B with a dormant ice stream C. Ice streams appear to require ice discharges above a certain threshold in order to maintain frictional heat generation and fast flow. Individual ice streams can therefore interact through changes in catchment-area size: a reduction in catchment area reduces the volume of ice entering a stream and can cause stagnation as the amount of frictional heating falls. Received: 22 July 1997/Accepted: 27 July 1998     

Studies of ice nucleation by Arctic aerosol on AGASP-II

Aerosol particles were collected on filters for studies of their ability to nucleate ice during the second Arctic Gas and Aerosol Sampling Program (AGASP-II) in April, 1986. The ice nuclei (IN) samples were collected from an aircraft at altitudes ranging from the surface to the vicinity of the tropopause in Arctic locations over Alaska, northern Canada and Greenland. Samples of other components of the aerosol were collected and measurements were made of other properties of the aerosol coincident in time with the IN samples. The IN filters were exposed to water saturation in a dynamic developing chamber at –15° C and –25° C. Ice crystals grew on the IN and were counted on the filters at discrete time intervals during the exposure period to determine the rate of ice nucleation and the final concentration of (IN). Results show that Arctic haze aerosol, identified by pollutant signatures, had lower IN concentrations, a lower IN to total aerosol fraction and slower ice nucleation rates than aerosol which had a chemical signature more indicative of the remote unpolluted troposphere. These observations suggest that the Arctic haze aerosol does not efficiently form ice in the arctic troposphere. This may be a factor contributing to the long-range transport of Arctic haze.     

Boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice

During the field experiment ARKTIS 1993 ten cases of boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice in the Spitsbergen region were observed by aircraft over a distance ranging from about 50 km over the ice to about 300 km over the water. The modification depends decisively on the initial conditions over the ice, the boundary conditions at the bottom and top of the boundary layer and on the conditions of the large-scale flow. The modification of the bulk boundary-layer characteristics in relation to these conditions is presented.Besides the air-sea temperature contrast, the most important role for the boundary-layer modification is played by the stability on top of the boundary layer and by the divergence of the large-scale flow. According to the high variability of these conditions the observed boundary-layer modifications were very variable ranging from 100 to 300 m thick boundary layers with air temperatures between -32 and -22 °C over the ice to thicknesses between 900 and 2200 m and air temperatures between -15 and -5 °C after 300 km fetch over the open water. In most cases the large-scale flow was anticyclonic and divergent over the ice and changed to cyclonic and convergent over the water and an ice-sea breeze was superimposed on it.The sensible and latent heat fluxes are the dominant terms in the surface energy budget over the open water and ranged between 200 and 700 W m^-2 whereas the net longwave radiation is the dominating term over the ice with the heat fluxes only about 10 W m^-2.     

Projections of Climate Change Effects on Water Temperature Characteristics of Small Lakes in the Contiguous U.S.

To simulate effects of projected climate change on water temperature characteristics of small lakes in the contiguous U.S., a deterministic, one-dimensional year-round water temperature model is applied. In cold regions the model simulates ice and snow cover on a lake. The lake parameters required as model input are surface area, maximum depth, and Secchi depth as a measure of radiation attenuation and trophic state. The model is driven by daily weather data. Weather records from 209 stations in the contiguous U.S. for the period 1961–1979 were used to represent present climate conditions. The projected climate change owing to a doubling of atmospheric CO₂ was obtained from the output of the Canadian Climate Center General Circulation Model. The simulated water temperature and ice characteristics are related to the geometric and trophic state lake characteristics and to geographic location. By interpolation, the sensitivity of lake water temperature characteristics to latitude, longitude, lake geometry and trophic status can therefore be quantified for small lakes in the contiguous U.S. The 2× CO₂ climate scenario is projected to increase maximum and minimum lake surface temperatures by up to 5.2°C. (Maximum surface water temperatures in lakes near the northern and the southern border of the contiguous U.S. currently differ by up to 13°C.) Maximum temperature differences between lake surface and lake bottom are projected to increase in average by only 1 to 2°C after climate warming. The duration of seasonal summer stratification is projected to be up to 66 days longer under a 2×CO₂ climate scenario. Water temperatures of less than 8°C are projected to occur on lake bottoms during a period which is on the order of 50 days shorter under a 2×CO₂ climate scenario. With water temperature change projected to be as high as 5.2°C, ecological impacts such as shifts in species distributions and in fish habitat are most likely. Ice covers on lakes of northern regions would also be changed strongly.     

On the role of high latitude ice,snow, and vegetation feedbacks in the climatic response to external forcing changes

A seasonal energy balance climate model containing a detailed treatment of surface and planetary albedo, and in which seasonally varying land snow and sea ice amounts are simulated in terms of a number of explicit physical processes, is used to investigate the role of high latitude ice, snow, and vegetation feedback processes. Feedback processes are quantified by computing changes in radiative forcing and feedback factors associated with individual processes. Global sea ice albedo feedback is 5–8 times stronger than global land snowcover albedo feedback for a 2% solar constant increase or decrease, with Southern Hemisphere cryosphere feedback being 2–5 times stronger than Northern Hemisphere cryosphere feedback.In the absence of changes in ice extent, changes in ice thickness in response to an increase in solar constant are associated with an increase in summer surface melting which is exactly balanced by increased basal winter freezing, and a reduction in the upward ocean-air flux in summer which is exactly balanced by an increased flux in winter, with no change in the annual mean ocean-air flux. Changes in the mean annual ocean-air heat flux require changes in mean annual ice extent, and are constrained to equal the change in meridional oceanic heat flux convergence in equilibrium. Feedback between ice extent and the meridional oceanic heat flux obtained by scaling the oceanic heat diffusion coefficient by the ice-free fraction regulates the feedback between ice extent and mean annual air-sea heat fluxes in polar regions, and has a modest effect on model-simulated high latitude temperature change.Accounting for the partial masking effect of vegetation on snow-covered land reduces the Northern Hemisphere mean temperature response to a 2% solar constant decrease or increase by 20% and 10%, respectively, even though the radiative forcing change caused by land snowcover changes is about 3 times larger in the absence of vegetational masking. Two parameterizations of the tundra fraction are tested: one based on mean annual land air temperature, and the other based on July land air temperature. The enhancement of the mean Northern Hemisphere temperature response to solar constant changes when the forest-tundra ecotone is allowed to shift with climate is only 1/3 to 1/2 that obtained by Otterman et al. (1984) when the mean annual parameterization is used here, and only 1/4 to 1/3 as large using the July parameterization.The parameterized temperature dependence of ice and snow albedo is found to enhance the global mean temperature response to a 2% solar constant increase by only 0.04 °C, in sharp contrast to the results of Washington and Meehl (1986) obtained with a mean annual model. However, there are significant differences in the method used here and in Washington and Meehl to estimate the importance of this feedback process. When their approach is used in a mean annual version of the present model, closer agreement to their results is obtained.     

General circulation model simulation of mild nuclear winter effects

The climatic effects of an elevated uniform global layer of purely absorbing smoke of absorption optical depth 0.2 have been simulated using a version of the 9-level spectral model of McAvaney et al. (1978). The model was run at rhomboidal wave number 21 with convective adjustment, prognostic precipitation and soil hydrology, but fixed zonally averaged climatological cloud and fixed sea surface temperature, for constant January and July conditions with and without smoke absorption. Results show a reduction in convective rainfall in the tropics and monsoonal regions of the order of 50%, with diurnal average soil surface coolings of several degrees C except in those locations where the reduction in soil moisture is sufficient to effectively stop evaporation at the surface. In that case, small increases in temperature may occur. Results over Australia are consistent with the zonal mean picture. Run in a diurnal cycle mode, the model shows that daily maximum temperatures are more strongly affected, with soil surface coolings of the order of 2°–3° C in summer (with some local warmings) and 4°–6° C in winter. Overninght minimum temperatures cool by only 1°–2° C in both summer and winter. Possible effects of a lowering of sea surface temperature, variations in cloud cover, neglect of scattering by smoke, and infrared absorption and emission by the smoke are discussed.     

Glacial meltwater cooling of the Gulf of Mexico: GCM implications for Holocene and present-day climates

An atmospheric general circulation model, the NCAR CCM, has been used to investigate the possible effects that reduced Gulf of Mexico sea surface temperatures (SST) could have on regional and hemispheric climates. ¹⁸O records and terrestrial evidence indicate at least two major glacial meltwater discharges into the Gulf of Mexico subsequent to the last glacial maximum. It is probable that these discharges reduced Gulf of Mexico SST. We have conducted three numerical experiments, with imposed gulf-wide SST coolings of 3°C, 6°C, and 12°C, and find in all three experiments significant reductions in the North Atlantic storm-track intensity, along with a strong decrease in transient eddy water vapor transport out of the Gulf of Mexico. Surface pressures are higher over the North Atlantic, indicating a reduction of the climatological Icelandic low. The region is generally cooler and drier, with a reduction in precipitation that agrees well with evidence from Greenland ice cores. Other statistically significant changes occur across the Northern Hemisphere, but vary between the three experiments. In particular, warmer, wetter conditions are found over Europe for both the 6°C and 12°C SST reductions, but cooler conditions are found for the 3°C reduction. This indicates a dependence, in both the sign and magnitude of the model response, on the magnitude of the imposed SST anomaly. The results suggest that the present-day North Atlantic storm track is dependent on warm Gulf of Mexico SST for much of its intensity. They also suggest that meltwater-induced coolings may help account, in part, for some of the climatic oscillations that occurred during the last glacial/interglacial transition.     

Ice-albedo feedback in a CO2-doubling simulation

We estimate the feedback of sea-ice change to the warming from CO₂-doubling according to the simulation of Washington and Meehl (1984). Without ice-snow albedo feedback, their global warming of 3.5 °C would have been 2.2. °C according to our estimate of the ice-snow feedback. About 80% of the albedo change from ice and snow occurred in the Southern Hemisphere. Whether this change was an overestimate will require further study.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A regional model for studies of atmosphere-ice-ocean interaction in the western Arctic

Summary Asa step in the development of a fully coupled regional model of the atmosphere-ice-ocean system, atmospheric and sea ice models have been adapted to a western Arctic domain centered on the Bering Strait. Lateral boundary conditions derived from operational analyses drive the models through simulations on grids having horizontal resolutions of 21 km and 7 km. Sensitivities to the presence of sea ice are large after only 48 hours, by which time the surface temperatures in the Bering and Chukchi Seas are 10–15°C higher without sea ice than with sea ice. The temperatures, in turn, modify the fields of sea level pressure, surface wind and precipitation. By influencing the surface wind stress through the static static stability, the surface state feeds back to the surface momentum exchange, ice/ocean transport, and the rate of formation of new ice. The results also show a resolution-dependence of the surface winds, precipitation rates and new ice formation rates, particularly in areas in which the coastal configuration and topography are spatially complex. The experiments will be augmented by the implementation of an ocean model on the same grids.With 12 Figures     

Simulating the Holocene climate evolution at northern high latitudes using a coupled atmosphere-sea ice-ocean-vegetation model

The response of the climate at high northern latitudes to slowly changing external forcings was studied in a 9,000-year long simulation with the coupled atmosphere-sea ice-ocean-vegetation model ECBilt-CLIO-VECODE. Only long-term changes in insolation and atmospheric CO₂ and CH₄ content were prescribed. The experiment reveals an early optimum (9–8 kyr BP) in most regions, followed by a 1–3°C decrease in mean annual temperatures, a reduction in summer precipitation and an expansion of sea-ice cover. These results are in general agreement with proxy data. Over the continents, the timing of the largest temperature response in summer coincides with the maximum insolation difference, while over the oceans, the maximum response is delayed by a few months due to the thermal inertia of the oceans, placing the strongest cooling in the winter half year. Sea ice is involved in two positive feedbacks (ice-albedo and sea-ice insulation) that lead regionally to an amplification of the thermal response in our model (7°C cooling in Canadian Arctic). In some areas, the tundra-taiga feedback results in intensified cooling during summer, most notably in northern North America. The simulated sea-ice expansion leads in the Nordic Seas to less deep convection and local weakening of the overturning circulation, producing a maximum winter temperature reduction of 7°C. The enhanced interaction between sea ice and deep convection is accompanied by increasing interannual variability, including two marked decadal-scale cooling events. Deep convection intensifies in the Labrador Sea, keeping the overall strength of the thermohaline circulation stable throughout the experiment.     

Mid latitude winter climate variability in the South Indian and southwest Pacific regions since 1300 AD

Mid-latitude winter atmospheric variability in the South Indian Ocean and southwest Pacific Ocean regions of the circum-Antarctic are reconstructed using sea-salt aerosol concentrations measured in the high resolution Law Dome (DSS) ice core from East Antarctica. The sea-salt aerosol concentration data, as sodium (Na), were measured at approximately monthly resolution spanning the past 700 years. Analyses of covariations between Na concentrations in Law Dome ice, and mean sea-level pressure (MSLP) and wind field data were conducted to define the mid-latitude and sub-Antarctic atmospheric circulation patterns associated with variations in Na delivery. High Na concentrations in Law Dome snow are associated with increased meridional aerosol transport from mid-latitude sources. The seasonal average Na concentration for early winter (May, June, July (MJJ)) is strongly correlated to the mid-latitude MSLP field in the South Indian and southwest Pacific Oceans, and southern Australian regions. In addition, the average MJJ Na concentrations display a strong association with the stationary Rossby wave number 3 circulation, and are anti-correlated to the Southern Annular Mode (SAM) index of climate variability: high (low) Na concentrations occurring during negative (positive) SAM phases. This observed relationship is used to derive a proxy record for early-winter MSLP anomalies and the SAM in the South Indian and southwest Pacific Ocean regions over the period 1300–1995 AD. The proxy SAM index from 1300 to 1995 AD shows pronounced decadal-scale variability throughout. The period after 1500 AD is marked by a tendency toward slower variations and a weakly-positive mean SAM (enhanced westerlies in the 50° to 65°S zone) compared to the early part of the record.     

Satellite-derived outgoing longwave radiation,surface temperature and sea ice concentration off the coast of Adélie-Land,Eastern Antarctica

Summary Outgoing longwave radiation (OLR) as seen from satellite (NOAA-series), and sea ice concentration in Mer Dumont d'Urville, Eastern Antarctica were analyzed. For the time period 1974–1990 available radiative flux data showed a slight increase of 1.8 W/m² or 1% for the period. If thistrend should continue — trends over a 16-year period in any geophysical data is a rather questionable concept in isolation — a 4°C warming would be observed from space for this polar region over a century. The observed increase is, however, in agreement with Dumont d'Urville, the only ground station within the study area, which displayed a similar temperature increase (Periard and Pettré, 1991). Further it is in agreement with the general temperature increase which has been observed for the high southern latitudes (Boden et al., 1990). In addition, models of climatic change due to increased CO₂ and other trace gases predict for polar regions values of similar size.Sea ice concentration showed a slight decrease for the time period 1974–1989, for which data were available. However, a relationship existed between the radiative flux and the ice concentration, not only for the actual data, but also for the deviation series with the annual cycles removed. A correlation factor of –0.74 was found; the sensitivity for an increase of 10% in ice concentration was –2.9 W/m². This represents a temperature change as seen from space of about 1°C for a 10% change in sea ice concentration.With 7 Figures     

Simulated Climate Change Effects on Year-Round Water Temperatures in Temperate Zone Lakes

A deterministic, one-dimensional model is presented to simulate daily water temperature profiles and associated ice and snow covers for dimictic and polymictic lakes of the temperate zone. The lake parameters required as model input are surface area (A_s), maximum depth (H_MAX), and Secchi depth (z_s), the latter, used as a measure of light attenuation and trophic state. The model is driven by daily weather data and operates year-round over multiple years. The model has been tested with extensive data (over 5,000 temperature points). Standard error between simulated and measured water temperatures is 1.4°C in the open water season and 0.5°C in the ice cover season. The model is applied to simulate the sensitivity of Minnesota lake water temperature characteristics to climate change. The projected climate changes due to a doubling of atmospheric CO₂ are obtained from the output of the Canadian Climate Center General Circulation Model (CCC GCM) and the Goddard Institute of Space Studies General Circulation Model (GISS GCM). Simulated lake temperature characteristics have been plotted in a coordinate system with a lake geometry ratio (A _s ^0.25 /H_MAX) on one axis and Secchi depth on the other. The lake geometry ratio expresses a lake's susceptibility to stratification. By interpolation, the sensitivity of lake temperature characteristics to changes of water depth and Secchi depth under the projected climate scenarios can therefore be obtained. Selected lake temperature characteristics simulated with past climate conditions (1961–1979) and with a projected 2 × CO₂ climate scenario as input are presented herein in graphical form. The simulation results show that under the 2 × CO₂ climate scenario ice formation is delayed and ice cover period is shortened. These changes cause water temperature modifications throughout the year.     

Simulation of lake ice and its effect on the late-Pleistocene evaporation rate of Lake Lahontan

A model of lake ice was coupled with a model of lake temperature and evaporation to assess the possible effect of ice cover on the late-Pleistocene evaporation rate of Lake Lahontan. The simulations were done using a data set based on proxy temperature indicators and features of the simulated late-Pleistocene atmospheric circulation over western North America. When a data set based on a mean-annual air temperature of 3° C (7° C colder than present) and reduced solar radiation from jet-stream induced cloud cover was used as input to the model, ice cover lasting 4 months was simulated. Simulated evaporation rates (490–527 mm a^–1) were 60% lower than the present-day evaporation rate (1300 mm a^–1) of Pyramid Lake. With this reduced rate of evaporation, water inputs similar to the 1983 historical maxima that occurred in the Lahontan basin would have been sufficient to maintain the 13.5 ka BP high stand of Lake Lahontan.