Recent advances in understanding ice sheet dynamics

Glaciers and ice sheets play a dynamic role in Earth's climate system, influencing regional- and global-scale climate and responding to climate change on time scales from years to millennia. They are also an integral part of Earth's landscape in alpine and polar regions, where they are an active agent in isostatic, tectonic, and Earth surface processes. This review paper summarizes recent progress in understanding and modelling ice sheet dynamics, from the microphysical processes of ice deformation in glaciers to continental-scale processes that influence ice dynamics. Based on recent insights and research directions, it can be expected that a new generation of ice sheet models will soon replace the current standard. Improvements that can be foreseen in the near future include: (i) the addition of internally-consistent evolutionary equations for ice crystal fabric (anisotropic flow laws), (ii) more generalized flow laws that include different deformation mechanisms under different stress regimes, (iii) explicit incorporation of the effects of chemical impurities and grain size (dynamic recrystallization) on ice deformation, (iv) higher-order stress solutions to the momentum balance (Stokes' equation) that governs ice sheet flow, and (v) the continued merger of ice sheet models with increasingly complex Earth systems models, which include fully-coupled subglacial hydrological and geological processes. Examples from the Greenland Ice Sheet and Vatnajökull Ice Cap, Iceland are used to illustrate several of these new directions and their importance to glacier dynamics.     

Representing Grounding Line Dynamics in Numerical Ice Sheet Models: Recent Advances and Outlook

Recent satellite observations of the Antarctic and Greenland ice sheets show accelerated ice flow and associated ice sheet thinning along coastal outlet glaciers in contact with the ocean. Both processes are the result of grounding line retreat due to melting at the grounding line (the grounding line is the contact of the ice sheet with the ocean, where it starts to float and forms an ice shelf or ice tongue). Such rapid ice loss is not yet included in large-scale ice sheet models used for IPCC projections, as most of the complex processes are poorly understood. Here we report on the state-of-the art of grounding line migration in marine ice sheets and address different ways in which grounding line migration can be attributed and represented in ice sheet models. Using one-dimensional ice flow models of the ice sheet/ice shelf system we carried out a number of sensitivity experiments with different spatial resolutions and stress approximations. These are verified with semi-analytical steady state solutions. Results show that, in large-scale finite-difference models, grounding line migration is dependent on the numerical treatment (e.g. staggered/non-staggered grid) and the level of physics involved (e.g. shallow-ice/shallow-shelf approximation).     

Present State and Prospects of Ice Sheet and Glacier Modelling

Since the late 1970s, numerical modelling has become established as an important technique for the understanding of ice sheet and glacier dynamics, and several models have been developed over the years. Ice sheet models are particularly relevant for predicting the possible response of ice sheets to climate change. Recent observations suggest that ice dynamics could play a crucial role for the contribution of ice sheets to future sea level rise under global warming conditions, and the need for further research into the matter was explicitly stated in the Fourth Assessment Report (AR4) of the United Nations Intergovernmental Panel on Climate Change (IPCC). In this paper, we review the state of the art and current problems of ice sheet and glacier modelling. An outline of the underlying theory is given, and crucial processes (basal sliding, calving, interaction with the solid Earth) are discussed. We summarise recent progress in the development of ice sheet and glacier system models and their coupling to climate models, and point out directions for future work.     

Response of a Coupled Ocean�CAtmosphere Model to Greenland Ice Melting

We investigate the transient response of the global coupled ocean?Catmosphere system to enhanced freshwater forcing representative of melting of the Greenland ice sheets. A 50-year long simulation by a coupled atmosphere?Cocean general circulation model (CGCM) is compared with another of the same length in which Greenland melting is prescribed. To highlight the importance of coupled atmosphere?Cocean processes, the CGCM results are compared with those of two other experiments carried out with the oceanic general circulation model (OGCM). In one of these OGCM experiments, the prescribed surface fluxes of heat, momentum and freshwater correspond to the unperturbed simulation by the CGCM; in the other experiment, Greenland melting is added to the freshwater flux. The responses by the CGCM and OGCM to the Greenland melting have similar patterns in the Atlantic, albeit the former having five times larger amplitudes in sea surface height anomalies. The CGCM shows likewise stronger variability in all state variables in all ocean basins because the impact of Greenland melting is quickly communicated to all ocean basins via atmospheric bridges. We conclude that the response of the global climate to Greenland ice melting is highly dependent on coupled atmosphere?Cocean processes. These lead to reduced latent heat flux into the atmosphere and an associated increase in net freshwater flux into the ocean, especially in the subpolar North Atlantic. The combined result is a stronger response of the coupled system to Greenland ice sheet melting.     

Sensitivity of sea ice and ocean simulations to sea ice salinity in a coupled global climate model

The impacts of the spatiotemporal variations of sea ice salinity on sea ice and ocean characteristics have not been studied in detail, as the existing climate models neglect or misrepresent this process. To address this issue, this paper formulated a parameterization with more realistic sea ice salinity budget, and examined the sensitivity of sea ice and ocean simulations to the ice salinity variations and associated salt flux into the ocean using a coupled global climate model. Results show that the inclus...     

Sea-level rise: A review of recent past and near-future trends

Global mean sea level is a potentially sensitive indicator of climate change. Global warming will contribute to worldwide sea-level rise (SLR) from thermal expansion of ocean water, melting of mountain glaciers and polar ice sheets. A number of studies, mostly using tide-gauge data from the Permanent Service for Mean Sea Level, Bidston Observatory, England, have obtained rates of global SLR within the last 100 years that range between 0·3 and 3 mm yr^?1, with most values concentrated between 1 and 2 mm yr^?1. However, the reliability of these results has been questioned because of problems with data quality and physical processes that introduce a high level of spatial and temporal variability. Sources of uncertainty in the sea-level data include variations in winds, ocean currents, river runoff, vertical earth movements, and geographically uneven distribution of long-term records. Crustal motions introduce a major source of error. To a large extent, these can be filtered by employing palaeo-sea-level proxies, and geophysical modelling to remove glacio-isostatic changes. Ultimately, satellite geodesy will help resolve the inherent ambiguity between the land and ocean level changes recorded by tide gauges. Future sea level is expected to rise by ～ 1 m, with a ‘best-guess’ value of 48 cm by the year 2100. Such rates represent an acceleration of four to seven times over present rates. Local land subsidence could substantially increase the apparent SLR. For example, Louisiana is currently experiencing SLR trends nearly 10 times the global mean rate. These recently reduced SLR estimates are based on climate models that predict a zero to negative contribution to SLR from Antarctica. Most global climate models (GCMs) indicate an ice accumulation over Antarctica, because in a warmer world, precipitation will exceed ablation/snow-melt. However, the impacts of attritional processes, such as thinning of the ice shelves, have been downplayed according to some experts. Furthermore, not all climate models are in agreement. Opposite conclusions may be drawn from the results of other GCMs. In addition, the West Antarctic Ice Sheet is potentially subject to dynamic and volcanic instabilities that are difficult to predict. Because of the great uncertainty in SLR projections, careful monitoring of future sea-level trends by upgraded tide-gauge networks and satellite geodesy will become essential. Finally, because of the high spatial variability in crustal subsidence rates, wave climates and tidal regimes, it will be the set of local conditions (especially the relative sea-level rise), rather than a single global mean sea-level trend, that will determine each locality's vulnerability to future SLR.     

基于远震接收函数的南极大陆冰盖厚度研究

冰盖厚度是研究南极冰盖质量、建立冰盖动力学模型的基本参数,对于冰川均衡调整、冰盖物质平衡及全球气候变化研究具有重要意义.基于地震学的远震接收函数和H-Kappa格网搜索方法可以用于地震台站下方冰盖厚度的可靠探测,不仅能与冰雷达获得的冰盖厚度进行独立对比,还可以与冰雷达方法相互补充,进一步填补南极大陆冰盖厚度探测空白区.本文利用布设于南极大陆冰盖上方的流动地震台阵记录到的远震波形数据,基于接收函数方法对台阵下方的冰盖厚度进行了研究.结果显示:基于远震接收函数方法的冰盖厚度与Bedmap2冰厚格网模型相比,二者差别大多在200 m以内;少数台站差值达到600 m左右,这一差别可能与Bedmap2测线分布空区、冰雷达测深不确定性以及冰盖内部复杂波速结构等因素有关.本文研究结果表明:利用南极大陆冰盖上方的流动地震台阵,基于远震接收函数方法可以获得比较可靠的南极冰盖厚度,为独立验证冰雷达的探测结果并弥补冰雷达探测空白区提供了有效方法.同时,部分台站接收函数波形的复杂性可能暗示了南极大陆数千米厚的冰盖内部结构不是均一的,仍然存在比较复杂的内部结构变化.因此,有必要进一步利用包括接收函数波形拟合、地震面波反演等方法对南极大陆冰盖厚度及其内部精细结构进行更为深入的研究.     

Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling

We utilise a global finite-element sea ice–ocean model (FESOM), focused on the Antarctic marginal seas, to analyse projections of ice shelf basal melting in a warmer climate. Ice shelf–ocean interaction is described using a three-equation system with a diagnostic computation of temperature and salinity at the ice–ocean interface. A tetrahedral mesh with a minimumhorizontal resolution of 4 km and hybrid vertical coordinates is used. Ice shelf draft, cavity geometry, and global ocean bathymetry have been derived from the RTopo-1 data set. The model is forced with the atmospheric output from two climate models: (1) the Hadley Centre Climate Model (HadCM3) and (2) Max Planck Institute’s ECHAM5/MPI-OM coupled climate model. Results from experiments forced with their twentieth century output are used to evaluate the modelled present-day ocean state. Sea ice coverage is largely realistic in both simulations; modelled ice shelf basal melt rates compare well with observations in both cases, but are consistently smaller for ECHAM5/MPI-OM. Projections for future ice shelf basal melting are computed using atmospheric output for the Intergovernmental Panel on Climate Change (IPCC) scenarios E1 and A1B. In simulations forced with ECHAM5 data, trends in ice shelf basal melting are small. In contrast, decreasing convection along the Antarctic coast in HadCM3 scenarios leads to a decreasing salinity on the continental shelf and to intrusions of warm deep water of open ocean origin. In the case of the Filchner–Ronne Ice Shelf (FRIS), this water reaches deep into the cavity, so that basal melting increases by a factor of 4 to 6 compared to the present value of about 90 Gt/year. By the middle of the twenty-second century, FRIS becomes the dominant contributor to total ice shelf basal mass loss in these simulations. Our results indicate that the surface freshwater fluxes on the continental shelves may be crucial for the future of especially the large cold water ice shelves in the Southern Ocean.     

Oblique flexural-gravity wave scattering by a submerged semi-circular ridge

This study examines the flexural-gravity wave scattering by a semi-circular ridge submerged in the ocean covered by a thin ice sheet. The ice sheet is treated as an elastic plate and its deflection is modelled using the Euler-Bernoulli beam equation. The fluid motion is described by linear potential theory. The series solutions of velocity potentials for obliquely and normally incident waves are both developed using multipole expansions. The unknowns in the series solutions are determined in terms of the impermeable condition on the ridge surface. Accurate results for the reflection and transmission coefficients of ridge and the deflection of ice sheet are presented. The effects of ridge radius, ice thickness, wave frequency and wave incident angle on the hydrodynamic quantities are shown.     

Climatic effects on ice phenology and ice-jam flooding of the Athabasca River in western Canada

ABSTRACT

In cold region environments, any alteration in the hydro-climatic regime can have profound impacts on river ice processes. This paper studies the implications of hydro-climatic trends on river ice processes, particularly on the freeze-up and ice-cover breakup along the Athabasca River in Fort McMurray in western Canada, which is an area very prone to ice-jam flooding. Using a stochastic approach in a one-dimensional hydrodynamic river ice model, a relationship between overbank flow and breakup discharge is established. Furthermore, the likelihood of ice-jam flooding in the future (2041–2070 period) is assessed by forcing a hydrological model with meteorological inputs from the Canadian regional climate model driven by two atmospheric–ocean general circulation climate models. Our results show that the probability of ice-jam flooding for the town of Fort McMurray in the future will be lower, but extreme ice-jam flood events are still probable.     

Landscape preservation under Fennoscandian ice sheets determined from in situ produced Be and Al

Some areas within ice sheet boundaries retain pre-existing landforms and thus either remained as ice free islands (nunataks) during glaciation, or were preserved under ice. Differentiating between these alternatives has significant implications for paleoenvironment, ice sheet surface elevation, and ice volume reconstructions. In the northern Swedish mountains, in situ cosmogenic ¹⁰Be and ²⁶Al concentrations from glacial erratics on relict surfaces as well as glacially eroded bedrock adjacent to these surfaces, provide consistent last deglaciation exposure ages (∼8-13 kyr), confirming ice sheet overriding as opposed to ice free conditions. However, these ages contrast with exposure ages of 34-61 kyr on bedrock surfaces in these same relict areas, demonstrating that relict areas were preserved with little erosion through multiple glacial cycles. Based on the difference in radioactive decay between ²⁶Al and ¹⁰Be, the measured nuclide concentration in one of these bedrock surfaces suggests that it remained largely unmodified for a minimum period of 845₋₄₁₈⁺⁴⁶¹ kyr. These results indicate that relict areas need to be accounted for as frozen bed patches in basal boundary conditions for ice sheet models, and in landscape development models. Subglacial preservation also implies that source areas for glacial sediments in ocean cores are considerably smaller than the total area covered by ice sheets. These relict areas also have significance as potential long-term subglacial biologic refugia.     

E) Organisations Nationales

Abstract

Sudden rapid advances or surges of glaciers and sections of smail ice caps are well known. After remaining dormant or in retreat over long periods of time these ice masses suddenly move forward rapidly with speeds about 2 orders of magnitude greater than usual. If such surges were to occur in large sections of the Antarctic ice sheet serious consequences could result. These include a significant rise of sea level, a substantial increase in the high-albedo ice cover around the continent especially in summer, and a cooling of the Antarctic ocean by the additional ice melting.

A numerical model has now been developed which simulates surging in certain glaciers and ice sheets in an apparently realistic manner. This model has been found to give close representations to a number of existing real surging temperate ice masses from small mountain glaciers to large sectors of ice caps. The model reproduces realistically many features of these ice masses such as the period of the surge, the duration, the velocity of advance, the magnitude of the advance, and the changes in ice thickness.

The application of the model to the Antarctic ice sheet is made more difficult by the problem caused by the temperature dependence of the flow properties of ice. This means that for a complete study the interaction with the environment needs to be considered. However, at this stage preliminary calculations indicate a number of features that are relevant to the effect of Antarctic ice surges on the global climate. These include the period between surges, the duration of the surge, the amount of ice advanced and the changes in thickness of the ice sheet.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Toward Improved Estimation of the Dynamic Topography and Ocean Circulation in the High Latitude and Arctic Ocean: The Importance of GOCE

The Arctic plays a fundamental role in the climate system and shows significant sensitivity to anthropogenic climate forcing and the ongoing climate change. Accelerated changes in the Arctic are already observed, including elevated air and ocean temperatures, declines of the summer sea ice extent and sea ice thickness influencing the albedo and CO₂ exchange, melting of the Greenland Ice Sheet and increased thawing of surrounding permafrost regions. In turn, the hydrological cycle in the high latitude and Arctic is expected to undergo changes although to date it is challenging to accurately quantify this. Moreover, changes in the temperature and salinity of surface waters in the Arctic Ocean and Nordic Seas may also influence the flow of dense water through the Denmark Strait, which are found to be a precursor for changes in the Atlantic meridional overturning circulation with a lead time of around 10 years (Hawkins and Sutton in Geophys Res Lett 35:L11603, 2008). Evidently changes in the Arctic and surrounding seas have far reaching influences on regional and global environment and climate variability, thus emphasizing the need for advanced quantitative understanding of the ocean circulation and transport variability in the high latitude and Arctic Ocean. In this respect, this study combines in situ hydrographical data, surface drifter data and direct current meter measurements, with coupled sea ice–ocean models, radar altimeter data and the latest GOCE-based geoid in order to estimate and assess the quality, usefulness and validity of the new GOCE-derived mean dynamic topography for studies of the ocean circulation and transport estimates in the Nordic Seas and Arctic Ocean.     

Coherent variations of the obliquity components in global ice volume and ocean carbon reservoir over the past 5 Ma

The Milankovi theory stresses that the summer insolation in the high northern latitudes that is dominated by the precession cycle controls the glacial/interglacial cycles in global climate change.If the climate system responds linearly to the external insolation forcing,the precession cycle of 23 or 19 ka should dominate the variations in the climatic proxy records.I performed spectral and evolutive cross spectral analyses on the high resolution benthic 18O and 13C records from the South China Sea and the North Atlantic,the proxies of global ice volume and ocean carbon reservoir respectively.I found that the obliquity instead of the eccentricity or the precession is the most marked cycle in the global ice volume and ocean carbon reservoir variations over the past 5 Ma.The analysis further reveals that only at the obliquity band instead of the eccentricity or the precession band does the global ice volume and ocean carbon reservoir display consistently high coherency and stable phase relationship over the past 5 Ma.The consistently positive or near-zero phases of the benthic 18O relative to the benthic13C at the obliquity band suggest that the global carbon cycle is involved in the polar ice sheet growth as an important internal feedback,not a determinative driving factor.The obliquity instead of the precession or the eccentricity takes the dominant role of driving the global climate change during the Pliocene and Pleistocene.     

利用ICESat数据解算南极冰盖冰雪质量变化

南极冰盖冰雪质量变化反映了全球气候变化,并且直接影响着全球海平面变化.ICESat测高卫星的主要任务之一就是要确定南北两极冰盖的质量变化情况并评估其对全球海平面变化的影响.本文利用2003年10月至2008年12月的ICESat测高数据,针对南极DEM分辨率有限的特殊性,通过求解坡度改正值,解决重复轨道地面脚点不重合的问题,计算了南极大陆(86°S以北区域,后文所述南极冰盖均不包括86°S以南区域)在这5年里的冰雪质量变化情况,得到东南极冰盖的质量变化为-18±20Gt/a,西南极-26±6Gt/a,南极冰盖的冰雪质量变化为-44±21Gt/a,对全球海平面上升的影响约为0.12mm·a~(-1).解算结果表明,南极冰盖质量亏损主要集中在西南极阿蒙森海岸附近冰川以及东南极波因塞特角区域.     

Response of the Greenland and Antarctic Ice Sheets to Multi-Millennial Greenhouse Warming in the Earth System Model of Intermediate Complexity LOVECLIM

Calculations were performed with the Earth system model of intermediate complexity LOVECLIM to study the response of the Greenland and Antarctic ice sheets to sustained multi-millennial greenhouse warming. Use was made of fully dynamic 3D thermomechanical ice-sheet models bidirectionally coupled to an atmosphere and an ocean model. Two 3,000-year experiments were evaluated following forcing scenarios with atmospheric CO₂ concentration increased to two and four times the pre-industrial value, and held constant thereafter. In the high concentration scenario the model shows a sustained mean annual warming of up to 10°C in both polar regions. This leads to an almost complete disintegration of the Greenland ice sheet after 3,000 years, almost entirely caused by increased surface melting. Significant volume loss of the Antarctic ice sheet takes many centuries to initiate due to the thermal inertia of the Southern Ocean but is equivalent to more than 4 m of global sea-level rise by the end of simulation period. By that time, surface conditions along the East Antarctic ice sheet margin take on characteristics of the present-day Greenland ice sheet. West Antarctic ice shelves have thinned considerably from subshelf melting and grounding lines have retreated over distances of several 100 km, especially for the Ross ice shelf. In the low concentration scenario, corresponding to a local warming of 3?C4°C, polar ice-sheet melting proceeds at a much lower rate. For the first 1,200 years, the Antarctic ice sheet is even slightly larger than today on account of increased accumulation rates but contributes positively to sea-level rise after that. The Greenland ice sheet loses mass at a rate equivalent to 35 cm of global sea level rise during the first 1,000 years increasing to 150 cm during the last 1,000 years. For both scenarios, ice loss from the Antarctic ice sheet is still accelerating after 3,000 years despite a constant greenhouse gas forcing after the first 70?C140 years of the simulation.     

LASG耦合气候系统模式FGCM-1.0

本文描述了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG)最新发展的一个耦合气候系统模式的基本性能. 该模式是在LASG灵活的全球耦合气候系统模式（英文缩写为FGCM）的初始版本FGCM-0的基础上发展而来的,是该系列耦合模式的第二个版本,即FGCM-1.0. FGCM-1.0通过一个通量耦合器将大气、海洋和海冰三个分量模式耦合在一起,其中海洋分量模式是LASG发展的一个涡相容分辨率（eddy-permitting）全球海洋环流模式,大气和海冰分量模式则为美国国家大气研究中心（NCAR）的大气环流模式CAM2和海冰模式CSIM4. 耦合模式完整地考虑了海气界面上的动量、热量和淡水通量交换,尽管在模式中没有使用任何形式的人为的通量调整或者通量距平方案,模式还是比较合理地模拟出基本的气候形态. 通过对该耦合模式长期积分结果的进一步分析发现,模式能够比较好地模拟出厄尔尼诺－南方涛动（ENSO）以及印度洋偶极子事件的基本特征;与FGCM系列耦合模式的最初版本FGCM-0相比,FGCM-1.0模拟的北赤道逆流（NECC）和ENSO循环更加真实.     

利用CryoSat-2波形数据建立南极Lambert冰川流域DEM

Lambert冰川-Amery冰架系统是南极冰盖最大的冰流系统之一,对南极冰盖物质平衡研究有着重要的作用.数字高程模型(DEM)是进行南极冰盖研究的基础.本文基于CryoSat-2 L1b波形数据,研究建立了Lambert冰川流域高分辨率DEM.测高卫星返回波形在冰盖区域存在变形,需进行波形重跟踪处理.利用交叉点分析方法对重心偏移法(OCOG)、阈值法和β参数法等常用的波形重跟踪方法对不同类型的CryoSat-2波形的适用性进行了研究.最后,利用克里金插值方法建立了500 m分辨率的Lambert冰川流域DEM——LAS DEM (Lambert Glacier-Amery Ice Shelf system DEM).利用ICESat卫星测高数据和GPS地面实测数据对LAS DEM进行精度验证,并与另外两种基于CryoSat-2数据的南极DEM进行了对比.结果表明:LAS DEM的整体精度约为0.295±2.7 m,优于另外两种CryoSat-2 DEM;在冰盖内陆地区,LAS DEM的高程误差在2 m之内;在Amery冰架上,LAS DEM的精度优于1 m.     

Sensitivity of physical lake processes to climate change within a large Precambrian Shield catchment

Potential future changes in lake physical processes (e.g. stratification and freezing) can be assessed through exploring their sensitivity to climate change, and assessing the current vulnerability of different lake types to plausible changes in meteorological drivers. This study quantifies the impacts of climate change and sensitivity of lake physical processes within a large (5100 km²) Precambrian Shield catchment in south‐central Ontario. Historic regional relationships are established between climate drivers, lake morphology, and lake physical changes through generalized linear modelling (GLM), and are used to quantify likely changes in timing of ice phenology and lake stratification across 72 lakes under a range of future climate models and scenarios. In response to projections of increased temperature (ensemble mean of +3.3 °C), both earlier ice‐off and onset of summer stratification were projected, with later ice‐on and fall turnover compared to the baseline. Process sensitivity to climate change varied by lake type; shallower lakes with a smaller volume (less than 15 m deep and less than 0.05 km³) were more sensitive to processes associated with lake heating (stratification onset and ice‐off), and deeper lakes with a larger surface area (greater than 30 m deep and greater than 1000 ha) were more sensitive to processes associated with lake cooling (fall turnover and ice‐on). These results indicate that whereas small lakes are vulnerable to climate warming because of changes that occur in spring and summer, larger lakes are particularly sensitive during the fall. The findings suggest that lake morphology and associated sensitivity should be considered in the development of sustainable lake management strategies. Copyright © 2016 John Wiley & Sons, Ltd.