Estimation of snow water equivalent using microwave radiometry over Arctic first‐year sea ice

The magnitude and spatial distribution of snow on sea ice are both integral components of the ocean–sea‐ice–atmosphere system. Although there exists a number of algorithms to estimate the snow water equivalent (SWE) on terrestrial surfaces, to date there is no precise method to estimate SWE on sea ice. Physical snow properties and in situ microwave radiometry at 19, 37 and 85 GHz, V and H polarization were collected for a 10‐day period over 20 first‐year sea ice sites. We present and compare the in situ physical, electrical and microwave emission properties of snow over smooth Arctic first‐year sea ice for 19 of the 20 sites sampled. Physical processes creating the observed vertical patterns in the physical and electrical properties are discussed. An algorithm is then developed from the relationship between the SWE and the brightness temperature measured at 37 GHz (55°) H polarization and the air temperature. The multiple regression between these variables is able to account for over 90% of the variability in the measured SWE. This algorithm is validated with a small in situ data set collected during the 1999 field experiment. We then compare our data against the NASA snow thickness algorithm, designed as part of the NASA Earth Enterprise Program. The results indicated a lack of agreement between the NASA algorithm and the algorithm developed here. This lack of agreement is attributed to differences in scale between the Special Sensor Microwave/Imager and surface radiometers and to differences in the Antarctic versus Arctic snow physical and electrical properties. Copyright © 2003 John Wiley & Sons, Ltd.     

On the winter evolution of snow thermophysical properties over land‐fast first‐year sea ice

The geophysical, thermodynamic and dielectric properties of snow are important state variables that are known to be sensitive to Arctic climate variability and change. Given recent observations of changes in the Arctic physical system (Arctic Climate Impact Assessment, 2004), it is important to focus on the processes that give rise to variability in the horizontal, vertical and temporal dimensions of the life‐history of snow on sea ice. The objectives in this study are to present these ‘state’ variables and to investigate the processes that govern variability in the vertical, horizontal and temporal dimension by using a case study over land‐fast first‐year sea ice for the period December 2003 to June 2004. Results from two sampling areas (thin and thick snowpacks) show that differences in snowpack thickness can substantially change the vertical and temporal evolution of snow properties. During the late fall and early winter (cooling period) we measured no significant changes in the physical properties, except for thin snow‐cover salinity, which decreased throughout the period. Fall‐snow desalination was only observed under thin snowpacks with a rate of ?0·12 ppt day^?1. Significant changes occurred in the late winter and early spring (warming period), especially for snow grain size. Snow grain kinetic growth of 0·25–0·48 mm·day^?1 was measured coincidently with increasing salinity and wetness for both thin and thick snowpacks. Copyright © 2006 John Wiley & Sons, Ltd.     

New satellite derived sea ice motion tracks Arctic contamination

Sea ice has been reported to contain contaminants from atmospheric and nearshore sediment resuspension processes. In this study successive passive microwave images from the 85.5 GHz channels on the Special Sensor Microwave Imager (SSM/I) were merged with drifting buoy trajectories from the International Arctic Buoy Program to compute Arctic sea ice motion in the Russian Arctic between 1988 and 1994. Smooth daily motion fields were averaged to prepare monthly maps making it possible to compute the 7-year mean and mean seasonal ice motions as well as principal components of directional variability of sea ice motion for the entire Arctic and surrounding basins. These mean motion vectors are used to simulate the advection of contaminants deposited on or contained within the sea ice and subsequently transported into the Arctic Ocean in order to predict both their mean trajectories and dispersal over time. The 3-year displacement of contaminants from a number of Russian sites and one American site display various behaviours from substantial displacement and dispersal to almost no movement. This computational procedure could be applied to realtime SSM/I and ice buoy data to provide detailed, all-weather, vector motion maps of ice circulation to predict the path and dispersal of any new substance introduced to the sea ice and transported into the Arctic or Antarctic ocean surface.     

Sea ice thickness measurement and its underside morphology analysis using radar penetration in the Arctic Ocean

Sea ice, as an important component of the Arctic climate system, has drawn significant sci-entific interest. Sea ice thickness and its morphology have dramatic impacts on ocean-atmos- phere-ice interactions[1—4], which directly affect the exchange proces…     

Hydrogeomorphic processes of thermokarst lakes with grounded‐ice and floating‐ice regimes on the Arctic coastal plain,Alaska

Thermokarst lakes cover > 20% of the landscape throughout much of the Alaskan Arctic Coastal Plain (ACP) with shallow lakes freezing solid (grounded ice) and deeper lakes maintaining perennial liquid water (floating ice). Thus, lake depth relative to maximum ice thickness (1·5–2·0 m) represents an important threshold that impacts permafrost, aquatic habitat, and potentially geomorphic and hydrologic behaviour. We studied coupled hydrogeomorphic processes of 13 lakes representing a depth gradient across this threshold of maximum ice thickness by analysing remotely sensed, water quality, and climatic data over a 35‐year period. Shoreline erosion rates due to permafrost degradation ranged from < 0·2 m/year in very shallow lakes (0·4 m) up to 1·8 m/year in the deepest lakes (2·6 m). This pattern of thermokarst expansion masked detection of lake hydrologic change using remotely sensed imagery except for the shallowest lakes with stable shorelines. Changes in the surface area of these shallow lakes tracked interannual variation in precipitation minus evaporation (P ? E_L) with periods of full and nearly dry basins. Shorter‐term (2004–2008) specific conductance data indicated a drying pattern across lakes of all depths consistent with the long‐term record for only shallow lakes. Our analysis suggests that grounded‐ice lakes are ice‐free on average 37 days longer than floating‐ice lakes resulting in a longer period of evaporative loss and more frequent negative P ? E_L. These results suggest divergent hydrogeomorphic responses to a changing Arctic climate depending on the threshold created by water depth relative to maximum ice thickness in ACP lakes. Copyright © 2011 John Wiley & Sons, Ltd.     

Re-calibration of Arctic sea ice extent datasets using Arctic surface air temperature records

A new seasonal and annual dataset describing Arctic sea ice extents for 1901–2015 was constructed by individually re-calibrating sea ice data sources from the three Arctic regions (North American, Nordic and Siberian) using the corresponding surface air temperature trends for the pre-satellite era (1901–1978), so that the strong relationship between seasonal sea ice extent and surface air temperature observed for the satellite era (1979-present) also applies to the pre-satellite era. According to this new dataset, the recent period of Arctic sea ice retreat since the 1970s followed a period of sea ice growth after the mid-1940s, which in turn followed a period of sea ice retreat after the 1910s. Arctic sea ice is a key component of the Arctic hydrological cycle, through both its freshwater storage role and its influence on oceanic and atmospheric circulation. Therefore, these new insights have significance for our understanding of Arctic hydrology.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR not assigned     

An examination of snow redistribution over smooth land‐fast sea ice

An understanding of temporal evolution of snow on sea ice at different spatial scales is essential for improvement of snow parameterization in sea ice models. One of the problems we face, however, is that long‐term climate data are routinely available for land and not for sea ice. In this paper, we examine the temporal evolution of snow over smooth land‐fast first‐year sea ice using observational and modelled data. Changes in probability density functions indicate that depositional and drifting events control the evolution of snow distribution. Geostatistical analysis suggests that snowdrifts increased over the study period, and the orientation was related to the meteorological conditions. At the microscale, the temporal evolution of the snowdrifts was a product of infilling in the valleys between drifts. Results using two shore‐based climate reporting stations (Paulatuk and Tuktoyuktuk, NWT) suggest that on‐ice air temperature and relative humidity can be estimated using air temperature recorded at either station. Wind speed, direction and precipitation on ice cannot be accurately estimated using meteorological data from either station. The temporal evolution of snow distribution over smooth land‐fast sea ice was modelled using SnowModel and four different forcing regimes. The results from these model runs indicate a lack of agreement between observed distribution and model outputs. The reasons for these results are lack of meteorological measurements prior to the end of January, lack of spatially adequate surface topography and discrepancies between meteorological variables on land and ice. Copyright © 2009 John Wiley & Sons, Ltd.     

First‐year sea ice spring melt transitions in the Canadian Arctic Archipelago from time‐series synthetic aperture radar data, 1992–2002

This paper synthesizes 10‐years' worth of interannual time‐series space‐borne ERS‐1 and RADARSAT‐1 synthetic aperture radar (SAR) data collected coincident with daily measurement of snow‐covered, land‐fast first‐year sea ice (FYI) geophysical and surface radiation data collected from the Seasonal Sea Ice Monitoring and Modeling Site, Collaborative‐Interdisciplinary Cryospheric Experiment and 1998 North Water Polynya study over the period 1992 to 2002. The objectives are to investigate the seasonal co‐relationship of the SAR time‐series dataset with selected surface mass (bulk snow thickness) and climate state variables (surface temperature and albedo) measured in situ for the purpose of measuring the interannual variability of sea ice spring melt transitions and validating a time‐series SAR methodology for sea ice surface mass and climate state parameter estimation. We begin with a review of the salient processes required for our interpretation of time‐series microwave backscatter from land‐fast FYI. Our results suggest that time‐series SAR data can reliably measure the timing and duration of surface albedo transitions at daily to weekly time‐scales and at a spatial scales that are on the order of hundreds of metres. Snow thickness on FYI immediately prior to melt onset explains a statistically significant portion of the variability in timing of SAR‐detected melt onset to pond onset for SAR time‐series that are made up of more than 25 images. Our results also show that the funicular regime of snowmelt, resolved in time‐series SAR data at a temporal resolution of approximately 2·5 images per week, is not detectable for snow covers less than 25 cm in thickness. Copyright © 2007 John Wiley & Sons, Ltd.     

秋季北极海冰与欧亚冬季气温在年代际和年际尺度上的不同联系

北极海冰的急剧消融在近年来欧亚大陆频发的低温事件中扮演着关键角色.秋季北极海冰的偏少对应着冬季欧亚大陆的低温天气,然而二者的联系在年代际和年际两种时间尺度上存在显著区别.本文运用1979—2012年哈德莱中心第一套海冰覆盖率(HadISST1)、欧洲中心(ERA_Interim)的2m温度、风场、海平面气压场、高度场等资料,分别研究了年代际和年际时间尺度上前期秋季北极海冰与欧亚冬季气温的联系.结果表明,欧亚和北极地区(0°—160°E,15°N—90°N)的冬季气温具有显著的年代际和年际变化.在年代际尺度上,温度异常分布在21世纪初由北极冷-大陆暖转为北极暖-大陆冷.这一年代际转折与前期秋季整个北极地区的海冰年代际减少联系密切.秋季北极全区海冰年代际偏少对应冬季欧亚大陆中高纬地区的高压异常,有利于北大西洋的暖湿气流北上和北极的冷空气南侵,造成北极暖-大陆冷的温度分布;在年际时间尺度上,温度异常分布主要由第一模态的年际变化部分和第二模态组成,且第一模态包含的年际变率信号也存在显著的年代际变化.年际尺度上全区北极海冰对欧亚冬季气温的影响远不及位于北冰洋西南边缘的巴伦支海、喀拉海和拉普捷夫海西部(30°E—120°E,75°N—85°N)的关键区海冰影响显著.关键区内海冰的偏少会引发冬季的北大西洋涛动负位相,导致北大西洋吹往欧亚大陆的暖湿气流减弱和欧亚大陆中高纬地区的气温偏低.     

Detection and classification of surface roughness in an Arctic marginal sea ice zone

Sea ice dynamic and thermodynamic processes are important and highly variable elements of the marginal ice zone (MIZ). This study examines the detection and classification of statistically separable sea ice classes in the MIZ through a range of temporal and spatial scales. A helicopter‐based laser system was used to obtain large‐scale and a ship‐based laser profiler to identify small‐scale roughness types, respectively. The analysis of variance of surface height data from helicopter‐ and ship‐based laser systems, active microwave (AMW) C‐band backscattering data and passive microwave (PMW) (37 and 89 GHz) brightness temperature data reveal different classes that statistically differ from one another. We found significant statistical difference in variances in AMW data with six classes that differ in VV polarization, three classes in VH polarization, and five classes in HH polarization in the MIZ (e.g. snow‐covered first‐year ice, ice rubble, pancake ice, frost flowers, melt pond, flooded ice, and ice edge) of southeastern Beaufort Sea. The PMW emission was not as effective at discrimination, yielding only one statistically separable class. The results can potentially be extended to satellite‐based investigations of the MIZ at regional scales. Copyright © 2012 John Wiley & Sons, Ltd.     

Seafloor evidence for ice shelf flow across the Alaska–Beaufort margin of the Arctic Ocean

New imagery of ～14 100 km² of seafloor along a 640 km stretch of the Alaska and Beaufort margins (ABM) in water depths from 250 to 2800 m depicts a repetitive association of glaciogenic bedforms (lineations and iceberg scours), broad erosional bathymetric features and adjacent downslope turbidite gullies. These bedforms have styles, depths and orientations similar to features discovered earlier on the Chukchi Borderland, up to 800 km northwest of the ABM. Lineations occur across the surface of a flattened bathymetric bench interpreted to have formed by an ice shelf sliding along the continental slope and scraping the seafloor at temporary grounding locations. The glacial geology of surrounding areas suggests that an ice shelf probably flowed from the mouths of overdeepened glacial troughs in the Canadian Arctic Archipelago westward along the ABM and across the Chukchi Borderland. This curved pathway indicates an obstruction to ice flow in the central Canada Basin, possibly caused by either a basin‐wide ice shelf or by a pile‐up of mega‐bergs originating from the Eurasian side of the Arctic Ocean. The ice shelf that affected the ABM may have formed between Oxygen Isotopic Stage 4 to 5b, possibly correlating to an inferred intra‐Stage 5 widespread Beringian glaciation. Evidence for glaciogenic features on the ABM corroborates suggestions that large ice volumes and extents existed in the Arctic during Pleistocene glacial periods. These findings have far‐reaching implications for Arctic climate studies, ocean circulation, sediment stratigraphy and the stability of circum‐Arctic continental ice masses. Copyright © 2007 John Wiley & Sons, Ltd.     

Dipole anomaly in the Arctic atmosphere and winter Arctic sea ice motion

Arctic sea ice, an important component of the cli- mate system, has received significant attention re- cently. Arctic sea ice variation is an important indica- tor of changes in the climate system, such as global change and polar amplification, and observation and climate modeling suggests that sea ice can itself be an agent of climate change[1―4]. Previous studies[3,5] have shown that sea ice influences the surface and atmos- pheric boundary layer temperature. Potentially more important is t…     

Sea ice as an evaluation target for HEM modelling and inversion

We investigate helicopter electromagnetic (HEM) inversion schemes applied to synthetic and measured HEM sea ice profiling data. Direct HEM data-to-ice-thickness inversion is compared to three different formal, least squares layered earth inversion algorithms.By making several approximations, it is possible to directly invert a single channel measurement (i.e., the in-phase or quadrature component of a single frequency measurement) to obtain an estimate of sea ice thickness. Measurements from multiple input channels, however, can be used in a layered earth inversion to simultaneously recover several model parameters such as sea ice thickness, sea ice conductivity and sub-ice bathymetry. Synthetic data sets for a particular two-frequency HEM system showed that simple least squares inversion algorithms produce reliable estimates of sea ice thickness in cases where the ice is thicker than 3 m. These methods could also recover acceptable estimates of sea ice thickness when a thin, conductive, partially melted sea ice layer was present, and could determine shallow, sub-ice bathymetry in brackish water. As expected, 1D transformations and inversions of synthetic data for a three-dimensional pressure ridge keel structure contained artifacts, notably broadening of the apparent width of the keel.Prior to inverting a field data set acquired over rather thin (~ 0.5 m) Antarctic sea ice, we found it necessary to recalibrate the phase angle of the measurements using a phasor diagram-based method. Direct transformation of a single channel from the recalibrated data set produced more accurate estimates of sea-ice thickness than formal inversion of multi-channel data. We suggest that the least squares inversion methods are inferior in this situation because of the particular characteristics of the two-frequency HEM system used in this evaluation; the extreme differences in sensitivity of high and low frequency data components, the overall low sensitivity to sea ice conductivity (especially for thin ice), and the partially low signal-to noise ratios of the measurements. The data sets used in this study will be made available to the public to allow alternate inversion approaches to be applied and evaluated. It is suggested that inclusion of parameter bounds and other forms of regularization could help to improve the inversion results.     

Ice processes on medium‐sized north‐temperate lakes

Lake ice supports a range of socio‐economic and cultural activities including transportation and winter recreational actives. The influence of weather patterns on ice‐cover dynamics of temperate lakes requires further understanding for determining how changes in ice composition will impact ice safety and the range of ecosystem services provided by seasonal ice cover. An investigation of lake ice formation and decay for three lakes in Central Ontario, Canada, took place over the course of two winters, 2015–2016 and 2016–2017, through the use of outdoor digital cameras, a Shallow Water Ice Profiler (upward‐looking sonar), and weekly field measurements. Temperature fluctuations across 0°C promoted substantial early season white ice growth, with lesser amounts of black ice forming later in the season. Ice thickening processes observed were mainly through meltwater, or midwinter rain, refreezing on the ice surface. Snow redistribution was limited, with frequent melt events limiting the duration of fresh snow on the ice, leading to a fairly uniform distribution of white ice across the lakes in 2015–2016 (standard deviations week to week ranging from 3 to 5 cm), but with slightly more variability in 2016–2017 when more snow accumulated over the season (5 to 11 cm). White ice dominated the end‐of‐season ice composition for both seasons representing more than 70% of the total ice thickness, which is a stark contrast to Arctic lake ice that is composed mainly of black ice. This research has provided the first detailed lake ice processes and conditions from medium‐sized north‐temperate lakes and provided important information on temperate region lake ice characteristics that will enhance the understanding of the response of temperate lake ice to climate and provide insight on potential changes to more northern ice regimes under continued climate warming.     

Spatio‐temporal analysis of glacial ice area distribution of Hunza River Basin,Karakoram region of Pakistan

There is a high degree of uncertainty about the state and fate of Pakistan's Karakoram glaciers due to data scarcity in high altitude regions. They are thought to be less vulnerable to climatic change because they behave differently as compared with eastern Himalayas. This study measures the decadal temporal changes in the glacial ice area of Karakoram's Hunza River Basin, one of the eight subbasins of Upper Indus Basin. An attempt has been made to investigate the relationship between glacial ice area changes and calculated values of precipitation, temperature and run‐off. A combination of satellite and field‐based approach is applied. Output includes maps of glacial ice hypsometries of eight glacial ice subregions of Hunza River Basin for 3 years (i.e., 1989, 2002, and 2010). The results show a decreasing trend in the glacial ice‐covered area signifying a reduction of 20.47% with the largest reduction being in the lower elevation bands. There is presently no conclusive answer as to why glacial ice in the Karakoram is acting differently from the near‐global indication of glacial ice changes. Climate data from high altitudes are needed to find answer for this anomalous behaviour.     

Debris entrainment by basal freeze‐on and thrusting during the 1995–1998 surge of Kuannersuit Glacier on Disko Island,west Greenland

Kuannersuit Glacier, a valley glacier on Disko Island in west Greenland, experienced a major surge from 1995 to 1998 where the glacier advanced 10·5 km and produced a ～65 m thick stacked sequence of debris‐rich basal ice and meteoric glacier ice. The aim of this study is to describe the tectonic evolution of large englacial thrusts and the processes of basal ice formation using a multiproxy approach including structural glaciology, stable isotope composition (δ¹⁸O and δD), sedimentology and ground‐penetrating radar. We argue that the major debris layers that can be traced in the terminal zone represent englacial thrusts that were formed early during the surge. Thrust overthrow was at least 200–300 m and this lead to a 30 m thick repetition of basal ice at the ice margin. It is assumed that the englacial thrusting was initiated at the transition between warm ice from the interior and the cold snout. The basal debris‐rich ice was mainly formed after the thrusting phase. Two sub‐facies of stratified basal ice have been identified; a lower massive ice facies (S_M) composed of frozen diamict enriched with heavy stable isotopes overlain by laminated ice facies (S_L) consisting of millimetre thick lamina of alternating debris‐poor and debris‐rich ice. We interpret the stratified basal ice as a continuum formed mainly by freeze‐on processes and localized regelation. First laminated basal ice is formed and as meltwater is depleted more sediment is entrained and finally the glacier freezes to the base and massive diamict is frozen‐on. The increased ability to entrain sediments may partly be associated with higher basal freezing rates enhanced by loss of frictional heat from cessation of fast flow and conductive cooling through a thin heavily crevassed ice during the final phase of the glacier surge. The dispersed basal ice facies (D) was mainly formed by secondary processes where fine‐grained sediment is mobilized in the vein system of ice. Our results have important implications for understanding the significance of basal ice formation and englacial thrusting beneath fast‐flowing glaciers and it provides new information about the development of landforms during a glacier surge. Copyright © 2010 John Wiley & Sons, Ltd.     

Improved sea ice parcel trajectories in the Arctic via data assimilation

An assimilated sea ice motion product is used to track ice parcels in several regions of the Arctic over time periods of one day to several weeks during 1992-1993. Motions simulated using a two-dimensional, dynamic-thermodynamic sea ice model are combined with motions derived from daily 85 GHz special sensor microwave/imager (SSM/I) imagery using an optimal interpolation method that minimizes error covariance. Assimilation attenuates the tracking error over the stand-alone model in comparison to buoy trajectories with the same starting location and time. The average 14-day assimilated trajectory's displacement error is as much as 34% lower than the model trajectory, while the RMS direction error is decreased by up to 10 degrees (24%). Assimilation can also yield an estimate of dispersion, which is not retrievable by point buoy observations. An assimilation approach improves estimates of ice drift and has the potential to further the understanding of ice mass flux, freshwater flux, and pollutant transport in the polar regions.     

Structure of summer atmospheric boundary layer in the center of Arctic Ocean and its relation with sea ice extent change

The atmospheric vertical structure and changed characteristics of boundary layer parameters, as well as their relations with sea ice and temperature changes in the center of Arctic Ocean(80°–88°N) are presented by adopting GPS sounding data obtained by the 4th–6th Arctic expeditions of China and NCEP(National Centre for Environmental Prediction) reanalysis data. Obvious differences are observed regarding the tropopause, boundary layer height, temperature inversion, and vertical structure of wind speed and direction in the center Arctic Ocean in the summer of 2012, 2010, and 2014. These differences can be explained by the relations between temperature and changes in sea ice extent in September from 1979 to 2014. In September 2012, the Arctic sea ice extent decreased by 44% an with obvious warming process. In September 2010 and 2014, it decreased by 22.6% and 17% with an obvious cooling process, respectively. A comparison of the two processes shows that sea ice change has a significant influence on the structure of the atmospheric boundary layer. In the recent 30 years, the temperature changes of 1000 and 850 h Pa in the center of the Arctic Ocean have displayed an obvious warming trend and negative correlation with sea ice extent. These changes indicate that the continuous reduction of Arctic sea ice will continue the warming of the troposphere middle layer.     

The role of GPR techniques in determining ice cave properties: Peña Castil ice cave,Picos de Europa

The structure and ice content of ice caves are poorly understood. Ground penetrating radar (GPR) can provide useful insights but has only rarely been applied to ice caves. This paper interprets GPR images (radargrams) in terms of internal structure, stratification, compaction, thickness and volume of the ice block in the Peña Castil ice cave (Central Massif of Picos de Europa, northern Spain), providing the endokarst geometry of the ice cave in GPR data reflections. Eight radargrams were obtained by applying a shielded ground‐coupled antenna with a nominal frequency of 400 MHz. Although the radargrams do not depict the ice–basal bedrock interface, they suggest that the ice block is at least 54 m deep and similarly thick. Some curved reflection signatures suggest a potential vertical displacement in the block of ice, and thus certain dynamics in the ice body. Other images show numerous interbedded clasts and thin sediment layers imaged as banded reflections. In this particular cave a direct visual inspection of the ice stratigraphy is a difficult task but GPR provides clear reflectivity patterns of some of its internal features, making GPR a suitable instrument for this and future studies to achieve a better and broader understanding of the internal behavior of ice caves.     

Macro‐scale bed roughness of the siple coast ice streams in West Antarctica

Subglacial bed conditions are key to understanding ice stream behaviour and evolution, with bed roughness re?ecting substrate composition and ?ow resistance. Here we present an analysis of bed roughness in the Siple Coast region of West Antarctica from airborne radio‐echo sounding data. The ice streams are associated generally with low bed roughness values, which decrease downstream. The bed of the slow‐?owing Ice Stream C (～10 000 km²) is also characterized by being smooth at all scales (wavelengths ranging from 5 km to in excess of 40 km). Furthermore, the bed is smooth either side of Ice Stream C. This suggests the location of the ice stream is controlled by internal ice sheet dynamics rather than by bed morphology. If the ice stream were encouraged to migrate laterally, when active, there would be little resistance offered by the subglacial morphology. Other inter‐ice stream regions are rough, however, indicating a subglacial topographic in?uence on ice stream position. Bed roughness increases up‐?ow of ice streams, which, unless the bed is modi?ed, may limit the inland migration of these systems. Copyright © 2004 John Wiley & Sons, Ltd.