To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5 × 10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.

We present twenty-three years (1993–2016) of automatic weather station (AWS) data, collected along the K-transect near Kangerlussuaq in west Greenland. The transect runs from east to west, roughly perpendicular to the ice sheet edge at about 67° N. The K-transect originated from the Greenland Ice Margin Experiments (GIMEX), held in the summers of 1990 and 1991. Until recently, surface mass balance and ice velocity measurements were performed at nine locations along the K-transect, of which four are equipped with AWS: two in the ablation zone at approximately 500 m and 1,000 m asl, one at the approximate equilibrium-line altitude (~1,500 m asl), and one in the lower accumulation zone (~1,850 m asl) at distances of 5, 38, 88, and 140 km from the ice edge, respectively. Here, we present an overview of the various AWS types and their data corrections, quality, and availability, including a preliminary trend analysis. Recent increases in temperature and radiation components are associated with the frequent occurrence of anti-cyclonic conditions in west Greenland, resulting in clear skies and relatively warm summers. Strong melt concurs with a decrease in winter accumulation, lowering the surface albedo of the ice sheet. The AWS situated at 1,500 m asl, the former equilibrium-line altitude (ELA), observed almost a doubling of the summertime net shortwave radiation since 2004; as a result, the ELA along the K-transect has been steadily increasing and is currently situated well above 1,700 m asl.

We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z _s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z ₀). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z ₀ below a threshold value of approximately 10⁻³ m, the model performs well and in accord with earlier studies. However, with growing hummock size, z ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z _s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z ₀ > 10⁻³ m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z _s /z ₀) = 1.5 − 0.2 ln (Re _*) − 0.11(ln (Re _*))².

Long term methane flux measurements have been mostly performed with plant or soil enclosure techniques on specific components of an ecosystem. New fast response methane analyzers make it possible to use the eddy covariance (EC) technique instead. The EC technique is advantageous because it allows continuous flux measurements integrating over a larger and more representative area including the complete ecosystem, and allows fluxes to be observed as environmental conditions change naturally without disturbance. We deployed the closed-path Fast Methane analyzer (FMA) from Los Gatos Research Ltd and demonstrate its performance for EC measurements at a Ponderosa pine plantation at the Blodgett Forest site in central California. The fluctuations of the CH₄ concentration measured at 10 Hz appear to be small and their standard deviation is comparable to the magnitude of the signal noise (±5 ppbv). Consequently, the power spectra typically have a white noise signature at the high frequency end (a slope of +1). Nevertheless, in the frequency range important for turbulent exchange, the cospectra of CH₄ compare very well with all other scalar cospectra confirming the quality of the FMA measurements are good for the EC technique. We furthermore evaluate the complications of combined open and closed-path measurements when applying the Webb-Pearman-Leuning (WPL) corrections (Webb et al., 1980) and the consequences of a phase lag between the water vapor and methane signal inside the closed path system. The results of diurnal variations of CH₄ concentrations and fluxes are summarized and compared to the monthly results of process-based model calculations.

The southward progression of ice shelf collapse in the Antarctic Peninsula is partially attributed to a strengthening of the circumpolar westerlies and the associated increase in föhn conditions over its eastern ice shelves. We used observations from an automatic weather station at Cabinet Inlet on the northern Larsen C ice shelf between 25 November 2014 and 31 December 2016 to describe föhn dynamics. Observed föhn frequency was compared to the latest version of the regional climate model RACMO2.3p2, run over the Antarctic Peninsula at 5.5-km horizontal resolution. A föhn identification scheme based on observed wind conditions was employed to check for model biases in föhn representation. Seasonal variation in total föhn event duration was resolved with sufficient skill. The analysis was extended to the model period (1979–2016) to obtain a multidecadal perspective of föhn occurrence over Larsen C ice shelf. Föhn occurrence at Cabinet Inlet strongly correlates with near-surface air temperature, and both are found to relate strongly to the location and strength of the Amundsen Sea Low. Furthermore, we demonstrated that föhn occurrence over Larsen C ice shelf shows high variability in space and time.

The ice surface temperature (IST) derived from thermal infrared remote sensing is crucial for accurately monitoring ice or snow surface temperatures in the polar region. Generally, the remote sensing IST needs to be validated by the in situ IST to ensure its accuracy. However, due to the limited availability of in situ IST measurements, previous studies in the validation of remote sensing ISTs are scarce in the Antarctic ice sheet. This study utilizes ISTs from eight broadband radiation stations to assess the accuracy of the latest-released Moderate Resolution Imaging Spectroradiometer (MODIS) IST and Visible Infrared Imager Radiometer Suite (VIIRS) IST products, which were derived from two different algorithms, the Split-Window (SW-based) algorithm and the Temperature–Emissivity Separation (TES-based) algorithm, respectively. This study also explores the sources of uncertainty in the validation process. The results reveal prominent errors when directly validating remote sensing ISTs with the in situ ISTs, which can be attributed to incorrect cloud detection due to the similar spectral characteristics of cloud and snow. Hence, cloud pixels are misclassified as clear pixels in the satellite cloud mask during IST validation, which emphasizes the severe cloud contamination of remote sensing IST products. By using a cloud index (n) to remove the cloud contamination pixels in the remote sensing IST products, the overall uncertainties for the four products are about 2 to 3 K, with the maximum uncertainty (RMSE) reduced by 3.51 K and the bias decreased by 1.26 K. Furthermore, a progressive cold bias in the validation process was observed with decreasing temperature, likely due to atmospheric radiation between the radiometer and the snow surface being neglected in previous studies. Lastly, this study found that the cloud mask errors of satellites are more pronounced during the winter compared to that in summer, highlighting the need for caution when directly using remote sensing IST products, particularly during the polar night.