ABSTRACT. Multi-annual records of glacier surface meteorology and energy balance are necessary to resolve glacier-climate interactions but remain sparse, especially in the Southern Hemisphere. To address this, we present a record from the ablation zone of Brewster Glacier, New Zealand, between October 2010 and September 2012. The mean air temperature was 1.2°C at 1760 m a.s.l., with only a moderate temperature difference between the warmest and coldest months (�8°C). Long-term annual precipitation was estimated to exceed 6000 mm a -1 , with the majority of precipitation falling within a few degrees of the freezing level. The main melt season was between November and March (83% of annual ablation), but melt events occurred during all months. Annually, net radiation was positive (a source of energy) and supplied 64% of the melt energy, driven primarily by net shortwave radiation. Net longwave radiation was often positive during cloudy conditions in summer, demonstrating the radiative importance of clouds during melt. Turbulent sensible and latent heat fluxes were directed towards the surface in the summer months, accounting for just over a third of the energy for melt (34%). The energy gain associated with rainfall was small except during heavy events in summer.

The turbulent sensible and latent heat fluxes are important components of the surface energy balance over glaciers in the Southern Alps of New Zealand, contributing over half the energy available for ablation during large melt events. To calculate these terms confidently in glacier mass-balance models it is essential to use appropriate parameterizations for surface roughness and atmospheric stability. Eddy covariance measurements at Brewster Glacier were obtained over an ice surface to help facilitate an assessment of the calculation of the turbulent heat fluxes. The roughness length for momentum was found to be 3.6 x 10−3m, while the roughness lengths for temperature and humidity were two orders of magnitude smaller, in agreement with surface renewal theory. A Monte Carlo approach was used to assess the uncertainty in turbulent heat fluxes calculated using the bulk aerodynamic method. It was found that input-data and roughness-length uncertainty could not explain underestimates of observed sensible heat fluxes during periods with low wind speed and large temperature gradients. During these periods a katabatic wind speed maximum alters the formulation of the turbulent exchange coefficient to that typically observed in a neutral atmosphere and this has implications for glacier mass-balance sensitivity.

ABSTRACT. Recognising the scarcity of glacier mass-balance data in the Southern Hemisphere, a massbalance measurement programme was started at Brewster Glacier in the Southern Alps of New Zealand in 2004. Evolution of the measurement regime over the 11 years of data recorded means there are differences in the spatial density of data obtained. To ensure the temporal integrity of the dataset a new geostatistical approach is developed to calculate mass balance. Spatial co-variance between elevation and snow depth allows a digital elevation model to be used in a co-kriging approach to develop a snow depth index (SDI). By capturing the observed spatial variability in snow depth, the SDI is a more reliable predictor than elevation and is used to adjust each year of measurements consistently despite variability in sampling spatial density. The SDI also resolves the spatial structure of summer balance better than elevation. Co-kriging is used again to spatially interpolate a derived mean summer balance index using SDI as a co-variate, which yields a spatial predictor for summer balance. The average glacier-wide surface winter, summer and annual balances over the period 2005-15 are 2484, −2586 and −102 mm w.e., respectively, with changes in summer balance explaining most of the variability in annual balance.

ABSTRACT:Clouds are important features of many high-altitude and glaciated areas, yet detecting their presence and specifying their effects on incoming shortwave (SW↓), longwave (LW↓) and net all-wave radiation (Rnet) remains challenging in these environments. These limitations hamper efforts to understand atmospheric controls on glacier surface mass balance (SMB) in the Southern Alps of New Zealand, as both cloud and airmass forcing accompanies key synoptic controls on SMB. Multi-year datasets of four-component broadband radiation from two sites at Brewster Glacier, Southern Alps of New Zealand, are used here to develop cloud metrics to account for the effects of clouds on SW↓, LW↓ and Rnet. On average 23% of top-of-atmosphere shortwave radiation (SW TOA ) is attenuated by the clear-sky atmosphere, while clouds attenuate a further 31%, resulting in <50% of SW TOA being received at the surface. The transmission of shortwave radiation by clouds (trc) during overcast conditions is found to vary with season and airmass characteristics. A simple parameterization is developed to account for lower trc observed during periods of higher water vapour pressure. Cloud metrics derived at the site show overcast conditions are frequent (45% of period) and strongly dependent on wind direction, highlighting the dominant role of orography in cloud formation and enhancement in the Southern Alps. The effect of clouds on Rnet exhibits a distinct seasonal variation; during summer when albedo and trc are lower, clouds decrease Rnet by 20-40 W m −2 , while during autumn, winter and spring, clouds enhance Rnet by approximately 20 W m −2 . Idealized modelling shows that these patterns are strongly dependent on albedo and extend across the elevation range of glaciers in the Southern Alps. Thus, overcast conditions appear to aid the extension of ablation into spring and autumn by increasing the energy available for snow and ice melt.