JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109, D16105, doi:10.1029/2003JD004451, 2004
Greenland ice sheet surface mass balance 1991––2000: Application of Polar MM5 mesoscale model and in situ data Jason E. Box and David H. Bromwich Polar Meteorology Group, Byrd Polar Research Center, Ohio State University, Columbus, Ohio, USA
Department of Geography, Atmospheric Sciences Program, Ohio State University, Columbus, Ohio, USA
Le-Sheng Bai Polar Meteorology Group, Byrd Polar Research Center, Ohio State University, Columbus, Ohio, USA Received 15 December 2003; revised 30 April 2004; accepted 12 May 2004; published 26 August 2004.
[1] The Polar Pennsylvania State University–National Center for Atmospheric Research Fifth-Generation Mesoscale Model (Polar MM5) regional climate model was run over the North Atlantic region for 1991–2000. We analyze 24-km output over the Greenland ice sheet to evaluate spatial and temporal variability of the surface mass balance and its subcomponents. The model output is compared with 3 years of automatic weather station (AWS) data from 17 sites to identify biases. Using the in situ data, we derive simple corrections for biases in melt energy and in water vapor fluxes from the surface and from blowing snow. The simulated accumulation rate is in agreement with AWS and snow pit observations. Estimates for runoff and the surface mass balance distribution over the ice sheet are produced using modeled melt volume and a meltwater retention scheme. From the decade investigated, the magnitude of interannual variability in surface mass balance components is tentatively established. The largest variability is concentrated along the ice sheet margin, where both accumulation and ablation rates are largest. The simulated interannual fluctuations suggest a large absolute variability, ±187 km3 yr 1 for total ice sheet surface mass balance. Variability in simulated equilibrium line altitude is suggestive of a dominance of thermal variability in the south with increasing importance of accumulation variability with increasing latitude. Empirical functions for the sensitivity of surface mass balance to temperature and precipitation anomalies are presented. The precise locations and regions of maximum and minimum surface energy and mass fluxes are suggested. Using an estimate for iceberg discharge and bottom melting, the total ice sheet mass balance is estimated be 78 km3 yr 1, producing 2.2 mm of eustatic sea level rise over the 1991–2000 decade and contributing 15% to the observed (1.5 mm yr 1) global sea level rise. The more negative mass balance is attributed to including blowing- snow sublimation loss and to regional warming in the 1990s. INDEX TERMS: 1620 Global Change: Climate dynamics (3309); 1827 Hydrology: Glaciology (1863); 3360 Meteorology and Atmospheric Dynamics: Remote sensing; 3349 Meteorology and Atmospheric Dynamics: Polar meteorology; 4556 Oceanography: Physical: Sea level variations; KEYWORDS: Greenland ice sheet, mass balance, runoff Citation: 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.
1. Scientific Background balance components are imprecisely known [van der Veen, 2002]. Therefore ranking the importance of the Greenland ice [2] Ice sheet mass balance plays an important role in global sheet in global ocean dynamics has been associated with sea level and thermohaline circulation changes. Greenland large uncertainty. Contributing to this uncertainty, surface ice sheet mass balance changes appear to have contributed mass balance estimates have been produced on the basis of a several meters to sea level fluctuations since the last climatic temporally nonuniform set of observations [e.g., Ohmura et optimum, 125,000 years ago [Cuffey and Marshall, 2000], al., 1999]. The interest in Greenland ice sheet mass balance and are expected to contribute to sea level rise under is particularly relevant given the dramatic thinning rates projected future global warming throughout this century observed in the 1990s at low elevations [Krabill et al., [Church et al., 2001]. However, Greenland ice sheet mass 2000]. The thinning is the combined result of ice flow changes and surface ablation. Apparently critical to this Copyright 2004 by the American Geophysical Union. thinning is the interplay between water and ice dynamics, 0148-0227/04/2003JD004451$09.00 providing a mechanism for accelerated ice sheet response to
D16105 1of21 D16105 BOX ET AL.: GREENLAND ICE SHEET SURFACE MASS BALANCE D16105 climate warming [Zwally et al., 2002]. Data and models both have tested a coupled RCM-snow model using a detailed set suggest that abrupt climate change during the last glaciation of observations from the ETH camp in western Greenland. originated through changes in the Atlantic thermohaline Using the Mode`le Atmosphe´rique Re´gional (MAR) RCM, circulation in response to small changes in the hydrological they produce a significantly better treatment of the near- cycle [Clark et al., 2002]. The other ice sheet, in Antarctica, is surface state variables than in the relatively coarse European 9 times larger in volume, yet only Greenland exhibits sub- Centre for Medium-Range Weather Forecasts (ECMWF) 15- stantial seasonal melting, with runoff representing approxi- year reanalysis (ERA-15). Galle´e and Duynkerke [1997] mately half the ablation total [Janssens and Huybrechts, 2000; employed an early version of MAR to investigate modeling Zwally and Giovinetto, 2000] and the remainder of ablation surface energy and mass balance over a transect of observa- coming from surface water vapor fluxes [Box and Steffen, tional stations for the Greenland Ice Margin Experiment 2001; De´ry and Yau, 2002] as well as iceberg discharge and (GIMEX) 1990–1991, citing limited success in simulating subglacial melting [Rignot et al., 1997; Reeh et al., 1999]. basic meteorological variables and surface melting and run- Surface water vapor mass fluxes are generally small in off. Lefebre et al. [2003] have validated and applied an comparison with precipitation rates; however, the associated evolving albedo model applicable by RCMs. Box and Rinke latent heat flux represents an important energy sink, without [2003] used the HIRHAM RCM to evaluate multiyear which surface melt rates would be significantly larger. simulations of meteorological and mass balance components [3] Global climate models (GCMs) have been used in over Greenland. The Polar MM5 has been run over the 1991– experiments to simulate ice sheet mass balance under future 2000 period to evaluate precipitation in an 8-km nested climate scenarios to relate climate fluctuations to mass domain over Iceland (D. H. Bromwich et al., High resolution balance variability [e.g., Ohmura et al., 1996; Thompson regional climate simulations over Iceland using Polar MM5: and Pollard, 1997; Glover, 1999; Wild and Ohmura, 2000; 1991–2000, submitted to Monthly Weather Review, 2003). In Church et al., 2001; van de Wal et al., 2001; Murphy et al., the present study, we take advantage of 10 years of data 2002; Wild et al., 2003]. Collectively, the GCM results availability over Greenland at 24-km horizontal resolution produce a large envelope of mass balance estimates as a generated in support of the Iceland study. result of the different physical parameterizations used. A [6] Hanna et al. [2002] and Mote [2003] have made major limitation of GCMs is coarse horizontal resolution, significant progress in incorporating RCM, in situ, and currently 100–250 km. At this spatial scale the steep ice remotely sensed data to estimate Greenland ice sheet surface sheet margin and ablation zone are not adequately resolved. mass balance and runoff. In the latter study, for a 12-year Further, terrain smoothing along steep ice margins can period (1988–1999), there was some evidence that larger produce net elevation lowering and results in systematic runoff is associated with lower accumulation. To derive overestimation of ablation rates from the resulting positive surface mass balance and runoff, Mote [2003] employed a bias in sensible heat flux [e.g., Glover, 1999; Hanna and methodology depending on the following: a constant thresh- Valdes, 2001]. Wild et al. [2003] have interpolated GCM old in daily average passive microwave brightness tempera- temperature fields to a fine (2-km) grid to estimate ablation ture to indicate melting; constant degree day factors for ice or rates using an empirical model. In a departure from previous snow melt; precipitation distributed evenly throughout the results, their projection implies a future dominance of year; constant spatial distribution of 1995–2000 average accumulation over ablation in a climate-warming scenario. sublimation/evaporation rate after Box and Steffen [2001]; [4] Regional climate models (RCMs), operating on hor- constant meltwater retention fraction; and liquid precipitation izontal scales of tens of kilometers, offer an attractive not added to the accumulation rate, but assumed as direct alternative to GCMs, automatic weather station networks, runoff. We pursue a similar methodology, however, with the and statistical climatologies. RCM applications to the following differences: calculation of meltwater production observed past rather than the future benefit from an assim- based on surface energy balance closure, variable meltwater ilated set of observations from surface weather stations, retention fraction, variable water vapor fluxes from the atmospheric soundings, and satellite remote sensing, i.e., surface and blowing snow, and liquid precipitation added to from atmospheric analyses. As such, RCMs act as physi- melt volume and accumulation. cally based interpolators that may be used to address persistent questions in better understanding the magnitude of climate variability over ice sheets, questions nearly 2. Objectives impossible to answer, even with much larger surface obser- [7] The first objective of this study is to reassess the vational networks. RCMs have the benefit of providing accuracy of the Polar MM5 RCM, under a different model climate data for explicit time periods, as opposed to tem- configuration, i.e., from that of Bromwich et al. [2001a] and porally nonuniform climatologies. Greenland is ideally Cassano et al. [2001], using in situ observations collected suited for atmospheric circulation modeling, with prevailing by a network of automatic weather stations (AWS) and large-scale flow across the domain, unlike Antarctica, where supplemented by other glaciological observations. Given the flow is quasi-concentric; thus less benefit comes from this assessment, estimates of surface mass balance compo- lateral constraints imposed by the atmospheric analyses. nents may be produced with better understood uncertainties. [5] Ice sheet studies with RCMs are increasingly common. Adjusting mesoscale model output based on comparisons Cassano et al. [2001] and Bromwich et al. [2001a] made a with AWS data is explored in an effort to produce more validation study of the Polar Pennsylvania State University– reliable surface mass balance estimates. By analyzing mul- National Center for Atmospheric Research Fifth-Generation tiple years of RCM output, the variability in surface mass Mesoscale Model (Polar MM5) RCM over Greenland using balance components may be gauged. Further, by incorpo- automatic weather station (AWS) data. Fettweis et al. [2003] rating interannual surface mass balance results with iceberg
2of21 D16105 BOX ET AL.: GREENLAND ICE SHEET SURFACE MASS BALANCE D16105
[11] Polar MM5 was configured to employ a K-theory bulk method for turbulent fluxes (medium-range forecast planetary boundary layer (MRF PBL) scheme), which corrects for missing local gradients in large PBL-scale eddies using a countergradient flux after Troen and Mahrt [1986]. The length scales for temperature and humidity were selected to be equal to that for momentum (0.1 mm), although they are probably some fraction of that for momentum [Garratt, 1994]. Grid and subgrid precipitation and cloud processes are represented by the same configu- ration as given by Bromwich et al. [2001a]. Cassano et al. [2001] includes details for the model options that were selected, including the following: the NCAR Community Climate Model, Version 2 (CCM2), radiation code [Hack et al., 1993] with the Meyers parameterization for ice nuclei concentration; the Reisner mixed-phase cloud parameteri- zation; the two-way nest that allows interaction of atmo- spheric motion and heat transfer across model domains; and the ‘‘multilayer’’ soil model. The rigid lid upper boundary condition was used. Surface temperature was calculated Figure 1. Polar MM5 model domains defined in a 10-year from the surface energy budget. Sea surface temperature simulation (1991–2000). Data from domain 2 are analyzed was invariant over each daily simulation. Solar radiation in this study. Domain 3, 8-km grid resolution is not shown. was computed hourly. Model grid cells are represented as homogeneous, i.e., no mixed land/ice points. discharge and basal melting estimates, we assess total ice [12] The Polar MM5 was initialized at the surface and sheet mass balance and provide estimates for the contribu- updated every 12 hours at the lateral boundaries with the tion of Greenland ice sheet mass balance to sea level. 2.5 horizontal resolution ECMWF operational analyses. [8] Sections 3.1–3.6 describe the atmospheric model, Although these analyses have the disadvantage of being techniques to apply the model data to questions of surface affected by inhomogeneities due to modifications in the mass balance, the observational data, and methods used for ECMWF assimilation/modeling system over time as well as model validation and correction. Maps of the mean surface changes in the global observing system, they were selected mass balance components and related meteorological param- because the ERA-15 does not cover when AWS measure- eters are then presented and discussed in light of previous ments are available and the relatively homogeneous ERA- results and specific regional controls. Integrated total ice 40 data were unavailable. sheet surface mass fluxes, equilibrium line variability and [13] The Polar MM5 was run at 24-km horizontal sensitivity, and the sensitivity of surface mass balance to resolution in a series of 30-hour forecasts. The first 6 hours temperature and precipitation anomalies are presented in were discarded for model equilibration, and the remaining succession. The paper concludes with recommendations 24 hours were concatenated into a 6-hourly output time about the state of mass balance of the Greenland ice sheet series spanning 1991–2000. With 6-hourly model output it and the challenge of reducing persistent uncertainties. is arguable that the diurnal cycle is inadequately resolved, e.g., for melting. Output data are available for 0600, 1200, 1800, and 2400 UTC (or 0400, 1000, 1600, and 2200 west 3. Methods Greenland standard time). Although solar noon is not 3.1. Atmospheric Model represented, data near the warmest and coldest times of [9] The Pennsylvania State University–National Center day are available. We integrate the 6-hourly data into for Atmospheric Research (NCAR) Fifth-Generation Meso- annual distributions of surface mass balance components, scale Model (MM5) has been modified for use in polar i.e., precipitation P, surface water vapor flux E, blowing- regions [Bromwich et al., 2001a; Cassano et al., 2001]. The snow sublimation QS, and runoff R. The surface mass major polar-specific modifications are as follows: explicit balance was calculated as P-E-QS -R. Model bias versus ice phase cloud microphysics, polar physical parameter- time into the simulation was investigated, as a short izations for cloud-radiation interactions, stable boundary (6-hour) spin-up is suspected to lead to model bias. layer turbulence parameterization, improved treatment of Earlier simulations used 24-hour spin-up time [i.e., heat transfer through snow and ice surfaces, and implemen- Cassano et al., 2001]. Shorter spin-up time was selected tation of a sea ice surface type. for computational efficiency for this 1.6-times-higher- [10] The model domains configured for this study are resolution and 10-times-longer-duration simulation includ- shown in Figure 1. Each domain had 28 vertical levels. The ing a high-resolution (8-km) model subdomain (Iceland) lowest two model levels correspond to roughly 14 m and not analyzed here. Of note is that the model configuration 0 m. The model is run in nonhydrostatic mode with a sixth- was not optimized for Greenland; however, this study was order finite-difference scheme after Chu and Fan [1997] pursued given the opportunity to extract useful informa- implemented to improve estimation of vertical motion tion from the simulation byproduct over Greenland. The significant for precipitation formation, particularly along ice sheet area in the 24-km model domain is 1.691 steep slopes. 106 km2 and includes 2936 grid cells.
3of21 D16105 BOX ET AL.: GREENLAND ICE SHEET SURFACE MASS BALANCE D16105
3.2. Surface Energy Balance and Meltwater tant in the surface mass budget over snow and ice Production surfaces [Pomeroy and Essery, 1999; Bintanja, 2001; [14] Changes in surface and subsurface temperatures are De´ry and Yau, 2002]. In the similar environment over directly linked to imbalances in the surface energy budget. Antarctica, wind-driven sublimation represents one of the In the case of a snow/ice surface, when the melting point major sources of uncertainty in surface mass balance has been reached during a period of energy surplus, the closure [e.g., Turner et al., 2002]. Here, we estimate excess energy goes into conversion of solid ice into liquid the blowing-snow sublimation component using a param- water. As such, meltwater production can be calculated eterization from De´ry and Yau [2001], equating QS with from the residual (QM) of the surface energy budget: 10-m wind speed, 2-m air temperature, humidity, and pressure. This parameterization performs similarly to an independent model developed for Antarctica [Mann et al., QM ¼ QN ðÞQH þ QE þ QG þ QR ; ð1Þ 2000], compares well with results over a Canadian Arctic land surface, and therefore appears to be applicable over where QN is the surface net radiative flux. Problematic for melt simulations, the albedo over the ice sheet in this Greenland. The parameterization of Bintanja [1998] pro- realization of Polar MM5 was fixed at 0.8. To account for duced similar results but was not selected owing to this, daily satellite-derived albedos [Key et al., 2002] were uncertainties stated by Bintanja [2001]. applied in surface energy balance calculations of absorbed [16] We evaluate the magnitude of the blowing-snow solar irradiance. Although the temperature-albedo feedback transport horizontal divergence as De´ry and Yau [2002] is not considered here, sensitivity studies using Polar MM5 and Guo et al. [2003] have done for Antarctica, but we indicate that albedo errors do not greatly affect near-surface explicitly account for snow availability. Potential blowing- air or ice temperatures [Cassano et al., 2001]. Year 2000 snow transport in the lowest 10 m (QTP) is estimated using albedo values were unavailable; thus monthly climatological the bulk formulation of Tabler [1991] based on Antarctic values compiled from the 1982–1999 data were used then. blowing-snow measurements [Budd et al., 1966] QH and QE are the turbulent sensible and latent heat fluxes, respectively. Q is the firn/ice conductive heat flux, u3:93 G Q ¼ 10 ; ð4Þ incorporating density-dependent effective thermal conduc- TP 290;951 tivity after Sturm et al. [1997]. A depth-density relationship based on Greenland snow pits is taken from Box [2001], and where u10 is the 10-m wind speed given by the model and 1 1 the conductivity estimates are applied only off-line to the QTP has units of kg m s perpendicular to the wind. model output to estimate QG. Although we expect potentially Snow transport is taken to occur when the wind speed large ‘‘soil’’ temperature errors owing to insufficient spin-up exceeds the temperature-dependent threshold [Li and time and biases in climatological soil temperature initializa- Pomeroy, 1997]. A constant threshold of 7 m s 1 was used tion, we have verified that QG represents a small component for temperatures below 27 C owing to unrealistic in the modeled surface energy budget and is of tertiary behavior of the quadratic function in that extreme. Actual importance to melting as compared to radiation and turbulent snow transport QTA is estimated using a reduction of QTP 0 fluxes. QR is the sensible heat flux from rain. QR may be accounting for snow availability A , which is assumed to be calculated using the rain Tr and surface temperature Ts unlimited when snow is freshly deposited and decays with time from snow events defined by the model. QR ¼ r cwRTðÞr Ts ; ð2Þ w ÀÁ 1 A0 ¼ 1:038 þ 0:03758t 0:00014349t2 þ 1:911315e 7t3 where rw is the density of water, cw is the specific heat of 1 1 1 ð5Þ water (4.2 kJ kg K ), and R is the rainfall rate (m s ). Tr is the rain temperature, and Ts is the surface skin temperature. 0 The contribution of QR to the ablation rate is assessed here Time t is in hours after a snowfall event. After 12 days, A since parts of the ice sheet, particularly in the south, have a decays to a value of 0.22 that is then held constant and is maritime character [Putnins, 1970]. Volumetric melt M in proportional to the force required to disaggregate bonds of 1 cubic meters or surface height variation (dz dt ) is related snow particles of that age measured by Jellinek [1957]. QTA 0 directly to QM when surface temperature is at or above the is estimated by multiplying values of QTP by A . Blowing- melting point, Ts >0 C, as snow divergence D is estimated by differencing snow transport at model grid points with the surrounding grid 1 points and represents the combined effect of blowing-snow M ¼ QM tLðÞr ; ð3Þ sublimation and snow redistribution; thus it cannot be used to distinguish blowing-snow sublimation alone. where t is time (s), r is ice density (917 kg m 3), and L is the latent heat of fusion (384 kJ kg 1). 3.4. Meltwater Retention and Runoff [17] To gauge runoff, we use the Pfeffer et al. [1991] model 3.3. Blowing-Snow Sublimation for meltwater retention, which includes both capillary tension [15] It is important to note that sublimation/evaporation and latent heating effects for a simple uniform wetting front. derived from bulk atmospheric profile methods, i.e., that This method considers meltwater retention for single annual based on QE, corresponds to surface water vapor flux (snow accumulation season) layers only and thus does not only and does not include physics for sublimation from consider the case of meltwater percolating into the previous blowing and drifting snow QS, now recognized as impor- accumulation layer. Consequently, in extremely warm years,
4of21 D16105 BOX ET AL.: GREENLAND ICE SHEET SURFACE MASS BALANCE D16105 runoff may be overestimated, as annual melt exceeding the accumulation rate is not retained. Janssens and Huybrechts [2000] evaluated this and simpler models for meltwater retention using climatological average accumulation and melt. Annual potential meltwater retention fraction pr, esti- mated using the annual mean temperature T,snowmelt volume M, precipitation P, surface sublimation/evaporation E, blowing-snow sublimation QS, specific dry snow accu- mulation C calculated as P-E-QS-liquid precipitation, and the depth of melting d, here is taken as the water equivalent depth of the annual accumulation rate c d C M re pr ¼ T þ 1 ; ð6Þ Lf P P ro