Impact of Greenland's topographic height on ... - Wiley Online Library

2 downloads 63 Views 4MB Size Report
Oct 25, 2011 - cending air over the central cold region. Following the above theory, the reduction of the GrIS height would place. Greenland more often in the ...
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117, D09107, doi:10.1029/2011JD017052, 2012

Impact of Greenland’s topographic height on precipitation and snow accumulation in idealized simulations M. Z. Hakuba,1 D. Folini,1 M. Wild,1 and C. Schär1 Received 25 October 2011; revised 19 March 2012; accepted 20 March 2012; published 3 May 2012.

[1] The Greenland Ice Sheet (GrIS) is one of the most dominant orographic obstacles for atmospheric flow in the Northern Hemisphere. Within an idealized framework, we investigate the potential impact of a reduced Greenland topography upon precipitation, snow accumulation, and atmospheric circulation over the GrIS. Using the global atmospheric model ECHAM5-HAM at about 1 spatial resolution (T106) and with present-day climatological mean conditions, we perform four 16-year sensitivity experiments that are identical except for the height of the GrIS topography: one control simulation with the present-day Greenland topography, and three simulations with topographies reduced to 75%, 50%, and 25% of the present-day height. This simple reduction of the GrIS topography (as compared to realistic melt dynamics) leads to an overall increase in annual total precipitation and snow accumulation, composed of significant increases in eastern, northern, and central Greenland and decreases on the western slopes. In principle, this gain in snow accumulation raises the possibility of a negative feedback that would stabilize the height of the GrIS. However, this feedback is likely overcompensated by enhanced ablation (positive feedback), as surface air temperatures strongly increase with reduced topographic height. The analysis of changes in circulation patterns indicates that flatter topographies allow the atmospheric flow to penetrate farther inland, enabling precipitation in regions that are presently desert-like. Prominent circulation features change, in particular the all-season Greenland Anticyclone and wintertime Icelandic Low become weaker with lower topography. Citation: Hakuba, M. Z., D. Folini, M. Wild, and C. Schär (2012), Impact of Greenland’s topographic height on precipitation and snow accumulation in idealized simulations, J. Geophys. Res., 117, D09107, doi:10.1029/2011JD017052.

1. Introduction [2] Greenland is the world’s largest island stretching over 2600 km meridionally, between 59 and 84 N, and 950 km zonally, between 11 and 74 W. More than 80% of the island’s area is covered by the Greenland Ice Sheet (GrIS) with a total volume of 2.65  106 km3 [Ohmura et al., 1996]. The GrIS largely determines the topography and peaks at approximately 3200 meters above sea level. The coastal slopes are very steep, most prominently in the southeast, the wettest region in Greenland, which is directly hit by onshore flow from the northern flank of the Icelandic Low [Ohmura and Reeh, 1991]. The evolution of the GrIS has received considerable attention in climate research, as its mass balance is particularly sensitive to climate change, and a full or partial meltdown of the GrIS would have major repercussions on sea level rise. 1

Institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland. Corresponding Author: M. Z. Hakuba, Institute for Atmospheric and Climate Science, ETH Zurich, Universitätsstrasse 16, CH-8092 Zurich, Switzerland. ([email protected]) Copyright 2012 by the American Geophysical Union. 0148-0227/12/2011JD017052

[3] The massive size of the present-day GrIS is known to impinge substantially on synoptic and planetary scale atmospheric flow [e.g., Broccoli and Manabe, 1992; Kristjánsson and Mcinnes, 1999; Petersen et al., 2004]. Moreover, a number of studies point out the relevance of atmospheric circulation for Greenland precipitation, which is mostly orographically induced over the coastal slopes by passing cyclones [Chen et al., 1997; Schuenemann et al., 2009]. The orographic uplift leads to immediate precipitation over the outer margins of the ice sheet, not affecting regions farther inland due to a lack of precipitable water [Ohmura and Reeh, 1991]. Most of Greenland’s interior is in the precipitation shadow and represents the so-called “dry slot” in annual mean precipitation amounts [Roe, 2005]. [4] Although the flow past Greenland is very complex, with a continuously changing upstream flow, some indication of the atmospheric response may be obtained from idealized studies of uniform flow past isolated topography. Sufficiently high orography, as embodied by the present-day GrIS, may induce upstream blocking and splitting of atmospheric flow, which in some cases is associated with vortex shedding in the lee of the obstacle [Smith, 1989; Schär and Durran, 1997]. Depending on the dimensionless obstacle height and the Rossby number, two main flow regimes may be distinguished [Trüb and Davies, 1995; Ólafsson and

D09107

1 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

Bougeault, 1997; Schär, 2002; Petersen et al., 2003]: either (1) the flow passes over the obstacle, remaining unblocked and quasi-stationary for low dimensionless heights, or (2) the flow is blocked and forced to go around with a potential for vortex shedding to the lee of the obstacle. In this framework, present-day Greenland is mainly placed in the “flow around” regime, but also partially in the “flow over” regime [Schär, 2002]. In the latter case, flow that is directed over Greenland may induce the formation of an anticyclone sitting at the ice sheet top, similar to quasigeostrophic flow solutions. Indeed, the so-called Greenland (glacial) anticyclone [Hobbs, 1945] is a well-known allseason circulation feature over the GrIS, induced by descending air over the central cold region. Following the above theory, the reduction of the GrIS height would place Greenland more often in the “flow over” regime, which in turn would lead to larger vertical and weaker horizontal deflections of the incident flow. As the flow field considerably determines the precipitation pattern over the ice sheet, a feedback upon snow accumulation and surface mass balance is to be expected as well. [5] The fourth IPCC assessment report 2007 [Lemke et al., 2007] summarizes observational estimates of the presentday GrIS mass balance from numerous studies, and finds that most results indicate accelerating mass loss from Greenland during the 1990s up to 2005. Various studies using GCM calculations support these findings and indicate that the mass loss of the ice sheet will become even more substantial for the coming century [Meehl et al., 2007]. Ridley et al. [2005] (HadCM3) find that under a 4  CO2 climate the GrIS undergoes continuous reduction. According to Ridley et al. [2005, Figure 4], the GrIS might shrink to 75% of its volume after approximately 300 years, to 50% after 750 years, and down to 25% after 1500 years. [6] Ablation is highly sensitive to temperature changes near the ice sheet margins, where runoff of meltwater is seen to increase in all models. Most studies indicate that the increase in ablation exceeds the increase in accumulation, leading to enhanced GrIS mass loss and a net positive contribution to global sea level in the 21st century. [7] The present study investigates in an idealized framework how a reduced Greenland topographic height, as may occur in a warmer climate, would interact with the local atmospheric flow and affect the amount and spatial distribution of precipitation and snow accumulation over the GrIS. The investigation of such atmosphere-topography interactions will add to our understanding of the dynamics of an eventual meltdown of the GrIS in a warmer climate. However, in assessing the potential of such feedbacks, we will assume an otherwise unperturbed climate and investigate the role of Greenland’s height under current climatic conditions, i.e., we prescribe global sea-surface temperatures (SST) and greenhouse gas concentrations at current levels. [8] Greenland precipitation is expected to change not only due to the topographic effects as discussed in this paper, but will also increase in response to the anticipated global mean warming. There is a scientific consensus that the global specific humidity will rise in a warmer climate by about 6–7%K 1, consistent with the Clausius-Clapeyron relation [e.g., Boer, 1993; Allen and Ingram, 2002; Held and Soden, 2006; Meehl et al., 2007]. The intensity of the hydrological cycle is, however, also constrained by the availability of

D09107

energy and not by moisture alone [Pall et al., 2007; Wild et al., 2008]. Thus, the increase in global mean precipitation is smaller and estimated to be near 3.4%K 1 (based on model simulations under equilibrium 2  CO2 [Allen and Ingram, 2002]) or 2%K 1 (based on climate change experiments (A1B) of the fourth IPCC assessment report [Held and Soden, 2006]). This additional effect is not accounted for in the present study. It is also evident that in the short term (i.e. the next 100 years) changes in Greenland precipitation are likely driven by global mean warming, while the potential precipitation changes addressed in the current paper would act on a slower timescale but become relevant as soon as there is a significant net mass loss over the GrIS. [9] In Section 2 we give a brief description of the GCM and experimental setup used in this study. A validation of total precipitation (TP) and snow accumulation for the current GrIS is presented in Section 3. Section 4 describes the topographically induced changes in amount and spatial distribution of TP, snow accumulation, and ablation over the flatter Greenland topographies. Changes in circulation patterns, the main reason for the redistribution of Greenland precipitation, are the topic of Section 5. Sections 6 and 7 discuss our results and present our conclusions.

2. Model and Experiments [10] We use the fifth generation general circulation model ECHAM5, which was developed at the Max Planck Institute for Meteorology, Hamburg (MPI) [Roeckner et al., 2003] and evolved from the weather forecasting model of the European Centre for Medium Range Weather Forecasts (ECMWF) [Simmons et al., 1989]. In our version, ECHAM5 incorporates the Hamburg Aerosol Model (HAM) [Stier et al., 2005] and sophisticated cloud microphysics [Lohmann et al., 2007], providing a detailed description of tropospheric aerosols and their interaction with clouds. We run the model with spectral resolution of T106 (about 1.1 grid spacing on a corresponding Gaussian grid) with 31 vertical levels from the surface up to 10 hPa (ca. 30 km). It is clear that a resolution of 1.1 is poor compared to the resolution used in RCMs [Box et al., 2006; Fettweis, 2007; Ettema et al., 2009] but only a GCM can evaluate the impact of a lower GrIS topography on the general circulation, which impacts in its turn the precipitation pattern. [11] We chose to prescribe climatological conditions in order to reduce the natural variability of the system. Aerosol emissions and greenhouse gas (GHG) concentrations are kept constant at their 1870 values. For the sea surface temperatures (SST) and sea ice coverage (SIC) we use a climatological data set based on 1961–1990 Hadley Centre data [Rayner et al., 2003]. [12] Four simulations were carried out, each of them covering 16 years, plus 3 months spin up time. The simulations are identical except for the Greenland topographies. We conducted a control simulation (in the following referred to as normGR) with a present-day Greenland topography, and three simulations with Greenland topographies reduced to 75%, 50%, and 25% of normGR (in the following referred to as 75GR, 50GR, and 25GR). The reduction was applied evenly over the entire Greenland topography by multiplying the input orography (bedrock plus ice sheet) at each spatial

2 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

D09107

only those land grid points that are permanently covered by ice. Based on this approach, the horizontal extent of the GrIS is 2 * 106 km2, covered by 493 grid points in our model resolution (T106). We further subdivided the GrIS into five regions: south (S), west (W), central (C), east (E), and north (N) Greenland. The subdivision follows earlier studies [e.g., Chen et al., 1997; Schuenemann and Cassano, 2010] and is shown in Figure 1, along with the Greenland topography as represented by the model resolution (T106).

3. Validation

Figure 1. Greenland topography (normGR) as represented by the model resolution T106 in meters, and subdivided into five regions: south (S), west (W), east (E), north (N), and central (C) Greenland. grid point in the GrIS with the corresponding factors. This simple reduction contrasts with more detailed modeling studies by, e.g., Ridley et al. [2005], who find that the GrIS melt in a 4  CO2 climate starts in the low-lying coastal regions (i.e., the north and southwest) and is strongly constrained by the eastern mountain ranges. Our approach has been selected as it represents a well controlled and more easily interpretable idealized framework for sensitivity studies about topography-precipitation feedbacks. [13] To define the spatial extent of the GrIS in the model data, we combined the model output variables “land sea mask” and “fraction of land covered by glaciers” to keep

[14] Previous studies already demonstrated the suitability of ECHAM5 for studying Greenland’s climate [Schuenemann et al., 2009; Mernild et al., 2010]. Comparing 15 GCMs, Walsh et al. [2008] found ECHAM5 to give realistic results with respect to output variables such as temperature, precipitation, and sea level pressure for Greenland, Alaska, panArctic, and the Northern Hemisphere, ranking it among the top two GCMs (MPI ECHAM5 and GFDL CM2.1) for such studies. [15] In line with these studies, Greenland precipitation from our normGR simulation is well within the range of published observational estimates, as listed in Table 1, and also compares well with a selection of previous modeling studies. The GrIS-averaged annual mean total precipitation (TP) for normGR is 353 mm yr 1. The published observational estimates (Table 1) range from 340 to 360 mm yr 1 with identical mean and median values of 350 mm yr 1. The modeling study results, as listed in Table 1, range from 336 to 494 mm yr 1, with a mean value of 395 mm yr 1 and median of 388 mm yr 1. [16] The spatial distribution of annual TP over normGR, shown in Figure 2 (left), is found to be in general agreement with the precipitation map by Ohmura and Reeh [1991] (Figure 2, right) based on observations made at 35 meteorological stations in the coastal regions of Greenland, and accumulation measurements of 251 pits and cores obtained from the upper accumulation zone. The resulting pattern also compares well with the distribution map obtained from calibrated regional climate model (RCM) output by Box et al. [2004] (not shown). In accordance with their findings, the

Table 1. Annual Total Precipitation (16-Year Mean) and Mass Balance Parameters Over the NormGR Ice Sheet and Published Estimates From Literaturea

NormGR Obs., in situ Obs., in situ Obs., in situ Obs., microwave Obs., ice cores ECMWF ERA-40 ECMWF ERA-40 ECHAM3 T106 ECHAM4 T106 Polar MM5 MAR RCM HIRHAM4 RCM ECHAM5 T213/T319

Reference

TP

Snowfall

Acc.

Abl.

SMB

This study Benson [1960] Ohmura and Reeh [1991] Reeh et al. [1999] Mote [2003] Bales et al. [2009] Hanna et al. [2005] Hanna et al. [2008] Ohmura et al. [1996] Wild et al. [2003] Box et al. [2006] Fettweis [2007] Mernild et al. [2010] Bengtsson et al. [2011]

353

303

290 340 317 327 293 300 311

75/379/212

215/ 89/78

200 165 151

117 162 142

152 158 146 152 219 165 218 161/205

159 175 280 179 100 167 60 193/153

340 350 360 494 383 377 336 354 408/426

a

337

366 363 354/358

426 331 319 332

Columns denote model or observational reference, total precipitation (TP), snowfall, snow accumulation (Acc.), ablation (Abl.), and surface mass balance (SMB) in mm yr 1 water equivalent (w.e.). For normGR, ablation and SMB are determined in three ways: (1) based on summer mean temperature [Ohmura, 2001], (2) based on energy balance [Braithwaite and Olesen, 1990], and (3) by combining methods 1 and 2. See text for details.

3 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

Figure 2. (left) Annual precipitation over normGR (mm yr 1 w.e.) with contour interval 100 mm yr and (right) observed precipitation by Ohmura and Reeh [1991]. normGR precipitation maxima are located in the southeastern coastal regions, exceeding 1900 mm yr 1, and on the western slopes with up to 600 mm yr 1. The smallest amount of precipitation reaches north and central Greenland, with a minimum around 70 mm yr 1. [17] The normGR seasonal cycles of TP over the GrIS and for the individual regions (see Section 2) compare quite well with corresponding ERA-Interim climatologies (1979–2011), except for the southern part, where the reanalysis shows a distinct seasonal cycle with a minimum during summer and a maximum during winter (not shown). We attribute the less pronounced seasonal cycle in normGR to the comparatively coarse resolution of the coastal margins [cf. Wild et al., 2003]. ERA-Interim has been found to reproduce the spatial pattern and temporal variability of Greenland TP reasonably well [Chen et al., 2011]. [18] The simulated annual snow accumulation (snow fall minus evaporation) over the normGR ice sheet is with 290 mm yr 1 in good agreement, although somewhat on the lower side, compared to the observation based estimates listed in Table 1, which have a median of 314 mm yr 1. The spatial pattern closely follows the TP pattern, with lowest rates found in central Greenland around the summit area, local peaks along the west coast (around 500 mm yr 1), and

D09107

1

maxima in the southeast exceeding 1400 mm yr 1. The spatial distribution of snow accumulation is in general agreement with the observation-based accumulation map by Ohmura et al. [1999] and the aforementioned RCM mean accumulation pattern for the period 1958–2007 by Burgess et al. [2010] (not shown).

4. Impact on Precipitation and Surface Mass Balance 4.1. Changes in Precipitation [19] The GrIS-averaged annual mean total precipitation (16-year mean) increases over a reduced Greenland topography. Reducing the GrIS height from 100% down to 25% of its present-day height, raises the annual mean TP by 36%. The increase is approximately linear with on average +11% per 25% topographic height reduction (Table 2). According to a Kolmogorov-Smirnov-Lilliefors test, the GrIS-averaged annual means are normally distributed at the 5% significance level, as is to be expected for natural variability. To examine whether these changes are statistically significant, we can thus apply the parametric two-sample t test on the GrISaveraged annual means. We find the differences between normGR and the reduced-GrIS simulations 50GR and 25GR

4 of 15

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

D09107

D09107

Table 2. Annual (16-Year Mean) Precipitation, Temperature, and Mass Balance Parameters Over the Entire Ice Sheet From our NormGR and Reduced-Greenland Simulationsa

NormGR 75GR 50GR 25GR

TP

Snowfall

snow TP

353 382 432 479

303 323 344 352

86 85 80 73

Temp2 21 19 17 15

Acc.

Abl. (1), (2), (3)

SMB (1), (2), (3)

290 304 319 328

75, 379, 212 114, 547, 293 206, 856, 534 476, 1091, 1091

215, 89, 78 190, 243, 11 113, 537, 215 148, 763, 763

a Units are mm yr 1 w.e. unless otherwise indicated. Columns denote simulation, total precipitation (TP), snowfall (absolute and % of TP), 2-meter air temperature (Temp2) ( C), snow accumulation (Acc.), ablation (Abl.) and surface mass balance (SMB), estimated in three different ways (see Table 1). Year-to-year standard deviations of TP are in the range 25–42 mm yr 1.

to be significant at the 95% significance level, but insignificant for 75GR. The result is confirmed by the nonparametric Wilcoxon-Mann-Whitney test. All mean values considered in the following passed a Kolmogorov-SmirnovLilliefors test and we always obtained very similar results for both the t test and the Wilcoxon-Mann-Whitney test. Only the results of the t test are given in the following. [20] The spatial pattern of annual TP changes significantly with reduced GrIS height (Figure 3, hatched areas indicate significance at the 95% significance level as determined for each individual grid box). Precipitation on the western slopes weakens in all three reduced-Greenland simulations, exceeding 50% around Thule in the northwest of 25GR. East of it, i.e., in most of north, (south-)eastern, and central Greenland, precipitation is enhanced. In terms of percentage,

most of the TP increase occurs in central Greenland, where the absolute topography decrease is largest and absolute TP values are smallest (not shown). Regions that are presently very dry receive enhanced precipitation, whereas rather wet regions, especially the western and partially the southeastern slopes (25GR), become drier. The resulting TP pattern of 25GR shows a more distinct north-south gradient instead of a desert in the central summit area and less extreme TP rates along the coastal margins than for normGR (Figure 3). [21] The multiyear seasonal cycles of monthly mean TP averaged over the entire GrIS (Figure 4, top row, left) show that TP increases almost uniformly throughout the year, seasonality does not change substantially when reducing the topography. In all simulations, most TP is generated during summer and autumn.

Figure 3. (top) Annual total precipitation (16-year mean) over (left to right) normGR, 75GR, 50GR, and 25GR (mm yr 1 w.e.), with contour interval 100 mm yr 1. (bottom) Difference fields (left to right) 75GRnormGR, 50GR-normGR, and 25GR-normGR. Hatching indicates regions where the difference in the TP annual means is significant at the 95% significance level (based on the two-sample t test). 5 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

D09107

Figure 4. Multiyear monthly mean total precipitation (mm month 1 w.e.): normGR (black line), 75GR (blue), 50GR (green), 25GR (red), averaged over the entire GrIS, as well as over south, west, east, north, and central Greenland. Included are the 95% confidence intervals for normGR (gray shading) and 25GR (reddish shading), based on the one-sample Student’s t test. [22] A look at individual regions (see Section 2) instead of the entire GrIS, reveals substantial differences in normGR TP seasonality (Figure 4, black line). Most precipitation reaches the west during summer, whereas in east, north, and central Greenland the seasonal cycles are overall flatter with weak maxima during summer and autumn. In south Greenland, the normGR TP shows no distinct seasonality. Reducing the GrIS height increases TP in the south throughout the year (75GR, 50GR) with a slight tendency to decrease during winter in 25GR. Precipitation in west Greenland declines during summer but barely changes during the colder months. In east and central Greenland, TP increases significantly (95% significance level of twosample t test applied to monthly and regional means) throughout the year (25GR-normGR), and the summer and autumn peaks become more distinct. In north Greenland, TP increases abruptly during most of the seasons (except summer) for a topographic reduction down to 25GR. 4.2. Changes in Snow Accumulation [23] The snow accumulation over the ice sheet is determined by the amount and spatial distribution of solid precipitation or snowfall, respectively, and to a lesser extent by evaporation. Over normGR, the total precipitation (TP) is mostly solid, exceeding 86% in the annual mean (see Table 2). [24] As the solid precipitation fraction decreases over a flatter and thus warmer Greenland, the increase in accumulation cannot be as high as the increase in TP. Reducing Greenland’s height increases the GrIS-averaged annual snow accumulation approximately linearly by on average

+4% per 25% topographic reduction (Table 2). As with total precipitation, these changes are significant at the 95% level for 50GR and 25GR, but not for 75GR. The overall increase in snow accumulation with decreasing height of topography is a composite of different tendencies in space and season. Snow accumulation decreases most prominently on the western slopes during the warmer months, but increases most distinctly in east and central Greenland during the colder months (Figures 5 and 6). [25] The spatial patterns of annual mean snow accumulation and TP are closely related to each other (see Section 3.1). Nonetheless, the seasonality differs because of the influence of surface air temperature on the solid fraction of precipitation. Unlike for TP, accumulation has strong minima during the summer months and is dominated by the colder months. This pattern gets more pronounced as the GrIS is reduced, especially in eastern and central Greenland (Figure 6). The warmer the summers, the lower the associated accumulation rates. In July, the GrIS-averaged multiyear monthly mean surface air temperature reaches close to 0 C in 25GR (not shown). The implied higher frequency in positive degree days not only results in a shift from solid to liquid precipitation, but also strongly increases ablation (see next Section). 4.3. Ablation Estimates [26] Ablation is strongly affected by surface air temperature. In the context of global warming it has been shown that beyond a critical temperature increase over Greenland (i.e., +4.5 C [Gregory and Huybrechts, 2006]) ablation accelerates, leading to a decline in surface mass balance (SMB)

6 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

Figure 5. Same as Figure 3 but for snow accumulation.

Figure 6. Same as Figure 4 but for snow accumulation. For comparison, TP of normGR (black dashed) and 25GR (red dashed) are included. 7 of 15

D09107

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

[e.g., Meehl et al., 2007]. Recent ECHAM5 T213 calculations by Bengtsson et al. [2011] project 200% more ablation for a 5.2 C warming over the GrIS (IPCC scenario A1B, +3 C globally), turning the initially positive SMB negative. In our simulations, the surface warming is due to the reduction of the topographic height. With a topographically induced warming of around +2 C per 25% topographic reduction, up to 6 C over 25GR (Table 2), we expect a comparable increase in ablation, as the melt zone will expand spatially into Greenland’s central regions, and temporarily into presently too cold seasons. [27] Even though the model resolution of T106 is higher than typically used in climate model simulations, it is still fairly coarse for an adequate determination of ablation [Glover, 1999; Wild et al., 2003]. To derive a rough estimate, we determined ablation in three ways: [28] (1) Following Wild et al. [2003], we estimated ablation from the mean summer temperature, based on an empirical formulation after Ohmura [2001], which we applied to each grid point over the GrIS. The computation is based on the observed relationship between the annual ablation and the collocated mean air temperature for the summer months (JJA) [Ohmura et al., 1996]. For normGR, this methodology yields an annual ablation of 75 mm yr 1, which is significantly below published estimates (Table 1), and may be considered as a conservative estimate. For the reduced Greenland topographies, this approach is not fully justified, as the empirical relationship stems from the current-day GrIS. Applying it nevertheless, suggests accelerated mass loss with reduced topography, the enhanced accumulation being overcompensated by the warming of up to +6 C (25GR) and the associated ablation increase (Table 2). [29] (2) Based on a simple energy balance model [Braithwaite and Olesen, 1990], we calculated the residual energy available for melting over the GrIS as the sum of the net radiation, sensible, and latent heat fluxes. This approach does not capture the refreezing of meltwater in the firn [Pfeffer et al., 1991; Ohmura et al., 1996] and overestimates ablation strongly, resulting in a negative SMB for normGR (Table 1). [30] (3) To overcome the shortcoming of method (2), we combined it with method (1) in calculating the melt based on the energy balance only in grid boxes, in which the temperature criterion of method (1) indicates that ablation takes place and, thus, refreezing does not occur. The GrISaveraged annual ablation for normGR lies then with 212 mm yr 1 well within the range of observational and model estimates (Table 1). [31] All three methods, regardless whether our normGR ablation agrees with published estimates for present-day conditions or not (Table 1), indicate that the ablation increase over a flatter GrIS sooner or later dominates over the gain in snow accumulation (Table 2). The SMB declines (all methods) and eventually turns negative in 50GR (method 3) or 25GR (method 1), respectively.

5. Circulation Changes and TP Patterns [32] Various earlier studies find a clear relationship between the circulation and precipitation patterns over Greenland, arguing that TP is mostly due to orographic effects acting upon

D09107

the mean flow and passing extratropical cyclones [e.g., Ohmura and Reeh, 1991; Chen et al., 1997; Schuenemann et al., 2009; Schuenemann and Cassano, 2010]. Here, we use two methodologies to assess this relationship and the associated changes as the height of Greenland is reduced. First, one common approach is to inspect the mean sea level pressure (SLP) patterns and wind fields in summer (JJA) and winter (DJF). This is well motivated as cyclonic activity highly depends on the season, and as snow accumulation and also TP show a seasonal cycle in at least some of the regions of Greenland (especially the west, Figures 4 and 6). Second, we will identify high precipitation months for each of the five regions introduced in Section 2, and determine the corresponding circulation patterns by composite analysis, following Chen et al. [1997]. This approach is particularly appropriate for regions where precipitation shows no distinct seasonality. 5.1. Winter and Summer Flow Patterns [33] The analysis is conducted using maps of precipitation, sea level pressure (SLP), 500 hPa geopotential height (gpm500), and low-level wind (Figures 7 and 8). In these plots we exclude SLP in the area of Greenland, as pressure reduction leads to substantial errors. [34] The spatial distributions of winter and summer precipitation are determined by characteristic flow patterns and differ substantially: In winter, TP maximizes on the southeastern slopes, while the other regions experience their seasonal minima (Figure 7, first row, left). Reducing Greenland’s height significantly increases TP over the higher plateau in north, east, and central Greenland (Figure 7, first row). The TP maxima in the southeast initially amplify (75GR and 50GR), but then decrease again (25GR). In summer, TP is more evenly distributed over Greenland than in winter, maximizing in the very south and on the slopes of western and southeastern normGR (Figure 8, top row, left). A flatter GrIS leads to a decrease in TP on the western slopes and an increase in all other regions, predominantly the summit area (Figure 8, top row). 5.1.1. Winter Circulation [35] In normGR, the winter circulation in terms of SLP and wind field (at the sixth terrain-following model level, approximately 860 hPa over sea) depicts a long-stretched, banana-shaped Icelandic Low with its center just east of the southern tip that produces strong onshore winds striking the southeastern coast (Figure 7, second and third rows, left). Once these southeasterlies reach the higher plateau, they are deflected north and transition into anticyclonic motion around the summit area. Commonly referred to as the Greenland Anticyclone, this is an expected flow feature over Greenland’s glacial surface [Hobbs, 1945], and in agreement with the quasi-geostrophic “flow over” theory [Schär, 2002]. [36] The reduction of the GrIS height leads to a weaker Icelandic Low, SLP increases over the Greenland and Norwegian Seas, and decreases to the southwest (50GR and 25GR) (Figure 7, second row, left). The weakening of the Icelandic Low is accompanied by the diversion of the North Atlantic storm track. The storm track difference pattern in terms of 2.5–8 days band-pass filtered variability of the 500 hPa geopotential height (gpm500), indicates enhanced storm activity above and to the west of a flatter Greenland, and weaker activity over the eastern North Atlantic (not

8 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

D09107

Figure 7. (top) Winter (DJF) (left to right) precipitation normGR and differences 75GR-normGR, 50GRnormGR, and 25GR-normGR (mm yr 1 w.e.). The hatching indicates regions where the difference is significant at the 95% significance level. (middle) Same as top plots but for sea level pressure (hPa). (bottom) DJF wind fields at the sixth model level (approx. 860 hPa over sea) plus gpm500 (red contours), over (left to right) normGR, 75GR, 50GR, and 25GR. The vector length scales with wind velocity. shown). Thus, in the case of a flatter Greenland topography, storms move more frequently north along the west coast instead of the east coast, which is their natural track under current conditions, and traverse the island without being forced to go around it. These findings are in line with Junge et al. [2005], who suggest that more cyclones pass over a flat-bottom Greenland topography. [37] Previous studies have shown that westerly flow crossing the southern tip of Greenland induces leecyclogenesis that significantly strengthens the Icelandic low [Chen et al., 1997; Kristjánsson and Mcinnes, 1999; Doyle and Shapiro, 1999], and that the “lee cyclone” declines or even vanishes in simulations without a Greenland topography [Kristjánsson and Mcinnes, 1999; Petersen et al., 2003; Junge et al., 2005; Tsukernik et al., 2007]. Indeed, the normGR wind field indicates the presence of a winter mean “lee cyclone” or “tip vortex” to the east of the southern tip, which becomes weaker as the height of the GrIS is reduced and essentially disappears in 25GR (Figure 7, second and third rows). In 75GR and 50GR the “lee cyclone” still exists and the TP maxima in the southeast amplify, as onshore winds are less deflected and hit the coast more perpendicular (Figure 7, first and second rows). In 50GR, the TP maxima are shifted inland, which is in accordance with theory, implying that orographic precipitation moves closer

to the mountain crest when the topographic height is reduced [Roe, 2005]. In 25GR, the missing “lee cyclone” and weaker Icelandic Low result in weaker easterly onshore winds and less precipitation at the southeastern margin. [38] The wind fields are furthermore consistent with a weakening of the Greenland Anticyclone, most prominently in 25GR (Figure 7, third row). The anticyclone weakens over a flatter (quasi-geostrophic “flow over” theory) and warmer (less cooling and descent of low-level air) Greenland. Its absence allows easterly onshore flow to penetrate farther inland (Figure 7, second row), which appears to be the main reason for the increase in TP over eastern and northern Greenland. The effect on west Greenland is reversed. The southerly flow to the west of Greenland (which brings warm moist air to the northern portion of the west coast) is diminished, in response to the weakening of the Greenland Anticyclone. As a result, the flow is more directed over (rather than around) the Greenland topography, yielding increased upstream and reduced downstream precipitation. At the same time, as the orography is reduced, it becomes easier for the southwesterly and -easterly flows to reach the central plateau with the potential for higher TP rates. Looking at the gpm500 contours, the geopotential height increases over and northwest of Greenland, indicating the weakening and spatial

9 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

D09107

Figure 8. Same as Figure 7 but for summer (JJA). Note the different color scale for normGR SLP. expansion of the winter polar vortex for a flatter GrIS (Figure 7, third row). 5.1.2. Summer Circulation [39] In summer, the normGR circulation is dominated by low pressure over Northern Canada and Baffin Bay, producing westerly onshore flow. The westerlies sharply deviate north over the western slopes and accelerate around the summit area, again inducing the formation of a Greenland Anticyclone (Figure 8, second and third rows, left). A reduced Greenland topography shifts the mean cyclone southeast into Baffin Bay and over west Greenland in 25GR (Figure 8, second row). The intensification of this Baffin Bay Low [Ohmura and Reeh, 1991] changes the flow pattern quite dramatically, as it produces offshore flow over northwestern Greenland. As in winter, the storm activity increases strongly over Greenland’s interior (not shown). The Greenland Anticyclone ultimately vanishes in 25GR (Figure 8, third row). 5.2. Flow Patterns for High TP Months [40] For all Greenland regions, composite maps are constructed that relate to rainy (or snowy) periods in the respective regions. To this end, for each simulation and Greenland region, we select the 25% of all months with the highest TP during the 16-year simulation period (on average 3 months per year). We then compute their composite mean SLP patterns, 500 hPa geopotential heights and 860 hPa wind fields. The seasonal distribution of the high TP months per region in normGR is as follows: central, north, and west Greenland show distinct summer maxima (with more than

50% of high TP months) and minima in winter (down to 4% in the west). In east and south Greenland, high TP months are more evenly distributed between summer, autumn, and winter (25–31%), the minimum with 17% lies in spring. The reduction in topographic height does not substantially change the seasonal distributions of high TP months. In general, the 25GR simulation exhibits a broadened seasonal peak of the high TP maxima in west, central, and east Greenland, spanning summer (on average 44%) and autumn (42%). [41] The circulation patterns responsible for high TP in south, west, and east Greenland differ substantially (Figure 9), whereas north and central Greenlands’ patterns are combinations of the eastern and western high TP flow regimes. Thus, we concentrate on the former three regions. [42] Over southern Greenland for normGR, high TP is related to a mean cyclone that engulfs the southern tip with the low pressure center just east of it (Figure 9, first row, left). This cyclone produces onshore flow approaching the southeastern slopes (Figure 9, second row, left). Chen et al. [1997] found that if a cyclone exists in the mean monthly sea level pressure in the Labrador sea (Labrador Low), heavy precipitation occurs in southern Greenland. The pressure pattern we found for high TP in the south appears to be a combination of such a Labrador Low and an additional “lee cyclone”, induced by westerly flow crossing the southern tip and causing lee cyclogenesis [Chen et al., 1997; Kristjánsson and Mcinnes, 1999; Doyle and Shapiro, 1999]. Reducing Greenland’s topography changes the pattern noticeably. The mean cyclone’s center weakens and relocates westward

10 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

Figure 9 11 of 15

D09107

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

around the southern tip. In 25GR, the cyclone looks just like the Labrador Low proposed by Chen et al. [1997] (Figure 9, first row, right). The absence of lee cyclogenesis may play a role and the formerly intense onshore winds (Figure 9, second row) and precipitation (Figure 3) decline. [43] Over western Greenland, high TP is associated with a summer-like low pressure structure centered over Baffin Island (Figure 9, third row, left). The westerly onshore winds accelerate around the summit area (Figure 9, fourth row, left). Conducting the analysis with reduced Greenland topography yields a pattern with an intensified low pressure system over the Baffin Bay (Figure 9, third row). A more uniform southwesterly wind field dominates over 25GR (Figure 9, fourth row), reduces TP on the western slopes and provides enhanced TP to the central region. These changes do not necessarily imply that the Baffin Bay low becomes stronger, but that such lows become more important in yielding precipitation over western Greenland. [44] Over eastern Greenland, high TP is produced in late summer by an Icelandic Low that, compared to the winter mean, stretches broader into the Norwegian sea with its center farther off the southeastern coast of Greenland (Figure 9, fifth row, left). The easterly onshore winds sharply deviate south, as they are blocked by steep orography and forthcoming westerly downslope winds, and run parallel to the shore (Figure 9, sixth row, left). The Icelandic Low in these circulation patterns becomes more compact in 50GR and 25GR (Figure 9, fifth and sixth rows), although the wintertime Icelandic Low itself becomes generally weaker with the reduction of the Greenland topography (see Figures 7 and 8). The result thus suggests that TP in eastern Greenland would become more sensitive to the location of the Icelandic Low in a scenario with reduced topographic height. 5.3. Flow and TP Pattern: Synthesis [45] Both analyses in terms of regional high TP and summer/winter circulation patterns yield the following major flow changes as the GrIS height is reduced: (1) The decay of the Greenland Anticyclone reduces the lateral diversion of the incoming flow and thereby increases the lifting and precipitation with the easterly onshore flow (TP: east ↑, north ↑, west ↓), (2) Reduced orographic flow blocking is associated with the shift of orographic TP maxima closer to the mountain crest and increased cross-orographic circulation (slopes ↓, interior ↑), (3) The diversion of the NA storm track contributes to a less intense Icelandic Low (southeast ↓) and to a more intense Baffin Bay Low that wells over a flatter Greenland, inducing offshore flow over the northwest and enhanced inflow into the central regions (interior ↑, west ↓), and (4) Lee-cyclogenesis occurs less frequently (south ↓).

6. Discussion [46] In this study, we use simplified topographic reductions as compared to realistic melt dynamics of the GrIS. Simulations using dynamic ice sheet models [e.g., Ridley

D09107

et al., 2005, 2010] show that the GrIS meltdown evolves non-linearly and exposes bedrock asymmetrically, starting in the north and southwest of the GrIS. Along the low-lying ice sheet margins, the initial warming is strongest and the retreat of ice is accompanied by enhanced surface air temperatures due to elevation- and albedo-temperature feedbacks [e.g., Toniazzo et al., 2004; Ridley et al., 2005]. The eastern mountains act as a barrier, restricting the ice flow to lower elevations, and a dome of ice remains in central-east Greenland long after most of the ice sheet has retreated. Thus, the maximum elevation height of the GrIS is barely affected for many hundred years, even for a climate warming up to +10 C [Greve, 2000]. Bedrock elevation maps by Letréguilly et al. [1991] have mean altitudes of 440 meters before and 812 meters after isostatic rebound and are characterized by two main mountain ranges in eastern Greenland and on the southern tip, reaching up to 1000 and 600 meters before, and up to 1400 and 800 meters after isostatic rebound [Letréguilly et al., 1991; Greve, 2000; Bamber et al., 2001; Toniazzo et al., 2004; Ridley et al., 2005]. The mean height of all our reduced-Greenland topographies (461 meters for 25GR) lies above the present-day bedrock mean altitude. However, our Greenland orographies are generally too smooth, i.e., too flat around the summit area, and too high in low-lying coastal areas (maximum height 830 meters for 25GR). Overall, we expect our results to remain basically valid, once the bedrock is fully exposed. Differences may be expected during the different phases of the complex melt [Ridley et al., 2005] or because of other reasons, e.g., globally changing atmospheric conditions in a 4  CO2 world. [47] Considering the impact of changes in Greenland topography upon the atmospheric circulation, both the height and width of the GrIS determine the flow regime [Trüb and Davies, 1995; Ólafsson and Bougeault, 1997; Schär, 2002; Petersen et al., 2003]. Changes in both these parameters associated with the melt of the GrIS can potentially affect the proportion of flows in the “flow over” and “flow around” regimes. If the melt of Greenland proceeds such as to primarily reduce the height of the GrIS (rather than the shape), then the flow will become more likely to flow over the GrIS and thereby affect precipitation as simulated in the presented numerical experiments. If on the other hand the melt proceeds such as to maintain the peak height of Greenland but to reduce its horizontal scale (in particular in the north-south direction), then the flow will change such as to be more often in the “flow around” regime, thereby potentially reducing precipitation. This effect occurs as a reduction of horizontal extent will reduce the impact of the Earth’s rotation on the selection of the flow regime (i.e. an increase in effective Rossby number). [48] We further introduced uncertainties through the computation of ablation, as we neglect several feedbacks that might either enhance or reduce ablation. On the one hand, we observe a shift from solid to liquid precipitation (see Table 2 and Figure 6) that might enhance ablation during summer through the direct contribution of rain to the

Figure 9. SLP patterns and wind fields at the sixth model level (approx. 860 hPa over sea level, vector length scales with wind velocity) with geopotential height at the 500 hPa pressure level (red contours) for high TP in south (S), west (W), and east (E) Greenland: (left to right) normGR, 50GR, and 25GR. 12 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

melt, which is estimated to be extremely small for presentday conditions [Box et al., 2004], and through the decrease in surface albedo due to less frequent snowfall [Ettema et al., 2009; Screen and Simmonds, 2011]. On the other hand, we have not taken into account the exposure of bare ice during summer that might be delayed due to increased snowfall during the colder months [Tedesco et al., 2011]. Indeed, the GrIS-averaged snow accumulation (Figure 6, top row, left) is enhanced from October to April over a reduced GrIS topography, but there are regional differences. In east, central, and north Greenland a thicker snowpack might build up, whereas in west and partially south Greenland, snow accumulation even decreases (see Figures 5 and 6). In the latter two regions, earlier bare ice appearance might enhance the mass loss of the ice sheet.

7. Conclusions [49] Within this study, we investigate how a reduced height of the Greenland topography would affect the flow and precipitation patterns over the Greenland Ice Sheet (GrIS), and whether there is a significant increase in snow accumulation tied to this, which might act as a negative feedback and slow down the ongoing GrIS mass loss. We use an idealized framework to quantify the feedback at question, rather than a scenario approach that would need to consider the transient melt of the GrIS and specifications of climate forcings. Our approach thus makes several simplifications as we use (1) simplified topographies, (2) prescribe present-day climatological mean conditions, and (3) estimate ablation based on rather vague methods. Nevertheless, the idealized framework has proven to be suitable in isolating the topography-precipitation feedback. [50] We conducted four climatological 16-year ECHAM5HAM T106 simulations that are identical except for the height of the GrIS: 100%, 75%, 50%, and 25% of the present-day height. In the control simulation normGR, total precipitation (TP) and snow accumulation are in reasonable agreement with observations, reanalyses, and previous modeling studies. High precipitation over Greenland occurs mainly in the coastal regions and maximizes in the southeast, predominantly during autumn and winter, and less pronounced on the western slopes, distinctly peaking during summer. [51] Reducing the GrIS height increases the total precipitation (+11% per 25% topographic reduction) and snow accumulation (+4% per 25% topographic reduction). The change in total amounts is accompanied by a regional redistribution (wetter central, east, north, and drier west). The seasonal distribution of TP remains largely unchanged, but snow accumulation becomes more concentrated in winter. The annual mean surface air temperature increases by +2 C per 25% topographic reduction. From the current analysis, there is strong evidence that the associated increase in ablation dominates over the small gain in snow accumulation. A flatter GrIS topography, as might occur in a warmer climate, thus may accelerate ablation, which represents a positive feedback on the GrIS mass loss. [52] We conclude that there must be a relationship between the changes in precipitation pattern – regional increases and decreases – and the partial transition from the “flow around” to the “flow over” regime as the Greenland

D09107

topography is reduced. The topographically induced changes in circulation are thus essential for the increase and redistribution in TP. [53] The present study demonstrates that feedback processes between Greenland’s topography, circulation, and precipitation/snow accumulation may play a relevant role for the evolution of the surface mass balance of the GrIS, potentially both in the past and in the future. [54] Based on the presented results and previous publications [e.g., Broccoli and Manabe, 1992; Kristjánsson and Mcinnes, 1999; Dethloff et al., 2004; Petersen et al., 2004; Junge et al., 2005; Tsukernik et al., 2007], it is to be expected that other prominent flow features in Greenland’s vicinity (e.g., Greenland tip jet [Doyle and Shapiro, 1999]) and in the far field are affected by a shrinking Greenland topography. Thus, it may be interesting to extend the spatial scope of this study into the Arctic region and the entire Northern Hemisphere. Indeed, preliminary analysis of our data at this spatial scale suggests a weakening and spatial expansion of the polar vortex (strongest in 25GR), accompanied by enhanced surface temperatures over the polar region and negative temperature anomalies over mid-latitudinal regions in North America and Eurasia. These findings resemble the features of the negative Arctic oscillation (AO) index [Wallace and Gutzler, 1981]. A flatter Greenland topography seems to shift the AO permanently toward its negative phase. Further research is needed to scrutinize the interaction between the melting GrIS in a changing climate and the Northern Hemispheric large-scale circulation. [55] Acknowledgments. Computations were performed at the Swiss National Supercomputing Center (CSCS), and the authors are most grateful for the continuous help and assistance of CSCS staff. We would also like to thank Atsumu Ohmura and Heinz Blatter for support and extensive discussions, and Ulrike Lohmann and Sylvaine Ferrachat (ETH) as well as MPI Hamburg and the Center for Climate System Modeling (C2SM) for their support in maintaining this model. ECMWF ERA-Interim data used in this study have been obtained from the ECMWF data server. The work was financially supported by the National Center for Competence in Climate Research (NCCR Climate).

References Allen, M. R., and W. J. Ingram (2002), Constraints on future changes in climate and the hydrologic cycle, Nature, 419, 224–232, doi:10.1038/ nature01092. Bales, R. C., Q. Guo, D. Shen, J. R. McConnell, G. Du, J. F. Burkhart, V. B. Spikes, E. Hanna, and J. Cappelen (2009), Annual accumulation for Greenland updated using ice core data developed during 2000–2006 and analysis of daily coastal meteorological data, J. Geophys. Res., 114, D06116, doi:10.1029/2008JD011208. Bamber, J. L., R. L. Layberry, and S. P. Gogineni (2001), A new ice thickness and bed data set for the Greenland Ice Sheet: 1. Measurement, data reduction, and errors, J. Geophys. Res., 106(D24), 33,773–33,780, doi:10.1029/2001JD900054. Bengtsson, L., S. Koumoutsaris, and K. Hodges (2011), Large-scale surface mass balance of ice sheets from a comprehensive atmospheric model, Surv. Geophys., 32, 459–474, doi:10.1007/s10712-011-9120-8. Benson, C. S. (1960), Stratigraphic studies in the snow and firn of the Greenland Ice Sheet, PhD dissertation, Calif. Inst. of Technol., Pasadena. Boer, G. J. (1993), Climate change and the regulation of the surface moisture and energy budgets, Clim. Dyn., 8(5), 225–239, doi:10.1007/ BF00198617. Box, J. E., D. H. Bromwich, and L.-S. Bai (2004), Greenland Ice Sheet surface mass balance 1991–2000: Application of Polar MM5 mesoscale model and in situ data, J. Geophys. Res., 109, D16105, doi:10.1029/ 2003JD004451. Box, J. E., D. H. Bromwich, B. A. Veenhuis, L.-S. Bai, J. C. Stroeve, J. C. Rogers, K. Steffen, T. Haran, and S.-H. Wang (2006), Greenland Ice Sheet surface mass balance variability (1988–2004) from calibrated

13 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

Polar MM5 output, J. Clim., 19(12), 2783–2800, doi:10.1175/ JCLI3738.1. Braithwaite, R. J., and O. B. Olesen (1990), A simple energy-balance model to calculate ice ablation at the margin of the Greenland Ice Sheet, J. Glaciol., 36(123), 222–228. Broccoli, A. J., and S. Manabe (1992), The effects of orography on midlatitude Northern Hemisphere dry climates, J. Clim., 5, 1181–1201, doi:10.1175/1520-0442(1992)0052.0.CO;2. Burgess, E. W., R. R. Forster, J. E. Box, E. Mosley-Thompson, D. H. Bromwich, R. C. Bales, and L. C. Smith (2010), A spatially calibrated model of annual accumulation rate on the Greenland Ice Sheet (1958–2007), J. Geophys. Res., 115, F02004, doi:10.1029/2009JF001293. Chen, L., O. M. Johannessen, H. Wang, and A. Ohmura (2011), Accumulation over the Greenland Ice Sheet as represented in reanalysis data, Adv. Atmos. Sci., 28, 1030–1038, doi:10.1007/s00376-010-0150-9. Chen, Q., D. H. Bromwich, and L. Bai (1997), Precipitation over Greenland retrieved by a dynamic method and its relation to cyclonic activity, J. Clim., 10, 839–870, doi:10.1175/1520-0442(1997)010 2.0.CO;2. Dethloff, K., W. Dorn, A. Rinke, K. Fraedrich, M. Junge, E. Roeckner, V. Gayler, U. Cubasch, and J. H. Christensen (2004), The impact of Greenland’s deglaciation on the Arctic circulation, Geophys. Res. Lett., 31, L19201, doi:10.1029/2004GL020714. Doyle, J. D., and M. A. Shapiro (1999), Flow response to large-scale topography: The Greenland tip jet, Tellus, Ser. A, 51(5), 728–748, doi:10.1034/ j.1600-0870.1996.00014.x. Ettema, J., M. R. van den Broeke, E. van Meijgaard, W. J. van de Berg, J. L. Bamber, J. E. Box, and R. C. Bales (2009), Higher surface mass balance of the Greenland Ice Sheet revealed by high-resolution climate modeling, Geophys. Res. Lett., 36, L12501, doi:10.1029/2009GL038110. Fettweis, X. (2007), Reconstruction of the 1979–2006 Greenland Ice Sheet surface mass balance using the regional climate model MAR, Cryosphere, 1, 21–40, doi:10.5194/tc-1-21-2007. Glover, R. W. (1999), Influence of spatial resolution and treatment of orography on GCM estimates of the surface mass balance of the Greenland Ice Sheet, J. Clim., 12, 551–563, doi:10.1175/15200442(1999)0122.0.CO;2. Gregory, J. M., and P. Huybrechts (2006), Ice-sheet contributions to future sea-level change, Philos. Trans. R. Soc. A, 364(1844), 1709–1731, doi:10.1098/rsta.2006.1796. Greve, R. (2000), On the response of the Greenland Ice Sheet to greenhouse climate change, Clim. Change, 46, 289–303, doi:10.1023/A:1005647226590. Hanna, E., P. Huybrechts, I. Janssens, J. Cappelen, K. Steffen, and A. Stephens (2005), Runoff and mass balance of the Greenland Ice Sheet: 1958–2003, J. Geophys. Res., 110, D13108, doi:10.1029/2004JD005641. Hanna, E., P. Huybrechts, K. Steffen, J. Cappelen, R. Huff, C. Shuman, T. Irvine-Fynn, S. Wise, and M. Griffiths (2008), Increased runoff from melt from the Greenland Ice Sheet: A response to global warming, J. Clim., 21, 331–341, doi:10.1175/2007JCLI1964.1. Held, I. M., and B. J. Soden (2006), Robust responses of the hydrological cycle to global warming, J. Clim., 19, 5686–5699. Hobbs, W. W. (1945), The Greenland glacial anticyclone, J. Meteorol., 2(3), 143–153, doi:10.1175/1520-0469(1945)0022.0.CO;2. Junge, M. M., R. Blender, V. Fraedrich, V. Gayler, U. Luksch, and F. Lunkeit (2005), A world without Greenland: Impacts on the Northern Hemisphere winter circulation in low- and high-resolution models, Clim. Dyn., 24, 297–307, doi:10.1007/s00382-004-0501-2. Kristjánsson, J. E., and H. Mcinnes (1999), The impact of Greenland on cyclone evolution in the North Atlantic, Q. J. R. Meteorol. Soc., 125, 2819–2834, doi:10.1002/qj.49712556003. Lemke, P., et al. (2007), Observations: Changes in snow, ice and frozen ground, in Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, pp. 339–383, Cambridge Univ. Press, New York. Letréguilly, A., P. Huybrechts, and N. Reeh (1991), Steady-state characteristics of the Greenland Ice Sheet under different climates, J. Glaciol., 37(125), 149–157. Lohmann, U., P. Stier, C. Hoose, S. Ferrachat, S. Kloster, E. Roeckner, and J. Zhang (2007), Cloud microphysics and aerosol indirect effects in the global climate model ECHAM5-HAM, Atmos. Chem. Phys., 7(13), 3425–3446, doi:10.5194/acp-7-3425-2007. Meehl, G. A., et al. (2007), Global climate projections, in Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, pp. 749–845, Cambridge Univ. Press, New York. Mernild, S. H., G. E. Liston, C. A. Hiemstra, and J. H. Christensen (2010), Greenland Ice Sheet surface mass-balance modeling in a 131-yr perspective, 1950–2080, J. Hydrometeorol., 11(1), 3–25, doi:10.1175/2009JHM1140.1.

D09107

Mote, T. L. (2003), Estimation of runoff rates, mass balance, and elevation changes on the Greenland Ice Sheet from passive microwave observations, J. Geophys. Res., 108(D2), 4056, doi:10.1029/2001JD002032. Ohmura, A. (2001), Physical basis for the temperature-based meltindex method, J. Appl. Meteorol., 40, 753–761, doi:10.1175/1520-0450 (2001)0402.0.CO;2. Ohmura, A., and N. Reeh (1991), New precipitation and accumulation maps for Greenland, J. Glaciol., 37(125), 140–148. Ohmura, A., M. Wild, and L. Bengtsson (1996), A possible change in mass balance of Greenland and Antarctic Ice Sheets in the coming century, J. Clim., 9, 2124–2135, doi:10.1175/1520-0442(1996)009 2.0.CO;2. Ohmura, A., P. Calanca, M. Wild, and M. Anklin (1999), Precipitation, accumulation and mass balance of the Greenland Ice Sheet, Z. Gletscherkd. Glazialgeol., 35, 1–20. Ólafsson, H., and P. Bougeault (1997), The effect of rotation and surface friction on orographic drag, J. Atmos. Sci., 54, 193–210, doi:10.1175/ 1520-0469(1997)0542.0.CO;2. Pall, P., M. R. Allen, and D. A. Stone (2007), Testing the ClausiusClapeyron constraint on changes in extreme precipitation under CO2 warming, Clim. Dyn., 28(4), 351–363, doi:10.1007/s00382-006-0180-2. Petersen, G. N., H. Ólafsson, and J. E. Kristjánsson (2003), Flow in the lee of idealized mountains and Greenland, J. Atmos. Sci., 60(17), 2183–2195, doi:10.1175/1520-0469(2003)0602.0.CO;2. Petersen, G. N., J. E. Kristjánsson, and H. Ólafsson (2004), Numerical simulations of Greenland’s impact on the Northern Hemisphere winter circulation, Tellus, Ser. A, 56(2), 102–111, doi:10.1111/j.1600-0870. 2004.00047.x. Pfeffer, W. T., M. F. Meier, and T. H. Illangasekare (1991), Retention of Greenland runoff by refreezing: Implications for projected future sea level change, J. Geophys. Res., 96(C12), 22,117–22,124, doi:10.1029/ 91JC02502. Rayner, N. A., D. E. Parker, E. B. Horton, C. K. Folland, L. V. Alexander, D. P. Rowell, E. C. Kent, and A. Kaplan (2003), Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res., 108(D14), 4407, doi:10.1029/ 2002JD002670. Reeh, N., C. Mayer, H. Miller, H. H. Thomsen, and A. Weidick (1999), Present and past climate control on fjord glaciations in Greenland: Implications for IRD-deposition in the sea, Geophys. Res. Lett., 26, 1039–1042, doi:10.1029/1999GL900065. Ridley, J. K., P. Huybrechts, J. M. Gregory, and J. A. Lowe (2005), Elimination of the Greenland Ice Sheet in a high CO2 climate, J. Clim., 18, 3409–3427, doi:10.1175/JCLI3482.1. Ridley, J. K., J. M. Gregory, P. Huybrechts, and J. Lowe (2010), Thresholds for irreversible decline of the Greenland Ice Sheet, Clim. Dyn., 35, 1065–1073, doi:10.1007/s00382-009-0646-0. Roe, G. H. (2005), Orographic precipitation, Annu. Rev. Earth Planet. Sci., 33, 645–671, doi:10.1146/annurev.earth.33.092203.122541. Roeckner, E., et al. (2003), The atmospheric general circulation model ECHAM5, Part I. Model description, MPI Rep. 349, Max Planck Inst. for Meteorol., Hamburg, Germany. Schär, C. (2002), Mesoscale mountains and the larger-scale atmospheric dynamics: A review, in Meteorology at the Millennium: 150th Anniversary Conference of the Royal Meteorological Society, Int. Geophys. Ser., vol. 83, pp. 29–41, Academic, London. Schär, C., and D. R. Durran (1997), Vortex formation and vortex shedding in continously stratified flows past isolated topography, J. Atmos. Sci., 54, 534–554. Schuenemann, K. C., and J. J. Cassano (2010), Changes in synoptic weather patterns and Greenland precipitation in the 20th and 21st centuries: 2. Analysis of 21st century atmospheric changes using self-organizing maps, J. Geophys. Res., 115, D05108, doi:10.1029/2009JD011706. Schuenemann, K. C., J. J. Cassano, and J. Finnis (2009), Synoptic forcing of precipitation over Greenland: Climatology for 1961–99, J. Hydrometeorol., 10, 60–78, doi:10.1175/2008JHM1014.1. Screen, J. A., and I. Simmonds (2011), Declining summer snowfall in the Arctic: Causes, impacts and feedbacks, Clim. Dyn., doi:10.1007/ s00382-011-1105-2, in press. Simmons, A. J., D. M. Burridge, M. Jarraud, C. Girard, and W. Wergen (1989), The ECMWF medium-range prediction models development of the numerical formulations and the impact of increased resolution, Meteorol. Atmos. Phys., 40, 28–60, doi:10.1007/BF01027467. Smith, R. B. (1989), Hydrostatic airflow over mountains, Adv. Geophys., 31, 1–41, doi:10.1016/S0065-2687(08)60052-7. Stier, P., et al. (2005), The aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 5, 1125–1156, doi:10.5194/acp-5-1125-2005. Tedesco, M., X. Fettweis, M. R. van den Broeke, R. S. W. van de Wal, C. J. P. P. Smeets, W. J. van de Berg, M. C. Serreze, and J. E. Box (2011),

14 of 15

D09107

HAKUBA ET AL.: IMPACT OF GREENLAND’S TOPOGRAPHIC HEIGHT

The role of albedo and accumulation in the 2010 melting record in Greenland, Environ. Res. Lett., 6(1), 014005, doi:10.1088/1748-9326/ 6/1/014005. Toniazzo, T., J. M. Gregory, and P. Huybrechts (2004), Climatic impact of a Greenland deglaciation and its possible irreversibility, J. Clim., 17, 21–33. Trüb, J., and H. C. Davies (1995), Flow over a mesoscale ridge: Pathways to regime transition, Tellus, Ser. A, 47, 502–524, doi:10.1034/j.16000870.1995.t01-2-00008.x. Tsukernik, M., D. N. Kindig, and M. C. Serreze (2007), Characteristics of winter cyclone activity in the northern North Atlantic: Insights from observations and regional modeling, J. Geophys. Res., 112, D03101, doi:10.1029/2006JD007184.

D09107

Wallace, J. M., and D. S. Gutzler (1981), Teleconnections in the geopotential height field during the Northern Hemisphere winter, Mon. Weather Rev., 109, 784–812, doi:10.1175/1520-0493(1981)1092.0.CO;2. Walsh, J. E., W. L. Chapman, V. Romanovsky, J. H. Christensen, and M. Stendel (2008), Global climate model performance over Alaska and Greenland, J. Clim., 21, 6156–6174, doi:10.1175/2008JCLI2163.1. Wild, M., P. Calanca, S. C. Scherrer, and A. Ohmura (2003), Effects of polar ice sheets on global sea level in high-resolution greenhouse scenarios, J. Geophys. Res., 108(D5), 4165, doi:10.1029/2002JD002451. Wild, M., J. Grieser, and C. Schär (2008), Combined surface solar brightening and increasing greenhouse effect support recent intensification of the global land-based hydrological cycle, Geophys. Res. Lett., 35, L17706, doi:10.1029/2008GL034842.

15 of 15